CN113687296A - 一种基于阵列信号数值计算的水下目标辐射噪声相干信源测向方法 - Google Patents

一种基于阵列信号数值计算的水下目标辐射噪声相干信源测向方法 Download PDF

Info

Publication number
CN113687296A
CN113687296A CN202110780405.1A CN202110780405A CN113687296A CN 113687296 A CN113687296 A CN 113687296A CN 202110780405 A CN202110780405 A CN 202110780405A CN 113687296 A CN113687296 A CN 113687296A
Authority
CN
China
Prior art keywords
arrival
coherent
radiation noise
array
matrix
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202110780405.1A
Other languages
English (en)
Other versions
CN113687296B (zh
Inventor
戴卫华
戴高千子
王彬
李天昀
孟玉婷
于沛东
李勇斌
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Information Engineering University of PLA Strategic Support Force
Original Assignee
Information Engineering University of PLA Strategic Support Force
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Information Engineering University of PLA Strategic Support Force filed Critical Information Engineering University of PLA Strategic Support Force
Priority to CN202110780405.1A priority Critical patent/CN113687296B/zh
Publication of CN113687296A publication Critical patent/CN113687296A/zh
Application granted granted Critical
Publication of CN113687296B publication Critical patent/CN113687296B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S3/00—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
    • G01S3/02—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using radio waves
    • G01S3/14—Systems for determining direction or deviation from predetermined direction

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)

Abstract

本发明公开一种基于阵列信号数值计算的水下目标辐射噪声相干信源测向方法，首先通过非实时计算，预先存储线谱分量在不同相干信源波达方向组合下的阵列信号的生成矩阵和转换矩阵；然后通过盖氏圆信源估计方法和常规波束形成方法，分别实时估计水下目标辐射噪声中的相干信源数目和波达方向分布范围，确定相干信源波达方向的备选组合；最后依据波达方向备选组合，搜索预先存储的生成矩阵和转换矩阵，实时求解相干信源与阵列信号模型误差，将模型误差最小时的波达方向组合确定为最优解，并将最优模型系数作为相干信源强度。本发明能够克服多途影响，实时并精确计算水下目标辐射噪声的相干信源波达方向，对阵元安装误差和噪声观测误差具有良好包容性。

Description

