发明名称 基于四阶累积量的单基地MIMO雷达分布式目标角度估计方法
摘要 本发明提供的是一种基于四阶累积量的单基地MIMO雷达分布式目标角度估计方法。M个收发共置的天线阵元,发射端发射相互正交的相位编码信号,接收端接收相互正交的相位编码信号;接收端每个接收阵元的匹配滤波器对接收到的正交信号进行匹配滤波进行分离,得到分布式目标接收信号矩阵;利用分布式目标接收信号矩阵,计算接收信号的四阶累积量矩阵;对四阶累积量矩阵进行特征值分解,求得相互正交的信号子空间与噪声子空间,利用多重信号分类算法构造空间谱函数;利用二维谱峰搜索计算分布式目标角度。本发明即使在色噪声环境下对中心波达方向角仍有很高的估计精度,且能够估计出分布式目标的扩展角度、实现中心波达方向角与扩展角的自动配对。
申请公布号 CN103926573B 申请公布日期 2017.02.22
申请号 CN201410153039.7 申请日期 2014.04.17
申请人 哈尔滨工程大学 发明人 王伟;王犇;刘琦;李欣;黄平;王冠男;司佳希;马跃华;王慧
分类号 G01S7/41(2006.01)I 主分类号 G01S7/41(2006.01)I
代理机构 代理人
主权项 一种基于四阶累积量的单基地MIMO雷达分布式目标角度估计方法,其特征是:(1)M个收发共置的天线阵元,发射端发射相互正交的相位编码信号,接收端接收相互正交的相位编码信号,每个收发阵元的间距均为半波长;(2)接收端每个接收阵元的匹配滤波器对接收到的正交信号进行匹配滤波进行分离,得到分布式目标接收信号矩阵;(3)利用得到的分布式目标接收信号矩阵,计算接收信号的四阶累积量矩阵;(4)对四阶累积量矩阵进行特征值分解,求得相互正交的信号子空间与噪声子空间,利用多重信号分类算法构造空间谱函数;(5)利用二维谱峰搜索计算分布式目标角度,即中心波达方向角角和扩展角;所述的分布式目标接收信号矩阵按如下方法得到:1)接收阵列接收到的第q个发射脉冲重复周期的回波信号表示为:<img file="FDA0001109908540000011.GIF" wi="1254" he="158" />式中P是目标个数,α<sub>pq</sub>是第p个目标在第q个发射脉冲重复周期的雷达截面反射系数,<img file="FDA0001109908540000012.GIF" wi="855" he="78" />为接收导向矢量,λ为载波波长,<img file="FDA0001109908540000013.GIF" wi="194" he="62" />是分布式目标确定性信号分布函数,s(t)是发射的正交编码信号,N(t)是高斯色噪声;2)对接收信号进行匹配滤波,<img file="FDA0001109908540000014.GIF" wi="1253" he="157" />式表示为:<img file="FDA0001109908540000015.GIF" wi="534" he="134" />式中<img file="FDA0001109908540000016.GIF" wi="955" he="86" />S(t)=[s<sub>1</sub>(t),s<sub>2</sub>(t),…,s<sub>M</sub>(t)]<sup>T</sup>,n(t)表示经匹配滤波后的高斯色噪声,在快拍个数是Q时,得到分布式目标信号接收矩阵:<img file="FDA0001109908540000017.GIF" wi="710" he="63" />式中<img file="FDA0001109908540000018.GIF" wi="1099" he="71" />所述的四阶累积量矩阵构造按如下进行:1)对于n×1维接收信号x,其四阶累积量由下式确定:<maths num="0001"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><msub><mi>C</mi><mrow><mn>4</mn><mi>x</mi></mrow></msub><mrow><mo>(</mo><msub><mi>k</mi><mn>1</mn></msub><mo>,</mo><msub><mi>k</mi><mn>2</mn></msub><mo>,</mo><msub><mi>k</mi><mn>3</mn></msub><mo>,</mo><msub><mi>k</mi><mn>4</mn></msub><mo>)</mo></mrow><mo>=</mo><mi>c</mi><mi>u</mi><mi>m</mi><mo>{</mo><msub><mi>x</mi><mrow><mi>k</mi><mn>1</mn></mrow></msub><mo>,</mo><msub><mi>x</mi><mrow><mi>k</mi><mn>2</mn></mrow></msub><mo>,</mo><msub><mi>x</mi><mrow><mi>k</mi><mn>3</mn></mrow></msub><mo>,</mo><msub><mi>x</mi><mrow><mi>k</mi><mn>4</mn></mrow></msub><mo>}</mo></mrow></mtd></mtr><mtr><mtd><mrow><mo>=</mo><mi>E</mi><mo>{</mo><msub><mi>x</mi><mrow><mi>k</mi><mn>1</mn></mrow></msub><msub><mi>x</mi><mrow><mi>k</mi><mn>2</mn></mrow></msub><msubsup><mi>x</mi><mrow><mi>k</mi><mn>3</mn></mrow><mo>*</mo></msubsup><msubsup><mi>x</mi><mrow><mi>k</mi><mn>4</mn></mrow><mo>*</mo></msubsup><mo>}</mo><mo>-</mo><mi>E</mi><mo>{</mo><msub><mi>x</mi><mrow><mi>k</mi><mn>1</mn></mrow></msub><msubsup><mi>x</mi><mrow><mi>k</mi><mn>3</mn></mrow><mo>*</mo></msubsup><mo>}</mo><mi>E</mi><mo>{</mo><msub><mi>x</mi><mrow><mi>k</mi><mn>2</mn></mrow></msub><msubsup><mi>x</mi><mrow><mi>k</mi><mn>4</mn></mrow><mo>*</mo></msubsup><mo>}</mo></mrow></mtd></mtr><mtr><mtd><mrow><mo>-</mo><mi>E</mi><mo>{</mo><msub><mi>x</mi><mrow><mi>k</mi><mn>1</mn></mrow></msub><msubsup><mi>x</mi><mrow><mi>k</mi><mn>4</mn></mrow><mo>*</mo></msubsup><mo>}</mo><mi>E</mi><mo>{</mo><msub><mi>x</mi><mrow><mi>k</mi><mn>2</mn></mrow></msub><msubsup><mi>x</mi><mrow><mi>k</mi><mn>4</mn></mrow><mo>*</mo></msubsup><mo>}</mo><mo>-</mo><mi>E</mi><mo>{</mo><msub><mi>x</mi><mrow><mi>k</mi><mn>1</mn></mrow></msub><msub><mi>x</mi><mrow><mi>k</mi><mn>2</mn></mrow></msub><mo>}</mo><mi>E</mi><mo>{</mo><msubsup><mi>x</mi><mrow><mi>k</mi><mn>3</mn></mrow><mo>*</mo></msubsup><msubsup><mi>x</mi><mrow><mi>k</mi><mn>4</mn></mrow><mo>*</mo></msubsup><mo>}</mo></mrow></mtd></mtr></mtable><mo>;</mo></mrow>]]></math><img file="FDA0001109908540000019.GIF" wi="1214" he="223" /></maths>2)根据上式,分布式目标接收信号矩阵的四阶累积量矩阵的计算方法为:<maths num="0002"><math><![CDATA[<mrow><msub><mi>R</mi><mn>4</mn></msub><mo>=</mo><mi>E</mi><mo>{</mo><mrow><mo>(</mo><mi>Y</mi><mo>&CircleTimes;</mo><msup><mi>Y</mi><mo>*</mo></msup><mo>)</mo></mrow><msup><mrow><mo>(</mo><mi>Y</mi><mo>&CircleTimes;</mo><msup><mi>Y</mi><mo>*</mo></msup><mo>)</mo></mrow><mi>H</mi></msup><mo>}</mo><mo>-</mo><mi>E</mi><mo>{</mo><mi>Y</mi><mo>&CircleTimes;</mo><msup><mi>Y</mi><mo>*</mo></msup><mo>}</mo><mi>E</mi><mo>{</mo><msup><mrow><mo>(</mo><mi>Y</mi><mo>&CircleTimes;</mo><msup><mi>Y</mi><mo>*</mo></msup><mo>)</mo></mrow><mi>H</mi></msup><mo>}</mo><mo>-</mo><mi>E</mi><mo>{</mo><mi>Y</mi><mo>&CenterDot;</mo><msup><mi>Y</mi><mi>H</mi></msup><mo>}</mo><mi>E</mi><mo>{</mo><msup><mrow><mo>(</mo><mi>Y</mi><mo>&CenterDot;</mo><msup><mi>Y</mi><mi>H</mi></msup><mo>)</mo></mrow><mo>*</mo></msup><mo>}</mo><mo>;</mo></mrow>]]></math><img file="FDA0001109908540000021.GIF" wi="1614" he="63" /></maths>所述对四阶累积量矩阵进行特征值分解为:<maths num="0003"><math><![CDATA[<mrow><msub><mi>R</mi><mn>4</mn></msub><mo>=</mo><msub><mi>U</mi><mi>s</mi></msub><msub><mi>&Sigma;</mi><mi>s</mi></msub><msubsup><mi>U</mi><mi>s</mi><mi>H</mi></msubsup><mo>+</mo><msub><mi>U</mi><mi>n</mi></msub><msub><mi>&Sigma;</mi><mi>n</mi></msub><msubsup><mi>U</mi><mi>n</mi><mi>H</mi></msubsup></mrow>]]></math><img file="FDA0001109908540000022.GIF" wi="491" he="64" /></maths>式中U<sub>s</sub>是R<sub>4</sub>最大的P<sup>2</sup>个特征值对应的特征向量所张成的信号子空间,Σ<sub>s</sub>为最大的P<sup>2</sup>个特征值组成的对角阵;U<sub>n</sub>是M<sup>4</sup>‑P<sup>2</sup>个最小的特征值对应的特征向量所张成的噪声子空间,Σ<sub>n</sub>为M<sup>4</sup>‑P<sup>2</sup>个最小的特征值组成的对角阵;所述的基于多重信号分类算法的空间谱函数构造方法为:1)信号子空间由导向矢量矩阵的列向量张成,即<img file="FDA0001109908540000023.GIF" wi="526" he="79" />2)信号子空间与噪声子空间相互正交<img file="FDA0001109908540000024.GIF" wi="198" he="70" />3)得到空间谱函数<img file="FDA0001109908540000025.GIF" wi="470" he="131" />通过对空间谱函数进行二维谱峰搜索计算得到分布式目标中心波达方向角角和扩展角。
地址 150001 黑龙江省哈尔滨市南岗区南通大街145号哈尔滨工程大学科技处知识产权办公室