发明名称 基于波束保形的阵列雷达自适应抗干扰方法
摘要 本发明公开了一种基于波束保形的阵列雷达自适应抗干扰方法,主要解决现有技术形成的波束主瓣不平坦以及不适用于多点约束的问题。其实现步骤是:1.采集阵列雷达检测目标的采样数据;2.设计优化的三个约束点的导向矢量;3.根据步骤2中的导向矢量构建数学模型;4.设计数学模型的约束矢量,并根据该约束矢量求解数学模型的自适应权矢量;5.用自适应权矢量对干扰信号进行抑制,得到输出数据。本发明通过在主波束多个点内采用线性约束,通过设置相位响应的值实现对目标信号相位的匹配,避免了干扰对检测性能的影响;同时通过设置主波束幅度响应为1,使主波束能均匀覆盖观测区域,提高了阵列雷达的目标检测性能,可用于多目标的检测。
申请公布号 CN105929374A 申请公布日期 2016.09.07
申请号 CN201610471771.8 申请日期 2016.06.24
申请人 西安电子科技大学 发明人 许京伟;廖桂生;冯阳;王成浩;兰岚
分类号 G01S7/36(2006.01)I;G01S7/41(2006.01)I 主分类号 G01S7/36(2006.01)I
代理机构 陕西电子工业专利中心 61205 代理人 王品华;朱红星
主权项 基于波束保形的阵列雷达自适应抗干扰方法,包括:1)在阵列信号处理机中输入阵列雷达检测目标的采样数据X,该数据X包含目标和干扰信号;2)设计优化的三个约束点的导向矢量:2a)设计三个约束点的角度,即第一个约束点的角度θ<sub>1</sub>是估计的目标到达角,第二个约束点的角度θ<sub>2</sub>=θ<sub>1</sub>‑3°,第三个约束点的角度θ<sub>3</sub>=θ<sub>1</sub>+3°;2b)根据步骤2a)中三个约束点的角度θ<sub>1</sub>、θ<sub>2</sub>、θ<sub>3</sub>,计算优化的三个约束点的导向矢量分别为a(θ<sub>1</sub>),a(θ<sub>2</sub>),a(θ<sub>3</sub>);3)根据步骤2b)中三个约束点的导向矢量a(θ<sub>1</sub>),a(θ<sub>2</sub>),a(θ<sub>3</sub>)构建数学模型,即<img file="FDA0001028254690000011.GIF" wi="253" he="87" />使得w<sup>H</sup>C=u<sup>T</sup>,其中,w为待求的自适应权矢量,(·)<sup>H</sup>表示共轭转置操作,(·)<sup>T</sup>表示转置操作,<img file="FDA0001028254690000012.GIF" wi="60" he="71" />为采样数据的协方差矩阵,C=[a(θ<sub>1</sub>)a(θ<sub>2</sub>)a(θ<sub>3</sub>)]为约束矩阵,<img file="FDA0001028254690000013.GIF" wi="412" he="93" />为约束矢量,<img file="FDA0001028254690000014.GIF" wi="158" he="47" />为待求的两个不同约束相位,j表示虚数;4)设计约束矢量:4a)根据步骤3)中的约束矩阵C和采样协方差矩阵<img file="FDA0001028254690000015.GIF" wi="62" he="74" />计算得到一个中间矩阵<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msup><mrow><mo>(</mo><msup><mi>C</mi><mi>H</mi></msup><msup><msub><mover><mi>R</mi><mo>^</mo></mover><mi>x</mi></msub><mrow><mo>-</mo><mn>1</mn></mrow></msup><mi>C</mi><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mover><mi>a</mi><mo>~</mo></mover><mn>11</mn></msub></mtd><mtd><mrow><msub><mover><mi>a</mi><mo>~</mo></mover><mn>12</mn></msub><msup><mi>e</mi><mrow><mi>j</mi><msub><mover><mi>b</mi><mo>~</mo></mover><mn>12</mn></msub></mrow></msup></mrow></mtd><mtd><mrow><msub><mover><mi>a</mi><mo>~</mo></mover><mn>13</mn></msub><msup><mi>e</mi><mrow><mi>j</mi><msub><mover><mi>b</mi><mo>~</mo></mover><mn>13</mn></msub></mrow></msup></mrow></mtd></mtr><mtr><mtd><mrow><msub><mover><mi>a</mi><mo>~</mo></mover><mn>12</mn></msub><msup><mi>e</mi><mrow><mo>-</mo><mi>j</mi><msub><mover><mi>b</mi><mo>~</mo></mover><mn>12</mn></msub></mrow></msup></mrow></mtd><mtd><msub><mover><mi>a</mi><mo>~</mo></mover><mn>22</mn></msub></mtd><mtd><mrow><msub><mover><mi>a</mi><mo>~</mo></mover><mn>23</mn></msub><msup><mi>e</mi><mrow><mi>j</mi><msub><mover><mi>b</mi><mo>~</mo></mover><mn>23</mn></msub></mrow></msup></mrow></mtd></mtr><mtr><mtd><mrow><msub><mover><mi>a</mi><mo>~</mo></mover><mn>13</mn></msub><msup><mi>e</mi><mrow><mo>-</mo><mi>j</mi><msub><mover><mi>b</mi><mo>~</mo></mover><mn>13</mn></msub></mrow></msup></mrow></mtd><mtd><mrow><msub><mover><mi>a</mi><mo>~</mo></mover><mn>23</mn></msub><msup><mi>e</mi><mrow><mo>-</mo><mi>j</mi><msub><mover><mi>b</mi><mo>~</mo></mover><mn>23</mn></msub></mrow></msup></mrow></mtd><mtd><msub><mover><mi>a</mi><mo>~</mo></mover><mn>33</mn></msub></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0001028254690000016.GIF" wi="878" he="270" /></maths>其中,(·)<sup>‑1</sup>表示求逆操作,<img file="FDA0001028254690000017.GIF" wi="998" he="71" />为上述中间矩阵中的九个数值不同的中间数;4b)根据步骤4a)中的第七个中间数<img file="FDA0001028254690000018.GIF" wi="51" he="71" />和第九个中间数<img file="FDA0001028254690000019.GIF" wi="59" he="75" />得到第二个约束点的约束相位<img file="FDA00010282546900000110.GIF" wi="179" he="71" />和第三个约束点的约束相位<img file="FDA00010282546900000111.GIF" wi="515" he="71" />4c)根据步骤4b)中的第二个约束点的约束相位<img file="FDA0001028254690000021.GIF" wi="51" he="46" />和第三个约束点的约束相位<img file="FDA0001028254690000022.GIF" wi="51" he="47" />得到约束矢量为<img file="FDA0001028254690000023.GIF" wi="611" he="118" />5)求解自适应权矢量:根据步骤3)中的约束矩阵C以及步骤4c)中的约束矢量u得到自适应权矢量为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>w</mi><mo>=</mo><msup><msub><mover><mi>R</mi><mo>^</mo></mover><mi>x</mi></msub><mrow><mo>-</mo><mn>1</mn></mrow></msup><mi>C</mi><msup><mrow><mo>(</mo><msup><mi>C</mi><mi>H</mi></msup><msup><msub><mover><mi>R</mi><mo>^</mo></mover><mi>x</mi></msub><mrow><mo>-</mo><mn>1</mn></mrow></msup><mi>C</mi><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><mi>u</mi><mo>;</mo></mrow>]]></math><img file="FDA0001028254690000024.GIF" wi="542" he="99" /></maths>6)根据步骤5)中的自适应权矢量w对干扰信号进行抑制,得到抑制干扰后的输出信号Y=w<sup>T</sup>X,其中,(·)<sup>T</sup>表示转置操作。
地址 710071 陕西省西安市太白南路2号