一种基于阵列信号数值计算的水下目标辐射噪声相干信源测 向方法
技术领域
本发明属于水下目标辐射噪声测向定位技术领域，尤其涉及一种基于阵列信号数值计算的水下目标辐射噪声相干信源测向方法，适用于多途信道中相干信源的测向。
背景技术
对水下目标辐射噪声进行测向定位具有重要意义。但由于水声信道的不均匀波导与有界空间结构特性，目标辐射噪声在海水中容易发生多途传播现象，从而形成多个相干信源。在众多测向定位算法中，常规波束形成(Conditional Beamforming,CBF)算法具有最好的稳健性，即对阵列安装误差和噪声特性变化包容性最好。但CBF算法受限于“傅里叶限”限制，无法分辨位于同一个波束宽度内的多个相干信源。若要提高CBF算法的方位分辨能力，需要扩大阵列有效孔径或阵元数目，但在水下机动平台上较难实施。其它测向算法，例如最小方差无畸变响应波束形成算法(Minimum Variance Distortionless Response，MVDR)、多重信号分类方法(MUltiple SIgnal Characteristics，MUSIC)和旋转不变子空间技术(Estimation of Signal Parameters via Rotational Invariance Techniques，ESPRIT)等虽然具有更高的方位分辨能力，但它们普遍要求海洋背景噪声与目标辐射噪声不相关、多个目标辐射噪声之间不相关，并且要求背景噪声严格服从某些确定的分布，以及依赖阵列高精度安装和标定。在多途水声信道中，阵列接收的目标辐射噪声可能由多个相干信源构成，而且海洋背景噪声统计特性也十分复杂，往往受到洋流、航运和气象等时变空变因素影响。因此高分辨测向算法在实际应用中面临较大的局限性，最直接表现就是稳健性不足，远不如CBF方法。为了对多个相干信源进行测向，在这些高分辨测向算法的基础上产生了一些改进算法，即在测向之前采用空间平滑等方法做解相干预处理，但预处理会损失阵列孔径和方位分辨能力。除此之外，还有一些折中算法在方位分辨率和稳健性方面取得一定的平衡，例如波束域MVDR算法、对角加载算法和基于最差性能优化的稳健波束形成算法等。这些折中方法的方位分辨率和稳健性均介于CBF和其它高分辨算法之间。
与上述算法相比，以极大似然(Maximum Likelihood，ML)算法和加权子空间拟合(Weight Subspace Fitting，WSF)为典型代表的空间拟合类算法性能更加突出：通过构造拟合关系得到一个需要进行多维求解的目标函数，目标方位估计精度高，尤其在低信噪比和小快拍数据的情况下估计性能优良，而且不需要进行解相干处理就可以对相关信号进行DOA估计(安春莲.独立信号与相干信号并存的测向算法研究[D].哈尔滨工程大学.2013)。从性能上，空间拟合类算法十分适合多途条件下的水下目标辐射噪声的测向定位。但这类算法的不足在于目标函数求解过程涉及多维搜索，计算量繁重，很难满足工程应用中实时计算要求。
发明内容
本发明针对空间拟合类算法存在的目标函数求解过程涉及多维搜索，计算量繁重，很难满足工程应用中实时计算要求的问题，提出一种基于阵列信号数值计算的水下目标辐射噪声相干信源测向方法。
为了实现上述目的，本发明采用以下技术方案：
一种基于阵列信号数值计算的水下目标辐射噪声相干信源测向方法，包括：
步骤1：生成矩阵和转换矩阵的非实时计算存储：预先将水下目标辐射噪声不同频率线谱的到达方向做数值化离散，然后依据不同信源数目分别做组合排列，最后计算每个组合下的生成矩阵和转换矩阵；
步骤2：目标辐射噪声相干信源数目的实时估计：实时计算各个阵元通道采集的目标辐射噪声中线谱分量的样本协方差矩阵，然后通过盖氏圆信源估计方法估计目标辐射噪声相干信源数目；
步骤3：目标辐射噪声相干信源到达方向分布范围的实时估计：计算各阵元通道采集的目标辐射噪声的线谱分量的DFT值，然后通过常规波束形成方法估计目标辐射噪声相干信源到达方向分布范围；
步骤4：目标辐射噪声相干信源到达方向的实时精确计算：依据目标辐射噪声相干信源数目估计结果和到达方向分布范围，确定目标辐射噪声相干信源到达方向的所有可能组合以及阵列信号数学模型；然后选取预先存储的生成矩阵、转换矩阵与频域阵列信号向量相乘，计算模型系数与模型偏差；最后将模型偏差最小时对应的到达方向组合确定为目标辐射噪声的相干信源真实到达方向，并将此时的模型系数作为相干信源强度。
进一步地，所述步骤1包括：
将水下目标辐射噪声的不同频率线谱到达方向做数值化离散，离散间隔为Δθ分辨率，离散波达方向的数目为M；
相干辐射噪声到达方向依据信源组合从1至L-3做组合排列，组合个数分别为
Figure BDA0003156566660000034
其中L为阵元数目；
对于每个频率的每个组合，计算转换矩阵：
D(fi，θ1，…，θm)＝[GH(fi，θ1，…，θm)·G(fi，θ1，…，θm)]-1·GH(fi，θ1，…，θm)
其中，1＜m≤M，G(fi，θ1，…，θm)表示在线谱频率fi以及波达方向组合(θ1，…，θm)条件下，阵列信号模型中的生成矩阵；阵列信号模型为：
Figure BDA0003156566660000031
记
Figure BDA0003156566660000032
S(fi)＝(S(fi，θ1) S2(fi，θ2)...Sm(fi，θm))T
X(fi)＝(X1(fi) X2(fi)...XL(fi))T
其中，S(fi)为待求的相干信源DFT值；X(fi)为阵列实际接收的目标辐射噪声的线谱DFT值；G(fi，θ1，…，θm)的第一行元素取值均为1，表示以第一个阵元为参考阵元，其它行元素表示其它阵元接收信号线谱DFT值相对参考阵元接收信号线谱DFT值的增量；G(fi，θ1，θ2，…，θm)的第l行第k列元素取值为：
Figure BDA0003156566660000033
其中，d为阵元间距，满足d≤λi/2，λi为线谱频率fi对应的载波波长；
将生成矩阵G(fi，θ1，θ2，…，θm)和转换矩阵D(fi，θ1，…，θm)进行非实时计算，并存储；
每个频率的线谱需要预先存储的生成矩阵、转换矩阵的元素个数总数均为：
Figure BDA0003156566660000041
其中，L为阵元数目。
进一步地，所述步骤2包括：
首先，对水下目标辐射噪声进行观测接收和实时采集，做DFT变换，实时计算L个阵元通道的基频线谱分量DFT值的样本协方差矩阵：
C＝X(fi)H·X(fi)
其中X(fi)＝(X1(fi)，X2(fi)，…，XL(fi))，为阵列实际接收的目标辐射噪声的线谱DFT值；L为阵元数目；fi为线谱频率；
然后，根据盖氏圆信源估计方法，由C的前(L-1)行和前(L-1)列构成新矩阵C1，对C1做特征值分解，计算C1的特征值和特征向量；将所述特征值从大到小排序，构造对角矩阵D1＝diag(λ′1，λ′2，…，λ′K，…，λ′L-1)，其中K是待估计的相干辐射噪声信源个数，λ′1、λ′2…λ′K…λ′L-11分别为从大到小排序的C1的特征值；根据所述特征向量构造酉矩阵U1＝[q′1，q′2，…，q′M，…，q′L-1]，q′1、q′2…q′M…q′L-1分别为λ′1、λ′2…λ′K…λ′L-1对应的特征向量，满足
Figure BDA0003156566660000042
再重构酉矩阵
Figure BDA0003156566660000043
用生成的酉矩阵U对样本协方差矩阵C做酉变换：
Figure BDA0003156566660000044
其中，ρ1、ρ2、ρ3…ρL-1、
Figure BDA0003156566660000045
cLL为UHCU计算后矩阵中的元素；
则按照下式估计相干信源数目：
Figure BDA0003156566660000046
其中，0≤D(N)≤1，是关于观测时间内共采集的噪声点数N的递减函数；GDE(k)是用来帮助判断信源数目的一个过程量，当k＝k0时，GDE(k0)为第一个比零小的值，则信源数的估计值取K＝k0-1。
进一步地，所述步骤3包括：
首先，依据X(fi)＝(X1(fi)，X2(fi)，…，XL(fi))做常规波束形成：
d(θm)＝a(θm)·XH(k)，m＝1，2，…，M.
其中，θm为波达方向离散后第m个波达方向，a(θm)为该波达方向的阵列导向矢量；
基于d(θm)计算功率，总共计算得到M个E(θm)：
E(θm)＝d*(θm)·d(θm)，m＝1，2，…，M.
然后，再根据步骤2估计得到的相干信源数目K，从M个E(θm)中搜索峰值和次峰值，将峰值和次峰值所在方位记作
Figure BDA0003156566660000051
当κ＝K时，以
Figure BDA0003156566660000052
的每个方位值为中心，取峰值或次峰值-3dB宽度区间以内的离散方位值集合作为相干辐射噪声到达方向分布范围，记作
Figure BDA0003156566660000053
当κ＜K时，以
Figure BDA0003156566660000054
的每个方位值为中心，取峰值或次峰值-6dB宽度区间以内的离散方位值集合作为相干辐射噪声到达方向分布范围，记作
Figure BDA0003156566660000055
进一步地，所述步骤4包括：
首先，依据步骤2估计的相干信源数目K以及步骤3估计的目标辐射噪声相干信源到达方向
Figure BDA0003156566660000056
确定相干信源到达方向的所有可能组合；记
Figure BDA0003156566660000057
的元素个数分别为m1，m2，…，mκ；当κ＝K时，从κ个集合
Figure BDA0003156566660000058
中各分配一个信源波达方向进行组合，组合数目为每个集合长度的乘积，即
Figure BDA0003156566660000059
当κ＜K时，重新合并集合
Figure BDA00031565666600000510
新集合元素个数为
Figure BDA00031565666600000511
在
Figure BDA00031565666600000512
中做K个相干信源的组合，组合数目为
Figure BDA00031565666600000513
然后，在
Figure BDA00031565666600000514
个波达方向组合中，对每一个相干信源波达方向的组合，记任意一个波达方向组合为
Figure BDA00031565666600000515
查找预先存储的生成矩阵
Figure BDA00031565666600000516
和转换矩阵
Figure BDA00031565666600000517
计算K个相干信源的DFT值：
Figure BDA0003156566660000061
以及模型偏差：
Figure BDA0003156566660000062
最后，对于
Figure BDA0003156566660000063
个波达方向组合，总共有
Figure BDA0003156566660000064
个模型偏差，从中搜索模型偏差最小值，将最小值对应的波达方向组合
Figure BDA0003156566660000065
作为K个相干信源的真实波达方向，并将对应的DFT值
Figure BDA0003156566660000066
中每个元素的幅度平方值作为相干信源强度。
与现有技术相比，本发明具有的有益效果：
本发明通过非实时计算，预先存储阵列信号的生成矩阵和转换矩阵，大幅度降低了线谱分量信号空间的多维搜索求解复杂度；并通过盖氏圆方法和常规波束形成方法分别估计相干信源数目和波达方向分布范围，进一步降低线谱分量信号空间的多维搜索求解范围。本发明充分利用了空间拟合类算法对相干信源波达方位估计精度高的优点，又大幅降低了实时计算复杂度，适用于多途条件下对水下目标辐射噪声的相干信源测向定位。
附图说明
图1为本发明实施例一种基于阵列信号数值计算的水下目标辐射噪声相干信源测向方法的流程图。
具体实施方式
下面结合附图和具体的实施例对本发明做进一步的解释说明：
如图1所示，一种基于阵列信号数值计算的水下目标辐射噪声相干信源测向方法，具体步骤如下：
S1：生成矩阵和转换矩阵的非实时计算存储：预先将水下目标辐射噪声不同频率线谱的到达方向做数值化离散，然后依据不同信源(多途)数目分别做组合排列，最后计算每个组合下的生成矩阵和转换矩阵；
S2：目标辐射噪声相干信源数目的实时估计：实时计算各个阵元通道采集的目标辐射噪声中线谱分量的样本协方差矩阵，然后通过盖氏圆信源估计方法估计目标辐射噪声相干信源数目；
S3：目标辐射噪声相干信源到达方向分布范围的实时估计：计算各阵元通道采集的目标辐射噪声的线谱分量的DFT值，然后通过常规波束形成方法估计目标辐射噪声相干信源到达方向分布范围；
S4：目标辐射噪声相干信源到达方向的实时精确计算：依据目标辐射噪声相干信源数目估计结果和到达方向分布范围，确定目标辐射噪声相干信源到达方向的所有可能组合以及阵列信号数学模型；然后选取预先存储的生成矩阵、转换矩阵与频域阵列信号向量相乘，计算模型系数与模型偏差；最后将模型偏差最小时对应的到达方向组合确定为目标辐射噪声的相干信源真实到达方向，并将此时的模型系数作为相干信源强度。
具体地，所述S1包括：
首先，将水下目标辐射噪声的不同频率线谱到达方向做数值化离散。当前，水下目标普遍采用螺旋桨推动的方式进行机动，在定速巡航状态下，目标的辐射噪声中一般包含了与螺旋桨转速成整倍数关系的多条线谱分量，线谱的强度往往与水下目标的排水量和航速成正比，因此可对目标辐射噪声的线谱分量进行观测。例如，对水下目标辐射噪声中的基频线谱分量进行观测。基频大小与水下目标的航速成正比，即由螺旋桨转速决定。当目标处于定速潜航状态时，该基频取值可以稳定分布在某个区间，记作fbasee∈(fa，fb)，按照一定分辨率Δf将其离散化，记作f∈{f1，f2，…，fγ}。进一步，对水下目标辐射噪声的波达方向进行离散。由于阵列接收具有指向性，能够观测的目标辐射噪声波达方向的范围有限，记作θ∈(θa，θb)，按照一定分辨率Δθ将其离散化，记作θ∈{θ1，θ2，…，θM}。
然后，对于每一个线谱频率f，按照信源数目从1增加至(L-3)，对其波达方向进行组合排列，组合个数分别为
Figure BDA0003156566660000071
其中L为阵元数目，M为离散波达方向的数目。以频率f1的线谱为例，在1个信源的情况下，波达方向组合为(θ1)，(θ2)，…，(θM)，共
Figure BDA0003156566660000072
个。在2个信源的情况下，波达方向组合为：
(θ1，θ2)，(θ1，θ3)，…，(θ1，θM)；(θ2，θ3)，(θ2，θ4)，…，(θ2，θM)；…；(θM-1，θM).共
Figure BDA0003156566660000073
最后，对于每个频率的每个组合，计算转换矩阵：
D(fi，θ1，…，θm)＝[GH(fi，θ1，…，θm)·G(fi，θ1，…，θm)]-1·GH(fi，θ1，…，θm)
其中，1＜m≤M，G(fi，θ1，…，θm)表示在线谱频率fi以及波达方向组合(θ1，…，θm)条件下，阵列信号模型中的生成矩阵。阵列信号模型为：
Figure BDA0003156566660000081
记
Figure BDA0003156566660000082
S(fi)＝(S(fi，θ1) S2(fi，θ2)…Sm(fi，θm))T
x(fi)＝(X1(fi) X2(fi)...XL(fi))T
其中，S(fi)为待求的相干信源DFT值，x(fi)为阵列实际接收的目标辐射噪声的线谱DFT值。G(fi，θ1，…，θm)的第一行元素取值均为1，表示以第一个阵元为参考阵元，其它行元素表示其它阵元接收信号线谱DFT值相对参考阵元接收信号线谱DFT值的增量。对于构型固定和通道一致性经过校准的小规模接收阵列，G(fi，θ1，…，θm)中的其它元素只与阵列构型、线谱频率fi以及波达方向组合(θ1，…，θm)有关，可以预先计算。以均匀线阵为例，阵元数目为L，阵元间距为d，满足d≤λi/2，λi为线谱频率fi对应的载波波长，G(fi，θ1，θ2，…，θm)的第l行第k列元素取值为：
Figure BDA0003156566660000083
因此，φl，k、生成矩阵G(fi，θ1，θ2，…，θm)和转换矩阵D(fi，θ1，…，θm)都可非实时计算，并存储，计算量可忽略不计。每个频率的线谱需要预先存储的生成矩阵、转换矩阵的元素个数总数均分别为：
Figure BDA0003156566660000084
具体地，所述S2包括：
首先，采用构型固定和通道一致性经过校准的小规模接收阵列，对水下目标辐射噪声进行观测接收和实时采集。假设阵元数目为L，采样频率fs，在观测时间内共采集噪声点数为N。第l个阵元通道采集的噪声样本记作xl＝(xl[1]，xl[2]，…，xL[N])，均为复数向量或IQ数据。做DFT变换，计算L个阵元通道的基频线谱分量DFT值的样本协方差矩阵：
C＝X(fi)H·X(fi)
其中X(fi)＝(X1(fi)，X2(fi)，…，XL(fi))，XL(fi)为第L个阵元实际接收的目标辐射噪声的线谱DFT值。
然后，参考盖氏圆信源估计方法[Hsien-Tsai Wu，Jar-Fen Yang and Fwu-KuenChen.Source number estimation using transformed Geschgorin radii[J].IEEETransactions on signal processing，1995.vol.43(6)：1325-1333x]，由C的前(L-1)行和前(L-1)列构成新矩阵C1；对C1做特征值分解，计算C1的特征值和特征向量；将特征值从大到小排序，构造对角矩阵，记作D1＝diag(λ′1，λ′2，…，λ′K，…，λ′L-1)，其中K是待估计的相干辐射噪声信源个数，λ′1、λ′2…λ′K…λ′L-1分别为从大到小排序的C1的特征值；用所述特征向量构造酉矩阵，记作U1＝[q′1，q′2，…，q′M，…，q′L-1]，q′1、q′2…q′M…q′L-1分别为λ′1、λ′2…λ′K…λ′L-1对应的特征向量，则满足
Figure BDA0003156566660000091
再重构酉矩阵
Figure BDA0003156566660000092
用生成的酉矩阵U对样本协方差矩阵C做酉变换：
Figure BDA0003156566660000093
其中，ρ1、ρ2、ρ3…ρL-1、
Figure BDA0003156566660000094
cLL为UHCU计算后矩阵中的元素，
Figure BDA0003156566660000095
分别为ρ1、ρ2、ρ3…ρL-1的共轭；
进一步，计算
Figure BDA0003156566660000096
其中，0≤D(N)≤1，是关于N的递减函数，例如取值D(N)＝1/N。GDE(k)是用来帮助判断信源数目的一个过程量，如果当k＝k0时，GDE(k0)为第一个比零小的值，则信源数的估计值取K＝k0-1。显然K≤L-3，即该盖氏圆方法适用于相干辐射噪声信源数目少于阵元数目3个以上的场合。例如阵元数为9，那么采用该盖氏圆方法可以估计信源数目不超过6个的相干辐射噪声。考虑到海洋中水声信道的多途现象具有随机性，所有多途噪声同时被阵列接收的概率较低，可近似为稀疏多径信道条件，因此该信源估计方法可以满足大部分水下目标测向定位需求。并且该盖氏圆信源估计方法在低信噪比和相干条件下具有良好的估计性能。
具体地，所述S3包括：
首先，依据所述S2计算的基频线谱分量的DFT值向量X(fi)＝(X1(fi)，X2(fi)，…，XL(fi))做常规波束形成：
d(θm)＝a(θm)·XH(k)，m＝1，2，…，M.
其中，θm为波达方向离散后第m个波达方向，a(θm)为该波达方向的阵列导向矢量。进一步计算功率E(θm)，总共计算得到M个E(θm)。
E(θm)＝d*(θm)·d(θm)，m＝1，2，…，M.
其中，d*(θm)表示d(θm)的共轭。
然后，再根据S2估计得到的相干信源数目K，从M个E(θi)中搜索峰值和其它次峰值，将峰值和次峰值所在方位记作
Figure BDA0003156566660000101
当κ＝K时，以
Figure BDA0003156566660000102
的每个方位值为中心，取峰值或次峰值-3dB宽度区间以内的离散方位值集合作为相干辐射噪声到达方向分布范围，记作
Figure BDA0003156566660000103
当κ＜K时，以
Figure BDA0003156566660000104
的每个方位值为中心，取峰值或次峰值-6dB宽度区间以内的离散方位值集合作为相干辐射噪声到达方向分布范围，记作
Figure BDA0003156566660000105
具体地，所述S4包括：
首先，依据S2估计的相干信源数目K以及S3估计的目标辐射噪声相干信源到达方向
Figure BDA0003156566660000106
确定相干信源到达方向的所有可能组合。记
Figure BDA0003156566660000107
的元素个数分别为m1，m2，…，mκ。当κ＝K时，从κ个集合
Figure BDA0003156566660000108
中各分配一个信源波达方向进行组合，组合数目为每个集合长度的乘积，即
Figure BDA0003156566660000109
当κ＜K时，重新合并集合
Figure BDA0003156566660000111
新集合元素个数为
Figure BDA0003156566660000112
在
Figure BDA0003156566660000113
中做K个相干信源的组合，组合数目为
Figure BDA0003156566660000114
然后，在
Figure BDA0003156566660000115
个波达方向组合中，对每一个相干信源波达方向的组合，记任意一个波达方向组合为
Figure BDA0003156566660000116
查找预先存储的生成矩阵
Figure BDA0003156566660000117
和转换矩阵
Figure BDA0003156566660000118
计算K个相干信源的DFT值：
Figure BDA0003156566660000119
以及模型偏差：
Figure BDA00031565666600001110
最后，对于
Figure BDA00031565666600001111
个波达方向组合，总共有
Figure BDA00031565666600001112
个模型偏差，从中搜索模型偏差最小值，将最小值对应的波达方向组合
Figure BDA00031565666600001113
作为K个相干信源的真实波达方向，并将对应的DFT值
Figure BDA00031565666600001114
中每个元素的幅度平方值作为相干信源强度。作为一种可实施方式，信源强度最大的波达方向可参考视为水下目标辐射噪声的最短路径到达方向。
综上，本发明通过非实时计算，预先存储阵列信号的生成矩阵和转换矩阵，大幅度降低了线谱分量信号空间的多维搜索求解复杂度；并通过盖氏圆方法和常规波束形成方法分别估计相干信源数目和波达方向分布范围，进一步降低线谱分量信号空间的多维搜索求解范围。本发明充分利用了空间拟合类算法对相干信源波达方位估计精度高的优点，又大幅降低了实时计算复杂度，适用于多途条件下对水下目标辐射噪声的相干信源测向定位。
以上所示仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。

