CN113238208B - 非规则轨迹水中运动目标前向声散射多普勒频移计算方法 - Google Patents

Abstract

本发明涉及海洋声学技术领域，尤其涉及水中运动目标前向声散射多普勒频移计算方法。非规则轨迹水中运动目标前向声散射多普勒频移计算方法，包括：(1)将发射端、目标轨迹和接收端的地球经纬度坐标转换为直角坐标，并依次标记为(x_S,y_S)、(x_T,y_T)和(x_R,y_R)；(2)将目标轨迹坐标序列记为(x_T(t_i),y_T(t_i)),i＝1,2,…,N，计算所有时段内的平均航速

(3)对目标轨迹进行平滑处理，获得虚拟直线轨迹，分别计算目标和接收端到其虚拟直线轨迹与基线交点的水平距离、虚拟直线轨迹与基线的夹角；(4)利用步骤(2)和(3)计算得到的参数，计算各时刻目标的前向声散射多普勒频移。本发明解决了现有技术无法计算目标运动轨迹非规则时的前向声散射多普勒频移的问题。

Description

非规则轨迹水中运动目标前向声散射多普勒频移计算方法

技术领域

本发明涉及海洋声学技术领域，尤其涉及水中运动目标前向声散射多普勒频移计算方法。

背景技术

在前向声散射探测系统中，由于发射端和接收端的位置均固定而只有目标运动，直达波的多普勒频移为零而目标具有时变的多普勒频移。从接收声信号中提取目标的前向散射多普勒频移信息，可以探测到入侵收发连线的水中运动目标。然而，目标前向散射多普勒频移提取结果的真实性和准确性，需要通过相应的理论计算结果作对比检验。

水中目标的前向声散射多普勒频移的理论计算方法有两种：

第一种方法主要利用双基地分置角、目标航向与双基地分置角平分线的夹角这两个参数，其表达式为：

其中v为目标运动速度，λ为声源信号的波长，β为双基地分置角，δ为目标航向与双基地分置角平分线的夹角。该方法本质上是一种几何描述，其优势是能够计算目标任意运动形式下的前向声散射多普勒频移，但也存在两个明显不足之处。一个是难以直观地建立目标多普勒频移与目标参数的对应关系，而且对多普勒频移关于目标运动时间的演变历程的描述不够清晰。另一个是要求记录目标航向角，否则在只给定目标轨迹坐标和航速的条件下，通常需要多层计算才能得到角度δ的值。显然，数值计算层数越多，误差积累越大，计算结果的精度也就越差。

第二种方法的计算原理与第一种方法的原理相同，但用代数描述替换了几何描述，直接建立了目标的前向声散射多普勒频移关于目标运动时间、穿越收发连线时刻、穿越点位置、航向角等参数的表达式，即：

式中t为目标运动时间、t_c为目标穿越收发连线的时刻、d为收发连线长度、d_TR为目标穿越点到接收端的水平距离、α为目标航迹与基线的夹角。该方法的突出优势是能够直接描述目标的前向散射多普勒频移与各参数的依赖关系，并能够直观地刻画多普勒频移关于目标运动时间的演变历程，且不需要目标的航向记录。然而，该方法的缺点是只适用于目标轨迹为直线的情形，并不适合非规则目标轨迹的情形。

目前尚未有针对非规则轨迹的水中目标前向声散射多普勒频移的计算方法。

发明内容

本发明的目的是为了解决现有技术存在的上述缺陷，提供一种针对非规则轨迹的水中目标前向声散射多普勒频移的计算方法，

本发明解决其技术问题采用的技术方案是：非规则轨迹水中运动目标前向声散射多普勒频移计算方法，包括：(1)将发射端、目标轨迹和接收端的地球经纬度坐标转换为直角坐标，并依次标记为(x_S,y_S)、(x_T,y_T)和(x_R,y_R)；

(2)将目标轨迹坐标序列记为(x_T(t_i),y_T(t_i)),i＝1,2,…,N，计算所有时段内的平均航速

(3)对目标轨迹进行平滑处理，获得虚拟直线轨迹，分别计算目标和接收端到其虚拟直线轨迹与基线交点的水平距离、虚拟直线轨迹与基线的夹角；

