发明名称 一种协作多点联合传输系统在多用户场景下的两步预编码方法
摘要 本发明公开了一种协作多点联合传输系统在多用户场景下的两步预编码方法,首先通过干扰对齐处理,在单个基站端建立基于最小化每个用户到其他各用户的干扰泄漏的优化目标,并采用交替最小化方法,分别得到针对每个用户的第一步预编码矩阵和接收矩阵;其次,将第一步得到的针对每个用户的接收矩阵、第一步预编码矩阵和有用的信道矩阵结合为等效信道矩阵,并在此基础上设置相应的相位旋转矩阵,将每个用户不同天线之间的干扰转化为有用信号,有效利用各基站到该用户的发射功率,并得到第二步预编码矩阵和最终的预编码矩阵。本发明的误码率下降比较明显,加快了误码率曲线的下降趋势,并且系统平均速率也有明显提高。
申请公布号 CN103209051B 申请公布日期 2015.12.02
申请号 CN201310074564.5 申请日期 2013.03.08
申请人 西安交通大学 发明人 任品毅;徐大同;张国梅;孙黎
分类号 H04L1/00(2006.01)I 主分类号 H04L1/00(2006.01)I
代理机构 西安通大专利代理有限责任公司 61200 代理人 汪人和
主权项 一种协作多点联合传输系统在多用户场景下的两步预编码方法,其特征在于,包括以下步骤:1)通过干扰对齐处理,将各用户所受的用户间干扰重叠放置在接收空间的正交空间内,并使用交替最小化方法在基站端迭代优化各用户的第一步预编码矩阵和接收矩阵,使用户能够接收无干扰的有用信号;2)在用户间干扰抑制的基础上,各基站得到针对不同用户的等效信道矩阵,设置相位旋转因子,将各用户所受到的数据流间干扰转化为有用信号,有效利用基站发射功率提高系统性能;所述的干扰对齐处理中,干扰对齐处理的目标表示为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msubsup><mi>&Phi;</mi><mi>k</mi><mi>H</mi></msubsup><msub><mi>H</mi><mrow><mi>k</mi><mo>,</mo><mi>i</mi></mrow></msub><msubsup><mi>W</mi><mrow><mi>l</mi><mo>,</mo><mi>i</mi></mrow><mrow><mo>(</mo><mi>a</mi><mo>)</mo></mrow></msubsup><mo>=</mo><mn>0</mn><mo>,</mo><mo>&ForAll;</mo><mi>l</mi><mo>&NotEqual;</mo><mi>k</mi><mo>,</mo><mi>i</mi><mo>;</mo></mrow>]]></math><img file="FDA0000801250210000011.GIF" wi="550" he="83" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>r</mi><mi>a</mi><mi>n</mi><mi>k</mi><mrow><mo>(</mo><msubsup><mi>&Phi;</mi><mi>k</mi><mi>H</mi></msubsup><msub><mi>H</mi><mrow><mi>k</mi><mo>,</mo><mi>i</mi></mrow></msub><msubsup><mi>W</mi><mrow><mi>k</mi><mo>,</mo><mi>i</mi></mrow><mrow><mo>(</mo><mi>a</mi><mo>)</mo></mrow></msubsup><mo>)</mo></mrow><mo>=</mo><msub><mi>L</mi><mi>k</mi></msub><mo>,</mo><mo>&ForAll;</mo><mi>k</mi><mo>,</mo><mi>i</mi><mo>;</mo></mrow>]]></math><img file="FDA0000801250210000012.GIF" wi="667" he="94" /></maths>其中<img file="FDA0000801250210000013.GIF" wi="106" he="89" />表示第i个基站到第k个用户的第一步预编码矩阵,rank表示获得矩阵的秩;Φ<sub>k</sub>为用户k的接收矩阵,H表示共轭转置;H<sub>k,i</sub>为第i个基站到第k个用户的信道矩阵,L<sub>k</sub>为第k个用户的传输数据流数;采用交替最小化方法得到干扰对齐处理的目标的次优解,建立以下优化目标:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>min</mi><mi>i</mi><mi>m</mi><mi>i</mi><mi>z</mi><mi>e</mi><mi> </mi><msub><mi>U</mi><mrow><mi>I</mi><mi>A</mi></mrow></msub><mo>{</mo><mrow><mo>(</mo><msubsup><mi>W</mi><mrow><mi>l</mi><mo>,</mo><mi>i</mi></mrow><mrow><mo>(</mo><mi>a</mi><mo>)</mo></mrow></msubsup><mo>)</mo></mrow><mo>,</mo><mrow><mo>(</mo><msub><mi>&Phi;</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>}</mo><mo>,</mo><mo>&ForAll;</mo><mi>k</mi><mo>,</mo><mi>i</mi><mo>,</mo><mi>l</mi></mrow>]]></math><img file="FDA0000801250210000014.GIF" wi="750" he="110" /></maths><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>s</mi><mo>.</mo><mi>t</mi><mo>.</mo><msubsup><mi>W</mi><mrow><mi>l</mi><mo>,</mo><mi>i</mi></mrow><mrow><mo>(</mo><mi>a</mi><mo>)</mo><mi>H</mi></mrow></msubsup><msubsup><mi>W</mi><mrow><mi>l</mi><mo>,</mo><mi>i</mi></mrow><mrow><mo>(</mo><mi>a</mi><mo>)</mo></mrow></msubsup><mo>=</mo><msub><mi>I</mi><msub><mi>L</mi><mi>l</mi></msub></msub><mo>,</mo><mi>l</mi><mo>&Element;</mo><mo>{</mo><mn>1</mn><mo>,</mo><mo>...</mo><mo>,</mo><mi>N</mi><mo>}</mo></mrow>]]></math><img file="FDA0000801250210000015.GIF" wi="666" he="92" /></maths><maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msubsup><mi>&Phi;</mi><mi>k</mi><mi>H</mi></msubsup><msub><mi>&Phi;</mi><mi>k</mi></msub><mo>=</mo><msub><mi>I</mi><msub><mi>L</mi><mi>k</mi></msub></msub><mo>,</mo><mi>k</mi><mo>&Element;</mo><mo>{</mo><mn>1</mn><mo>,</mo><mo>...</mo><mo>,</mo><mi>N</mi><mo>}</mo></mrow>]]></math><img file="FDA0000801250210000016.GIF" wi="531" he="86" /></maths>其中<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mi>U</mi><mrow><mi>I</mi><mi>A</mi></mrow></msub><mrow><mo>{</mo><mrow><mrow><mo>(</mo><msubsup><mi>W</mi><mrow><mi>l</mi><mo>,</mo><mi>i</mi></mrow><mrow><mo>(</mo><mi>a</mi><mo>)</mo></mrow></msubsup><mo>)</mo></mrow><mo>,</mo><mrow><mo>(</mo><msub><mi>&Phi;</mi><mi>k</mi></msub><mo>)</mo></mrow></mrow><mo>}</mo></mrow><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msubsup><mrow><mo>||</mo><mrow><msubsup><mi>&Phi;</mi><mi>k</mi><mi>H</mi></msubsup><mrow><mo>(</mo><mrow><munderover><mo>&Sigma;</mo><mrow><mi>l</mi><mo>=</mo><mn>1</mn><mo>,</mo><mi>l</mi><mo>&NotEqual;</mo><mi>k</mi></mrow><mi>N</mi></munderover><mrow><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><mrow><msub><mi>H</mi><mrow><mi>k</mi><mo>,</mo><mi>i</mi></mrow></msub><msubsup><mi>W</mi><mrow><mi>l</mi><mo>,</mo><mi>i</mi></mrow><mrow><mo>(</mo><mi>a</mi><mo>)</mo></mrow></msubsup></mrow></mrow></mrow><mo>)</mo></mrow></mrow><mo>||</mo></mrow><mi>F</mi><mn>2</mn></msubsup><mo>;</mo></mrow>]]></math><img file="FDA0000801250210000017.GIF" wi="1057" he="182" /></maths>优化目标表示最小化各用户泄漏到其他用户的干扰的总和;所述的步骤1)的执行过程如下:(a)首先根据优化目标初始化各接收矩阵Φ<sub>k</sub>,<img file="FDA0000801250210000021.GIF" wi="110" he="63" />Φ<sub>k</sub>表示第k个用户的接收矩阵;(b)在每个基站逐一优化各用户的第一步预编码矩阵<img file="FDA0000801250210000022.GIF" wi="292" he="84" />将优化目标改写为:<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msubsup><mi>W</mi><mrow><mi>k</mi><mo>,</mo><mi>i</mi></mrow><mrow><mo>(</mo><mi>a</mi><mo>)</mo></mrow></msubsup><mo>=</mo><mi>arg</mi><munder><mi>min</mi><mrow><msubsup><mi>&Phi;</mi><mi>k</mi><mi>H</mi></msubsup><msub><mi>&Phi;</mi><mi>k</mi></msub><mo>=</mo><msub><mi>I</mi><msub><mi>L</mi><mi>k</mi></msub></msub><mo>,</mo><mo>&ForAll;</mo><mi>k</mi></mrow></munder><mi>t</mi><mi>r</mi><mrow><mo>(</mo><mrow><msubsup><mi>W</mi><mrow><mi>k</mi><mo>,</mo><mi>i</mi></mrow><mrow><mrow><mo>(</mo><mi>a</mi><mo>)</mo></mrow><mi>H</mi></mrow></msubsup><mrow><mo>(</mo><mrow><munderover><mo>&Sigma;</mo><mrow><mi>l</mi><mo>=</mo><mn>1</mn><mo>,</mo><mi>l</mi><mo>&NotEqual;</mo><mi>k</mi></mrow><mi>N</mi></munderover><mrow><msubsup><mi>H</mi><mrow><mi>l</mi><mo>,</mo><mi>i</mi></mrow><mi>H</mi></msubsup><msub><mi>&Phi;</mi><mi>l</mi></msub><msubsup><mi>&Phi;</mi><mi>l</mi><mi>H</mi></msubsup><msub><mi>H</mi><mrow><mi>l</mi><mo>,</mo><mi>i</mi></mrow></msub></mrow></mrow><mo>)</mo></mrow><msubsup><mi>W</mi><mrow><mi>k</mi><mo>,</mo><mi>i</mi></mrow><mrow><mo>(</mo><mi>a</mi><mo>)</mo></mrow></msubsup></mrow><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000801250210000023.GIF" wi="1311" he="178" /></maths>其中,arg min(A)表示A能够使该式最小化,求解结果中<img file="FDA0000801250210000024.GIF" wi="113" he="84" />由矩阵<img file="FDA0000801250210000025.GIF" wi="383" he="148" />最小的若干个特征值对应的L<sub>k</sub>特征向量组成,其中tr表示矩阵求迹运算;(c)在每个基站逐一优化用户的接收矩阵<img file="FDA0000801250210000026.GIF" wi="246" he="79" />将优化目标改写为<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><msub><mi>&Phi;</mi><mi>k</mi></msub><mo>=</mo><mi>arg</mi><munder><mi>min</mi><mrow><msubsup><mi>W</mi><mrow><mi>l</mi><mo>,</mo><mi>i</mi></mrow><mrow><mrow><mo>(</mo><mi>a</mi><mo>)</mo></mrow><mi>H</mi></mrow></msubsup><msubsup><mi>W</mi><mrow><mi>l</mi><mo>,</mo><mi>i</mi></mrow><mrow><mo>(</mo><mi>a</mi><mo>)</mo></mrow></msubsup><mo>=</mo><msub><mi>I</mi><msub><mi>L</mi><mi>l</mi></msub></msub><mo>,</mo><mo>&ForAll;</mo><mi>l</mi><mo>,</mo><mi>i</mi></mrow></munder><mi>t</mi><mi>r</mi><mrow><mo>(</mo><mrow><msubsup><mi>&Phi;</mi><mi>k</mi><mi>H</mi></msubsup><mrow><mo>(</mo><mrow><munderover><mo>&Sigma;</mo><mrow><mi>l</mi><mo>=</mo><mn>1</mn><mo>,</mo><mi>l</mi><mo>&NotEqual;</mo><mi>k</mi></mrow><mi>N</mi></munderover><mrow><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><mrow><msub><mi>H</mi><mrow><mi>k</mi><mo>,</mo><mi>i</mi></mrow></msub><msubsup><mi>W</mi><mrow><mi>l</mi><mo>,</mo><mi>i</mi></mrow><mrow><mo>(</mo><mi>a</mi><mo>)</mo></mrow></msubsup><msubsup><mi>W</mi><mrow><mi>l</mi><mo>,</mo><mi>i</mi></mrow><mrow><mrow><mo>(</mo><mi>a</mi><mo>)</mo></mrow><mi>H</mi></mrow></msubsup><msubsup><mi>H</mi><mrow><mi>k</mi><mo>,</mo><mi>i</mi></mrow><mi>H</mi></msubsup></mrow></mrow></mrow><mo>)</mo></mrow><msub><mi>&Phi;</mi><mi>k</mi></msub></mrow><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000801250210000027.GIF" wi="1343" he="172" /></maths>求解结果中Φ<sub>k</sub>由矩阵<img file="FDA0000801250210000028.GIF" wi="642" he="144" />最小的若干个特征值对应的L<sub>k</sub>个特征向量组成;(d)重复过程(b)、(c),直至相邻两次迭代后的目标函数<img file="FDA0000801250210000029.GIF" wi="499" he="173" />的差值小于预先设定的阈值η或达到最大迭代次数规定;在步骤(c)中求解接收矩阵时,第i个基站可以预先计算出<img file="FDA00008012502100000210.GIF" wi="486" he="146" />的值,并将该值作为等效信道状态信息在协作簇内共享;在优化各用户的第一步预编码矩阵和接收矩阵后,步骤2)的执行过程如下:(a)基于第k个用户的第一步预编码矩阵<img file="FDA0000801250210000031.GIF" wi="135" he="90" />有用信道矩阵H<sub>k,i</sub>和接收矩阵Φ<sub>k</sub>,生成第i个基站到第k个用户的等效信道矩阵D<sub>k,i</sub>:<maths num="0009" id="cmaths0009"><math><![CDATA[<mrow><msub><mi>D</mi><mrow><mi>k</mi><mo>,</mo><mi>i</mi></mrow></msub><mo>=</mo><msubsup><mi>&Phi;</mi><mi>k</mi><mi>H</mi></msubsup><msub><mi>H</mi><mrow><mi>k</mi><mo>,</mo><mi>i</mi></mrow></msub><msubsup><mi>W</mi><mrow><mi>k</mi><mo>,</mo><mi>i</mi></mrow><mrow><mo>(</mo><mi>a</mi><mo>)</mo></mrow></msubsup><mo>,</mo><mo>&ForAll;</mo><mi>k</mi><mo>,</mo><mi>i</mi></mrow>]]></math><img file="FDA0000801250210000032.GIF" wi="501" he="84" /></maths>D<sub>k,i</sub>表示第i个基站到第k个用户的等效信道矩阵;(b)根据等效信道矩阵D<sub>k,i</sub>和信号向量s<sub>k</sub>的相关信息,在基站端计算每个用户的各个数据流的相位旋转因子<img file="FDA0000801250210000033.GIF" wi="307" he="93" />为降低相位旋转处理带来的复杂度,需要根据实际信道情况确定相位旋转因子的取值;定义等效信道矩阵D<sub>k,i</sub>中元素的比例因子<img file="FDA0000801250210000034.GIF" wi="134" he="89" /><maths num="0010" id="cmaths0010"><math><![CDATA[<mrow><msubsup><mi>&mu;</mi><mrow><mi>p</mi><mo>,</mo><mi>q</mi></mrow><mrow><mo>(</mo><mi>k</mi><mo>,</mo><mi>i</mi><mo>)</mo></mrow></msubsup><mo>=</mo><mfrac><mrow><mo>|</mo><msubsup><mi>d</mi><mrow><mi>p</mi><mo>,</mo><mi>q</mi></mrow><mrow><mo>(</mo><mi>k</mi><mo>,</mo><mi>i</mi><mo>)</mo></mrow></msubsup><mi>c</mi><mi>o</mi><mi>n</mi><mi>j</mi><mrow><mo>(</mo><msubsup><mi>d</mi><mrow><mi>p</mi><mo>,</mo><mi>q</mi></mrow><mrow><mo>(</mo><mi>k</mi><mo>,</mo><mi>i</mi><mo>)</mo></mrow></msubsup><mo>)</mo></mrow><mo>|</mo></mrow><mrow><msubsup><mi>&Sigma;</mi><mrow><mi>r</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>L</mi><mi>k</mi></msub></msubsup><mo>|</mo><msubsup><mi>d</mi><mrow><mi>p</mi><mo>,</mo><mi>r</mi></mrow><mrow><mo>(</mo><mi>k</mi><mo>,</mo><mi>i</mi><mo>)</mo></mrow></msubsup><mi>c</mi><mi>o</mi><mi>n</mi><mi>j</mi><mrow><mo>(</mo><msubsup><mi>d</mi><mrow><mi>p</mi><mo>,</mo><mi>r</mi></mrow><mrow><mo>(</mo><mi>k</mi><mo>,</mo><mi>i</mi><mo>)</mo></mrow></msubsup><mo>)</mo></mrow><mo>|</mo></mrow></mfrac><mo>,</mo><mi>p</mi><mo>,</mo><mi>q</mi><mo>&Element;</mo><mo>{</mo><mn>1</mn><mo>,</mo><mo>...</mo><mo>,</mo><msub><mi>L</mi><mi>k</mi></msub><mo>}</mo></mrow>]]></math><img file="FDA0000801250210000035.GIF" wi="950" he="219" /></maths>其中<img file="FDA0000801250210000036.GIF" wi="108" he="95" />表示D<sub>k,i</sub>的第p行第q列的元素,<img file="FDA0000801250210000037.GIF" wi="106" he="96" />表示D<sub>k,i</sub>的第p行第r列的元素,<img file="FDA0000801250210000038.GIF" wi="80" he="92" />和<img file="FDA0000801250210000039.GIF" wi="84" he="90" />分别表示s<sub>k</sub>的第q个和第p个元素,conj表示取共轭运算,| |表示取模值运算;设置第p行的阈值δ,将<img file="FDA00008012502100000310.GIF" wi="100" he="93" />与δ比较:当<maths num="0011" id="cmaths0011"><math><![CDATA[<mrow><msubsup><mi>&mu;</mi><mrow><mi>p</mi><mo>,</mo><mi>q</mi></mrow><mrow><mo>(</mo><mi>k</mi><mo>,</mo><mi>i</mi><mo>)</mo></mrow></msubsup><mo>&lt;</mo><mi>&delta;</mi></mrow>]]></math><img file="FDA00008012502100000311.GIF" wi="200" he="92" /></maths>时<img file="FDA00008012502100000312.GIF" wi="214" he="96" />当<maths num="0012" id="cmaths0012"><math><![CDATA[<mrow><msubsup><mi>&mu;</mi><mrow><mi>p</mi><mo>,</mo><mi>q</mi></mrow><mrow><mo>(</mo><mi>k</mi><mo>,</mo><mi>i</mi><mo>)</mo></mrow></msubsup><mo>&GreaterEqual;</mo><mi>&delta;</mi></mrow>]]></math><img file="FDA00008012502100000313.GIF" wi="188" he="92" /></maths>时<img file="FDA00008012502100000314.GIF" wi="540" he="105" /><img file="FDA00008012502100000315.GIF" wi="106" he="93" />代表数据<img file="FDA00008012502100000316.GIF" wi="84" he="94" />相对于<img file="FDA00008012502100000317.GIF" wi="82" he="93" />所做的相位旋转;(c)在各基站计算针对每个用户的相位旋转矩阵<img file="FDA00008012502100000318.GIF" wi="316" he="87" /><img file="FDA00008012502100000319.GIF" wi="772" he="272" /><img file="FDA0000801250210000041.GIF" wi="94" he="85" />表示第i个基站和第k个用户之间的相位旋转矩阵;(d)利用最小均方误差准则,得到第i个基站到第k个用户的MMSE生成矩阵M<sub>k,i</sub>:<maths num="0013" id="cmaths0013"><math><![CDATA[<mrow><msub><mi>M</mi><mrow><mi>k</mi><mo>,</mo><mi>i</mi></mrow></msub><mo>=</mo><msubsup><mi>D</mi><mrow><mi>k</mi><mo>,</mo><mi>i</mi></mrow><mi>H</mi></msubsup><msup><mrow><mo>(</mo><msub><mi>D</mi><mrow><mi>k</mi><mo>,</mo><mi>i</mi></mrow></msub><msubsup><mi>D</mi><mrow><mi>k</mi><mo>,</mo><mi>i</mi></mrow><mi>H</mi></msubsup><mo>+</mo><mfrac><mrow><msub><mi>L</mi><mi>k</mi></msub><msub><mi>N</mi><mn>0</mn></msub></mrow><msubsup><mi>P</mi><mi>t</mi><mrow><mo>(</mo><mi>k</mi><mo>,</mo><mi>i</mi><mo>)</mo></mrow></msubsup></mfrac><msub><mi>I</mi><msub><mi>L</mi><mi>k</mi></msub></msub><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>,</mo><mo>&ForAll;</mo><mi>k</mi><mo>,</mo><mi>i</mi></mrow>]]></math><img file="FDA0000801250210000042.GIF" wi="813" he="178" /></maths>其中N<sub>0</sub>表示加性高斯噪声的功率谱密度,<img file="FDA0000801250210000043.GIF" wi="112" he="87" />表示第i个基站到第k个用户的发射功率,符号<sub>‑1</sub>表示矩阵求逆运算;(e)计算第i个基站到第k个用户的第二步预编码矩阵<img file="FDA0000801250210000044.GIF" wi="315" he="85" />其中<maths num="0014" id="cmaths0014"><math><![CDATA[<mrow><msubsup><mi>W</mi><mrow><mi>k</mi><mo>,</mo><mi>i</mi></mrow><mrow><mo>(</mo><mi>b</mi><mo>)</mo></mrow></msubsup><mo>=</mo><msub><mi>M</mi><mrow><mi>k</mi><mo>,</mo><mi>i</mi></mrow></msub><msubsup><mi>R</mi><mrow><mi>k</mi><mo>,</mo><mi>i</mi></mrow><mi>&phi;</mi></msubsup><mo>;</mo></mrow>]]></math><img file="FDA0000801250210000045.GIF" wi="344" he="93" /></maths>(f)计算最终的预编码矩阵W<sub>k,i</sub>,<img file="FDA0000801250210000046.GIF" wi="528" he="90" />W<sub>k,i</sub>表示第i个基站到第k个用户的预编码矩阵。
地址 710049 陕西省西安市咸宁西路28号