CN1181808A - 确定地磁场方向的方法 - Google Patents

Abstract

本发明给出一种采用一个包括三个磁场传感器和两个倾斜采集设备的电子磁针罗盘确定可以受到固定在仪器上的磁性材料和若干电流的磁场干扰的地磁场方向的方法,其中,电子磁针罗盘置于不同的空间位置N,其中,在此N位置的每一个位置上测量倾斜传感器信号和磁场传感器信号并且由这些信号确定倾斜值和磁场值。在这些倾斜值和磁场值的基础上按照矢量方程确定地磁场矢量值。const(恒定值)=b_g=b_Esin(i)=g^Tb_E=g^Tm(b_gem－b_o)其中b_g=与万有引力矢量g平行的磁场分量b_E=磁场矢量b_E的长度b_gem=所测量的磁场矢量b_E=在测量地点实际的地磁场矢量△M=软磁场失真,即通过地磁场感应的磁性M=I+△M;I=单位矩阵b_o=硬磁干扰场矢量i=磁偏向。其中,N至少应与矢量方程待确定的参数的数量相等。

Description

确定地磁场方向的方法

本发明涉及一种采用一个包括三个磁场传感器和两个倾斜采集设备的电子磁针罗盘确定可以受到与仪表固定的磁性材料和若干电流磁场干扰的地磁场方向的方法。

在US4 686 772中公开了一种例如能够确定装甲车行驶方向的电子磁针罗盘。装甲车具有两个铁制物体，即炮塔和车体，它们相互可以绕垂直轴转动。磁针罗盘安装在车体上，并且包括一个不摆动的三轴磁力计。此磁力计给出表示磁力计所在部件的三个磁场分量的电磁信号。配置了两个倾角传感器，它们给出表示具有倾角传感器的车体的抑角和滚动角的电信号。此外，角度测量设备给出一个表示两个铁制物体之间的角度的转角信号。在一个存储中存储着多个预先校准过的修正系数，这些修正系数应修正由车辆感应的磁场对磁力计对两个铁制物体之间的转角测量的影响。借助于一个计算机由电磁场信号和倾角传感器信号实时计算装甲车行驶方向，这些信号通过储存的修正系数对相应的转角已进行了修正。

在能够应用所述磁针罗盘之前有必要进行校准测量。这里此车辆置于至少两个不同的、相互之间不平行的平面上。在这些部位上在车位不同的排列但是在装甲炮塔相对相同的排列下进行方位角、抑角和滚动角的测量。为了能够测量这些角度，例如需要一个经纬仪、在校准时假定在测量部位的地磁场是已知的。这里由卡片材料可以获得相应测量部位地磁场的磁倾角和磁偏角的数值。

在US4 539 760中描述了一个具有三个磁传感器的电子磁针罗盘。这些传感器对包括地磁场和一个附加的与车辆有关的干扰场的磁场的三个正交分量作出反应。这些磁传感器产生与这些分量相应的电信号。此外，若干倾角传感器对车辆与水平平面的倾斜角作出反应。采用一个数据处理设备和一个存储器存储这些作为校准修正值的、由测量传感器导出的信号，如果车辆为了磁针罗盘的校准而转动了一个圆。为了避免干扰磁场的影响，数据处理设备在校准过程结束后，在考虑到校准修正值的情况下计算在车辆所在部位地磁场的修正值，然后由这些修正值出发在应用用于确定地平线的倾角传感器数值的情况下计算车辆行驶方向的方位。在计算时由修正矩阵是对称的这一点出发，尽管这只是在非常的情况下才能出现的。

本发明的任务在于给出一种采用电子磁针罗盘确定地磁场方向的比较简单的方法。

按照本发明此项任务是按照权利要求1所述特征加以解决的。

本发明的进一步发展在各项从属权利要求中作了阐述。

以有利的方式按照本发明此方法不要求采用附加的测量仪器进行特定的校准测量。此外，在电子磁针罗盘中也无需特定数据的输入，以便于在应用中进行推广。唯一需要的是将此磁针罗盘置于若干不同的任何空间位置。在每一个空间位置优先确定三个磁场分量。这些磁场分量可以相互正交，如果提出这种要求的话。但是这并不是一定要这么作的。借助于在每个空间位置附加测量的倾斜值、即倾角和滚动角、可以由相应的磁场分量中确定地场实际磁场矢量的方向。

本发明方法不要求电子磁针罗盘在安装前必须进行校准。采用本发明方法不仅考虑到了干扰磁场，而且还考虑到了制造容差以及各个传感器具有不同的灵敏度等。此外，有必要采用一个在出厂前已由厂家校准的磁针罗盘。

