基于卡尔曼滤波的磁偶极子目标张量定位
摘 要：针对现有单点磁梯度张量定位方法存在受测量噪声影响容易出现定位结果发散的不足，根据磁偶极子磁感应强度与磁场梯度张量的基本关系式，建立了张量定位的状态空间模型，并采用卡尔曼滤波器对状态向量进行估计，实现了目标定位。设计了数值算例对定位方法进行了检验，定位结果表明基于卡尔曼滤波的定位方法具有较高的精度和稳定性，可用于航空磁探测和水下小目标磁场探测。
关键词：张量定位；状态空间模型；卡尔曼滤波
水中铁磁物体，如钢铁航行体、水雷等，在地磁场的作用下会引起局部地磁异常，磁定位技术就是利用该磁异常信号来确定这些磁性目标的位置等参数[1-3]。基于磁偶极子模型的磁性目标定位常见于航空探潜、水下地磁反演等领域。在航空探潜领域，通常先采用声探测技术来确定可疑海域，再用装载于飞机上的磁探仪接收到的潜艇磁场信号推算目标潜艇的运动要素及其与飞机的相对位置等参数，实现目标精确定位。
对适用于航空探潜的各种磁性目标定位方法，按采用的物理量来可以分为基于磁场总量[4]、
磁场分量[5-7]和单点梯度张量[8,9]等方法，其中单点梯度张量定位方法基于磁偶极子目标的线性探测方程，具有计算量小、实时性好的优点，一直以来是研究的热点[8-10]。文献[8]对单点梯度张量定位方法进行了较为深入的研究，指出目标磁矩特征、测点与目标距离、磁探仪分辨率和磁梯度测量系统的基线距离是影响定位结果的四个主要因素，但却忽略了测量噪声对定位精度的影响分析。在航空探潜的实际工作中，磁场测量噪声不可避免地存在，其主要来源为测量系统噪声和飞机磁场干扰[11]，这些测量噪声对定位结果具有重要影响。为此，本文以航空探潜和水下小目标磁场探测为应用背景，在考虑磁场测量噪声情况下，采用卡尔曼滤波理论对磁偶极子目标梯度张量定位进行研究，并设计了相应的仿真实验对方法的有效性进行验证。
1 目标定位数学模型
在卡尔曼滤波框架下，磁性目标张量定位的数学模型称为状态空间模型，状态空间模型由用于描述系统动态特性的状态模型和描述系统测量信息的量测模型两部分组成。
1.1 状态模型
状态模型的三个要素为状态变量，状态转移矩阵和状态噪声。构建状态模型的基本原则是所建立的状态模型既要符合系统动态特性实际，又要便于数据处理。在航空探潜过程中或小目标磁场探测中，传感器平台的位置、速度和姿态一般由平台的导航系统提供，为已知量，在进行张量定位方法分析时可假定测量平台静止而目标运动，如图1所示，其中坐标系原点o为探测系统的中心，x轴，y轴和z轴分别平行于地磁坐标系相应的三轴。
由于探测过程持续时间较短，一般假定在定位过程中潜艇目标作深度不变的匀速直线运动[10,12-13]。此时，可采用目标在水平方面航速vx，vy以及目标与探测系统的相对位置r来描述其运动过程的动态特性，即可定义状态变量x=(x,滴滴哒哒y,z,vx,vy)T，则：
（1）
目字的成语式中T为测量设备采样周期，k=0,1,2,…N为采样时刻。上式写成矩阵的形式为xk+1=Ak,k-1xk+wk，其中Ak,k-1称为状态转移矩阵，状态噪声wk=(wx,k,wy邓紫棋的歌曲,k,wz,k,wvx,k,wvy,k)T为相互独立的高斯白噪声，其方差为R。
图1 目标张量定位示意图
1.2 量测模型
磁偶极子是磁性目标定位中使用最为广泛且最有效的数学模型。对位于源点S的磁偶极子，其空间中任意场点P处的磁感应强度B(Bx,By挨打的作文,Bz)可按下式计算：
（2）
式中：r为由源点S指向场点P的位置矢量，m为磁偶极子磁矩，μ为介质磁导率。
磁感应强度的三个分量(Bx,By,Bz)在空间三个方向(娃娃鱼养殖x,y,z)的变化率构成一个张量矩阵G=[Gij]3×3，Gij的表达式如下：
（3）
式中：i,j=1,2,3；B1，B2和B3分别表示磁感应强度的三个分量Bx，By和Bz；我的偶像英语作文r1，r2和r3分别表示矢径r的三个分量rx，r总结大会主持词y和rz；δij为Kronecker-Delta函数。
由文献[8,9]，磁偶极子磁场梯度张量与磁场的关系为：
（4）
在磁偶极子目标定位实际工程中，目标磁场一般采用精度较高的超导量子磁力仪进行测量，其测量精度可达到pT级，磁场测量值中包含目标磁场真值、地磁噪声和载体干扰磁场等。磁场梯度张量可采用5个三轴磁传感器[9]或10个单轴磁传感器[14]的输出经差分后得到，该测量方式极大地降低了同模测量噪声的干扰，因此可认为梯度张量测量值即为目标梯度张量真值。为此，观测方程可写为：
（5）
其中Yk和xk分别表示k时刻的观测值（磁场测量值）和状态向量；vk为k鱼皮怎么做好吃又简单时刻的测量噪声，
一般认为其为方差为Q的高斯白噪声；Hk为观测矩阵，具体表达式为：
（6）
此外，在估计得到目标与探测系统的相对位置后，可采用式（7）计算目标磁矩参数，该参数可用于目标的分类识别。
（7）
1.3 可观测性分析
采用卡尔曼滤波求解式（1）和式（5）描述的磁性目标定位模型的前提是定位模型可观测。由线性系统理论，可采用如下的观测性矩阵Qo来分析其可观测性：
（8）
若Qo矩阵的秩大于等于待观测向量维数，则系统可观测。将H，A代入上式后，取Qo阵的某一6阶子阵如下：
， （9）
明显，若||G||≠0，则TG至少存在一个行列式不为零的2阶子阵；从而Qo1至少存在一个行列式不为零的5阶子阵，即若||G||≠0，系统可观测。
