一种RIS辅助无线供电通信网络的优化方法

一种ris辅助无线供电通信网络的优化方法
技术领域
1.本发明涉及无线供电系统优化技术领域，更具体地，涉及一种ris辅助无线供电通信网络的优化方法。


背景技术：

2.到2030年，全球物联网设备将达到5000亿台。能源受限已经成为物联网发展的一个主要问题，能源收集技术的最新进展，特别是无线射频技术，为可持续物联网设备从专用的环境射频源收集能源开辟了一条新途径，这促使了无线供电通信网络(wireless powered communication networks,wpcns)的出现。在这种网络中，物联网设备可以从专用的混合接入点收集能量，然后使用收集的能量将数据传输到接入点。但是在广泛应用于实际应用之前，wpcns还面临着一些挑战。特别是在大规模工业物联网中，设备的部署是随机的，无法控制它们所在的位置，射频能量信号会由于距离遭受衰减，并且出现障碍物时，影响通信质量。因此，需要更高效、更经济的解决方案来提高wpcns的能量效率和通信质量，以保证wpcns能够令人满意的性能适配到物联网环境中。
3.公开号cn112272384a，公开日为2021-01-26的发明专利一种基于可重构智能表面的通信系统吞吐量优化方法，应用于具有可重构智能反射面的通信系统，并对通信系统进行优化，但该专利技术方案未考虑ris元件功耗，考虑的约束条件不够全面，因此与实际情况不够相符。


技术实现要素：

