发明名称 基于子带处理的宽带信号目标检测多普勒色散消除方法
摘要 本发明涉及一种基于子带处理的宽带信号目标检测多普勒色散消除方法,属于雷达目标探测技术领域。本发明方法将直达波分解为多个子带,在各子带内分别进行相参积累,然后将各子带的相参积累结果进行复数叠加,即得到最终的CAF输出结果。为使各子带积累出的目标峰值在同一个多普勒单元内,需要根据各子带的载波频率,选择不同的积累时间。与传统的Keystone变换方法相比,子带方法不存在信号波形不变的假设,避免了因此而带来的信噪比损失。本发明通过子带方法明显消除了目标多普勒色散,有效提高了相参积累增益,进而提高了系统的检测性能,特别适合于采用长时间积累探测高速目标的应用场合。
申请公布号 CN103278806B 申请公布日期 2015.04.22
申请号 CN201310227099.4 申请日期 2013.06.08
申请人 北京理工大学 发明人 单涛;刘升恒;陶然;张果;冯远
分类号 G01S7/285(2006.01)I;G01S7/32(2006.01)I 主分类号 G01S7/285(2006.01)I
代理机构 代理人
主权项 基于子带处理的宽带信号目标检测多普勒色散消除方法,其特征在于,具体包括如下步骤:步骤一,子带分解将天线接收的宽带直达波信号通过子带分解滤波器组,分解为K个子带:r<sub>k</sub>(n)=r(n)*h<sub>k</sub>(n)    (4)其中,r<sub>k</sub>(n)为直达波信号矢量r(n)经过子带分解得到的第k个子带的信号;所述子带分解滤波器组的数学描述为:h<sub>k</sub>(n)=h(n)e<sup>‑j2πkn/K</sup>    (5)式中k=0,1,2,...K‑1,K为子带数量;h<sub>k</sub>(n)为第k个子带分解滤波器;h(n)是子带分解滤波器的原型滤波器;步骤二,共轭复乘将步骤一得到的各个子带中的直达波信号,分别和另一天线接收的宽带回波信号做共轭复乘,第k个子带共轭复乘公式为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>x</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>&tau;</mi><mo>,</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><mi>e</mi><mrow><mo>(</mo><mi>n</mi><mo>+</mo><mi>&tau;</mi><mo>)</mo></mrow><msubsup><mi>r</mi><mi>k</mi><mo>*</mo></msubsup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000616297630000011.GIF" wi="1255" he="85" /></maths>其中,e(n)为目标回波信号,τ为时延;步骤三,信号下抽将步骤二共轭复乘后得到的各路信号依次进行抗混叠低通滤波、M倍下抽取;x<sub>k</sub>(τ,n)下抽后所得信号为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>y</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>&tau;</mi><mo>,</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>l</mi><mo>=</mo><mo>-</mo><mo>&infin;</mo></mrow><mo>&infin;</mo></munderover><mi>g</mi><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow><msub><mi>x</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>&tau;</mi><mo>,</mo><mi>Mn</mi><mo>-</mo><mi>l</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000616297630000012.GIF" wi="1309" he="157" /></maths>其中M为下抽倍数;g(l)为理想低通滤波器,其DTFT形式为<img file="FDA0000616297630000013.GIF" wi="1216" he="180" />步骤四,相参积累时间调整对步骤三下抽后的各路子带信号分别进行相参积累时间调整,使得不同子带k得到的目标峰值处于相同的多普勒单元内;具体方法为:步骤4.1,求得各路子带信号的相参积累的点数多普勒频率为f<sub>d</sub>的目标,所处的多普勒单元位置为:N<sub>Dop_bin</sub>=f<sub>d</sub>/Δf=f<sub>d</sub>T<sub>int</sub>    (9)其中,目标检测的频率分辨率Δf=1/T<sub>int</sub>,T<sub>int</sub>为相参积累时间;根据目标回波信号的多普勒频率公式<img file="FDA0000616297630000021.GIF" wi="430" he="150" />其中v<sub>r</sub>为目标的径向运动速度,c为光速,f<sub>c</sub>为载波频率,λ为载波波长,进一步得到目标所在N<sub>Dop_bin</sub>为:N<sub>Dop_bin</sub>=2v<sub>r</sub>f<sub>c</sub>T<sub>int</sub>/c    (10)每个子带的相参积累时间T<sub>int</sub>正比于1/f<sub>c</sub>,相参积累的点数为:<img file="FDA0000616297630000022.GIF" wi="1148" he="198" />其中f<sub>s</sub>为采样率,f<sub>c_k</sub>是第k个子带的载波频率,f<sub>center</sub>是整个数字电视信号载波的中心频率;步骤4.2,根据各路子带信号的相参积累的点数,截取点数对应长度的各路子带信号,作为参与子带互模糊函数计算的子带信号;其中,每路子带信号的起始截取位置相同;步骤五,计算各子带CAF根据步骤四调整得到的各子带相参积累时间,采用CAF计算得到各子带的相参积累结果为:χ<sub>k</sub>(τ,f<sub>d</sub>)=DFT{y<sub>k</sub>(τ,n)}    (12)步骤六,子带相参积累综合将各个子带所得的相参积累结果复数值相加,得到已消除多普勒色散的目标检测结果:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>&chi;</mi><mi>&Sigma;</mi></msub><mrow><mo>(</mo><mi>&tau;</mi><mo>,</mo><msub><mi>f</mi><mi>d</mi></msub><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mi>&chi;</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>&tau;</mi><mo>,</mo><msub><mi>f</mi><mi>d</mi></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>13</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000616297630000023.GIF" wi="1144" he="164" /></maths>至此,完成宽带信号目标检测。
地址 100081 北京市海淀区中关村南大街5号