发明名称 基于特征值重构的相干信号到达角估计方法
摘要 基于特征值重构的相干信号到达角估计方法利用信号协方差矩阵的大特征值对应的特征向量进行重新构造新的协方差矩阵,通过对新的协方差矩阵进行特征值分解,估计信号的波达方向;当信号相干时,解除由于多径等因素影响导致信号的相关性,使得新信号协方差矩阵的秩得到恢复,同时在同样的天线阵元数下能够估计更多的相干信号的到达方向,本发明相较于传统的方法有两大优点:1、能够估计空间相干信号的到达方向,2、估计同样数目的信号到达角,本发明只需更少的天线阵元,因此大大降低了运行成本。
申请公布号 CN101592721B 申请公布日期 2011.09.14
申请号 CN200910032757.8 申请日期 2009.06.19
申请人 中国人民解放军理工大学 发明人 郭艳;李宁;钱祖平;吴启辉;王金龙;刘学亮
分类号 G01S3/14(2006.01)I;G01S7/00(2006.01)I;G01S7/36(2006.01)I 主分类号 G01S3/14(2006.01)I
代理机构 南京经纬专利商标代理有限公司 32200 代理人 叶连生
主权项 1.一种基于特征值重构的相干信号到达角估计方法,其特征在于该方法利用信号协方差矩阵的大特征值对应的特征向量进行重新构造新的协方差矩阵,通过对新的协方差矩阵进行特征值分解,估计信号的波达方向;当信号相干时,解除由于多径因素影响导致信号的相关性,使得新信号协方差矩阵的秩得到恢复,同时在同样的天线阵元数下能够估计更多的相干信号的到达方向,该方法的具体步骤为:1.)由天线阵列接收数据向量X(n),根据估计式<img file="FSB00000534525900011.GIF" wi="487" he="137" />计算数据协方差矩阵R<sub>x</sub>的估计值<img file="FSB00000534525900012.GIF" wi="88" he="58" />其中N为数据快拍数,“H”为矩阵共轭转置运算,n表示时刻;2.)对协方差矩阵<img file="FSB00000534525900013.GIF" wi="64" he="60" />进行特征值分解,对<img file="FSB00000534525900014.GIF" wi="63" he="56" />进行特征值分解,得到<img file="FSB00000534525900015.GIF" wi="335" he="81" />U<sub>1</sub>=[u<sub>1</sub>,…,u<sub>M</sub>]表示由特征向量u<sub>1</sub>,…,u<sub>M</sub>构成的矩阵,∑<sub>1</sub>=diag(λ<sub>1</sub>,…,λ<sub>M</sub>)表示以特征值λ<sub>1</sub>,…,λ<sub>M</sub>为对角线元素的对角矩阵;3.)将特征值以降序排列,q个大特征值对应的特征向量构成信号子空间<img file="FSB00000534525900016.GIF" wi="87" he="73" />即<img file="FSB00000534525900017.GIF" wi="378" he="80" />u<sub>1</sub>,…,u<sub>q</sub>表示q个大特征值对应的特征向量,其中,u<sub>1</sub>=e<sub>1</sub>=[e<sub>1</sub>(1) e<sub>1</sub>(2) … e<sub>1</sub>(M)]<sup>T</sup>,…,u<sub>q</sub>=e<sub>q</sub>=[e<sub>q</sub>(1) e<sub>q</sub>(2) … e<sub>q</sub>(M)]<sup>T</sup>;4.)取<img file="FSB00000534525900019.GIF" wi="1469" he="395" /><img file="FSB000005345259000110.GIF" wi="1309" he="394" />5.)重新构造协方差矩阵得到<img file="FSB000005345259000111.GIF" wi="412" he="63" />6.)对<img file="FSB000005345259000112.GIF" wi="43" he="52" />进行特征分解,得到<img file="FSB000005345259000113.GIF" wi="262" he="54" />其中U=[u<sub>1</sub>,…,u<sub>M</sub>]表示由特征向量u<sub>1</sub>,…,u<sub>M</sub>构成的矩阵,∑=diag(λ<sub>1</sub>,…,λ<sub>M</sub>)表示以特征值λ<sub>1</sub>,…,λ<sub>M</sub>为对角线元素的对角矩阵,7.)将特征值以降序排列,q个大特征值对应的特征向量构成信号子空间,即 <img file="FSB00000534525900021.GIF" wi="394" he="113" />u<sub>1</sub>,…,u<sub>q</sub>表示q个大特征值对应的特征向量,而M-q个小特征值对应的特征向量构成噪声子空间,即<img file="FSB00000534525900022.GIF" wi="417" he="119" />u<sub>q+1</sub>,…,u<sub>M</sub>表示M-q个小特征值对应的特征向量,8.)根据噪声子空间<img file="FSB00000534525900023.GIF" wi="95" he="73" />构造谱估计公式:<img file="FSB00000534525900024.GIF" wi="579" he="185" />其中P(ω)为功率谱值,a(ω)为导引向量,9.)根据谱估计公式计算结果,极大值点对应的角度就是信号入射方向θ<sub>i</sub>,i=1,…,q的估计值。 
地址 210007 江苏省南京市御道街标营2号