发明名称 一种U型阵波束形成加权方法
摘要 本发明涉及声纳信号处理领域,主要是一种U型阵波束形成加权方法,根据“U”型阵阵元在与波束预成方向垂直的波阵面方向的投影坐标位置关系确定插值点,通过对窗函数在插值点的插值,计算出与“U”型阵阵元位置相匹配的修正窗函数权值。本发明提供了一种新型的,计算方法简单,计算量小,易于实现的“U”型阵波束形成加权方法,把基于均匀阵的窗函数加权推广到非均匀阵加权,有效地抑制了旁瓣,提高了“U”型阵波束形成的性能。
申请公布号 CN101149435A 申请公布日期 2008.03.26
申请号 CN200710156454.8 申请日期 2007.10.23
申请人 中国船舶重工集团公司第七一五研究所 发明人 胡鹏;丁烽;李然威
分类号 G01S7/52(2006.01);G01S7/539(2006.01);G01S15/00(2006.01) 主分类号 G01S7/52(2006.01)
代理机构 杭州九洲专利事务所有限公司 代理人 陈继亮
主权项 1.一种U型阵波束形成加权方法,其特征在于:该方法包含以下步骤:步骤一:采用M元“U”型阵接收信号,对基阵接收信号进行前置预处理,令阵元编号分别为1,2,...,M,则阵元输出信号为:x(n)=[x<sub>1</sub>(n),x<sub>2</sub>(n),...,x<sub>M</sub>(n)]<sup>T</sup>    (1);步骤二:信号处理机对阵元输出信号进行快速傅立叶变换(FFT),将信号从时域变换到频域,频域输出为:X(f)=[x<sub>1</sub>(f),x<sub>2</sub>(f),...,x<sub>m</sub>(f)]<sup>T</sup>    (2)其中x<sub>m</sub>(f)为第m个阵元在频率f的值,m=1,2,...,M,设接收信号为带限信号,则f<sub>l</sub>≤f≤f<sub>h</sub>,f<sub>l</sub>,f<sub>h</sub>分别为信号频率的下限和上限;步骤三:通过计算斜率为tgθ的直线ll与“U”型阵的切点位置P,形成θ方位的波束,根据P点位置,确定第1,2,...,N号阵元参加波束形成,由(2)式选出这N个阵元的频域输出为:<img file="A2007101564540002C1.GIF" wi="1207" he="65" />则波束形成的驾驶向量为:W(f,θ)=[w<sub>1</sub>(f,θ),w<sub>2</sub>(f,θ),...,w<sub>N</sub>(f,θ)]    (4)其中w<sub>n</sub>(f,θ)=exp(j2πfτ<sub>n</sub>(θ)),τ<sub>n</sub>(θ)为第n号阵元到参考点(坐标原点0点)的时延补偿,n=1,2,...,N;步骤四:把“U”型阵中参加波束形成的N个阵元投影到与θ垂直的波阵面,即直线qq方向,得到投影点A<sub>1</sub>,A<sub>2</sub>,...,A<sub>N</sub>;假定在投影点左右两端最边缘的两点,即点A<sub>1</sub>和点A<sub>N</sub>之间形成了一个均匀阵,对该均匀阵进行常规窗函数加权,再对形成的“窗口”进行插值;插值点由各投影点在投影连线,即线A<sub>1</sub>A<sub>N</sub>的位置确定,投影点A<sub>n</sub>(其中,n=1,2,...,N)对应的插值点计算公式如下:<maths num="0001"><math>&lt;mrow&gt;&lt;msub&gt;&lt;mover&gt;&lt;mi&gt;A&lt;/mi&gt;&lt;mo&gt;^&lt;/mo&gt;&lt;/mover&gt;&lt;mi&gt;n&lt;/mi&gt;&lt;/msub&gt;&lt;mo&gt;=&lt;/mo&gt;&lt;mover&gt;&lt;mrow&gt;&lt;msub&gt;&lt;mi&gt;A&lt;/mi&gt;&lt;mn&gt;1&lt;/mn&gt;&lt;/msub&gt;&lt;msub&gt;&lt;mi&gt;A&lt;/mi&gt;&lt;mi&gt;n&lt;/mi&gt;&lt;/msub&gt;&lt;/mrow&gt;&lt;mo&gt;&amp;OverBar;&lt;/mo&gt;&lt;/mover&gt;&lt;mo&gt;*&lt;/mo&gt;&lt;mrow&gt;&lt;mo&gt;(&lt;/mo&gt;&lt;mi&gt;N&lt;/mi&gt;&lt;mo&gt;-&lt;/mo&gt;&lt;mn&gt;1&lt;/mn&gt;&lt;mo&gt;)&lt;/mo&gt;&lt;/mrow&gt;&lt;mo&gt;/&lt;/mo&gt;&lt;mover&gt;&lt;mrow&gt;&lt;msub&gt;&lt;mi&gt;A&lt;/mi&gt;&lt;mn&gt;1&lt;/mn&gt;&lt;/msub&gt;&lt;msub&gt;&lt;mi&gt;A&lt;/mi&gt;&lt;mi&gt;N&lt;/mi&gt;&lt;/msub&gt;&lt;/mrow&gt;&lt;mo&gt;&amp;OverBar;&lt;/mo&gt;&lt;/mover&gt;&lt;mo&gt;,&lt;/mo&gt;&lt;mo&gt;-&lt;/mo&gt;&lt;mo&gt;-&lt;/mo&gt;&lt;mo&gt;-&lt;/mo&gt;&lt;mrow&gt;&lt;mo&gt;(&lt;/mo&gt;&lt;mn&gt;5&lt;/mn&gt;&lt;mo&gt;)&lt;/mo&gt;&lt;/mrow&gt;&lt;/mrow&gt;</math></maths>其中<img file="A2007101564540003C2.GIF" wi="92" he="63" />表示线段A<sub>1</sub>A<sub>n</sub>的长度,<img file="A2007101564540003C3.GIF" wi="101" he="64" />表示线段A<sub>1</sub>A<sub>N</sub>的长度;根据插值点对常规窗函数在[0,N-1]范围内进行插值,插值后的值即为窗函数的修正权值;设所得到的修正窗函数权值为:B=[b<sub>1</sub>,b<sub>2</sub>,...,b<sub>N</sub>]    (6)其中,b<sub>n</sub>为对应第n号阵元的修正窗函数权值,n=1,2,...,N;步骤五:结合步骤三和步骤四的结果,进行频域波束形成得到θ方位的频域波束输出为:<img file="A2007101564540003C4.GIF" wi="1254" he="57" />
地址 310012浙江省杭州市华星路96号
您可能感兴趣的专利