发明名称 前视阵列雷达的俯仰-慢时间空时二维自适应处理方法
摘要 发明公开了一种前视阵列雷达的俯仰-慢时间空时二维自适应处理方法,主要解决现有技术对强杂波抑制效果差的问题。其实现过程是:1)把整个平面阵列划分成若干行子阵,经过行合成等效得到一个沿俯仰方向的等距线阵;2)对等效线阵的回波数据进行采样;3)将处理器接收的单个距离单元空时采样数据排成矢量形式;4)计算该距离单元的多普勒补偿矩阵,对该距离单元的数据矢量进行多普勒频移补偿;5)用多普勒频移补偿后的数据估计待检测单元杂波协方差矩阵,计算空时二维自适应处理最优权,得到处理后系统的输出数据。本发明滤波后获得的改善因子优于传统方位—慢时间空时二维自适应处理,可用于抑制前视阵列雷达的强杂波。
申请公布号 CN102798842B 申请公布日期 2014.03.12
申请号 CN201210245519.7 申请日期 2012.07.16
申请人 西安电子科技大学 发明人 刘峥;陈熠;谢荣;张磊;严湘南;张伦
分类号 G01S7/36(2006.01)I 主分类号 G01S7/36(2006.01)I
代理机构 陕西电子工业专利中心 61205 代理人 王品华;朱红星
主权项 1.一种前视阵列雷达的俯仰-慢时间空时二维自适应处理方法,包括如下步骤:(1)在机载前视平面相控阵雷达中,设计天线阵面使其与水平地面有一倾角α,0<α<90°,阵面经过行合成后形成沿俯仰方向的N元等距离线阵接收地面回波信号;(2)对上述阵列接收的地面回波信号进行采样,获得第p个阵元第q个脉冲在第l距离单元的采样数据x<sub>pq</sub>(l),其中,p=1,2,…N,N为俯仰方向阵元数,q=1,2,…K,K为脉冲数;(3)将上述采样数据x<sub>pq</sub>(l)排列成NK×1维空时二维采样数据矢量X(l);(4)对空时二维采样数据矢量X(l)进行多普勒频移补偿,补偿后的数据矢量Y(l)为:Y(l)=T<sub>l</sub><sup>H</sup>X(l),其中[]<sup>H</sup>为共轭转置操作,T<sub>l</sub>为补偿矩阵:<img file="FDA0000380904960000011.GIF" wi="682" he="305" />式中,<img file="FDA0000380904960000017.GIF" wi="169" he="64" />为虚数单元符号,I<sub>N</sub>是N×N的单位矩阵,<img file="FDA0000380904960000016.GIF" wi="56" he="51" />为Kronecker积,Δ为中间变量,<img file="FDA0000380904960000012.GIF" wi="642" he="141" />l<sub>0</sub>为待检测单元距离单元数,<img file="FDA0000380904960000013.GIF" wi="64" he="67" />为待检测单元的俯仰角,<img file="FDA0000380904960000014.GIF" wi="54" he="54" />为第l采样距离单元的俯仰角,v为载机速度,λ为波长,f<sub>r</sub>为脉冲重复频率;(5)用上述多普勒频移补偿后的数据Y(l)进行空时二维自适应处理,求得空时二维自适应处理的最优权值W<sub>opt</sub>:<![CDATA[<math><mrow><msub><mi>W</mi><mi>opt</mi></msub><mo>=</mo><mfrac><mrow><msup><mover><msub><mi>R</mi><mn>0</mn></msub><mo>^</mo></mover><mrow><mo>-</mo><mn>1</mn></mrow></msup><mi>S</mi><mrow><mo>(</mo><msub><mi>&psi;</mi><mn>0</mn></msub><mo>,</mo><msub><mi>f</mi><mi>d</mi></msub><mo>)</mo></mrow></mrow><mrow><msup><mi>S</mi><mi>H</mi></msup><mrow><mo>(</mo><msub><mi>&psi;</mi><mn>0</mn></msub><mo>,</mo><msub><mi>f</mi><mi>d</mi></msub><mo>)</mo></mrow><msup><mover><msub><mi>R</mi><mn>0</mn></msub><mo>^</mo></mover><mrow><mo>-</mo><mn>1</mn></mrow></msup><mi>S</mi><mrow><mo>(</mo><msub><mi>&psi;</mi><mn>0</mn></msub><mo>,</mo><msub><mi>f</mi><mi>d</mi></msub><mo>)</mo></mrow></mrow></mfrac><mo>,</mo></mrow></math>]]></maths>式中,<img file="FDA0000380904960000022.GIF" wi="68" he="82" />为待检测单元的杂波协方差矩阵,[]<sup>-1</sup>为矩阵求逆运算,ψ<sub>0</sub>为目标方向锥角,f<sub>d</sub>为目标多普勒频率,S(ψ<sub>0</sub>,f<sub>d</sub>)为目标导向矢量:<![CDATA[<math><mrow><mi>S</mi><mrow><mo>(</mo><msub><mi>&psi;</mi><mn>0</mn></msub><mo>,</mo><msub><mi>f</mi><mi>d</mi></msub><mo>)</mo></mrow><mo>=</mo><msub><mi>S</mi><mi>s</mi></msub><mrow><mo>(</mo><msub><mi>&psi;</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>&CircleTimes;</mo><msub><mi>S</mi><mi>t</mi></msub><mrow><mo>(</mo><msub><mi>f</mi><mi>d</mi></msub><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths>S<sub>s</sub>(ψ<sub>0</sub>)和S<sub>t</sub>(f<sub>d</sub>)分别为目标信号的空间导向矢量和时间导向矢量:<![CDATA[<math><mrow><msub><mi>S</mi><mi>s</mi></msub><mrow><mo>(</mo><msub><mi>&psi;</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>=</mo><msup><mfenced open='[' close=']'><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mi>exp</mi><mrow><mo>(</mo><mi>j</mi><mn>2</mn><mi>&pi;</mi><mfrac><mi>d</mi><mi>&lambda;</mi></mfrac><mi>cos</mi><msub><mi>&psi;</mi><mn>0</mn></msub><mo>)</mo></mrow></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><mi>exp</mi><mrow><mo>(</mo><mi>j</mi><mn>2</mn><mi>&pi;</mi><mfrac><mi>d</mi><mi>&lambda;</mi></mfrac><mrow><mo>(</mo><mi>N</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mi>cos</mi><msub><mi>&psi;</mi><mn>0</mn></msub><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mi>T</mi></msup><mo>,</mo></mrow></math>]]></maths><![CDATA[<math><mrow><msub><mi>S</mi><mi>t</mi></msub><mrow><mo>(</mo><msub><mi>f</mi><mi>d</mi></msub><mo>)</mo></mrow><mo>=</mo><msup><mfenced open='[' close=']'><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mi>exp</mi><mrow><mo>(</mo><mi>j</mi><mn>2</mn><mi>&pi;</mi><msub><mi>f</mi><mi>d</mi></msub><mo>)</mo></mrow></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><mi>exp</mi><mrow><mo>(</mo><mi>j</mi><mn>2</mn><mi>&pi;</mi><mfrac><mi>d</mi><mi>&lambda;</mi></mfrac><mrow><mo>(</mo><mi>K</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><msub><mi>f</mi><mi>d</mi></msub><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mi>T</mi></msup><mo>,</mo></mrow></math>]]></maths><img file="FDA0000380904960000028.GIF" wi="48" he="49" />为Kronecker积,N为俯仰方向阵元数,K为脉冲数,λ为波长,d为阵元间距,[]<sup>T</sup>为转置操作;(6)用上述获得的最优权W<sub>opt</sub>滤除待检测单元的杂波,得到空时二维自适应处理后的输出数据Z<sub>l0</sub>:<![CDATA[<math><mrow><msub><mi>Z</mi><mrow><mi>l</mi><mn>0</mn></mrow></msub><mo>=</mo><msubsup><mi>W</mi><mi>opt</mi><mi>H</mi></msubsup><mi>X</mi><mrow><mo>(</mo><msub><mi>l</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths>其中X(l<sub>0</sub>)为待检测单元的数据矢量。
地址 710071 陕西省西安市太白南路2号