发明名称 任意阵列流形的单基地MIMO雷达波达方向估计方法
摘要 本发明公开了一种基于任意阵列流形的单基地MIMO雷达波达方向估计方法。主要解决现有技术只适用于线性阵列和计算量大的问题。其实现步骤是:1)根据阵列流形写出单基地MIMO雷达的导向矢量;2)对MIMO雷达的导向矢量流形展开,获得范得蒙导向矢量,采用空间平滑思想进行导向矢量变换,求出从MIMO导向矢量到变换后导向矢量的转换矩阵;3)利用接收阵列与发射波形进行匹配滤波,形成自相关矩阵;4)对自相关矩阵进行特征分解,得到特征值和特征向量,选取特征向量形成噪声子空间;5)利用噪声子空间形成空间零谱函数,采用多项式求根方法得到方位角。本发明可实现对任意阵列流形的单基地MIMO雷达波达方向快速估计,且运算量小,可用于目标定位和雷达追踪。
申请公布号 CN103364762A 申请公布日期 2013.10.23
申请号 CN201310296467.0 申请日期 2013.07.13
申请人 西安电子科技大学 发明人 曹运合;陈婷;姜婕;王敏;王胜华;彭志刚
分类号 G01S7/02(2006.01)I 主分类号 G01S7/02(2006.01)I
代理机构 陕西电子工业专利中心 61205 代理人 王品华;朱红星
主权项 1.一种任意阵列流形的单基地MIMO雷达波达方向估计方法,包括如下步骤:1)根据阵列流形写出单基地MIMO雷达的导向矢量A(θ),其中θ为方位角;2)对MIMO雷达的导向矢量A(θ)流形展开,获得范得蒙导向矢量D(θ),对D(θ)进行空间平滑处理,处理后的导向矢量记为D<sub>S</sub>(θ),利用最小二乘法得出从MIMO雷达导向矢量A(θ)到平滑处理后的导向矢量D<sub>S</sub>(θ)的转换矩阵G;3)对雷达的接收数据与发射波形进行匹配滤波,将匹配滤波结果记为X(t),再利用滤波结果X(t)形成自相关矩阵R<sub>x</sub>;4)对自相关矩阵R<sub>x</sub>进行特征分解,得到一系列特征值λ<sub>m</sub>及其对应的特征向量e<sub>m</sub>,m=1,2,…,N<sup>2</sup>,其中N为阵元数;将特征值λ<sub>m</sub>进行从小到大排序,依次选取N<sup>2</sup>-P个小特征值,P为同一距离门内的目标个数,取这些特征值对应的特征向量<img file="FDA00003507573800017.GIF" wi="355" he="64" />形成噪声子空间<maths num="0001"><![CDATA[<math><mrow><msub><mover><mi>E</mi><mo>^</mo></mover><mi>n</mi></msub><mo>=</mo><mo>[</mo><msub><mi>e</mi><mrow><mi>P</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><msub><mi>e</mi><mi>m</mi></msub><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><msub><mi>e</mi><msup><mi>N</mi><mn>2</mn></msup></msub><mo>]</mo><mo>;</mo></mrow></math>]]></maths>5)利用噪声子空间<img file="FDA00003507573800012.GIF" wi="71" he="86" />形成MUSIC空间零谱函数:<maths num="0002"><![CDATA[<math><mrow><mi>f</mi><mrow><mo>(</mo><mi>&theta;</mi><mo>)</mo></mrow><mo>=</mo><msup><mrow><mo>|</mo><mo>|</mo><msup><msub><mover><mi>E</mi><mo>^</mo></mover><mi>n</mi></msub><mi>H</mi></msup><mi>A</mi><mrow><mo>(</mo><mi>&theta;</mi><mo>)</mo></mrow><mo>|</mo><mo>|</mo></mrow><mn>2</mn></msup><mo>=</mo><msup><mi>A</mi><mi>H</mi></msup><mrow><mo>(</mo><mi>&theta;</mi><mo>)</mo></mrow><msub><mover><mi>E</mi><mo>^</mo></mover><mi>n</mi></msub><msup><msub><mover><mi>E</mi><mo>^</mo></mover><mi>n</mi></msub><mi>H</mi></msup><mi>A</mi><mrow><mo>(</mo><mi>&theta;</mi><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths>其中|| ||<sup>2</sup>表示2-范数,<img file="FDA00003507573800014.GIF" wi="90" he="85" />表示矩阵<img file="FDA00003507573800015.GIF" wi="64" he="87" />的共轭转置;6)把步骤5)中的A(θ)用步骤2)中的GD<sub>s</sub>(θ)代替,并令z=e<sup>jθ</sup>,则f(θ)可转化为:<img file="FDA00003507573800016.GIF" wi="662" he="91" />令f(z)=0,采用多项式求根的方法求得P个方位角的复指数形式z<sub>k</sub>,其中k=1,2…P;7)利用接收角的复指数形式z<sub>k</sub>求得第k个接收方位角θ<sub>k</sub>:θ<sub>k</sub>=angle(z<sub>k</sub>),其中angle()表示求相位角。
地址 710071 陕西省西安市太白南路2号