一种基于傅里叶积分法的水下目标辐射噪声测量方法技术

本发明专利技术涉及一种基于傅里叶积分法的水下目标辐射噪声测量方法，利用阵元之间的互谱组合出更多不同相位的输出，并将这些不同相位的输出等效为具有不同坐标的虚拟阵元并组成大孔径虚拟阵列，从而将垂直阵的孔径拓展至大于原物理孔径，获得更高的自由度以抑制噪声。在此基础上，对孔径扩展后的虚拟阵列进行空间平滑处理和自适应波束形成，获得优于传统单声压水听器测量方法的水下目标辐射噪声测量结果。本发明专利技术与传统单声压水听器测量法相比，本发明专利技术所提方法利用傅里叶积分法的互相关处理的噪声抑制能力，又获得了虚拟阵列的孔径扩展效果，提高信噪比，可以获得比已有方法更精准的测量结果。

【技术实现步骤摘要】
一种基于傅里叶积分法的水下目标辐射噪声测量方法
本专利技术属于噪声测量领域，特别涉及一种基于傅里叶积分法的水下目标辐射噪声测量方法。
技术介绍
水下目标辐射噪声是影响水下目标隐蔽性的重要指标之一，对其精确测量是进行减振降噪以降低被发现概率的前提。因此需要研究高精度的水下目标辐射噪声测量方法。目前，对于水下目标辐射噪声的测量方法主要是单声压水听器测量法。该方法原理简单，布放时只需要将测量水听器布放得与水下目标几何中心等深度即可，水下目标从距离测量水听器一千米处起航，以固定航向、恒定深度，匀速通过测量水听器。在此过程中，记录测量水听器数据，将记录数据进行分段时间平均来估计水下目标辐射噪声级。但是，随着水下目标辐射噪声级越来越低，测量海区海洋环境背景噪声因航道大吨位航船的增多而变得越来越高，使得测量信噪比变低，导致所测量的水下目标辐射噪声存在较大的误差。
技术实现思路
本专利技术解决的技术问题是：为了克服传统单声压水听器测量法因低信噪比而导致的水下目标辐射噪声源级测量误差较大的问题，本专利技术提出一种基于傅里叶积分法(FIM：Fourierintegralmethod)的辐射噪声测量方法。该方法利用阵元之间的互谱组合出更多不同相位的输出，并将这些不同相位的输出等效为具有不同坐标的虚拟阵元并组成大孔径虚拟阵列，从而将垂直阵的孔径拓展至大于原物理孔径，获得更高的自由度以抑制噪声。在此基础上，对孔径扩展后的虚拟阵列进行空间平滑处理和自适应波束形成，获得优于传统单声压水听器测量方法的水下目标辐射噪声测量结果。本专利技术的技术方案是：一种基于傅里叶积分法的水下目标辐射噪声测量方法，包括以下步骤：步骤一：设水下目标深度不变，目标与水面之间的距离为D1，目标与水底之间的距离为D2；建立N元垂直均匀直线阵ULA(ULA：uniformlineararray)，且定义N元均匀垂直直线阵的几何中心深度为D，水下目标和垂直阵几何中心的水平距离为r，海深h为(D1+D2)，从而得到辐射噪声信号直达波与垂直阵几何中心水平方向形成的夹角即垂直到达角θdw：步骤二：利用垂直阵进行辐射噪声的采集以及获得孔径拓展后的虚拟阵列，包括以下子步骤：子步骤一：利用N元均匀垂直直线阵采集辐射噪声信号，得到N×L维的时域数据采样矩阵x，L为快拍数；任意两个阵元上接收信号的相关输出yi，j(τ)可表示为：yi，j(τ)＝E[xi*(t)xj(t+τ)](2)其中，xi(t)为第i(i为N元接收阵列阵元编号，i＝1，2，...，N)个接收阵元上信号，xj(t)为第j(j为N元接收阵列阵元编号，j＝1，2，...，N)个接收阵元上信号，[]*表示求共轭，E[]表示求时间平均，t表示时间，τ表示时延；根据式(2)，对所有N个接收阵元上的信号进行两两互相关处理和自相关处理，N元ULA可获得N2个互相关输出，且N2个互相关输出可看做N2个虚拟阵元上的输出yi，j(τ)＝{yi，j(τ)}(i＝1，2，...，N，j＝1，2，...，N)(共有N2项)；子步骤二：对N2个虚拟阵元上的输出yi，j(τ)进行去冗余处理，得到2N-1个虚拟阵元上的输出，组成2N-1元虚拟ULA，阵元间距为d；去冗余处理后的虚拟阵元的输出为：其中，i为N元接收阵列阵元编号，i＝1，2，...，N，j为N元接收阵列阵元编号，j＝1，2，...，N，m为编号i和j之差，m＝1-N，2-N，...，0，...，N-2，N-1；表示将满足i-j＝m的所有互相关输出求和；根据式(3)，对N2个虚拟阵元上的输出yi，j(τ)进行去冗余处理，最终可得到2N-1元虚拟ULA的输出即：其中，为去冗余处理后的虚拟阵元上的输出，m＝1-N，2-N，...，0，...，N-2，N-1子步骤三：对2N-1元虚拟ULA的采样协方差矩阵做子阵平滑处理，得到处理后的协方差矩阵RFBSS2N-1元虚拟阵列ULA阵列采样协方差矩阵为：其中：[]H表示求共轭转置；为2N-1元虚拟ULA的输出；划分2N-Q个子阵，每个子阵阵元数为Q，则子阵平滑后的协方差矩阵(表示为RFBSS)为：步骤三：使用最小方差无畸变相应MVDR(MVDR：minimumvariancedistortionlessresponse)波束形成器获得加权向量w并反推出目标辐射噪声源级：子步骤一：利用MVDR波束形成器获得加权向量wMVDR波束形成器可表示为：其中，a(θm)为Q×1维扫描向量：θm∈[-90°，90°]，是波束形成器的指向角度，指向辐射噪声信号的直达波与垂直阵几何中心水平方向形成的夹角即垂直到达角θdw，f为阵列设计频率，d为阵元间距，c为声速；对(7)式进行求解，获得加权向量w。子步骤二：通过加权求和得到Q元虚拟ULA的波束形成输出，并反推出目标辐射噪声源级利用子步骤一得到的加权向量w加权求和并做傅里叶变换得到频域输出RSSL：其中，FFT表示傅里叶变换，[]H表示求共轭转置，为子阵平滑处理后Q元虚拟ULA的输出，为去冗余处理后虚拟阵元上的输出，m＝1-N，2-N，...Q-N；再利用该测量结果反推水下目标的辐射噪声级SSL，即：SSL＝RSSL+TL(9)TL为传播损失，按球面扩展定律计算，即：TL＝20log(r)(10)r为垂直阵几何中心与水下目标的水平距离，由于Q元虚拟ULA阵列长度(Q-1)×d与r相比较小，各虚拟阵元上的传播损失可近似为水下目标到垂直阵几何中心的传播损失。专利技术效果本专利技术的技术效果在于：本专利技术的基本原理和实施方案经过了计算机数值仿真的验证，其结果表明：与传统单声压水听器测量法相比，本专利技术所提方法利用傅里叶积分法的互相关处理的噪声抑制能力，又获得了虚拟阵列的孔径扩展效果，提高信噪比，可以获得比已有方法更精准的测量结果。附图说明图1为本专利技术提出的测量方法的坐标示意图；图2为本专利技术中主要步骤的流程；图3为处理采集的信号获得辐射噪声测量结果的流程；图4(a)为仿真浅海波导环境；(b)为实施实例中利用几何射线对辐射噪声信号直达波的垂直到达角估计结果；图5为实施实例中传统单声压水听器测量辐射噪声的结果；图6为本专利技术提出的测量方法测量辐射噪声的结果；具体实施方式参见图1-图6，本专利技术的主要内容有：1)在浅海波导环境中采用垂直直线阵接收目标辐射噪声信号。垂直阵采用等间隔布阵，阵列长度大于等于1米、小于等于测量场水深的一半。根据测量场的浅海波导环境、声速剖面、水下目标所在深度、距离，对水下目标的垂直到达角进行估计、跟踪。方法是利用几何射线直接预估出辐射噪声信号直达波的垂直到达角；2)对垂直直线阵N个元上的接收信号做互相关处理，可以等效为N2个虚拟阵元上的输出，再对其进行去冗余处理，获得2N-1个虚拟阵元上的输出。3)由于多途效应，垂直直线阵上接收到的直达波信号与反射信号之间具有较高的相关系数，需要采用子阵平滑处理以解相干。再利用子阵平滑后的协方差矩阵进行自适应波束形成，获得虚拟阵列的加权向量，加权求和获得入射信号，反推出水下目标辐射噪声级。4)通过计算机数值仿真给出了传统单声压水听器测量法和本专利技术所提方法的辐射噪声源级测量结果，以此证明了本专利技术所提方法可以获得更高精度的辐射噪声源级测量结果。本专利技术解决现存问题所采用的技术方案可分为以下4个步本文档来自技高网...