在本发明方法中和用了以下事实，即在万有引力矢量和地磁场矢量之间的偏角在相应的固定部位均与系统的瞬时位置无关而保持恒定值。

所给出的附图是为了便于理解本发明：

图1示出了电子磁针罗盘和软磁和磁磁干扰场发生器的一个装置的示意图；

图2示出与测量地磁场相关的矢量。

波松早在19世纪末就探讨了这个问题，即如何确定实际存在的磁场，如果在系统中放置了一个磁力针，而它本身具有磁场分量。描述这一情况的波松方程指出，在这样一种情况下所测量的磁场是实际存在场的一个线性函数，它涉及到一个仿射性图象。在这方面请参阅上述US4686 772第2部分第26至30行。

对于一般情况可以说所测量的磁场是由测量位置被软磁失真了的地磁场和一个硬磁分量构成的。在被软磁失真了的地磁场中涉及到一个通过地磁场感应的磁性。硬磁分量包括例如通过永磁铁或在系统中的若干电流产生在磁力计部位恒定的磁场。硬磁分量不受外部场改变的影响。

在数学上上述波松方程稍加变更后可以表示为：

b_gem＝Mb_E+b_o＝b_E+ΔMb_E+b_o    (1)其中

b_gem＝所测量的磁场矢量

b_E＝在测量部位实际的地磁场分量

ΔM＝软磁场失真，就是说通过地磁场感应的磁性，其中

M＝I+ΔM；I＝单位矩阵

b_o＝硬磁分量

以上公式以及下面的表达式中矢量和矩阵用黑体符号表示。这些矢量一般是参照一个笛卡儿系统。矩阵一般涉及至3×3矩阵。

在图1中的示意图表示在仪器中安装了磁针罗盘DMC，在此仪器中作为硬磁干扰器的永磁铁与软磁材料相互作用并且影响DMC的测量值。在图2中示出与地磁场测量有关的矢量及其在水平面上的投影。在固定安装在车辆上的磁针罗盘中目标方向相当于行驶方向。

电子磁针罗盘具有三个磁场传感器和两个倾斜传感器，其中，前者确定三个磁场分量，这些磁场分量不必相互正交。

永磁铁和若干电流在磁场传感器部位产生一个固定场，此固定场表示由磁场分量形成的坐标系统的零点偏移。

软磁材料由已存在的场在行驶方向产生一个减弱或加强了的场以及在与其垂直方向产生附加的场分量。这可以看作是场方向X、Y和Z的“串扰”。

如果观察一个磁力针，此磁力计具有三个不同灵敏度的三个非正交磁场传感器和零点偏移，例如一个由厂家提供的未经校准的“原始”磁力计，则也可以得出上述相同的方程式。

采用这样一个“原始”磁力针进行磁场测量在数学上可以以如下方式表示：

b_gem，i＝f_ie_ib_Ei十o_i；i＝x，y，z其中f_i为增强，e_i为测量方向，即单位矢量，o_i为I-传感器的分支。

如果写成f_ie_i＝(M_i1，M_i2，M_i3)和o_i＝bo_i，则可以重新得到以上方程(1)。

为了确定实际地磁场的大小，有必要从上面给出的矢量式(1)中解出b_E。

通过转换和减法得出：

b_E＝M^-1(b_oem-b_o)＝m(b_gem-b_o)                         (2)其中m＝M^-1。

为由方程(1)、(2)中确定未知量M或m＝M^-1以及b_o可以采用上述US4 686 772中公开了的求解途径。在此文献中的出发点是在测量的磁场矢量b_gem中每一次还清楚的知道地磁场矢量b_E。矢量方程(1)和(2)相应表示未知的M_ij以及m_ij和b_oi的一个线性方程线。通过至少在4个不同几何位置的测量在已知地磁场情况下可以借助于求解线性方程组的基本方程直接求得M_ij和b_oi。

在上述US4 539 760中给出的结果是从以下事实出发的，即地磁场的值与磁力计所在部位不相关。如果正确的确定了m_ij和b_o，则可以得出地磁场矢量b_E，此矢量在磁力计的每一个部位均具有相同的长度。于是求得：

const＝b_E ²＝b_E ^Tb_E＝(b_gem-b_o)^Tm^Tm(b_oem-b_o)

          ＝(b_gem-b_o)^TU(b_gem-b_o)

                                                     (3)其中U^T＝U，这意味着U是一个对称矩阵。

const为恒定值。

由此不难看出，由此方程(3)只能计算所求矩阵m的乘积U＝m^Tm。只有假定此矩阵是对称的并且假定对角线为正时才能计算此矩阵的各项。第一种假定只有在很少的情况下才能出现，因为有可能意味着是一个软磁对称，这种对称在技术设备中，例如在车辆或汽车上是极其不大可能出现的。

