一种多参量水下磁目标体定位方法技术

本发明专利技术提供一种多参量水下磁目标体定位方法，将水下磁性目标体视为静止的磁偶极子，其中心视为笛卡尔坐标系原点。拖曳系统沿测线工作时同时测得三分量磁异常数据、两个总场磁异常数据和三轴姿态角数据，利用三轴姿态角数据对三分量磁异常数据进行误差补偿。将三分量磁异常数据分别用最近距离点的坐标、水下磁目标体的磁矩矢量、以及磁测系统速度矢量和时间t0来表示，建立关于时间的函数。利用两个总场数据形成对三分量函数表达式的约束条件，结合粒子群算法粗略求解函数表达式的系数矩阵，再用L

【技术实现步骤摘要】
一种多参量水下磁目标体定位方法


[0001]本专利技术属于水下磁目标定位领域，尤其是指一种基于粒子群算法和L
‑
M混合算法的多参量水下磁目标体。

技术介绍

[0002]水下磁目标探测为军事探浅，定位未爆炸物以及海底管道检测等提供了手段手段。随着高精度磁力仪的发展，水下拖曳平台所获得的磁传感器数据也越来越丰富，配合广泛应用的智能算法，对水下磁目标体的定位精度也逐渐提高。
[0003]现有的水下拖曳平台通常配备两个总场磁力仪在水平或者竖直方向构成总场梯度测量，或者仅采用单个三分量磁力仪完成测量任务。现有的方法的优点在于采用总场磁力仪测量时，对拖曳运动时产生的姿态旋转与震动不敏感，磁测数据较为可靠。缺点在于总场磁力仪只能获得单个磁场数据，想要准确定位水下磁目标体必须增加磁力仪数量或者沿着不同测线反复测量。然而三分量磁力仪能够同时测得更加丰富的磁场数据，但由于它对旋转震动较为敏感，目前很少有研究将其布置于水下拖曳平台上。定位算法方面，L
‑
M算法在各种领域被广泛应用，但是缺点在于对初值的要求较高，一旦给定远离理想值的初值，便会产生迭代发散现象，影响磁目标体的定位精度，而粒子群算法具有收敛速度快，搜索全剧最优解强的优点。

技术实现思路

