一种基于信道匹配的阵列阵元位置校准方法及系统

1.本发明属于水声信号处理领域，具体涉及一种基于信道匹配的阵列阵元位置校准方法及系统。


背景技术：

2.海底柔性水平阵列在布放过程中，受海浪海流、布放船机动等因素影响，各阵元位置较预设位置往往存在偏差。阵元位置坐标的准确与否对目标探测与定位性能有着显著影响，相关研究表明，如果想控制阵形参数误差带来的阵增益损失在一分贝以内，则需要阵元位置偏差不大于所处理信号波长的十分之一。为提高后续阵列信号处理算法的探测性能，有必要进行阵形校准工作。
3.阵形校准方法大致可以分为两类：一类为无源校准，此类方法通过噪声源进行估计校准，利用海洋环境噪声或者过往船只噪声即可，但是只能反演出阵列的形状，即两两阵元间的距离信息，还需要其他信息辅助对初步结果进行平移、旋转、翻转等操作才能获取准确位置。另一类为有源校准，例如基于时延估计的阵形校准方法，通过配置校准声源，分析各阵元间信号的时延关系，并结合声源与参考阵元的位置信息，估计所有阵元相对参考阵元的位置，但局部最优解会受信道起伏影响，导致测量结果存在误差、并且有误差累积效应，阵形估计结果的准确性和稳定性不高；全局最优解则要求初始位置在真实值附近，否则校准结果也容易陷入局部最优值、导致较大误差；子空间分解的拟合算法对声源的信噪比与快拍数有一定要求，否则会引起阵列协方差矩阵估计值的误差，从而影响校准精度。
4.上述无论无源校准还是有源校准方法，都涉及到参考声速c0的选取问题，通常假设c0为频率不变的常数，并选择为接收阵深度处的介质声速。然而根据简正波理论可知，浅海波导环境属于频散介质，阵列接收到的声信号是不同阶模态的叠加，各阶模态的相速度、群速度不同，且会随频率变化，因此准确的声速参考值往往难以选择。因此在浅海水声信道下，会由于参考声速与实际声传播速度的偏差，给阵元位置校准结果带来系统误差，且阵列孔径越大影响越显著。


技术实现要素：

