CN108770007B - 基于noma的无线携能通信系统多目标优化方法 - Google Patents

基于noma的无线携能通信系统多目标优化方法 Download PDF

Info

Publication number
CN108770007B
CN108770007B CN201810493017.3A CN201810493017A CN108770007B CN 108770007 B CN108770007 B CN 108770007B CN 201810493017 A CN201810493017 A CN 201810493017A CN 108770007 B CN108770007 B CN 108770007B
Authority
CN
China
Prior art keywords
energy
user
power
data rate
optimization
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
CN201810493017.3A
Other languages
English (en)
Other versions
CN108770007A (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.)
South China University of Technology SCUT
Original Assignee
South China University of Technology SCUT
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 South China University of Technology SCUT filed Critical South China University of Technology SCUT
Priority to CN201810493017.3A priority Critical patent/CN108770007B/zh
Publication of CN108770007A publication Critical patent/CN108770007A/zh
Application granted granted Critical
Publication of CN108770007B publication Critical patent/CN108770007B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04WWIRELESS COMMUNICATION NETWORKS
    • H04W24/00Supervisory, monitoring or testing arrangements
    • H04W24/02Arrangements for optimising operational condition
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04WWIRELESS COMMUNICATION NETWORKS
    • H04W72/00Local resource management
    • H04W72/04Wireless resource allocation
    • H04W72/044Wireless resource allocation based on the type of the allocated resource
    • H04W72/0473Wireless resource allocation based on the type of the allocated resource the resource being transmission power
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02DCLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D30/00Reducing energy consumption in communication networks
    • Y02D30/70Reducing energy consumption in communication networks in wireless communication networks

Abstract

本发明公开了一种基于NOMA的无线携能通信系统多目标优化方法,首先建立一个基于NOMA的无线携能通信系统模型并提出优化问题,由于问题为一个复杂的多目标优化问题,利用香农定理将收集到的能量转化为等效数据率,并将优化目标重新定义为系统吞吐量和等效数据率的加权和,将原问题转化为单目标优化问题。由于目标函数对于功率分配和功率分割系数不是联合凹的,优化问题非凸,为此将其拆分成两个子问题并通过对两个子问题迭代求解寻求优化问题的次优解。所述方法在保证满足每个用户对最小数据率和最小收集能量要求的前提下,通过合理设置系统吞吐量和总收集能量的偏好系数,在保证两者公平性的基础上进行同时最大化,从而优化了系统的整体性能。

Description

