发明名称 双向全双工MIMO中继通信系统中的自干扰消除方法
摘要 本发明公开一种全双工MIMO中继自干扰消除方法,主要解决现有技术中只考虑最小化自干扰而忽视系统容量性能的问题。其实现步骤是:1)估计信道参数;2)设计第一源节点的发射,接收波束成形矩阵以消除第一源节点的自干扰;3)设计第二源节点的发射,接收波束成形矩阵,以消除第二源节点的自干扰;4)设计中继节点的发射,接收波束成形矩阵,以消除中继节点的自干扰,并同时选择符合条件的发射,接收波束成形矩阵最大化系统容量。本发明不仅消除了中继节点的自干扰,又消除了两个源节点的自干扰,且能在消除自干扰的同时最大化系统的容量,可用于未来第五代无线移动通信的分布式协作传输系统。
申请公布号 CN106301457A 申请公布日期 2017.01.04
申请号 CN201610876815.5 申请日期 2016.10.08
申请人 西安电子科技大学 发明人 刘毅;朱芳美;吴炯;陈甜甜;张海林
分类号 H04B1/52(2015.01)I;H04L25/03(2006.01)I 主分类号 H04B1/52(2015.01)I
代理机构 陕西电子工业专利中心 61205 代理人 王品华
主权项 一种最大化系统容量的双向全双工MIMO中继自干扰消除方法,包括:(1)采用最小均方误差信道估计方法估计如下信道参数:第一源节点到中继节点的信道参数H<sub>1r</sub>,第二源节点到中继节点的信道参数H<sub>2r</sub>,中继节点到第一源节点的信道参数H<sub>r1</sub>,中继节点到第二源节点的信道参数H<sub>r2</sub>,中继节点的自干扰信道参数H<sub>rr</sub>,第一源节点自干扰信道参数H<sub>11</sub>,第二源节点的自干扰信道参数H<sub>22</sub>;(2)采用零空间投影法设计第一源节点的接收波束成形矩阵U<sub>1</sub>和发射波束成形矩阵V<sub>1</sub>,消除第一源节点的自干扰;(3)采用零空间投影法设计第二源节点的接收波束成形矩阵U<sub>2</sub>和发射波束成形矩阵V<sub>2</sub>,消除第二源节点的自干扰;(4)设计中继节点的发射波束成形矩阵W<sub>t</sub>和接收波束成形矩阵W<sub>r</sub>,消除中继节点的自干扰:(4a)对中继节点的自干扰信道H<sub>rr</sub>进行如下奇异值分解:<maths num="0001"><math><![CDATA[<mrow><msub><mi>H</mi><mrow><mi>r</mi><mi>r</mi></mrow></msub><mo>=</mo><msub><mi>U</mi><mrow><mi>r</mi><mi>r</mi></mrow></msub><msub><mi>&Sigma;</mi><mrow><mi>r</mi><mi>r</mi></mrow></msub><msubsup><mi>V</mi><mrow><mi>r</mi><mi>r</mi></mrow><mi>H</mi></msubsup></mrow>]]></math><img file="FDA0001126347360000011.GIF" wi="334" he="72" /></maths>其中,U<sub>rr</sub>是H<sub>rr</sub>的左奇异矩阵,∑<sub>rr</sub>是由H<sub>rr</sub>的奇异值组成的对角矩阵,V<sub>rr</sub>是H<sub>rr</sub>右奇异矩阵;(4b)利用(4a)中的U<sub>rr</sub>,V<sub>rr</sub>得到中继节点的发射波束成形矩阵W<sub>t</sub>和接收波束成形矩阵W<sub>r</sub>:<maths num="0002"><math><![CDATA[<mrow><msub><mi>W</mi><mi>r</mi></msub><mo>=</mo><msubsup><mi>D</mi><mrow><mi>R</mi><mi>r</mi></mrow><mi>T</mi></msubsup><msubsup><mi>U</mi><mrow><mi>r</mi><mi>r</mi></mrow><mi>H</mi></msubsup><mo>,</mo><msub><mi>W</mi><mi>t</mi></msub><mo>=</mo><msub><mi>V</mi><mrow><mi>r</mi><mi>r</mi></mrow></msub><msub><mi>D</mi><mrow><mi>R</mi><mi>t</mi></mrow></msub><mo>,</mo></mrow>]]></math><img file="FDA0001126347360000012.GIF" wi="726" he="78" /></maths>(4c)根据(4a)和(4b)得到中继节点的自干扰项表达式:<maths num="0003"><math><![CDATA[<mrow><msub><mi>W</mi><mi>r</mi></msub><msub><mi>H</mi><mrow><mi>r</mi><mi>r</mi></mrow></msub><msub><mi>W</mi><mi>t</mi></msub><mo>=</mo><msubsup><mi>D</mi><mrow><mi>R</mi><mi>r</mi></mrow><mi>T</mi></msubsup><msubsup><mi>U</mi><mrow><mi>r</mi><mi>r</mi></mrow><mi>H</mi></msubsup><msub><mi>U</mi><mrow><mi>r</mi><mi>r</mi></mrow></msub><msub><mi>&Sigma;</mi><mrow><mi>r</mi><mi>r</mi></mrow></msub><msubsup><mi>V</mi><mrow><mi>r</mi><mi>r</mi></mrow><mi>H</mi></msubsup><msub><mi>V</mi><mrow><mi>r</mi><mi>r</mi></mrow></msub><msub><mi>D</mi><mrow><mi>R</mi><mi>t</mi></mrow></msub><mo>=</mo><msubsup><mi>D</mi><mrow><mi>R</mi><mi>r</mi></mrow><mi>T</mi></msubsup><msub><mi>&Sigma;</mi><mrow><mi>r</mi><mi>r</mi></mrow></msub><msub><mi>D</mi><mrow><mi>R</mi><mi>t</mi></mrow></msub><mo>=</mo><msubsup><mi>&Sigma;</mi><mrow><mi>n</mi><mo>=</mo><msub><mi>N</mi><mrow><mi>R</mi><mi>r</mi></mrow></msub><mo>+</mo><msub><mi>N</mi><mrow><mi>R</mi><mi>t</mi></mrow></msub><mo>-</mo><mrow><mo>(</mo><mrow><msub><mover><mi>N</mi><mo>^</mo></mover><mrow><mi>R</mi><mi>r</mi></mrow></msub><mo>&times;</mo><msub><mover><mi>N</mi><mo>^</mo></mover><mrow><mi>R</mi><mi>t</mi></mrow></msub></mrow><mo>)</mo></mrow><mo>+</mo><mn>1</mn></mrow><mrow><mi>min</mi><mrow><mo>{</mo><mrow><msub><mi>N</mi><mrow><mi>R</mi><mi>r</mi></mrow></msub><mo>,</mo><msub><mi>N</mi><mrow><mi>R</mi><mi>t</mi></mrow></msub></mrow><mo>}</mo></mrow></mrow></msubsup><msubsup><mi>&sigma;</mi><mrow><mi>R</mi><mi>R</mi></mrow><mn>2</mn></msubsup><mrow><mo>&lsqb;</mo><mi>n</mi><mo>&rsqb;</mo></mrow></mrow>]]></math><img file="FDA0001126347360000013.GIF" wi="1575" he="108" /></maths>其中<img file="FDA0001126347360000014.GIF" wi="334" he="86" />是行子集选择矩阵,<img file="FDA0001126347360000015.GIF" wi="318" he="87" />是列子集选择矩阵,N<sub>Rt</sub>是中继节点的发射天线数,N<sub>Rr</sub>是中继节点的接收天线数,<img file="FDA0001126347360000021.GIF" wi="83" he="94" />是行子集选择矩阵的秩,<img file="FDA0001126347360000022.GIF" wi="83" he="94" />是列子集选择矩阵的秩;<img file="FDA0001126347360000023.GIF" wi="133" he="70" />表示中继自干扰矩阵H<sub>rr</sub>的第n个奇异值的平方,(.)<sup>H</sup>表示共轭转置矩阵,(.)<sup>T</sup>表示转置矩阵;(4d)通过(4b)中的行子集选择矩阵<img file="FDA0001126347360000024.GIF" wi="75" he="68" />和列子集选择矩阵D<sub>Rt</sub>,选择出使中继节点的自干扰项W<sub>r</sub>H<sub>rr</sub>W<sub>t</sub>=0的接收波束成形矩阵W<sub>r</sub>,发射波束成形矩阵W<sub>t</sub>:W<sub>r</sub>=[u<sub>r,k1</sub>u<sub>r,k2</sub>...u<sub>r,ki</sub>...u<sub>r,kn</sub>]<sup>H</sup>,W<sub>t</sub>=[v<sub>r,j1</sub>v<sub>r,j2</sub>...v<sub>r,ji</sub>...v<sub>r,jn</sub>]其中u<sub>r,k1</sub>u<sub>r,k2</sub>...u<sub>r,kn</sub>是H<sub>rr</sub>的左奇异向量,v<sub>r,j1</sub>v<sub>r,j2</sub>...v<sub>r,jn</sub>是H<sub>rr</sub>的右奇异向量,下标ki∈[k1,k2···kn],ji∈[j1,j2···jn],且ki≠ji;(4e)计算双向全双工系统总的信道容量:R=(R<sub>1</sub>+R<sub>2</sub>),其中R<sub>1</sub>是第一源节点端的信道容量,R<sub>1</sub>=log<sub>2</sub>(1+SNIR<sub>1</sub>),R<sub>2</sub>是第二源节点端的信道容量R<sub>2</sub>=log<sub>2</sub>(1+SNIR<sub>2</sub>),SNIR<sub>1</sub>第一源节点的信干扰比,SNIR<sub>2</sub>是第二源节点的信干扰比。其中SNIR<sub>2</sub>,SNIR<sub>1</sub>分别如下式:<maths num="0004"><math><![CDATA[<mrow><msub><mi>SNIR</mi><mn>1</mn></msub><mo>=</mo><mfrac><mrow><msub><mi>P</mi><mn>2</mn></msub><mo>|</mo><mo>|</mo><msubsup><mi>U</mi><mn>1</mn><mi>H</mi></msubsup><msubsup><mi>H</mi><mrow><mn>1</mn><mi>r</mi></mrow><mi>H</mi></msubsup><msub><mi>&beta;W</mi><mi>t</mi></msub><msub><mi>W</mi><mi>r</mi></msub><msub><mi>H</mi><mrow><mn>2</mn><mi>r</mi></mrow></msub><msub><mi>V</mi><mn>2</mn></msub><mo>|</mo><msup><mo>|</mo><mn>2</mn></msup></mrow><mrow><msub><mi>P</mi><mn>1</mn></msub><mo>|</mo><mo>|</mo><msubsup><mi>U</mi><mn>1</mn><mi>H</mi></msubsup><msub><mi>H</mi><mn>11</mn></msub><msub><mi>V</mi><mn>1</mn></msub><mo>|</mo><msup><mo>|</mo><mn>2</mn></msup><mo>+</mo><mo>|</mo><mo>|</mo><msubsup><mi>U</mi><mn>1</mn><mi>H</mi></msubsup><msubsup><mi>H</mi><mrow><mn>1</mn><mi>r</mi></mrow><mi>H</mi></msubsup><msub><mi>&beta;W</mi><mi>t</mi></msub><msub><mi>W</mi><mi>r</mi></msub><mo>|</mo><msup><mo>|</mo><mn>2</mn></msup><msubsup><mi>&sigma;</mi><mi>r</mi><mn>2</mn></msubsup><mo>+</mo><mo>|</mo><mo>|</mo><msubsup><mi>U</mi><mn>1</mn><mi>H</mi></msubsup><mo>|</mo><msup><mo>|</mo><mn>2</mn></msup><msubsup><mi>&sigma;</mi><mn>1</mn><mn>2</mn></msubsup></mrow></mfrac></mrow>]]></math><img file="FDA0001126347360000025.GIF" wi="1221" he="143" /></maths><maths num="0005"><math><![CDATA[<mrow><msub><mi>SNIR</mi><mn>2</mn></msub><mo>=</mo><mfrac><mrow><msub><mi>P</mi><mn>1</mn></msub><mo>|</mo><mo>|</mo><msubsup><mi>U</mi><mn>2</mn><mi>H</mi></msubsup><msubsup><mi>H</mi><mrow><mn>2</mn><mi>r</mi></mrow><mi>H</mi></msubsup><msub><mi>&beta;W</mi><mi>t</mi></msub><msub><mi>W</mi><mi>r</mi></msub><msub><mi>H</mi><mrow><mn>1</mn><mi>r</mi></mrow></msub><msub><mi>V</mi><mn>1</mn></msub><mo>|</mo><msup><mo>|</mo><mn>2</mn></msup></mrow><mrow><msub><mi>P</mi><mn>2</mn></msub><mo>|</mo><mo>|</mo><msubsup><mi>U</mi><mn>2</mn><mi>H</mi></msubsup><msub><mi>H</mi><mn>22</mn></msub><msub><mi>V</mi><mn>2</mn></msub><mo>|</mo><msup><mo>|</mo><mn>2</mn></msup><mo>+</mo><mo>|</mo><mo>|</mo><msubsup><mi>U</mi><mn>2</mn><mi>H</mi></msubsup><msubsup><mi>H</mi><mrow><mn>2</mn><mi>r</mi></mrow><mi>H</mi></msubsup><msub><mi>&beta;W</mi><mi>t</mi></msub><msub><mi>W</mi><mi>r</mi></msub><mo>|</mo><msup><mo>|</mo><mn>2</mn></msup><msubsup><mi>&sigma;</mi><mi>r</mi><mn>2</mn></msubsup><mo>+</mo><mo>|</mo><mo>|</mo><msubsup><mi>U</mi><mn>2</mn><mi>H</mi></msubsup><mo>|</mo><msup><mo>|</mo><mn>2</mn></msup><msubsup><mi>&sigma;</mi><mn>2</mn><mn>2</mn></msubsup></mrow></mfrac></mrow>]]></math><img file="FDA0001126347360000026.GIF" wi="1261" he="143" /></maths>其中P<sub>2</sub>是第二源节点的发射功率,||.||<sup>2</sup>表示F范数的平方,β是放大转发协议系数,P<sub>1</sub>是第一源节点的发射功率,<img file="FDA0001126347360000027.GIF" wi="61" he="66" />是中继节点的噪声功率,<img file="FDA0001126347360000028.GIF" wi="61" he="69" />是第一源节点的噪声功率,<img file="FDA0001126347360000029.GIF" wi="61" he="70" />是第二源节点的噪声功率;(4f)根据(4d)和(4e),利用穷举法得到(4d)中满足中继节点自干扰项W<sub>r</sub>H<sub>rr</sub>W<sub>t</sub>=0的所有W<sub>r</sub>,W<sub>t</sub>的解集合,再按照容量最大化准则从解集合里选择一组最优发射波束成形矩阵<img file="FDA00011263473600000210.GIF" wi="106" he="63" />和最优接收波束成形矩阵<img file="FDA00011263473600000211.GIF" wi="123" he="61" />使系统的容量最大化,即:<maths num="0006"><math><![CDATA[<mrow><mo>(</mo><msubsup><mi>W</mi><mi>t</mi><mrow><mi>o</mi><mi>p</mi><mi>t</mi></mrow></msubsup><mo>,</mo><msubsup><mi>W</mi><mi>r</mi><mrow><mi>o</mi><mi>p</mi><mi>t</mi></mrow></msubsup><mo>)</mo><mo>=</mo><mi>arg</mi><munder><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow><mrow><msub><mi>W</mi><mi>t</mi></msub><mo>,</mo><msub><mi>W</mi><mi>r</mi></msub></mrow></munder><mo>(</mo><mi>R</mi><mo>)</mo><mo>.</mo></mrow>]]></math><img file="FDA00011263473600000212.GIF" wi="557" he="86" /></maths>
地址 710071 陕西省西安市雁塔区太白南路2号