[0004]针对上述现有水下磁目标体定位方法的局限性，尤其是现有方法采用总场磁力仪所测磁异常数据信息较少，且需要沿不同测线反复测量导致的效率低下问题，提供一种多参量水下磁目标体定位方法。
[0005]本专利技术是这样实现的，r/>[0006]一种多参量水下磁目标体定位方法，该方法包括：
[0007]a.拖曳式磁测系统配备有两个光总场磁力仪布置在水平方向，三分量磁力仪、惯性测量单元和定位系统布置于总场磁力仪轴线中心后方，假设磁性目标体位于系统下方某处，并以其中心点为系统坐标原点，拖曳式磁测系统沿着某条测线以速度V直线行驶，测线上每点坐标由定位系统获得，建立磁测系统工作时的测量模型，并保存测量得到的三分量磁场数据；
[0008]b.利用惯性测量单元所测得的数据对三分量磁异常数据进行误差补偿，拖曳式磁测系统沿测线运行时在x、y、z轴上产生的姿态角变化分别为α、β、γ，将地理坐标系上的磁场数据B
g
转化为载体坐标系上的磁场数据B
c
，则在每个数据时间点上其补偿公式如下所示：
[0009][0010]c.假设在磁目标体与测线距离最近的点处的最短距离为r0，磁测系统在XOY平面
上的速度分量为v
x
和v
y
，则建立磁测系统测线上任何一点(x，y，z0)与磁性目标体的距离r的表达式，该表达式是关于时间t的函数，如下所示：
[0011]r(t)＝(x+v
x
t，y+v
y
t，z0)
[0012]d.磁性目标体产生的磁异常在测线上某点处的数据值B表示为：
[0013][0014]式中M为磁性目标体的磁矩矢量，R为磁性目标体与测线上测量点的距离矢量，r为磁性目标体与测线上测量点的距离标量，μ0＝4π
×
10
‑7H/m为真空磁导率；
[0015]e.结合步骤c、步骤d分别建立三分量磁常数据B
x
,B
y
,B
z
关于最近距离点时刻t0、最近距离点坐标(x0,y0,z0)、磁目标体的磁矩矢量(M
x
,M
y
,M
z
)、磁测系统的速度矢量(v
x
,v
y
)的表达式，该表达式包含一个系数矩阵C并且是关于时间t的函数，表示为如下:
[0016][0017]式中B
x
,B
y
,B
z
为磁场数值的三个分量，(M
x
,M
y
,M
z
)为磁目标体的磁矩矢量，系数矩阵C中各元素表示为下式：
[0018][0019]式中(x，y，z0)为测线上任意一点的坐标值，v
x
和v
y
磁测系统在XOY平面上的速度分量，μ0＝4π
×
10
‑7H/m为真空磁导率；
[0020]f.由总场磁力仪记录的磁偶极子信号，近似地表示为B在周围地球磁场B
E
上的投影,三分量磁异常数据通过下式转化为总场磁异常数据T，表示为如下：
[0021][0022]式中λ1、λ2、λ3可由当地地磁倾角I和地磁偏角D求得，具体求解公式为：
[0023]λ1＝cosIcosD
[0024]λ2＝cosIsinD
[0025]λ3＝sinI
[0026]将其整理可得：
[0027]T＝b0M
x
+b1M
y
+b2M
z
[0028]其中(M
x
,M
y
,M
z
)为磁目标体的磁矩矢量，(b0，b1，b2)为
系数向量；
[0029]g.根据磁测系统中传感器的布置方式，一个总场磁力仪的数据表达式中系数b0、b1、b2和另一个总场磁力仪的磁异常数据表达式中的b'0、b'1、b'2通过曲线拟合确定，并且系数的大小与距离的高次方成反比，通过两个总场磁力仪数据确定三分量磁异常数据中系数矩阵部分元素的取值范围，如下式所示：
[0030][0031]h.经过步骤g取值范围的约束后，通过三分量磁异常数据的测量值和计算值之间的差异，建如下目标函数：
[0032][0033]式中B
x
,B
y
,B
z
是三分量磁异常数据的测量值，B
xcal
,B
ycal
,B
zcal
是三分量磁异常数据的计算值，
[0034]i.利用粒子群算法最小化目标函数，获得三分量磁异常数据的系数矩阵的粗略解。
[0035]j.将粒子群算法求解出的系数矩阵作为L
‑
M算法的初值，精确求解系数矩阵，其中将所求系数矩阵内的9个元素定义为9维变量x,通过下式迭代公式求解最优解
[0036]x
k+1
≈x
k
‑
[J
T
(x
k
)J(x
k
)+μ
k
I]‑1J
T
(x
k
)v(x
k
)
[0037]式中x
k
为当前位置，x
k+1
为下一步位置，J为雅可比矩阵，μ
k
I为调整因子，当调整因子大时相当于梯度下降法，小时相当于高斯牛顿法；
[0038]k.获得系数矩阵的精确解后，通过解超定非线性方程组的方本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种多参量水下磁目标体定位方法，其特征在于，该方法包括：a.拖曳式磁测系统配备有两个总场磁力仪布置在水平方向，三分量磁力仪、惯性测量单元和定位系统布置于总场磁力仪轴线中心后方，假设磁性目标体位于系统下方某处，并以其中心点为系统坐标原点，拖曳式磁测系统沿着某条测线以速度V直线行驶，测线上每点坐标由定位系统获得，建立磁测系统工作时的测量模型，并保存测量得到的三分量磁场数据；b.利用惯性测量单元所测得的数据对三分量磁异常数据进行误差补偿，拖曳式磁测系统沿测线运行时在x、y、z轴上产生的姿态角变化分别为α、β、γ，将地理坐标系上的磁场数据B
g
转化为载体坐标系上的磁场数据B
c
，则在每个数据时间点上其补偿公式如下所示：c.假设在磁目标体与测线距离最近的点处的最短距离为r0，磁测系统在XOY平面上的速度分量为v
x
和v
y
，则建立磁测系统测线上任何一点(x，y，z0)与磁性目标体的距离r的表达式，该表达式是关于时间t的函数，如下所示：r(t)＝(x+v
x
t，y+v
y
t，z0)d.磁性目标体产生的磁异常在测线上某点处的数据值B表示为：式中M为磁性目标体的磁矩矢量，R为磁性目标体与测线上测量点的距离矢量，r为磁性目标体与测线上测量点的距离标量，μ0＝4π
×
10
‑7H/m为真空磁导率；e.结合步骤c、步骤d分别建立三分量磁常数据B
x
,B
y
,B
z
关于最近距离点时刻t0、最近距离点坐标(x0,y0,z0)、磁目标体的磁矩矢量(M
x
,M
y
,M
z
)、磁测系统的速度矢量(v
x
,v
y
)的表达式，该表达式包含一个系数矩阵C并且是关于时间t的函数，表示为如下:式中B
x
,B
y
,B
z
为磁场数值的三个分量，(M
x
,M
y
,M
z
)为磁目标体的磁矩矢量，系数矩阵C中各元素表示为下式：
式中(x，y，z0)为测线上任意一点的坐标值，v
x
和v
y
磁测系统在XOY平面上的速度分量，μ0＝4π
×
10
‑7H/m为真空磁导率；f.由总场磁力仪记录的磁偶极子信号，近似地表示为B在周围地球磁场B
E
上的投影,三分量磁异常数据通过下式转化为总场磁异常数据T，表示为如下：式中λ1、λ2、λ3可由当地地磁倾角I和地磁偏角D求得，具体求解公式为：λ1＝cosIcosDλ2＝cosIsinDλ3＝sinI将其整理可得：T＝b0M
x
+b1M
y
+b2M
z
其中(M
x
,M
y
,M
z
)为磁目标...

【专利技术属性】
技术研发人员：赵静，王震，林君，郭言国，刘生威，吴永鹏，
申请(专利权)人：南方海洋科学与工程广东省实验室湛江，
类型：发明
国别省市：