(4)利用步骤(2)和(3)计算得到的参数，计算各时刻目标的前向声散射多普勒频移。

进一步优选地，所述平滑处理的方法为：采用长度为l_w点的时间窗并以步长t_w点依次截取目标轨迹坐标序列，将每个窗口内轨迹坐标的算术平均值作为该时段内目标的等效轨迹坐标，平滑后的目标轨迹坐标为：

进一步优选地，目标和接收端到其虚拟直线轨迹与基线交点的水平距离的计算方法为：

在任意t_i时刻，目标到其虚拟直线轨迹与基线交点

的水平距离为：

目标虚拟直线轨迹与基线的交点到接收端的水平距离为：

目标当前位置到接收端的水平距离为：

进一步优选地，虚拟直线轨迹与基线的夹角的计算方法为：应用余弦定理，得到t_i时刻目标虚拟直线轨迹与基线的夹角：

进一步优选地，目标的前向声散射多普勒频移的计算方法包括：

在t_i时刻，将计算得到的参数：

d_TC(t_i)、d_TR(t_i)、α(t_i)及信号中心频率f_c、环境声速c，一并代入线性轨迹的前向声散射多普勒频移参数公式：

遍历平滑后的所有目标轨迹点，获得多普勒频移关于目标运动时间的变化曲线。

本发明的非规则轨迹水中运动目标前向声散射多普勒频移计算方法，具有的有益效果是：

解决了现有技术无法计算目标运动轨迹非规则时的前向声散射多普勒频移的问题；在对目标轨迹作平滑处理时采用等时间间隔方位采样并作算术平均计算，过程严谨；对估计目标虚拟直线航迹与基线的夹角时依据余弦定理，并利用几何关系计算目标、轨迹交点和接收端两两之间的水平距离，具有严格理论依据；在计算前向声散射多普勒频移时依据严格的参数方程，得到的结果准确。

附图说明

图1为本发明实施例提供的非规则轨迹水中运动目标前向声散射多普勒频移计算方法流程图；

图2目标轨迹分段线性近似示意图；

图3试验中某航程的发射端、接收端位置以及目标轨迹；

图4目标运动轨迹分段估计示意图；

图5平滑前后的目标轨迹对比曲线，其中，(a)为平滑前；(b)为平滑后；

图6目标虚拟直线轨迹与基线的交点分布图；

图7目标到其虚拟直线轨迹与基线交点的水平距离；

图8目标在基线的虚拟穿越点到接收端的水平距离；

图9目标的虚拟直线轨迹与基线的夹角信息；

图10目标前向声散射多普勒频移的估算结果；

图11前向声散射多普勒频移的近似计算与试验数据结果对比。

具体实施方式

为了便于理解本发明，下面结合附图和具体实施例，对本发明进行更详细的说明。附图中给出了本发明的较佳的实施例。但是，本发明可以以许多不同的形式来实现，并不限于本说明书所描述的实施例。相反地，提供这些实施例的目的是使对本发明公开内容的理解更加透彻全面。

实施例1本发明提供的非规则轨迹水中运动目标前向声散射多普勒频移计算方法，主要流程如图1所示，具体包括以下步骤：

(1)目标航迹坐标变换

在二维平面内选择合适的参考位置为原点建立直角坐标系，将发射端、接收端和目标轨迹的地球经纬度坐标转换为直角坐标，并依次标记为(x_S,y_S)、(x_T,y_T)和(x_R,y_R)，其中的基线长度(收发连线长度)可表示为：

(2)估算目标运动速度

将目标轨迹坐标序列记为(x_T(t_i),y_T(t_i)),i＝1,2,…,N，则目标在任意两个相邻时刻t_i和t_i+1之间的运动速度可直接求得：

对所有时段内的航速求算术平均，得

将该平均速度作为计算目标前向声散射多普勒频移时的目标运动速度。

(3)计算平滑处理后的目标轨迹坐标

