地磁背景下基于传感器阵列的磁偶极子目标跟踪方法

陈路昭; 冯永强; 郭瑞杰; 朱万华; 方广有

doi:10.11999/JEIT190236

地磁背景下基于传感器阵列的磁偶极子目标跟踪方法

doi: 10.11999/JEIT190236

陈路昭^{1, 2, 3},
冯永强^{1, 2},
郭瑞杰³,
朱万华^1, ,,
方广有¹

1.
中国科学院电磁辐射与探测技术重点实验室北京 100190
2.
中国科学院大学北京 100049
3.
北京自动化控制设备研究所北京 100074

基金项目: 国家自然科学基金青年基金(41704177)和国家重点研发计划“深地资源勘查开采”重点专项(2018YFC0603201)

详细信息

作者简介:
陈路昭：男，1992年生，博士，研究方向为传感器阵列的磁信号处理方法、运动平台磁干扰补偿技术

冯永强：男，1985年生，博士生，研究方向为航磁系统硬件设计与航磁数据处理方法

郭瑞杰：女，1992年生，硕士，研究方向为组合导航与系统控制技术

朱万华：男，1982年生，副研究员，主要研究方向为高灵敏度磁场传感器关键技术、高性能磁屏蔽关键技术

方广有：男，1963年生，研究员，主要研究方向为超宽带电磁场理论及工程应用、超宽带雷达成像技术、微波成像新方法和新技术

通讯作者:
朱万华　whzhu@mail.ie.ac.cn

中图分类号: O411.5; TN911.73
计量
- 文章访问数: 3705
- HTML全文浏览量: 1228
- PDF下载量: 134
- 被引次数: 14
出版历程
- 收稿日期: 2019-04-10
- 修回日期: 2019-12-01
- 网络出版日期: 2019-12-09
- 刊出日期: 2020-03-19

Magnetic Dipole Object Tracking Algorithm Based on Magnetometer Array in Geomagnetic Background

Luzhao CHEN^{1, 2, 3},
Yongqiang FENG^{1, 2},
Ruijie GUO³,
Wanhua ZHU^{1
, ,},
Guangyou FANG¹

1.
Key Laboratory of Electromagnetic Radiation and Sensing Technology, Chinese Academy of Sciences, Beijing 100190, China
2.
University of Chinese Academy of Sciences, Beijing 100049, China
3.
Beijing Automation Control Equipment Institute, Beijing 100074, China

Funds: The Youth Program of National Natural Science Foundation of China (41704177), The National Key R&D Program of China (2018YFC0603201)

摘要

摘要:
针对地磁背景下磁偶极子目标跟踪过程中存在的地磁干扰与模型非线性的问题，该文提出一种基于差量磁异常的蒙特卡洛卡尔曼滤波(MCKF)跟踪方法。新的跟踪方法以传感器阵列测量磁场的差量作为观测信号，并利用蒙特卡洛卡尔曼滤波算法解决模型的非线性问题，实现磁偶极子目标的实时跟踪。通过仿真跟踪实验，结果表明该文算法较传统的扩展或无迹卡尔曼滤波算法在稳定跟踪过程中对目标特征参数的估计更精确；通过地磁背景跟踪实验，结果验证了该文算法较传统算法在低信噪比下的性能优势。
- 磁偶极子 /
- 蒙特卡洛卡尔曼滤波 /
- 传感器阵列 /
- 地磁背景
Abstract:
In order to solve the problem of geomagnetic interference and model nonlinearity in the tracking process of magnetic dipole under geomagnetic background, Monte Carlo Kalman Filter (MCKF) tracking method based on differential magnetic anomaly is proposed in this paper. The new tracking method takes the difference of magnetic field measured by sensor array as the observation signal, and uses Monte Carlo Kalman Filtering (MCKF) algorithm to solve the nonlinear problem of the model to realize the real-time tracking of magnetic dipole targets. The simulation results show that the proposed method is more accurate than the traditional Extended Kalman Filter (EKF) or Untracked Kalman Filter (UKF) in the stable tracking process. The results of real geomagnetic background tracking experiments show that the proposed algorithm has better tracking performance under low SNR.
- Magnetic dipole /
- Monte Carlo Kalman Filter (MCKF) /
- Sensor array /
- Geomagnetic background

HTML全文

1. 引言

随着浅层陆地资源的日益匮乏与开发难度的增大，世界各国对海洋资源的开发与海域权益的争夺日趋激烈。面对部分邻国蓄意挑起的海洋争端及来自境外军事侦查、非法入侵和资源掠夺的敌方潜艇及无人潜航器等对我国的临海主权，尤其是对港口、码头及军事设施等的安全构成严重威胁^[1,2]。磁探测技术是一种被动探测技术，因其能够有效弥补声学近海、浅海探测的不足受到国内外的广泛关注^[3]。

