发明名称 频分复用多用户MIMO能效优化方法
摘要 本发明提出了频分复用多用户MIMO能效优化方法,能够运用于多用户MIMO系统。优化的目的是为每一个用户设计预编码方案,使总能效效用函数最大。优化方法按如下步骤进行:首先将各个用户的能效函数进行优化预处理,具体分为三个步骤:噪声白化与信道对角化,求出链路功率分配方案的注水解,获取优化后的能效函数;然后利用预优化的能效函数,确定各个用户的发送功率以最大化多用户的总能效效用;最后使用用户组能效效用最大化算法解决这个问题。本发明定义的多用户能效效用函数的定义具有一般性,能够涵盖最大化各用户能效的加权平均、幂运算的乘积、加权调和平均、加权最小值等多种指标,且随后提出的高效算法对这些指标均适用。
申请公布号 CN103997775B 申请公布日期 2017.03.15
申请号 CN201410243380.1 申请日期 2014.06.03
申请人 东南大学 发明人 王家恒;赵立成;梁霄
分类号 H04W52/02(2009.01)I;H04B7/0413(2017.01)I;H04Q11/02(2006.01)I 主分类号 H04W52/02(2009.01)I
代理机构 南京瑞弘专利商标事务所(普通合伙) 32249 代理人 杨晓玲
主权项 一种频分复用多用户MIMO能效优化方法,其特征在于,该方法包括如下步骤:1)噪声白化与预编码矩阵换元:根据下式进行信道噪声白化,得到所有用户的等效白化信道矩阵:<maths num="0001"><math><![CDATA[<mrow><msub><mover><mi>H</mi><mo>~</mo></mover><mi>i</mi></msub><mo>=</mo><msubsup><mi>R</mi><mi>i</mi><mrow><mo>-</mo><mfrac><mn>1</mn><mn>2</mn></mfrac></mrow></msubsup><msub><mi>H</mi><mi>i</mi></msub><mo>;</mo></mrow>]]></math><img file="FDA0001172961310000011.GIF" wi="262" he="71" /></maths>其中,i表示用户编号,i∈{1,2,…,K},<img file="FDA0001172961310000012.GIF" wi="54" he="69" />表示用户i的等效白化信道矩阵,K为用户总数,R<sub>i</sub>为用户i的噪声协方差矩阵,H<sub>i</sub>为用户i的信道矩阵;然后进行信道对角化,即根据下式对等效白化信道矩阵<img file="FDA0001172961310000013.GIF" wi="56" he="68" />做奇异值分解:<maths num="0002"><math><![CDATA[<mrow><msub><mover><mi>H</mi><mo>~</mo></mover><mi>i</mi></msub><mover><mo>=</mo><mi>&Delta;</mi></mover><msub><mi>U</mi><mi>i</mi></msub><msub><mi>&Sigma;</mi><mi>i</mi></msub><msubsup><mi>V</mi><mi>i</mi><mi>H</mi></msubsup></mrow>]]></math><img file="FDA0001172961310000014.GIF" wi="274" he="71" /></maths>其中U<sub>i</sub>为用户i的等效白化信道矩阵<img file="FDA0001172961310000015.GIF" wi="57" he="71" />的左奇异矩阵,Σ<sub>i</sub>为用户i的等效白化信道矩阵<img file="FDA0001172961310000016.GIF" wi="57" he="70" />的奇异值矩阵,<img file="FDA0001172961310000017.GIF" wi="70" he="62" />为用户i的等效白化信道矩阵<img file="FDA0001172961310000018.GIF" wi="57" he="70" />的右奇异矩阵的共轭转置,<img file="FDA0001172961310000019.GIF" wi="38" he="46" />为数学中用于定义的符号;最后将修正预编码矩阵<img file="FDA00011729613100000110.GIF" wi="165" he="63" />换元为紧凑预编码矩阵<img file="FDA00011729613100000111.GIF" wi="76" he="70" />其中Q<sub>i</sub>为用户i的简单预编码矩阵,V<sub>i</sub>为用户i的等效白化信道矩阵<img file="FDA00011729613100000112.GIF" wi="57" he="70" />的右奇异矩阵,<img file="FDA00011729613100000113.GIF" wi="70" he="62" />为用户i的等效白化信道矩阵<img file="FDA00011729613100000114.GIF" wi="58" he="70" />的右奇异矩阵的共轭转置;2)运用拉格朗日乘值法求解链路功率分配问题,得到最优功率分配方案为:<maths num="0003"><math><![CDATA[<mrow><msubsup><mi>p</mi><mrow><mi>i</mi><mi>j</mi></mrow><mo>*</mo></msubsup><mo>=</mo><mi>m</mi><mi>a</mi><mi>x</mi><mo>{</mo><mfrac><mn>1</mn><msub><mi>&upsi;</mi><mi>i</mi></msub></mfrac><mo>-</mo><mfrac><mn>1</mn><msubsup><mi>&gamma;</mi><mrow><mi>i</mi><mi>j</mi></mrow><mn>2</mn></msubsup></mfrac><mo>,</mo><mn>0</mn><mo>}</mo></mrow>]]></math><img file="FDA00011729613100000115.GIF" wi="430" he="142" /></maths>其中,<img file="FDA00011729613100000116.GIF" wi="54" he="78" />为用户i的第j个链路的最优功率分配方案,j为链路标号,max表示取最大值,υ<sub>i</sub>为注水水平,υ<sub>i</sub>的取值须满足<img file="FDA00011729613100000117.GIF" wi="235" he="135" />P′<sub>i</sub>是链路功率的总和,γ<sub>ij</sub>为Σ<sub>i</sub>对角线上的元素,即等效白化信道矩阵<img file="FDA00011729613100000118.GIF" wi="58" he="71" />的奇异值,且<img file="FDA00011729613100000119.GIF" wi="541" he="85" />L<sub>i</sub>为等效白化信道矩阵<img file="FDA00011729613100000120.GIF" wi="59" he="67" />的秩;3)根据下式获取预优化能效函数<img file="FDA00011729613100000121.GIF" wi="147" he="63" /><maths num="0004"><math><![CDATA[<mrow><msubsup><mi>u</mi><mi>i</mi><mo>*</mo></msubsup><mrow><mo>(</mo><msub><mi>P</mi><mi>i</mi></msub><mo>)</mo></mrow><mover><mo>=</mo><mi>&Delta;</mi></mover><munder><mi>sup</mi><mrow><msubsup><mi>P</mi><mi>i</mi><mo>&prime;</mo></msubsup><mo>&Element;</mo><mrow><mo>&lsqb;</mo><mrow><mn>0</mn><mo>,</mo><msub><mi>P</mi><mi>i</mi></msub></mrow><mo>&rsqb;</mo></mrow></mrow></munder><msub><mi>u</mi><mi>i</mi></msub><mrow><mo>(</mo><msubsup><mi>P</mi><mi>i</mi><mo>&prime;</mo></msubsup><mo>)</mo></mrow><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>u</mi><mi>i</mi></msub><mrow><mo>(</mo><msub><mi>P</mi><mi>i</mi></msub><mo>)</mo></mrow></mrow></mtd><mtd><mrow><mn>0</mn><mo>&le;</mo><msub><mi>P</mi><mi>i</mi></msub><mo>&lt;</mo><msubsup><mi>P</mi><mi>i</mi><mo>*</mo></msubsup></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>u</mi><mi>i</mi></msub><mrow><mo>(</mo><msubsup><mi>P</mi><mi>i</mi><mo>*</mo></msubsup><mo>)</mo></mrow></mrow></mtd><mtd><mrow><msub><mi>P</mi><mi>i</mi></msub><mo>&GreaterEqual;</mo><msubsup><mi>P</mi><mi>i</mi><mo>*</mo></msubsup></mrow></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0001172961310000021.GIF" wi="913" he="158" /></maths>其中,P<sub>i</sub>为用户i的功率资源,u<sub>i</sub>(P<sub>i</sub>)为用户能效效用,<img file="FDA0001172961310000022.GIF" wi="59" he="70" />为用户i的最优功率资源,<img file="FDA0001172961310000023.GIF" wi="132" he="66" />为最优用户能效效用;4)用户组能效效用的数据初始化:初始化λ<sub>min</sub>=0,<img file="FDA0001172961310000024.GIF" wi="621" he="103" />ε&gt;0;其中λ<sub>min</sub>为拉格朗日乘子λ的下界,λ<sub>max</sub>为拉格朗日乘子λ的上界,ε为迭代精确度,s<sub>i</sub>为用户i对应的权重,γ<sub>i1</sub>为用户i的等效白化信道矩阵<img file="FDA0001172961310000025.GIF" wi="59" he="71" />的最大奇异值,β为基站的静态功率;5)首先根据二分法的折半查找原理,按照下式计算拉格朗日乘子:λ=(λ<sub>min</sub>+λ<sub>max</sub>)/2;然后根据下式求取所有用户的当前最佳功率分配方案为:<maths num="0005"><math><![CDATA[<mrow><msubsup><mi>P</mi><mi>i</mi><mrow><mi>o</mi><mi>p</mi><mi>t</mi></mrow></msubsup><mrow><mo>(</mo><mi>&lambda;</mi><mo>)</mo></mrow><mo>=</mo><mi>m</mi><mi>i</mi><mi>n</mi><mo>{</mo><msub><mi>P</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>&lambda;</mi><mo>)</mo></mrow><mo>,</mo><msub><mi>P</mi><mrow><mi>i</mi><mi>u</mi></mrow></msub><mo>}</mo></mrow>]]></math><img file="FDA0001172961310000026.GIF" wi="500" he="76" /></maths>其中,<img file="FDA0001172961310000027.GIF" wi="149" he="62" />为用户i的的当前最佳功率分配方案,P<sub>iu</sub>为用户i的上限功率,P<sub>i</sub>为用户i的无功率限制情况下的功率分配方案,其为关于λ的函数,P<sub>i</sub>的值满足<img file="FDA0001172961310000028.GIF" wi="435" he="127" /><img file="FDA0001172961310000029.GIF" wi="134" he="71" />为预优化能效函数<img file="FDA00011729613100000210.GIF" wi="123" he="70" />的一阶导数,a为描述指标的参数;6)更新拉格朗日乘子的上限或下限:如果<img file="FDA00011729613100000211.GIF" wi="325" he="129" />则不更新λ<sub>min</sub>,并根据λ<sub>max</sub>=λ更新λ<sub>max</sub>;如果<img file="FDA00011729613100000212.GIF" wi="326" he="134" />则不更新λ<sub>max</sub>,并根据λ<sub>min</sub>=λ更新λ<sub>min</sub>;其中P<sub>s</sub>为基站的总功率;7)获得最优的功率分配方案:如果λ<sub>max</sub>‑λ<sub>min</sub>≤ε,则将最后计算得到的P<sub>i</sub><sup>opt</sup>(λ)作为用户i的最优功率分配方案,从而得到最优功率分配方案集合,否则返回步骤5)。
地址 211189 江苏省南京市江宁区东南大学路2号