发明名称 基于能效最优的多输入多输出中继系统联合功率分配方法
摘要 本发明公开了一种基于能效最优的多输入多输出中继系统联合功率分配方法。该通信系统由一个信源节点,一个信宿节点,一个放大转发中继节点所组成,且三个节点均配置多根收发天线。该方法是以最大化系统能效为目标,以指定的最小频谱效率为约束条件,首次建立了以多天线信源节点和多天线中继节点对多条数据流上的发射功率为变量的联合优化数学模型。通过高信噪比区间下的近似转换以及Jensen不等式,将原始的非凸优化问题转化为凸优化问题。再利用拉格朗日对偶函数凸优化算法,并借助于Lambert W函数,首次得到信源节点和中继节点功率变量的闭合形式最优解。
申请公布号 CN105142209A 申请公布日期 2015.12.09
申请号 CN201510591850.8 申请日期 2015.09.17
申请人 东南大学 发明人 李春国;王毅;杨绿溪;王东明;郑福春
分类号 H04W52/24(2009.01)I;H04W52/46(2009.01)I;H04W72/04(2009.01)I 主分类号 H04W52/24(2009.01)I
代理机构 江苏永衡昭辉律师事务所 32250 代理人 王斌
主权项 基于能效最大化的多输入多输出中继系统功率分配方法,其特征在于,所述方法包括以下步骤:1).信源节点通过信道估计获取信道系数矩阵H,中继节点通过信道估计获取信道系数矩阵H和G,信宿节点通过信道估计获取信道系数矩阵G。在三个节点处分别对各自获取的信道系数矩阵进行奇异值分解,即<img file="FDA0000804656910000011.GIF" wi="432" he="102" />和<img file="FDA0000804656910000012.GIF" wi="459" he="107" />其中,<img file="FDA0000804656910000013.GIF" wi="558" he="116" />和<img file="FDA0000804656910000014.GIF" wi="564" he="113" />分别对应于H和G的奇异值;2).信源节点对L条数据流所组成的信号向量z=[z<sub>1</sub>,z<sub>2</sub>,...,z<sub>L</sub>]<sup>T</sup>采用预编码矩阵<img file="FDA0000804656910000015.GIF" wi="358" he="105" />进行发射,即信源节点的发射信号向量为Fz,其中,p<sub>s</sub>=[p<sub>s,1</sub>,p<sub>s,2</sub>,...,p<sub>s,L</sub>]<sup>T</sup>表示信源节点对于L条数据流上分配的功率所组成的向量。信源发射信号经过信道H到达中继节点的信号为y=HFz+n<sub>1</sub>,并采用预编码矩阵<img file="FDA0000804656910000016.GIF" wi="911" he="189" /><img file="FDA0000804656910000017.GIF" wi="75" he="80" />进行发射,即中继转发信号为Wy,其中,p<sub>r</sub>=[p<sub>r,1</sub>,p<sub>r,2</sub>,...,p<sub>r,L</sub>]<sup>T</sup>表示中继节点对于L条数据流上分配的功率所组成的向量。最后,中继转发的信号经过信道G后到达信宿节点。在信宿节点处,利用解码处理矩阵为<img file="FDA0000804656910000019.GIF" wi="202" he="84" />对接收信号进行解码操作,获得解码后的数据流信号向量d=QGWHFz+QGWn<sub>1</sub>+Qn<sub>2</sub>。3).在中继节点处,以最大化系统能效为准则,以指定的最小系统频谱效率为约束,联合求解以信源节点和中继节点对多条数据流上发射功率为参量的功率分配优化问题,如下:<img file="FDA0000804656910000018.GIF" wi="1444" he="309" />其中,p<sub>s</sub>=[p<sub>s,1</sub>,...,p<sub>s,L</sub>]<sup>T</sup>,p<sub>r</sub>=[p<sub>r,1</sub>,...,p<sub>r,L</sub>]<sup>T</sup><img file="FDA0000804656910000021.GIF" wi="1469" he="308" />4).由于3)中优化问题非凸问题,此处利用大信噪比条件,舍去3)中目标函数分子中所包含的噪声功率项<img file="FDA0000804656910000022.GIF" wi="147" he="82" />并利用变量代换x<sub>q</sub>=1/p<sub>r,q</sub>和y<sub>q</sub>=1/p<sub>s,q</sub>,将原最大化目标函数问题,近似等价的变换为一种最小化目标函数问题,如下:<img file="FDA0000804656910000023.GIF" wi="878" he="331" /><img file="FDA0000804656910000024.GIF" wi="887" he="187" />5).利用拉格朗日对偶优化方法,获得4)中的优化问题的拉格朗日对偶函数<img file="FDA0000804656910000025.GIF" wi="79" he="83" />如下:<img file="FDA0000804656910000026.GIF" wi="988" he="346" />6).将<img file="FDA0000804656910000027.GIF" wi="50" he="85" />对x<sub>l</sub>(l=1,2,...,L)取一阶偏导,并令其为零,得到表达式,如下:<img file="FDA0000804656910000028.GIF" wi="1692" he="332" />7).由于6)中等式右侧是与x<sub>l</sub>无关的,可以得到如下等式关系:<img file="FDA0000804656910000029.GIF" wi="1718" he="165" />8).将<img file="FDA00008046569100000210.GIF" wi="52" he="75" />对y<sub>l</sub>(l=1,2,...,L)取一阶偏导,并令其为零,得到表达式,如下:<img file="FDA0000804656910000031.GIF" wi="1690" he="330" />9).合并6)与8)中等式,可以化简得到如下表达式:<img file="FDA0000804656910000032.GIF" wi="624" he="170" />10).将9)中表达式代入7),化简后可以得到y<sub>l</sub>表达式,如下:<img file="FDA0000804656910000033.GIF" wi="1864" he="357" />11).对8)中等式右侧分母项的第二项进行化简,可以得到如下表达式:<img file="FDA0000804656910000034.GIF" wi="707" he="149" />其中,μ<sub>i</sub>(i=1,2,3,4)为常数参量。12).为了求得功率分配变量的闭合形式解,利用Jensen不等式,即<img file="FDA0000804656910000035.GIF" wi="1000" he="155" />对8)中等式右边分子项进行转化,如下:<img file="FDA0000804656910000036.GIF" wi="1383" he="572" />13).将9),11)和12)代入8)中表达式,化简可以得到如下等式:<img file="FDA0000804656910000037.GIF" wi="1751" he="320" />14).由于<img file="FDA0000804656910000041.GIF" wi="658" he="165" />将13)中等式两边同乘以该项<img file="FDA0000804656910000042.GIF" wi="552" he="169" />并化简得到如下表达式:<img file="FDA0000804656910000043.GIF" wi="1723" he="181" />其中,μ<sub>7</sub>和μ<sub>8</sub>为常数参数。15).对14)中等式两边取以自然常数为底的指数运算,并化简得到如下等式:<img file="FDA0000804656910000044.GIF" wi="1691" he="180" />其中,π<sub>1</sub>为常数参数。16).利用Lambert W函数,可以直接求得y<sub>1</sub>的闭合表达式如下:<img file="FDA0000804656910000045.GIF" wi="641" he="241" />其中,<img file="FDA0000804656910000046.GIF" wi="111" he="68" />表示Lambert W函数,其定义为:关于变量x的方程θ=xe<sup>x</sup>,则关于x的解可以用Lambert W函数表示,即<img file="FDA0000804656910000047.GIF" wi="244" he="72" />17).将16)中表达式代入10)式,可以得到信源节点对每条数据流的最优发射功率闭合解,如下:<img file="FDA0000804656910000048.GIF" wi="911" he="182" />其中,π<sub>2</sub>为常数参数18).将16)中表达式代入9)式,可以得到中继节点对每条数据流的最优发射功率闭合形式解,如下:<img file="FDA0000804656910000049.GIF" wi="1214" he="204" />其中:(·)<sup>T</sup>—表示矩阵的转置运算,(·)<sup>H</sup>—表示矩阵的共轭转置运算, diag{x}—以向量x为对角元素的对角阵,M—信源节点发射机天线数,K—中继节点收发机天线数,N—信宿接受机天线数,EE—系统总能效,R—系统总频谱效率,P<sub>total</sub>—系统总功率消耗,K—信源节点发射端输入数据流数且满足L=min{M,N,K},H—信源节点到中继节点间信道系数矩阵,G—中继节点到信宿节点间信道系数矩阵,U<sub>i</sub>|<sub>i</sub><sub>=</sub><sub>1,2</sub>—列正交矩阵,V<sub>i</sub>|<sub>i</sub><sub>=</sub><sub>1,2</sub>—列正交矩阵,<img file="FDA0000804656910000057.GIF" wi="82" he="84" />—信道矩阵H的奇异值组成的列向量,<img file="FDA0000804656910000058.GIF" wi="81" he="91" />—信道矩阵G的奇异值组成的列向量,<img file="FDA0000804656910000059.GIF" wi="71" he="80" />—中继节点处的加性高斯白噪声功率,<img file="FDA00008046569100000510.GIF" wi="79" he="75" />—信宿节点处的加性高斯白噪声功率,p<sub>s,l</sub>—信源节点发射机在第l条数据流上分配的发射功率,p<sub>r,l</sub>—中继节点处在第l条数据流上分配的发射功率,α—信源节点发射机功率放大器功率转换系数,δ<sub>1</sub>—信源节点发射机每根天线上的静态功耗,φ—信源节点发射机的固定电路功耗,β—中继节点发射机功率放大器功率转换系数,δ<sub>1</sub>—中继节点收发机每根天线上的静态功耗,<img file="FDA00008046569100000511.GIF" wi="49" he="63" />—中继节点发射机的固定电路功耗,r<sub>0</sub>—系统要求的最小频谱效率值,λ—拉格朗日乘子,算法中的常量参数包括:<img file="FDA0000804656910000051.GIF" wi="1110" he="197" /><img file="FDA0000804656910000052.GIF" wi="1046" he="205" /><img file="FDA0000804656910000053.GIF" wi="1208" he="363" /><img file="FDA0000804656910000054.GIF" wi="1171" he="178" /><img file="FDA0000804656910000055.GIF" wi="1756" he="178" /><img file="FDA0000804656910000056.GIF" wi="759" he="71" /><img file="FDA0000804656910000061.GIF" wi="616" he="140" /><img file="FDA0000804656910000062.GIF" wi="960" he="182" /><img file="FDA0000804656910000063.GIF" wi="1798" he="1313" /><img file="FDA0000804656910000064.GIF" wi="1120" he="545" />。
地址 210096 江苏省南京市四牌楼2号