4.本发明为克服上述技术问题，提供对能量传输距离较长，与实际情况更符合的一种ris辅助无线供电通信网络的优化方法。
5.本发明技术方案如下：
6.一种ris辅助无线供电通信网络的优化方法，包括以下步骤：
7.s1：建立ris辅助无线供电通信网络系统模型；
8.s2：确定所述ris辅助无线供电通信网络系统模型优化的原问题，所述原问题是指，在ris辅助无线供电通信网络系统模型的系统参数满足预设的初始约束条件的情况下，使系统吞吐量最大；
9.s3：将所述原问题简化为若干个子问题，所述子问题需满足对应的子约束条件，若所述子问题的子约束条件存在非凸约束，则该子问题为非凸子问题，将所述非凸子问题分类为可解非凸子问题和不可解非凸子问题；
10.s4：对不可解非凸子问题进行凸优化处理，使其子约束条件的非凸约束全部转化为凸约束，从而将所有不可解非凸子问题转化为可解的凸子问题；
11.s5：建立迭代算法，通过所述迭代算法对所述子问题进行交替迭代求解，得到使ris辅助无线供电通信网络系统模型的系统吞吐量最大时，所述ris辅助无线供电通信网络系统模型模型的各系统参数。
12.本技术方案提出了一种ris辅助无线供电通信网络的优化方法，首先构建包括可重构智能反射面的ris辅助无线供电通信网络系统模型，其中，可重构智能反射面集成有大量低成本的无源反射元件，能自适应调整反射信号的相位，从而获得更高的反射波束赋形增益，实现更长的无线能量传输距离。并且通过更全面的考虑系统优化的约束条件，得到较好的系统优化效果，获得最优的系统吞吐量。
13.进一步地，步骤s1所述ris辅助无线供电通信网络系统模型包括：具有n个反射元件的ris、天线数为t的多天线基站、多个单天线用户；
14.其中ris辅助无线供电通信网络系统模型中一共有个用户，被分为k簇，每一簇有mk个用户，u
k,m
表示第k簇的第m个用户，m∈ω＝{1,...,mk}；系统的工作被分成两个部分，包括下行传输能量和上行传输信息，基站通过下行链路给用户传输能量，用户通过上行链路向基站传输信息。
15.进一步地，步骤s2所述系统参数包括：w,τ0,{τk},φ0,{φk},βe，其中w表示基站的波束赋形矢量；τ0为基站通过下行链路传输能量的下行时隙；{τk}为用户簇中的用户终端通过上行链路向基站上行传输的上行时隙集合；φ0为是下行时隙τ0的ris反射系数矩阵；{φk}为上行时隙集合{τk}中的上行时隙τk的ris反射系数矩阵；βe为ris的反射系数。
16.进一步地，步骤s2所述原问题的目标函数p1为所述原问题的目标函数用于在满足所述预设的初始约束条件的情况下，计算使网络通信系统吞吐量最大的各系统参数的取值；目标函数的推导过程为：
17.所述ris辅助无线供电通信网络系统模型为非线性能量收集模型，因此ris和用户u
k,m
所获得的能量分别表示为：
[0018][0019]ek，m
＝min(η|wh(ghβeφ0g
i，(k，m)
g
d，（k，m）
|2，p
(k，m)，sat)
t0(2)
[0020]
其中η表示能量转换效率，φ0＝diag(φ
1,0
,...,φn，0)是下行时隙τ0的ris反射系数矩阵，表示基站的波束赋形矢量，p
ris,sat
和p
(k,m),sat
分别表示ris和用户u
k,m
的饱和功率；
[0021]
根据noma协议，基站使用连续干扰消除来解码每个用户群的信息；具体地，当解码第k个簇的用户信息时，为了检测第m个用户的消息，基站首先解码第i个用户的消息，然后从接收的消息中移除该消息；则其它用户的信号则视为干扰信号；假设用户将获得的能量全部用与上行信息传输，则基站接收的第k簇用户的信号表示为：
[0022][0023]
其中，s
k,m
表示用户u
k,m
的传输数据符号，其均值为0和方差为1；φk表示τk时隙的ris反射系数矩阵；nk表示零均值和方差σ2的加性噪声；第k个簇中第m个用户的可实现吞吐量表示为：
[0024][0025]
其中，p
k,m
表示用户u
k,m
的发射功率；因此，系统总吞吐量由下式给出：
[0026][0027]
优化的原问题表示为最大化系统吞吐量，通过优化上行信息传输时间、下行能量传输时间，基站波束赋形矢量以及不同时隙中ris的反射系数来最大化系统吞吐量；考虑的原问题的目标函数p1表述如下：
[0028][0029]
进一步地，步骤s2中原问题的目标函数p1预设的初始约束条件包括：
[0030][0031]
|φ
n，0
|＝1,n∈n；(9)
[0032]
|φ
n,k
|＝1,n∈n；(10)
[0033]
τ0≥0，τk≥0；(11)
[0034]
w||2≤p
max
；(12)
[0035][0036]
0≤βe≤1；(14)
[0037]
其中w表示基站的波束赋形矢量；τ0为基站通过下行链路传输能量的下行时隙；{τk
}
为用户终端通过上行链路向基站上行传输的上行时隙集合；φ0为是下行时隙τ0的ris反射系数矩阵，φ0＝diag(φ
1,0
,...,φ
n，0
)；{φk
}
为上行时隙集合{τk}中的上行时隙τk的ris反射系数矩阵，φk＝diag(φ
1,k
,...,φ
n，k
)；βe为ris的反射系数；p
max
表示基站最大发射功率；e
ris
表示ris获得的能量；μ表示单位时间ris每个反射元件所消耗的能量，n为ris中反射单元的数量。
[0038]
进一步地，步骤s3将所述原问题简化为四个子问题，分别是：
[0039]
第一子问题p2：给定波束赋形矢量w,时隙τ0、{τk}以及{φk}，优化φ0，令
[0040]
第一子问题p2表示为：
[0041][0042]
s.t.|φ
n，0
|＝1,n∈n；(17)
[0043]
第二子问题p3：给定其他变量优化{φk}，第二子问题p3表示为：
[0044][0045]
s.t.|φ
n,k
|＝1,n∈n；(33)
[0046]
第三子问题p4：给定φ0，{φk}优化波束赋形矢量w，第三子问题p4表示为：
[0047]
[0048][0049]
第四子问题p5：给定φ0，{φk}以及波束赋形矢量w,优化时隙τ0、{τk}，第四子问题p5表示为：
[0050][0051][0052]
其中，第二子问题为可解非凸子问题，第一子问题、第三子问题和第四子问题为不可解非凸子问题，经过步骤s4将第一子问题、第三子问题和第四子问题凸优化处理为可解的凸子问题。
[0053]
进一步地，经过步骤s4得到第一子问题、第二子问题、第三子问题和第四子问题均为可解的子问题，步骤s5所述迭代算法包括以下步骤：
[0054]
s51.初始化w、τ0、{τk}、{φk}、φ0，目标函数为r
(0)
，差阈值ε＝10-4
；
[0055]
s52.求解第一子问题，获得令
[0056]
s53.求解第二子问题，获得令
[0057]
s54.将代入第三子问题，获得w
*
，令w＝w
*
；
[0058]
s55.在使用一维搜索给定令将w
*
，代入第四子问题，获得搜索出使得第四子问题的目标函数最优的和其中p
ris,sat
表示ris的饱和功率；
[0059]
s56.使用w
*
，计算系统吞吐量r
(r 1)
；
[0060]
s57.如果时，可得到系统最大吞吐量，否则返回执行步骤s52-步骤s56。
[0061]
进一步地，得到步骤s52所述第一个子问题p2的凸优化过程为：
[0062]
通过引入松弛变量e，令e≤η|wh(ghβeφ0g
i,(k,m)
g
d,(k,m)
)|2，e≤p
(k,m),sat
；因此，问题p2等价为：
[0063]
p2.1
[0064][0065][0066]
|φ
n，0
|＝1,n∈ν(21)
[0067]
将η|wh(ghβeφ0g
i，(k,m)
g
d,(k,m)
)|2等价表示为：
[0068]
[0069]
其中，q0＝[φ
1,0
,
…
,φ
n,0
]
t
,(.)'表示对元素进行共轭操作，
[0070]a0,(k,m
)＝diag(g
′
i，(k,m)
)g
′w′wtgt
diag(g
i，(k，m)
)(23)
[0071][0072]c0，(k，m)
＝g
′
d，(k，m)w′wt
diag(g
i，(k，m)
)(25)
[0073]
因此，子问题p2.1等价为
[0074]
p2.2
[0075][0076][0077]
rank(x0)＝1(29)
[0078]
tr(e0x0)＝1(30)
[0079]
其中，约束表达式(30)是约束表达式(9)的等价表达式，en满足在第(n,n)个元素为1，其它元素为0的矩阵，即：
[0080][0081]
子问题p2.2是关于e
k,m
的凹函数，约束表达式(27)、(28)、(30)均是仿射的，在不考虑秩1约束表达式(29)的情况下，此时该问题是标准的凸优化问题。
[0082]
进一步地，得到步骤s54所述第三个子问题的凸优化过程为：
[0083]
给定φ0，{φk}优化波束赋形矢量w，时隙τ0、{τk}，则第三个子问题表示为
[0084]
p4p4
[0085]
引入松弛变量ω
k，m
，令ω≤η|wh(ghβeφ0g
i，(k，m)
g
d，(k，m)
)|2，ω≤p
(k，m)，sat
；令hhβeφ
khi，(k，m)
h
d，(k，m)
＝δh；定义w＝wwh，ghβeφ0g
i，(k，m)
g
d，(k，m)
＝δ9，因此问题p4等价为：
[0086]
p4.2
[0087]
s.t.ω≤ηtr(wθg)(42)
[0088]
tr(w)≤p
max
(43)
[0089]
rank(w)＝1(44)
[0090]
ω≤p
(k，m)，sat
(46)
[0091][0092]
在不考虑秩1约束表达式(44)的情况下，此时该问题是标准的半正定规划问题。
[0093]
进一步地，得到步骤s55所述第四个子问题的凸优化过程为：
[0094]
给定φ0，{φk}以及波束赋形矢量w，优化时隙τ0、{τk}，则第四个子问题表示为
[0095]
p5p5
[0096]
为了最大化目标函数，时隙τ0，{τk}的最优解必须满足(8)式等式成立，即即最优解由于所以通过一维搜索，给定τ0，求解{τk}；对于给定τ0，优化{τk}的优化问题表示为：
[0097]
p5.1
[0098][0099]
目标函数(50)是关于{τk}的凹函数，约束表达式(51)都是凸约束，所以问题p5.1是标准的凸优化问题。
[0100]
本技术方案提出了一种ris辅助无线供电通信网络的优化方法，与现有技术相比，本发明技术方案的有益效果是：本技术方案提出了一种ris辅助无线供电通信网络的优化方法，首先构建包括可重构智能反射面的ris辅助无线供电通信网络系统模型，其中，可重构智能反射面集成有大量低成本的无源反射元件，能自适应调整反射信号的相位，从而获得更高的反射波束赋形增益，实现更长的无线能量传输距离。并且通过更全面的考虑系统优化的约束条件，得到较好的系统优化效果，获得最优的系统吞吐量。
附图说明
[0101]
图1为本发明优化方法步骤示意图
[0102]
图2为ris辅助无线供电通信网络系统模型示意图；
[0103]
图3为ris反射元件数量对系统吞吐量性能影响示意图；
[0104]
图4为系统吞吐量与基站最大发射功率关系示意图；
[0105]
图5为迭代算法流程图。
具体实施方式
[0106]
为清楚地说明本发明一种ris辅助无线供电通信网络的优化方法，结合实施例和附图对本发明作进一步说明，但不应以此限制本发明的保护范围。
[0107]
实施例1
[0108]
一种ris辅助无线供电通信网络的优化方法，如图1所示，包括以下步骤：
[0109]
s1：建立ris辅助无线供电通信网络系统模型；
[0110]
s2：确定所述ris辅助无线供电通信网络系统模型优化的原问题，所述原问题是指，在ris辅助无线供电通信网络系统模型的系统参数满足预设的初始约束条件的情况下，
使系统吞吐量最大；
[0111]
s3：将所述原问题简化为若干个子问题，所述子问题需满足对应的子约束条件，若所述子问题的子约束条件存在非凸约束，则该子问题为非凸子问题，将所述非凸子问题分类为可解非凸子问题和不可解非凸子问题；
[0112]
s4：对不可解非凸子问题进行凸优化处理，使其子约束条件的非凸约束全部转化为凸约束，从而将所有不可解非凸子问题转化为可解的凸子问题；
[0113]
s5：建立迭代算法，通过所述迭代算法对所述子问题进行交替迭代求解，得到使ris辅助无线供电通信网络系统模型的系统吞吐量最大时，所述ris辅助无线供电通信网络系统模型模型的各系统参数。
[0114]
本实施例首先构建包括可重构智能反射面的ris辅助无线供电通信网络系统模型，其中，可重构智能反射面集成有大量低成本的无源反射元件，能自适应调整反射信号的相位，从而获得更高的反射波束赋形增益，实现更长的无线能量传输距离。并且通过更全面的考虑系统优化的约束条件，得到较好的系统优化效果，获得最优的系统吞吐量。
[0115]
实施例2
[0116]
一种ris辅助无线供电通信网络的优化方法，如图1所示，包括以下步骤：
[0117]
s1：建立ris辅助无线供电通信网络系统模型；
[0118]
所述ris辅助无线供电通信网络系统模型如图2所示，包括：基站、用户簇、ris；
[0119]
所述用户簇中包括至少一个用户终端，所述基站和用户簇连接通信，基站和可重构智能反射面ris连接通信，ris和用户簇连接通信，实现基站通过下行链路向用户终端传输能量，用户终端通过上行链路向基站传输信息。
[0120]
s2：确定所述ris辅助无线供电通信网络系统模型优化的原问题，所述原问题是指，在ris辅助无线供电通信网络系统模型的系统参数满足预设的初始约束条件的情况下，使系统吞吐量最大；
[0121]
系统参数包括：w,τ0,{τk},φ0,{φk},βe，其中w表示基站的波束赋形矢量；τ0为基站通过下行链路传输能量的下行时隙；{τk}为用户簇中的用户终端通过上行链路向基站上行传输的上行时隙集合；φ0为是下行时隙τ0的ris反射系数矩阵；{φk}为上行时隙集合{τk}中的上行时隙τk的ris反射系数矩阵；βe为ris的反射系数。
[0122]
所述原问题的计算公式为所述原问题的计算公式用于在满足所述预设的初始约束条件的情况下，计算使网络通信系统吞吐量最大的各系统参数的取值。
[0123]
所述预设的初始约束条件包括：
[0124][0125]
|φ
n，0
|＝1,n∈n；(9)
[0126]
|φ
n,k
|＝1,n∈n；(10)
[0127]
τ0≥0，τk≥0；(11)
[0128]
w||2≤p
max
；(12)
[0129]
[0130]
0≤βe≤1； (14)
[0131]
其中w表示基站的波束赋形矢量；τ0为基站通过下行链路传输能量的下行时隙；{τk}为用户终端通过上行链路向基站上行传输的上行时隙集合；φ0为是下行时隙τ0的ris反射系数矩阵，φ0＝diag(φ
1,0
,...,φ
n，0
)；{φk}为上行时隙集合{τk}中的上行时隙τk的ris反射系数矩阵，φk＝diag(φ
1,k
,...,φ
n，k
)；βe为ris的反射系数；p
max
表示基站最大发射功率；e
ris
表示ris获得的能量；μ表示单位时间ris每个反射元件所消耗的能量，n为ris中反射单元的数量。
[0132]
s3：使用块坐标下降法，将所述原问题简化为若干个子问题，所述子问题需满足对应的子约束条件，若所述子问题的子约束条件存在非凸约束，则该子问题为非凸子问题，将所述非凸子问题分类为可解非凸子问题和不可解非凸子问题；
[0133]
步骤s3将所述原问题简化为四个子问题，分别是第一子问题：给定波束赋形矢量w,时隙τ0、{τk}以及{φk}，优化φ0；第二子问题：给定其他变量优化{φk}；第三子问题：给定其他变量优化波束赋形矢量w；第四子问题：给定φ0，{φk}以及波束赋形矢量w,优化时隙τ0、{τk}；其中第二子问题为可解非凸子问题，第一子问题、第三子问题和第四子问题为不可解非凸子问题，经过步骤s4将第一子问题、第三子问题和第四子问题凸优化处理为可解的凸子问题。
[0134]
s4：对不可解非凸子问题进行凸优化处理，使其子约束条件的非凸约束全部转化为凸约束，从而将所有不可解非凸子问题转化为可解的凸子问题；
[0135]
经过步骤s4得到第一子问题、第二子问题、第三子问题和第四子问题均为可解的子问题，步骤s5所述迭代算法如图5所示，包括以下步骤：
[0136]
s51.初始化w、τ0、{τk}、{φk}、φ0，目标函数为r
(0)
，差阈值ε＝10-4
；
[0137]
s52.求解第一子问题，获得令
[0138]
s53.求解第二子问题，获得令
[0139]
s54.将代入第三子问题，获得w
*
，令w＝w
*
；
[0140]
s55.在使用一维搜索给定令将w
*
，代入第四子问题，获得搜索出使得第四子问题的目标函数最优的和其中p
ris,sat
表示ris的饱和功率；
[0141]
s56.使用w
*
，计算系统吞吐量r
(r 1)
；
[0142]
s57.如果时，可得到系统最大吞吐量，否则返回执行步骤s52-步骤s56。
[0143]
s5：建立迭代算法，通过所述迭代算法对所述子问题进行交替迭代求解，得到使ris辅助无线供电通信网络系统模型的系统吞吐量最大时，所述ris辅助无线供电通信网络系统模型模型的各系统参数。
[0144]
其中s52所述第一个子问题和s54所述第三个子问题均是通过引入松弛变量的方法进行凸优化得到的。第一个子问题和第三个子问题均使用内点法和标准的凸优化工具箱
如cvx进行求解。
[0145]
本实施例首先构建包括可重构智能反射面的ris辅助无线供电通信网络系统模型，其中，可重构智能反射面集成有大量低成本的无源反射元件，能自适应调整反射信号的相位，从而获得更高的反射波束赋形增益，实现更长的无线能量传输距离。并且通过更全面的考虑系统优化的约束条件，得到较好的系统优化效果，获得最优的系统吞吐量。
[0146]
实施例3
[0147]
可重构智能反射面(renconfigurable intelligent surface,ris)通过集成大量低成本的无源反射元件，能自适应调整反射信号的相位，从而获得更高的反射波束赋形增益。另外，ris可以实现信号补盲作用，在通信死角处，使用ris可以实现信号精准定位。因此，利用ris的智能反射被认为是提高未来无线通信网络频谱和能量效率的一种有前景的方法，ris的高反射波束赋形增益可以实现提高无线功率传输和信息传输效率，解决无线供电通信网络所遇到的效率较低的瓶颈性难题。在自我可持续的物联网中，设备被期望以不间断的方式运行，并在理论上具有永久寿命。考虑实际ris元件不可忽略的功耗，提出有效的策略来保持ris长时间工作是很重要的。
[0148]
本实施例的一种ris辅助无线供电通信网络的优化方法，包括以下步骤：
[0149]
s1：建立ris辅助无线供电通信网络系统模型；
[0150]
s2：确定所述ris辅助无线供电通信网络系统模型优化的原问题，所述原问题是指，在ris辅助无线供电通信网络系统模型的系统参数满足预设的初始约束条件的情况下，使系统吞吐量最大；
[0151]
s3：将所述原问题简化为若干个子问题，所述子问题需满足对应的子约束条件，若所述子问题的子约束条件存在非凸约束，则该子问题为非凸子问题，将所述非凸子问题分类为可解非凸子问题和不可解非凸子问题；
[0152]
s4：对不可解非凸子问题进行凸优化处理，使其子约束条件的非凸约束全部转化为凸约束，从而将所有不可解非凸子问题转化为可解的凸子问题；
[0153]
s5：建立迭代算法，通过所述迭代算法对所述子问题进行交替迭代求解，得到使ris辅助无线供电通信网络系统模型吞吐量最大时，所述ris辅助无线供电通信网络系统模型模型的各系统参数。
[0154]
其中，步骤s2确定所述原问题的方法为：
[0155]
考虑一个ris辅助的无线供电通信网络，其中ris具有n个反射元件，天线数为t的多天线基站和多个单天线用户，其中一共有个用户，被分为k簇，每一簇有mk个用户。u
k,m
表示第k簇的第m个用户，m∈ω＝{1,...,mk}。系统的工作被分成两个部分，其中包括下行传输能量和上行传输信息。如图1所示，红线表示基站通过下行链路给用户传输能量，蓝线表示用户通过上行链路往基站传输信息。其中，基站与irs之间传输能量的信道表示为从ris到用户u
k,m
传输能量的信道表示为从基站到用户u
k,m
传输能量的信道表示为从用户u
k,m
到irs传输信息的信道表示为从ris到基站传输信息的反射信道表示为从用户到基站传输信息的
信道表示为这里，表示x
×
y维度的复值矩阵的集合。图2是系统时隙图，在时隙τ0中ris和所有用户收集能量，ris将收集的能量一部分作为自身的能耗，另一部分反射给用户，定义βe为ris的反射系数。本实施例考虑更符合实际的非线性能量收集模型，因此ris和用户u
k,m
所获得的能量分别表示为
[0156][0157]ek，m
＝min(η|wh(ghβeφ0g
i，(k，m)
g
d，(km)
)|2，p
(k，m)，sat
)τ0(2)
[0158]
其中η表示能量转换效率，φ0＝diag(φ
1,0
,...,φn，0)是下行时隙τ0的ris反射系数矩阵，表示基站的波束赋形矢量，p
ris,sat
和p
(k,m),sat
分别表示ris和用户u
k,m
的饱和功率。
[0159]
根据noma协议，基站使用连续干扰消除来解码每个用户群的信息。具体地，当解码第k个簇的用户信息时，为了检测第m个用户的消息，基站首先解码第i个用户的消息，然后从接收的消息中移除该消息。则其它用户的信号则视为干扰信号。假设用户将获得的能量全部用与上行信息传输，则基站接收的第k簇用户的信号可以表示为
[0160][0161]
其中，s
k,m
表示用户u
k,m
的传输数据符号，其均值为0和方差为1。φk的定义类似于φ0，表示τk时隙的ris反射系数矩阵。nk表示零均值和方差σ2的加性噪声。第k个簇中第m个用户的可实现吞吐量(比特/赫兹)可以表示为
[0162][0163]
其中，p
k,m
表示用户u
k,m
的发射功率。因此，系统总吞吐量由下式给出
[0164][0165]
优化问题可以表示为最大化系统吞吐量，通过优化上行信息传输时间、下行能量传输时间,基站波束赋形矢量以及不同时隙中ris的反射系数来最大化系统吞吐量。考虑的问题表述如下。
[0166]
(p1)
[0167][0168]
|φ
n，0
|＝1,n∈n；(9)
[0169]
|φ
n,k
|＝1,n∈n；(10)
[0170]
τ0≥0，τk≥0；(11)
[0171]
w||2≤p
max
；(12)
[0172][0173]
0≤βe≤1；(14)
[0174]
问题(p1)为本实施例步骤s2所述原问题，在问题(p1)中，(8)是总时间约束，总时间为1s,(12)是基站最大发射功率约束，(13)是ris能耗约束，其中μ表示单位时间ris每个反射元件所消耗的能量。由于问题(p1)目标函数是非凹的，约束(9)和(10)都是非凸的约束，使得问题(p1)难以求解。下面将给出算法解决问题的近似解。
[0175]
步骤s2将原问题简化为若干个子问题并进行步骤s4凸优化处理的方法为：
[0176]
首先分析ris的反射系数βe，为了满足系统下行能量传输和上行信息传输，根据约束(13)可知，必须满足如果则ris的能量无法满足它正常的工作。在最优情况下，ris能量收集的接收功率不得大于饱和功率。这是因为如果大于饱和功率，可以通过增加ris反射元件的幅度来增加反射功率，而不影响收集功率，这是不实际的。因此，必须满足必须满足所以根据式(13)可以得出因此，必须满足
[0177][0178]
为了确保的可行性，必须满足由于在最优情况下，反射系数必须为上限，这样才能使得反射功率最大化。因此，并且
[0179]
本实施例使用块坐标法将问题(p1)简化为四个子问题进行求解，分别是第一子问题：给定波束赋形矢量w,时隙τ0、{τk}以及{φk},优化φ0；第二子问题：给定其他变量优化{φk}；第三子问题：给定其他变量优化波束赋形矢量w；第四子问题：给定φ0，{φk}以及波束赋形矢量w,优化时隙τ0、{τk}。
[0180]
第一子问题，波束赋形矢量w，时隙τ0、{τk}以及{φk}，优化φ0，令，令则问题(p1)简化为第一子问题为：
[0181]
(p2)
[0182]
s.t.|φ
n，0
|＝1，n∈n；(17)
[0183]
由式(2)可知，目标函数是一个非凹函数，通过引入松弛变量e，令e≤η|wh(ghβeφ0g
i，(k，m)
g
d，(k，m)
)|2，e≤p
(k，m)，sat
。因此，问题(p2)可以等价为
[0184]
(p2.1)
[0185][0186][0187]
|φ
n，0
|＝1，n∈n；(21)
[0188]
对约束表达式(19)不等号右侧的表达式可以等价表示为
[0189][0190]
其中，q0＝[φ
1，0
，
…
，φ
n，0
]
t
，(.)
′
表示对元素进行共轭操作，
[0191]a0，(k，m)
＝diag(g
′
i，(k，m)
)g
′w′wtgt
diag(g
i，(k，m)
)(23)
[0192][0193]c0，(k，m)
＝g
′
d，(k，m)w′wt
diag(g
i，(k，m)
)(25)
[0194]
因此，子问题(p2.1)可以等价为
[0195]
(p2.2)
[0196][0197][0198]
rank(x0)＝1(29)
[0199]
tr(enx0)＝1(30)
[0200]
其中，约束(30)是约束(9)的等价表达式，en满足在第(n，n)个元素为1，其它元素为0的矩阵，即
[0201][0202]
子问题(p2.2)是关于e
k，m
的凹函数，约束(27)、(28)、(30)均是仿射的，但是由于秩1约束的存在，使得子问题(p2.2)仍然不是标准的半正定规划问题，先不考虑秩1约束(29)，此时该问题是标准的凸优化问题，可以使用内点法和标准的凸优化工具箱如cvx进行求解，得到的解使用高斯随机恢复秩为1的解。
[0203]
第二子问题，给定其他变量优化{φk}，则问题(p1)简化为第二子问题为：
[0204]
(p3)
[0205]
s.t.|φ
n，k
|＝1，n∈n；(33)
[0206]
虽然问题(p3)不是凸问题，但是由于目标函数的特殊性，可以得到闭式解，由三角不等式可知
[0207]
|wh(hhβeφ
khi，(k，m)
h
d，(k，m)
)|＝|whhhβeφ
khi，(k，m)
whh
d，(k，m)
|
[0208]
≤|whhhβeφ
khi，(k，m)
| |whh
d，(k，m)
|(34)
[0209]
当且仅当arg(whhhβeφ
khi，(k，m)
)＝arg(whh
d，(k，m)
)＝ψ0时等式成立，arg(z)表示一个向量，每一个元素都是向量z对应元素的相位。令向量，每一个元素都是向量z对应元素的相位。令f
k，m
＝whhhdiag(βeh
i，(k，m)
)。因此，问题(p3)可以等价为
[0210]
(p3.1)
[0211]
s.t.arg(ξhf)＝ψ0(36)
[0212]
|φ
n，k
|＝1，n∈n；(37)
[0213]
不难看出，问题(p3.1)的最优解为因此，第k个时隙的ris第n个反射元件的相位为
[0214][0215]
其中，表示向量h
i，(k，m)
第n个元素，hn表示矩阵h的第n行矢量。
[0216]
第三子问题：给定φ0，{φk}优化波束赋形矢量w，时隙τ0、{τk}，则第三子问题可以表示为
[0217]
(p4)
[0218]
s.t.(12)(13)(40)
[0219]
与第一子问题的求解方法类似，引入松弛变量ω
k，m
，令ω≤η|wh(ghβeφ0g
i，(k，m)
g
d，(k，m)
)|2，ω≤p
(k，m)，sat
。令hhβeφ
khi，(k，m)
h
d，(k，m)
＝δh。定义w＝wwh，ghβeφ0g
i，(k，m)
g
d，(k，m)
＝δg，因此问题(p4)可以等价为
[0220]
(p4.2)
[0221]
s.t.ω≤ηtr(wθg)(42)
[0222]
tr(w)≤p
max
(43)
[0223]
rank(w)＝1(44)
[0224]
ω≤p
(k，m)，sat
(46)
[0225][0226]
与问题(p2.2)的解法一致，先不考虑秩1约束(44)，此时该问题是标准的半正定规划问题，可以使用内点法或者通过标准的凸优化问题求解器如cvx工具箱有效地解决，得到的解再进行高斯随机化求得w的秩为1的解。
[0227]
第四子问题：给定φ0，{φk}以及波束赋形矢量w，优化时隙τ0、{τk}，则第四子问题可以表示为
[0228]
(p5)
[0229]
s.t.(8)(11)(13)(49)为了最大化目标函数，时隙τ0，{τk}的最优解必须满足(8)式等式成立，即即假设最优解可以找到一个可行解和使得目标函数的值比最优解更大，这与最优解相矛盾，所以目标函数的值比最优解更大，这与最优解相矛盾，所以由之前的分析可知所以可以通过一维搜索，给定τ0，求解{τk}。对于给定τ0，优化{τk}的优化问题可以表示为
[0230]
(p5.1)
s.t.(8)(13)(51)
[0231]
目标函数(50)是关于{τk}的凹函数，约束表达式(51)都是凸约束，所以问题(p5.1)是标准的凸优化问题，可以使用内点法或者通过标准的凸优化问题求解器有效地解决。
[0232]
步骤s5建立所述迭代算法
[0233]
(1)初始化w、τ0、{τk}、{φk}、φ0，目标函数为r
(0)
，迭代次数r＝0，差阈值ε＝10-4
；
[0234]
(2)r＝r 1，求解问题p2.2，获得令
[0235]
(3)求解问题p3.1，获得令
[0236]
(4)将代入问题p4.2，获得w
*
，令w＝w
*
；
[0237]
(5)在使用一维搜索给定令将w
*
，代入问题p5.1，获得搜索出使得p5的目标函数(48)最优的和
[0238]
(6)将w
*
，代入(6)式，求得r
(r 1)
[0239]
(7)如果时，可得到系统最大吞吐量，否则重复步骤(2)，(3)，(4)，(5)，(6)。
[0240]
本实施例通过块坐标法将目标问题简化为四个子问题，在求解下行能量传输时隙时ris的反射相位矩阵使用了半正定松弛技术，求解上行信息传输时的ris反射相位，得到了闭式解，降低了算法优化的复杂度。优化时隙时，使用一维搜索可以得到最优解。
[0241]
本实施例提出的算法适用范围广。首先，所考虑的模型采用了一种通用的多址接入方式，tdma和noma都是该方式的特例。其次，本实施例的算法对用户数量、用户簇的数量没有限制，适用于任意用户簇数量及簇内任意用户数量的情况。再次，对于ris和用户的能量收集，与传统研究考虑的线性能量收集模型相比。本实施例考虑了更符合实际的非线性能量收集模型，适用范围更广。
[0242]
可以增加无源反射元件的数量来提高系统的吞吐量性能，提高系统的信息传输效率、能量使用效率和能量传输效率。
[0243]
在上行信息传输过程，使用混合多址接入的方法能更好的满足用户多的情况的通信，即可以提高频谱效率，也能降低在基站处做串行干扰消除的复杂度，做到性能和复杂度的折中。
[0244]
考虑更符合实际的ris和用户的非线性能量收集模型，算法适用范围更广。
[0245]
进行仿真实验得到实验结果，ris反射元件数量对系统吞吐量性能的影响示意图如图3所示，系统吞吐量随基站最大发射功率变化趋势图如图4所示，图中：
[0246]
optimized ris w/ta：优化时隙和每个时隙的不同的ris反射矩阵以及基站波束成形矢量
[0247]
optimized ris w/o ta：只优化每个时隙ris反射矩阵和基站波束成形矢量，固定时隙
[0248]
random ris w/ta：优化时隙和基站波束成形矢量，ris的相位时随机相位
[0249]
random ris w/o ta:随机相位情况下，固定时隙，只优化基站波束成形矢量
[0250]
without ris:没有ris情况。
[0251]
可以从图中看出，提出的算法在性能上与其他算法的对比，性能得到一个提高，并且当反射元件数量增大时，系统的吞吐量也随着增大。
[0252]
显然，本发明的上述实施例仅仅是为清楚地说明本发明所作的举例，而并非是对本发明的实施方式的限定。对于所属领域的普通技术人员来说，在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明权利要求的保护范围之内。