Claims (5)

1.一种基于阵列信号数值计算的水下目标辐射噪声相干信源测向方法，其特征在于，包括：
步骤1：生成矩阵和转换矩阵的非实时计算存储：预先将水下目标辐射噪声不同频率线谱的到达方向做数值化离散，然后依据不同信源数目分别做组合排列，最后计算每个组合下的生成矩阵和转换矩阵；
步骤2：目标辐射噪声相干信源数目的实时估计：实时计算各个阵元通道采集的目标辐射噪声中线谱分量的样本协方差矩阵，然后通过盖氏圆信源估计方法估计目标辐射噪声相干信源数目；
步骤3：目标辐射噪声相干信源到达方向分布范围的实时估计：计算各阵元通道采集的目标辐射噪声的线谱分量的DFT值，然后通过常规波束形成方法估计目标辐射噪声相干信源到达方向分布范围；
步骤4：目标辐射噪声相干信源到达方向的实时精确计算：依据目标辐射噪声相干信源数目估计结果和到达方向分布范围，确定目标辐射噪声相干信源到达方向的所有可能组合以及阵列信号数学模型；然后选取预先存储的生成矩阵、转换矩阵与频域阵列信号向量相乘，计算模型系数与模型偏差；最后将模型偏差最小时对应的到达方向组合确定为目标辐射噪声的相干信源真实到达方向，并将此时的模型系数作为相干信源强度。
2.根据权利要求1所述的一种基于阵列信号数值计算的水下目标辐射噪声相干信源测向方法，其特征在于，所述步骤1包括：
将水下目标辐射噪声的不同频率线谱到达方向做数值化离散，离散间隔为Δθ分辨率，离散波达方向的数目为M；
相干辐射噪声到达方向依据信源组合从1至L-3做组合排列，组合个数分别为
Figure FDA0003156566650000011
其中L为阵元数目；
对于每个频率的每个组合，计算转换矩阵：
D(fi，θ1，…，θm)＝[GH(fi，θ1，…，θm)·G(fi，θ1，…，θm)]-1·GH(fi，θ1，…，θm)
其中，1＜m≤M，G(fi，θ1，…，θm)表示在线谱频率fi以及波达方向组合(θ1，…，θm)条件下，阵列信号模型中的生成矩阵；阵列信号模型为：
Figure FDA0003156566650000021
记
Figure FDA0003156566650000022
S(fi)＝(S(fi，θ1) S2(fi，θ2) … Sm(fi，θm))T
X(fi)＝(X1(fi) X2(fi) … XL(fi))T
其中，s(fi)为待求的相干信源DFT值；X(fi)为阵列实际接收的目标辐射噪声的线谱DFT值；G(fi，θ1，…，θm)的第一行元素取值均为1，表示以第一个阵元为参考阵元，其它行元素表示其它阵元接收信号线谱DFT值相对参考阵元接收信号线谱DFT值的增量；G(fi，θ1，θ2，…，θm)的第l行第k列元素取值为：
Figure FDA0003156566650000023
其中，d为阵元间距，满足d≤λi/2，λi为线谱频率fi对应的载波波长；
将生成矩阵G(fi，θ1，θ2，…，θm)和转换矩阵D(fi，θ1，…，θm)进行非实时计算，并存储；
每个频率的线谱需要预先存储的生成矩阵、转换矩阵的元素个数总数均为：
Figure FDA0003156566650000024
其中，L为阵元数目。
3.根据权利要求1所述的一种基于阵列信号数值计算的水下目标辐射噪声相干信源测向方法，其特征在于，所述步骤2包括：
首先，对水下目标辐射噪声进行观测接收和实时采集，做DFT变换，实时计算L个阵元通道的基频线谱分量DFT值的样本协方差矩阵：
C＝X(fi)H·X(fi)
其中X(fi)＝(X1(fi)，X2(fi)，…，XL(fi))，为阵列实际接收的目标辐射噪声的线谱DFT值；L为阵元数目；fi为线谱频率；
然后，根据盖氏圆信源估计方法，由C的前(L-1)行和前(L-1)列构成新矩阵C1，对C1做特征值分解，计算C1的特征值和特征向量；将所述特征值从大到小排序，构造对角矩阵D1＝diag(λ′1，λ′2，…，λ′K，…，λ′L-1)，其中K是待估计的相干辐射噪声信源个数，λ′1、λ′2…λ′K…λ′L-1分别为从大到小排序的C1的特征值；根据所述特征向量构造酉矩阵U1＝[q′1，q′2，…，q′M，…，q′L-1]，q′1、q′2…q′M…q′L-1分别为λ′1、λ′2…λ′K…λ′L-1对应的特征向量，满足
Figure FDA0003156566650000031
再重构酉矩阵
Figure FDA0003156566650000032
用生成的酉矩阵U对样本协方差矩阵C做酉变换：
Figure FDA0003156566650000033
其中，ρ1、ρ2、ρ3…ρL-1、ρ1 *、ρ2 *、ρ3 *…ρL-1 *、cLL为UHCU计算后矩阵中的元素；
则按照下式估计相干信源数目：
Figure FDA0003156566650000034
其中，0≤D(N)≤1，是关于观测时间内共采集的噪声点数N的递减函数；GDE(k)是用来帮助判断信源数目的一个过程量，当k＝k0时，GDE(k0)为第一个比零小的值，则信源数的估计值取K＝k0-1。
4.根据权利要求3所述的一种基于阵列信号数值计算的水下目标辐射噪声相干信源测向方法，其特征在于，所述步骤3包括：
首先，依据X(fi)＝(X1(fi)，X2(fi)，…，XL(fi))做常规波束形成：
d(θm)＝a(θm)·XH(k)，m＝1，2，…，M.
其中，θm为波达方向离散后第m个波达方向，a(θm)为该波达方向的阵列导向矢量；
基于d(θm)计算功率，总共计算得到M个E(θm)：
E(θm)＝d*(θm)·d(θm)，m＝1，2，…，M.
然后，再根据步骤2估计得到的相干信源数目K，从M个E(θm)中搜索峰值和次峰值，将峰值和次峰值所在方位记作
Figure FDA0003156566650000041
当κ＝K时，以
Figure FDA0003156566650000042
的每个方位值为中心，取峰值或次峰值-3dB宽度区间以内的离散方位值集合作为相干辐射噪声到达方向分布范围，记作
Figure FDA0003156566650000043
当κ＜K时，以
Figure FDA0003156566650000044
的每个方位值为中心，取峰值或次峰值-6dB宽度区间以内的离散方位值集合作为相干辐射噪声到达方向分布范围，记作
Figure FDA0003156566650000045
5.根据权利要求4所述的一种基于阵列信号数值计算的水下目标辐射噪声相干信源测向方法，其特征在于，所述步骤4包括：
首先，依据步骤2估计的相干信源数目K以及步骤3估计的目标辐射噪声相干信源到达方向
Figure FDA0003156566650000046
确定相干信源到达方向的所有可能组合；记
Figure FDA0003156566650000047
的元素个数分别为m1，m2，…，mκ；当κ＝K时，从κ个集合
Figure FDA0003156566650000048
中各分配一个信源波达方向进行组合，组合数目为每个集合长度的乘积，即
Figure FDA0003156566650000049
当κ＜K时，重新合并集合
Figure FDA00031565666500000410
新集合元素个数为
Figure FDA00031565666500000411
在
Figure FDA00031565666500000412
中做K个相干信源的组合，组合数目为
Figure FDA00031565666500000413
然后，在
Figure FDA00031565666500000414
个波达方向组合中，对每一个相干信源波达方向的组合，记任意一个波达方向组合为
Figure FDA00031565666500000415
查找预先存储的生成矩阵
Figure FDA00031565666500000416
和转换矩阵
Figure FDA00031565666500000417
计算K个相干信源的DFT值：
Figure FDA00031565666500000418
以及模型偏差：
Figure FDA00031565666500000419
最后，对于
Figure FDA00031565666500000420
个波达方向组合，总共有
Figure FDA00031565666500000421
个模型偏差，从中搜索模型偏差最小值，将最小值对应的波达方向组合
Figure FDA00031565666500000422
作为K个相干信源的真实波达方向，并将对应的DFT值
Figure FDA00031565666500000423
中每个元素的幅度平方值作为相干信源强度。
CN202110780405.1A 2021-07-09 2021-07-09 一种基于阵列信号数值计算的水下目标辐射噪声相干信源测向方法 Active CN113687296B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110780405.1A CN113687296B (zh) 2021-07-09 2021-07-09 一种基于阵列信号数值计算的水下目标辐射噪声相干信源测向方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110780405.1A CN113687296B (zh) 2021-07-09 2021-07-09 一种基于阵列信号数值计算的水下目标辐射噪声相干信源测向方法

