发明名称 针对雷达主瓣干扰的极化抑制处理方法
摘要 本发明提供一种针对雷达主瓣干扰的极化抑制处理方法。技术方案包括五个处理步骤:第一步,建立关于雷达天线极化特性的函数库;第二步,干扰和目标混合信号的预处理;第三步,干扰和目标混合信号的正交极化分解处理;第四步,干扰和目标混合信号的极化滤波矢量估计;第五步,干扰和目标混合信号的空域多凹口极化滤波。本发明利用最小二乘算法对干扰和目标混合信号的电压矩阵进行处理,实现了对干扰和目标混合信号的正交极化分解与估计;在不能精确获知俯仰角的情况下,利用空域多凹口虚拟极化滤波器抑制干扰信号。本发明对硬件系统无额外要求,易于实现,工程适用性强。
申请公布号 CN102495394A 申请公布日期 2012.06.13
申请号 CN201110400197.4 申请日期 2011.12.06
申请人 中国人民解放军国防科学技术大学 发明人 王雪松;李永祯;戴幻尧;常宇亮;董健;代大海;李棉全
分类号 G01S7/36(2006.01)I 主分类号 G01S7/36(2006.01)I
代理机构 国防科技大学专利服务中心 43202 代理人 王文惠
主权项 1.一种针对雷达主瓣干扰的极化抑制处理方法,应用于单极化的机械扫描体制雷达,其特征在于,包括下述步骤:第一步、建立关于雷达天线极化特性的函数库;已知雷达天线主极化方式为水平极化H,交叉极化方式为垂直极化V;用现有的天线特性测量方法测量在不同方位角<img file="FDA0000116449360000011.GIF" wi="30" he="35" />和俯仰角θ下的雷达天线增益<img file="FDA0000116449360000012.GIF" wi="445" he="136" />其中<img file="FDA0000116449360000013.GIF" wi="165" he="50" />表示在方位角<img file="FDA0000116449360000014.GIF" wi="30" he="35" />和俯仰角θ下雷达天线的水平极化增益,<img file="FDA0000116449360000015.GIF" wi="161" he="50" />表示在方位角<img file="FDA0000116449360000016.GIF" wi="30" he="35" />和仰角θ下雷达天线的垂直极化增益;其中,俯仰角θ的分别取值为θ<sub>1</sub>,θ<sub>2</sub>,L,θ<sub>N</sub>,俯仰角θ的取值是根据雷达天线的俯仰波束宽度和测量精度决定;俯仰波束宽度越宽,俯仰角θ的取值数越多;测量精度要求越高,俯仰角θ的取值数越多;方位角<img file="FDA0000116449360000017.GIF" wi="30" he="35" />在雷达天线方位扫描区间<img file="FDA0000116449360000018.GIF" wi="192" he="61" />内取值,即方位角<img file="FDA0000116449360000019.GIF" wi="667" he="52" /><img file="FDA00001164493600000110.GIF" wi="223" he="38" />方位角<img file="FDA00001164493600000111.GIF" wi="30" he="35" />的取值根据测量精度的需要决定;第二步、干扰和目标混合信号的预处理;雷达天线在工作过程中俯仰角为θ′,接收的干扰和目标混合信号的电压值<img file="FDA00001164493600000112.GIF" wi="158" he="61" />是采样时刻t和方位角<img file="FDA00001164493600000113.GIF" wi="30" he="35" />的函数,则在采样时刻t=t<sub>1</sub>,t<sub>2,L</sub>,t<sub>M</sub>和方位角<img file="FDA00001164493600000114.GIF" wi="420" he="48" />下接收的干扰和目标混合信号的电压矩阵V<sub>θ′</sub>为:<img file="FDA00001164493600000115.GIF" wi="1263" he="286" />第三步、干扰和目标混合信号的正交极化分解处理;通过第一步测量可知在俯仰角θ=θ<sub>1</sub>,θ<sub>2</sub>,L,θ<sub>N</sub>下,雷达天线极化特性矩阵G<sub>θ</sub>为:<img file="FDA00001164493600000116.GIF" wi="1168" he="136" />由电压矩阵V<sub>θ′</sub>可知在采样时刻t的电压列矢量<img file="FDA00001164493600000117.GIF" wi="118" he="49" /><img file="FDA00001164493600000118.GIF" wi="779" he="69" />则利用下式计算在采样时刻t=t<sub>1</sub>,t<sub>2,L</sub>,t<sub>M</sub>和俯仰角θ=θ<sub>1</sub>,θ<sub>2</sub>,L,θ<sub>N</sub>下正交极化分量R<sub>t,θ</sub>:<maths num="0001"><![CDATA[<math><mrow><msub><mi>R</mi><mrow><mi>t</mi><mo>,</mo><mi>&theta;</mi></mrow></msub><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>r</mi><msub><mi>H</mi><mrow><mi>t</mi><mo>,</mo><mi>&theta;</mi></mrow></msub></msub></mtd></mtr><mtr><mtd><msub><mi>r</mi><msub><mi>V</mi><mrow><mi>t</mi><mo>,</mo><mi>&theta;</mi></mrow></msub></msub></mtd></mtr></mtable></mfenced><mo>=</mo><msup><mfenced open='[' close=']'><mtable><mtr><mtd><msup><msub><mi>G</mi><mi>&theta;</mi></msub><mi>H</mi></msup></mtd><mtd><msub><mi>G</mi><mi>&theta;</mi></msub></mtd></mtr></mtable></mfenced><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>&CenterDot;</mo><msup><msub><mi>G</mi><mi>&theta;</mi></msub><mi>H</mi></msup><mo>&CenterDot;</mo><msub><mi>V</mi><mi>t</mi></msub></mrow></math>]]></maths>(公式一)上式中,G<sub>θ</sub><sup>H</sup>表示矩阵G<sub>θ</sub>的共轭转置;将采样时刻t=t<sub>1</sub>,t<sub>2,L</sub>,t<sub>M</sub>代入公式一即可得到关于采样时刻的正交极化信号序列R<sub>θ</sub>:<maths num="0002"><![CDATA[<math><mrow><msub><mi>R</mi><mi>&theta;</mi></msub><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>R</mi><mrow><msub><mi>t</mi><mn>1</mn></msub><mo>,</mo><mi>&theta;</mi></mrow></msub><mo>,</mo></mtd><mtd><msub><mi>R</mi><mrow><msub><mi>t</mi><mn>2</mn></msub><mo>,</mo><mi>&theta;</mi></mrow></msub><mo>,</mo></mtd><mtd><mi>L</mi></mtd><mtd><msub><mi>R</mi><mrow><msub><mi>t</mi><mi>M</mi></msub><mo>,</mo><mi>&theta;</mi></mrow></msub></mtd></mtr></mtable></mfenced><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>r</mi><msub><mi>H</mi><mi>&theta;</mi></msub></msub></mtd></mtr><mtr><mtd><msub><mi>r</mi><msub><mi>V</mi><mi>&theta;</mi></msub></msub></mtd></mtr></mtable></mfenced><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>r</mi><msub><mi>H</mi><mrow><msub><mi>t</mi><mn>1</mn></msub><mo>,</mo><mi>&theta;</mi></mrow></msub></msub><mo>,</mo></mtd><mtd><msub><mi>r</mi><msub><mi>H</mi><mrow><msub><mi>t</mi><mn>2</mn></msub><mo>,</mo><mi>&theta;</mi></mrow></msub></msub><mo>,</mo></mtd><mtd><mi>L</mi><mo>,</mo></mtd><mtd><msub><mi>r</mi><msub><mi>H</mi><mrow><msub><mi>t</mi><mi>M</mi></msub><mo>,</mo><mi>&theta;</mi></mrow></msub></msub></mtd></mtr><mtr><mtd><msub><mi>r</mi><msub><mi>V</mi><mrow><msub><mi>t</mi><mn>1</mn></msub><mo>,</mo><mi>&theta;</mi></mrow></msub></msub><mo>,</mo></mtd><mtd><msub><mi>r</mi><msub><mi>V</mi><mrow><msub><mi>t</mi><mn>2</mn></msub><mo>,</mo><mi>&theta;</mi></mrow></msub></msub><mo>,</mo></mtd><mtd><mi>L</mi><mo>,</mo></mtd><mtd><msub><mi>r</mi><msub><mi>V</mi><mrow><msub><mi>t</mi><mi>M</mi></msub><mo>,</mo><mi>&theta;</mi></mrow></msub></msub></mtd></mtr></mtable></mfenced></mrow></math>]]></maths>(公式二)将俯仰角θ=θ<sub>1</sub>,θ<sub>2</sub>,L,θ<sub>N</sub>代入公式二,可得到N路的正交极化信号<maths num="0003"><![CDATA[<math><mrow><msub><mi>R</mi><msub><mi>&theta;</mi><mn>1</mn></msub></msub><mo>,</mo><msub><mi>R</mi><msub><mi>&theta;</mi><mn>2</mn></msub></msub><mo>,</mo><mi>L</mi><mo>,</mo><msub><mi>R</mi><msub><mi>&theta;</mi><mi>N</mi></msub></msub><mo>;</mo></mrow></math>]]></maths>第四步、干扰和目标混合信号的极化滤波矢量估计;在俯仰角θ=θ<sub>1</sub>,θ<sub>2</sub>,L,θ<sub>N</sub>下,用下式求均值,估计出信号的极化状态<img file="FDA0000116449360000024.GIF" wi="84" he="64" /><maths num="0004"><![CDATA[<math><mrow><msub><mover><mi>J</mi><mo>^</mo></mover><mi>&theta;</mi></msub><mo>=</mo><mfrac><mn>1</mn><mi>M</mi></mfrac><mfenced open='[' close=']'><mtable><mtr><mtd><munderover><mi>&Sigma;</mi><mrow><mi>m</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mi>r</mi><msub><mi>H</mi><mrow><msub><mi>t</mi><mi>m</mi></msub><mo>,</mo><mi>&theta;</mi></mrow></msub></msub></mtd></mtr><mtr><mtd><munderover><mi>&Sigma;</mi><mrow><mi>m</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mi>r</mi><msub><mi>V</mi><mrow><msub><mi>t</mi><mi>m</mi></msub><mo>,</mo><mi>&theta;</mi></mrow></msub></msub></mtd></mtr></mtable></mfenced><mo>,</mo></mrow></math>]]></maths>求出满足下面表达式的极化滤波矢量H<sub>q</sub>:<maths num="0005"><![CDATA[<math><mfenced open='{' close=''><mtable><mtr><mtd><msubsup><mi>H</mi><mi>&theta;</mi><mi>T</mi></msubsup><msub><mover><mi>J</mi><mo>^</mo></mover><mi>&theta;</mi></msub><mo>=</mo><mn>0</mn></mtd></mtr><mtr><mtd><msup><mrow><mo>|</mo><mo>|</mo><msub><mi>H</mi><mi>&theta;</mi></msub><mo>|</mo><mo>|</mo></mrow><mn>2</mn></msup><mo>=</mo><mn>1</mn></mtd></mtr></mtable></mfenced></math>]]></maths>由俯仰角θ的取值为θ<sub>1</sub>,θ<sub>2</sub>,L,θ<sub>N</sub>,可以得到由N个极化滤波矢量H<sub>q</sub>组成的极化滤波器<maths num="0006"><![CDATA[<math><mrow><msub><mi>H</mi><msub><mi>&theta;</mi><mn>1</mn></msub></msub><mo>,</mo><msub><mi>H</mi><msub><mi>&theta;</mi><mn>2</mn></msub></msub><mo>,</mo><mi>L</mi><mo>,</mo><msub><mi>H</mi><msub><mi>&theta;</mi><mi>N</mi></msub></msub><mo>;</mo></mrow></math>]]></maths>第五步、干扰和目标混合信号的空域多凹口极化滤波;利用下式得到对干扰和目标混合信号进行抑制干扰后的输出信号O<sub>f</sub>:<maths num="0007"><![CDATA[<math><mrow><msub><mi>O</mi><mi>f</mi></msub><mo>=</mo><munderover><mi>min</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mo>{</mo><msup><msub><mi>H</mi><msub><mi>&theta;</mi><mi>k</mi></msub></msub><mi>T</mi></msup><mo>&CenterDot;</mo><msub><mi>R</mi><msub><mi>&theta;</mi><mi>k</mi></msub></msub><mo>}</mo><mo>.</mo></mrow></math>]]></maths>
地址 410073 湖南省长沙市开福区德雅路109号