磁干扰模型补偿系数的复共线性研究.docxVIP

磁干扰模型补偿系数的复共线性研究 0 干扰磁场模型求解方法国外发展现状 航空磁体勘探广泛应用于考古、矿物勘探和水磁体勘探领域。地磁场和飞机磁体的干扰将重叠在磁体目标信号的测量结果中。为了准确探测弱磁体的磁体信号，有必要补偿飞机自身和运动引起的干扰变量。这就要对飞机磁干扰场的组成和特性进行研究,建立准确的数学模型,用合适的方法求解出模型系数,从而对磁干扰进行补偿,提高磁探质量。 20世纪40年代,美国学者最早开始对飞机干扰磁场的数学模型进行研究,参考文献中给出了飞机干扰磁场的数学建模。由于模型具有较强的复共线性,他们未对模型的求解和干扰磁场的实时补偿提出有效的方法。在其后几十年时间,国内外一些专家学者相继提出各种模型系数求解方法。参考文献提出了小幅度机动磁干扰补偿方法,但这种方法是建立在飞机做小幅度机动的条件下,具有一定的局限性且求解的精度不高。参考文献研究了截断奇异值分解的飞机干扰磁场模型求解算法,此方法对奇异值截断位置比较敏感,如截断位置选择不合理,会使估计误差迅速增大。针对上述问题,本文提出了基于遗传算法的磁干扰系数求解方法。 1 飞机磁体模型的解决方案和遗产算法 1.1 补偿系数求解算法 飞机干扰磁场包括硬磁恒定干扰场、软磁感应干扰场和涡流干扰场,本文采用文献中的模型来描述干扰磁场。建立坐标系如图1所示,x轴与飞机纵轴平行,机首方向为正,y轴与飞机横轴平行,向右方向为正,z轴与平面xOy正交,向下方向为正。α,β,γ为H和3个坐标轴之间的夹角,其方向余弦为cosα,cosβ,cosγ。 在实际测量中,可用高精度光泵磁强计测量含干扰外部磁场的总场H,三轴磁强计测量含干扰外部磁场的3个方向余弦cosα,cosβ,cosγ,从而可得三轴方向上含干扰磁场的输出Hx,Hy,Hz。 作为一种算法研究,本文暂不考虑涡流场干扰,则由参考文献可知磁场模型为: ?????Hx=x0+a11Tx+a12Ty+a13TzHy=y0+a21Tx+a22Ty+a23TzHz=z0+a31Tx+a32Ty+a33Tz(1){Ηx=x0+a11Τx+a12Τy+a13ΤzΗy=y0+a21Τx+a22Τy+a23ΤzΗz=z0+a31Τx+a32Τy+a33Τz(1) 为求解补偿系数,可选取一处地磁场较为稳定的区域让飞机飞行,以测量三轴上的磁场输出Hx,Hy,Hz,由式(1)可得: ???TxTyTz???=???a11a21a31a12a22a32a13a23a33????1???Hx?x0Hy?y0Hz?z0???(2)[ΤxΤyΤz]=[a11a12a13a21a22a23a31a32a33]-1[Ηx-x0Ηy-y0Ηz-z0](2) 则总场输出: O1=Tx2+Ty2+Tz2?????????????√(3)Ο1=Τx2+Τy2+Τz2(3) 此时总场输出O1为一个固定的当地地磁总场值,此式可作为遗传算法适应度函数的输入,用于求解12个补偿系数。 假设存在目标信号,目标信号在三轴上的分量分别为Sx,Sy,Sz,由式(1)得: ?????Hx=x0+a11(Tx+Sx)+a12(Ty+Sy)+a13(Tz+Sz)Hy=y0+a21(Tx+Sx)+a22(Ty+Sy)+a23(Tz+Sz)Hz=z0+a31(Tx+Sx)+a32(Ty+Sy)+a33(Tz+Sz)(4)???Tx+SxTy+SyTz+Sz???=???a11a21a31a12a22a32a13a23a33????1???Hx?x0Hy?y0Hz?z0???(5){Ηx=x0+a11(Τx+Sx)+a12(Τy+Sy)+a13(Τz+Sz)Ηy=y0+a21(Τx+Sx)+a22(Τy+Sy)+a23(Τz+Sz)Ηz=z0+a31(Τx+Sx)+a32(Τy+Sy)+a33(Τz+Sz)(4)[Τx+SxΤy+SyΤz+Sz]=[a11a12a13a21a22a23a31a32a33]-1[Ηx-x0Ηy-y0Ηz-z0](5) 此时总场输出: O2=(Tx+Sx)2+(Ty+Sy)2+(Tz+Sz)2??????????????????????????????√(6)Ο2=(Τx+Sx)2+(Τy+Sy)2+(Τz+Sz)2(6) 由式(6)可知,利用求解出的补偿系数,可得到总场输出O2,其包含地磁场与目标信号,而我们只关心目标信号的大小,故可采用式(7)作为系统输出O。 O=O2?M(7)Ο=Ο2-Μ(7) 式(7)中, M=limSx,Sy,Sz→0(Tx+Sx)2+(Ty+Sy)2+(Tz+Sz)2??????????????????????????????√Μ=limSx,Sy,Sz→0(Τx+Sx)2+(Τy+Sy)2+(Τz+Sz)2 为一个总场常数,即飞机远离目