拖曳线列阵声纳子阵误差失配估计方法技术

本发明专利技术公开的一种拖曳线列阵声纳子阵误差失配估计方法，提供一种能够减小定位误差，并能得到准确的方位估计值，更高角度分辨率的估计方法，本发明专利技术通过下述技术方案予以实现：在阵列流形矩阵模型中，将真实的全阵列流形矩阵表示为每个子阵的位移误差量和每个子阵误差量对全阵列流形矩阵贡献的线性组合；将子阵间的位置误差引入测向模型，对含有子阵间位移失配的全阵列模型运用贝叶斯法则建立数据融合模型，对子阵位置误差向量和目标真实方位进行同时求解，根据融合传感器节点采集到的数据使用贝叶斯算法同时对子阵位移误差和波达方向进行估计，得到多快拍观测值的似然函数；用后验函数得到均方根误差随信噪比变化的方位估计值和位移误差估计值。

【技术实现步骤摘要】
拖曳线列阵声纳子阵误差失配估计方法
本专利技术涉及一种将水听器镶嵌在电缆上形成线列阵，由拖曳电缆拖在舰尾后水中探测目标的声纳。具体涉及一种声纳拖曳线列阵进行被动合成孔径时存在子阵位移误差情况下，同时对子阵位阵列移误差和波达方向估计方法。
技术介绍
拖曳声纳(towedsonar)是将换能器基阵拖曳在运载平台尾后水中探侧目标的声纳。拖线阵声纳基阵是柔性的。拖曳过程中由于舰艇机动和海流影响，以及自身抖动共振，其阵形难以保持稳定，阵型畸变将使得拖线阵声纳难以达到理论性能，这个问题尤其在采用自适应信号处理，空间谱估计等现代阵处理方法中更为严重。因为常规阵处理方法只对阵元能量进行累加，而现代阵处理方法还对阵元信号相关性，信号协方差矩阵特征值等进行计算，虽然现代阵处理方法可以大大提高精度和目标分辨能力，但对阵型稳定更为严格。此外，拖线阵没有垂直孔径，不能区分水面水下目标，在航运密集海域或编队作战时大量目标对其性能产生严重干扰，会产生到处是目标的“满天星”现象。低频主动拖曳线列阵声呐是目前探测安静型潜艇最有效的手段。由于拖曳线列阵声呐基阵远离本舰,且为柔性阵,探测的方位存在较大误差。现代拖曳线列阵声纳系统逐渐趋向于工作在低频段，探测距离不断提高，向着更高的作用距离和探测精度发展。在这样环境下，若要获得较好的空间分辨能力，需要使用孔径较大的阵列，这在实际中，通常意味着更大的系统复杂度和更高的设备成本。在不改变拖曳阵列参数的情况下，利用拖曳线列阵进行被动合成孔径(PSA)，增加了有效孔径的可能性。PSA技术是利用信号时间和空间上的相关性，以及阵列的运动信息、位置信息，构造一个含有多个虚拟子阵的阵列来获得更高的方位分辨力。目前，常见的运动直线阵PSA方法有：Yen、Carey被动合成孔径方法、基于快速傅立叶变换的被动合成孔径，以及扩展拖曳阵测量方法等。虽然具有多个虚拟子阵的PSA已经被广泛地应用到了声纳波束形成，以及波达方位(DOA)估计等领域，但在目标数目较多，以及目标舰船与本舰运动参数与真实条件不一致的情况下，性能会急剧下降。具体地，只存在单目标时，PSA方法可以在相对速度出现误差情况下工作，直接从数据中得到真实的相位校正因子，合成阵元位置不重叠时，仍可准确地进行DOA估计；但在多目标时，若相对速度无误差，PSA方法可以工作，若相对速度存在误差，则此时的DOA结果，甚至会差于单阵常规波束形成的DOA结果，阵列扩展不但没有带来好处，反而增加了定位误差。除此之外，还有其它因素影响DOA的估计精度,例如阵元间距(通常取工作波长的一半)、目标来波之间的角度差、信号的相干性等等。在PSA的基础上，通常利用高分辨DOA方法可以得到目标的精确角度估计。许多高分辨方法，例如多信号分类MUSIC(MultipleSignalClassification)多重信号分类方法、旋转不变子空间方法、稀疏空间谱估计等可以在阵列参数精确已知的情况，提高阵列的空间分辨率。MUSIC算法是一种基于矩阵特征空间分解的方法，从几何角度讲，从几何角度，信号处理的观测空间可以分解为两个空间是正交的信号子空间和噪声子空间。MUSIC算法利用这两个互补空间之间的正交特性来估计空间信号的方位。噪声子空间的所有向量被用来构造谱，所有空间方位谱中的峰值位置对应信号的来波方位，其基本思想是对任意阵列输出数据的协方差矩阵进行特征分解，从而得到与信号分类相对应的信号子空间和与信号分量相正交的噪声。MUSIC算法处理任务是设法估计出入射到阵列的空间信号的个数D以及空间信号源的强度及其来波方向。在实际处理中，观测到的阵列输出数据复向量Y得到的数据是有限时间段内的有限次数的样本，也称快拍或快摄，但是原型MUSIC算法要求来波信号是不相干的。以MUSIC为代表的算法存在一个缺点，即对相干信号处理的不理想。在针对相干信号源的一系列处理方案中，比较经典的是空间平滑技术，如空间平滑(SS)和修正的空间平滑(MSS)算法。然而，空间平滑技术是以损失阵列有效孔径为代价的，而且只适用于等距均匀线阵(ULA)。当存在阵元误差时，假定的杂波模型与实际接收数据不匹配，在模型参数失配的情况下，原本属于信号子空间的目标信号会被错误划分到噪声子空间内，破坏MUSIC算法目标函数的正交性，导致定位失效和。同理，其余各类高分辨DOA算法也对阵列流形的误差很敏感，且在低信噪比、少快拍数、声源相干情况出现时，高分辨算法的性能将会大为退化。由于重构杂波协方差矩阵依赖于杂波模型的准确性，若在假定的杂波模型与实际的数据不匹配时，将导致算法性能下降甚至失效。因此在拖曳线列阵PSA应用中，阵元误差，信号相干、阵元互耦、通道失配等均会使理想模型不再成立。如何更有效地利用接收数据来估计实际的真实阵列流行向量，并进行准确方位估计，成为一个十分重要的课题。目前，国内外在解决这类问题中已进行了一系列的探索。这些文献根据各自应用领域中存在的实际误差来源，从不同侧面对阵列误差进行校正。例如Weiss和Friedlander提出了一种基于最大似然估计阵列自校正的方法，该方法对目标方位和位移失配量进行分块交替优化迭代。由于此类方法是对所有阵元进行建模，在应用到长基线、大规模、长时间PSA中时，会导致维数急剧增多，目标函数收敛极为缓慢、计算效率急剧下降。因此，需要发展针对长基线、大规模、长时间情况拖曳线列阵PSA的阵列流行自校正和方位估计方法。拖曳线列阵声呐也称“拖曳阵声呐”(可简称拖线阵)。它是将水听器镶嵌在电缆上形成线列阵，由拖曳电缆拖在舰艇尾后水中探测目标的声呐。主要用于听测潜艇辐射噪声，进行远程监视、测向和识别，有的也可用于测距。由线列基阵、拖曳电缆、收放装置和绞盘、电子机柜等组成。拖曳线列阵又由前导段、仪器段、基阵段、后导段和尾段构成，阵长数十米至数百米，工作深度可变。具有基阵尺寸大、工作频率低、利于线谱检测，能远距离隐蔽地发现目标等优点；但对拖带舰艇的旋回和倒车等机动有不利影响。声源、海洋信道与水听器阵是水声学研究中的三个基本要素。声源在水中辐射声信号，是形成声场的源头；海洋信道则决定了声波在海洋中的传播特点；水听器阵用以接收声信号，对水中的声场分布进行采样。这三者之间相互密切联系，构成不可分割的统一整体。已知其中两者，就可以推断第三者，这就是水声匹配场处理的基本依据。如果已知水听器阵接收信号和海洋信道信息，待求解的是包括声源位置在内的声源信息，这就是匹配场被动定位。如果已知的是水听器阵接收信号和包括声源位置在内的声源信息，待求解的是海洋信道信息，这就是匹配场反演(MFI：MatchedFieldInversion)。它们都是水声匹配场处理研究的重要内容。匹配场处理技术近年来在水下目标检测、被动定位、海洋环境参数反演等方面的应用受到广泛的关注。水声阵列信号的匹配场匹配场处理(MFP：MatchedFieldProcessing)是利用海洋环境参数和声传播信道特性，通过水下声场模型计算得到接收基阵的声场幅度和相位，形成拷贝场向量，并与基阵接收数据进行匹配，从而实现水下目标的被动定位和海洋环境参数的精确估计。基于远近场声传播特性的拖线阵声纳平台，远场目标信号的平面波传播特性将匹配场定位技术和平面波目标方位估计技术结合，在本文档来自技高网...

