发明名称 一种基于空域滤波的米波雷达自适应频率选择方法
摘要 本发明涉及一种基于空域滤波的米波雷达自适应频率选择方法,跟据米波雷达同样具有较低频段的特点,对其进行频率分段,数字波束形成,时域FFT处理,最后,通过频谱比较,建立最佳工作频率库,按照用户的不同要求选择工作频率,本发明采用不预置最佳工作频段的处理流程,对所有频段的数据都进行DBF空域滤波处理;同时考虑到常规目标在米波波段谐振散射的特性,在判别准则上,不以输出平均功率最小为准则选择工作频率,而是建立最佳工作频率库,结合谐振频率数据库进行频率选择,从而保证雷达最佳的检测性能。该方法适用于米波雷达在复杂电磁环境下的自适应变频对抗外部干扰。
申请公布号 CN103033797B 申请公布日期 2014.08.13
申请号 CN201210545575.2 申请日期 2012.12.17
申请人 西安电子工程研究所 发明人 游俊;强勇;黄金杰;黄志忠;罗丁利
分类号 G01S7/36(2006.01)I 主分类号 G01S7/36(2006.01)I
代理机构 西北工业大学专利中心 61204 代理人 王鲜凯
主权项 一种基于空域滤波的米波雷达自适应频率选择方法,其特征在于步骤如下:步骤1频率分段:在雷达休止期对M个天线通道分别同时以F<sub>s</sub>采样率采集数据,然后进行数字滤波,将接收信号划分为N个频率段;所述第n个数字滤波器通带范围为f<sub>n,b</sub>~f<sub>n,e</sub>,第m个通道采集的数据经过滤波器之后的输出为<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>y</mi><mrow><mi>m</mi><mo>,</mo><mi>n</mi></mrow></msub><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><msub><mi>l</mi><mi>n</mi></msub><mo>=</mo><mn>0</mn></mrow><mrow><msub><mi>L</mi><mi>n</mi></msub><mo>-</mo><mn>1</mn></mrow></munderover><msub><mi>a</mi><mi>n</mi></msub><mrow><mo>(</mo><msub><mi>l</mi><mi>n</mi></msub><mo>)</mo></mrow><msub><mi>x</mi><mi>m</mi></msub><mrow><mo>(</mo><mi>l</mi><mo>-</mo><msub><mi>l</mi><mi>n</mi></msub><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000481474920000011.GIF" wi="567" he="153" /></maths>其中:a<sub>n</sub>(l<sub>n</sub>)是根据通带范围设计的滤波器系数,L<sub>n</sub>表示滤波器的阶数,x<sub>m</sub>(l)(l=0,1,2,…)表示第m(m=1,2,…,M)个天线通道经过低噪放之后采集的数据;步骤2数字波束形成:对第n个数字滤波器的输出按照空域划分进行DBF,<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>z</mi><mrow><mi>p</mi><mo>,</mo><mi>n</mi></mrow></msub><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>m</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msubsup><mi>w</mi><mrow><mi>m</mi><mo>,</mo><mi>n</mi></mrow><mo>*</mo></msubsup><mrow><mo>(</mo><msub><mi>&theta;</mi><mi>p</mi></msub><mo>)</mo></mrow><msub><mi>y</mi><mrow><mi>m</mi><mo>,</mo><mi>n</mi></mrow></msub><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000481474920000012.GIF" wi="559" he="141" /></maths>其中w<sub>m,n</sub>(θ<sub>p</sub>),m=1,2,…,M表示天线阵在θ<sub>p</sub>方向上形成波束f<sub>n,0</sub>频率所对应的权值,f<sub>n,0</sub>∈[f<sub>n,b</sub>,f<sub>n,e</sub>];非自适应权值<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>w</mi><mrow><mi>m</mi><mo>,</mo><mi>n</mi></mrow></msub><mrow><mo>(</mo><msub><mi>&theta;</mi><mi>p</mi></msub><mo>)</mo></mrow><mo>=</mo><mi>exp</mi><mo>[</mo><mo>-</mo><mi>j</mi><mn>2</mn><mi>&pi;</mi><mrow><mo>(</mo><mi>m</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mfrac><mi>d</mi><mi>c</mi></mfrac><msub><mi>f</mi><mrow><mi>n</mi><mo>,</mo><mn>0</mn></mrow></msub><msub><mrow><mi>sin</mi><mi>&theta;</mi></mrow><mi>p</mi></msub><mo>]</mo><mo>,</mo></mrow>]]></math><img file="FDA0000481474920000013.GIF" wi="872" he="145" /></maths>其中d是阵元间距;c是光速;步骤3时域FFT处理:按照每D点抽取对z<sub>p,n</sub>(l)进行K点的FFT<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>Z</mi><mrow><mi>p</mi><mo>,</mo><mi>n</mi></mrow></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>l</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>K</mi><mo>-</mo><mn>1</mn></mrow></munderover><msub><mi>z</mi><mrow><mi>p</mi><mo>,</mo><mi>n</mi></mrow></msub><mrow><mo>(</mo><mi>lD</mi><mo>)</mo></mrow><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mi>j</mi><mn>2</mn><mi>&pi;</mi><mfrac><mi>kl</mi><mi>K</mi></mfrac><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000481474920000014.GIF" wi="738" he="145" /></maths>其中D满足<img file="FDA0000481474920000015.GIF" wi="507" he="126" />k=0,1,…,K‑1,|Z<sub>p,n</sub>(k)|用来描述在θ<sub>p</sub>方向上干扰噪声在频率<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mrow><mo>(</mo><mi>i</mi><mo>+</mo><mfrac><mi>k</mi><mi>K</mi></mfrac><mo>)</mo></mrow><mfrac><msub><mi>F</mi><mi>s</mi></msub><mi>D</mi></mfrac><mo>&Element;</mo><mo>[</mo><msub><mi>f</mi><mrow><mi>n</mi><mo>,</mo><mi>b</mi></mrow></msub><mo>,</mo><msub><mi>f</mi><mrow><mi>n</mi><mo>,</mo><mi>e</mi></mrow></msub><mo>]</mo></mrow>]]></math><img file="FDA0000481474920000016.GIF" wi="481" he="145" /></maths>的幅度谱,i=i<sub>b</sub>或i<sub>e</sub>,<img file="FDA0000481474920000017.GIF" wi="695" he="92" /><img file="FDA00004814749200000110.GIF" wi="79" he="81" />是向下取整;步骤4,频谱比较,得到干扰噪声最小的频率:对FFT处理的输出进行比较,得到θ<sub>p</sub>方向上f<sub>n,e</sub>~f<sub>n,b</sub>频段内的干扰噪声最小的频率,记为f<sub>n,p</sub>;若k<sub>0</sub>满足条件<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><munder><mi>min</mi><mrow><mi>k</mi><mo>=</mo><mn>0</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>K</mi><mo>-</mo><mn>1</mn></mrow></munder><mi>E</mi><mo>{</mo><msup><mrow><mo>|</mo><msub><mi>Z</mi><mrow><mi>p</mi><mo>,</mo><mi>n</mi></mrow></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup><mo>}</mo><mo>,</mo></mrow>]]></math><img file="FDA0000481474920000018.GIF" wi="454" he="121" /></maths>则<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msub><mi>f</mi><mrow><mi>n</mi><mo>,</mo><mi>p</mi></mrow></msub><mo>=</mo><mrow><mo>(</mo><msub><mi>i</mi><mn>0</mn></msub><mo>+</mo><mfrac><msub><mi>k</mi><mn>0</mn></msub><mi>K</mi></mfrac><mo>)</mo></mrow><mfrac><msub><mi>F</mi><mi>s</mi></msub><mi>D</mi></mfrac><mo>,</mo></mrow>]]></math><img file="FDA0000481474920000019.GIF" wi="395" he="146" /></maths>假设k<sub>b</sub>=[DKf<sub>n,b</sub>/F<sub>s</sub>]‑Ki<sub>b</sub>,k<sub>e</sub>=[DKf<sub>n,e</sub>/F<sub>s</sub>]‑Ki<sub>e</sub>,[·]是四舍五入式的取整;当k<sub>e</sub>>k<sub>0</sub>>k<sub>b</sub>时,i<sub>0</sub>=i<sub>b</sub>;当k<sub>0</sub>>k<sub>b</sub>>k<sub>e</sub>时,i<sub>0</sub>=i<sub>b</sub>;当k<sub>0</sub><k<sub>e</sub><k<sub>b</sub>时,i<sub>0</sub>=i<sub>e</sub>。
地址 710100 陕西省西安市长安区凤栖东路