主权项 |
一种利用发射栅瓣的成像方法,其特征在于包括下述步骤:1)设计所需的成像阵列,发射阵列为M元ULA,接收阵列为N元ULA,所述的ULA表示均匀直线阵;发射ULA和接收ULA位于同一条直线上,且几何中心重叠;发射ULA为稀疏阵列,接收ULA为密布阵列,发射阵元间距dt和接收阵元间距dr满足dt=Ndr;接收主瓣指向其中某个发射栅瓣或主瓣时,接收波束图的零点或低旁瓣同时指向其它发射栅瓣;2)使用步骤1)设计的成像阵列进行信号发射、接收波束形成和成像处理,其中,目标被建模为P个散射点,第p个散射点的角度为θ<sub>p</sub>,即第p个散射点相对于接收ULA法线方向的夹角,p=1,2,…,P;在一个发射与接收周期中,M个发射阵元同时发射经过加权的脉冲信号s(t)=[a<sub>1</sub>(θ<sub>0</sub>) a<sub>2</sub>(θ<sub>0</sub>) … a<sub>M</sub>(θ<sub>0</sub>)]<sup>H</sup>s<sub>0</sub>(t)<sub>,</sub>其中,t代表时间,s<sub>0</sub>(t)为原始发射信号,a<sub>m</sub>(θ<sub>0</sub>)=A<sub>m</sub>exp[‑j2πf<sub>0</sub>(m‑1)d<sub>t</sub>sinθ<sub>0</sub>/c]为M元发射ULA中第m个发射阵元上的复加权,A<sub>m</sub>为对应的幅度加权,θ<sub>0</sub>代表发射主瓣指向,f<sub>0</sub>为窄带发射信号的中心频率,c为信号传播速度;此时,M元发射ULA指向θ<sub>0</sub>的波束<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>B</mi><mi>T</mi></msub><mrow><mo>(</mo><mi>θ</mi><mo>;</mo><msub><mi>θ</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mi>M</mi></mfrac><mo>|</mo><mfrac><mrow><mi>sin</mi><mo>[</mo><mi>M</mi><mfrac><msub><mi>πd</mi><mi>t</mi></msub><mi>λ</mi></mfrac><mrow><mo>(</mo><mi>sin</mi><mi>θ</mi><mo>-</mo><msub><mrow><mi>sin</mi><mi>θ</mi></mrow><mn>0</mn></msub><mo>)</mo></mrow><mo>]</mo></mrow><mrow><mi>sin</mi><mo>[</mo><mfrac><msub><mi>πd</mi><mi>t</mi></msub><mi>λ</mi></mfrac><mrow><mo>(</mo><mi>sin</mi><mi>θ</mi><mo>-</mo><msub><mrow><mi>sin</mi><mi>θ</mi></mrow><mn>0</mn></msub><mo>)</mo></mrow><mo>]</mo></mrow></mfrac><mo>|</mo><mo>,</mo></mrow>]]></math><img file="FDA0000707138440000011.GIF" wi="956" he="320" /></maths>其中λ是与f<sub>0</sub>对应的信号波长;利用接收ULA采集回波,第n个接收阵元上的回波<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>x</mi><mi>n</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>Σ</mi><mrow><mi>p</mi><mo>=</mo><mn>1</mn></mrow><mi>P</mi></munderover><msub><mi>σ</mi><mi>p</mi></msub><munderover><mi>Σ</mi><mrow><mi>m</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mi>s</mi><mi>m</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>-</mo><msubsup><mi>τ</mi><mrow><mi>t</mi><mo>,</mo><mi>m</mi></mrow><mi>p</mi></msubsup><mo>-</mo><msubsup><mi>τ</mi><mrow><mi>r</mi><mo>,</mo><mi>n</mi></mrow><mi>p</mi></msubsup><mo>)</mo></mrow><mo>+</mo><msub><mi>z</mi><mi>n</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000707138440000012.GIF" wi="956" he="181" /></maths>其中,n=1,2,…,N,σ<sub>p</sub>为第p个散射点的散射强度,<img file="FDA0000707138440000013.GIF" wi="87" he="91" />是第m个发射阵元到第p个散射点的时延,<img file="FDA0000707138440000014.GIF" wi="74" he="79" />是第p个散射点到第n个接收阵元的时延,z<sub>n</sub>(t)代表噪声项;在进行接收波束形成时,接收波束图的主瓣指向位于<img file="FDA0000707138440000015.GIF" wi="588" he="187" />的发射波束图主瓣和栅瓣,第q个波束的输出<img file="FDA0000707138440000016.GIF" wi="534" he="179" />其中,q=1,2,…,Q,Q代表波束数;第n个接收阵元上的复加权<img file="FDA0000707138440000021.GIF" wi="1102" he="196" />A<sub>n</sub>为幅度加权,θ<sub>q</sub>为第q个波束的指向角。 |