发明名称 一种多径环境下的到来角估计算法
摘要 本发明涉及一种多径环境下的到来角估计算法,属于无线通信技术领域。本发明针对移动通信中存在的多径问题,基于多径信号的衰落特性,提出了将时间平均与合成空间谱方法结合的算法,该方法解决了多径环境下信号源相干、天线数相对少时的到来角估计问题,不存在阵列孔径的损失;同时也有效的解决了多径环境下的低信噪比、来波方位间隔较近时的到来角估计问题。
申请公布号 CN103558581B 申请公布日期 2016.05.04
申请号 CN201310340255.8 申请日期 2013.08.06
申请人 大连工业大学 发明人 王智森;申莹;唐艺鸣;于志远
分类号 G01S3/46(2006.01)I 主分类号 G01S3/46(2006.01)I
代理机构 大连东方专利代理有限责任公司 21212 代理人 曲永祚;李洪福
主权项 一种多径环境下的到来角估计算法,其特征在于:接收一个移动的移动用户的m个时间间隔的信号,在这些时间间隔,到来角的所有波前认为是不变的,在每个时间间隔,接收信号的协方差矩阵为:<maths num="0001"><math><![CDATA[<mfenced open = "" close = ""><mtable><mtr><mtd><mrow><msub><mi>R</mi><msup><mi>X</mi><mi>j</mi></msup></msub><mo>=</mo><mi>E</mi><mrow><mo>&lsqb;</mo><mrow><msup><mi>X</mi><mi>j</mi></msup><msup><mrow><mo>(</mo><msup><mi>X</mi><mi>j</mi></msup><mo>)</mo></mrow><mi>H</mi></msup></mrow><mo>&rsqb;</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mo>=</mo><msup><mi>AC</mi><mi>j</mi></msup><msub><mi>R</mi><mi>s</mi></msub><msup><mrow><mo>(</mo><msup><mi>C</mi><mi>j</mi></msup><mo>)</mo></mrow><mi>H</mi></msup><msup><mi>A</mi><mi>H</mi></msup><mo>+</mo><msubsup><mi>&sigma;</mi><mi>n</mi><mn>2</mn></msubsup><mi>I</mi><mo>,</mo></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000851770580000011.GIF" wi="700" he="263" /></maths>经推导,<img file="FDA0000851770580000017.GIF" wi="379" he="87" />而信号子空间的秩为P,因此协方差矩阵<img file="FDA0000851770580000018.GIF" wi="95" he="83" />不能包含信号子空间,此时取m个时间间隔的协方差矩阵,做平均:<maths num="0002"><math><![CDATA[<mfenced open = "" close = ""><mtable><mtr><mtd><mrow><msub><mi>R</mi><mrow><mi>T</mi><mo>-</mo><mi>X</mi></mrow></msub><mo>=</mo><mfrac><mn>1</mn><mi>m</mi></mfrac><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msub><mi>R</mi><msup><mi>X</mi><mi>j</mi></msup></msub></mrow></mtd></mtr><mtr><mtd><mrow><mo>=</mo><mi>A</mi><mrow><mo>(</mo><mrow><mfrac><mn>1</mn><mi>m</mi></mfrac><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><mrow><msup><mi>C</mi><mi>j</mi></msup><msub><mi>R</mi><mi>s</mi></msub><msup><mrow><mo>(</mo><msup><mi>C</mi><mi>j</mi></msup><mo>)</mo></mrow><mi>H</mi></msup></mrow></mrow><mo>)</mo></mrow><msup><mi>A</mi><mi>H</mi></msup><mo>+</mo><msubsup><mi>&sigma;</mi><mi>n</mi><mn>2</mn></msubsup><mi>I</mi><mo>,</mo></mrow></mtd></mtr><mtr><mtd><mrow><mo>=</mo><mi>A</mi><mover><msub><mi>R</mi><mi>s</mi></msub><mo>&OverBar;</mo></mover><msup><mi>A</mi><mi>H</mi></msup><mo>+</mo><msubsup><mi>&sigma;</mi><mi>n</mi><mn>2</mn></msubsup><mi>I</mi><mo>,</mo></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000851770580000012.GIF" wi="854" he="447" /></maths>其中,<img file="FDA0000851770580000013.GIF" wi="661" he="157" />由于各个时间间隔之间是快变的,所以<img file="FDA0000851770580000014.GIF" wi="333" he="77" />又因为<img file="FDA0000851770580000015.GIF" wi="670" he="86" />rank(A)),rank(A)=P,所以,平均后的协方差矩阵R<sub>T‑X</sub>也满秩,即rank(R<sub>T‑X</sub>)=P,其中,A是阵列导向矢量矩阵,由信号的方向向量a(θ)组成,θ为信号的到来角;R<sub>s</sub>是Q个独立信号源的Q×Q的协方差矩阵;H是共轭转置符号;<img file="FDA0000851770580000019.GIF" wi="66" he="71" />是噪声功率;I是M×M的单位矩阵,C<sup>j</sup>是P×Q的矩阵;所述阵列输出功率与到来角的关系为:<maths num="0003"><math><![CDATA[<mrow><msub><mi>P</mi><mrow><mi>T</mi><mo>-</mo><mi>S</mi><mi>S</mi><mi>M</mi></mrow></msub><mrow><mo>(</mo><mi>&theta;</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>P</mi></munderover><msubsup><mi>&lambda;</mi><mi>j</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mo>|</mo><msup><mi>a</mi><mi>H</mi></msup><mrow><mo>(</mo><mi>&theta;</mi><mo>)</mo></mrow><msub><mi>q</mi><mi>j</mi></msub><msup><mo>|</mo><mn>2</mn></msup></mrow><mrow><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mi>P</mi><mo>+</mo><mn>1</mn></mrow><mi>M</mi></munderover><mo>|</mo><msup><mi>a</mi><mi>H</mi></msup><mrow><mo>(</mo><mi>&theta;</mi><mo>)</mo></mrow><msub><mi>q</mi><mi>j</mi></msub><msup><mo>|</mo><mn>2</mn></msup></mrow></mfrac><mo>,</mo></mrow>]]></math><img file="FDA0000851770580000016.GIF" wi="813" he="382" /></maths>其中,λ<sub>j</sub>,q<sub>j</sub>分别为平均后的R<sub>T‑X</sub>的第j个特征值和对应的特征向量,a(θ)是信号的方向向量,H是共轭转置符号。
地址 116034 辽宁省大连市甘井子区轻工苑1号