采用长度为l_w点的时间窗并以步长t_w点依次截取目标轨迹坐标序列，将每个窗口内轨迹坐标的算术平均值作为该时段内目标的等效轨迹坐标，即可实现对目标原始轨迹的平滑，以减弱局部时段内目标轨迹剧烈摇摆对计算结果的影响。平滑后的目标轨迹坐标为：

(4)计算目标轨迹与基线的交点坐标

对于平滑后的目标轨迹采用分段线性近似处理，即假设目标在任意两个相邻时刻t_i和t_i+1之间的运动轨迹为直线，如图2所示，这一虚拟直线轨迹必然会与基线交于某点。显然，在目标运动历程的不同时刻，目标的虚拟直线轨迹会与基线相交于多个不同的点。

利用目标在t_i和t_i+1两个时刻的位置坐标

和

求得目标在此时段内的虚拟直线轨迹方程：

同时，由发射端和接收端的位置坐标可得基线所在的直线方程：

联立(5)、(6)两式，得到t_i时刻目标虚拟直线轨迹与基线的交点坐标为：

(5)计算目标、轨迹交点和接收端两两之间的水平距离

在任意t_i时刻，目标到其虚拟直线轨迹与基线交点

的水平距离为：

目标虚拟直线轨迹与基线的交点到接收端的水平距离为：

目标当前位置到接收端的水平距离为：

(6)计算目标虚拟直线轨迹与基线的夹角

应用余弦定理并联立(9)、(10)和(11)式，得到t_i时刻目标虚拟直线轨迹与基线的夹角：

(7)计算目标的前向声散射多普勒频移

在t_i时刻，将对目标真实轨迹记录处理后的参数：

d_TC(t_i)、d_TR(t_i)、α(t_i)连同信号中心频率f_c、环境声速c，一并代入线性轨迹的前向声散射多普勒频移参数公式：

遍历平滑后的所有目标轨迹点，依次经过步骤(4)-(7)的计算，得到各时刻目标的前向声散射多普勒频移，最终得到对应的多普勒频移关于目标运动时间的变化曲线。

实施例2利用2011年千岛湖缩尺目标前向声散射探测试验的目标航迹数据和目标前向声散射多普勒频移提取结果对本发明的算法进行验证。

试验中的发射端和接收端分置、间距约为1100米。目标为含泡沫夹层的双层铝板组合体(无自动力)，由拖船通过长度为10米的柔性绳拖拽穿越收发连线。发射端、接收端位置以及目标航迹如图3所示。试验中声源发射信号的中心频率f_c＝10千赫兹，试验水层的声速近似恒定，声速值c≈1485米/秒。受风和水流影响，拖船的运动轨迹并非直线且存在较显著的偏转和摇摆。试验中采用便携式GPS设备记录拖船的真实轨迹位置信息，采样间隔为1秒。

(1)估算目标航速

将(2)式和(3)式依次应用于图3的目标轨迹序列，得到航速估计结果示于图4，其中点划线表示实测的拖船航速，速度值的起伏振荡源于拖船轨迹的摇摆。实线为求算术平均后的结果，约为0.5米/秒，将该值作为目标的运动速度。

(2)目标轨迹平滑

采用时间长度为5秒的矩形窗以3秒为步长对拖船位置点序列滑动截取，并以每次截取的五个位置点的算术平均值作为目标在该时刻的位置，从而实现对拖船原始轨迹的平滑处理。图5(a)和图5(b)分别示出了平滑前、后的拖船轨迹，可见平滑处理能够减弱局部时段的剧烈震荡。如前所述，因拖船与目标通过柔性绳连接，故可把平滑后的轨迹视为目标的真实轨迹。

(3)目标轨迹与基线的交点坐标分布

图6给出了利用图5(b)的目标轨迹信息经(7)式和(8)式计算得到的目标虚拟直线轨迹与基线的交点分布(*号)。结果表明：由于目标轨迹并非直线形式，经分段线性近似处理后各时刻的虚拟轨迹与基线的交点位置不同于目标的真实穿越位置，但也大都聚集在真实位置附近区域。

(4)目标到轨迹交点的水平距离

