发明名称 一种基于相位滤波的波束形成方法
摘要 本发明涉及声纳信号处理领域,主要是一种基于相位滤波的波束形成方法,根据分裂波束的相位差序列,对波束形成的结果进行相位滤波,得到基于相位滤波的波束形成输出。本发明有益的效果:本发明提供了一种新型的,计算方法简单,计算量小,易于实现的基于相位滤波的波束形成方法,充分利用了相位信息,使波束的主瓣变窄,尽可能的抑制了波束主轴以外方向的混响回波信号,提高后续波达时间(TOA)估计的精度。
申请公布号 CN101504458B 申请公布日期 2011.07.20
申请号 CN200910096476.9 申请日期 2009.03.10
申请人 中国船舶重工集团公司第七一五研究所 发明人 丁烽;胡鹏;祝恒年;汪素萍;高博
分类号 G01S7/52(2006.01)I 主分类号 G01S7/52(2006.01)I
代理机构 杭州九洲专利事务所有限公司 33101 代理人 陈继亮
主权项 1.一种基于相位滤波的波束形成方法,其特征在于:该方法包含以下步骤:步骤一:采用M元任意阵型接收阵的接收信号,前置预处理后利用信号的希尔伯特变换,得到阵元输出信号的解析信号为:<maths num="0001"><![CDATA[<math><mrow><mover><mi>x</mi><mo>~</mo></mover><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><msup><mrow><mo>[</mo><msub><mover><mi>x</mi><mo>~</mo></mover><mn>1</mn></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>,</mo><msub><mover><mi>x</mi><mo>~</mo></mover><mn>2</mn></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msub><mover><mi>x</mi><mo>~</mo></mover><mi>M</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>]</mo></mrow><mi>T</mi></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,<img file="FSB00000443897900012.GIF" wi="425" he="52" />为第1,2...,M号阵元输出信号的解析信号,上标“T”表示转置;步骤二:把基阵分成两个子阵,第1,2,...,N号阵元组成阵A,第L,L+1,...,M号阵元组成阵B,那么两个子阵各自输出信号的解析信号为:<maths num="0002"><![CDATA[<math><mrow><msub><mover><mi>x</mi><mo>~</mo></mover><mi>A</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><msup><mrow><mo>[</mo><msub><mover><mi>x</mi><mo>~</mo></mover><mn>1</mn></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>,</mo><msub><mover><mi>x</mi><mo>~</mo></mover><mn>2</mn></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msub><mover><mi>x</mi><mo>~</mo></mover><mi>N</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>]</mo></mrow><mi>T</mi></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0003"><![CDATA[<math><mrow><msub><mover><mi>x</mi><mo>~</mo></mover><mi>B</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><msup><mrow><mo>[</mo><msub><mover><mi>x</mi><mo>~</mo></mover><mi>L</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>,</mo><msub><mover><mi>x</mi><mo>~</mo></mover><mrow><mi>L</mi><mo>+</mo><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msub><mover><mi>x</mi><mo>~</mo></mover><mi>M</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>]</mo></mrow><mi>T</mi></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow></math>]]></maths>步骤三:假设要形成θ方位的波束,则整阵和A、B两个子阵波束形成的驾驶向量分别为:W(θ)=[w<sub>1</sub>(θ),w<sub>2</sub>(θ),…,w<sub>M</sub>(θ)]        (4)W<sub>A</sub>(θ)=[w<sub>1</sub>(θ),w<sub>2</sub>(θ),…,w<sub>N</sub>(θ)]       (5)W<sub>B</sub>(θ)=[w<sub>L</sub>(θ),w<sub>L+1</sub>(θ),…,w<sub>M</sub>(θ)]     (6)其中,w<sub>j</sub>(f,θ)=exp(j2πfτ<sub>j</sub>(θ))为第j号阵元的驾驶向量,τ<sub>j</sub>(θ)为第j号阵元到参考点的时延补偿,j=1,2,…,N,则θ方位的整阵和A、B两个子阵波束输出分别为:<maths num="0004"><![CDATA[<math><mrow><mi>y</mi><mrow><mo>(</mo><mi>&theta;</mi><mo>.</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><mi>W</mi><msup><mrow><mo>(</mo><mi>&theta;</mi><mo>)</mo></mrow><mo>&prime;</mo></msup><mo>*</mo><mover><mi>x</mi><mo>~</mo></mover><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0005"><![CDATA[<math><mrow><msub><mi>y</mi><mi>A</mi></msub><mrow><mo>(</mo><mi>&theta;</mi><mo>,</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>W</mi><mi>A</mi></msub><msup><mrow><mo>(</mo><mi>&theta;</mi><mo>)</mo></mrow><mo>&prime;</mo></msup><mo>*</mo><msub><mover><mi>x</mi><mo>~</mo></mover><mi>A</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0006"><![CDATA[<math><mrow><msub><mi>y</mi><mi>B</mi></msub><mrow><mo>(</mo><mi>&theta;</mi><mo>,</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>W</mi><mi>B</mi></msub><msup><mrow><mo>(</mo><mi>&theta;</mi><mo>)</mo></mrow><mo>&prime;</mo></msup><mo>*</mo><msub><mover><mi>x</mi><mo>~</mo></mover><mi>B</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow></math>]]></maths>步骤四:计算子阵A和B的相位差序列为:<img file="FSB00000443897900018.GIF" wi="1424" he="120" />其中,angle[□]表示对序列求相位,上标“*”表示求共轭,当接收到海底回波信号时,两子阵的相位差序列的值处于零点附近,设定一个相位差门限值δ,使当相位差序列<img file="FSB00000443897900019.GIF" wi="213" he="78" />大于δ时,y(θ,n)为零,则θ方位的所要求得的波束输出结果为:<img file="FSB00000443897900021.GIF" wi="1361" he="141" />
地址 310012 浙江省杭州市华星路96号