舰船下方磁场的信号重构门限值参数优选方法

摘要

本申请属于舰船磁场信号重构方法技术领域，尤其涉及一种舰船下方磁场的信号重构门限值参数优选方法。包括如下步骤：测量获得舰船在不同深度下的磁场时域分布曲线；获取不同截止频率对应的不同深度最佳采样间隔，取同一深度下磁场三份量对应采样间隔的最小值作为Nyquist间隔，通过深度换算得出换算精度最佳的间隔，并确定其对应的最佳截止频率门限值；并以精度最高的重构曲线对应的截止频率门限值Q为基准确定最佳采样频率。该方法对目标形状鲁棒性好，水上和水下舰船目标皆适用，在较大采样间隔条件下能完成较高精度的信号重构。

说明书

技术领域

本申请属于舰船磁场信号重构方法技术领域，尤其涉及一种舰船下方磁场的信号重构门限值参数优选方法。

背景技术

通过测量得到有限的舰船磁场数据反演其等效源强度，或推算不同深度的磁场数据是研究舰船磁场特征分布规律的重要手段。目前国内对于舰船目标的反演建模都是建立在已知其较为完整的磁场通过特性的前提下进行，虽能实现对舰船等效源的模拟，但建模采用的磁场信号采样间隔较小，包含了完整的磁场信息，这对实际采样过程中采样点数目有限、有效采样间隔范围的选取等工作的参考性不足，随着现代舰船磁性隐身技术的提高，需要一种更好的磁场信号重构方法，以完成在最大采样间隔条件下对舰船磁场反演。

发明内容

本申请的目的在于，提供一种目标形状鲁棒性好，水上和水下舰船目标皆适用，在较大采样间隔条件下能完成较高精度的信号重构，当由浅向深进行换算时，换算误差小的舰船下方磁场的信号重构门限值参数优选方法。

为实现上述目的，本申请采用如下技术方案。

一种舰船下方磁场的信号重构门限值参数优选方法，包括如下步骤：

步骤一、测量获得舰船在不同深度下的磁场时域分布曲线；

步骤二、获取不同截止频率对应的不同深度最佳采样间隔，具体包括：

B1.定义Q为截止频率门限值，根据能谱幅值B(f)与最大能谱B(f)

步骤三、基于不同截止频率门限值对舰船磁场信号进行重构及深度换算，确定最佳截止频率门限值；

C1.取同一深度下磁场三份量对应采样间隔的最小值作为Nyquist间隔，计算不同Q取值对应的不同深度最佳采样间隔；针对不同的间隔对舰船磁场信号进行重构，通过深度换算得出换算精度最佳的间隔，并确定其对应的最佳截止频率门限值；

C2.分别计算不同截止频率门限值Q下不同深度对应的磁场分布；将磁场分布与对应深度的原始磁场采集数据进行换算误差的计算，对重构曲线进行深度换算的精度进行验证，并以精度最高的重构曲线对应的截止频率门限值Q为基准确定最佳采样频率。

对前述舰船下方磁场的信号重构门限值参数优选方法的进一步改进或者优选方案，所述步骤一具体包括：

A1.采集舰船磁场三分量B

A2.求平均值得到舰船磁场三分量的平均值

A3.两者相减得到去除直流分量的舰船磁场分布

A4.对B

A5.获得不同深度下的最大能谱

对前述舰船下方磁场的信号重构门限值参数优选方法的进一步改进或者优选方案，所述步骤B1具体是指，对于满足

对前述舰船下方磁场的信号重构门限值参数优选方法的进一步改进或者优选方案，所述步骤C1具体包括：

4.1.根据不同截止频率门限值Q对应的不同深度最大采样间隔在原始采集磁场中等间隔取样，对完成截取的能量谱进行离散傅里叶逆变换，得到不同截止频率门限值Q对应的不同深度的重构曲线；

4.2.基于前述重构曲线，得到不同截止频率门限值Q对应的磁场分布；

采取磁体模拟法，将船模简化为均匀磁化的旋转椭球体和磁偶极子混合模型进行深度换算；均匀磁化的旋转椭球体的几何尺寸与舰船相当，将N个磁偶极子以相同的距离d沿着椭球体的纵向进行排列，长半轴为a，短半轴为b，椭球体的几何中心和舰船的几何中心重合，从左到右分别记为1，2...N；第i个磁偶极子的坐标为：

当1≤i≤N时，

