一种飞机航磁干扰的补偿方法
【技术领域】
[0001] 本发明涉及属于航磁探测领域，尤其涉及一种飞机航磁干扰的补偿方法。
【背景技术】
[0002] 航磁干扰补偿技术是在航磁探测过程中去除航空运动平台自身磁干扰的一种技 术。通过分析航空运动平台自身磁干扰的类型和性质，建立航空运动平台磁干扰数学模型， 然后在校准飞行过程中按照规定的方法测得磁总场及三分量数据，并将其用来计算航空运 动平台磁干扰数学模型的系数。在实际进行航磁探测时，利用估计出的模型系数及飞机姿 态数据估计航空运动平台产生的磁干扰并将其从磁总场中去除，从而得到不含航空运动平 台磁干扰的磁场数据。现有航磁干扰补偿系数计算方法多是基于T-L模型，该模型将航空 运动平台磁干扰分为恒定场、感应场和涡流场三种类型，其中恒定场系数有3项，感应场系 数和涡流场系数各有9项，且两者分别与地磁场的大小及变化率有关。由于可以把补偿系 数的估计过程看作线性回归，因此系数求解矩阵的性质势必会影响估计结果的精度。其中， 系数求解矩阵的病态性是影响补偿系数估计精度的关键因素。
【发明内容】
[0003] 为了降低系数求解矩阵的病态性，本发明提供了一种降低系数求解矩阵病态性的 预处理方法，该方法可应用于补偿系数估计算法并有助于提高其求解精度。
[0004] 本方法解决其技术问题所采用的技术原理是：
[0005] 在校准飞行过程中总场磁力仪测得的磁场数据为NX 1列向量Ητ，实际地磁场值为 NX 1列向量HE，而飞机产生的磁干扰为NX 1列向量托，则有：
[0006] HT= He+H: (1)
[0007] 其中，只有Ητ通过直接测量得到；
[0008] 航磁干扰补偿的最终目的是确定HE，观察上式可以发现，只要计算出氏并将其从 Ητ中减掉就可以得到He。
[0009] 根据T-L模型，飞机产生的磁干扰可以表示为：
[0010]
(2)
[0011] 其中，Ppaybu为需要估计出的航磁干扰补偿系数，ujpu,为由地磁场与飞机轴 向所成夹角的方向余弦，〃为^的导数；
[0012] 三分量磁力仪测得的磁场三分量为Vl、v2、v3,Ul可根据下式计算得到：
[0013]
⑶
[0014] 根据几何对称性，公式（2)可以转化成更为紧凑的形式：
[0015]
[0016]式中Θ是由待定系数构成的列向量，δ是由Ui、ua和〃〃构成的行向量。
[0017]根据公式（4)，氏可表示为：
[0018] Η：= Δ Θ (5)
[0019] 其中，A= 茗f为NX16矩阵；飞机产生的磁干扰与校准机动动作的频 率有关，而地磁场与校准机动动作的频率无关，设bpf为FIR带通滤波器，结合式（1)和（2) 则有：
[0020] bpf(HT) =bpf(HE)+bpf(H:) =bpf(HE)+bpf(A)Θ (6)
[0021] 其中bpf()表示对矩阵的每一列进行带通滤波；
[0022] 由于bpf的通带范围是根据飞机校准机动动作的频率设置的，故bpf(HE)~0,则 公式（6)转化为：
[0023] bpf (H:) = bpf ( A ) θ (7)
[0024] 补偿系数构成的列向量可以根据式（7)计算得到。
[0025] 容易看出bpf (Δ)的性质将直接影响对Θ的估计。由于Δ是确知的，因此通过 调整bpf的相关参数可以影响bpf(A)的性质。
[0026] 本发明的主要目标是降低bpf(A)的病态性程度，利用小波分解技术使得 bpf(A)的病态性最低；
[0027] 本发明为解决现有航磁干扰补偿系数计算方法中由于系数求解矩阵的病态性而 影响补偿系数估计精度的问题，而提出一种飞机航磁干扰的补偿方法。
[0028] -种飞机航磁干扰的补偿方法，按以下步骤进行：
[0029]一、构造矩阵Δ :
[0030] 在校准飞行过程中总场磁力仪测得的磁场数据为NX1列向量Ητ，实际地磁场值为 NX1列向量HE，而飞机产生的磁干扰为NX1列向量托，则有：
[0031] HT= He+H: (1)
[0032] 其中，只有Ητ通过直接测量得到；
[0033] 航磁干扰补偿的最终目的是确定HE，观察上式可以发现，只要计算出氏并将其从 Ητ中减掉就可以得到Η E;
[0034] 根据T-L模型，飞机产生的磁干扰可以表示为：
[0035]
(2)
[0036]其中，为需要估计出的航磁干扰补偿系数，uJPu,为由地磁场与飞机轴 向所成夹角的方向余弦，4为^的导数；
[0037] 三分量磁力仪测得的磁场三分量为Vl、v2、v3, Ul可根据下式计算得到：
[0038]
⑶
[0039] 根据几何对称性，公式（2)可以转化成更为紧凑的形式：
[0040]
[0041] 式中Θ是由待定系数构成的列向量，δ是由Ui、ua和构成的行向量；
[0042] 根据公式（4)，氏可表示为：
[0043] Η：=ΔΘ (5)
[0044] 其中，Δ= 为NX 16矩阵；
[0045] 二、利用小波分解矩阵Δ:
[0046] 矩阵Λ的每一列为fi;(i= 1，2,…16)，则利用小波分解可将仁分解成如下形式：
[0047]fi=f〇+w1+w2+···+Wj (8)
[0048] 上式中％，（1 = 1，2, ···，j)表示t中的不同频率成分，且1越大w:的频率越高， Wj频率最高高表示fi中的噪声；f。表示fi中的极低频率成分，即地磁场分量；
[0049] 三、对步骤二的结果实现带通滤波：
[0050] 定义U1= {Wll|l彡1彡j-Ι}为对矩阵Δ中的第i列进行小波分解的结果，其 中fi。和Wi_j均已忽略；对于任意的ui需要确定相同的连续下标集合L= {s,S+1,…t} (1 <s<t<j-1)，则针对匕的相应带通滤波结果为：
[0051] bpf^f,) =wls+wl(s+1)+---wlt (9)
[0052] 进而有：
[0053] bpfL(A) = (bpf^fi),bpfL(f2), ···,bpfL(f16)) (10)
[0054] 四、通过确定下标集合L使得bpfV(A)的病态性最弱：
[0055] 利用条件数来度量矩阵bpf； (Δ)的病态性程度：
[0056]
(11)
[0057] 其中，λ_和λ_分别表示bpf\( △)的最大特征值和最小特征值， C〇nd(bpf；(A))取值越大，矩阵的病态性越严重；
[0058] 由于需要选取的下标是连续的，可知L= {s,s+1,…t} (1彡s彡t彡j-Ι)的取 值共有（j_l)j/2种不同情况，j的取值不超过10,本发明中利用穷举法来确定，使得
5小。
[0059] 五、确定补偿系数方程并求解航磁干扰补偿系数：
[0060] 根据步骤四得到补偿系数方程如下：
[0061]
02)
[0062] 由上式求解航磁干扰补偿系数。
[0063] 本发明包括以下有益效果：
[0064] 1、提出的方法有效降低了系数求解矩阵的病态性，进而提高了补偿系数估计算法 的求解精度；
[0065] 2、本发明方法，为了实现带通滤波，合理地忽略了Wj分量和f。分量，既不影响计算 结果，又简化了计算过程；
[0066] 3、本发明中利用穷举法来确定，使得(Δ))最小，最大程度地降低了 系数求解矩阵的病态性。
【