现在按照本发明为了确定b_E值，就是说解矢量方程由列式出发，其中利用了如下事实，即在测量系统的每一个位置在相同的地理位置在水平面和地磁场之间的角度、即磁倾角是恒定的。这一点显然对万有引力矢量g的方向和地磁场矢量b_E之间的角度也成立。这样可以写为：const＝b_g＝b_Esin(i)＝g^Tb_E＝g^Tm(b_gem-b_o)

                                                     (4)其中，const为恒定值，

b_g＝平行于万有引力矢量g的磁场分量

b_E＝磁场矢量b_E的长度

i表示磁倾角。

此方程表明，在地重力场方向的磁场分量，就是说垂直于水平面，在所有的系统部位均相等。

在此方程中m值是线性的并且不是一个乘积。万有引力矢量g可以借助倾斜传感器确定。所以m值可以直接计算出来而无需测量场，如上述US4 686 772所指出的，或者应假定特定的对称条件，如上述US 4 539760所指出的。

待确定参数的数量为：

m＝3×3＝9

b_o＝3

b_g＝1

这样我们就总共得出了9+3+1＝13个参数，确定这些参数至少需要13个方程。

还可以选择比例系数，例如在US4 539 760中(第4栏第三行)公开的那样。可能的一些比例确定例如可参考下表：

b_g＝const(恒定值)，

m₁₁＝const(恒定值)；

m₁₁+m₂₂+m₃₃＝const(恒定值)

m₁₁ ²+m₂₂ ²+m₃₃ ²＝const(恒定值)

m₁₁ ²+m₁₂ ²+…+m₃₃ ²＝const(恒定值)

det m＝const(恒定值)

或者选择其它适当的表格，其中恒定值也可以选择为1。通过确定比例参数的数量可减少1个，即现在只还剩12个参数，即还需要相同数量的方程。在磁针罗盘初始校准时需占有12个不同的几何位置j＝1，……12，在这些位置上测量三个磁场分量和两个倾斜角。

现在得到一个线性的方程组：

b_g＝g_i ^Tm(b_gemj-b_o)＝g_j ^Tmb_gemj-g_j ^Tu_o                 (5)其中U_o＝mb_o并且例如b_g＝1，其中明确得出：

1＝g_j1b_gem，1m₁₁+g_j1b_gem，2m₁₂+...

...+g_j1u_o1+g_j2u_o2+g_j3u_o3

如果提供了12个以上的方程式，则这里例如可以与最小平方误差最佳的适配。这是大家知道的。

在本发明范围内还提供了利用在具有磁场矢量长度为恒定值(b_E＝|b_E|)的同一固定位置上磁偏角以及b_g为恒定值的事实的可能性。换言之，这就意味着不仅上面的方程(5)，而且方程(3)也可以应用在本方法中。然后还需要附加的确定参数b_E，这样参数的数量并且从而所必要的测量提高了1次、即N＝13次。

这里的优点在于通过应用方程(3)和方程(5)可以更好的利用所提供的数据，因为在两个方程中要用到每一次测量值。因此必要的测量数量只有一半。为了实现此方法要求进行统计平衡计算，为此例如可以最小平方误差方法作为基础。

在下述方法中上面的方程(3)和(5)以统计精确方法加以利用。

下面将进一步阐述磁场传感器信号M_ij；i＝1，2，3和地磁场矢量b_Ej之间形成的关系。此关系可以通过以下方程表示：μ_j＝μ_o+Mb_Ej+ε_j；j＝1，.N       (测量部位)

b_Ej＝地磁场矢量

ε_j＝表示传感器噪声的偶然矢量

<ε_j>＝0；<ε_jε_j ^T>＝σ_μ ²l₃；<ε_jε_k ^T>＝O j＝k；

<>＝统计预期值

地磁场矢量b_Ej的各分量在磁针罗盘的不同部位是未知的，但是可以借助于倾斜传感器部分地加以表示。其中，两个倾斜传感器表示归一化的万有引力场矢量的各一个分量。

b_Ej＝sg_j+cn_j

s＝sin(i)；c＝cos(i)

|b_Ej|＝1  地磁场强度＝1

i＝地磁场的磁偏角

η_j＝cos(α_j)e_j+sin(α_j)f_j；其中α_j为方位角，就是说，在水平面上传感器自身的坐标系统与一个空间固定的坐标系统的转角。

g_j＝万有引力矢量