磁性目标的特性识别主要包括优化定位与实时滤波跟踪两大类。其中，优化定位算法通常依赖大量观测数据，比如传感器网络、磁梯度系统或航空测线数据等，通过最优化的方法获取目标位置、磁矩等参量。优化定位算法研究的典型代表有，以色列的Sheinker等人^[4,5]在2005～2008年提出依次基于遗传算法、模拟退火算法及种群增量学习算法实现磁偶极子目标的定位。2010年Yang等人^[6]将粒子群与(Levenberg-Marquardt, LM)优化算法结合实现了3个磁偶极子目标的定位与跟踪，2018年Gao等人^[7]将LM算法用于运动磁性目标的定位。此外，利用磁梯度张量系统对UXO目标的定位与识别UXO同样吸引了大量研究人员的关注^[8-10]，其通常将搭载在移动平台上，如无人机等。在目标的实时滤波跟踪方面，通常利用静止平台的磁传感器阵列测量目标磁异常，实现对目标体的实时跟踪。滤波跟踪算法研究的典型代表有，Birsan等人^[11,12]在2003～2006年利用扩展/无迹卡尔曼滤波方法和无迹粒子滤波算法等迭代贝叶斯算法对磁偶极子目标进行实时跟踪；2008年，Kozick等人^[13]提出磁传感器阵列目标跟踪的矢量与标量模型，并进行仿真实验研究；2012年日本学者Alhmiedat等人^[14]提出将无线磁传感器网络用于实现军事目标的识别与跟踪应用。国内研究学者在近几年针对磁性目标的跟踪已开展大量工作，主要的研究成果有：2010年，张朝阳等人^[15]对磁性目标跟踪中磁偶极子模型适用性进行研究；2012年，于振涛等人^[16]利用海底磁力仪阵列对水中磁目标定位与跟踪方法进行研究；2013年，吴志东等人^[17,18]提出用于实现磁偶极子目标的非线性跟踪算法；2014年，张宏欣等人^[19]提出改进的粗糙化粒子滤波的算法改善磁偶极子目标跟踪中的粒子贫化；2015年，周穗华等人^[20]提出混合卡尔曼滤波算法解决经典粒子滤波在磁偶极子跟踪中高维模型估计时维数灾问题；2017年吴垣甫等人^[21]提出一种递推更新的卡尔曼滤波方法对磁偶极子目标进行跟踪，有效抑制磁偶极子跟踪中初始误差过大导致的滤波发散问题；张宏欣等人^[22]提出渐近贝叶斯估计算法和多模型自适应滤波算法对磁偶极子目标进行跟踪，有效抑制目标跟踪中初始误差导致的性能下降和发散等问题；在2018年张宏欣^[23]提出多模型自适应滤波算法对磁偶极子目标进行跟踪，算法可在目标源及位置等先验信息完全未知条件下准确估计出真实目标状态，提升真实目标跟踪中的适应力。

分析当前国内外研究现状，磁目标跟踪的主要研究对象在于非线性跟踪模型的优化与不同体制跟踪系统的研究，其中很少考虑到复杂跟踪环境对目标跟踪性能的影响，比如海水中目标在运动过程中存在的涡流磁场、尾流磁场^[18]，以及目标跟踪过程中存在地磁背景干扰等，都会对磁目标的跟踪性能产生影响。本文主要针对磁目标跟踪过程中存在的地磁背景干扰的问题开展相关的研究，并在传统磁跟踪模型的基础上提出利用差量磁场对磁偶极子进行跟踪的模型，并考虑到差量磁场模型的高度非线性度，提出基于蒙特卡洛卡尔曼滤波(Monte Carlo Kalman Filter, MCKF)跟踪方法提升跟踪性能。

2. 磁偶极子运动模型

铁磁性目标周围产生磁场叠加原背景磁场导致周围的磁场发生畸变，称为磁异常。当目标探测距离较近时，磁性目标可等效为若干个偶极子源叠加；当探测距离为目标尺寸3倍以上时，可将目标视为磁偶极子，其矢量磁场如式(1)，式中 ${\mu _0}$ 为自由空间磁导率， ${{m}}$ 为目标磁矩矢量， ${{r}}$ 为距离矢量，其幅度为 $r$ 。

$B\left( {{{m}},{{r}}} \right) = \frac{{{\mu _0}}}{{4\pi }}\frac{{3\left( {{{r}} \cdot {{m}}} \right){{r}} - {r^2}{{m}}}}{{{r^5}}}$

(1)

确立直角坐标系，式(1)可以写为如式(2)形式。式中 $\left( {{x_0},{y_0},{z_0}} \right)$ 代表了在3维坐标系中矢量磁力仪的坐标， $\left( {x,y,z} \right)$ 代表坐标系中磁偶极子的坐标； $\left( {{B_x},{B_y},{B_z}} \right)$ 代表了磁场矢量在各坐标轴上的投影分量，而 $\left( {{m_x},{m_y},{m_z}} \right)$ 代表了磁矩矢量在3维坐标系统的投影分量。

$\left[ {\begin{array}{*{20}{c}} {{B_x}} \\ {{B_y}} \\ {{B_z}} \end{array}} \right] = \frac{{{\mu _0}}}{{4\pi {r^5}}} \cdot \left[ {\begin{array}{*{20}{c}} {3{{\left( {x - {x_0}} \right)}^2} - {r^2}}&{3\left( {x - {x_0}} \right)\left( {y - {y_0}} \right)}&{3\left( {x - {x_0}} \right)\left( {z - {z_0}} \right)} \\ {3\left( {y - {y_0}} \right)\left( {x - {x_0}} \right)}&{3{{\left( {y - {y_0}} \right)}^2} - {r^2}}&{3\left( {y - {y_0}} \right)\left( {z - {z_0}} \right)} \\ {3\left( {z - {z_0}} \right)\left( {x - {x_0}} \right)}&{3\left( {z - {z_0}} \right)\left( {y - {y_0}} \right)}&{3{{\left( {z - {z_0}} \right)}^2} - {r^2}} \end{array}} \right] \cdot \left[ {\begin{array}{*{20}{c}} {{m_x}} \\ {{m_y}} \\ {{m_z}} \end{array}} \right]$

(2)

由式(2)可见，任意磁偶极子目标产生的磁场仅依赖于目标与传感器的相对位置及磁偶极子的磁矩。因此，通过对目标位置与磁矩参量的估计，即可完成对运动目标的实时跟踪。

