发明名称 一种快速宽带频域扩展拖曳阵波束形成方法
摘要 本发明提出一种快速宽带频域扩展拖曳阵波束形成方法,该方法将ETAM被动合成孔径技术和宽带频域波束形成方法结合在一起,包含如下步骤:1)用线阵接收空间信号,得到M个阵元的时域信号;取重叠阵元数目为N;2)对第k快拍和第k+1快拍数据在时间域上做快速傅里叶变换;3)确定目标辐射信号的频带范围,根据相邻任意两快拍的重叠阵元计算相移参数:4)对不同频率分量分别进行相位补偿得到虚拟的合成阵元;5)重复以上步骤,进行K次的孔径合成;6)对合成的阵元信号在空间域上补零做FFT变换,然后将获得数据的零频分量移至谱中心;7)在每一个频点进行校正,得到不同频点的波束输出;8)合成不同频率的波束形成结果,从而实现目标的精确定位。
申请公布号 CN101813772A 申请公布日期 2010.08.25
申请号 CN201010163261.7 申请日期 2010.04.30
申请人 中国科学院声学研究所;北京中科海跃科技有限责任公司 发明人 李峥;薛山花;田彪;李宇;黄海宁
分类号 G01S7/539(2006.01)I 主分类号 G01S7/539(2006.01)I
代理机构 北京法思腾知识产权代理有限公司 11318 代理人 杨小蓉;高宇
主权项 一种快速宽带频域扩展拖曳阵波束形成方法,该方法将ETAM被动合成孔径技术和宽带频域波束形成方法结合在一起,可用于对微弱目标的检测和定位,同时速度较快,便于实时处理;所述的方法包含如下步骤:1)用线阵接收空间信号,得到M个阵元的时域信号;取重叠阵元数目为N,其中N大于零且小于阵元数目,计算两次快拍数据的时间间隔阵τ,即阵运动M-N个阵元距离需要的时间;2)对第k快拍和第k+1快拍数据在时间域上做快速傅里叶变换;3)确定目标辐射信号的频带范围[fmin,fmax],根据第k快拍和第k+1快拍的重叠阵元计算相移参数ψk+1(fi): <mrow> <msub> <mi>&psi;</mi> <mrow> <mi>k</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mi>arg</mi> <mo>{</mo> <mrow> <mfrac> <mrow> <munderover> <mi>&Sigma;</mi> <mrow> <mi>n</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <msub> <mi>X</mi> <mrow> <mi>k</mi> <mo>,</mo> <mi>M</mi> <mo>-</mo> <mi>N</mi> <mo>+</mo> <mi>n</mi> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>&times;</mo> <msubsup> <mi>X</mi> <mrow> <mi>k</mi> <mo>+</mo> <mn>1</mn> <mo>,</mo> <mi>n</mi> </mrow> <mo>*</mo> </msubsup> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>&times;</mo> <msub> <mi>&rho;</mi> <mi>n</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> </mrow> <mrow> <munderover> <mi>&Sigma;</mi> <mrow> <mi>n</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <msub> <mi>&rho;</mi> <mi>n</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> </mrow> </mfrac> </mrow> <mo>}</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow>k表示快拍序号,n表示阵元序号,fi表示频率,ρn(fi)表示第n阵元频率fi的补偿因子: <mrow> <msub> <mi>&rho;</mi> <mi>n</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <mo>|</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mo>-</mo> <mi>Q</mi> <mo>/</mo> <mn>2</mn> </mrow> <mrow> <mi>Q</mi> <mo>/</mo> <mn>2</mn> </mrow> </munderover> <msub> <mi>X</mi> <mrow> <mi>k</mi> <mo>,</mo> <mi>M</mi> <mo>-</mo> <mi>N</mi> <mo>+</mo> <mi>n</mi> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>&times;</mo> <msubsup> <mi>X</mi> <mrow> <mi>k</mi> <mo>+</mo> <mn>1</mn> <mo>,</mo> <mi>n</mi> </mrow> <mo>*</mo> </msubsup> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>|</mo> </mrow> <msqrt> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mo>-</mo> <mi>Q</mi> <mo>/</mo> <mn>2</mn> </mrow> <mrow> <mi>Q</mi> <mo>/</mo> <mn>2</mn> </mrow> </munderover> <msup> <mrow> <mo>|</mo> <msub> <mi>X</mi> <mrow> <mi>k</mi> <mo>,</mo> <mi>M</mi> <mo>-</mo> <mi>N</mi> <mo>+</mo> <mi>n</mi> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>|</mo> </mrow> <mn>2</mn> </msup> <mo>&times;</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mo>-</mo> <mi>Q</mi> <mo>/</mo> <mn>2</mn> </mrow> <mrow> <mi>Q</mi> <mo>/</mo> <mn>2</mn> </mrow> </munderover> <msup> <mrow> <mo>|</mo> <msubsup> <mi>X</mi> <mrow> <mi>k</mi> <mo>+</mo> <mn>1</mn> <mo>,</mo> <mi>n</mi> </mrow> <mo>*</mo> </msubsup> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>|</mo> </mrow> <mn>2</mn> </msup> </msqrt> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow>4)对第k+1快拍、第M-N+1~M阵元信号的不同频率分量分别进行相位补偿,将其虚拟成第k快拍的第M+1~M+N阵元;依次类推,这N个阵元信号可以虚拟成为第1快拍的第M+N*(k-1)+1~M+N*K阵元;5)重复以上步骤,进行K次的孔径合成,可以虚拟得到MM=M+N*K个阵元接收的信号在[fmin,fmax]范围内的频域响应;6)对MM个阵元信号在空间域上补零,再对其在空间域上做FFT变换,然后将获得数据的零频分量移至谱中心;7)在处理带宽内的每一个频点上依据频率——波束网格对其进行校正,得到不同频点的波束输出;8)合成不同频率的波束形成结果,得到MM个阵元在不同方位上的功率输出,从而实现目标的精确定位。
地址 100190 北京市海淀区北四环西路21号