【技术保护点】
1.一种基于傅里叶积分法的水下目标辐射噪声测量方法，其特征在于，包括以下步骤：步骤一：设水下目标深度不变，目标与水面之间的距离为D1，目标与水底之间的距离为D2；建立N元垂直均匀直线阵ULA(ULA:uniform linear array)，且定义N元均匀垂直直线阵的几何中心深度为D，水下目标和垂直阵几何中心的水平距离为r，海深h为(D1+D2)，从而得到辐射噪声信号直达波与垂直阵几何中心水平方向形成的夹角即垂直到达角θdw：

【技术特征摘要】
1.一种基于傅里叶积分法的水下目标辐射噪声测量方法，其特征在于，包括以下步骤：步骤一：设水下目标深度不变，目标与水面之间的距离为D1，目标与水底之间的距离为D2；建立N元垂直均匀直线阵ULA(ULA:uniformlineararray)，且定义N元均匀垂直直线阵的几何中心深度为D，水下目标和垂直阵几何中心的水平距离为r，海深h为(D1+D2)，从而得到辐射噪声信号直达波与垂直阵几何中心水平方向形成的夹角即垂直到达角θdw：步骤二：利用垂直阵进行辐射噪声的采集以及获得孔径拓展后的虚拟阵列，包括以下子步骤：子步骤一：利用N元均匀垂直直线阵采集辐射噪声信号，得到N×L维的时域数据采样矩阵x，L为快拍数；任意两个阵元上接收信号的相关输出yi,j(τ)可表示为：yi,j(τ)＝E[xi*(t)xj(t+τ)](2)其中，xi(t)为第i(i为N元接收阵列阵元编号，i＝1,2,…,N)个接收阵元上信号，xj(t)为第j(j为N元接收阵列阵元编号，j＝1,2,…,N)个接收阵元上信号，[]*表示求共轭，E[]表示求时间平均，t表示时间，τ表示时延；根据式(2)，对所有N个接收阵元上的信号进行两两互相关处理和自相关处理，N元ULA可获得N2个互相关输出，且N2个互相关输出可看做N2个虚拟阵元上的输出(共有N2项)；子步骤二：对N2个虚拟阵元上的输出进行去冗余处理，得到2N-1个虚拟阵元上的输出，组成2N-1元虚拟ULA，阵元间距为d；去冗余处理后的虚拟阵元的输出为：其中，i为N元接收阵列阵元编号，i＝1,2,…,N，j为N元接收阵列阵元编号，j＝1,2,…,N，m为编号i和j之差，m＝1-N,2-N,…,0,…,N-2,N-1；表示将满足i-j＝m的所有互相关输出求和；根据式(3)，对N2个虚拟阵元上的输出进行去冗...

【专利技术属性】
技术研发人员：刘雄厚，李晓彬，孙超，蒋国庆，张少东，
申请(专利权)人：西北工业大学，
类型：发明
国别省市：陕西,61