发明名称 一种适用于大规模MIMO系统的二维波达角的测定方法
摘要 本发明公开了一种适用于大规模MIMO系统的二维波达角的测定方法,属于无限移动通信技术领域。本发明基于M<sub>x</sub>×M<sub>y</sub>的二维均匀矩形天线阵列执行下列步骤:基于接收信号建立接收信号的协方差矩阵,通过协方差矩阵获取信号子空间,将天线阵列划分为4个子阵,并得到相应的选择信号子空间,获得旋转矩阵,最后分别确定信号源的俯仰角与方位角。该发明在大规模MIMO场景下具有较小的测定误差,避免了角度搜多与非线性优化。因此,本发明具有测定精度高,数据处理复杂度低等优点,应用于大规模MIMO场景可有效提高系统的性能和实用价值。
申请公布号 CN105510872A 申请公布日期 2016.04.20
申请号 CN201610043980.2 申请日期 2016.01.22
申请人 电子科技大学 发明人 孟会鹏;郑植;杨雨轩;刘柯宏;葛琰;杨娇
分类号 G01S3/14(2006.01)I 主分类号 G01S3/14(2006.01)I
代理机构 电子科技大学专利中心 51203 代理人 周刘英
主权项 一种适用于大规模MIMO系统的二维波达角的测定方法,其特征在于,基于M<sub>x</sub>×M<sub>y</sub>的二维均匀矩形天线阵列执行下列步骤:步骤1:基于接收信号计算接收信号的协方差矩阵;步骤2:对协方差矩阵进行特征值分解,基于前K个最大特征值所对应的特征向量e<sub>1</sub>,e<sub>2</sub>,…e<sub>K</sub>构建信号子空间E<sub>S</sub>=[e<sub>1</sub>,e<sub>2</sub>,...,e<sub>K</sub>],其中K为信号源个数;步骤3:把M<sub>x</sub>×M<sub>y</sub>的二维均匀矩形天线阵划分为4个子阵,各子阵的选择矩阵J<sub>1</sub>、J<sub>2</sub>、J<sub>3</sub>、J<sub>4</sub>分别为:<img file="FDA0000912622900000011.GIF" wi="966" he="245" /><img file="FDA0000912622900000012.GIF" wi="1029" he="242" /><math><![CDATA[<mrow><msub><mrow><mo>&lsqb;</mo><msub><mi>J</mi><mn>3</mn></msub><mo>&rsqb;</mo></mrow><mrow><mover><mi>m</mi><mo>~</mo></mover><mo>,</mo><mover><mi>n</mi><mo>~</mo></mover></mrow></msub><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mn>1</mn><mo>,</mo></mrow></mtd><mtd><mrow><mover><mi>n</mi><mo>~</mo></mover><mo>=</mo><mover><mi>m</mi><mo>~</mo></mover><mo>,</mo><mover><mi>m</mi><mo>~</mo></mover><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>...</mo><mo>,</mo><msub><mover><mi>M</mi><mo>~</mo></mover><mn>2</mn></msub><mo>,</mo></mrow></mtd></mtr><mtr><mtd><mrow><mn>0</mn><mo>,</mo></mrow></mtd><mtd><mrow><mi>o</mi><mi>t</mi><mi>h</mi><mi>e</mi><mi>r</mi><mi>w</mi><mi>i</mi><mi>s</mi><mi>e</mi><mo>,</mo></mrow></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000912622900000013.GIF" wi="747" he="166" /></maths><math><![CDATA[<mrow><msub><mrow><mo>&lsqb;</mo><msub><mi>J</mi><mn>4</mn></msub><mo>&rsqb;</mo></mrow><mrow><mover><mi>m</mi><mo>~</mo></mover><mo>,</mo><mover><mi>n</mi><mo>~</mo></mover></mrow></msub><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mn>1</mn><mo>,</mo></mrow></mtd><mtd><mrow><mover><mi>n</mi><mo>~</mo></mover><mo>=</mo><mover><mi>m</mi><mo>~</mo></mover><mo>+</mo><msub><mi>M</mi><mi>x</mi></msub><mo>,</mo><mover><mi>m</mi><mo>~</mo></mover><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>...</mo><mo>,</mo><msub><mover><mi>M</mi><mo>~</mo></mover><mn>2</mn></msub><mo>,</mo></mrow></mtd></mtr><mtr><mtd><mrow><mn>0</mn><mo>,</mo></mrow></mtd><mtd><mrow><mi>o</mi><mi>t</mi><mi>h</mi><mi>e</mi><mi>r</mi><mi>w</mi><mi>i</mi><mi>s</mi><mi>e</mi><mo>,</mo></mrow></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000912622900000014.GIF" wi="860" he="166" /></maths>其中J<sub>1</sub>和J<sub>2</sub>为<img file="FDA0000912622900000015.GIF" wi="165" he="70" />矩阵,J<sub>3</sub>和J<sub>4</sub>为<img file="FDA0000912622900000016.GIF" wi="173" he="71" />矩阵,<img file="FDA0000912622900000017.GIF" wi="397" he="95" /><math><![CDATA[<mrow><msub><mover><mi>M</mi><mo>~</mo></mover><mn>2</mn></msub><mo>=</mo><msub><mi>M</mi><mi>x</mi></msub><mrow><mo>(</mo><msub><mi>M</mi><mi>y</mi></msub><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000912622900000018.GIF" wi="402" he="87" /></maths>M=M<sub>x</sub>×M<sub>y</sub>;基于各个子阵的选择矩阵J<sub>q</sub>将信号子空间E<sub>S</sub>划分为4个选择信号子空间<img file="FDA0000912622900000019.GIF" wi="254" he="71" />其中q=1,2,3,4;步骤4:通过总体最小二乘法构建选择信号子空间E<sub>S1</sub>和E<sub>S2</sub>的旋转矩阵ψ<sub>1</sub>,选择信号子空间E<sub>S3</sub>和E<sub>S4</sub>的旋转矩阵ψ<sub>2</sub>;步骤5:确定方位角和俯仰角:对旋转矩阵ψ<sub>1</sub>进行特征值分解得到对角矩阵Λ<sub>1</sub>,其对角元素为ψ<sub>1</sub>的特征值,对ψ<sub>2</sub>进行特征值分解得到对角矩阵Λ<sub>2</sub>,其对角元素为ψ<sub>2</sub>的特征值,对Λ<sub>1</sub>和Λ<sub>2</sub>的对角元素进行配对对齐,得到对齐的对角矩阵Λ<sub>1</sub>和<img file="FDA00009126229000000110.GIF" wi="95" he="85" />定义ξ<sub>1,k</sub>和ξ<sub>2,k</sub>分别为Λ<sub>1</sub>和<img file="FDA00009126229000000111.GIF" wi="75" he="78" />的第k个对角线元素,则方位角<img file="FDA0000912622900000021.GIF" wi="62" he="78" />和俯仰角<img file="FDA00009126229000000215.GIF" wi="58" he="78" />分别为:<math><![CDATA[<mrow><msub><mover><mi>&theta;</mi><mo>^</mo></mover><mi>k</mi></msub><mo>=</mo><mi>arctan</mi><mrow><mo>(</mo><mfrac><mrow><msub><mi>ln&xi;</mi><mrow><mn>2</mn><mo>,</mo><mi>k</mi></mrow></msub></mrow><mrow><msub><mi>ln&xi;</mi><mrow><mn>1</mn><mo>,</mo><mi>k</mi></mrow></msub></mrow></mfrac><mo>)</mo></mrow><mo>,</mo><msub><mover><mi>&phi;</mi><mo>^</mo></mover><mi>k</mi></msub><mo>=</mo><mi>arcsin</mi><mrow><mo>(</mo><mrow><mfrac><mn>1</mn><mi>u</mi></mfrac><msqrt><mrow><mo>-</mo><munderover><mo>&Sigma;</mo><mrow><mi>l</mi><mo>=</mo><mn>1</mn></mrow><mn>2</mn></munderover><msup><mrow><mo>(</mo><mrow><msub><mi>ln&xi;</mi><mrow><mi>l</mi><mo>,</mo><mi>k</mi></mrow></msub></mrow><mo>)</mo></mrow><mn>2</mn></msup></mrow></msqrt></mrow><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA00009126229000000216.GIF" wi="1134" he="173" /></maths>其中k=1,2,...,K,u=2πd/λ,d为相邻阵元之间的间距,λ为载波波长。
地址 611731 四川省成都市高新区(西区)西源大道2006号