一种基于畸变拖曳阵的保真增强波束形成的方法技术

本发明专利技术公开了一种基于畸变拖曳阵的保真增强波束形成的方法，包括如下步骤：(1)模拟水声目标辐射噪声s(t)；(2)模拟观测阵列信号xi(t)，i＝1,2,...,M，M为拖曳阵中阵元数目；(3)基于理想波束形成粗略估计目标方位

A method to enhance the fidelity of towed array beamforming based on distortion

The invention discloses a method for enhancing the fidelity of beamforming distortion based on towed array, which comprises the following steps: (1) simulation of underwater target radiated noise s (T); (2) the simulation observation array signal Xi (T), I = 1,2, M,..., M is the number of elements in the array (3; a rough estimate of target azimuth) form an ideal beam based on

【技术实现步骤摘要】
一种基于畸变拖曳阵的保真增强波束形成的方法
本专利技术属于声纳信号处理领域，具体涉及一种基于畸变拖曳阵的保真增强波束形成的方法。
技术介绍
当前的声纳系统主要分为水听器岸基阵声纳和水听器拖曳线列阵声纳。水听器拖曳线列阵声纳简称拖曳阵，是拖曳在舰船尾部一定距离的声探测系统，通过接收航行目标自身的辐射噪声或反射的主动信号，来检测目标的有无并估计目标的有关参数。它具有探测能力强，探测频率低，水文适应能力强和无盲区的特点。水听器阵列流形是水听器阵列的一个重要参数，利用水听器阵列接收信号作波束形成等声纳信号处理时，一般都要求阵形己知。然而水听器阵列布放到水中之后，其阵形不好控制，一般认为当阵形畸变超过λ/10(λ为接收信号波长)时，就应在波束形成时对阵形畸变进行补偿，否则会对拖曳线列阵声呐性能产生重大影响，从而影响方位估计性能。现有的拖曳阵阵形估计方法主要可以分为两类：一类是声学计算的方法，它是利用水听器阵列的接收信号来反推出阵列流形；另一类是非声辅助测量的方法，主要是在水听器拖曳阵上安装若干个深度或航向传感器，利用这些传感器的测量值来估计出阵形，实现畸变线列阵的保真波束形成。声学计算的方法主要有锐度提取法和特征矢量法。锐度提取法在阵元较多的情况下搜索算法过于复杂，随后很少有人研究，之后该方法被应用于阵形未知的自适应波束形成，但是相比于该方法，特征矢量法性能较好。该方法只需要一个声源就可以估计出阵形。但一般要求声源方位精确已知，此外，它还要求首阵元的坐标已知，阵列形变不是很大，相邻阵元间距己知且固定。除了这两种方法外，还可以用最大似然估计方法和隐藏马尔科夫模型的方法估计畸变，实现畸变线列阵的较精确的波束形成。非声辅助测量方法也主要有两种，流体力学法和插值拟合法。流体力学法将辅助传感器的测量值作为边界条件，通过求解拖缆的流体力学方程来估计畸变。这种方法为了得到更加精确的估计结果，往往在阵列尾部安装另一个航向传感器以做修正。另外，在方程中，阵列所受到的牵引力只来自阵列前端的拖点处。然而在实际情况中，当海况不好时，阵列会受到海浪涌流等因素的影响，这种方法估计出来的波束形成结果不太可信。插值拟合法是用辅助传感器测量阵列上的多点状态，然后用样条插值拟合出畸变的阵形实现精确的波束形成。这种方法假设各航向传感器的横坐标一直保持不变，即仍然等于阵列没有形变时的横坐标，因此估计出来的阵形与实际的阵形会有差别，阵列的总长度会拉长。当阵形变化较小时，这种差别影响不大，但当阵形变化较大时就需要对估计出的阵形做修正，从而使波束形成的结果更接近于真实值，提高阵列的方位估计性能。自从1917年第一套拖曳阵声纳系统被专利技术以来，拖曳阵的发展已经经历了100年的时间。如今，拖曳阵在军事目标定向及海洋地层、石油勘探等商业领域都得到广泛地应用。其中，大多数应用都是基于阵列流形己知的前提，然而，由于海洋环境的复杂多变，水听器阵列在拖曳过程中，海浪、涌流、拖船机动都会改变阵列的形状，即阵列流形在拖曳过程中是不断变化的。常用的解决办法是假设拖曳阵始终保持为一条直线阵，但是在军事安全和商业利益的双重推动下，拖曳阵不断地向多基元、长距离、多维化的大规模方向发展，导致阵列在拖曳过程中阵形越来越难控制，这种假设阵形不变的做法不能满足实际应用需求。关于畸变的拖曳线列阵的目标方位估计准确度问题日显突出。如何在拖曳线列阵发生畸变的情况下，更为准确的进行目标信号的方位估计将成为提高阵列探测性能、促使拖曳阵走向高精度、高分辨率应用方向的关键。目前已有多种针对畸变的拖曳线列阵进行估计的方法，一种最常用的途径是在拖线阵上放置多个水平深度测量仪，以获得阵元的位置信息，进而估计阵形，实现较精确的波束形成和目标信号方位估计。这种方法应用起来比较直接，但是经济代价太高。另外一种方法是通过估计误差参数来实现对阵形的校正，该类方法首先对阵列误差进行建模，将阵列误差校正转化为参数估计的问题。该类阵形校正方法通常可以分为有源校正类和自校正类。对于有源校正方法，该方法对辅助信号源有着较高的精确方位信息的要求，所以当辅助信号源的方位信息有偏差时，这类校正算法会带来较大的偏差。自校正算法由于阵元位置误差与方位参数之间的耦合和某些病态的阵列结构，参数估计的唯一辨识往往无法保证，更为重要的是参数联合估计对应的高维、多模非线性优化问题带来了庞大的运算量，估计的全局收敛性往往无法保证。
技术实现思路
专利技术目的：针对现有技术中存在的问题，本专利技术公开了一种基于畸变拖曳阵的保真增强波束形成的方法，该方法通过时延估计来校正拖曳线列阵的畸变对波束形成影响，获得保真增强的目标辐射噪声跟踪波束。技术方案：本专利技术采用如下技术方案：一种基于畸变拖曳阵的保真增强波束形成的方法，包括如下步骤：(1)模拟水声目标辐射噪声s(t)；(2)模拟观测阵列信号xi(t)，i＝1,2,...,M，M为拖曳阵中阵元数目；(3)基于理想波束形成粗略估计目标方位为目标信号波束能量最大时的引导角；(4)检测目标信号的L个功率最大的线谱位置l＝1,2,...,L；(5)由强线谱位置处的相位差估计拖曳阵中每个阵元的平均时延差Δτi，i＝1,2,...,M；(6)基于估计时延获取保真增强的目标跟踪波束具体地，水声目标辐射噪声s(t)包括平稳连续谱分量sc(t)和线谱分量sl(t)；所述平稳连续谱分量sc(t)的获取步骤如下：(A.1)采用三参数模型法模拟平稳连续谱的功率谱Gxf(ωt)：其中ωm，ωc和λ为三参数模型的三个参数，决定了该连续谱的形状；ωt为频率，ωm为尖锐度因子，决定谱锋的尖锐程度和高度，ωc决定谱锋的位置，λ决定功率谱高、低频端幅度的相对比例，σ表示平稳连续谱信号的能量；(A.2)建立p阶AR滤波器，其Yule-Walker方程为：其中a[l]，l∈{1,2,…,p}和b0为p阶AR滤波器系数，δ[k]为冲击函数；rx[k]为Gxf(ωt)的自相关函数rc(τ)的采样值；(A.3)采用Levison-Durbin算法求解式(2)方程，得到p阶AR滤波器系数；高斯白噪声通过该AR滤波器后得到的信号，即为水声目标辐射噪声中的平稳连续谱分量sc(t)；所述线谱分量sl(t)的获取步骤如下：(B.1)采用K个正弦信号来模拟目标信号的线谱分量，其中Ak为正弦信号幅度，fk为正弦信号的频率，t∈[0,T]为观测时间；(B.2)在线谱位置fk处计算平稳连续谱分量sc(t)的能量PIk，k＝1,2,...,K；(B.3)根据已知的信号干扰比计算出各个正弦信号幅度Ak，即得到水声目标辐射噪声中的线谱分量sl(t)。具体地，步骤(2)包括如下步骤：(2.1)将拖曳阵中第一个阵元设置为参考阵元，其阵元数据为：s1(t)＝s(t)；(2.2)拖曳阵中其余M-1个阵元的阵元数据为：si(t)＝s[t-timeDelay(i)]，i＝2,…,M；其中timeDelay(i)为第i个阵元相对于参考阵元的时延：tarDis(i)为声源与第i个阵元之间的距离，v是声音在水中的传播速度。(2.3)根据已知的信噪比计算能量Pn，并生成能量为Pn的M路高斯白噪声ni(t)，其中i＝1,…,M，sl(t)为水声目标辐射噪声线谱分量；(2.4)观测阵列信号xi(t)为：xi本文档来自技高网...