图7给出了由(9)式得到的整个航程中目标到其虚拟直线轨迹与基线交点的水平距离计算结果。显然，目标在基线的虚拟穿越点与实际穿越点的位置偏差越大，则d_TC的取值也越大。但是总体的轨迹演变趋势是清晰的，尤其当目标位于基线附近时，距离波动较小。

(5)轨迹交点和接收端的水平距离

由(10)式得到的计算结果示于图8，目标在基线的虚拟穿越点到接收端的水平距离变化趋势与图6中的虚拟穿越点在基线上的分布形式保持一致。对比图8和图7可知：目标位置越是靠近基线，各个距离参量的起伏程度也就越小。

(6)目标虚拟直线轨迹与基线的夹角

图9给出了(12)式计算的目标虚拟直线轨迹与基线的夹角信息。当目标刚好穿越收发连线时，目标航迹与基线的夹角会变为其补角。因而，在目标穿越基线的时刻航迹夹角会存在一定程度的跃变。目标航向越是偏离正橫方向(α＝90°)，跃变量就越大。但是，这里的角度跃变并不影响对目标前向散射多普勒频移的计算结果。

(7)计算目标的前向声散射多普勒频移

将从目标真实轨迹记录中提取的参数：

d_TC(t)、d_TR(t)、α(t)连同信号中心频率f_c、环境声速c，一并代入(13)式，得到的前向声散射多普勒频移曲线示于图10。

(8)数值计算与数据分析结果比对

由于发射端和接收端的位置均固定，直达波的多普勒频移为零。水中运动目标穿越收发连线时目标的前向声散射波存在时变的多普勒频移。但由于前向声散射波通常与直达波干涉在一起，最终效果是使得接收声场发生低频调制现象。提取接收信号包络并对该包络作时频分析得到了图11所示的多普勒频移条纹。穿越收发连线的目标散射波对应了图中“V”型多普勒频移条纹，直达波则对应了多普勒频移近似为零的水平条纹。

本方法得到的多普勒频移曲线如图11的点划线所示，该曲线所表征的多普勒频移范围与变化趋势均同数据分析得到的多普勒频移条纹一致，结果证实了本发明方法的可行性。

Claims (1)

1.非规则轨迹水中运动目标前向声散射多普勒频移计算方法，其特征在于，包括：

(1)将发射端、目标轨迹和接收端的地球经纬度坐标转换为直角坐标，并依次标记为(x_S,y_S)、(x_T,y_T)和(x_R,y_R)；

(2)将目标轨迹坐标序列记为(x_T(t_i),y_T(t_i)),i＝1,2,…,N，计算所有时段内的平均航速

目标在任意两个相邻时刻t_i和t_i+1之间的运动速度可直接求得：

对所有时段内的航速求算术平均，得

将该平均速度作为计算目标前向声散射多普勒频移时的目标运动速度；

(3)对目标轨迹进行平滑处理，获得虚拟直线轨迹，分别计算目标和接收端到其虚拟直线轨迹与基线交点的水平距离、虚拟直线轨迹与基线的夹角；

所述平滑处理的方法为：采用长度为l_w点的时间窗并以步长t_w点依次截取目标轨迹坐标序列，将每个窗口内轨迹坐标的算术平均值作为该时段内目标的等效轨迹坐标，平滑后的目标轨迹坐标为：

目标和接收端到其虚拟直线轨迹与基线交点的水平距离的计算方法为：

在任意t_i时刻，目标到其虚拟直线轨迹与基线交点

的水平距离为：

目标虚拟直线轨迹与基线的交点到接收端的水平距离为：

目标当前位置到接收端的水平距离为：

虚拟直线轨迹与基线的夹角的计算方法为：应用余弦定理，得到t_i时刻目标虚拟直线轨迹与基线的夹角：

(4)利用步骤(2)和(3)计算得到的参数，计算各时刻目标的前向声散射多普勒频移：

在t_i时刻，将计算得到的参数：

d_TC(t_i)、d_TR(t_i)、α(t_i)及信号中心频率f_c、环境声速c，一并代入线性轨迹的前向声散射多普勒频移参数公式：

其中，d_SR表示基线长度；

遍历平滑后的所有目标轨迹点，获得多普勒频移关于目标运动时间的变化曲线。

Images