发明名称 基于超指向性小孔径圆柱阵的舰船辐射噪声源分辨方法
摘要 本发明提供了一种基于超指向性小孔径圆柱阵的舰船辐射噪声源分辨方法,提出了采用超指向性小孔径圆柱阵实现舰船辐射噪声源分辨;利用测量基阵各阵元接收环境噪声信号之间的相关特性设计出超指向性加权向量,以实现窄波束宽度和高增益处理系统;在六亮点目标模型的假设下,通过计算机仿真给出了这种基于超指向性小孔径圆柱阵舰船辐射噪声源分辨的结果;与常规波束形成进行对比仿真研究,证明了该发明在舰船辐射噪声源分辨问题上的优越性。
申请公布号 CN103438987A 申请公布日期 2013.12.11
申请号 CN201310390863.X 申请日期 2013.08.30
申请人 西北工业大学 发明人 孙超;向龙凤;马远良;杨益新;杨坤德
分类号 G01H17/00(2006.01)I;G01S5/28(2006.01)I 主分类号 G01H17/00(2006.01)I
代理机构 西北工业大学专利中心 61204 代理人 顾潮琪
主权项 1.一种基于超指向性小孔径圆柱阵的舰船辐射噪声源分辨方法,其特征在于包括下述步骤: (1)建立测量水听器阵输出信号模型 均匀圆柱阵由N个阵元个数为M的均匀圆环阵等间隔排列而成,均匀圆环阵的半径为a,相邻两圆环阵的间隔为d;当入射信号为平面波信号时,任意均匀圆环阵阵元相对于坐标X轴的角度为γ<sub>k</sub>=0,β,2β,...,(M-1)β,β=2π/M,则阵元相对于坐标原点的位置矢量P<sub>k,l</sub>=a[cos(γ<sub>k</sub>),sin(γ<sub>k</sub>),(l-1)d]<sup>T</sup>,式中,k=1,2,...,M表示单个均匀圆环阵阵元编号,l=1,2,...,N为圆环阵编号,符号T表示求转置;归一化信号方向矢量<img file="FDA0000375022530000011.GIF" wi="940" he="78" />其中,信号方位角为θ,俯仰角为<img file="FDA0000375022530000012.GIF" wi="55" he="37" />则声程差Δd<sub>k,l</sub>=P<sub>s</sub><sup>T</sup>P<sub>k,l</sub>;当入射信号为CW信号,基阵输出信号<img file="FDA0000375022530000013.GIF" wi="571" he="131" />式中,<img file="FDA0000375022530000014.GIF" wi="228" he="78" />分别表示第p=kM+l,(p=0,1,...,MN-1)号阵元接收到的时域信号和权值,其中,上标*为求共轭,t为时间,W=[w<sub>0</sub>,w<sub>1</sub>,...,w<sub>MN-1</sub>]<sup>T</sup>为加权向量,<img file="FDA0000375022530000015.GIF" wi="633" he="78" />为基阵接收信号矩阵,当入射信号频率为f<sub>0</sub>时,阵元输出信号可进一步写为<img file="FDA0000375022530000016.GIF" wi="309" he="50" />式中,A为信号幅值,U为信号导向矢量,<img file="FDA0000375022530000017.GIF" wi="581" he="84" />其中,<img file="FDA0000375022530000018.GIF" wi="277" he="72" />为相位,c为声速,Δd<sub>p</sub>=Δd<sub>k,l</sub>;可得基阵输出信号<img file="FDA0000375022530000019.GIF" wi="462" he="66" />(2)求解最优加权向量 定义基阵波束幅度响应B(θ,φ)=|W<sup>H</sup>U|,当信号入射方向一定时,基阵波束输出的幅度响应由加权向量W确定;在声纳信号处理中,假定圆柱阵接收数据向量X=X<sub>d</sub>+X<sub>n</sub>,式中,X<sub>d</sub>和X<sub>n</sub>分别为期望信号和噪声信号,则基阵输出功率谱B<sub>f</sub>(θ,φ)=W<sup>H</sup>R<sub>X</sub>W=W<sup>H</sup>R<sub>d</sub>W+W<sup>H</sup>R<sub>n</sub>W,式中,R<sub>X</sub>,R<sub>d</sub>和R<sub>n</sub>分别为接收数据向量、信号向量和噪声向量的协方差矩阵,R<sub>X</sub>=E{X<sup>H</sup>X},<img file="FDA00003750225300000110.GIF" wi="333" he="78" /><img file="FDA00003750225300000111.GIF" wi="357" he="78" />定义基阵的空间增益为基阵的输出信噪比和参考阵元输出信噪比,则空间增益<img file="FDA0000375022530000021.GIF" wi="634" he="133" />式中,P为归一化信号协方差矩阵,P=UU<sup>H</sup>,Q为归一化噪声协方差矩阵,令增益G为最大值,可以得到<img file="FDA0000375022530000022.GIF" wi="327" he="148" />当加权向量W满足上式时,基阵输出信噪比最大,即可获得最大的空间处理增益;基阵的常规处理增益<img file="FDA0000375022530000023.GIF" wi="276" he="125" />在空间均匀噪声场中,当阵元间隔等于λ/2时,λ为信号波长,各阵元接收到的噪声信号互不相关,Q退化成单位阵,这时,G<sub>0</sub>=G<sub>c</sub>;当阵元间隔小于λ/2时,G<sub>0</sub>&gt;G<sub>c</sub>,基阵各阵元接收噪声信号是相关的,在各向同性噪声场中,其相关系数<img file="FDA0000375022530000024.GIF" wi="389" he="153" />式中,k=ω/c为波数,ω为角频率,Δl<sub>p,q</sub>表示编号为p和q阵元之间的距离,p,q=0,1,2,...,MN-1;p和q阵元的位置坐标可表示为<img file="FDA0000375022530000025.GIF" wi="929" he="150" /><img file="FDA0000375022530000026.GIF" wi="913" he="150" />式中,θ<sub>p</sub>和θ<sub>q</sub>分别为第p号、第q号阵元与坐标轴X轴的水平夹角,ceil(·)为向上取整。则有<img file="FDA0000375022530000027.GIF" wi="1216" he="163" />可得噪声相关矩阵为<img file="FDA0000375022530000028.GIF" wi="836" he="300" />式中,ρ<sub>p,q</sub>=ρ<sub>q,p</sub>,ρ<sub>p,p</sub>=1,则Q为对称正定矩阵,可得具有超指向性基阵的最优权及最大阵增益值; (3)实现舰船辐射噪声多源分辨 将步骤2计算出的最优加权向量代入基阵输出功率谱,可得舰船辐射噪声源方位谱,从而实现舰船辐射噪声多源分辨。 
地址 710072 陕西省西安市友谊西路127号