基于线阵综合声速补偿的近场波达方向估计方法及装置的制造方法
【技术领域】
[0001] 本发明涉及目标定位的技术领域，尤其是指一种基于线阵综合声速补偿的近场波 达方向估计方法及装置。
【背景技术】
[0002] 目标定位的方法有很多，针对不同应用场景有不同的定位方法。在陆地上，一般采 用卫星定位系统、雷达等技术来定位目标。但是对于广阔的江河、海洋等水领域下，由于电 磁波在水下衰减严重，而声波是水下传播最远的辐射形式。所以一般采用声呐技术来定位 目标。根据收发方式的不同，水声定位系统分为主动式探测和被动式探测。目前主要有三种 方法：一、基于波束形成的定位方法，这类方法对阵列接收的声波信号滤波、加权求和，调节 接收阵列的方向，并在整个接收空间扫描，具有最大能量的方向就是目标的方位角；二、基 于时延估计的定位方法，这类方法先计算声源分别到达阵列上两个阵元的时间差，然后再 结合阵列的几何结构估计出声源的位置；三、基于高分辨率的空间谱估计方法，这类算法应 用广，分辨率高，它基本思想是对阵列接收数据协方差矩阵进行特征分解，划分为两个相互 正交的子空间：信号子空间和噪声子空间，最后进行谱峰搜索确定方位角。但是现有的近场 波达方向（简称：D0A)估计算法中，如经典MUSIC算法，均采用声速在每条声传播路径上都 是相同的，即具有相同的波长的假设，忽略了水流等实际环境因素对声音传播速度的影响， 在复杂的水流环境下易产生较大的探测误差。
[0003] 为了实现高精度探测定位，本专利考虑每一条传播路径上的水流速分量的影响， 将感知的水流速用于补偿静止水流时的声速，形成空变综合声速，研宄一种基于线阵综合 声速补偿的近场D0A估计方法。
[0004] 因此，人们需要一种基于线阵综合声速补偿的近场波达方向估计方法及装置来 满足需求。
【发明内容】
[0005] 本发明的目的在于克服现有技术的不足，提供一种基于线阵综合声速补偿的近场 波达方向估计方法及装置，可实现精准、快速定位，实用性高。
[0006] 为实现上述目的，本发明所提供的技术方案其基于线阵综合声速补偿的近场波达 方向估计方法，包括以下步骤：
[0007] 1)推导综合声速补偿的近场波达方向估计模型
[0008] 模型场景中有一均匀直线阵，有M个接收阵元，阵元间间距为d，被动式探测，有 K个目标声源Si，S2，…，SK为点源，频率为f，目标对应的方位角及距离参数分别为（0p A)，（02，r2)，…，（0K，rK)，其中0i是声源阵列的参考点之间的连线与阵元所在直线 之间的夹角，9mi是声源Si与第m个阵元之间的连线与阵元所在直线之间的夹角，dm是第 m个阵元与参考点之间的距离，ri是声源Si与阵列参考点之间的距离，i= 1，2,…，K;设水 流方向与阵元所在直线之间的夹角为a，水流速度为vwatOT，声音在静止的水流中的速度设 为V ;对模型详细推导如下：
[0009] 1. 1)第i个声源到第m个阵元的距离为：
[0011] 其中，dm= [m_(M+l)/2]d ;
[0012] 1. 2)声源Si与第m个阵元之间的连线与阵元所在直线之间的夹角0 mi为：
[0014] 1. 3)根据几何关系，可以求得第i个声源Si到达阵元参考点与到达第m个阵元所 需的时间之差Tmi，求得：
[0017]因此，时间差Tmi可以这样表达：
[0019] 1.4)对于K个声源信号，第m个阵元所接收到的全部信号为：
[0021] 其中，叫（〇是接收第i个声源时的噪声，a时是幅度衰减参数，且
[0023] 1. 5)由上述四步可以得出整个阵列第t个时刻的接收信号为：
[0024] X(t) = AS(t)+N(t)
[0025] 其中，阵列流形矩阵i
[0026] 2)在步骤1)推导的模型下具体近场波达方向估计方法估计出方位角和距离位置 参数，定位目标，详细步骤如下：
[0027] 2. 1)对接收M阵源的阵列信号做L次的快拍采样，得到MX L维的数据矩阵X;
[0028] 2. 2)利用最大似然估计，计算出数据X协方差矩阵为：
[0029]Rx=XXh/L
[0030] 2. 3)对数据协方差矩阵艮进行特征值分解，得到M个特征值入= 1，...，M)和 相对应的特征矢量ei:
[0031] Rxei=入々，i = 1，2,…，M
[0032] 2. 4)划分信号子空间和噪声子空间，对艮的特征值进行升序排列 {44，…，4,4+1，…，4}，即…>4，其前k个小特征值允12，…，4对应的特征 矢量组成信号子空间％ =叵,§2，…,^]，后m-k个大特征值if+15...对应的特征矢 量组成噪声子空间％ =[^+1，…，U;
[0033] 2. 5)构造新模型下的新的空间谱函数，即构造步骤1)推导的近场波达方向估计 模型的空间谱函数，因为新模型是有考虑到水流速度的影响，所以求出的时间差Tmi与传 统模型有所不同，也就使得新模型下所求得的阵列流形矩阵A自然与传统理想场景的A有 所不同的，也即新模型下会有新的空间谱函数，对于第i个声源，新的空间谱函数为：
[0035] 其中a( 0 i，ri)是导向矢量，即阵列流形矩阵A的列矢量；
[0037] 2. 6)进行方向和距离联合二维谱峰搜索，确定第K个声源的方向和距离信息：
[0038]
[0039] 本发明所述的基于线阵综合声速补偿的近场波达方向估计装置，包括控制模块、 发射模块、接收模块、处理模块、显示模块、无线通信模块，其中：
[0040] 所述控制模块，分别与发射模块、接收模块、处理模块、显示模块，无线通信模块相 连，用于对各个模块进行控制；
[0041] 所述发射模块，与控制模块相连，根据控制模块的指令，进行超声发射；
[0042] 所述接收模块，分别与控制模块、处理模块相连，根据控制模块的指令，接收目标 声源信号并传送给处理模块；
[0043] 所述处理模块，分别与控制模块、接收模块、显示模块相连，根据控制模块的指令 进行数据处理，其通过对接收信号进行分析，运用相关算法计算出目标的距离和方位角位 置参数；
[0044] 所述显示模块，分别与控制模块、处理模块相连，提供人机交互界面，供用户输入 设定参数，并根据控制模块的指令，将处理模块计算出的目标位置参数显示出来；
[0045] 所述无线通信模块，与控制模块相连，提供水下与船舶或陆地之间数据传送的功 能。
[0046] 本发明与现有技术相比，具有如下优点与有益效果：
[0047] 1、本发明相对传统的水下目标定位方法更具有实用性，精度高。传统的方法在每 一条传播路径上都使用了相同的声速，忽略了水流影响导致的综合声速的改变，在复杂的 水流环境下会产生较大的误差，定位精度不高，不具有实用性，而本发明方法采用的是综合 的空变声速，考虑到水流速的补偿，具有实用性和高精度的定位。
[0048] 2、本发明是在改进传统方法的前提下，算法复杂度与传统方法相当。因为估计算 法仍旧可以采用各种经典的波达方向估计算法，算法的复杂度与传统的波达方向估计方法 相当。
[0049] 3、本发明装置可行性强、安装简单。利用测水流速的应用已经很广泛，技术相当成 熟，而且，本发明使用的处理器等芯片，集成度高，计算能力强大，保证了本发明的可行性。
【附图说明】
[0050] 图1为本发明所述装置的硬件结构框图。
[0051] 图2为本发明所述方法的总体流程图。
[0052] 图3为本发明所述水下近场波达方向估计D0A模型示意图。
[0053] 图4为基于本发明所述模型的具体D0A方法流程图。
[0054] 图5为本发明所述装置的总体结构示意框图。
[0055] 图6为方向和距离联合二维谱峰搜索流程图。
【具体实施方式】
[0056] 下面结合具体实施例对本发明作进一步说明。
[0057] 本实施例所述的基于线阵综合声速补偿的近场波达方向估计方法，总体流程如图 2所示，包含以下两大步骤，而每一大步骤均分子步骤实现，下面分别对这两大步骤进行详 细说明：
[0058] 1)推导综合声速补偿的近场波达方向估计模型
[0059] 在水下目标定位实际应用中，大多数情况都是基于近场模型的，此时声源达到接 收阵列的信号应该是球面波，考虑到有水流的影响，则近场D0A模型场景示意图如图3所 不〇
[0060] 模型场景中有一均匀直线阵，有M个接收阵元，阵元间间距为d，被动式探测，有 K个目标声源Si，S 2，…，SK为点源，频率为f，目标对应的方位角及距离参数分别为（0 p A)，（92，r2)，…，（0K，rK)，其中0i是声源阵列的参考点（这里选定为阵列中心）之 间的连线与阵元所在直线之间的夹角，0 mi是声源S ,与第m个阵元之间的连线与阵元所在 直线之间的夹角，dm是第m个阵元与参考点之间的距离，ri是声源Si与阵列参考点之间的 距离，i = 1，2,…，K ;设水流方向与阵元所在直线之间的夹角为a，水流速度为VwatOT，声音 在静止的水流中的速度设为V。对模型详细推导如下：
[0061] 1. 1)第i个声源到第m个阵元的距离为：
[0063] 其中，dm= [m_(M+l)/2]d ;
[0064] 1. 2)声源Si与第m个阵元之间的连线与阵元所在直线之间的夹角0 mi为：
[0066] 1. 3)根据几何关系，可以求得第i个声源Si到达阵元参考点与到达第m个阵元所 需的时间之差Tmi，求得：
[0069] 因此，时间差tmi可以这样表达：
[0071] 1.4)对于K个声源信号，第m个阵元所接收到的全部信号为：
[0073] 其中，叫⑴是接收第i个声源时的噪声，a时是幅度衰减参数，且
[0075] 1. 5)由上述四步可以得出整个阵列第t个时刻的接收信号为：
[0076] X(t) = AS(t)+N(t)
[0077] 其中，阵列流形矩阵
[0078] 上述步骤详细地推导了一种考虑水流影响的近场波达方向估计D0A模型，可以发 现新场景下阵列流形矩阵A与传统理想场景的A是有所不同的，原因是对于特定某个入射 角度9 i，考虑到水流速度在声音在阵元参考点和阵元m两条传播路径的投影分量的不同， 导致时间差Tmi的求解有所不同，新模型下所求得的阵列流形矩阵A自然与传统理想场景 的A是有所不同的。因此，新建立的模型相对于传统