磁偶极子目标跟踪过程中状态变量由两部分组成：在t时刻，描述目标状态的位置参数 $\left[ {x\left( t \right),y\left( t \right),z\left( t \right)} \right]$ 和速度参数 $\left[ {{v_x}\left( t \right),{v_y}\left( t \right),{v_z}\left( t \right)} \right]$ 以及目标的磁矩参数 $\left[ {{m_x}\left( t \right),{m_y}\left( t \right),{m_z}\left( t \right)} \right]$ 。因此，目标状态向量 ${{{x}}_{{t}}}$ 如式(3)。

$\begin{split} {{{x}}_{{t}}} =\,& \left[ x\left( t \right),y\left( t \right),z\left( t \right),{v_x}\left( t \right),{v_y}\left( t \right),{v_z}\left( t \right),\right. \\ & \left.{m_x}\left( t \right),{m_y}\left( t \right),{m_z}\left( t \right) \right] \end{split}$

(3)

考虑到水面舰船或航行器等磁性目标运动状态相对稳定，目标运动状态可以用恒速模型和恒磁矩模型描述，因此目标运动状态方程如式(4)。

$\left. \begin{array}{l} x\left( t \right) = x\left( {t - \Delta t} \right) + \Delta t \cdot {v_x}\left( {t - \Delta t} \right) \\ y\left( t \right) = y\left( {t - \Delta t} \right) + \Delta t \cdot {v_y}\left( {t - \Delta t} \right) \\ z\left( t \right) = z\left( {t - \Delta t} \right) + \Delta t \cdot {v_z}\left( {t - \Delta t} \right) \\ {v_x}\left( t \right) = {v_x}\left( {t - \Delta t} \right) \\ {v_y}\left( t \right) = {v_y}\left( {t - \Delta t} \right) \\ {v_z}\left( t \right) = {v_z}\left( {t - \Delta t} \right) \\ {m_x}\left( t \right) = {m_x}\left( {t - \Delta t} \right) \\ {m_y}\left( t \right) = {m_y}\left( {t - \Delta t} \right) \\ {m_z}\left( t \right) = {m_z}\left( {t - \Delta t} \right) \\ \end{array} \right\}$

(4)

根据目标的位置参量与磁矩参量，通过式(2)，可以获得不同位置传感器的磁异常，当使用3个矢量传感器阵列时，各传感器测量场表示为 $\left[ b_x^i\left( t \right),\right.$ $\left.b_y^i\left( t \right),b_z^i\left( t \right),i = 1,2,3 \right]$ 。因此，t时刻的观测场 ${{{z}}_{{t}}}$ 可表示如式(5)。

$\begin{split} {{{z}}_{{t}}} =\,& \left[ b_x^1\left( t \right),b_y^1\left( t \right),b_z^1\left( t \right),b_x^2\left( t \right),b_y^2\left( t \right),b_z^2\left( t \right),\right.\\ & \left.b_x^3\left( t \right),b_y^3\left( t \right),b_z^3\left( t \right) \right] \end{split}$

(5)

综上，偶极子目标跟踪过程中的状态方程式(4)为线性的；而量测方程式(2)、式(5)为非线性方程。

3. 基于差量磁场的MCKF跟踪算法

3.1 差量磁场信号

传感器测量磁场中同时包含目标磁异常与背景地磁场。通常，采用远参考传感器来消除背景地磁场的影响，但基于参考传感器的方法存在两个明显缺陷：其一，参考传感器与传感器阵列的距离很难把握，距离过近则容易受目标磁异常影响；距离过远，则由于相关性的下降致使背景场消除不彻底。其二，参考传感器设置处需要可靠的背景环境，要保证工作中不受外界磁干扰影响。本文针对上述缺陷，提出传感器阵列测量磁场差量实现磁偶极子目标跟踪，在消除背景地磁场干扰的同时，又避免参考传感器的设置难题。

通过仿真示例，阐述磁偶极子差量磁场与测量场及梯度磁场的特性。仿真磁偶极子运动轨迹如图1所示，其轨迹中类似波浪式的运动用于模拟真实目标运动过程中的随机性，“*”代表测量磁传感器，仿真磁偶极子磁矩为 $\left[ {10,\;20,\;30} \right]{\rm{A}} \cdot {{\rm{m}}^2}$ 。

图 1 仿真磁偶极子运动轨迹与传感器分布位置

下载: 全尺寸图片幻灯片

传感器测量磁场如图2所示，其中图2(a)代表传感器1测量的磁场3分量；图2(b)代表传感器1处的磁梯度场中3个典型量，依次是x分量场在x方向梯度、y分量场在y方向梯度及z分量场在z方向的梯度；图2(c)代表传感器1, 2的差量场。

图 2 传感器测量结果

下载: 全尺寸图片幻灯片

从图2可知，差量场与传感器原始测量磁场在磁异常幅度值上相当，而梯度磁场最小。因此，与磁梯度场跟踪相比，采用差量磁场进行目标跟踪时可以获得更远的探测距离。

3.2 基于MCKF的跟踪算法

常用的非线性滤波跟踪算法有扩展卡尔曼滤波(Extended Kalman Filter, EKF)和无迹卡尔曼滤波(Unscented Kalman Filter, UKF)等。EKF算法利用泰勒级数展开，对非线性的系统方程作1阶线性化近似处理，仅适用于弱非线性系统。UKF以确定性的方式选择一组近似高斯分布的样本点，通过非线性方程变换，获得均值和方差的估计，是最常用的非线性卡尔曼滤波算法。本文提出的蒙特卡洛卡尔曼滤波方法，与UKF原理相似，区别在于通过蒙特卡洛的方法生成随机sigma样本点，通过非线性方程变换，获取非线性系统均值与方差估计，完成卡尔曼滤波增益计算与状态的更新，算法流程如图3所示。