|g_j|＝1 $g_j=\left[\begin{array}{c}{g}_{1j}\\ g_{2j}\\ g_{3j}\end{array}\right];e_j={(1-{g_{1j}}^2)}^{-1/2}\left[\begin{array}{c}1-g_{1j}^2\\ -g_{1j}{g}_{2j}\\ -g_{1j}{g}_{3j}\end{array}\right];f_j={(1-{g_{1j}}^2)}^{-1/2}\left[\begin{array}{c}0\\ -g_{3j}\\ g_{2j}\end{array}\right]$

g_1j，g_2j：倾斜传感器测量值，

参数M₁₁…M₃₃，μ₀₁，μ₀₂，μ₀₃，α₁…α_N和I可以应用已知的最佳化计算确定，其中，要求统计和 $L(M,{\mathrm{\&mu;}}_0,{\mathrm{\&alpha;}}_1\&CenterDot;\&CenterDot;\&CenterDot;{\mathrm{\&alpha;}}_{N,},i)=\underset{j=1}{\overset{N}{\mathrm{\&Sigma;}}}{({\mathrm{\&mu;}}_j-{\mathrm{\&mu;}}_o-{\mathrm{Mb}}_{\mathrm{Ej}})}^2$ 具有最小值。

在本发明范围内也可以采用其它的各种求解方法，例如通过卡尔曼滤波器，模糊逻辑方法或神经网络。

不难看出，当ΔM＝0时在数学方法上将显著简化。这种情况是不考虑软磁干扰场的存在。这时M相当于单位矩阵。

上面在说明中参照了倾斜传感器。除此之外还可以应用两个正交安装的编码器。采用这些编码器有可能与一个参照点相应地进行角度度测量。在实际实施中可能应采用一个固定装配的支架，与此支架相比有可能磁针罗盘和干扰系统转动一个角度。

还有可能利用两个在系统上配置的转动率传感器，就是说旋磁确定所期望的倾斜角。通过此转动率信息的积分可以推导出转角。

Claims (8)

1.采用一个包括三个磁场传感器和两个倾斜采集设备的电子磁针罗盘确定可以确定受到固定在仪器上的磁性材料和若干电流的磁场干扰的地磁场方向的方法，其中，电子磁针罗盘置于不同的空间位置N，

在此N位置的每一个位置上测量倾斜采集设备的信号和磁场传感信号并且由这些信号确定倾斜值和磁场值，并且

在此倾斜值和磁场值的基础上计算地磁场矢量的大小，其中采用矢量方程const＝b_g＝b_Esin(i)＝g^Tb_E＝g^Tm(b_oem-b_o)其中，

b_g＝与万有引力矢量g平行的磁场分量

b_E＝磁场矢量b_E的长度

b_gem＝所测量的磁场矢量

b_E＝在测量地点实际的地磁场矢量

ΔM＝软磁场失真，即通过地磁场感应的磁性

M＝I+ΔM；I＝单位矩阵

m＝M^-1

b_o＝硬磁干扰场矢量

i＝磁偏向其中，N至少应与待确定矢量方程的参数数量相等。

2.根据权利要求1所述方法，其特征在于，选择如下的比例确定

b_g＝const(恒定值)，m₁₁＝const(恒定值)；

m₁₁+m₂₂+m₃₃＝const(恒定值)

m₁₁ ²+m₂₂ ²+m₃₃ ²＝const(恒定值)

m₁₁ ²+m₁₂ ²+…+m₃₃ ²＝const(恒定值)

detm＝const(恒定值)

3.根据权利要求1所述方法，其特征在于，如果不考虑软磁干扰场，则ΔM＝0。

4.根据权利要求1所述方法，其特征在于，如果在不同部位的测量数量N多于所提供的方程数量，则磁场矢量借助于最佳化计算确定。

5.根据权利要求1所述方法，其特征在于，地磁场矢量值在确定了倾斜值和磁场值的基础上，在附加的应用方程

const(恒定值)＝b_E ²＝b_E ^Tb_E＝(b_gem-b_o)^Tm^Tm(b_gem-b_o)

b_gem＝所测量的磁场矢量

b_E＝在测量部位实际的地磁场矢量

ΔM＝软磁场失真，就是说通过地磁场感应的磁性

M＝I+ΔM，I＝单位矩阵

m＝M^-1

b_o＝硬磁干扰场矢量借助于统计优化方法加以确定。

6.根据权利要求1至5其中之一所述方法，其特征在于，倾斜采集设备是用的倾斜传感器。

7.根据权利要求1至5其中之一所述方法，其特征在于，倾斜采集设备是用的转动率传感器(旋磁)，其中通过转动率信息的积分推导出倾斜角。

8.根据权利要求1至5其中之一所述方法，其特征在于，倾斜采集设备是用的编码器，用此编码器测量与一个参照点的角度。

Images