发明名称 基于空时自适应处理雷达的和差波束形成方法
摘要 本发明公开一种基于空时自适应处理雷达的和差波束形成方法,主要解决现有技术主瓣响应失真的问题。其实现步骤是:1.在信号处理机中输入空时自适应处理雷达的采样数据;2.根据采样数据估计得到采样协方差矩阵;3.设计五个约束点的角度和多普勒;4.根据采样协方差矩阵和五个约束点,设计和波束权矢量;5.根据采样协方差矩阵和五个约束点,设计空域差波束权矢量;6.根据采样协方差矩阵和五个约束点,设计时域差波束权矢量;7.根据采样数据和波束权矢量,得到和波束、空域差波束和时域差波束的输出信号。本发明通过幅度相位联合约束、导数约束和多个零点约束,实现了空时自适应处理雷达的和差波束主瓣保形,可用于目标跟踪。
申请公布号 CN106226749A 申请公布日期 2016.12.14
申请号 CN201610512752.5 申请日期 2016.07.01
申请人 西安电子科技大学 发明人 许京伟;廖桂生;王成浩;张玉洪
分类号 G01S7/41(2006.01)I;G01S13/66(2006.01)I 主分类号 G01S7/41(2006.01)I
代理机构 陕西电子工业专利中心 61205 代理人 王品华;朱红星
主权项 基于空时自适应处理STAP雷达的和差波束形成方法,包括:1)在信号处理机中输入空时自适应处理STAP雷达检测目标的采样数据X,该数据X包含目标、杂波和噪声;2)根据采样数据X估计得到采样协方差矩阵<img file="FDA0001038759020000011.GIF" wi="67" he="62" />3)设计五个约束点的角度和多普勒,即五个约束点分别为(θ<sub>i</sub>,f<sub>di</sub>),i=0,1,…,4,其中θ<sub>0</sub>和f<sub>d0</sub>分别为假设目标的角度和多普勒,其余四个约束点的角度分别为θ<sub>1</sub>=θ<sub>0</sub>,θ<sub>2</sub>=θ<sub>0</sub>,θ<sub>3</sub>=θ<sub>0</sub>+3°,θ<sub>4</sub>=θ<sub>0</sub>‑3°,其余多普勒分别为f<sub>d1</sub>=f<sub>d0</sub>‑30Hz,f<sub>d2</sub>=f<sub>d0</sub>+30Hz,f<sub>d3</sub>=f<sub>d0</sub>,f<sub>d4</sub>=f<sub>d0</sub>;4)设计和波束权矢量:4a)根据步骤3)中五个约束点的角度及多普勒(θ<sub>i</sub>,f<sub>di</sub>),计算五个约束点的导向矢量s(θ<sub>i</sub>,f<sub>di</sub>);4b)根据步骤4a)中五个约束点的导向矢量s(θ<sub>i</sub>,f<sub>di</sub>),构建和波束约束矩阵:C<sub>Σ</sub>=[s(θ<sub>0</sub>,f<sub>d0</sub>),s(θ<sub>1</sub>,f<sub>d1</sub>),s(θ<sub>2</sub>,f<sub>d2</sub>),s(θ<sub>3</sub>,f<sub>d3</sub>),s(θ<sub>4</sub>,f<sub>d4</sub>)];4c)根据步骤4a)中第一约束点的导向矢量s(θ<sub>0</sub>,f<sub>d0</sub>)与步骤4b)中和波束约束矩阵C<sub>Σ</sub>,构建和波束响应矢量:<maths num="0001"><math><![CDATA[<mrow><msub><mi>f</mi><mi>&Sigma;</mi></msub><mo>=</mo><msubsup><mi>C</mi><mi>&Sigma;</mi><mi>H</mi></msubsup><mi>s</mi><mrow><mo>(</mo><msub><mi>&theta;</mi><mn>0</mn></msub><mo>,</mo><msub><mi>f</mi><mrow><mi>d</mi><mn>0</mn></mrow></msub><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001038759020000012.GIF" wi="370" he="71" /></maths>其中(·)<sup>H</sup>表示共轭转置操作;4d)根据步骤2)、步骤4b)和步骤4c)的结果,计算和波束权矢量:<maths num="0002"><math><![CDATA[<mrow><msub><mi>w</mi><mi>&Sigma;</mi></msub><mo>=</mo><msup><mover><mi>R</mi><mo>^</mo></mover><mrow><mo>-</mo><mn>1</mn></mrow></msup><msub><mi>C</mi><mi>&Sigma;</mi></msub><msup><mrow><mo>(</mo><msubsup><mi>C</mi><mi>&Sigma;</mi><mi>H</mi></msubsup><msup><mover><mi>R</mi><mo>^</mo></mover><mrow><mo>-</mo><mn>1</mn></mrow></msup><msub><mi>C</mi><mi>&Sigma;</mi></msub><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><msub><mi>f</mi><mi>&Sigma;</mi></msub></mrow>]]></math><img file="FDA0001038759020000013.GIF" wi="558" he="102" /></maths>其中(·)<sup>‑1</sup>表示求逆操作;5)设计空域差波束权矢量:5a)根据步骤4a)中第一约束点的导向矢量s(θ<sub>0</sub>,f<sub>d0</sub>),计算第一约束点导向矢量的角度偏导数s<sub>θ</sub>(θ<sub>0</sub>,f<sub>d0</sub>)和多普勒偏导数<img file="FDA0001038759020000021.GIF" wi="259" he="77" />5b)根据步骤4a)中五个约束点的导向矢量s(θ<sub>i</sub>,f<sub>di</sub>)与步骤5a)中第一约束点导向矢量的角度偏导数s<sub>θ</sub>(θ<sub>0</sub>,f<sub>d0</sub>)和多普勒偏导数<img file="FDA0001038759020000022.GIF" wi="259" he="70" />构建空域差波束约束矩阵:<maths num="0003"><math><![CDATA[<mrow><msub><mi>C</mi><mrow><mi>&Delta;</mi><mi>s</mi></mrow></msub><mo>=</mo><mo>&lsqb;</mo><msub><mi>s</mi><mi>&theta;</mi></msub><mrow><mo>(</mo><msub><mi>&theta;</mi><mn>0</mn></msub><mo>,</mo><msub><mi>f</mi><mrow><mi>d</mi><mn>0</mn></mrow></msub><mo>)</mo></mrow><mo>,</mo><mi>s</mi><mrow><mo>(</mo><msub><mi>&theta;</mi><mn>0</mn></msub><mo>,</mo><msub><mi>f</mi><mrow><mi>d</mi><mn>0</mn></mrow></msub><mo>)</mo></mrow><mo>,</mo><mi>s</mi><mrow><mo>(</mo><msub><mi>&theta;</mi><mn>1</mn></msub><mo>,</mo><msub><mi>f</mi><mrow><mi>d</mi><mn>1</mn></mrow></msub><mo>)</mo></mrow><mo>,</mo><mi>s</mi><mrow><mo>(</mo><msub><mi>&theta;</mi><mn>2</mn></msub><mo>,</mo><msub><mi>f</mi><mrow><mi>d</mi><mn>2</mn></mrow></msub><mo>)</mo></mrow><mo>,</mo><msub><mi>s</mi><msub><mi>f</mi><mi>d</mi></msub></msub><mrow><mo>(</mo><msub><mi>&theta;</mi><mn>0</mn></msub><mo>,</mo><msub><mi>f</mi><mrow><mi>d</mi><mn>0</mn></mrow></msub><mo>)</mo></mrow><mo>&rsqb;</mo><mo>;</mo></mrow>]]></math><img file="FDA0001038759020000023.GIF" wi="1302" he="86" /></maths>5c)设定空域差波束响应矢量为:f<sub>Δs</sub>=[1,0,0,0,0]<sup>T</sup>,其中(·)<sup>T</sup>表示转置操作;5d)根据步骤2)、步骤5b)与步骤5c)的结果,计算空域差波束权矢量:<maths num="0004"><math><![CDATA[<mrow><msub><mi>w</mi><mrow><mi>&Delta;</mi><mi>s</mi></mrow></msub><mo>=</mo><msup><mover><mi>R</mi><mo>^</mo></mover><mrow><mo>-</mo><mn>1</mn></mrow></msup><msub><mi>C</mi><mrow><mi>&Delta;</mi><mi>s</mi></mrow></msub><msup><mrow><mo>(</mo><msubsup><mi>C</mi><mrow><mi>&Delta;</mi><mi>s</mi></mrow><mi>H</mi></msubsup><msup><mover><mi>R</mi><mo>^</mo></mover><mrow><mo>-</mo><mn>1</mn></mrow></msup><msub><mi>C</mi><mrow><mi>&Delta;</mi><mi>s</mi></mrow></msub><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><msub><mi>f</mi><mrow><mi>&Delta;</mi><mi>s</mi></mrow></msub><mo>;</mo></mrow>]]></math><img file="FDA0001038759020000024.GIF" wi="638" he="101" /></maths>6)设计时域差波束权矢量:6a)根据步骤4a)中五个约束点的导向矢量s(θ<sub>i</sub>,f<sub>di</sub>)与步骤5a)中第一约束点导向矢量的空域偏导数s<sub>θ</sub>(θ<sub>0</sub>,f<sub>d0</sub>)和时域偏导数s<sub>fd</sub>(θ<sub>0</sub>,f<sub>d0</sub>),构建时域差波束约束矩阵:<maths num="0005"><math><![CDATA[<mrow><msub><mi>C</mi><mrow><mi>&Delta;</mi><mi>t</mi></mrow></msub><mo>=</mo><mo>&lsqb;</mo><msub><mi>s</mi><msub><mi>f</mi><mi>d</mi></msub></msub><mrow><mo>(</mo><msub><mi>&theta;</mi><mn>0</mn></msub><mo>,</mo><msub><mi>f</mi><mrow><mi>d</mi><mn>0</mn></mrow></msub><mo>)</mo></mrow><mo>,</mo><mi>s</mi><mrow><mo>(</mo><msub><mi>&theta;</mi><mn>0</mn></msub><mo>,</mo><msub><mi>f</mi><mrow><mi>d</mi><mn>0</mn></mrow></msub><mo>)</mo></mrow><mo>,</mo><mi>s</mi><mrow><mo>(</mo><msub><mi>&theta;</mi><mn>3</mn></msub><mo>,</mo><msub><mi>f</mi><mrow><mi>d</mi><mn>3</mn></mrow></msub><mo>)</mo></mrow><mo>,</mo><mi>s</mi><mrow><mo>(</mo><msub><mi>&theta;</mi><mn>4</mn></msub><mo>,</mo><msub><mi>f</mi><mrow><mi>d</mi><mn>4</mn></mrow></msub><mo>)</mo></mrow><mo>,</mo><msub><mi>s</mi><mi>&theta;</mi></msub><mrow><mo>(</mo><msub><mi>&theta;</mi><mn>0</mn></msub><mo>,</mo><msub><mi>f</mi><mrow><mi>d</mi><mn>0</mn></mrow></msub><mo>)</mo></mrow><mo>&rsqb;</mo><mo>;</mo></mrow>]]></math><img file="FDA0001038759020000025.GIF" wi="1309" he="87" /></maths>6b)设定时域差波束响应矢量为:f<sub>Δt</sub>=[1,0,0,0,0]<sup>T</sup>;6c)根据步骤2)、步骤6a)和步骤6b)的结果,计算时域差波束权矢量:<maths num="0006"><math><![CDATA[<mrow><msub><mi>w</mi><mrow><mi>&Delta;</mi><mi>t</mi></mrow></msub><mo>=</mo><msup><mover><mi>R</mi><mo>^</mo></mover><mrow><mo>-</mo><mn>1</mn></mrow></msup><msub><mi>C</mi><mrow><mi>&Delta;</mi><mi>t</mi></mrow></msub><msup><mrow><mo>(</mo><msubsup><mi>C</mi><mrow><mi>&Delta;</mi><mi>t</mi></mrow><mi>H</mi></msubsup><msup><mover><mi>R</mi><mo>^</mo></mover><mrow><mo>-</mo><mn>1</mn></mrow></msup><msub><mi>C</mi><mrow><mi>&Delta;</mi><mi>t</mi></mrow></msub><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><msub><mi>f</mi><mrow><mi>&Delta;</mi><mi>t</mi></mrow></msub><mo>;</mo></mrow>]]></math><img file="FDA0001038759020000026.GIF" wi="630" he="99" /></maths>7)根据步骤4d)、步骤5d)和步骤6c)的结果,得到和波束的输出信号Y<sub>Σ</sub>、空域差波束的输出信号Y<sub>Δs</sub>和时域差波束的输出信号Y<sub>Δt</sub>:Y<sub>Σ</sub>=w<sub>Σ</sub><sup>H</sup>X,Y<sub>Δs</sub>=w<sub>Δs</sub><sup>H</sup>X,Y<sub>Δt</sub>=w<sub>Δt</sub><sup>H</sup>X。
地址 710071 陕西省西安市太白南路2号