图 3 基于MCKF跟踪算法的流程图

下载: 全尺寸图片幻灯片

基于图3的算法流程，目标状态变量 ${{x}}$ 及其状态方程与观测方程分别用 $f\left( {{x}}\right)$ 和 $h\left( {{x}} \right)$ 表示，则磁偶极子的蒙特卡洛卡尔曼滤波跟踪算法的步骤可归结如下：

(1) 滤波器初始化：定义初始状态变量为 $9 \times 1$ 维列向量，记为x₀；初始化协方差矩阵为 $9 \times 9$ 矩阵，标记为P₀^xx。

(2) 状态变量预测：根据蒙特卡洛规则随机生成N个样本点，记 ${{{x}}^{(i)}}_{n - 1|n - 1} = {{{x}}_{n - 1|n - 1}} +$ ${{{D}}_{n - 1|n - 1}}c(i),i =1,2, ··· ,N$ ，样本点权重记为 ${w_i} = 1/N,i = 1,2, ··· ,N$ 。根据随机样本点获取状态量均值预测记为式(6)，其协方差矩阵预测如式(7)，式中 $c\left( i \right)$ ～ $N\left( {0,{{I}}} \right),i = 1,2, ··· ,N$ ; ${{{D}}_{n - 1|n - 1}} =$ ${\left[ {{{P}}_{n - 1|n - 1}^{{xx}}} \right]^{ - 1}},n = 2,3,4, ···$ ; Q为变量预测协方差矩阵。

${{{x}}_{n|n - 1}} = \mathop \sum \limits_{i = 1}^N {w_i}f\left( {{{x}}_{n - 1|n - 1}^{\left( i \right)}} \right) \hspace{58pt}$

(6)

$\begin{split} {{P}}_{n|n - 1}^{{xx}} =\,& \mathop \sum \limits_{i = 1}^N {w_i}\left[ {f\left( {{{x}}_{n - 1|n - 1}^{\left( i \right)}} \right) - {{{x}}_{n|n - 1}}} \right]\\ &\cdot {\left[ {f\left( {{{X}}_{n - 1|n - 1}^{\left( i \right)}} \right) - {{{x}}_{n|n - 1}}} \right]^{\rm{T}}} + {{Q}} \end{split}$

(7)

(3) 观测值及其协方差矩阵估计：根据蒙特卡洛规则生成样本点 ${{x}}_{n|n - 1}^{\left( i \right)} = {{{x}}_{n|n - 1}} + {{{D}}_{n|n - 1}}c\left( i \right),$ $i = 1,2, ··· ,N$ ，通过量测方程估计观测均值如式(8)。利用输出样本点及输入样本点获取观测值协方差及观测量与状态变量互协方差矩阵如式(9)，式中 ${{{D}}_{n|n - 1}} = {\left[ {{{P}}_{n|n - 1}^{{xx}}} \right]^{ - 1}},n = 2,3,4, ···$ ，R为观测噪声协方差矩阵。

${{{z}}_{n|n - 1}} = \mathop \sum \limits_{i = 1}^N {w_i}h\left( {{{X}}_{n|n - 1}^{\left( i \right)}} \right)\hspace{73pt}$

(8)

$\left. \begin{split} {{P}}_{n|n - 1}^{{zz}} =\,& \mathop \sum \limits_{i = 1}^N {w_i}\left[ {h\left( {{{X}}_{n - 1|n - 1}^{\left( i \right)}} \right) - {{{z}}_{n|n - 1}}} \right] \\ &\cdot {\left[ {h\left( {{{X}}_{n - 1|n - 1}^{\left( i \right)}} \right) - {{{z}}_{n|n - 1}}} \right]^{\rm{T}}} + {{R}} \\ {{P}}_{n|n - 1}^{{xz}} =\, & \mathop \sum \limits_{i = 1}^N {w_i}\left[ {f\left( {{{X}}_{n - 1|n - 1}^{\left( i \right)}} \right) - {{{x}}_{n|n - 1}}} \right]\\ & \cdot {\left[ {h\left( {{{X}}_{n - 1|n - 1}^{\left( i \right)}} \right) - {{{z}}_{n|n - 1}}} \right]^{\rm{T}}} \end{split} \right\}$

(9)

(4) 增益计算与状态变量更新：利用第(3)步得到的协方差矩阵求得滤波的增益估计，如式(10)所示。利用卡尔曼增益与实时观测值，更新状态变量及其协方差矩阵，如式(11)和式(12)。

${K_n} = {{P}}_{n|n - 1}^{{{xz}}}{\left( {{{P}}_{n|n - 1}^{{{zz}}}} \right)^{ - 1}}\hspace{40pt}$

(10)

${{{x}}_{n|n}} = {{{x}}_{n|n - 1}} + {K_n}\left( {{{{z}}_n} - {{{z}}_{n|n - 1}}} \right)$

(11)

${{P}}_{n|n}^{{{xx}}} = {{P}}_{n|n - 1}^{{{xx}}} - {K_n}{{P}}_{n|n - 1}^{{{zz}}}K_n^{\rm{T}} \hspace{20pt}$

(12)

(5) 存储与迭代：将该时刻状态变量及其协方差存储至跟踪文件，返回第(2)步进行下一时刻迭代。

4. 实验与算法验证

4.1 仿真实验