基于NOMA的无线携能通信系统多目标优化方法
技术领域
本发明涉及无线传输领域,具体涉及一种基于NOMA的无线携能通信系统同时最大化系统吞吐量和总收集能量的多目标优化方法。
背景技术
随着移动互联网、社交网络和物联网的蓬勃发展,网络数据流量和接入设备数量呈爆炸式增长,对无线通信网络提出了更高的服务质量要求,如何提高频谱效率(SpectralEfficiency,SE)和系统能效(Energy Efficiency,EE)为第五代移动通信网络亟需解决的问题。非正交多址接入(Non-orthogonal Multiple Access,NOMA)和无线携能通信(Simultaneous Wireless Information and Power Transfer,SWIPT)技术由于在提高频谱效率和系统能效上具有强大潜力成为5G的关键候选技术,引起广大科研工作者的关注。
NOMA在传输端功率域实现多址接入和多信号分离的基础上,通过采用先进的接收机设计和串行干扰消除(Successive Interference Cancellation,SIC)技术,实现了多个用户共享同一带宽资源,提高了频谱效率。无线携能通信技术利用无线电波在传输信息的同时运输能量使能源有限的系统拥有源源不断的能量来源,保障了无线通信系统的耗电需求,无线通信设备通过接收无线电波的能量为自身充电提高电池续航能力,延长了电池使用寿命,提高了无线通信供电系统的稳定性。两者的结合应用对于无线通信系统性能的改善有望满足5G物联网时代对无线通信系统的关键性能要求,具有广阔的应用前景。
当前已有科研工作者专注于NOMA和SWIPT结合应用的相关研究,主要关注系统吞吐量或用户收集能量的最大化以及对系统通信可靠性的探索。系统吞吐量和总收集能量是两个相互对立的目标变量,不仅如何保证两者之间的公平性值得深入探索,在此基础上对两者进行联合优化更是一个极具实际应用意义的科研课题。
发明内容
本发明的目的是针对5G物联网时代对于无线通信网络高频谱效率和高系统能效的需求,提供一种基于NOMA的无线携能通信系统同时最大化系统吞吐量和总收集能量的多目标优化方法,结合功率分割技术,通过优化功率分配方案和功率分割系数,在保证公平性的基础上,寻求系统能达到的最优性能。所述方法在保证所有用户满足最低数据率和最低收集能量要求的基础上,实现了系统吞吐量和总收集能量的同时最大化,一方面通过多个用户复用频谱提高了频谱利用率,另一方面利用无线携能通信技术对基站发射能量做最大限度的利用,提高了系统能效。
本发明的目的可以通过如下技术方案实现:
一种基于NOMA的无线携能通信系统多目标优化方法,所述方法以同时最大化系统吞吐量和总收集能量为目标,包括以下步骤:
1)建立基于NOMA的无线携能通信系统:假设系统包括一个小区基站BS和N个用户User,所有用户User都配备有信息接收机和能量接收机,具有信息解码(InformationDecoding,ID)和能量收集(Energy Harvesting,EH)能力,另外假设基站BS与每个用户User间的信道状态信息在基站处已知,同时信道增益按升序排列:|h1|2<|h2|2<…|hn|2…<|hN|2,其中|hn|2表示第n个用户的信道增益,在发送端,基站BS对要发送给各个用户User的信号在功率域上进行叠加,pn表示基站BS分配给第n个用户的功率,由于基站的总发射功率P是有限的,因此有:
Figure BDA0001668389880000021
在接收端,用户User采用功率分割技术对从基站BS接收到的能量进行分割,ρnpn用于信息解码,(1-ρn)pn用于能量收集,其中ρn表示功率分割系数;
2)利用功率分割和串行干扰消除技术,根据香农定理得出每个用户通过信息解码获得的数据率以及通过能量收集获得的能量,进而定义系统吞吐量和总收集能量;两个目标变量度量单位不同,不能直接相加,利用香农定理把用户收集到的能量转化为等效数据率以统一两个目标变量的单位,并通过设置合适的偏好系数在保证两目标变量公平性的基础上把优化目标转化为用户通过信息解码获得的系统吞吐量和通过能量收集获得的等效数据率的加权和,将原本的多目标优化问题转化单目标优化问题;
3)由于系统吞吐量和总收集能量关于基站发射功率都是单调递增的,假定基站总以其最大可用总发射功率P发射信号,考虑系统对用户的服务质量要求,将用户需达到的最低数据率Rmin和最小收集能量Emin作为约束条件,提出优化问题;
4)对功率分配和功率分割系数的联合优化使优化问题非凸,难以求解,鉴于穷尽搜索算法的高运算量,提出一个高效的资源分配方案:首先将原两目标变量优化问题拆分成两个单变量子问题,然后通过对两个单变量子问题迭代求解寻求优化问题的次优解。
进一步地,步骤2)中,在接收端,第n个用户通过信息解码获得的数据率为:
Figure BDA0001668389880000031
其中σ2表示信道加性高斯噪声的功率,pi表示基站BS分配给第i个用户的功率,利用香农定理把用户收集到的能量转化为等效数据率,此时第n个用户获得的等效数据率为:
Figure BDA0001668389880000032
其中η表示能量接收效率,
Figure BDA0001668389880000033
表示能量转化为数据率的转化效率,分别用RID、REH表示系统从信息解码和能量收集获得的数据率:
Figure BDA0001668389880000034
将系统通过信息解码获得的数据率的权重归一化设置为1,通过能量收集获得的等效数据率的权重设置为β,此时系统等效总吞吐量定义为:
R=RID+βREH
进一步地,步骤3)中,将用户需达到的最低数据率Rmin和最小收集能量Emin作为约束条件后提出的优化问题如下;
P1:
Figure BDA0001668389880000036
Figure BDA0001668389880000037
Figure BDA0001668389880000038
Figure BDA0001668389880000039
p>0, (1-4)
0<ρ<1, (I-5)
其中约束条件(1-1)、(1-2)分别对用户的最低数据率和最小收集能量做出了限制,约束条件(1-3)表明基站以其最大功率发射信号,p=(p1,p2,…,pN)T、ρ=(ρ1,ρ2,…,ρN)T分别表示各用户的功率分配方案和功率分割系数。
进一步地,步骤4)中将原两目标变量优化问题拆分成如下两个单变量子问题:
P2:
Figure BDA0001668389880000041
Figure BDA0001668389880000042
Figure BDA0001668389880000043
0<ρ<1. (2-3)
P3:
Figure BDA0001668389880000044
Figure BDA0001668389880000045
Figure BDA0001668389880000046
p>0. (3-3)
其中子问题P2假设用户间功率分配方案确定,以功率分割系数为优化变量进行系统优化,P3假设用户功率分割系数确定,以功率分配方案为优化变量进行系统优化,最终原问题通过对以上两个子问题迭代求解取得次优解。
进一步地,对子问题P2的具体求解步骤如下:
a、证明目标函数关于ρ是凹的,P2是凸优化问题,存在最优解;
b、将目标函数重新表示为:
Figure BDA0001668389880000047
其中Rn表示用户n通过信息解码和能量收集获得的等效总数据率,基于各个用户Rn相互独立,将子问题P2拆分成N个具有相同形式的独立的子问题,表示如下:
P4:
Figure BDA0001668389880000048
Figure BDA0001668389880000049
Figure BDA00016683898800000410
0<ρn<1. (4-3)
c、P4为一系列凹最大化问题,根据凹函数的性质,令Rn关于ρn的一阶导数等于0求出最大化Rn的ρn值,求解方程表示如下:
Figure BDA0001668389880000051
Figure BDA0001668389880000052
表示上述方程的唯一正实数根,根据约束条件(4-1)(4-2),子问题P4的最优解最终给出如下:
Figure BDA0001668389880000053
其中
Figure BDA00016683898800000514
确保流向信息接收机的能量使用户达到最低数据率要求,
Figure BDA00016683898800000515
确保流向能量接收机的能量使用户达到最小收集能量要求。
进一步地,对子问题P3的具体求解步骤如下:
a、证明目标函数在当且仅当ρ1|h1|2≤ρ2|h2|2≤…≤ρN|hN|2的前提条件下关于p是凹的,P3为凸优化问题,存在最优解;
b、利用Lagrange函数和KKT最优性条件对子问题P3进行求解,解得最大化系统等效总吞吐量的最优功率分配方案为:
Figure BDA0001668389880000054
进一步地,子问题P2的最优解可能出现的三种情况划分如下:
(1)、所有
Figure BDA0001668389880000055
都落在区间
Figure BDA0001668389880000056
内;
(2)、至少存在一个
Figure BDA0001668389880000057
Figure BDA0001668389880000058
(3)、至少存在一个
Figure BDA0001668389880000059
Figure BDA00016683898800000510
为确保基于优化的功率分割系数子问题P3是凸的,证明前两种情况均能确保
Figure BDA00016683898800000511
对于第三种情况,由于
Figure BDA00016683898800000512
受多个参数影响,难以确定各用户最小功率分割系数间的关系,将求得的各个用户的最优功率分割系数与相对应的信道增益相乘并将乘积结果按升序排列,在确保迭代求解过程中两个子问题始终为凸优化问题的前提下,对两个子问题进行迭代求解直至目标函数收敛;具体算法如下:
I.随机初始化用户的功率分割系数为ρ(1)
Figure BDA00016683898800000513
Figure BDA0001668389880000061
δ=10-5为精确度;
II.基于ρ(1)对子问题P3进行求解,记求得的最优功率分配方案为p(1),获得的系统等效总吞吐量为
Figure BDA0001668389880000062
III.初始化收敛标志为C=0,迭代计数参数t=2,重复以下过程直至C=1:
a、基于p(t-1)对子问题P2进行求解,记求得的最优功率分割系数为ρ(t),系统等效总吞吐量为
Figure BDA0001668389880000063
判断是否存在
Figure BDA0001668389880000064
若存在,则将功率分割系数与相对应的信道增益相乘并按升序排列;
b、判断
Figure BDA0001668389880000065
是否成立,若不成立,令C=1,得出可达的最大系统等效总吞吐量为
Figure BDA0001668389880000066
相对应的功率分配方案和功率分割系数分别为p*=p(t-1),ρ(*)=ρ(t-1)
c、若
Figure BDA0001668389880000067
成立,基于ρ(t)对子问题P3进行求解,记求得的最优功率分配方案为p(t),获得的系统等效总吞吐量为
Figure BDA0001668389880000068
d、判断
Figure BDA0001668389880000069
是否成立,若不成立,令C=1,得出可达的最大系统等效总吞吐量为
Figure BDA00016683898800000610
相对应的的功率分配方案和功率分割系数分别为p*=p(t-1),ρ(*)=ρ(t)
e、若
Figure BDA00016683898800000611
成立,t=t+1。
本发明与现有技术相比,具有如下优点和有益效果:
本发明提供的基于NOMA的无线携能通信系统同时最大化系统吞吐量和总收集能量的多目标优化方法,利用NOMA和SWIPT技术的优势,结合功率分割技术,不局限于单方面地最大化系统吞吐量或总收集能量,而是通过合理设置系统吞吐量和总收集能量的偏好系数来保证两者公平性,在保证所有用户满足最低数据率和最低收集能量要求的基础上,寻找最优的功率分配方案和功率分割控制方案同时最大化系统吞吐量和总收集能量。一方面通过利用NOMA,实现多个用户复用同一频谱资源,提高了频谱利用率;另一方面利用无线携能通信技术实现信息和能量的并行传输,提高了系统能效。
附图说明
图1为本发明实施例各用户功率分割系数ρn和等效总吞吐量R的关系示意图。
图2为本发明实施例系统用户数量不同情况下最低数据率Rmin和等效总吞吐量R的关系示意图。
图3为本发明实施例最低数据率Rmin和系统吞吐量RID、总收集能量E的关系示意图。
图4为本发明实施例基站发射功率P和等效总吞吐量R的关系示意图。
图5为本发明实施例偏好系数β和系统吞吐量RID、总收集能量E的关系示意图。
具体实施方式
下面结合实施例及附图对本发明作进一步详细的描述,但本发明的实施方式不限于此。
实施例:
本实施例提供了一种基于NOMA的无线携能通信系统同时最大化系统吞吐量和总收集能量的多目标优化方法,所述方法包括以下步骤:
步骤一、假设我们所讨论的基于NOMA的无线携能通信系统有1个单天线基站BS部署在小区中间,6个单天线用户User随机部署在基站覆盖范围内,所有用户均配备有信息接收机和能量接收机。基站对要发送给用户的信号在功率域上进行叠加,用户采用功率分割技术对从下行链路中接收到的信号进行功率分割,用ρn表示用户n用于ID的功率分割系数,即流向信息接收机的功率为ρnpn,流向能量接收机的功率为(1-ρn)pn,其中pn表示基站分配给用户n的功率。本实施例采用高斯信道模型,用|hn|2表示基站与用户间的信道增益,信道增益服从均值为0,方差为1的高斯分布,并按升序排列|h1|2<|h2|2<…<|hN|2
步骤二、在接收端用户通过SIC进行信息解码获得的数据率为:
Figure BDA0001668389880000071
其中σ2=0.001W为信道噪声功率。用户收集到的能量为:
Figure BDA0001668389880000072
其中η=0.1为功率接收效率,最终得出系统的吞吐量和总收集能量表示如下:
Figure BDA0001668389880000073
Figure BDA0001668389880000074
吞吐量的单位为bit/s/Hz,总收集能量的单位是Watt,两者度量单位不同,不能直接相加,利用香农定理把用户收集到的能量转化为等效数据率:
Figure BDA0001668389880000075
其中
Figure BDA0001668389880000076
表示能量转化为数据率的转化效率。系统通过EH获得的等效总数据率:
Figure BDA00016683898800000811
此时RID和REH的单位一致,将前者的权重归一化设置为1,后者的权重设置为β,系统等效总吞吐量定义为:
R=RID+βREH
不失一般性,在求解优化问题过程中我们假设偏好系数β为一固定的实数值,这里假设β=0.1。
步骤三、假定基站以其最大可用功率P=40W发射信号,用户需达到的最低数据率Rmin=2bit/s/Hz,最小收集能量Emin=0.01W,优化问题最终定义如下:
P1:
Figure BDA0001668389880000081
Figure BDA0001668389880000082
Figure BDA0001668389880000083
Figure BDA0001668389880000084
p>0, (1-4)
0<ρ<1, (1-5)
步骤四、将优化问题具体拆分如下为两个单变量子问题:
P2:
Figure BDA0001668389880000085
Figure BDA0001668389880000086
Figure BDA0001668389880000087
0<ρ<1. (2-3)
P3:
Figure BDA0001668389880000088
Figure BDA0001668389880000089
Figure BDA00016683898800000810
p>0. (3-3)
其中子问题P2基于固定的功率分配方案寻求最优的功率分割系数,P3基于固定的功率分割系数寻求最优的功率分配方案。
1)对子问题P2进行求解:
根据各用户间的功率分割系数相互独立可得:
Figure BDA0001668389880000091
据此将目标函数关于ρ的海森矩阵简化为:
Figure BDA0001668389880000092
其中:
Figure BDA0001668389880000093
Figure BDA0001668389880000099
因此
Figure BDA00016683898800000910
为负定矩阵,目标函数关于ρ是凹的。同时由于约束条件(2-1)是凸的,(2-2)(2-3)是仿射的,判定P2为一个凸优化问题,存在最优解。
将目标函数重新表示为:
Figure BDA0001668389880000094
其中Rn表示用户n通过ID和EH获得的个体等效总数据率。基于各个用户Rn相互独立,将P2拆分成N个具有相同形式的独立的子问题:
P4:
Figure BDA0001668389880000095
Figure BDA0001668389880000096
Figure BDA0001668389880000097
0<ρn<1. (4-3)
P4为一系列凹最大化问题,目标函数Rn关于ρn是严格凹的,根据凹函数的性质,令Rn关于ρn的一阶导数等于0求得最大化Rn的ρn值:
Figure BDA0001668389880000098
Figure BDA00016683898800000911
表示上述方程的唯一正实数根,由于最优功率分割系数需满足最小数据率和最小收集能量要求,根据约束条件(4-1)(4-2),子问题P4的最优解最终给出如下:
Figure BDA0001668389880000101
其中
Figure BDA0001668389880000106
2)对子问题P3进行求解:
证明子问题P3目标函数在当且仅当ρ1|h1|2≤ρ2|h2|2≤…≤ρN|hN|2的前提条件下关于p是凹的:
目标函数R关于功率分配因子的二阶导数表示如下:
Figure BDA0001668389880000102
由以上式子可得:
Figure BDA0001668389880000103
用wn表示
Figure BDA0001668389880000107
H表示目标函数关于p的海森矩阵,H的n阶顺序主子式为:
Figure BDA0001668389880000104
由w1<0可得H的一阶顺序主子式是负的,当n=2,3,…,N时,wi-wi-1表示如下:
Figure BDA0001668389880000105
当且仅当ρ1|h1|2≤ρ2|h2|2≤…≤ρN|hN|2时H的奇数阶顺序主子式非正,偶数阶顺序主子式非负,此时H为半负定矩阵,因此目标函数R关于p是凹的,P3为凸优化问题。
假设功率分割系数满足ρ1|h1|2<ρ2|h2|2<…<ρN|hN|2,P3的Lagrange函数为:
Figure BDA0001668389880000111
其中λ为对应约束条件(3-1)的N阶Lagrange乘子向量,μ为对应约束条件(3-2)的Lagrange乘子,fn以及h定义如下:
Figure BDA0001668389880000112
Figure BDA0001668389880000113
KKT最优性条件为:
Figure BDA0001668389880000114
Figure BDA0001668389880000115
h=0,
Figure BDA0001668389880000116
μh=0,
Figure BDA0001668389880000117
μ≥0,
Figure BDA0001668389880000118
表示任意一个用户,可得
Figure BDA0001668389880000119
因此有
Figure BDA00016683898800001110
Figure BDA00016683898800001111
可得
Figure BDA00016683898800001112
由于
Figure BDA0001668389880000124
当n=1,2,…,N-1时λn>0,又因为λnfn=0,因此
Figure BDA0001668389880000122
这意味着约束条件(3-1)对n=1,2,…,N-1取等号。当λN>0时,
Figure BDA0001668389880000123
计算出基站此时的发射功率,记为Pmin,Pmin即为所有用户均能达到最低数据率要求时基站所需的最小发射功率,为保证我们的讨论有实际意义,我们假设P≥Pmin
利用KKT最优性条件获得子问题P3的最优解给出如下:
当P=Pmin时,对于所有的用户,令
Figure BDA0001668389880000125
此时分配给各个用户的功率可按n=N,N-1,…,1的顺序得出。
当P>Pmin时,对于信道条件较差的前N-1个的用户,令
Figure BDA0001668389880000126
得出分配给前N-1个的用户的功率,剩余的功率全部分配给信道条件最好的用户。
即最大化系统等效总吞吐量的最优功率分配方案为:
Figure BDA0001668389880000121
3)对子问题P2的求解结果进行讨论,最优解可能出现的三种情况划分如下:
(1)所有
Figure BDA0001668389880000127
都落在区间
Figure BDA0001668389880000128
内;
(2)至少存在一个
Figure BDA0001668389880000129
Figure BDA00016683898800001210
(3)至少存在一个
Figure BDA00016683898800001211
Figure BDA00016683898800001212
基于最优的功率分配方案,
Figure BDA00016683898800001213
通过令Rn关于ρn的一阶导数等于0求得,并且求解结果有
Figure BDA00016683898800001214
对于情况(1),所有
Figure BDA00016683898800001215
都落在区间
Figure BDA00016683898800001216
内,对所有用户有
Figure BDA00016683898800001217
因而有
Figure BDA00016683898800001218
成立。对于情况(2),不失一般性,假设存在用户n,
Figure BDA00016683898800001219
此时
Figure BDA00016683898800001221
Figure BDA00016683898800001220
时,
Figure BDA00016683898800001222
考虑另一种情况,
Figure BDA00016683898800001223
此时
Figure BDA00016683898800001224
因此对于第二种情景,依然有
Figure BDA00016683898800001225
成立。
对于情况(3),由于
Figure BDA00016683898800001226
受多个参数影响,难以确定各用户最小功率分割系数间的关系,无法确保
Figure BDA00016683898800001227
成立。为此,当存在
Figure BDA00016683898800001228
时,将求得的用户最优功率分割系数与相对应的信道增益相乘并将乘积结果按升序排列。在确保迭代求解过程中两个子问题始终为凸优化问题的前提下,对两个子问题进行迭代求解直至目标函数收敛,具体算法如下:
I.随机初始化用户的功率分割系数为ρ(1)
Figure BDA00016683898800001229
Figure BDA00016683898800001230
δ=10-5为精确度;
II.基于ρ(1)对子问题P3进行求解,记求得的最优功率分配方案为p(1),获得的系统等效总吞吐量为
Figure BDA0001668389880000131
III.初始化收敛标志为C=0,迭代计数参数t=2,重复以下过程直至C=1:
a、基于p(t-1)对子问题P2进行求解,记求得的最优功率分割系数为ρ(t),系统等效总吞吐量为
Figure BDA0001668389880000132
判断是否存在
Figure BDA0001668389880000133
若存在,则将功率分割系数与相对应的信道增益相乘并按升序排列;
b、判断
Figure BDA0001668389880000134
是否成立,若不成立,令C=1,得出可达的最大系统等效总吞吐量为
Figure BDA0001668389880000135
相对应的的功率分配方案和功率分割系数分别为p*=p(t-1),ρ(*)=ρ(t-1)
c、若
Figure BDA0001668389880000136
成立,基于ρ(t)对子问题P3进行求解,记求得的最优功率分配方案为p(t),获得的系统等效总吞吐量为
Figure BDA0001668389880000137
d、判断
Figure BDA0001668389880000138
是否成立,若不成立,令C=1,得出可达的最大系统等效总吞吐量为
Figure BDA0001668389880000139
相对应的的功率分配方案和功率分割系数分别为p*=p(t-1),ρ(*)=ρ(t)
e、若
Figure BDA00016683898800001310
成立,t=t+1。
图1-图5为本实施例提供的基于NOMA的无线携能通信系统同时最大化系统吞吐量和总收集能量的多目标优化方案的仿真效果图。
图1证明了基于优化的功率分配方案,各用户等效数据率Rn关于功率分割系数ρn是凹的,对于每一个用户,存在一个最优功率分割系数使等效数据率Rn取得最大值,并且这一最大值随着信道条件的改善而增大。
图2的其他参数:用户数量分别为N=4,5,6。从图2可以看出,对于固定用户数量的系统,系统的等效总吞吐量随最低数据率要求的提高而降低;另一方面,随着用户数量的增加,系统可以获得的最大等效总吞吐量增大,增大的速度随着用户数量的增大而降低。这意味着尽管用户数量的增加有利于提高系统等效总吞吐量,但是由于基站发射功率有限,用户数量的增加会加大整个系统的负担,限制系统性能的提升。
图3对图2中系统的等效总吞吐量随最低数据率要求提高而降低的现象进行解释。从图3可以看出,随着最低数据率要求的提高,系统吞吐量和总收集能量均呈下降趋势。一方面更多的能量流向信息接收机以提高用户的数据率使其达到服务质量要求,这导致了更少的能量流向能量接收机,系统收集到的能量减少;另一方面,为达到提高了的数据率要求,信道条件较差的前N-1个用户需要更多的功率,这导致信道条件最好的用户分配到的功率减小,数据率减少,前N-1个用户数据率的提升不足以弥补信道条件最好的用户数据率的降低,最终导致系统整体吞吐量降低。系统吞吐量和总收集能量的减小最终导致等效总吞吐量的降低。
图4加入基于正交多址接入(OMA,orthogonal multiple access)的无线携能通信系统做对比,在基于OMA的无线携能通信系统中,基站发射功率平均分配给每一个用户,用户功率分割系数取最大化用户个体等效总数据率的最优值。从图4可以看出,随着基站发射功率的增大,两个系统的等效总吞吐量均增大,同时,我们所提系统的等效总吞吐量明显高于基于OMA的系统,这表明NOMA和SWIPT的结合应用能显著提升系统性能。
从图5可以看出,随着偏好系数β的增大系统吞吐量减小,总收集能量增大。偏好系数β的增大表明我们更看重用户收集能量的能力,更多的能量流向用户能量接收机,导致流向信息接收机的能量减少,最终使系统吞吐量减小,总收集能量增大。不存在绝对意义上最优的偏好系数β使系统性能达到最优,β的取值有赖于决策者对系统性能的具体要求,因此根据具体需求合理确定β的取值具有重要意义。
以上所述,仅为本发明专利较佳的实施例,但本发明专利的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明专利所公开的范围内,根据本发明专利的技术方案及其发明专利构思加以等同替换或改变,都属于本发明专利的保护范围。