通过已知的舰船磁场分布，就可以计算出每一个磁偶极子的磁矩和位置得到椭球体的磁性分布，计算舰船周围空间的磁场；根据浅深度的磁场信号的得到条件系数矩阵F和磁偶极子的磁矩M；根据系数矩阵F计算处磁偶极子的位置坐标，结合磁矩M得到模拟体的磁性分布，进行深度换算，得到不同深度的磁场分布B

4.3.将不同截止频率门限值Q对应的磁场分布与对应深度的原始磁场采集数据进行比较，确定不同深度下最佳的截止频率门限值Q。

对前述舰船下方磁场的信号重构门限值参数优选方法的进一步改进或者优选方案，所述C2中，取磁场三分量中相对残差最大值max{RRE

其有益效果在于：

本申请的舰船下方磁场的信号重构门限值参数优选方法，为了保证对舰船磁场信号的准确重构及换算的准确性，在缩比模型试验数据的基础上，对磁性处理后的船模磁场时频特性进行分析，基于起截止频率与最大能谱的关系，以及采样间隔与船长间的数学模型，提出适用于磁场的Nyquist间隔计算方法，结合磁体模拟法进行磁场深度换算完成对本文方法的验证。该方法对目标形状鲁棒性好，水上和水下舰船目标皆适用，在较大采样间隔条件下能完成较高精度的信号重构，当由浅向深进行换算时，换算误差皆小于7％，当换算的基准深度不小于3倍船宽时，能完成由浅向深和由深向浅双向的磁场换算，换算误差均小于10％。

附图说明

图1是船模实验布置示意图；

图2是不同深度下的磁场分布；

图3是不同深度下的磁场能量谱分布(x方向)；

图4是不同深度下的磁场能量谱分布(y方向)；

图5是不同深度下的磁场能量谱分布(z方向)；

图6是磁场三分量的累积能量与总能量之比随频率变化曲线图(x方向)

图7是磁场三分量的累积能量与总能量之比随频率变化曲线图(y方向)

图8是磁场三分量的累积能量与总能量之比随频率变化曲线图(z方向)

图9是均匀磁化的椭球体与磁偶极子混合模型

图10是算例不同深度下的磁场分布(x方向)

图11是算例不同深度下的磁场分布(y方向)

图12是算例不同深度下的磁场分布(z方向)

图13是算例不同深度下的磁场能量谱分布(x方向)；

图14是算例不同深度下的磁场能量谱分布(y方向)；

图15是算例不同深度下的磁场能量谱分布(z方向)。

具体实施方式

以下结合具体实施例对本申请作详细说明。

步骤一、测量获得舰船在不同深度下的磁场时域分布曲线；

为了明确舰船磁场的谱特征，得到磁性分布曲线与能量谱分布的关系，基于采样定理，对不同深度下舰船磁场信号三分量进行能量谱分析，提取频域特征量，得出舰船目标的磁性特征分布与频谱能量之间的关系

具体步骤包括：

采集舰船磁场三分量B

本实施例中，以典型潜艇模型作为试验目标(长l＝6M，宽B＝0.6m)，船模在传感器正上方以10cm/s匀速通过，采取动态连续采样，测量距离18m，设置采样频率f

表1不同深度下舰船磁场三分量最大能谱对应频率表

由前述内容可知，当测量深度小于3B时，舰船磁场三分量的最大能谱频率相同，当测量深度大于3B时，舰船磁场纵向分量B

步骤二、获取不同截止频率对应的不同深度最佳采样间隔；

由图6～8可知，舰船匀速通过时，舰船磁场95％的能量集中在0.05Hz以下，近98％的能量集中在0.1Hz以下。为了有效选取合适的截止频率，舍弃高频成分，使剩余能谱能重构完整原始信号的同时减少运算量和高频噪声的引入，建立起截止频率f

基于上述分析，设初始采样频率为f

表2截止频率门限值Q取不同值时，最小采样频率、采样距离与船长之间的关系

取同一深度下，磁场三份量对应采样间隔的最小值作为Nyquist间隔，得出不同Q取值对应的不同深度最佳采样间隔如表3所示。

表3不同截止频率门限值Q对应的不同深度最大采样间隔表

步骤三、基于不同截止频率门限值对舰船磁场信号进行重构及深度换算，确定最佳截止频率门限值；基于前述步骤，获得了多个基于不同截止频率门限值对应的不同深度最大采样间隔表，为确定不同间隔对磁场换算精度的影响，针对不同的间隔对舰船磁场信号进行重构，并通过深度换算得出换算精度最佳的间隔，并确定其对应的最佳截止频率门限值；具体步骤包括：

3.1.根据不同截止频率门限值Q对应的不同深度最大采样间隔在原始采集磁场中等间隔取样，对完成截取的能量谱进行离散傅里叶逆变换，得到不同截止频率门限值Q对应的不同深度的重构曲线，即得到不同截止频率门限值Q对应的不同深度磁场分布的基准值；

3.2.基于前述重构曲线，得到不同截止频率门限值Q对应的磁场分布；

为了实现对舰船磁场信号进行重构，需要通过深度换算得出换算精度最佳的磁场Nyquist间隔。本实施例中，采取磁体模拟法，将船模简化为均匀磁化的旋转椭球体和磁偶极子混合模型进行深度换算。假设均匀磁化的旋转椭球体的几何尺寸与舰船相当，如图9所示，将N个(N为奇数)磁偶极子以相同的距离d沿着椭球体的纵向进行排列，长半轴为a，短半轴为b，椭球体的几何中心和舰船的几何中心重合，从左到右分别记为1，2...N；那么第i个磁偶极子的坐标为：

当1≤i≤N时，M

通过已知的舰船磁场分布，就可以计算出每一个磁偶极子的磁矩和位置，从而得到椭球体的磁性分布，从而计算舰船周围空间的磁场，实现深度换算。

由式(9)可知，要获得磁偶极子的位置坐标，需要求解系数矩阵F，采用遗传算法对稀疏矩阵的条件数进行求解，矩阵条件数的优化问题表达式及其约束方程为：

算法中具体参数设定如表4所示。

表4算法参数设定表

根据浅深度的磁场信号的得到条件系数矩阵F和磁偶极子的磁矩M；根据系数矩阵F计算处磁偶极子的位置坐标，结合磁矩M得到模拟体的磁性分布，进行深度换算，得到不同深度的磁场分布B

3.3.将不同截止频率门限值Q对应的磁场分布与对应深度的原始磁场采集数据进行比较，确定不同深度下最佳的截止频率门限值Q；

基于前述步骤，分别计算不同截止频率门限值Q下不同深度对应的磁场分布B

具体实施时，采用相对残差来表示换算误差，通过二范数来定义相对残差RRE：

表5不同门限值Q对应的不同深度之间换算误差表

由计算结果可得：

(1)换算误差随着Q的取值减小而减小，当门限值Q＝0.0001时，换算误差最小；

(2)门限值Q＝0.001对应的采样间隔为Q＝0.0001的三倍，但Q＝0.0001时换算误差无显著改善，因此当门限值取Q＝0.001能在较少测量点的条件下实现较高的换算精度；

(3)对于同一基准深度来说，换算误差随着换算深度的增加而增大。这是由于信号重构后，基准深度本身存在一定误差，在进行深度换算时也会产生进一步的误差叠加造成的；

(4)从浅向深换算的误差远小于从深向浅，这是由于浅深度的磁场局部信息更复杂导致的；

(5)当门限值Q＝0.001、基准深度为3B至4B范围时(数据已标蓝)，能实现此范围内任意深度之间的高拟合度换算，换算误差皆小于10％。

为了验证结果的可靠性，本文在对潜艇模型进行分析后，采用相同方法对消磁后某型水面舰艇进行试验，(模型长度l＝4.82m、船宽B＝0.53m，其余实验布置与潜艇试验相同)，试验对所提方法的有效性进行了检验，试验过程中，获取消磁后船模在1.5B、2.5B、3B、3.5B、4B和4.5B深度下的原始磁场分布，进行信号的特征提取和重构，并采用本文反演算法得到磁源分布，最后通过深度换算对方法进行验证。

值得注意的是，由于在实际舰船探测中，无法控制船体的航向，无法进行舰船固定磁性的分离，因此，试验将基于实际采集的数据进行。试验中磁通门传感器分别置于1.5B、2.5B、3B、3.5B、4B和4.5B六个不同深度，电机拖动装置控制船模以10cm/s的速度匀速通过传感器上方，测量距离为18m，实时测量舰船磁性分布，采样频率f＝10Hz。

图10～12为船模通过时测得的不同深度磁场信号分布，图13～15为船模磁场的能量谱。为了便于观察曲线的变化，磁场分布数据采用归一化处理，能量谱采用局部放大处理。

最后应当说明的是，以上实施例仅用以说明本申请的技术方案，而非对本申请保护范围的限制，尽管参照较佳实施例对本申请作了详细地说明，本领域的普通技术人员应当理解，可以对本申请的技术方案进行修改或者等同替换，而不脱离本申请技术方案的实质和范围。