CN104954055B - 一种低复杂度的多用户信能同传系统能效优化方法 - Google Patents

一种低复杂度的多用户信能同传系统能效优化方法 Download PDF

Info

Publication number
CN104954055B
CN104954055B CN201510198147.0A CN201510198147A CN104954055B CN 104954055 B CN104954055 B CN 104954055B CN 201510198147 A CN201510198147 A CN 201510198147A CN 104954055 B CN104954055 B CN 104954055B
Authority
CN
China
Prior art keywords
msubsup
msub
mrow
msup
eta
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201510198147.0A
Other languages
English (en)
Other versions
CN104954055A (zh
Inventor
史清江
彭成
徐伟强
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhejiang Sci Tech University ZSTU
Original Assignee
Zhejiang Sci Tech University ZSTU
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Zhejiang Sci Tech University ZSTU filed Critical Zhejiang Sci Tech University ZSTU
Priority to CN201510198147.0A priority Critical patent/CN104954055B/zh
Publication of CN104954055A publication Critical patent/CN104954055A/zh
Application granted granted Critical
Publication of CN104954055B publication Critical patent/CN104954055B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B7/00—Radio transmission systems, i.e. using radiation field
    • H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
    • H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
    • H04B7/0413—MIMO systems
    • H04B7/0456—Selection of precoding matrices or codebooks, e.g. using matrices antenna weighting
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00—Baseband systems
    • H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/03—Shaping networks in transmitter or receiver, e.g. adaptive shaping networks
    • H04L25/03006—Arrangements for removing intersymbol interference
    • H04L25/03178—Arrangements involving sequence estimation techniques
    • H04L25/03248—Arrangements for operating in conjunction with other apparatus
    • H04L25/03292—Arrangements for operating in conjunction with other apparatus with channel estimation circuitry
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04W—WIRELESS COMMUNICATION NETWORKS
    • H04W52/00—Power management, e.g. TPC [Transmission Power Control], power saving or power classes
    • H04W52/04—TPC
    • H04W52/30—TPC using constraints in the total amount of available transmission power
    • H04W52/34—TPC management, i.e. sharing limited amount of power among users or channels or data types, e.g. cell loading
    • H04W52/346—TPC management, i.e. sharing limited amount of power among users or channels or data types, e.g. cell loading distributing total power among users or channels

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Power Engineering (AREA)
  • Pharmaceuticals Containing Other Organic And Inorganic Compounds (AREA)
  • Feedback Control In General (AREA)

Abstract

本发明公开了一种低复杂度的多用户信能同传系统能效优化方法，包括以下步骤：首先求出每个用户在满足服务质量QoS下的最小传输功率，然后根据最小传输功率值计算总功率分配比例得到每个用户独立的功率约束值，并将用户独立功率约束的信能同传系统能效最大化问题分解为一系列存在闭式解的子问题；再根据Dinkelbach方法，迭代求解这些子问题直至得到最大能效值；最后基站得到最大能效所对应的预编码向量并通过控制信道将接收功率分裂因子发送给用户，从而实现信能同传系统的设计。本发明通过巧妙构思引入用户独立功率约束从而将一个总功率约束的复杂问题分解为存在闭式解的独立子问题，从而降低系统的计算复杂度。

Description

一种低复杂度的多用户信能同传系统能效优化方法
技术领域
本发明涉及绿色无线通信传输技术领域，具体为多用户MISO(Multiple InputSingle Output)无线信能同传系统中的低复杂度能效最大化优化方案设计。
背景技术
随着绿色通信技术在当今通信界的地位越来越重要，寻求一种更加节能和环保的资源变得尤为迫切，因此无线信能同传系统，即能够同时传输信息和能量的系统，受到广大研究者的亲睐，并成为绿色通信中的重点研究领域。
在传统无线设备中通常使用有限电量的电池，这使得系统工作效率大大降低。因此，为了使无线设备能够更好地为人们提供服务，在绿色通信的大背景下，无线信能同传系统技术开始进入人们的视野，这种技术的出现使得无线通信设备能够更好地发挥其作用，同时也能缓解当今环境污染与能源匮乏问题。但在研究开始，大多研究工作多以提升频谱效率为目标，却忽略了绿色通信系统设计中一个重要的性能指标——能量效率，简称能效，在数学上可以表示为系统和速率(总的信道容量)与系统总功率消耗之比，通过文献检索可以发现，文献[D.W.K.Ng,E.S.Lo,and R.Schober,“Robust beamforming forsecurecommunication in systems with wireless information and power transfer,”IEEE Trans.Wireless Commun.,vol.13,no.8,pp.4599-4615,Aug.2013.]虽然针对能效问题进行研究并且设计了相关方法最终实现能效最大化，但是该方法却是一种双层迭代算法，增加了计算复杂度，不适合实际应用。因此，本发明针对这个问题设计了一种更实际的方案，即首先求出每个用户在满足服务质量QoS下的最小传输功率，然后根据最小传输功率值计算总功率分配比例得到每个用户独立的功率约束值，并将用户独立功率约束的信能同传系统能效最大化问题分解为一系列存在闭式解的子问题；再根据Dinkelbach方法，迭代求解这些子问题直至得到最大能效值；最后基站得到最大能效所对应的预编码向量并通过控制信道将接收功率分裂因子发送给用户，从而实现信能同传系统的设计。本发明通过巧妙构思引入用户独立功率约束从而将一个总功率约束的复杂问题分解为存在闭式解的独立子问题，从而降低系统的计算复杂度。
发明内容
本发明的目的在于针对现有技术的不足，提供一种低复杂度的多用户信能同传系统能效优化方法，包括以下步骤：
(1)初始化如下变量：总功率约束值Ptotal，采集功率约束值ek，信干噪比约束值γk，采集电路单元的能量转化效率ζk，k＝1,2,...,K；
(2)求解MU-MISO无线信能同传系统下所有用户在满足服务质量QoS时所需最小的传输功率，具体为：第k个用户根据方程计算其在满足服务质量QoS时的功率分裂因子ρk，即可求得最小的传输功率其中hk为信道估计得到的基站到用户k的信道向量，表示基站到用户k的最优迫零预编码方向向量；同理，其余用户均求得pi,min，i＝1,2,...,k-1,k+1,...,K，K为总的用户数；从而得到所有用户独立的功率约束目标Pk，即其中： 为射频信号转变为基带信号进行信号处理时引起的加性噪声方差；然后根据Dinkelbach思想，将MU-MISO无线信能同传系统下的能效最大化问题变为减式问题，并分解为独立子问题，第k个独立子问题如下：
其中：Rk(pk,ρk)为用户k的信道容量，Ek(pk,ρk)为用户k的采集功率，pk为用户k的发射功率，PC为系统固定总消耗功率，η表示能效值，表示功率放大器效率；
(3)初始化迭代次数n＝0、可行解并根据能效函数计算出对应的能效值η(n)，其中分别表示用户在第n次迭代所求得的传输功率和功率分裂因子，所述能效函数为系统和速率(总的信道容量)与系统总功率消耗之比，即:
(4)更新迭代次数n＝n+1，求解第k个子问题，即求解该子问题所对应的四个一元多次方程，即：
方程一：
其中：
方程二：
其中：
方程三：
其中：
方程四：
其中：
求出对应闭式解集合{pk,ρk}，找出满足可行域0≤ρk≤1和max(l1(ρk),l2(ρk))≤pk≤Pk并且使得目标函数取得最大值时对应的即为用户k所求的解，同理，按照此方式余下用户都找到对应解i＝1,2,...,k-1,k+1,...,K，并计算能效值η(n)；
(5)判断是否满足其中ε为判定阈值，如果满足则输出即为K个子问题的最优解，即最优发射功率和功率分裂因子否则执行步骤(3)；
(6)计算传输预编码向量k＝1,2,...K，基站利用vk对传输信号进行预编码，同时通过控制信道将每个功率分裂因子k＝1,2,...K发送到相应的用户，从而每个用户设定功率分裂因子，实现信息与能量的同时接收，即完成信能同传系统的收发机设计。
本发明的有益效果是，本发明针对多用户信能同传系统能效优化方法，通过巧妙地构思引入用户独立功率约束从而将一个总功率约束的复杂问题分解为存在闭式解的独立子问题，然后根据Dinkelbach方法，迭代求解这些子问题直至得到最大能效值；最后基站得到最大能效所对应的预编码向量并通过控制信道将接收功率分裂因子发送给用户，从而实现信能同传系统的设计。本发明通过巧妙构思引入用户独立功率约束从而将一个总功率约束的复杂问题分解为存在闭式解的独立子问题，从而降低系统的计算复杂度。
附图说明
图1是本发明一个实施例的系统模型图；
图2是本发明一个实施例采用该方法的收敛性仿真；
图3是本发明一个实施例的平均能效与发射功率比较图；
图4是本发明一个实施例的平均能效与发射天线数比较图。
具体实施方式
下面结合附图详细描述本发明，本发明的目的和效果将变得更加明显。
如图1所示，在该系统模型中假设基站(BS)(或接入点)天线数为Nt(Nt＞0)，基站利用传输预编码(或称传输波束赋形向量)vk，k＝1,...,K，向K个单天线接收机传输符号sk，k＝1,...,K，然而不同于传统的多用户MISO系统，用户接受到的信号将分裂两部分，其中一部分用于信息译码，另一部分则进行能量采集。由此信道模型可知，在功率分裂前第k个用户接收到的信号yk为：
其中：hk为基于导频的信道估计方法得到的基站到用户k的信道向量，sk表示用户k的符号，nk为天线引入的加性噪声，服从均值为0方差为的循环对称复高斯分布。
在经过功率分裂以后，信号分成两部分，其中信息译码部分可表示为：
上式中zk为射频信号转变为基带信号进行信号处理时引起的加性噪声，服从均值为0方差为的循环对称复高斯分布，ρk表示接收机k的功率分裂因子。。
第二部分能量采集函数表达式可写为：
因此，信干噪比SINRk(v,ρk)可表示为：
相应的速率函数Rk可表示为：
Rk＝W log(1+SINRk) (5)
其中W表示系统带宽。
而接收机采集功率Ek(v,ρk)，即单位时间内所采集的能量，可表示为：
其中ζk表示采集电路单元的能量转化效率。
再根据能效函数的定义，能效函数η可以表示为：
这里表示系统消耗的总功率，其中表示功率放大器效率，PC表示收发机总固定电路消耗(包括数模转换、频率合成等)，。
对于信能同传系统，通信服务质量包括两种，一种是为保证正常通信需要信干噪比SINRk达到一定的要求，另一种为了维持系统正常工作需要采集功率Ek满足一定的条件。特别地，为了简化发射机，本发明考虑常用的波束成形方案——迫零预编码来消除多用户之间的信号干扰。因此，基于迫零波束形成的多用户信能同传系统能效最大化设计问题可描述为：
其中：γk和ek分别为SINRk和Ek的设计目标，Hk＝[h1,...,hk-1,hk+1,...,hK],k＝1,2,...,K。
定义根据矩阵零空间知识，最优迫零波束成形向量的方向可以表示为： 其中：Uk表示的零空间正交基。
则问题(8)可等价地转化成下面问题
为使问题(9)容易求解，考虑将总功率约束条件变为每个用户自己独立约束，即pk≤PK,那么则原问题就可轻易分解为子问题然后只需要采用Dinkelbach迭代方法求解即可。从这个角度出发，我们将按照一定比例分配每个用户的最大的功率约束，即其中pk,min为第k个用户在满足Qos时所需最小的传输功率，如此一来，问题的关键就转向求解pk,min。
分析原问题，设则可以知道存在因此令l1(ρk)＝l2(ρk)，即可得到关于ρk的2次方程：
其中：所以通过方程(10)就可求出pk,min，那么问题(9)转化为求解下列问题“
并且第k个子问题可以表示为
再根据Dinkelbach迭代方法，首先引入符号η，将问题(12)中的分式问题转换为减式问题，即
然后求解问题(13)，根据已知知识，假设由ρk表示pk的表达式pk(ρk)可导且已知，那么问题(13)关于ρk的最优解必满足下列式子：
其中：0＜ρk＜1，
通过分析问题(13)，最优解pk只能在边界点或者目标函数的静态点取到，即最优值pk(ρk)只能满足以下4种情况：1)pk(ρk)＝Pk；2)pk(ρk)＝l1(ρk)；3)pk(ρk)＝l2(ρk)；4)其中表示问题(13)的目标函数的静态点，即所以只需结合式子(14)依次分析上面所述4种情形，即求得ρk的闭式解，从而求解问题(13)。
1)情形1：根据pk(ρk)＝Pk和式子(14)可得到关于ρk的一元二次方程：
其中：
2)情形2：当pk(ρk)＝l1(ρk)，则有关于ρk的一元二次方程：
其中:
3)情形3：当pk(ρk)＝l2(ρk)，同理有关于ρk的一元4次方程：
其中：
4)情形4：因为由于为目标函数静态点，那么有
可以推出：
则得到关于ρk的一元4次方程：
其中：
因此，所有子问题可以通过以上4个等式，即(15)、(16)、(17)和(20)求出{ρk,pk}，同时检查它们是否满足的可行域，即0≤ρk≤1和max(l1(ρk),l21(ρk))≤pk≤Pk，并从这些解中找到使取得最大值所对应的即问题(13)的解，同理，其余用户都按照此方式求解出对应的传输功率和功率分裂因子，最后迭代求解这些子问题直至得到最大能效值。
因此，一种低复杂度的多用户信能同传系统能效优化方法，包括以下步骤：
(1)初始化如下变量：总功率约束值Ptotal，采集功率约束值ek，信干噪比约束值γk，采集电路单元的能量转化效率ζk，k＝1,2,...,K；
(2)求解MU-MISO无线信能同传系统下所有用户在满足服务质量QoS时所需最小的传输功率，具体为：第k个用户根据方程(10)计算出在满足服务质量QoS时所需最小的传输功率pk,min；同理，其余用户均求得pi,min，i＝1,2,...,k-1,k+1,...,K，从而采集所有用户各自独立的功率约束目标Pk，即根据Dinkelbach迭代方法将MU-MISO无线信能同传系统下的能效最大化问题变为减式问题，并分解为独立子问题，第k个独立子问题如式子(13)；
(3)初始化迭代次数n＝0、可行解并根据能效函数计算出对应的能效值η(n)，其中分别表示用户在第n次迭代所求得的传输功率和功率分裂因子，所述能效函数为系统和速率(总的信道容量)与系统总功率消耗之比，即
(4)更新迭代次数n＝n+1，求解第k个子问题，即求解该子问题所对应的四个一元多次方程，即方式(15)、(16)、(17)和(20)，求出对应闭式解集合{pk,ρk}，找出满足可行域0≤ρk≤1和max(l1(ρk),l2(ρk))≤pk≤Pk并且使得问题(13)目标函数取得最大值时对应的即为用户k所求的解，同理，直到所有子问题都找到对应解i＝1,2,...,k-1,k+1,...,K，并计算能效值η(n)；
(5)判断是否满足其中ε为判定阈值，如果满足则输出即为K个子问题的最优解，即得到最优发射功率和功率分裂因子否则执行步骤(3)；
(6)计算传输预编码向量k＝1,2,...K，基站利用vk对传输信号进行预编码，同时通过控制信道将每个功率分裂因子k＝1,2,...K发送到相应的用户，从而每个用户设定功率分裂因子，实现信息与能量的同时接收，即完成信能同传系统的收发机设计。
图2、图3和图4是本发明通过Matlab对所设计方案的仿真验证。参数具体设置为：发送端天线数Nt＝4，信息接受者的数量K＝4，能量转换系数ξ＝0.65，天线噪声功率传输噪声功率功率放大器效率θ＝5，带宽W＝15KHz，此外,假设所有的信息接收者具有相同的SINRk和Ek的门限值,即γ1＝…＝γK＝γ和e1＝…＝eK＝e，如果不作特殊说明，在仿真中设置总传输功率Ptotal＝30dBm，γ＝20dB，e＝-20dBm，本文的收发机总固定电路消耗PC按如下方式取值：
其中：PDAC，Pmix，Pfilt，Psyn，PLNA，PIFA，PADC分别表示数模转化、混合器、发射机端的滤波器、频率混合器、低噪声放大器、中频放大器、接收机端滤波器以及模数转化所消耗的功率。在仿真中，每个参数的取值如表格1：
表1收发机总固定电路消耗中各参数取值
混合器Pmix 30.3mW
滤波器Pfilt＝Pfilr 2.5mW
频率混合器Psyn 50mW
低噪声放大器PLNA 20mW
中频放大器PIFA 3mW
数模转化PDAC 15.44mW
模数转化PADC 6.76mW
并设蒙特卡洛仿真次数为500000，比较结果分析如下：
图2验证了算法2的收敛性，图中的10条曲线分别对应十个随机信道得到的收敛曲线。从图可以看出该算法可以保证能效值的单调不减直至收敛，并且一般在3～4次取到最优解。体现出该算法的高效性。
图3给出了所设计方案与总传输功率的比较图。图中设置的参数为：γ＝15dB，e＝-25dBm，由图可以看出随着总传输功率的不断增大，平均能效值不断增大，特别地，当总传输功率Ptotal≥26dBm的时候，平均能效值不再增大而保持在同一水平线上，意味着此时总传输功率不再是影响平均能效的重要参数。
图4描绘出了所设计方法与天线数的比较图，图中设置的参数为：γ＝15dB，e＝-20dBm，当然收发机总固定电路消耗PC将随着天线数的增大而增大。由图同样可以看出上界值与本发明所设计方法得到值重合，并且当天线数Nt从4增大到12的时候，平均能效值不断增大，而从12到60，平均能效值却逐渐减小，意味着当发射机端要采用大规模天线阵列的时候应进行天线选择才能获得性能良好的能效。
通过前面的分析与性能仿真比较，本发明的方法不仅计算和通信复杂度低，而且能够优化能效性能，可以预见本发明方法能很好地适应未来绿色移动通信技术。
本发明不仅局限于上述具体实施方式，本领域一般技术人员根据本发明公开的内容，可以采用其它多种具体实施方案实施本发明。因此，凡是采用本发明的设计结构和思路，做一些简单的变化或更改的设计，都落入本发明保护范围。