【技术保护点】
1.一种拖曳线列阵声纳子阵误差失配估计方法，其特征在于包括如下步骤：在存在子阵间位移失配的阵列流形矩阵模型中，将真实的全阵列流形矩阵表示为每个子阵的位移误差量和每个子阵误差量对全阵列流形矩阵贡献的线性组合；将子阵间的位置误差引入测向模型，对含有子阵间位移失配的全阵列模型运用贝叶斯法则建立数据融合模型，对子阵位置误差向量β和目标真实方位角α‑1进行同时求解，根据融合传感器节点采集到的数据，使用贝叶斯算法同时对子阵位移误差和波达方向进行估计，计算得到多快拍观测值的似然函数；用后验函数得到均方根误差随信噪比变化的方位估计值和位移误差估计值。

【技术特征摘要】
1.一种拖曳线列阵声纳子阵误差失配估计方法，其特征在于包括如下步骤：在存在子阵间位移失配的阵列流形矩阵模型中，将真实的全阵列流形矩阵表示为每个子阵的位移误差量和每个子阵误差量对全阵列流形矩阵贡献的线性组合；将子阵间的位置误差引入测向模型，对含有子阵间位移失配的全阵列模型运用贝叶斯法则建立数据融合模型，对子阵位置误差向量β和目标真实方位角α-1进行同时求解，根据融合传感器节点采集到的数据，使用贝叶斯算法同时对子阵位移误差和波达方向进行估计，计算得到多快拍观测值的似然函数；用后验函数得到均方根误差随信噪比变化的方位估计值和位移误差估计值。2.如权利要求1所述的拖曳线列阵声纳子阵误差失配估计方法，其特征在于：失配后的真实子阵间阵列流行向量在真实子阵相对位置向量的预设子阵相对位置向量r处进行一阶泰勒展开逼近，得到子阵间近似阵列流行向量式中，向量h为子阵间近似阵列流行向量，e为欧拉常数，j为虚数单位常数，k为波数，上标T表示转置，P为子阵数，θ为表示扫描来波的方位角，向量为真实的子阵相对位置向量，为第p个子阵的真实相对位置，向量r＝[r1,...,rP]T为预设的子阵相对位置向量，rp为第p个子阵的预设相对位置，β＝[β1,...,βP]T为子阵位置误差向量，为第p个子阵位移误差，diag(β)表示以向量β的元素作为对角线元素的对角矩阵。3.如权利要求2所述的拖曳线列阵声纳子阵误差失配估计方法，其特征在于：根据一阶泰勒展开逼近得到的子阵间近似阵列流行向量，将真实的全阵阵列流形向量近似为子阵间阵列流行矩阵式中，是来波方位为θ时第p列为vp(θ)其余列为零向量的矩阵，向量vp为子阵间阵列流行矩阵V(θ)的第p列，ap(θ)为第p个子阵的阵列流行向量。4.如权利要求3所述的拖曳线列阵声纳子阵误差失配估计方法，其特征在于：在全阵列模型中，将全阵的阵列流形向量沿N个(θ1到θN)的来波方位进行扫描，组合成全阵的近似阵列流行矩阵A1，再将每个子阵的位移误差量和每个子阵误差量对全阵列流形矩阵的贡献的线性组合为:且无子阵间位置误差时预设的全阵阵列流形矩阵Aw＝[aw(θ1),...,aw(θN)]，第p个子阵误差投影矩阵Bp＝-jk[Vp(θ1)h(θ1,r),...,Vp(θN)h(θN,r)]，式中，N为扫描方位网格数量，aw为预设的全阵阵列流形向量，βpBp是第p个子阵位置误差引起的全阵阵列流形矩阵误差的一阶逼近乘积。第p个子阵误差投影矩阵Bp＝-jk[Vp(θ1)h(θ1,r),...,Vp(θN)h(θN,r)]，式中，N为扫描方位网格数量，aw为预设的全阵阵列流形向量，Aw＝[aw(θ1),...,aw(θN)]是无子阵间位置误差时预设的全阵阵列流形矩阵，βpBp是第p个子阵位置误差引起的全阵阵列流形矩阵误差的一阶逼近乘积。5.如权利要求1所述的拖曳线列阵声纳子阵误差失配估计方法，其特征在于：在基于子阵间位移失配阵列模型的贝叶斯定位中，使用贝叶斯算法对全阵列模型子阵间含子阵位置误差向量β和目标真实方位α-1进行同时求解，得到表示近似的全阵阵列流行矩阵和近似的全阵的阵列模型全阵接收信号向量x(t)＝Φ(β)s(t)+e(t)，根据全阵接收信号向量x(t)＝[x1(t),x2(t),...,xN(...

【专利技术属性】
技术研发人员：乔文昇，王立，雷志雄，李维科，李明兵，
申请(专利权)人：西南电子技术研究所中国电子科技集团公司第十研究所，
类型：发明
国别省市：四川,51