Publications (2)

Publication Number Publication Date
CN113687296A true CN113687296A (zh) 2021-11-23
CN113687296B CN113687296B (zh) 2023-05-16

Family

ID=78577178

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110780405.1A Active CN113687296B (zh) 2021-07-09 2021-07-09 一种基于阵列信号数值计算的水下目标辐射噪声相干信源测向方法

Country Status (1)

Country Link
CN (1) CN113687296B (zh)

Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0689178A1 (de) * 1994-06-21 1995-12-27 ANT Nachrichtentechnik GmbH Anordnung zum Empfangen von Signalen, die in einem vorgegebenen Gebiet sich bewegende Objekte aussenden
EP1253434A1 (en) * 2001-04-27 2002-10-30 Mitsubishi Electric Information Technology Centre Europe B.V. Method for estimating a direction of arrival
CN105785322A (zh) * 2014-12-26 2016-07-20 中国科学院声学研究所 一种波束形成方法
CN106093921A (zh) * 2016-07-25 2016-11-09 中国电子科技集团公司第五十四研究所 基于稀疏分解理论的声矢量阵宽带测向方法
CN106249244A (zh) * 2016-07-07 2016-12-21 哈尔滨工程大学 一种声矢量圆阵宽带相干信号源方位估计方法
CN106569180A (zh) * 2016-11-10 2017-04-19 中国人民解放军理工大学 一种基于Prony方法的方位估计算法
CN108375751A (zh) * 2018-01-31 2018-08-07 中国人民解放军战略支援部队信息工程大学 多信源波达方向估计方法
CN108845307A (zh) * 2018-08-02 2018-11-20 西北工业大学 一种基于傅里叶积分法的水下目标辐射噪声测量方法
CN109375154A (zh) * 2018-10-23 2019-02-22 哈尔滨工程大学 一种冲击噪声环境下基于均匀圆阵的相干信号参数估计方法
CN109412982A (zh) * 2018-09-28 2019-03-01 西北工业大学 一种基于信道观测冲激响应模型的多径数目估计方法
CN110007266A (zh) * 2019-04-22 2019-07-12 哈尔滨工程大学 一种冲击噪声下的任意阵列相干源测向方法
WO2019204231A1 (en) * 2018-04-16 2019-10-24 The Regents Of The University Of California Interferometric technique for measuring cerebral blood flow using inexpensive cmos sensors
CN110673085A (zh) * 2019-09-25 2020-01-10 南京航空航天大学 一种均匀面阵下基于快速收敛平行因子的相干信源测向方法