Claims (1)

1.一种低复杂度的多用户信能同传系统能效优化方法，其特征在于，包括以下步骤：
(1)初始化如下变量：总功率约束值Ptotal，采集功率约束值ek，信干噪比约束值γk，采集电路单元的能量转化效率ζk，k＝1,2,...,K，K为总的用户数；
(2)求解MU-MISO无线信能同传系统下所有用户在满足服务质量QoS时所需最小的传输功率，具体为：第k个用户根据方程计算其在满足服务质量QoS时的功率分裂因子ρk，即可求得最小的传输功率其中hk为信道估计得到的基站到用户k的信道向量，表示基站到用户k的最优迫零预编码方向向量，为天线引入的加性噪声方差；同理，其余用户均求得pi,min，i＝1,2,...,k-1,k+1,...,K；从而得到所有用户独立的功率约束目标Pk，即其中： 为射频信号转变为基带信号进行信号处理时引起的加性噪声方差；然后根据Dinkelbach思想，将MU-MISO无线信能同传系统下的能效最大化问题变为减式问题，并分解为独立子问题，第k个独立子问题如下：
<mrow> <munder> <mi>max</mi> <mrow> <msub> <mi>p</mi> <mi>k</mi> </msub> <mo>,</mo> <msub> <mi>&amp;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>&amp;rho;</mi> <mi>k</mi> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mi>&amp;eta;</mi> <mrow> <mo>(</mo> <msub> <mi>&amp;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> <mo>(</mo> <mrow> <msub> <mi>p</mi> <mi>k</mi> </msub> <mo>,</mo> <msub> <mi>&amp;rho;</mi> <mi>k</mi> </msub> </mrow> <mo>)</mo> <mo>)</mo> </mrow> </mrow>
<mrow> <mtable> <mtr> <mtd> <mrow> <mi>s</mi> <mo>.</mo> <mi>t</mi> <mo>.</mo> </mrow> </mtd> <mtd> <mrow> <mfrac> <mrow> <msub> <mi>&amp;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>&amp;rho;</mi> <mi>k</mi> </msub> <msubsup> <mi>&amp;sigma;</mi> <mi>k</mi> <mn>2</mn> </msubsup> <mo>+</mo> <msubsup> <mi>&amp;delta;</mi> <mi>k</mi> <mn>2</mn> </msubsup> </mrow> </mfrac> <mo>&amp;GreaterEqual;</mo> <msub> <mi>&amp;gamma;</mi> <mi>k</mi> </msub> </mrow> </mtd> </mtr> </mtable> <mo>,</mo> </mrow>
<mrow> <msub> <mi>&amp;zeta;</mi> <mi>k</mi> </msub> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msub> <mi>&amp;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>&amp;sigma;</mi> <mi>k</mi> <mn>2</mn> </msubsup> <mo>)</mo> </mrow> <mo>&amp;GreaterEqual;</mo> <msub> <mi>e</mi> <mi>k</mi> </msub> <mo>,</mo> </mrow>
pk≤Pk,
0≤ρk≤1.
其中：Rk(pk,ρk)为用户k的信道容量，Ek(pk,ρk)为用户k的采集功率，pk为用户k的发射功率，PC为系统固定总消耗功率，η表示能效值，表示功率放大器效率；
(3)初始化迭代次数n＝0、可行解并根据能效函数计算出对应的能效值η(n)，其中分别表示用户在第n次迭代所求得的传输功率和功率分裂因子，所述能效函数为系统和速率(总的信道容量)与系统总功率消耗之比，即:
<mrow> <msup> <mi>&amp;eta;</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </msup> <mo>=</mo> <mfrac> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>K</mi> </munderover> <msub> <mi>R</mi> <mi>k</mi> </msub> <mrow> <mo>(</mo> <msubsup> <mi>p</mi> <mi>k</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </msubsup> <mo>,</mo> <msubsup> <mi>&amp;rho;</mi> <mi>k</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </msubsup> <mo>)</mo> </mrow> </mrow> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>K</mi> </munderover> <msubsup> <mi>&amp;theta;p</mi> <mi>k</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </msubsup> <mo>+</mo> <msub> <mi>P</mi> <mi>C</mi> </msub> <mo>-</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>K</mi> </munderover> <msub> <mi>E</mi> <mi>k</mi> </msub> <mrow> <mo>(</mo> <msubsup> <mi>p</mi> <mi>k</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </msubsup> <mo>,</mo> <msubsup> <mi>&amp;rho;</mi> <mi>k</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </msubsup> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>;</mo> </mrow>
(4)更新迭代次数n＝n+1，求解第k个子问题，即求解该子问题所对应的四个一元多次方程，即：
方程一：
其中：
方程二：
其中：
方程三：
其中：
<mrow> <msub> <mi>C</mi> <mn>3</mn> </msub> <mo>=</mo> <mrow> <mo>(</mo> <mn>6</mn> <msub> <mi>g</mi> <mi>k</mi> </msub> <msubsup> <mi>&amp;zeta;</mi> <mi>k</mi> <mn>2</mn> </msubsup> <msubsup> <mi>&amp;sigma;</mi> <mi>k</mi> <mn>2</mn> </msubsup> <mo>-</mo> <mn>3</mn> <msub> <mi>g</mi> <mi>k</mi> </msub> <msub> <mi>e</mi> <mi>k</mi> </msub> <msub> <mi>&amp;zeta;</mi> <mi>k</mi> </msub> <mo>-</mo> <mn>3</mn> <msub> <mi>g</mi> <mi>k</mi> </msub> <msubsup> <mi>&amp;zeta;</mi> <mi>k</mi> <mn>2</mn> </msubsup> <msubsup> <mi>&amp;delta;</mi> <mi>k</mi> <mn>2</mn> </msubsup> <mo>-</mo> <msubsup> <mi>e</mi> <mi>k</mi> <mn>2</mn> </msubsup> <mi>&amp;eta;</mi> <mi>&amp;theta;</mi> <mo>+</mo> <msub> <mi>e</mi> <mi>k</mi> </msub> <msub> <mi>&amp;theta;&amp;zeta;</mi> <mi>k</mi> </msub> <msubsup> <mi>&amp;delta;</mi> <mi>k</mi> <mn>2</mn> </msubsup> <mi>&amp;eta;</mi> <mo>)</mo> </mrow> <msubsup> <mi>&amp;sigma;</mi> <mi>k</mi> <mn>2</mn> </msubsup> <mo>,</mo> </mrow>
<mrow> <msub> <mi>D</mi> <mn>3</mn> </msub> <mo>=</mo> <msub> <mi>e</mi> <mi>k</mi> </msub> <msub> <mi>g</mi> <mi>k</mi> </msub> <msub> <mi>&amp;zeta;</mi> <mi>k</mi> </msub> <mrow> <mo>(</mo> <mn>2</mn> <msubsup> <mi>&amp;sigma;</mi> <mi>k</mi> <mn>2</mn> </msubsup> <mo>-</mo> <msubsup> <mi>&amp;delta;</mi> <mi>k</mi> <mn>2</mn> </msubsup> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>g</mi> <mi>k</mi> </msub> <msubsup> <mi>&amp;zeta;</mi> <mi>k</mi> <mn>2</mn> </msubsup> <msubsup> <mi>&amp;sigma;</mi> <mi>k</mi> <mn>2</mn> </msubsup> <mrow> <mo>(</mo> <mn>3</mn> <msubsup> <mi>&amp;delta;</mi> <mi>k</mi> <mn>2</mn> </msubsup> <mo>-</mo> <mn>2</mn> <msubsup> <mi>&amp;sigma;</mi> <mi>k</mi> <mn>2</mn> </msubsup> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>e</mi> <mi>k</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>&amp;theta;&amp;zeta;</mi> <mi>k</mi> </msub> <msubsup> <mi>&amp;delta;</mi> <mi>k</mi> <mn>4</mn> </msubsup> <mi>&amp;eta;</mi> <mo>-</mo> <msub> <mi>e</mi> <mi>k</mi> </msub> <msubsup> <mi>&amp;eta;&amp;theta;&amp;delta;</mi> <mi>k</mi> <mn>2</mn> </msubsup> <mo>-</mo> <msub> <mi>&amp;theta;&amp;zeta;</mi> <mi>k</mi> </msub> <msubsup> <mi>&amp;delta;</mi> <mi>k</mi> <mn>2</mn> </msubsup> <msubsup> <mi>&amp;eta;&amp;sigma;</mi> <mi>k</mi> <mn>2</mn> </msubsup> <mo>)</mo> </mrow> <mo>,</mo> </mrow>
<mrow> <msub> <mi>E</mi> <mn>3</mn> </msub> <mo>=</mo> <mrow> <mo>(</mo> <msub> <mi>g</mi> <mi>k</mi> </msub> <msub> <mi>e</mi> <mi>k</mi> </msub> <mo>-</mo> <msub> <mi>g</mi> <mi>k</mi> </msub> <msub> <mi>&amp;zeta;</mi> <mi>k</mi> </msub> <msubsup> <mi>&amp;sigma;</mi> <mi>k</mi> <mn>2</mn> </msubsup> <mo>-</mo> <msub> <mi>e</mi> <mi>k</mi> </msub> <msubsup> <mi>&amp;theta;&amp;delta;</mi> <mi>k</mi> <mn>2</mn> </msubsup> <mi>&amp;eta;</mi> <mo>)</mo> </mrow> <msub> <mi>&amp;zeta;</mi> <mi>k</mi> </msub> <msubsup> <mi>&amp;delta;</mi> <mi>k</mi> <mn>2</mn> </msubsup> <mo>;</mo> </mrow>
方程四：
其中：
<mrow> <msub> <mi>B</mi> <mn>4</mn> </msub> <mo>=</mo> <mrow> <mo>(</mo> <msubsup> <mi>&amp;delta;</mi> <mi>k</mi> <mn>2</mn> </msubsup> <msub> <mi>&amp;zeta;</mi> <mi>k</mi> </msub> <mi>&amp;eta;</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> <msubsup> <mi>g</mi> <mi>k</mi> <mn>4</mn> </msubsup> <msubsup> <mi>&amp;zeta;</mi> <mi>k</mi> <mn>2</mn> </msubsup> <msup> <mi>&amp;eta;</mi> <mn>2</mn> </msup> <msubsup> <mi>&amp;sigma;</mi> <mi>k</mi> <mn>2</mn> </msubsup> <mo>-</mo> <mrow> <mo>(</mo> <msubsup> <mi>&amp;delta;</mi> <mi>k</mi> <mn>2</mn> </msubsup> <msub> <mi>&amp;theta;&amp;zeta;</mi> <mi>k</mi> </msub> <msup> <mi>&amp;eta;</mi> <mn>2</mn> </msup> <mo>+</mo> <mi>&amp;theta;</mi> <mi>&amp;eta;</mi> <mo>)</mo> </mrow> <msubsup> <mi>g</mi> <mi>k</mi> <mn>3</mn> </msubsup> <msubsup> <mi>&amp;sigma;</mi> <mi>k</mi> <mn>2</mn> </msubsup> <msub> <mi>&amp;zeta;</mi> <mi>k</mi> </msub> <mi>&amp;eta;</mi> <mo>,</mo> </mrow>
<mrow> <mtable> <mtr> <mtd> <mrow> <msub> <mi>C</mi> <mn>4</mn> </msub> <mo>=</mo> <mrow> <mo>(</mo> <msub> <mi>&amp;zeta;</mi> <mi>k</mi> </msub> <msubsup> <mi>&amp;delta;</mi> <mi>k</mi> <mn>2</mn> </msubsup> <mi>&amp;eta;</mi> <mo>-</mo> <mn>2</mn> <msubsup> <mi>&amp;sigma;</mi> <mi>k</mi> <mn>2</mn> </msubsup> <msub> <mi>&amp;zeta;</mi> <mi>k</mi> </msub> <mi>&amp;eta;</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <msubsup> <mi>&amp;zeta;</mi> <mi>k</mi> <mn>2</mn> </msubsup> <msubsup> <mi>g</mi> <mi>k</mi> <mn>4</mn> </msubsup> <msup> <mi>&amp;eta;</mi> <mn>2</mn> </msup> <msubsup> <mi>&amp;delta;</mi> <mi>k</mi> <mn>2</mn> </msubsup> <mo>+</mo> <mrow> <mo>(</mo> <mn>4</mn> <msubsup> <mi>&amp;sigma;</mi> <mi>k</mi> <mn>2</mn> </msubsup> <msub> <mi>&amp;theta;&amp;zeta;</mi> <mi>k</mi> </msub> <msup> <mi>&amp;eta;</mi> <mn>2</mn> </msup> <mo>+</mo> <mi>&amp;eta;</mi> <mi>&amp;theta;</mi> <mo>-</mo> <msub> <mi>&amp;zeta;</mi> <mi>k</mi> </msub> <msup> <mi>&amp;eta;</mi> <mn>2</mn> </msup> <msubsup> <mi>&amp;theta;&amp;delta;</mi> <mi>k</mi> <mn>2</mn> </msubsup> <mo>-</mo> <msub> <mi>&amp;zeta;</mi> <mi>k</mi> </msub> <msubsup> <mi>&amp;eta;&amp;theta;&amp;delta;</mi> <mi>k</mi> <mn>2</mn> </msubsup> <mo>)</mo> </mrow> <msubsup> <mi>g</mi> <mi>k</mi> <mn>3</mn> </msubsup> <msubsup> <mi>&amp;delta;</mi> <mi>k</mi> <mn>2</mn> </msubsup> <msub> <mi>&amp;zeta;</mi> <mi>k</mi> </msub> <mi>&amp;eta;</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <mn>2</mn> <msup> <mi>&amp;eta;</mi> <mn>2</mn> </msup> <msup> <mi>&amp;theta;</mi> <mn>2</mn> </msup> <msubsup> <mi>g</mi> <mi>k</mi> <mn>2</mn> </msubsup> <msubsup> <mi>&amp;sigma;</mi> <mi>k</mi> <mn>2</mn> </msubsup> <msub> <mi>&amp;zeta;</mi> <mi>k</mi> </msub> <msubsup> <mi>&amp;eta;&amp;delta;</mi> <mi>k</mi> <mn>2</mn> </msubsup> </mrow> </mtd> </mtr> </mtable> <mo>,</mo> </mrow>
<mrow> <mtable> <mtr> <mtd> <mrow> <msub> <mi>D</mi> <mn>4</mn> </msub> <mo>=</mo> <mrow> <mo>(</mo> <msubsup> <mi>&amp;sigma;</mi> <mi>k</mi> <mn>2</mn> </msubsup> <msub> <mi>&amp;zeta;</mi> <mi>k</mi> </msub> <mi>&amp;eta;</mi> <mo>-</mo> <mn>2</mn> <msubsup> <mi>&amp;delta;</mi> <mi>k</mi> <mn>2</mn> </msubsup> <msub> <mi>&amp;zeta;</mi> <mi>k</mi> </msub> <mi>&amp;eta;</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> <msubsup> <mi>&amp;zeta;</mi> <mi>k</mi> <mn>2</mn> </msubsup> <msubsup> <mi>g</mi> <mi>k</mi> <mn>4</mn> </msubsup> <msup> <mi>&amp;eta;</mi> <mn>2</mn> </msup> <msubsup> <mi>&amp;delta;</mi> <mi>k</mi> <mn>2</mn> </msubsup> <mo>+</mo> <mrow> <mo>(</mo> <mn>4</mn> <msubsup> <mi>&amp;delta;</mi> <mi>k</mi> <mn>2</mn> </msubsup> <msub> <mi>&amp;theta;&amp;zeta;</mi> <mi>k</mi> </msub> <msup> <mi>&amp;eta;</mi> <mn>2</mn> </msup> <mo>-</mo> <mn>3</mn> <msubsup> <mi>&amp;sigma;</mi> <mi>k</mi> <mn>2</mn> </msubsup> <msub> <mi>&amp;theta;&amp;zeta;</mi> <mi>k</mi> </msub> <msup> <mi>&amp;eta;</mi> <mn>2</mn> </msup> <mo>-</mo> <mn>3</mn> <msubsup> <mi>&amp;sigma;</mi> <mi>k</mi> <mn>2</mn> </msubsup> <msub> <mi>&amp;theta;&amp;zeta;</mi> <mi>k</mi> </msub> <mi>&amp;eta;</mi> <mo>-</mo> <mn>2</mn> <mi>&amp;eta;</mi> <mi>&amp;theta;</mi> <mo>)</mo> </mrow> <msubsup> <mi>g</mi> <mi>k</mi> <mn>3</mn> </msubsup> <msubsup> <mi>&amp;delta;</mi> <mi>k</mi> <mn>2</mn> </msubsup> <msub> <mi>&amp;zeta;</mi> <mi>k</mi> </msub> <mi>&amp;eta;</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>+</mo> <mrow> <mo>(</mo> <mn>3</mn> <msubsup> <mi>&amp;sigma;</mi> <mi>k</mi> <mn>2</mn> </msubsup> <msub> <mi>&amp;zeta;</mi> <mi>k</mi> </msub> <msup> <mi>&amp;eta;</mi> <mn>3</mn> </msup> <msup> <mi>&amp;theta;</mi> <mn>2</mn> </msup> <mo>-</mo> <mn>2</mn> <msub> <mi>&amp;zeta;</mi> <mi>k</mi> </msub> <msubsup> <mi>&amp;delta;</mi> <mi>k</mi> <mn>2</mn> </msubsup> <msup> <mi>&amp;eta;</mi> <mn>3</mn> </msup> <msup> <mi>&amp;theta;</mi> <mn>2</mn> </msup> <mo>+</mo> <msup> <mi>&amp;eta;</mi> <mn>2</mn> </msup> <msup> <mi>&amp;theta;</mi> <mn>2</mn> </msup> <mo>)</mo> </mrow> <msubsup> <mi>g</mi> <mi>k</mi> <mn>2</mn> </msubsup> <msubsup> <mi>&amp;delta;</mi> <mi>k</mi> <mn>2</mn> </msubsup> <mo>-</mo> <msup> <mi>&amp;eta;</mi> <mn>3</mn> </msup> <msup> <mi>&amp;theta;</mi> <mn>3</mn> </msup> <msub> <mi>g</mi> <mi>k</mi> </msub> <msubsup> <mi>&amp;sigma;</mi> <mi>k</mi> <mn>2</mn> </msubsup> <msubsup> <mi>&amp;delta;</mi> <mi>k</mi> <mn>2</mn> </msubsup> </mrow> </mtd> </mtr> </mtable> <mo>,</mo> </mrow>
<mrow> <msub> <mi>E</mi> <mn>4</mn> </msub> <mo>=</mo> <mo>-</mo> <msubsup> <mi>g</mi> <mi>k</mi> <mn>4</mn> </msubsup> <msubsup> <mi>&amp;zeta;</mi> <mi>k</mi> <mn>3</mn> </msubsup> <msup> <mi>&amp;eta;</mi> <mn>3</mn> </msup> <msubsup> <mi>&amp;delta;</mi> <mi>k</mi> <mn>4</mn> </msubsup> <mo>-</mo> <mn>3</mn> <msubsup> <mi>&amp;eta;&amp;theta;g</mi> <mi>k</mi> <mn>3</mn> </msubsup> <msubsup> <mi>&amp;zeta;</mi> <mi>k</mi> <mn>2</mn> </msubsup> <msup> <mi>&amp;eta;</mi> <mn>2</mn> </msup> <msubsup> <mi>&amp;delta;</mi> <mi>k</mi> <mn>4</mn> </msubsup> <mo>+</mo> <mn>3</mn> <msup> <mi>&amp;eta;</mi> <mn>2</mn> </msup> <msup> <mi>&amp;theta;</mi> <mn>2</mn> </msup> <msubsup> <mi>g</mi> <mi>k</mi> <mn>2</mn> </msubsup> <msub> <mi>&amp;zeta;</mi> <mi>k</mi> </msub> <msubsup> <mi>&amp;eta;&amp;delta;</mi> <mi>k</mi> <mn>4</mn> </msubsup> <mo>-</mo> <msub> <mi>g</mi> <mi>k</mi> </msub> <msubsup> <mi>&amp;delta;</mi> <mi>k</mi> <mn>4</mn> </msubsup> <msup> <mi>&amp;eta;</mi> <mn>3</mn> </msup> <msup> <mi>&amp;theta;</mi> <mn>3</mn> </msup> <mo>;</mo> </mrow>
求出对应闭式解集合{pk,ρk}，找出满足可行域0≤ρk≤1和max(l1(ρk),l2(ρk))≤pk≤Pk并且使得目标函数取得最大值时对应的即为用户k所求的解，同理，按照此方式余下用户都找到对应解并计算能效值η(n)；
(5)判断是否满足其中ε为判定阈值，如果满足则输出即为K个子问题的最优解，即最优发射功率和功率分裂因子否则执行步骤(3)；
(6)计算传输预编码向量基站利用vk对传输信号进行预编码，同时通过控制信道将每个功率分裂因子发送到相应的用户，从而每个用户设定功率分裂因子，实现信息与能量的同时接收，即完成信能同传系统的收发机设计。
CN201510198147.0A 2015-04-24 2015-04-24 一种低复杂度的多用户信能同传系统能效优化方法 Active CN104954055B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510198147.0A CN104954055B (zh) 2015-04-24 2015-04-24 一种低复杂度的多用户信能同传系统能效优化方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510198147.0A CN104954055B (zh) 2015-04-24 2015-04-24 一种低复杂度的多用户信能同传系统能效优化方法

Publications (2)

Publication Number Publication Date
CN104954055A CN104954055A (zh) 2015-09-30
CN104954055B true CN104954055B (zh) 2018-01-12

Family

ID=54168444

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510198147.0A Active CN104954055B (zh) 2015-04-24 2015-04-24 一种低复杂度的多用户信能同传系统能效优化方法

Country Status (1)

Country Link
CN (1) CN104954055B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106413059B (zh) * 2016-10-31 2019-11-08 北京邮电大学 一种发送功率的确定方法及基站
CN106533526B (zh) * 2016-12-06 2019-11-12 深圳大学 一种受独立功率约束的开关模拟波束成形系统
CN107396380B (zh) * 2017-06-26 2020-05-22 西安交通大学 一种合作认知网络中认知用户能量效率最优化的方法
CN108337024B (zh) * 2018-02-06 2021-02-09 重庆邮电大学 一种基于能量采集的大规模mimo系统能效优化方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103889061A (zh) * 2014-02-18 2014-06-25 南京邮电大学 一种基于多载波通信的d2d用户资源分配方法
CN103997775A (zh) * 2014-06-03 2014-08-20 东南大学 频分复用多用户mimo能效优化方法
CN104066098A (zh) * 2013-03-18 2014-09-24 上海贝尔股份有限公司 用于分布式无线网络的方法和设备

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104066098A (zh) * 2013-03-18 2014-09-24 上海贝尔股份有限公司 用于分布式无线网络的方法和设备
CN103889061A (zh) * 2014-02-18 2014-06-25 南京邮电大学 一种基于多载波通信的d2d用户资源分配方法
CN103997775A (zh) * 2014-06-03 2014-08-20 东南大学 频分复用多用户mimo能效优化方法

Also Published As

Publication number Publication date
CN104954055A (zh) 2015-09-30

Similar Documents

Publication Publication Date Title
CN105827273B (zh) 多小区大规模mimo系统用户双天线导频干扰消除方法
CN104821838B (zh) 基于能效最大化的多用户信能同传系统收发机设计方法
CN103516407B (zh) 传输模式选择方法、天线收发组合确定方法、装置及系统
CN103209051B (zh) 一种协作多点联合传输系统在多用户场景下的两步预编码方法
CN110808765B (zh) 一种基于不完全信道信息的大规模mimo系统频谱效率优化的功率分配方法
CN105634571A (zh) 大规模mimo系统中基于部分导频复用的导频污染减轻方法
CN104954055B (zh) 一种低复杂度的多用户信能同传系统能效优化方法
CN109168197B (zh) 基于异构小蜂窝网络能量采集效率最大化的功率分配方法
CN104702395B (zh) 一种协作认知网络中公平且能效高的联合资源分配方法
CN105450275A (zh) 基于能效最优的多用户大规模天线中继系统天线选择方法
CN101925070A (zh) 一种基于空间复用的认知系统资源分配方法
CN104168574B (zh) 一种混合蜂窝系统中基于干扰适变选择的上行传输方法
CN105871434A (zh) Mimo干扰信道的拓扑干扰对齐方法
CN105450274A (zh) 基于能效最优的大规模多天线中继系统用户数优化方法
CN104617996A (zh) 大规模mimo系统中最大化最小信噪比的预编码设计方法
CN102572864A (zh) 最大化吞吐量的多小区联合波束成形设计方法
CN102355294B (zh) 单基站功率约束的多点协作波束成型和功率分配方法
CN104092519A (zh) 一种基于加权和速率最大化的多用户mimo协作传输方法
CN103905106A (zh) 一种多天线多播波束赋形向量计算方法
CN112054830A (zh) 一种针对硬件损伤的Massive MIMO信能同传系统优化方法
CN103051581B (zh) 一种基于有效容量的mimo-ofdm系统能效优化方法
CN102684765B (zh) 一种基于最大化虚拟信噪比的分布式协作波束赋形方法
CN104320170B (zh) 大规模mimo系统中导频污染抑制波束赋形方法
CN106209188A (zh) 大规模mimo系统中基于部分导频交替复用的导频污染减轻方法
CN104010372B (zh) 大规模mu‑miso系统低复杂度用户调度方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant