主权项 |
一种基于空域滤波的米波雷达自适应频率选择方法,其特征在于步骤如下:步骤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>Σ</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>Σ</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>θ</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>θ</mi><mi>p</mi></msub><mo>)</mo></mrow><mo>=</mo><mi>exp</mi><mo>[</mo><mo>-</mo><mi>j</mi><mn>2</mn><mi>π</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>θ</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>Σ</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>π</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>∈</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>。 |