Claims (7)

1.一种基于NOMA的无线携能通信系统多目标优化方法,其特征在于,所述方法以同时最大化系统吞吐量和总收集能量为目标,包括以下步骤:
1)建立基于NOMA的无线携能通信系统:假设系统包括一个小区基站和N个用户,所有用户都配备有信息接收机和能量接收机,具有信息解码和能量收集能力,另外假设基站与每个用户间的信道状态信息在基站处已知,同时信道增益按升序排列:|h1|2<|h2|2<…|hn|2…<|hN|2,其中|hn|2表示第n个用户的信道增益,在发送端,基站对要发送给各个用户的信号在功率域上进行叠加,pn表示基站分配给第n个用户的功率,由于基站的总发射功率P是有限的,因此有:
Figure FDA0002390636110000011
在接收端,用户采用功率分割技术对从基站接收到的能量进行分割,ρnpn用于信息解码,(1-ρn)pn用于能量收集,其中ρn表示功率分割系数;
2)利用功率分割和串行干扰消除技术,根据香农定理得出每个用户通过信息解码获得的数据率以及通过能量收集获得的能量,进而计算得出系统吞吐量和总收集能量;两个目标变量度量单位不同,不能直接相加,利用香农定理把用户收集到的能量转化为等效数据率以统一两个目标变量的单位,并通过设置合适的偏好系数在保证两目标变量公平性的基础上把优化目标转化为用户通过信息解码获得的系统吞吐量和通过能量收集获得的等效数据率的加权和,将原本的多目标优化问题转化单目标优化问题;
3)由于系统吞吐量和总收集能量关于基站发射功率都是单调递增的,假定基站总以其最大可用总发射功率P发射信号,考虑系统对用户的服务质量要求,将用户需达到的最低数据率Rmin和最小收集能量Emin作为约束条件,提出优化问题;
4)对功率分配和功率分割系数的联合优化使优化问题非凸,难以求解,鉴于穷尽搜索算法的高运算量,提出一个高效的资源分配方案:首先将原两目标变量优化问题拆分成两个单变量子问题,然后通过对两个单变量子问题迭代求解寻求优化问题的次优解。
2.根据权利要求1所述的一种基于NOMA的无线携能通信系统多目标优化方法,其特征在于,步骤2)中,在接收端,第n个用户通过信息解码获得的数据率为:
Figure FDA0002390636110000012
其中σ2表示信道加性高斯噪声的功率,pi表示基站BS分配给第i个用户的功率,利用香农定理把用户收集到的能量转化为等效数据率,此时第n个用户获得的等效数据率为:
Figure FDA0002390636110000021
其中η表示能量接收效率,
Figure FDA0002390636110000028
表示能量转化为数据率的转化效率,分别用RID、REH表示系统从信息解码和能量收集获得的数据率:
Figure FDA0002390636110000022
Figure FDA0002390636110000023
将系统通过信息解码获得的数据率的权重归一化设置为1,通过能量收集获得的等效数据率的权重设置为β,此时系统等效总吞吐量定义为:
R=RID+βREH
3.根据权利要求2所述的一种基于NOMA的无线携能通信系统多目标优化方法,其特征在于,步骤3)中,将用户需达到的最低数据率Rmin和最小收集能量Emin作为约束条件后提出的优化问题如下;
P1:
Figure FDA0002390636110000024
Figure FDA0002390636110000025
Figure FDA0002390636110000026
Figure FDA0002390636110000027
p>0, (1-4)
0<ρ<1, (1-5)
其中约束条件(1-1)、(1-2)分别对用户的最低数据率和最小收集能量做出了限制,约束条件(1-3)表明基站以其最大功率发射信号,p=(p1,p2,…,pN)T、ρ=(ρ1,ρ2,…,ρN)T分别表示各用户的功率分配和功率分割系数。
4.根据权利要求3所述的一种基于NOMA的无线携能通信系统多目标优化方法,其特征在于,步骤4)中将原两目标变量优化问题拆分成如下两个单变量子问题:
P2:
Figure FDA0002390636110000031
Figure FDA0002390636110000032
Figure FDA0002390636110000033
0<ρ<1. (2-3)
P3:
Figure FDA0002390636110000034
Figure FDA0002390636110000035
Figure FDA0002390636110000036
p>0. (3-3)
其中子问题P2假设用户间功率分配方案确定,以功率分割系数为优化变量进行系统优化,P3假设用户功率分割系数确定,以功率分配方案为优化变量进行系统优化,最终原问题通过对以上两个子问题迭代求解取得次优解。
5.根据权利要求4所述的一种基于NOMA的无线携能通信系统多目标优化方法,其特征在于,对子问题P2的具体求解步骤如下:
a、证明目标函数关于ρ是凹的,P2是凸优化问题,存在最优解;
b、将目标函数重新表示为:
Figure FDA0002390636110000037
其中Rn表示用户n通过信息解码和能量收集获得的个体等效总数据率,基于各个用户Rn相互独立,将子问题P2拆分成N个具有相同形式的独立的子问题,表示如下:
P4:
Figure FDA0002390636110000038
Figure FDA0002390636110000039
Figure FDA00023906361100000310
0<ρn<1. (4-3)
c、P4为一系列凹最大化问题,根据凹函数的性质,令Rn关于ρn的一阶导数等于0求出最大化Rn的ρn值,求解方程表示如下:
Figure FDA0002390636110000041
Figure FDA0002390636110000042
表示上述方程的唯一正实数根,根据约束条件(4-1)(4-2),子问题P4的最优解最终给出如下:
Figure FDA0002390636110000043
其中
Figure FDA00023906361100000413
确保流向信息接收机的能量使用户达到最低数据率要求,
Figure FDA00023906361100000414
确保流向能量接收机的能量使用户达到最小收集能量要求。
6.根据权利要求5所述的一种基于NOMA的无线携能通信系统多目标优化方法,其特征在于,对子问题P3的具体求解步骤如下:
a、证明目标函数在当且仅当ρ1|h1|2≤ρ2|h2|2≤…≤ρN|hN|2的前提条件下关于p是凹的,P3为凸优化问题,存在最优解;
b、利用Lagrange函数和KKT最优性条件对子问题P3进行求解,解得最大化系统等效总吞吐量的最优功率分配方案为:
Figure FDA0002390636110000044
7.根据权利要求6所述的一种基于NOMA的无线携能通信系统多目标优化方法,其特征在于,子问题P2的最优解可能出现的三种情况划分如下:
(1)、所有
Figure FDA0002390636110000045
都落在区间
Figure FDA0002390636110000046
内;
(2)、至少存在一个
Figure FDA0002390636110000047
Figure FDA0002390636110000048
(3)、至少存在一个
Figure FDA0002390636110000049
Figure FDA00023906361100000410
为确保基于优化的功率分割系数子问题P3是凸的,证明前两种情况均能确保
Figure FDA00023906361100000411
对于第三种情况,由于
Figure FDA00023906361100000412
受多个参数影响,难以确定各用户最小功率分割系数间的关系,将求得的各个用户的最优功率分割系数与相对应的信道增益相乘并将乘积结果按升序排列,在确保迭代求解过程中两个子问题始终为凸优化问题的前提下,对两个子问题进行迭代求解直至目标函数收敛;具体算法如下:
I.随机初始化用户的功率分割系数为ρ(1)
Figure FDA0002390636110000051
Figure FDA0002390636110000052
δ=10-5为精确度;
II.基于ρ(1)对子问题P3进行求解,记求得的最优功率分配方案为p(1),获得的系统等效总吞吐量为
Figure FDA0002390636110000053
III.初始化收敛标志为C=0,迭代计数参数t=2,重复以下过程直至C=1:
a、基于p(t-1)对子问题P2进行求解,记求得的最优功率分割系数为ρ(t),系统等效总吞吐量为
Figure FDA0002390636110000054
判断是否存在
Figure FDA0002390636110000055
若存在,则将功率分割系数与相对应的信道增益相乘并按升序排列;
b、判断
Figure FDA0002390636110000056
是否成立,若不成立,令C=1,得出可达的最大系统等效总吞吐量为
Figure FDA0002390636110000057
相对应的功率分配方案和功率分割系数分别为p*=p(t-1),ρ(*)=ρ(t-1)
c、若
Figure FDA0002390636110000058
成立,基于ρ(t)对子问题P3进行求解,记求得的最优功率分配方案为p(t),获得的系统等效总吞吐量为
Figure FDA00023906361100000512
d、判断
Figure FDA0002390636110000059
是否成立,若不成立,令C=1,得出可达的最大系统等效总吞吐量为
Figure FDA00023906361100000510
相对应的的功率分配方案和功率分割系数分别为p*=p(t-1),ρ(*)=ρ(t)
e、若
Figure FDA00023906361100000511
成立,t=t+1。
CN201810493017.3A 2018-05-22 2018-05-22 基于noma的无线携能通信系统多目标优化方法 Active CN108770007B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810493017.3A CN108770007B (zh) 2018-05-22 2018-05-22 基于noma的无线携能通信系统多目标优化方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810493017.3A CN108770007B (zh) 2018-05-22 2018-05-22 基于noma的无线携能通信系统多目标优化方法