5.本发明的目的在于克服传统阵形校准方法由于参考声速与实际声传播速度的偏差，给阵元位置校准结果带来系统误差的缺陷。
6.为了实现上述目的，本发明提出了一种基于信道匹配的阵列阵元位置校准方法，所述方法包括：布设呈一定张角且分时工作的校准声源s1与s2，根据浅海信道环境计算可得第一阶、第二阶模态的截至频率，选择信号的频带范围介于第一阶模态与第二阶模态的截至频率间，即低频取值需要大于第一阶模态的截止频率，高频取值需要小于第二阶模态的截止频率，此时信号仅含有第一阶模态的信息，此步骤有效消除了浅海信道的模间频散的不利影响。带通滤波后将各阵元接收时域信号进行傅里叶变换后转为频域，再根据信道环境计算得到该波导条件下模态对应的相速度，然后基于随频率变化的相速度开展待校准阵
元与参考阵元间的距离差搜索，此步骤有效消除了浅海信道的模内频散造成的不利影响。再结合辅助声源、参考阵元位置建立代价函数，最后通过优化算法得到各个阵元相对参考阵元的位置信息。
7.作为上述方法的一种改进，所述方法具体包括：
8.步骤1：根据波导环境计算第一阶模态的截止频率f1与第二阶模态的截止频率f2，选择频带范围介于第一阶模态与第二阶模态的截至频率间，即低频取值需要大于第一阶模态的截止频率，高频取值需要小于第二阶模态的截止频率，此时信号仅含有第一阶模态的信息，此步骤有效消除了浅海信道的模间频散的不利影响。对参考阵元a1与待校准阵元ai接收到的时域信号进行带通滤波后得到s
ik
(t),k＝1,2，通过傅里叶变换转换得到频域信号
9.其中，k取值1或2代表声源s1、s2，f为频率，t为时间；
10.步骤2：计算p1(f)、pi(f)最大相关系数对应的即为ai较a1与声源sk的距离差估计值：
[0011][0012]
其中，δr为距离值，c
p
(f)为根据信道环境计算得到的该波导条件下第一阶模态对应的相速度，f
l
表示选取频段的最低频率，fh表示选取频段的最高频率，i表示阵元序号，i表示待校准阵元个数，上标*表示取共轭；
[0013]
步骤3：建立坐标(xi,yi)搜索代价函数如下：
[0014][0015]
其中，(x1,y1)为校准声源s1坐标位置，(x2,y2)为校准声源s2坐标位置；
[0016]
步骤4：通过优化搜索得到ai点坐标估计值：
[0017][0018]
作为上述方法的一种改进，所述优化搜索为模拟退火法。
[0019]
本发明还提供一种基于信道匹配的阵列阵元位置校准系统，所述系统包括：
[0020]
第一计算模块，用于根据波导环境计算第一阶模态的截止频率f1与第二阶模态的截止频率f2，选择频带范围介于第一阶模态与第二阶模态的截至频率间，即低频取值需要大于第一阶模态的截止频率，高频取值需要小于第二阶模态的截止频率，此时信号仅含有第一阶模态的信息，对参考阵元a1与待校准阵元ai接收到的时域信号进行带通滤波后得到时域信号通过傅里叶变换转换得到频域信号
[0021]
第二计算模块，用于计算ai较a1与声源sk的距离差估计值：
[0022]
第三计算模块，用于建立坐标(xi,yi)搜索代价函数；
[0023]
第四计算模块，用于通过优化搜索得到ai点坐标估计值。
[0024]
作为上述系统的一种改进，所述优化搜索为模拟退火法。
[0025]
本发明还提供一种计算机设备，包括存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现如上述任一项所述的方法。
[0026]
本发明还提供一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序当被处理器执行时使所述处理器执行如上述任一项所述的方法。
[0027]
与现有技术相比，本发明的优势在于：
[0028]
1、克服了传统阵形校准方法参考声速与实际声传播速度的偏差，给阵元位置校准结果带来系统误差的问题，校准所得阵元位置结果更加准确；
[0029]
2、当阵元间距很大时，本方法依然能够准确计算阵元位置。
附图说明
[0030]
图1所示为基于信道匹配的阵列阵元位置校准方法流程图；
[0031]
图2所示为阵形校准方法配置图；
[0032]
图3所示为校准声源布放位置示意图；
[0033]
图4所示为仿真波导环境示意图；
[0034]
图5所示为仿真信道下的第一阶模态相速度随频率的变化图；
[0035]
图6所示为接收信号波形图；(a)：s1声源，(b)：s2声源；
[0036]
图7所示为阵元位置校准结果图。
具体实施方式
[0037]
下面结合附图对本发明的技术方案进行详细的说明。
[0038]
现有的阵元位置校准方法，由于参考声速c0的选取往往与实际声传播速度存在偏差，导致阵元位置校准结果存在系统误差，对于大孔径阵列这种误差会更加显著。本发明通过充分利用浅海信道特征以克服上述技术缺陷，提出了一种基于信道匹配的阵列阵元位置校准方法。
[0039]
阵形校准方法配置如图2所示，s1、s2分别为校准声源的位置，a1为布设于海底的绝对同步水听器阵列参考阵元的位置，不失一般性，设定a1位置为直角坐标系原点，ai为第i个待校准阵元的位置。
[0040]
已知校准声源s1(x1,y1)、校准声源s2(x2,y2)与参考阵元a1的位置信息，声源s1、声源s2分时工作，参考阵元a1接收到信号待校准阵元ai接收到信号k取值1或2代表不同声源，阵元ai与阵元a1的水听器记录系统时钟绝对同步，求解ai点坐标(xi,yi)以实现阵元位置校准。
[0041]
如图1所示，具体方法如下：
[0042]
1、浅海波导环境具有复杂的频散效应，对于波导内的声传播，有模内频散与模间频散现象，这些都会导致信号波形的失真，再进行信号处理时需要结合信道特征予以补偿。首先，为了消除模间频散的影响，需要选取低频段仅激发一阶模态的信号，即根据波导环境计算第一阶模态的截止频率f1与第二阶模态的截止频率f2，选择频带范围介于第一阶模态与第二阶模态的截至频率间，即低频取值需要大于第一阶模态的截止频率，高频取值需要小于第二阶模态的截止频率，此时信号仅含有第一阶模态的信息，有效地消除了模间频散
的不利因素。对参考阵元a1与待校准阵元ai接收到的时域信号进行带通滤波后得到通过傅里叶变换转换得到频域信号
[0043]
其中，k取值1或2代表声源s1、s2，为第i号水听器与声源k的水平距离，f为频率，t为时间。
[0044]
2、为了消除模内频散的影响，需要选取与信道特征相匹配的相速度作为参考声速，按照如下公式计算p1(f)、pi(f)最大相关系数对应的即为ai较a1与声源sk的距离差估计值：
[0045][0046]
其中，δr为距离值，c
p
(f)为根据信道环境计算得到的该波导条件下第一阶模态对应的相速度，此参考声速的选取是与信道特征相匹配的，可以有效克服信道的频散特征，f
l
表示选取频段的最低频率，fh表示选取频段的最高频率，i表示阵元序号，i表示待校准阵元个数，上标*表示取共轭。
[0047]
3、进一步建立坐标(xi,yi)搜索代价函数如下：
[0048][0049]
其中，(x1,y1)为校准声源s1坐标位置，(x2,y2)为校准声源s2坐标位置。
[0050]
4、通过模拟退火等优化搜索可得ai点坐标估计值：
[0051][0052]
下面结合附图对校准方法的理论进行详细说明。
[0053]
1、校准声源布设图
[0054]
阵形校准方法配置如图2所示，s1、s2分别为校准声源的位置，a1为布设于海底的绝对同步水听器阵列参考阵元的位置，不失一般性，设定a1位置为直角坐标系原点，ai为某待校准阵元的位置，阵形校准问题描述如下：
[0055]
已知校准声源s1(x1,y1)、校准声源s2(x2,y2)与参考阵元a1的位置信息，声源s1、声源s2分时工作，参考阵元a1接收到信号待校准阵元ai接收到信号k取值1或2代表不同声源，阵元ai与阵元a1的水听器记录系统时钟绝对同步，求解ai点坐标(xi,yi)以实现阵元位置校准。
[0056]
2、信号模型
[0057]
根据简正波理论，水平不变的分层海洋波导环境中，深度为zs的点声源激发的远场声场可近似表示为m阶模态的线性叠加：
[0058][0059]
其中,s(f)为声源频谱，km(f)与βm(f)分别为第m阶模态的水平波数与指数衰减因
子，ψm为第m阶模态的深度函数，ρ(zs)为声源处水介质密度，z为水听器深度，r为水听器与声源的水平距离，m为波导中传播的模态个数。
[0060]
模态个数m受波导环境与声源频率制约，以液态海底均匀浅海波导为例，当声源频率f＞fm时，第m阶及其以下各阶次的模态才会被激发：
[0061][0062]
其中，c2为海底声速，c1为海水声速，h为水深。
[0063]
根据信道环境，通过公式(2)可以计算得到f1，f2。选择频带范围[f
l
,fh]，满足f1＜f
l
＜fh＜f2时，公式(1)可简化为：
[0064][0065]
其中,k(f)与β(f)分别为第一阶模态的水平波数与指数衰减因子，ψ为第一阶模态的深度函数。
[0066]
表征等相位面传播速度的相速度c
pm
(f)与模态的阶次和频率有关，第一阶模态的相速度c
p
(f)为：
[0067][0068]
设定δti为第i号阵元相对参考阵元的时钟偏移，则第i号水听器接收到的频域信号可以表示为：
[0069][0070]
其中，i为阵元个数，ri为第i号水听器与声源的水平距离。
[0071]
由于阵元ai与阵元a1的水听器记录系统时钟绝对同步，时域信号可以表示为：
[0072][0073]
结合公式(3)与公式(5)，阵元ai接收到k声源信号的频域形式可以表示为：
[0074][0075]
其中，为阵元ai较参考阵元a1距第k声源的水平距离差。
[0076]ai
较a1与声源sk的水平距离差估计值可使用相速度开展频域补偿估计得到：
[0077][0078]
其中，上标*表示取共轭。
[0079]
令δr
ik
为第i个阵元较参考阵元与第k声源的水平声程差，有如下关系：
[0080][0081]
参考阵元a1与待校准阵元ai接收到的时域信号通过傅里叶变换转换
到频域得到按照如下公式计算p1(f)、pi(f)最大相关系数对应的即为ai较a1与声源sk的水平距离差估计值：
[0082][0083]
进一步建立坐标(xi,yi)搜索代价函数如下：
[0084][0085]
通过模拟退火等优化搜索可得ai点坐标估计值：
[0086][0087]
下面通过实例验证基于信道匹配的阵列阵元位置校准方法与传统阵元位置校准方法的结果差异。
[0088]
为了量化所提方法与传统基于时延估计的阵形校准方法的优劣，定义阵元位置校准结果的均方根误差为：
[0089][0090]
其中，为校准估计得到的各阵元位置坐标，(x
0i
,y
0i
)为各阵元位置坐标的真实值。
[0091]
如图3所示，设定水平阵参考阵元(首阵元)为坐标系原点，自首阵元至末阵元方向为x轴正方向，逆时针旋转90度为y轴正方向。校准声源布放位置分别在s1[0
°
，30km]、s2[90
°
，20km]。
[0092]
如图4所示为仿真波导环境，参数设置具体如下：水层介质密度为1g/cm3,海深为96m，海水的声吸收衰减系数设为0db/λ；海底假设为平坦半无限液态基底,声速为1630m/s，密度为1.70g/cm3，衰减系数为0.1db/λ。声源深度为50m，接收阵为布放于海底的等间距60m的16元水平阵。分析信号的频带范围为15-30hz，频率分辨率为0.2hz，此频带在上述参数的仿真信道下远距离传播后仅存在一阶模态，使用kraken程序计算声场，得到该波导条件下模态对应的相速度如图5所示。
[0093]
校准声源分时工作后，各阵元接收信号的时域波形如图6所示。
[0094]
阵形坐标校准结果与真实位置坐标分别如图7所示。其中，信道匹配阵形校准方法选择参考声速为图5所示的随频率变化的相速度，而时延估计阵形校准方法选择参考声速为接收阵列处声速1526m/s。
[0095]
统计上述两种方法的阵形校准误差，如表1所示，表格第二行所示为全部阵元的校准误差
circuit，asic)、现场可编程门阵列(field programmable gate array，fpga)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件。可以实现或者执行上述公开的各方法、步骤及逻辑框图。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。结合上述公开的方法的步骤可以直接体现为硬件译码处理器执行完成，或者用译码处理器中的硬件及软件模块组合执行完成。软件模块可以位于随机存储器，闪存、只读存储器，可编程只读存储器或者电可擦写可编程存储器、寄存器等本领域成熟的存储介质中。该存储介质位于存储器，处理器读取存储器中的信息，结合其硬件完成上述方法的步骤。
[0107]
可以理解的是，本发明描述的这些实施例可以用硬件、软件、固件、中间件、微码或其组合来实现。对于硬件实现，处理单元可以实现在一个或多个专用集成电路(application specific integrated circuits，asic)、数字信号处理器(digital signal processing，dsp)、数字信号处理设备(dsp device，dspd)、可编程逻辑设备(programmable logic device，pld)、现场可编程门阵列(field-programmable gate array，fpga)、通用处理器、控制器、微控制器、微处理器、用于执行本技术所述功能的其它电子单元或其组合中。
[0108]
对于软件实现，可通过执行本发明的功能模块(例如过程、函数等)来实现本发明技术。软件代码可存储在存储器中并通过处理器执行。存储器可以在处理器中或在处理器外部实现。
[0109]
本发明还可提供一种非易失性存储介质，用于存储计算机程序。当该计算机程序被处理器执行时可以实现上述方法实施例中的各个步骤。
[0110]
最后所应说明的是，以上实施例仅用以说明本发明的技术方案而非限制。尽管参照实施例对本发明进行了详细说明，本领域的普通技术人员应当理解，对本发明的技术方案进行修改或者等同替换，都不脱离本发明技术方案的精神和范围，其均应涵盖在本发明的权利要求范围当中。