【技术保护点】
一种基于畸变拖曳阵的保真增强波束形成的方法，其特征在于，包括如下步骤：(1)模拟水声目标辐射噪声s(t)；(2)模拟观测阵列信号xi(t)，i＝1,2,...,M，M为拖曳阵中阵元数目；(3)基于理想波束形成粗略估计目标方位

【技术特征摘要】
1.一种基于畸变拖曳阵的保真增强波束形成的方法，其特征在于，包括如下步骤：(1)模拟水声目标辐射噪声s(t)；(2)模拟观测阵列信号xi(t)，i＝1,2,...,M，M为拖曳阵中阵元数目；(3)基于理想波束形成粗略估计目标方位为目标信号波束能量最大时的引导角；(4)检测目标信号的L个功率最大的线谱位置(5)由强线谱位置处的相位差估计拖曳阵中每个阵元的平均时延差△τi，i＝1,2,...,M；(6)基于估计时延获取保真增强的目标跟踪波束2.根据权利要求1所述的基于畸变拖曳阵的保真增强波束形成的方法，其特征在于，水声目标辐射噪声s(t)包括平稳连续谱分量sc(t)和线谱分量sl(t)；所述平稳连续谱分量sc(t)的获取步骤如下：(A.1)采用三参数模型法模拟平稳连续谱的功率谱Gxf(ωt)：其中ωm，ωc和λ为三参数模型的三个参数，决定了该连续谱的形状；ωt为频率，ωm为尖锐度因子，决定谱锋的尖锐程度和高度，ωc决定谱锋的位置，λ决定功率谱高、低频端幅度的相对比例，σ表示平稳连续谱信号的能量；(A.2)建立p阶AR滤波器，其Yule-Walker方程为：其中a[l]，l∈{1,2,…,p}和b0为p阶AR滤波器系数，δ[k]为冲击函数；rx[k]为Gxf(ωt)的自相关函数rc(τ)的采样值；(A.3)采用Levison-Durbin算法求解式(2)方程，得到p阶AR滤波器系数；高斯白噪声通过该AR滤波器后得到的信号，即为水声目标辐射噪声中的平稳连续谱分量sc(t)；所述线谱分量sl(t)的获取步骤如下：(B.1)采用K个正弦信号来模拟目标信号的线谱分量，其中Ak为正弦信号幅度，fk为正弦信号的频率，t∈[0,T]为观测时间；(B.2)在线谱位置fk处计算平稳连续谱分量sc(t)的能量PIk，k＝1,2,...,K；(B.3)根据已知的信号干扰比计算出各个正弦信号幅度Ak，即得到水声目标辐射噪声中的线谱分量sl(t)。3.根据权利要求1所述的基于畸变拖曳阵的保真增强波束形成的方法，其特征在于，步骤(2)包括如下步骤：(2.1)将拖曳阵中第一个阵元设置为参考阵元，其阵元数据为：s1(t)＝s(t)；(2.2)拖曳阵中其余M-1个阵元的阵元数据为：si(t)＝s[t-t...

【专利技术属性】
技术研发人员：武其松，徐萍，李腾飞，罗昕炜，方世良，
申请(专利权)人：东南大学，
类型：发明
国别省市：江苏,32