发明名称 一种基于共轭扩展的波达方向估计方法
摘要 本发明公开了一种基于共轭扩展的波达方向估计方法,该方法先计算接收信号之间的二阶统计量,利用延迟自相关函数的性质得到导向矩阵的共扼矩阵,然后利用扩展矢量的延迟采样构造出一个新的伪快拍矩阵,最后对这个矩阵用四阶累积量的方法来估计信号的波达方向。通过阵列扩展的信号波达方向估计方法,本发明能够大大增加估计的信号波达方向数目,并且能够提高信号波达方向的估计精度。
申请公布号 CN102694588B 申请公布日期 2014.12.03
申请号 CN201210199310.1 申请日期 2012.06.15
申请人 华南师范大学 发明人 单志龙;刘方伟;史景伦
分类号 H04L5/00(2006.01)I;H04B7/08(2006.01)I 主分类号 H04L5/00(2006.01)I
代理机构 广州新诺专利商标事务所有限公司 44100 代理人 张玲春
主权项 一种基于共轭扩展的波达方向估计方法,其特征在于,该方法首先对接收的阵列信号进行二阶统计量处理,进行延迟采样构造出一伪快拍矩阵,最后对该伪快拍矩阵用四阶累积量的方法来估计信号的波达方向; 该方法具体步骤为: 1)假设天线阵列为均匀线阵,每个阵元都是一致全向,阵元数为M,阵元间距为d,假如D个远场窄带非相干信号分别以θ<sub>i</sub>,i=1,2,...D投射到M个天线阵列上,根据第m个天线阵列的输出信号X<sub>m</sub>(t),其中 <img file="FDA0000563318760000011.GIF" wi="1140" he="146" />其中X<sub>m</sub>(t)为第m个天线的输出信号;S<sub>d</sub>(t)为第d个信号;n<sub>m</sub>(t)为高斯白噪声;λ表示波长;θ为第d个信号的波达方向; 利用延迟自相关的性质,按照公式 <img file="FDA0000563318760000012.GIF" wi="1130" he="140" />构造矩阵R<sub>x</sub>(τ),R'<sub>x</sub>(τ);其中<img file="dest_path_BDA00001769650900052.GIF" wi="159" he="55" />表示第m个阵元输出X<sub>m</sub>(t)与第一个阵元输出X<sub>1</sub>(t)间的延迟互相关函数,<img file="FDA0000563318760000014.GIF" wi="139" he="86" />表示s<sub>d</sub>(t)的延迟自相关函数,<img file="FDA0000563318760000015.GIF" wi="126" he="70" />表示噪声的延迟互相关函数;2)根据<img file="FDA0000563318760000016.GIF" wi="378" he="160" />对R(τ)在伪快拍数N<sub>p</sub>个延迟τ=T<sub>s</sub>,2T<sub>s</sub>,...N<sub>p</sub>T<sub>s</sub>上采样,得到一个伪快拍矩阵<img file="FDA0000563318760000017.GIF" wi="76" he="81" /><img file="FDA0000563318760000018.GIF" wi="651" he="97" />3)对所述伪快拍矩阵<img file="FDA0000563318760000019.GIF" wi="50" he="70" />做四阶累积量计算得到一个四阶累积量矩阵C,其中:<img file="FDA00005633187600000110.GIF" wi="1547" he="109" />4)将所述四阶累积量矩阵C演化成以下形式: <img file="FDA00005633187600000111.GIF" wi="310" he="85" />其中, ∑<sub>s</sub>为对角线上的元素为R<sub>s</sub>(τ)的峰值的对角矩阵; 得到导向矩阵<img file="FDA0000563318760000021.GIF" wi="75" he="74" /><img file="FDA0000563318760000022.GIF" wi="419" he="92" />其中<img file="FDA0000563318760000023.GIF" wi="924" he="176" />所述步骤4)中的导向矩阵<img file="FDA0000563318760000024.GIF" wi="52" he="65" />通过以下方式得到:设该方法的初始条件为:信号数D,信号入射方向为θ<sub>i</sub>,阵元数为M,其中一阵元的坐标为(a<sub>m</sub>,b<sub>m</sub>),得到: <img file="FDA0000563318760000025.GIF" wi="969" he="379" />以及,<img file="FDA0000563318760000026.GIF" wi="946" he="463" />根据公式: <img file="FDA0000563318760000027.GIF" wi="909" he="166" />得到导向矩阵<img file="FDA0000563318760000028.GIF" wi="75" he="66" /><img file="FDA0000563318760000029.GIF" wi="417" he="102" />5)对所述四阶累积量矩阵C进行特征值分解,得到特征值<img file="FDA00005633187600000210.GIF" wi="75" he="71" /><img file="FDA00005633187600000211.GIF" wi="608" he="88" />其中,U<sub>S</sub>是由大特征值对应的特征矢量张成的子空间,即信号子空间; U<sub>N</sub>是由小特征值对应的特征矢量张成的子空间,即噪声子空间; 6)对谱估计式P<sub>MUSIC</sub>(θ)进行谱峰搜索,其中,所述谱估计式P<sub>MUSIC</sub>(θ)为: <img file="FDA0000563318760000031.GIF" wi="620" he="142" />该谱估计式的P<sub>MUSIC</sub>(θ)谱峰对应的横坐标值即为信号的波达方向。 
地址 510631 广东省广州市中山大道西55号华南师范大学计算机学院