Patent Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0689178A1 (de) * 1994-06-21 1995-12-27 ANT Nachrichtentechnik GmbH Anordnung zum Empfangen von Signalen, die in einem vorgegebenen Gebiet sich bewegende Objekte aussenden
EP1253434A1 (en) * 2001-04-27 2002-10-30 Mitsubishi Electric Information Technology Centre Europe B.V. Method for estimating a direction of arrival
CN105785322A (zh) * 2014-12-26 2016-07-20 中国科学院声学研究所 一种波束形成方法
CN106249244A (zh) * 2016-07-07 2016-12-21 哈尔滨工程大学 一种声矢量圆阵宽带相干信号源方位估计方法
CN106093921A (zh) * 2016-07-25 2016-11-09 中国电子科技集团公司第五十四研究所 基于稀疏分解理论的声矢量阵宽带测向方法
CN106569180A (zh) * 2016-11-10 2017-04-19 中国人民解放军理工大学 一种基于Prony方法的方位估计算法
CN108375751A (zh) * 2018-01-31 2018-08-07 中国人民解放军战略支援部队信息工程大学 多信源波达方向估计方法
WO2019204231A1 (en) * 2018-04-16 2019-10-24 The Regents Of The University Of California Interferometric technique for measuring cerebral blood flow using inexpensive cmos sensors
CN108845307A (zh) * 2018-08-02 2018-11-20 西北工业大学 一种基于傅里叶积分法的水下目标辐射噪声测量方法
CN109412982A (zh) * 2018-09-28 2019-03-01 西北工业大学 一种基于信道观测冲激响应模型的多径数目估计方法
CN109375154A (zh) * 2018-10-23 2019-02-22 哈尔滨工程大学 一种冲击噪声环境下基于均匀圆阵的相干信号参数估计方法
CN110007266A (zh) * 2019-04-22 2019-07-12 哈尔滨工程大学 一种冲击噪声下的任意阵列相干源测向方法
CN110673085A (zh) * 2019-09-25 2020-01-10 南京航空航天大学 一种均匀面阵下基于快速收敛平行因子的相干信源测向方法