Publications (2)

Publication Number Publication Date
CN108770007A CN108770007A (zh) 2018-11-06
CN108770007B true CN108770007B (zh) 2020-07-28

Family

ID=64007598

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810493017.3A Active CN108770007B (zh) 2018-05-22 2018-05-22 基于noma的无线携能通信系统多目标优化方法

Country Status (1)

Country Link
CN (1) CN108770007B (zh)

Families Citing this family (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110191476B (zh) * 2019-04-18 2020-07-14 浙江大学 一种基于可重构天线阵列的非正交多址接入方法
CN110337148B (zh) * 2019-04-26 2022-04-05 重庆邮电大学 基于非正交多址接入的认知无线电能效资源分配方法
CN110278576B (zh) * 2019-04-26 2021-11-16 南京邮电大学 一种无线能量采集非正交多址接入系统资源分配方法
CN110166107A (zh) * 2019-05-17 2019-08-23 武汉大学 基于无线携能通信网络的无人机中继系统资源分配方法
CN110381589B (zh) * 2019-06-12 2023-08-08 广西智能交通科技有限公司 基于无线能量收集的协作无线网络子信道功率分配方法
CN110602722B (zh) * 2019-09-17 2021-07-02 北京邮电大学 一种基于noma的联合内容推送和传输的设计方法
WO2021068808A1 (en) * 2019-10-09 2021-04-15 Tcl Communication Limited Non-orthogonal multiple access method, base station, and user equipment
CN110730025B (zh) * 2019-10-22 2022-03-25 西北师范大学 一种适用于携能非正交多址通信系统的增量中继方法
CN110944375B (zh) * 2019-11-22 2021-01-12 北京交通大学 无线信息与能量同传辅助的雾计算网络资源配置的方法
CN110933757B (zh) * 2019-11-26 2022-04-05 重庆邮电大学 基于时间反演的wpcn系统的抗干扰资源分配方法
CN110868756B (zh) * 2019-11-26 2022-05-03 重庆邮电大学 一种noma系统中的功率分配方法及装置
CN111314959A (zh) * 2020-02-19 2020-06-19 西南交通大学 一种面向缓存辅助非正交多址网络的动态资源分配方法
CN111556577B (zh) * 2020-05-20 2022-06-07 西安电子科技大学 Noma系统中基于引力搜索算法的吞吐量最大化方法
CN111918320B (zh) * 2020-07-14 2022-05-24 吉林大学 时分双工下非正交多址接入的无线通信资源优化分配方法
CN113038510B (zh) * 2021-03-12 2022-11-15 上海海事大学 基于agm不等式变换的swipt与noma系统能量效率最优化方法
CN113114322B (zh) * 2021-04-19 2022-07-26 华南理工大学 无线携能通信系统中基于moead的波束赋形方法
CN113347728B (zh) * 2021-06-23 2022-12-06 西安邮电大学 反向散射网络吞吐量优化方法、装置及存储介质
CN113965245A (zh) * 2021-09-30 2022-01-21 广西电网有限责任公司柳州供电局 基于opgw接头盒的智能反射面通信系统资源优化方法
CN114520989B (zh) * 2022-01-21 2023-05-26 重庆邮电大学 多载波数能同传noma网络能效最大化方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101428919B1 (ko) * 2013-04-09 2014-08-11 한국과학기술원 다중 무선 접속 시스템에서의 상향 링크 다중 접속 방법 및 장치
CN106304362A (zh) * 2016-08-14 2017-01-04 辛建芳 一种基于ofdm的中继系统能效优化方法
CN107396380A (zh) * 2017-06-26 2017-11-24 西安交通大学 一种合作认知网络中认知用户能量效率最优化的方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101428919B1 (ko) * 2013-04-09 2014-08-11 한국과학기술원 다중 무선 접속 시스템에서의 상향 링크 다중 접속 방법 및 장치
CN106304362A (zh) * 2016-08-14 2017-01-04 辛建芳 一种基于ofdm的中继系统能效优化方法
CN107396380A (zh) * 2017-06-26 2017-11-24 西安交通大学 一种合作认知网络中认知用户能量效率最优化的方法

Also Published As

Publication number Publication date
CN108770007A (zh) 2018-11-06

Similar Documents

Publication Publication Date Title
CN108770007B (zh) 基于noma的无线携能通信系统多目标优化方法
CN109005551B (zh) 一种非理想信道状态信息的多用户noma下行功率分配方法
CN108495337B (zh) 基于noma的无线携能通信系统最大安全速率优化方法
CN108601087B (zh) 一种基于网络切片的无线通信资源分配方法
CN111314894B (zh) 一种面向noma与携能d2d融合网络的鲁棒资源分配方法
CN105680920B (zh) 一种多用户多天线数能一体化通信网络吞吐量优化方法
CN110519848B (zh) 一种认知中继无线传感器网络的联合资源分配方法
CN101951307A (zh) 在CoMP下的小区协作集的选择方法
CN111918320B (zh) 时分双工下非正交多址接入的无线通信资源优化分配方法
CN110933757B (zh) 基于时间反演的wpcn系统的抗干扰资源分配方法
CN113194492B (zh) 一种基于α公平性的安全D2D通信资源分配方法
CN107241180B (zh) 一种支持信息与能量同传的高效资源分配方法
CN109788540A (zh) D2d系统中基于能量采集的功率控制和信道分配方法
CN109451569A (zh) 一种无线携能异构网络中的资源分配方法
CN103944618B (zh) 大规模miso协同能效发送方法
CN105451315A (zh) 吞吐量最大化的串行能量采集方法
CN112702792B (zh) 基于gfdm的无线携能网络上下行链路资源联合分配方法
CN110418360A (zh) 无线携能网络多用户子载波比特联合分配方法
CN112105077B (zh) 一种基于swipt技术的大规模mimo系统uav中继通信方法
CN113507716A (zh) 一种基于swipt的cr-noma网络中断与能效的优化方法
CN107613565B (zh) 一种全双工超密集网络中的无线资源管理方法
Long et al. Throughput optimization in energy harvesting based cognitive iot with cooperative sensing
CN107733488B (zh) 一种大规模mimo系统中注水功率分配改进方法及系统
CN114520989B (zh) 多载波数能同传noma网络能效最大化方法
CN110381590B (zh) 一种基于无线能量收集的无线网络子信道功率分配方法

Legal Events

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