仿真实验如3.1节中描述，磁偶极子目标运动轨迹如图1。用于实现目标跟踪的传感器阵列由3个矢量磁力仪构成，坐标分别为(–2 m, 0, 0), (0, 0, 0), (2 m, 0, 0)，依次标记为传感器1、传感器2、传感器3。通过传感器1、传感器2和传感器2、传感器3分别获得差量磁场，并在原始差量场上叠加高斯白噪声模拟工作的噪声。利用观测磁场，分别采用EKF, UKF与本文提出的基于MCKF跟踪算法，实现目标实时跟踪。为评估不同算法的跟踪性能，采用100次随机实验，将跟踪误差的统计平均值作为评判指标。不同仿真算法误差结果如图4所示，其中图4(a)为x轴投影轨迹的误差，图4(b)为y轴投影轨迹的误差，图4(c)为z轴投影轨迹的误差。

图 4 仿真磁偶极子跟踪结果对比

下载: 全尺寸图片幻灯片

从图4可以看出，EKF与UKF具有类似的跟踪效果，MCKF的初期跟踪效果较差，而后期跟踪效果明显优于EKF/UKF的跟踪结果。定性来讲，MCKF在收敛速度上较EKF/UKF更慢，但稳定后在跟踪效果上跟踪精度最高。为定量分析不同算法跟踪性能，计算3种跟踪方法跟踪误差结果如表1。

表 1 不同跟踪算法各方向投影轨迹跟踪误差(m)

	跟踪算法	时间区间
	跟踪算法	1～40点	41～80点	81～120点	120～160点	160～200点
x方向	EKF	0.0651	0.0127	0.0023	0.0057	0.0409
	UKF	0.0689	0.0127	0.0024	0.0057	0.0410
	MCKF	0.2143	0.0153	0.0024	0.0039	0.0218
y方向	EKF	0.0491	0.0113	0.0059	0.0100	0.0319
	UKF	0.0512	0.0112	0.0061	0.0096	0.0318
	MCKF	0.1237	0.0115	0.0042	0.0061	0.0158
z方向	EKF	0.0430	0.0089	0.0021	0.0042	0.0163
	UKF	0.0436	0.0088	0.0022	0.0041	0.0164
	MCKF	0.0456	0.0086	0.0022	0.0043	0.0198

下载: 导出CSV

| 显示表格

表1中列举不同时间区间，x, y, z等3个方向上投影轨迹跟踪误差。从表1中看出MCKF算法与EKF/UKF算法相比，在1～40点内跟踪误差较大，随后在40～80点内跟踪误差相当，而在稳定的80～200点内，MCKF的跟踪精度明显优于其余两种跟踪算法；表1结果表明，在z方向3种算法跟踪误差相当。

依据MCKF算法原理，蒙特卡洛样本点数会影响算法跟踪性能。为评估不同蒙特卡洛样本点数的影响，依次采用50, 100, 200, 400样本点进行仿真，结果如表2。从表2中得出随粒子数由50增至200，算法的跟踪误差逐步减小，说明增加粒子数可有效提升算法跟踪性能；而当粒子数由200增至400时，跟踪性能的提升并不明显。由于粒子数多少决定计算量，因此目标跟踪过程中应在满足性能时尽量减少粒子数。

表 2 不同蒙特卡洛样本点数的算法各方向投影轨迹跟踪误差(m)

表2(a) x方向
蒙特卡洛样本点数	时间区间
蒙特卡洛样本点数	1～40	41～80	81～120	120～160	160～200
50	0.2817	0.0242	0.0035	0.0052	0.0264
100	0.2481	0.0194	0.0030	0.0046	0.0231
200	0.2143	0.0153	0.0024	0.0039	0.0218
400	0.2316	0.0159	0.0023	0.0039	0.0212
表2(b) y方向
蒙特卡洛样本点数	时间区间
蒙特卡洛样本点数	1～40	41～80	81～120	120～160	160～200
50	0.1658	0.0168	0.0049	0.0069	0.0192
100	0.1442	0.0141	0.0047	0.0063	0.0166
200	0.1237	0.0115	0.0042	0.0061	0.0158
400	0.1346	0.0118	0.0043	0.0059	0.0151
表2(c) z方向
蒙特卡洛样本点数	时间区间
蒙特卡洛样本点数	1～40	41～80	81～120	120～160	160～200
50	0.0569	0.0124	0.0031	0.0053	0.0242
100	0.0496	0.0105	0.0026	0.0046	0.0202
200	0.0456	0.0086	0.0022	0.0043	0.0198
400	0.0466	0.0090	0.0021	0.0042	0.0195

下载: 导出CSV

| 显示表格

4.2 地磁背景跟踪实验

地磁背景跟踪实验在张家口张北县某地草原进行，选用Mag-03磁通门传感器，传感器分布与磁性目标运动轨迹示意如图5，其中模拟跟踪磁性目标体为通电螺旋圈如图6，磁矩约为60 A·m²。根据图5(b)中运动轨迹，磁性目标在运动过程中的拐点位置由于磁性源的姿态发生变化，此时目标分量磁矩会产生突变。

图 5 真实目标跟踪实验

下载: 全尺寸图片幻灯片

图 6 跟踪实验模拟磁性目标体

下载: 全尺寸图片幻灯片

传感器1、传感器2、传感器3测量的磁异常如图7(a)、图7(b)、图7(c)所示，传感器1、传感器2的差量磁场如图7(d)，传感器2、传感器3的差量磁场如图7(e)。差量磁异常信号与原始磁异常信号在幅度上相近，甚至大于原始测量磁异常。利用本节中提出的基于MCKF跟踪方法与传统的基于EKF/UKF的滤波跟踪方法对目标进行跟踪，目标位置跟踪结果如图8所示，其中图8(a)、图8(b)、图8(c)分别代表x, y, z的跟踪结果，黑色实线为参考轨迹；目标磁矩的跟踪结果如图9所示，其中图9(a)、图9(b)、图9(c)分别代表磁矩x, y, z分量磁矩的跟踪结果，黑色实线为参考磁矩。

图 7 阵列传感器测量磁场与差量磁场

下载: 全尺寸图片幻灯片

图 8 模拟磁性目标的位置跟踪结果

下载: 全尺寸图片幻灯片

图 9 模拟磁性目标的磁矩跟踪结果

下载: 全尺寸图片幻灯片

由图8和图9的跟踪结果可以看出：(1)在跟踪初始阶段及拐点位置，EKF与UKF算法跟踪响应速度要优于MCKF跟踪方法，与仿真结果一致；(2)当目标运动状态稳定后(图中3500点后)跟踪过程中，MCKF的跟踪性能明显优于其余两种算法。其中，EKF与UKF算法在4300点后，由于信噪比下降算法的跟踪轨迹发生明显偏离，而MCKF直到4800点后轨迹才逐渐发生偏移。综上所述，真实目标的跟踪结果与仿真结果表现一致，说明本文提出的算法具有较好的抗地磁干扰性能，在低信噪比时跟踪性能优于传统滤波方法。

5. 结束语

本文分析了地磁场背景下磁偶极子目标跟踪过程中存在背景干扰问题，提出利用传感器差量磁场信号实现目标跟踪的模型。新模型的优势有：(1)传感器阵列测量磁场差量与原始测量场相比磁场强度相似，因此可以保证足够远的跟踪距离；(2)测量磁场差量不易受到背景地磁场的干扰。同时，考虑差量磁场模型的高度非线性，本文提出基于MCKF的滤波跟踪算法，从而有效解决减小跟踪模型的非线性误差。通过仿真实验，结果表明基于新模型的跟踪算法较经典的UKF、EKF算法具有更高的跟踪精度；地磁背景跟踪实验，同样证实算法的有效性。

图 1 仿真磁偶极子运动轨迹与传感器分布位置

下载: 全尺寸图片幻灯片

图 2 传感器测量结果

下载: 全尺寸图片幻灯片

图 3 基于MCKF跟踪算法的流程图

下载: 全尺寸图片幻灯片

图 4 仿真磁偶极子跟踪结果对比

下载: 全尺寸图片幻灯片

图 5 真实目标跟踪实验

下载: 全尺寸图片幻灯片

图 6 跟踪实验模拟磁性目标体

下载: 全尺寸图片幻灯片

图 7 阵列传感器测量磁场与差量磁场

下载: 全尺寸图片幻灯片

图 8 模拟磁性目标的位置跟踪结果

下载: 全尺寸图片幻灯片

图 9 模拟磁性目标的磁矩跟踪结果

下载: 全尺寸图片幻灯片

表 1 不同跟踪算法各方向投影轨迹跟踪误差(m)

	跟踪算法	时间区间
	跟踪算法	1～40点	41～80点	81～120点	120～160点	160～200点
x方向	EKF	0.0651	0.0127	0.0023	0.0057	0.0409
	UKF	0.0689	0.0127	0.0024	0.0057	0.0410
	MCKF	0.2143	0.0153	0.0024	0.0039	0.0218
y方向	EKF	0.0491	0.0113	0.0059	0.0100	0.0319
	UKF	0.0512	0.0112	0.0061	0.0096	0.0318
	MCKF	0.1237	0.0115	0.0042	0.0061	0.0158
z方向	EKF	0.0430	0.0089	0.0021	0.0042	0.0163
	UKF	0.0436	0.0088	0.0022	0.0041	0.0164
	MCKF	0.0456	0.0086	0.0022	0.0043	0.0198

下载: 导出CSV

表 2 不同蒙特卡洛样本点数的算法各方向投影轨迹跟踪误差(m)

表2(a) x方向
蒙特卡洛样本点数	时间区间
蒙特卡洛样本点数	1～40	41～80	81～120	120～160	160～200
50	0.2817	0.0242	0.0035	0.0052	0.0264
100	0.2481	0.0194	0.0030	0.0046	0.0231
200	0.2143	0.0153	0.0024	0.0039	0.0218
400	0.2316	0.0159	0.0023	0.0039	0.0212
表2(b) y方向
蒙特卡洛样本点数	时间区间
蒙特卡洛样本点数	1～40	41～80	81～120	120～160	160～200
50	0.1658	0.0168	0.0049	0.0069	0.0192
100	0.1442	0.0141	0.0047	0.0063	0.0166
200	0.1237	0.0115	0.0042	0.0061	0.0158
400	0.1346	0.0118	0.0043	0.0059	0.0151
表2(c) z方向
蒙特卡洛样本点数	时间区间
蒙特卡洛样本点数	1～40	41～80	81～120	120～160	160～200
50	0.0569	0.0124	0.0031	0.0053	0.0242
100	0.0496	0.0105	0.0026	0.0046	0.0202
200	0.0456	0.0086	0.0022	0.0043	0.0198
400	0.0466	0.0090	0.0021	0.0042	0.0195

下载: 导出CSV

参考文献(23)

IOANNIDIS G. Identification of a ship or submarine from its magnetic signature[J]. IEEE Transactions on Aerospace and Electronic Systems, 1977, AES-13(3): 327–329. doi: 10.1109/TAES.1977.308404

FAGGIONI O, SOLDANI M, GABELLONE A, et al. Undersea harbour defence: A new choice in magnetic networks[J]. Journal of Applied Geophysics, 2010, 72(1): 46–56. doi: 10.1016/j.jappgeo.2010.07.001

ZHANG Mengying, WANG Hua, GE Lin, et al. Automatic search algorithms for near-field ferromagnetic targets based on magnetic anomaly detection[J]. Mathematical Problems in Engineering, 2018, 2018: 2130236. doi: 10.1155/2018/2130236

SHEINKER A, LERNER B, SALOMONSKI N, et al. Localization and magnetic moment estimation of a ferromagnetic target by simulated annealing[J]. Measurement Science and Technology, 2007, 18(11): 3451–3457. doi: 10.1088/0957-0233/18/11/027

SHEINKER A, SALOMONSKI N, GINZBURG B, et al. Remote sensing of a magnetic target utilizing population based incremental learning[J]. Sensors and Actuators A: Physical, 2008, 143(2): 215–223. doi: 10.1016/j.sna.2007.10.064

YANG Wan’an, HU Chao, LI Mao, et al. A new tracking system for three magnetic objectives[J]. IEEE Transactions on Magnetics, 2010, 46(12): 4023–4029. doi: 10.1109/tmag.2010.2076823

GAO Xiang, YAN Shenggang, and LI Bi. A novel method of localization for moving objects with an alternating magnetic field[J]. Sensors, 2017, 17(4): 923(1–12). doi: 10.3390/s17040923

WANG Chen, QU Xiaodong, ZHANG Xiaojuan, et al. A fast calibration method for magnetometer array and the application of ferromagnetic target localization[J]. IEEE Transactions on Instrumentation and Measurement, 2017, 66(7): 1743–1750. doi: 10.1109/TIM.2017.2668558

贾文抖, 林春生, 陈春行, 等. 针对磁目标定位失效的改进欧拉定位方法[J]. 海军工程大学学报, 2018, 30(3): 37–42.

JIA Wendou, LIN Chunsheng, CHEN Chunxing, et al. Improved Euler method for preventing failure of positioning magnetic target[J]. Journal of Naval University of Engineering, 2018, 30(3): 37–42.

WIEGERT R and OESCHGER J. Generalized magnetic gradient contraction based method for detection, localization and discrimination of underwater mines and unexploded ordnance[C]. The OCEANS 2005 MTS/IEEE, Washington, USA, 2005: 1325–1332.

BIRSAN M. Non-linear Kalman filters for tracking a magnetic dipole[R]. Defence R&D Canada, Atlantic, 2005.

BIRSAN M. Unscented particle filter for tracking a magnetic dipole target[C]. The OCEANS 2005 MTS/IEEE, Washington, USA, 2005, 1656–1659.

KOZICK R J and SADLER B M. Algorithms for tracking with an array of magnetic sensors[C]. The 5th IEEE Sensor Array and Multichannel Signal Processing Workshop, Darmstadt, Germany, 2008: 423–427.

ALHMIEDAT T, ABU TALEB A, and BSOUL M. A study on threads detection and tracking systems for military applications using WSNs[J]. International Journal of Computer Applications, 2012, 40(15): 12–18. doi: 10.5120/5055-7347

张朝阳, 肖昌汉, 高俊吉, 等. 磁性物体磁偶极子模型适用性的试验研究[J]. 应用基础与工程科学学报, 2010, 18(5): 862–868. doi: 10.3969/j.issn.1005-0930.2010.05.016

ZHANG Zhaoyang, XIAO Changhan, GAO Junji, et al. Experiment research of magnetic dipole model applicability for a magnetic object[J]. Journal of Basic Science and Engineering, 2010, 18(5): 862–868. doi: 10.3969/j.issn.1005-0930.2010.05.016

于振涛, 吕俊伟, 张本涛. 基于海底磁力仪阵列的磁性目标定位方法[J]. 武汉理工大学学报, 2012, 34(6): 131–135. doi: 10.3963/j.issn.1671-4431.2012.06.028

YU Zhentao, LÜ Junwei, and ZHANG Bentao. A method to localize magnetic target based on a seabed array of magnetometers[J]. Journal of Wuhan University of Technology, 2012, 34(6): 131–135. doi: 10.3963/j.issn.1671-4431.2012.06.028

吴志东, 周穗华, 陈志毅. 基于非线性滤波算法的磁偶极子跟踪[J]. 鱼雷技术, 2013, 21(4): 262–267. doi: 10.3969/j.issn.1673-1948.2013.04.006

WU Zhidong, ZHOU Suihua, and CHEN Zhiyi. Magnetic dipole tracking based on nonlinear filtering algorithm[J]. Torpedo Technology, 2013, 21(4): 262–267. doi: 10.3969/j.issn.1673-1948.2013.04.006

高俊吉, 刘大明, 周国华. 水中非合作运动磁性目标跟踪及参数估计[J]. 哈尔滨工程大学学报, 2013, 34(9): 1124–1130. doi: 10.3969/j.issn.1006-7043.201209033

GAO Junji, LIU Daming, and ZHOU Guohua. Study on the tracking and parameter estimating of unknown moving magnetism objects[J]. Journal of Harbin Engineering University, 2013, 34(9): 1124–1130. doi: 10.3969/j.issn.1006-7043.201209033

张宏欣, 周穗华, 吴志东, 等. 基于改进粗糙化粒子滤波的磁偶极子跟踪[J]. 华中科技大学学报: 自然科学版, 2014, 42(9): 76–80. doi: 10.13245/j.hust.140917

ZHANG Hongxin, ZHOU Suihua, WU Zhidong, et al. Magnetic dipole localization based on improved roughening particle filter[J]. Journal of Huazhong University of Science and Technology:Natural Science Edition, 2014, 42(9): 76–80. doi: 10.13245/j.hust.140917

周穗华, 张文成, 张宏欣. 基于混合卡尔曼滤波的磁偶极子目标跟踪[J]. 水雷战与舰船防护, 2015, 23(4): 7–11.

ZHOU Suihua, ZHANG Wencheng, and ZHANG Hongxin. Magnetic dipole target tracking based on mixed kalman filter[J]. Mine Warfare &Ship Self-defence, 2015, 23(4): 7–11.

吴垣甫, 孙跃. 基于递推更新卡尔曼滤波的磁偶极子目标跟踪[J]. 北京航空航天大学学报, 2017, 43(9): 1805–1812. doi: 10.13700/j.bh.1001-5965.2016.0694

WU Yuanfu and SUN Yue. Magnetic dipole target tracking based on recursive update Kalman filter[J]. Journal of Beijing University of Aeronautics and Astronautics, 2017, 43(9): 1805–1812. doi: 10.13700/j.bh.1001-5965.2016.0694

张宏欣, 周穗华, 张伽伟. 磁偶极子跟踪的渐进贝叶斯滤波方法[J]. 自动化学报, 2017, 43(5): 822–834. doi: 10.16383/j.aas.2017.c160052

ZHANG Hongxin, ZHOU Suihua, and ZHANG Jiawei. A progressive bayesian filtering approach to magnetic dipole tracking[J]. Acta Automatica Sinica, 2017, 43(5): 822–834. doi: 10.16383/j.aas.2017.c160052

张宏欣. 磁性目标跟踪的多模型自适应滤波方法[J]. 数字海洋与水下攻防, 2018, 1(2): 57–62.

ZHANG Hongxin. Multiple-model adaptive filtering method for magnetic target tracking[J]. Digital Ocean &Underwater Warfare, 2018, 1(2): 57–62.

施引文献

期刊类型引用(7)

1.	陆斌杰，张晓兵，戴忠华. 基于椭球体阵列的磁目标跟踪系统能观性研究. 兵器装备工程学报. 2025(01): 157-164 . 百度学术
2.	朱明慧，谷鹏，高东 . 一种基于递推最小二乘法的运动平台剩磁差分测量与估计方法. 空间科学学报. 2024(03): 570-584 . 百度学术
3.	陆斌杰，张晓兵，戴忠华. 基于电流反演的螺线管磁场模拟能力分析. 兵器装备工程学报. 2024(09): 138-144 . 百度学术
4.	王莉君，唐骞，李鹏. 多物体遮挡下基于深度学习的目标跟踪仿真. 计算机仿真. 2023(09): 176-179+309 . 百度学术
5.	胡寅，张文英，朱荣. 埋地金属管线主动磁定位探测器设计与研究. 现代电子技术. 2022(24): 48-53 . 百度学术
6.	汪兴海，于力，周镇宇，王玉明. 基于矢量磁力仪阵列的水下磁性目标警戒系统. 兵工自动化. 2021(07): 78-83 . 百度学术
7.	戴忠华，周穗华，张晓兵. 多目标优化的舰船磁场建模方法. 物理学报. 2021(16): 147-159 . 百度学术

其他类型引用(7)

资源附件(0)

访问统计

图(9) / 表(2)

计量

文章访问数: 3705
HTML全文浏览量: 1228
PDF下载量: 134
被引次数: 14

1. 引言
2. 磁偶极子运动模型
3. 基于差量磁场的MCKF跟踪算法
3.1 差量磁场信号
3.2 基于MCKF的跟踪算法
4. 实验与算法验证
4.1 仿真实验
4.2 地磁背景跟踪实验
5. 结束语

姓名
邮箱
手机号码
标题
留言内容
验证码

留言板

地磁背景下基于传感器阵列的磁偶极子目标跟踪方法

doi: 10.11999/JEIT190236

通讯作者:
朱万华　whzhu@mail.ie.ac.cn

计量

Magnetic Dipole Object Tracking Algorithm Based on Magnetometer Array in Geomagnetic Background

1. 引言

2. 磁偶极子运动模型

3. 基于差量磁场的MCKF跟踪算法

3.1 差量磁场信号

3.2 基于MCKF的跟踪算法

4. 实验与算法验证

4.1 仿真实验

4.2 地磁背景跟踪实验

5. 结束语

期刊类型引用(7)

其他类型引用(7)

计量

目录

1. 引言

2. 磁偶极子运动模型

3. 基于差量磁场的MCKF跟踪算法

3.1 差量磁场信号

3.2 基于MCKF的跟踪算法

4. 实验与算法验证

4.1 仿真实验

4.2 地磁背景跟踪实验

5. 结束语

留言板

地磁背景下基于传感器阵列的磁偶极子目标跟踪方法

doi: 10.11999/JEIT190236

通讯作者: 朱万华 whzhu@mail.ie.ac.cn

计量

出版历程

Magnetic Dipole Object Tracking Algorithm Based on Magnetometer Array in Geomagnetic Background

1. 引言

2. 磁偶极子运动模型

3. 基于差量磁场的MCKF跟踪算法

3.1 差量磁场信号

3.2 基于MCKF的跟踪算法

4. 实验与算法验证

4.1 仿真实验

4.2 地磁背景跟踪实验

5. 结束语

期刊类型引用(7)

其他类型引用(7)

计量

出版历程

目录

1. 引言

2. 磁偶极子运动模型

3. 基于差量磁场的MCKF跟踪算法

3.1 差量磁场信号

3.2 基于MCKF的跟踪算法

4. 实验与算法验证

4.1 仿真实验

4.2 地磁背景跟踪实验

5. 结束语

通讯作者:
朱万华　whzhu@mail.ie.ac.cn