Non-Patent Citations (9)

* Cited by examiner, † Cited by third party
Title
HSIEN-TSAI WU 等: ""Source Number Estimators Using Transformed Gerschgorin Radii "", 《IEEE TRANSACTIONS ON SIGNAL PROCESSING》 *
JIE ZHUO 等: ""Joint detection and DOA estimation in multipath channels via multistage constant modulus array"", 《OCEANS \'04 MTS/IEEE TECHNO-OCEAN \'04 (IEEE CAT. NO.04CH37600)》 *
NAN JIANG 等: ""Modulation Recognition of Underwater Acoustic Communication Signals Based on Data Transfer"", 《2019 IEEE 8TH JOINT INTERNATIONAL INFORMATION TECHNOLOGY AND ARTIFICIAL INTELLIGENCE CONFERENCE (ITAIC) 》 *
R. ROSS HAACKE 等: ""Layer stripping of shear-wave splitting in marine PS waves"", 《GEOPHYSICAL JOURNAL INTERNATIONAL》 *
YINGHUA HAN: ""A Method for the Central Direction of Arrival Estimation of Coherently Distributed Source"", 《2009 THIRD INTERNATIONAL SYMPOSIUM ON INTELLIGENT INFORMATION TECHNOLOGY APPLICATION》 *
付淑娟 等: ""基于虚拟阵列改进MUSIC算法的相干信源DoA估计"", 《电讯技术》 *
杨龙: ""水下目标方位估计与稳健自适应波束形成方法研究"", 《中国博士论文全文数据库》 *
林志斌 等: ""基于声传播算子的声源定位实验模拟方法研"", 《应用声学》 *
钱诚: ""相干信源波达方向估计中的若干问题研究"", 《中国博士论文全文数据库》 *

