发明名称 基于能效最大化的多用户信能同传系统收发机设计方法
摘要 本发明公开了一种基于能效最大化的多用户信能同传系统收发机设计方法,该方法的核心思想为:利用拉格朗日松弛技术将复杂的基于能效最大化的多用户信能同传系统设计问题分解为一系列具有闭式解的子问题,即对传输总功率约束项进行等价变化,对等价的功率约束引入拉格朗日乘子得到部分拉格朗日函数和相应的对偶问题;然后利用二分法和Dinkelbach方法求解对偶问题得到最优能效;最后基站得到最大能效所对应的预编码向量并通过控制信道将接收功率分裂因子发送给用户,完成信能同传系统的收发机设计,在实现信息与能量传输的同时使得系统的能效达到最大。
申请公布号 CN104821838A 申请公布日期 2015.08.05
申请号 CN201510200542.8 申请日期 2015.04.24
申请人 浙江理工大学 发明人 史清江;彭成;徐伟强
分类号 H04B7/02(2006.01)I;H04B7/04(2006.01)I 主分类号 H04B7/02(2006.01)I
代理机构 杭州求是专利事务所有限公司 33200 代理人 邱启旺
主权项 一种基于能效最大化的多用户信能同传系统收发机设计方法,其特征在于,包括以下步骤:(1)初始化如下变量:总功率约束值P<sub>total</sub>,采集功率约束值e<sub>k</sub>,信干噪比约束值γ<sub>k</sub>,采集电路单元的能量转化效率ζ<sub>k</sub>,k=1,2,...,K;(2)对MU‑MISO无线信能同传系统下的能效最大化问题中的传输总功率约束项<img file="FDA0000705143940000011.GIF" wi="258" he="143" />进行等价变化,即<img file="FDA0000705143940000012.GIF" wi="562" he="213" />然后对等价的功率约束项引入拉格朗日乘子λ得到部分拉格朗日函数:<img file="FDA0000705143940000013.GIF" wi="1771" he="283" />和相应的对偶问题,即:min d(λ)s.t.λ>0其中,对偶函数d(λ)表示为<img file="FDA0000705143940000015.GIF" wi="611" he="110" /><maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>s</mi><mo>.</mo><mi>t</mi><mo>.</mo><mfrac><mrow><msub><mi>&rho;</mi><mi>k</mi></msub><msub><mi>p</mi><mi>k</mi></msub><msub><mi>g</mi><mi>k</mi></msub></mrow><mrow><msub><mi>&rho;</mi><mi>k</mi></msub><msubsup><mi>&sigma;</mi><mi>k</mi><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>&delta;</mi><mi>k</mi><mn>2</mn></msubsup></mrow></mfrac><mo>&GreaterEqual;</mo><msub><mi>&gamma;</mi><mi>k</mi></msub><mo>,</mo><mi>k</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mi>K</mi><mo>,</mo></mrow>]]></math><img file="FDA0000705143940000016.GIF" wi="701" he="137" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>&zeta;</mi><mi>k</mi></msub><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>&rho;</mi><mi>k</mi></msub><mo>)</mo></mrow><mrow><mo>(</mo><msub><mi>p</mi><mi>k</mi></msub><msub><mi>g</mi><mi>k</mi></msub><mo>+</mo><msubsup><mi>&sigma;</mi><mi>k</mi><mn>2</mn></msubsup><mo>)</mo></mrow><mo>&GreaterEqual;</mo><msub><mi>e</mi><mi>k</mi></msub><mo>,</mo><mi>k</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mi>K</mi><mo>,</mo></mrow>]]></math><img file="FDA0000705143940000017.GIF" wi="875" he="96" /></maths>0≤ρ<sub>k</sub>≤1,k=1,2,...K.根据Dinkelbach思想,将d(λ)转化为减式形式,并将相应问题分解为K个独立子问题,第k个独立子问题如下:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><munder><mi>max</mi><mrow><msub><mi>p</mi><mi>k</mi></msub><mo>,</mo><msub><mi>&rho;</mi><mi>k</mi></msub></mrow></munder><msub><mi>R</mi><mi>k</mi></msub><mrow><mo>(</mo><msub><mi>p</mi><mi>k</mi></msub><mo>,</mo><msub><mi>&rho;</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>-</mo><mi>&lambda;</mi><mrow><mo>(</mo><msub><mi>p</mi><mi>k</mi></msub><mo>-</mo><mfrac><msub><mi>P</mi><mi>total</mi></msub><mi>K</mi></mfrac><mo>)</mo></mrow><mo>-</mo><mi>&eta;</mi><mrow><mo>(</mo><msub><mi>&theta;p</mi><mi>k</mi></msub><mo>+</mo><mfrac><msub><mi>P</mi><mi>C</mi></msub><mi>K</mi></mfrac><mo>-</mo><msub><mi>E</mi><mi>k</mi></msub><mrow><mo>(</mo><msub><mi>p</mi><mi>k</mi></msub><mo>,</mo><msub><mi>&rho;</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000705143940000021.GIF" wi="1189" he="149" /></maths><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>s</mi><mo>.</mo><mi>t</mi><mo>.</mo><mfrac><mrow><msub><mi>&rho;</mi><mi>k</mi></msub><msub><mi>p</mi><mi>k</mi></msub><msub><mi>g</mi><mi>k</mi></msub></mrow><mrow><msub><mi>&rho;</mi><mi>k</mi></msub><msubsup><mi>&sigma;</mi><mi>k</mi><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>&delta;</mi><mi>k</mi><mn>2</mn></msubsup></mrow></mfrac><mo>&GreaterEqual;</mo><msub><mi>&gamma;</mi><mi>k</mi></msub></mrow>]]></math><img file="FDA0000705143940000022.GIF" wi="401" he="133" /></maths><maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>&zeta;</mi><mi>k</mi></msub><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>&rho;</mi><mi>k</mi></msub><mo>)</mo></mrow><mrow><mo>(</mo><msub><mi>p</mi><mi>k</mi></msub><msub><mi>g</mi><mi>k</mi></msub><mo>+</mo><msubsup><mi>&sigma;</mi><mi>k</mi><mn>2</mn></msubsup><mo>)</mo></mrow><mo>&GreaterEqual;</mo><msub><mi>e</mi><mi>k</mi></msub></mrow>]]></math><img file="FDA0000705143940000023.GIF" wi="562" he="95" /></maths>0<p<sub>k</sub><P<sub>total</sub>0≤ρ<sub>k</sub>≤1其中:R<sub>k</sub>(p<sub>k</sub>,ρ<sub>k</sub>)为用户k的信道容量,E<sub>k</sub>(p<sub>k</sub>,ρ<sub>k</sub>)为用户k的采集功率,ρ<sub>k</sub>为用户k的功率分裂因子,K为总的用户数,P<sub>C</sub>为系统固定总消耗功率,<img file="FDA0000705143940000024.GIF" wi="69" he="82" />为天线引入的加性噪声方差,<img file="FDA0000705143940000025.GIF" wi="69" he="81" />为射频信号转变为基带信号进行信号处理时引起的加性噪声方差,并定义<img file="FDA0000705143940000026.GIF" wi="292" he="105" />h<sub>k</sub>表示基站到用户k的信道向量,<img file="FDA0000705143940000027.GIF" wi="65" he="73" />表示基站到用户k的最优迫零预编码的方向向量,η表示能效值,ζ<sub>k</sub>表示采集电路单元的能量转化效率,<img file="FDA0000705143940000028.GIF" wi="61" he="137" />表示功率放大器效率;(3)设λ=0,求解K个独立子问题,得到对应的{p<sub>k</sub>(λ),ρ<sub>k</sub>(λ)},判断是否满足总功率约束,即<img file="FDA0000705143940000029.GIF" wi="353" he="140" />如果满足,输出{p<sub>k</sub>(λ),ρ<sub>k</sub>(λ)}求得最大能效值;否则执行下一步;(4)设λ=λ+L,其中L为步长,求解K个独立子问题输出{p<sub>k</sub>(λ),ρ<sub>k</sub>(λ)},重复该步骤直至找到满足总功率约束条件的拉格朗日乘子λ,输出拉格朗日乘子上界λ<sub>u</sub>=λ;(5)利用二分法思想求解拉格朗日乘子,即令<img file="FDA00007051439400000210.GIF" wi="263" he="133" />其中λ<sub>l</sub>=0为拉格朗日乘子下界,求解K个独立子问题得到{p<sub>k</sub>(λ),ρ<sub>k</sub>(λ)},判断是否满足总功率约束条件,如满足则令λ<sub>u</sub>=λ,否则令λ<sub>l</sub>=λ,重复该步骤直至<img file="FDA00007051439400000211.GIF" wi="273" he="161" />其中ε为判定阈值,输出{p<sub>k</sub>(λ),ρ<sub>k</sub>(λ)},如果{p<sub>k</sub>(λ),ρ<sub>k</sub>(λ)}满足总功率约束条件,该解即为问题的最优解,因此可求出无线信能同传系统下的最大能效;(6)如果步骤5最终得到的{p<sub>k</sub>(λ),ρ<sub>k</sub>(λ)}不满足总功率约束条件,那么令λ=λ<sub>u</sub>,求解K个独立子问题得到{p<sub>k</sub>(λ),ρ<sub>k</sub>(λ)},该解即为问题的一个可行解;(7)计算传输预编码向量<img file="FDA00007051439400000212.GIF" wi="341" he="91" />k=1,2,...K,基站利用v<sub>k</sub>对传输信号进行预编码,同时通过控制信道将每个功率分裂因子ρ<sub>k</sub>(λ)发送到相应的用户,每个用户设定功率分裂因子,完成信能同传系统的收发机设计,可以进行信息与能量的同时接收。
地址 310018 浙江省杭州市杭州经济技术开发区白杨街道2号大街928号