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

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

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为一个总场常数,即飞机远离目标达到一定距离,目标信号趋近于0时的总场输出。由式(7)可知,存在目标信号时,系统输出即为目标信号。1.2基于遗传算法的搜索空间生成技术遗传算法是一种模拟自然界生物进化过程,采用人工进化的方式对目标空间进行随机优化搜索的算法,它的适应度函数不受连续可微的约束,能很好地处理约束条件。遗传算法在搜索过程中能很好地跳出局部最优解,得到全局最优解。遗传算法的主要运算过程包括:1)编码。即把一个问题的可行解从其解空间转换到遗传算法所能处理的搜索空间。常用的编码方法有二进制编码与浮点数编码。2)生成初始群体Chrom(0)。设置群体大小为M,采用随机初始化的方式产生包含M个个体的初始群体。3)适应度值评价检测。适应度值函数表明个体或解的优劣性。对于不同的问题,适应度函数的定义方式不同。根据具体问题,计算群体中各个个体的适应度。4)选择。根据每个个体适应度值大小选择,使得适应度较高的个体被遗传到下一代群体中的概率较大,适应度低的个体被遗传到下一代群体中的概率较小。5)交叉。按一定的概率从群体中选择两个个体,交换两个个体的某个或某些位,从而形成两个新个体。交叉运算决定了遗传算法的全局搜索能力。6)变异。将个体编码串中的某些基因座上的基因值用该基因座的其他等位基因来替换,从而形成一个新的个体。变异运算决定了遗传算法的局部搜索能力。重复3)～6)至设定的最大遗传代数时停止,输出最优个体的解。2计算机制码采用遗传算法求解飞机磁场模型补偿系数的主要步骤为:1)编码。将12个待求解系数按浮点数编码方式进行编码,每个参数编码的取值范围如表1所示。这里采用浮点数编码方式是因为该编码方式适用于求解精度要求较高的场合,同时它也能提高运算效率。2)初始化种群。经过编码后,会产生一个初始种群,这里选择初始种群的大小为200。较大数目的初始种群可在搜索空间中的更多解搜索,容易找到全局最优解。3)根据适应度函数计算每个个体的适应度。按式(3)计算每个个体对所测得的每组外界磁场数据还原后的总场值,还原后的所有总场值标准差越小则适应度值越大。4)选择。以轮盘赌的方式进行正比例选择,这里选择概率取0.8。适应度值越高的个体,被选中的可能性就越大,进入下一代的概率就越大。5)交叉。按两两配对原则将群体中个体进行配对并进行交叉运算,这里交叉概率选为0.4,采用普通单点交叉方式进行。6)变异。根据设定的概率进行变异操作,得到新一代个体,这里选择变异概率为0.05。重复3)～6)至设定的最大遗传代数时停止,最后输出的最优个体解即为磁干扰模型补偿系数。3模拟测试3.1仿真结果与分析为获取相关磁场数据,模拟飞机从地磁正北方向开始按顺时针方向飞行一周,测量飞机运动时光泵的总场输出和三轴磁力仪三轴输出数据。假设所选取试验地区地磁场总场大小为51200nT,3个硬磁干扰参数为[-694.02,831.32,229.87],9个软磁干扰参数分别为a11=1.0379,a12=0.0099,a13=0.0216,a21=-0.0062,a22=1.0403,a23=-0.0079,a31=-0.0082,a32=0.0301,a33=1.0313,考虑当前三轴磁通门磁力仪静态噪声可达到0.1nT的水平,光泵磁力仪静态噪声水平可达到10pT左右,仿真中给三轴磁力仪加入0.1nT白噪声,给光泵加入10pT白噪声,利用这些假设条件,仿真生成的飞机飞行一周时光泵输出的磁场总场数据如图2所示,三轴磁力仪输出的三轴磁场数据如图3所示。利用三轴磁力仪测得的数据计算3个方向余弦cosα,cosβ,cosγ,对光泵总场数据进行三轴方向上的分解,得到包含外界硬磁、软磁干扰的三轴分量Hx,Hy,Hz,将此组数据代入遗传算法程序中,以式(2)、式(3)为适应度函数,求解出最优补偿参数,利用补偿参数对磁干扰进行补偿后的系统输出如图4所示。从仿真结果可知,利用光泵与磁通门数据代入遗传算法可稳定求出有效的系统补偿参数,补偿后系统输出峰峰值优于20pT,效果较为理想。3.2飞行坐标系仿真现在探测区域内加入一个模拟舰艇磁性目标,飞机模拟直线探测飞行,以验证系统还原微弱舰艇磁性信号的能力。首先建立目标舰艇磁矩坐标系和飞机飞行坐标系,如图5所示。假定这两个坐标系原点重合,目标舰艇磁矩坐标系中x1轴与舰艇纵轴方向平行,指向船首为正,y1轴与舰艇横轴平行,指向船体右侧为正,z1轴与平面x1Oy1垂直,向下方向为正,构成右手坐标系;飞机飞行坐标系中x2轴与飞机直线探测飞行方向重合,前进的方向为正,y2轴与x2轴垂直指向左侧机翼为正,z2轴与平面x2Oy2垂直,向上方向为正,也构成一个右手坐标系。目标舰船相对飞机为低速运动,可认为其在短时间内相对静止,假定其磁航向角为135°。根据舰艇磁偶极子模型,假定目标舰艇三轴磁矩分别为Mx=150000A·m2,My=3000A·m2,Mz=8000A·m2。在飞行坐标系中,飞机的初始坐标假定为(-3000,500,260),即飞机初始纵向距离为-3000m,正横距离为500m,距离目标舰艇高度为260m。假定飞机沿磁航向角为315°方向作直线探测飞行。该探测状态下,真实磁偶极子目标信号如图6所示,未进行干扰补偿时的系统输出如图7所示,采用基于遗传算法的磁干扰补偿方法进行补偿后,系统输出信号如图8所示。从这3个图中,我们可以很明显地看出,未经磁干扰补偿,微弱