Also Published As

Publication number Publication date
CN113687296B (zh) 2023-05-16

Similar Documents

Publication Publication Date Title
Zheng et al. DOA estimation for coprime linear arrays: An ambiguity-free method involving full DOFs
Zhu et al. Impaired sensor diagnosis, beamforming, and DOA estimation with difference co-array processing
Aghababaiyan et al. High-precision OMP-based direction of arrival estimation scheme for hybrid non-uniform array
CN111123192B (zh) 一种基于圆形阵列和虚拟扩展的二维doa定位方法
CN109633522B (zh) 基于改进的music算法的波达方向估计方法
Wong et al. Direction-finding with sparse rectangular dual-size spatial invariance array
WO1992008149A1 (en) Radar apparatus
CN110197112B (zh) 一种基于协方差修正的波束域Root-MUSIC方法
Baig et al. Comparison of direction of arrival (DOA) estimation techniques for closely spaced targets
CN112130111B (zh) 一种大规模均匀十字阵列中单快拍二维doa估计方法
CN106443571A (zh) 一种保角变换旋转不变子空间波达方向估计方法
CN108828586B (zh) 一种基于波束域的双基地mimo雷达测角优化方法
CN111368256B (zh) 一种基于均匀圆阵的单快拍测向方法
Wu et al. Unambiguous directions of arrival estimation of coherent sources using acoustic vector sensor linear arrays
Muhammad et al. Adaptive diagonal loading technique to improve direction of arrival estimation accuracy for linear antenna array sensors
Zhao et al. Improved capon estimators for DOD and DOA estimation in large array MIMO radar: A random matrix method
Yu et al. A robust minimum variance beamformer with new constraint on uncertainty of steering vector
CN113687296A (zh) 一种基于阵列信号数值计算的水下目标辐射噪声相干信源测向方法
CN114563760B (zh) 一种基于sca阵型的二阶超波束形成方法、设备及介质
CN113391266B (zh) 基于非圆多嵌套阵降维子空间数据融合的直接定位方法
CN112255625B (zh) 基于深度学习的二维角度依赖误差下的一维线阵测向方法
CN111366891B (zh) 一种基于伪协方差矩阵的均匀圆阵单快拍测向方法
CN114265004A (zh) 一种基于子空间对消的干扰下的目标角度估计方法
CN114325565A (zh) 一种基于子空间关系的阵列超分辨测向方法
Gu et al. Sparse linear arrays for estimating and tracking DOAs of signals with known waveforms

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant