发明名称 基于MIMO的双基地雷达多通道联合降维杂波抑制方法
摘要 本发明公开了一种基于MIMO的双基地雷达多通道联合降维杂波抑制方法,主要解决现有双基地MIMO雷达距离依赖杂波抑制方法运算量大,且需要的独立同分布样本数目高的问题。其实现步骤为:(1)用发射波形对雷达的回波数据进行匹配滤波;(2)构造多通道联合降维矩阵并对接收数据进行降维处理;(3)用降维后的数据估计出杂波协方差矩阵;(4)根据空时自适应处理原理得到最优权矢量;(5)用最优权对降维后的数据加权,抑制背景杂波,检测出目标信号。本发明与现有全维处理方法相比,具有计算复杂度低,对独立同分布样本数目的要求低和抑制杂波性能好的优点,可应用于MIMO的双基地雷达地面目标检测。
申请公布号 CN103353592B 申请公布日期 2015.03.04
申请号 CN201310244382.8 申请日期 2013.06.19
申请人 西安电子科技大学 发明人 李军;刘长赞;党博;伊孟磊;廖桂生
分类号 G01S7/36(2006.01)I 主分类号 G01S7/36(2006.01)I
代理机构 陕西电子工业专利中心 61205 代理人 王品华;朱红星
主权项 一种基于MIMO的双基地雷达多通道联合降维杂波抑制方法,包括如下步骤:(1)利用正侧视双基地MIMO雷达模式,对雷达接收机端N个天线在一次相干处理间隔内接收到的回波数据分别用M个发射天线阵元的发射波形<img file="FDA0000587234500000011.GIF" wi="59" he="78" />进行匹配滤波;将每个接收天线匹配滤波后的回波数据首尾相连,得到MNK×1维的空时数据矢量y:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>y</mi><mo>=</mo><msub><mi>&rho;</mi><mi>T</mi></msub><msub><mi>b</mi><mi>T</mi></msub><mrow><mo>(</mo><msub><mi>f</mi><mrow><mi>t</mi><mo>,</mo><mi>T</mi></mrow></msub><mo>,</mo><msub><mi>f</mi><mrow><mi>r</mi><mo>,</mo><mi>T</mi></mrow></msub><mo>,</mo><msub><mi>f</mi><mrow><mi>d</mi><mo>,</mo><mi>T</mi></mrow></msub><mo>)</mo></mrow><mo>+</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mrow><msub><mi>N</mi><mi>c</mi></msub><mo>-</mo><mn>1</mn></mrow></munderover><msub><mi>&rho;</mi><mi>i</mi></msub><msub><mi>b</mi><mi>i</mi></msub><mrow><mo>(</mo><msub><mi>f</mi><mrow><mi>t</mi><mo>,</mo><mi>i</mi></mrow></msub><mo>,</mo><msub><mi>f</mi><mrow><mi>r</mi><mo>,</mo><mi>i</mi></mrow></msub><mo>,</mo><msub><mi>f</mi><mrow><mi>d</mi><mo>,</mo><mi>i</mi></mrow></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>y</mi><mi>w</mi></msub><mo>,</mo></mrow>]]></math><img file="FDA0000587234500000012.GIF" wi="925" he="124" /></maths>其中,m=1,2,...,M,上标*表示共轭,K为一次相干处理间隔内的脉冲数,ρ<sub>T</sub>为目标的反射系数,<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>b</mi><mi>T</mi></msub><mrow><mo>(</mo><msub><mi>f</mi><mrow><mi>t</mi><mo>,</mo><mi>T</mi></mrow></msub><mo>,</mo><msub><mi>f</mi><mrow><mi>r</mi><mo>,</mo><mi>T</mi></mrow></msub><mo>,</mo><msub><mi>f</mi><mrow><mi>d</mi><mo>,</mo><mi>T</mi></mrow></msub><mo>)</mo></mrow><mo>=</mo><msub><mi>a</mi><mi>t</mi></msub><mrow><mo>(</mo><msub><mi>f</mi><mrow><mi>t</mi><mo>,</mo><mi>T</mi></mrow></msub><mo>)</mo></mrow><mo>&CircleTimes;</mo><msub><mi>a</mi><mi>r</mi></msub><mrow><mo>(</mo><msub><mi>f</mi><mrow><mi>r</mi><mo>,</mo><mi>T</mi></mrow></msub><mo>)</mo></mrow><mo>&CircleTimes;</mo><msub><mi>a</mi><mi>d</mi></msub><mrow><mo>(</mo><msub><mi>f</mi><mrow><mi>d</mi><mo>,</mo><mi>T</mi></mrow></msub><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000587234500000016.GIF" wi="846" he="66" /></maths>为目标的导向矢量,a<sub>i</sub>(f<sub>t,T</sub>)为目标发射阵列导向矢量,a<sub>r</sub>(f<sub>r,T</sub>)为目标接受阵列导向矢量,a<sub>d</sub>(f<sub>d,T</sub>)为目标多普勒导向矢量,符号<img file="FDA0000587234500000017.GIF" wi="118" he="55" />表示Kronecker积,N<sub>c</sub>为杂波点源的个数,ρ<sub>i</sub>为第i个杂波点的反射系数,<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>b</mi><mi>i</mi></msub><mrow><mo>(</mo><msub><mi>f</mi><mrow><mi>t</mi><mo>,</mo><mi>i</mi></mrow></msub><mo>,</mo><msub><mi>f</mi><mrow><mi>r</mi><mo>,</mo><mi>i</mi></mrow></msub><mo>,</mo><msub><mi>f</mi><mrow><mi>d</mi><mo>,</mo><mi>i</mi></mrow></msub><mo>)</mo></mrow><mo>=</mo><msub><mi>a</mi><mi>t</mi></msub><mrow><mo>(</mo><msub><mi>f</mi><mrow><mi>t</mi><mo>,</mo><mi>i</mi></mrow></msub><mo>)</mo></mrow><mo>&CircleTimes;</mo><msub><mi>a</mi><mi>r</mi></msub><mrow><mo>(</mo><msub><mi>f</mi><mrow><mi>r</mi><mo>,</mo><mi>i</mi></mrow></msub><mo>)</mo></mrow><mo>&CircleTimes;</mo><msub><mi>a</mi><mi>d</mi></msub><mrow><mo>(</mo><msub><mi>f</mi><mrow><mi>d</mi><mo>,</mo><mi>i</mi></mrow></msub><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000587234500000018.GIF" wi="789" he="66" /></maths>为第i个杂波点的空时二维导向矢量,a<sub>t</sub>(f<sub>t,i</sub>)为目标的发射阵列导向矢量,a<sub>r</sub>(f<sub>r,i</sub>)为目标的接受阵列导向矢量,a<sub>d</sub>(f<sub>d,i</sub>)为目标的多普勒导向矢量,y<sub>w</sub>为噪声的空时数据矢量;所述的正侧视双基地MIMO模式,是指雷达的发射机和接收机分置于不同地点,且飞行方向与各自天线法线方向垂直,在发射端由多个发射天线发射相互正交的信号产生多个发射通道,在接收端用多个天线接收目标的回波信号,雷达杂波谱位于发射空间频率‑接收空间频率‑多普勒频率的三维空间的一个平面上;(2)构造多通道联合降维方法的多普勒域滤波器G<sub>K</sub>和MNK×MNR<sub>K</sub>阶多通道联合降维矩阵T:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>G</mi><mi>K</mi></msub><mo>=</mo><mo>[</mo><msub><mi>a</mi><mi>d</mi></msub><mrow><mo>(</mo><msub><mi>f</mi><mrow><mi>d</mi><mo>,</mo><mi>T</mi><mo>-</mo><msub><mover><mi>r</mi><mo>~</mo></mover><mi>K</mi></msub></mrow></msub><mo>)</mo></mrow><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><msub><mi>a</mi><mi>d</mi></msub><mrow><mo>(</mo><msub><mi>f</mi><mrow><mi>d</mi><mo>,</mo><mi>T</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow><mo>,</mo><msub><mi>a</mi><mi>d</mi></msub><mrow><mo>(</mo><msub><mi>f</mi><mrow><mi>d</mi><mo>,</mo><mi>T</mi></mrow></msub><mo>)</mo></mrow><mo>,</mo><msub><mi>a</mi><mi>d</mi></msub><mrow><mo>(</mo><msub><mi>f</mi><mrow><mi>d</mi><mo>,</mo><mi>T</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><msub><mi>a</mi><mi>d</mi></msub><mrow><mo>(</mo><msub><mi>f</mi><mrow><mi>d</mi><mo>,</mo><mi>T</mi><mo>+</mo><msub><mover><mi>r</mi><mo>~</mo></mover><mi>K</mi></msub></mrow></msub><mo>)</mo></mrow><mo>]</mo><mo>,</mo></mrow>]]></math><img file="FDA0000587234500000013.GIF" wi="1367" he="101" /></maths><maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mi>T</mi><mo>=</mo><msub><mi>I</mi><mi>MN</mi></msub><mo>&CircleTimes;</mo><msub><mi>G</mi><mi>K</mi></msub><mo>,</mo></mrow>]]></math><img file="FDA0000587234500000019.GIF" wi="308" he="67" /></maths>其中,R<sub>K</sub>为构造降维矩阵时选取的多普勒通道的数目,I<sub>MN</sub>为MN阶的单位矩阵,<img file="FDA0000587234500000014.GIF" wi="517" he="79" />为目标多普勒通道附近第<img file="FDA0000587234500000015.GIF" wi="56" he="75" />个通道的多普勒导向矢量,<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mover><mi>r</mi><mo>~</mo></mover><mi>K</mi></msub><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><mrow><mo>(</mo><msub><mi>R</mi><mi>K</mi></msub><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>/</mo><mn>2</mn><mo>;</mo></mrow>]]></math><img file="FDA0000587234500000021.GIF" wi="482" he="77" /></maths>(3)用上述降维矩阵T乘以空时数据矢量y和目标导向矢量b<sub>T</sub>(f<sub>t,T</sub>,f<sub>r,T</sub>,f<sub>d,T</sub>),得到降维后的数据矢量z和降维后的目标导向矢量c<sub>T</sub>(f<sub>t,T,</sub>f<sub>r,T</sub>,f<sub>d,T</sub>):z=T<sup>H</sup>y,c<sub>T</sub>(f<sub>t,T</sub>,f<sub>r,T</sub>,f<sub>d,T</sub>)=T<sup>H</sup>b<sub>T</sub>(f<sub>t,T</sub>,f<sub>r,T</sub>,f<sub>d,T</sub>),其中,上标<sup>H</sup>表示共轭转置;(4)利用一个距离门的数据矢量z计算回波协方差矩阵R:R=zz<sup>H</sup>,并对2×MNR<sub>K</sub>个距离门的协方差矩阵取平均,得到其估计值<img file="FDA0000587234500000022.GIF" wi="76" he="76" /><maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><mover><mi>R</mi><mo>^</mo></mover><mo>=</mo><mfrac><mn>1</mn><mrow><mn>2</mn><mo>&times;</mo><msub><mi>MNR</mi><mi>K</mi></msub></mrow></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>0</mn></mrow><mrow><mn>2</mn><mo>&times;</mo><msub><mi>MNR</mi><mi>K</mi></msub><mo>-</mo><mn>1</mn></mrow></munderover><msub><mi>z</mi><mi>k</mi></msub><msubsup><mi>z</mi><mi>k</mi><mi>H</mi></msubsup><mo>,</mo></mrow>]]></math><img file="FDA0000587234500000023.GIF" wi="587" he="154" /></maths>其中,z<sub>k</sub>表示第k个距离门的降维后的数据矢量,距离门是指以接收机和发射机在地面的投影点为焦点,以地面杂波点到两个双基地的距离之和为固定值的一组椭圆环;(5)根据空时自适应处理原理,由上述回波协方差矩阵估计值<img file="FDA0000587234500000024.GIF" wi="50" he="70" />得到最优权矢量:<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><mi>w</mi><mo>=</mo><mi>&mu;</mi><msup><mover><mi>R</mi><mo>^</mo></mover><mrow><mo>-</mo><mn>1</mn></mrow></msup><msub><mi>c</mi><mi>T</mi></msub><mrow><mo>(</mo><msub><mi>f</mi><mrow><mi>t</mi><mo>,</mo><mi>T</mi></mrow></msub><mo>,</mo><msub><mi>f</mi><mrow><mi>r</mi><mo>,</mo><mi>T</mi></mrow></msub><mo>,</mo><msub><mi>f</mi><mrow><mi>d</mi><mo>,</mo><mi>T</mi></mrow></msub><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000587234500000025.GIF" wi="582" he="94" /></maths>其中μ为一个标量,<img file="FDA0000587234500000026.GIF" wi="87" he="76" />为回波协方差矩阵估计值<img file="FDA0000587234500000027.GIF" wi="52" he="77" />的逆矩阵;(6)利用上述最优权矢量w对降维后的数据矢量z进行加权,得到针对目标位置抑制杂波的回波数据,检测出目标。
地址 710071 陕西省西安市太白南路2号