CN113114317A - 一种irs辅助下行多用户通信***的相移优化方法 - Google Patents

一种irs辅助下行多用户通信***的相移优化方法 Download PDF

Info

Publication number
CN113114317A
CN113114317A CN202110396539.3A CN202110396539A CN113114317A CN 113114317 A CN113114317 A CN 113114317A CN 202110396539 A CN202110396539 A CN 202110396539A CN 113114317 A CN113114317 A CN 113114317A
Authority
CN
China
Prior art keywords
matrix
irs
phase shift
user
channel capacity
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.)
Granted
Application number
CN202110396539.3A
Other languages
English (en)
Other versions
CN113114317B (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.)
Chongqing University of Post and Telecommunications
Original Assignee
Chongqing University of Post and Telecommunications
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 Chongqing University of Post and Telecommunications filed Critical Chongqing University of Post and Telecommunications
Priority to CN202110396539.3A priority Critical patent/CN113114317B/zh
Publication of CN113114317A publication Critical patent/CN113114317A/zh
Application granted granted Critical
Publication of CN113114317B publication Critical patent/CN113114317B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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
    • H04B—TRANSMISSION
    • H04B17/00—Monitoring; Testing
    • H04B17/30—Monitoring; Testing of propagation channels
    • H04B17/391—Modelling the propagation channel
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04W—WIRELESS COMMUNICATION NETWORKS
    • H04W24/00—Supervisory, monitoring or testing arrangements
    • H04W24/02—Arrangements for optimising operational condition
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04W—WIRELESS COMMUNICATION NETWORKS
    • H04W24/00—Supervisory, monitoring or testing arrangements
    • H04W24/06—Testing, supervising or monitoring using simulated traffic
    • Y—GENERAL 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
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02D—CLIMATE 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/00—Reducing energy consumption in communication networks
    • Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Mobile Radio Communication Systems (AREA)

Abstract

本发明属于移动通信技术领域，具体涉及到一种IRS辅助下行多用户通信***的相移优化方法包括建立移动用户的IRS辅助通信***模型，并基于最优解码矩阵的均方误差矩阵和信道容量之间的关系，以信道容量最大化为目标构建IRS辅助通信***的信道容量优化模型；求解信道容量优化模型中最优解码矩阵、辅助矩阵以及预编码矩阵，根据信道容量优化模型及其已知矩阵构建IRS相移优化模型；采用黎曼信赖域方法对IRS相移优化模型进行求解，获得最佳相移，并将其作为IRS相移矩阵的主对角元素；本发明在未损失吞吐量的情况下，极大降低复杂度，并且本发明获得的信道容量均大于同等情况下传统方法获得的信道容量。

Description

一种IRS辅助下行多用户通信***的相移优化方法
技术领域
本发明属于移动通信技术领域，具体涉及到一种使智能反射表面(IntelligentReflecting Surfaces，IRS)辅助下行多用户通信***的相移优化方法。
背景技术
随着第五代(5G)无线通信网络的商业化，意味着未来的无限网络的频谱效率更高、传输速率更快、***容量更大，理论上，其频谱效率能够比4G高出5到15倍，其峰值传输速度能够达到10Gbit/s。大规模天线阵列(M-MIMO)和毫米波(mm-Wave)通信等5G关键技术为实现此功能做出了有效贡献。然而，额外的高硬件成本，巨大的功耗或者能源消耗以及5G基站的选址是它们在实践中实施的主要障碍。因此，为了在5G和无线网络之外实现绿色和可持续发展，寻找频谱和高能效技术的研究对于可持续容量增长仍然至关重要。为了解决上述挑战，在5G和下一代移动通信***中，智能反射面(Intelligent ReflectingSurfaces，IRS)被认为是一种有前途的绿色，高成本效益，节能和频谱高效的技术。为了满足新一代无线网络的要求，能够达到更高频谱速率、连接更多的用户数、实现更大的***容量，传统的5G技术无法完美的满足这些要求，智能反射表面技术的提出为这一方向找到了突破口。
IRS是由电磁材料组成的人造表面，可以通过在表面上集成大量低成本的无源反射元件来智能地调整无线传播环境。此外，IRS外形小巧，重量轻，为实际实施提供了很高的灵活性。例如，IRS可以轻松地从墙壁，天花板，广告板甚至衣服上安装/拆卸。并且，与传统的多用户通信***相比，多天线基站使用不同的线性预编码器同时服务不同的用户，这比单天线的基站极大地提高了频谱效率。但是当相邻的用户靠得比较近时，会产生较严重的共信道干扰，那么传统的线性预编码器对消除用户间的干扰无效。但是通过有效地布置IRS可以解决用户间干扰问题，因为IRS通过优化其反射系数来提供额外的控制信号路径，这可以有效的消除相邻用户间的信道干扰。
在IRS辅助通信***中，基站信号在经过IRS时，IRS通过在表面施加外部刺激(外部电流变化)进而从特定方向反射入射电磁波。IRS作为未来移动通信的关键技术，为了使其在实际应用中发挥出更优越的性能，得到更高的频率效率，更高的***容量以及更好的能量效率，对IRS辅助通信***中IRS相移优化研究尤为重要。现有的IRS辅助通信***相移优化算法大都采用半定松弛技术(Semi-Definite Relaxation，SDR)，这种算法虽然可以最大化***吞吐量，但是不仅算法实现复杂度高而且获得的只是一个近似解，在具体实现时可能产生不可避免的误差。
发明内容
为了在总功率限制和IRS相移所引起的单位模量约束的条件下最大化信道容量，本发明提出一种IRS辅助下行多用户通信***的相移优化方法，具体包括以下步骤：
建立移动用户的IRS辅助通信***模型，并基于最优解码矩阵的均方误差矩阵和信道容量之间的关系，以信道容量最大化为目标构建IRS辅助通信***的信道容量优化模型；
求解信道容量优化模型中最优解码矩阵、辅助矩阵以及预编码矩阵，根据信道容量优化模型及其已知矩阵构建IRS相移优化模型；
采用黎曼信赖域方法对IRS相移优化模型进行求解，获得最佳相移，并将其作为IRS相移矩阵的主对角元素。
进一步的，以信道容量最大化为目标构建IRS辅助通信***的信道容量优化模型表示为：
Figure BDA0003018791110000031
约束条件：
Figure BDA0003018791110000032
Figure BDA0003018791110000033
其中，G是预编码矩阵的集合；θ＝[θ1,…,θn,…,θN]是IRS反射相移收集向量；θn是IRS处第n个反射元件的相移；K为移动用户的数量；I为单位矩阵；σ2表示噪声的功率；Pmax表示发射功率阈值；N表示IRS反射元件数量；Gk是用户k的预编码矩阵；HIU,k是IRS与用户k间的信道矩阵；HAI,k是IRS与基站之间的信道矩阵；HAU,k是基站与用户k之间的信道矩阵；(·)H表示求矩阵的共轭转置。
进一步的，利用最优解码矩阵的均方误差和信道容量之间的关系对IRS辅助通信***的信道容量优化模型进行优化，优化后模型表示为：
Figure BDA0003018791110000034
约束条件：
Figure BDA0003018791110000035
Figure BDA0003018791110000036
其中，W是引入辅助矩阵的集合；D是用户解码矩阵的集合；rk(W,D,G,θ)是基于最优解码矩阵的均方误差矩阵和信道容量之间的关系所形成的目标函数；Wk是用户k的辅助矩阵；Ek是用户k的最小均方误差矩阵；q为用户的数据流；Tr( )表示求矩阵的迹；| |表示求矩阵的行列式；|| ||2表示求矩阵的二范数。
进一步的，IRS相移优化模型表示为：
Figure BDA0003018791110000037
约束条件：
Figure BDA0003018791110000038
其中，θ是相移收集向量；
Figure BDA0003018791110000039
是IRS处的相移对角矩阵；
Figure BDA00030187911100000310
K为移动用户的数量，
Figure BDA0003018791110000041
HIU,k是IRS到用户k之间的信道矩阵，Dk是用户k的解码矩阵，Wk用户k的辅助矩阵；
Figure BDA0003018791110000042
HAI,k是用户k的IRS到用户间的信道矩阵，
Figure BDA0003018791110000043
Gi是第i个用户的预编码矩阵；
Figure BDA0003018791110000044
HAU,k是基站与用户k之间的信道矩阵，
Figure BDA0003018791110000045
Gk是用户k的预编码矩阵；Tr()表示求矩阵的迹，上标H表示矩阵的共轭转置。
进一步的，采用黎曼信赖域方法对IRS相移优化模型进行求解包括：
S100：将IRS相移优化模型中相移约束转换为单位模量约束，获取IRS相移优化模型的流形优化，并将单位模量约束条件构成复圆流形嵌入到搜索空间，迭代求解搜索方向；
S101、采用函数二阶泰勒展开构造信赖域子问题，求解该子问题获取第i次迭代的搜索方向；
S102、根据搜索方向计算在第i次迭代时的目标函数实际下降量与根据二次模型函数的预测下降量的比值；
S103、根据获取的比值对信赖域进行调整；
S104、判断步骤S102计算得到的比值是否大于设定的阈值，若大于，则令ζi+1＝Rζi(ηi)，否则令ζi+1＝ζi；
S105、判断||gradf(ζi+1)||2≤ε是否成立，若不成立，则令i＝i+1并返回步骤S101；若成立，则IRS局部最优相移为
Figure BDA0003018791110000046
将其作为IRS相移矩阵
Figure BDA0003018791110000047
的主对角元素，重复以上步骤求取IRS相移矩阵
Figure BDA0003018791110000048
的其他对角元素；
其中，
Figure BDA0003018791110000049
为目标函数梯度为零时的相移；
Figure BDA00030187911100000410
为在点ζi的收缩；gradf(ζi+1)为黎曼梯度，|| ||2表示求矩阵的二范数；ε为迭代停止门限值。
进一步的，IRS相移优化模型的流形优化表示为：
Figure BDA00030187911100000411
其中，f(ζ)为相移优化目标函数；ζ为优化对象，是IRS相移优化模型中相移矩阵
Figure BDA0003018791110000051
对角元素的收集向量，表示为
Figure BDA0003018791110000052
θn为第n个反射元件的相移；U＝B⊙CT为一个半定矩阵，B⊙CT表示矩阵B与矩阵CT的哈达玛乘积；α为一个正数用来控制收敛；IN为N*N元单位阵；v*为矩阵V对角元素的收集向量；Re{}表示取实部。
进一步的，采用函数二阶泰勒展开构造信赖域子问题包括：
Figure BDA0003018791110000053
约束条件：||ηi||≤Δi，
其中，ζi表示第i次的迭代点，即优化对象的第i次迭代；ηi为第i次迭代的搜索方向；
Figure BDA0003018791110000054
为在SN上点ζi的切空间，其中
Figure BDA0003018791110000055
是由N个单位模量约束条件所形成的乘积复圆流形；
Figure BDA0003018791110000056
为二次函数模型；f(ζi)为第i次迭代时的相移优化目标函数；gradf(ζi)为流形上的黎曼梯度；Hessf(ζi)[ηi]为f(ζi)的黎曼海塞矩阵；<a,b>表示a与b的内积，Δi为信赖域半径；||ηi||为第i次迭代搜索方向的二范数。
进一步的，在第i次迭代时的实际下降量与根据二次模型函数的预测下降量的比值表示为：
Figure BDA0003018791110000057
其中，f(ζi)为点ζi的目标函数值；
Figure BDA0003018791110000058
为点ζi经过收缩后的目标函数值；
Figure BDA0003018791110000059
为0向量的二次模型函数值；
Figure BDA00030187911100000510
为搜索方向ηi的二次模型函数值。
进一步的，根据获取的比值对信赖域进行调整包括：
Figure BDA0003018791110000061
其中，Δi+1为第i+1次迭代的信赖域半径；
Figure BDA0003018791110000062
为最大信赖域半径，设定为
Figure BDA0003018791110000063
Figure BDA0003018791110000064
根据搜索方向计算在第i次迭代时的目标函数实际下降量与根据二次模型函数的预测下降量的比值；||ηi||为第i次迭代搜索方向的二范数。
本发明的有益效果包括：
(1)、与传统虽然能获得较大的***总吞吐量、但复杂度较高且为近似解的SDR相移优化算法相比，本发明在未损失吞吐量的情况下，杂度仅为O(N3)而SDR算法复杂度为O(N6)，在极大降低复杂度的情况下本发明仍保持优异性能；
(2)、本发明优化相移所获得的信道容量均达大于相同情况下的SDR算法、随机相移以及without-IRS***；
(3)、本发明方法相比without-IRS***的信道容量高出25.65％；比随机相移方法的信道容量依旧高出20.09％；比SDR方法的***容量高出2.5％。
附图说明
图1为本发明采用的黎曼信赖域算法流程图；
图2为本发明一种IRS辅助下行多用户通信***的相移优化方法流程图；
图3为本发明以及现有技术的仿真图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
本发明提出一种IRS辅助下行多用户通信***的相移优化方法，如图1和图2，具体包括以下步骤：
建立移动用户的IRS辅助通信***模型，并基于最优解码矩阵的均方误差矩阵和信道容量之间的关系，以信道容量最大化为目标构建IRS辅助通信***的信道容量优化模型；
求解信道容量优化模型中最优解码矩阵、辅助矩阵以及预编码矩阵，根据信道容量优化模型及其已知矩阵构建IRS相移优化模型；
采用黎曼信赖域方法对IRS相移优化模型进行求解，获得最佳相移，并将其作为IRS相移矩阵的主对角元素。
在本实施例中为了方便分析，假设基站功率Pmax恒为1W，小区用户数为5人，建立具有M1＝10根天线的基站，K＝5个M2＝2根天线的移动用户的IRS辅助通信***模型，IRS由N＝20个反射元件构成，其中水平方向上的反射元件数为5，垂直方向上的反射元件数为4；
设置
Figure BDA0003018791110000071
表示基站(Base Station,BS)和IRS之间的信道矩阵，
Figure BDA0003018791110000072
表示的IRS和用户之间的信道矩阵，
Figure BDA0003018791110000073
表示的BS和用户之间的信道矩阵；
Figure BDA0003018791110000074
是加性高斯白噪声；
Figure BDA0003018791110000075
是IRS处的对角矩阵，其中
Figure BDA0003018791110000076
α∈(0,1]是振幅系数，设置为最大化信号反射即α＝1，θn∈[0,2π)是相移；
Figure BDA0003018791110000077
表示在AP处的发射信号，
Figure BDA0003018791110000078
表示用户k的发送信号，
Figure BDA0003018791110000079
和
Figure BDA00030187911100000710
分别表示用户k的线性预编码矩阵和信息承载符号数据，其中q是每个用户的数据流。同时，由于所有用户发射信号之和应小于最大发送功率阈值Pmax，即
Figure BDA00030187911100000711
假定所有涉及的信道都是独立的瑞利衰落，并且在10m的参考距离处将路径损耗指数设置为3，基站与IRS之间距离为设置为50米，IRS坐标为(50,0)，基站坐标为(0,0)，ε＝10-4；以上述参数为基础，建立IRS辅助通信***的信道容量优化问题模型：
Figure BDA0003018791110000081
Figure BDA0003018791110000082
Figure BDA0003018791110000083
其中，式(1)表示信道容量最大化问题模型，式(2)表示有用户发射信号之和必须小于最大发射功率Pmax，式(3)表示相移约束；
G为所有用户预编码矩阵的集合，表示为
Figure BDA0003018791110000084
θ为相移的集合，表示为
Figure BDA0003018791110000085
I为单位矩阵；σ2表示噪声的功率。
上述优化问题(1)是一个NP难问题，需要将原问题转化为一个容易处理的形式。利用最优解码矩阵的均方误差(Mean-Square Error，MSE)矩阵和信道容量之间的关系，将将式(1)改成：
Figure BDA0003018791110000086
约束条件：
Figure BDA0003018791110000087
Figure BDA0003018791110000088
其中，
Figure BDA0003018791110000089
是引入的辅助矩阵的集合，≥0表示矩阵W为一个正半定矩阵；
Figure BDA00030187911100000810
是的解码矩阵的集合，θ＝[θ1,...,θN]是一个相移收集向量，E＝{Ek},
Figure BDA00030187911100000811
是所有用户MSE矩阵的集合，q是每个用户的数据流，第k个用户的MSE矩阵Ek如下式所示：
Figure BDA00030187911100000812
其中，
Figure BDA00030187911100000813
是基站到第k个用户的等效信道；yk表示第k个用户的基带接收信号；Gi为第i个用户的预编码矩阵；E[]表示求期望。
式(4)中的W与D只与rk(W,D,G,θ)有关，可以通过固定其余三个变量，对W或D的rk(W,D,G,θ)求一阶偏导并令其为零时获得最优解。那么第k个用户的最优解码矩阵与辅助矩阵分别是：
Figure BDA0003018791110000091
Figure BDA0003018791110000092
对于式(4)中的预编码矩阵G，可以通过复杂的矩阵变换，将问题转化为凸优化问题，再将其转化为二阶锥规划问题，可以直接使用标准的优化包，例如CVX进行求解。
通过固定矩阵W,D,G，将式(7)中Ek代入式(4)，并忽略常数项，那么相移优化问题如下：
Figure BDA0003018791110000093
Figure BDA0003018791110000094
其中，Tr()表示求矩阵的迹；令
Figure BDA0003018791110000095
将其和
Figure BDA0003018791110000096
并代入式(10)展开获得：
Figure BDA0003018791110000097
Figure BDA0003018791110000098
根据式(12)、(13)，定义
Figure BDA0003018791110000099
Figure BDA00030187911100000910
并将这些等式代入式(10)，并忽视与相移无关的常数项，则IRS相移优化模型可以变形为：
Figure BDA0003018791110000101
Figure BDA0003018791110000102
其中，
Figure BDA0003018791110000103
Figure BDA0003018791110000104
是一个对角矩阵，定义
Figure BDA0003018791110000105
同时令v＝[[V]1,1,...,[V]N,N]H为矩阵V对角元素的收集向量。通过利用矩阵性质可获得：
Figure BDA0003018791110000106
并且由
Figure BDA0003018791110000107
将相移约束转换为单位模量约束，即|ζn|＝1。那么问题式(14)可以变形为：
Figure BDA0003018791110000108
Figure BDA0003018791110000109
其中，f(ζ)为目标优化函数；v*为矩阵V对角元素的收集向量；Re{}为求实部；U＝B⊙CT是由B和C矩阵是半定矩阵，那么U也是一个半定矩阵，
Figure BDA00030187911100001010
是一个正数用来控制算法收敛性；λU为矩阵U的特征值。
基于上述问题，本文采用黎曼信赖域方法对其求解。包括：
构建关于式(16)的复圆流形：
Figure BDA00030187911100001011
其中
Figure BDA00030187911100001012
是一个复圆流形，该集合S可以看作是
Figure BDA00030187911100001013
的子流形，因此，这N个圆的乘积SN是
Figure BDA00030187911100001014
的子流形。将单位模量约束条件构成复圆流形嵌入到搜索空间，那么该优化问题与在欧氏空间的信赖域方法求解相似。初始化：
Figure BDA00030187911100001015
以及ζi∈SN，其中
Figure BDA00030187911100001016
为最大信赖域半径，通常设定为
Figure BDA00030187911100001017
Δ0为初始信赖域半径。
对于该最小化问题，在第i次迭代时首先要找到搜索方向。最常见的寻找搜索方向的方法是沿与f(ζi)的欧氏梯度相反的方向移动，但这会导致收敛速度较慢。本文利用二次逼近构造信赖域子问题求搜索方向，由下式给出：
Figure BDA0003018791110000111
s.t.||ηi||≤Δi. (20)
其中，
Figure BDA0003018791110000112
是流形上的黎曼梯度，它在当前迭代点ζi的切空间
Figure BDA0003018791110000113
上，
Figure BDA0003018791110000114
为在SN上点ζi的切空间，其中
Figure BDA0003018791110000115
是由N个单位模量约束条件所形成的乘积复圆流形；
Figure BDA0003018791110000116
为欧氏梯度，
Figure BDA0003018791110000117
欧氏梯度的共轭，⊙表示哈达玛乘积；Hessf(ζi)＝2(U+αIN)是f(ζi)的海塞矩阵；<a,b>表示a与b的内积，Δi为第i次迭代的信赖域半径，在每次迭代过程中会进行更新。对于该信赖域子问题，可以采用预条件共轭梯度法来求解。
Figure BDA0003018791110000118
的定义是在第i次迭代时的实际下降量
Figure BDA0003018791110000119
与二次模型函数
Figure BDA00030187911100001110
的预测下降量
Figure BDA00030187911100001111
的比值。具体如下式：
Figure BDA00030187911100001112
其中，
Figure BDA00030187911100001113
表示指数收缩,在黎曼流形优化过程中用来生成下一个迭代变量ζi+1。
计算所得的比值
Figure BDA00030187911100001114
根据以下原则调整信赖域：
Figure BDA00030187911100001115
根据上述步骤计算的
Figure BDA00030187911100001116
给定一个区间
Figure BDA00030187911100001117
若
Figure BDA00030187911100001118
那么
Figure BDA00030187911100001119
否则ζi+1＝ζi。
判断||gradf(ζi+1)||2≤ε，若不成立，则i＝i+1并返回采用函数二阶泰勒展开构造信赖域子问题，求解该子问题获取第i+1次迭代的搜索方向；若成立，则可以获得的IRS局部最优相移
Figure BDA0003018791110000121
其作为IRS相移矩阵
Figure BDA0003018791110000122
的主对角元素。按照该方法依次求解相移矩阵
Figure BDA0003018791110000123
的所有主对角元素。
如图3，根据TR算法，当用户均在水平距离基站60米的单位圆内，即在坐标(60,0)附近的单位圆内。可以计算出此时的信道容量约为28.9bps/Hz，相同条件下TR算法相比without-IRS***的信道容量高出25.65％；比随机相移方法的信道容量依旧高出20.09％；比SDR方法的***容量高出2.5％，在极大降低复杂度的情况下仍提高了IRS辅助通信***的信道容量。
尽管已经示出和描述了本发明的实施例，对于本领域的普通技术人员而言，可以理解在不脱离本发明的原理和精神的情况下可以对这些实施例进行多种变化、修改、替换和变型，本发明的范围由所附权利要求及其等同物限定。

Claims (9)

1.一种IRS辅助下行多用户通信***的相移优化方法，其特征在于，包括以下步骤：
建立移动用户的IRS辅助通信***模型，并基于最优解码矩阵的均方误差矩阵和信道容量之间的关系，以信道容量最大化为目标构建IRS辅助通信***的信道容量优化模型；
求解信道容量优化模型中最优解码矩阵、辅助矩阵以及预编码矩阵，根据信道容量优化模型及其已知矩阵构建IRS相移优化模型；
采用黎曼信赖域方法对IRS相移优化模型进行求解，获得最佳相移，并将其作为IRS相移矩阵的主对角元素。
2.根据权利要求1所述的一种IRS辅助下行多用户通信***的相移优化方法，其特征在于，以信道容量最大化为目标构建IRS辅助通信***的信道容量优化模型表示为：
Figure FDA0003018791100000011
约束条件：
Figure FDA0003018791100000012
Figure FDA0003018791100000013
其中，G是预编码矩阵的集合；θ＝[θ1,...,θn,...,θN]是IRS反射相移收集向量；θn是IRS处第n个反射元件的相移；K为移动用户的数量；I为单位矩阵；σ2表示噪声的功率；
Figure FDA0003018791100000014
是IRS处的反射相移对角矩阵；Pmax表示发射功率阈值；N表示IRS反射元件数量；Gk是用户k的预编码矩阵；HIU,k是IRS与用户k之间的信道矩阵；HAI,k是IRS与基站之间的信道矩阵；HAU,k是基站与用户k之间的信道矩阵；(·)H表示求矩阵的共轭转置。
3.根据权利要求2所述的一种IRS辅助下行多用户通信***的相移优化方法，其特征在于，利用最优解码矩阵的均方误差和信道容量之间的关系对IRS辅助通信***的信道容量优化模型进行优化，优化后模型表示为：
Figure FDA0003018791100000021
约束条件：
Figure FDA0003018791100000022
Figure FDA0003018791100000023
其中，W是引入辅助矩阵的集合；D是用户解码矩阵的集合；rk(W,D,G,θ)是基于最优解码矩阵的均方误差矩阵和信道容量之间的关系所形成的目标函数；Wk是用户k的辅助矩阵；Ek是用户k的最小均方误差矩阵；q为用户的数据流；Tr()表示求矩阵的迹；| |表示求矩阵的行列式；|| ||表示求矩阵的二范数。
4.根据权利要求1所述的一种IRS辅助下行多用户通信***的相移优化方法，其特征在于，IRS相移优化模型表示为：
Figure FDA0003018791100000024
约束条件：
Figure FDA0003018791100000025
其中，θ是IRS相移收集向量；
Figure FDA0003018791100000026
是IRS处的反射相移对角矩阵；
Figure FDA0003018791100000027
K为移动用户的数量，
Figure FDA0003018791100000028
HIU,k是用户k的IRS到用户间的信道矩阵，Dk是用户k的解码矩阵，Wk用户k的辅助矩阵；
Figure FDA0003018791100000029
HAI,k是用户k的IRS到用户间的信道矩阵，
Figure FDA00030187911000000210
Gi是第i个用户的预编码矩阵；
Figure FDA00030187911000000211
HAU,k是用户k的AP到用户间的信道矩阵，
Figure FDA00030187911000000212
Gk是用户k的预编码矩阵；Tr()表示求矩阵的迹，上标H表示矩阵的共轭转置。
5.根据权利要求1所述的一种IRS辅助下行多用户通信***的相移优化方法，其特征在于，采用黎曼信赖域方法对IRS相移优化模型进行求解包括：
S100：将IRS相移优化模型中相移约束转换为单位模量约束，获取IRS相移优化模型的流形优化，并将单位模量约束条件构成复圆流形嵌入到搜索空间，迭代求解搜索方向；
S101、采用函数二阶泰勒展开构造信赖域子问题，求解该子问题获取第i次迭代的搜索方向；
S102、根据搜索方向计算在第i次迭代时的相移优化目标函数实际下降量与根据二次模型函数的预测下降量的比值；
S103、根据获取的比值对信赖域进行调整；
S104、判断步骤S102计算得到的比值是否大于设定的阈值，若大于，则令
Figure FDA0003018791100000031
否则令ζi+1＝ζi；
S105、判断||gradf(ζi+1)||2≤ε是否成立，若不成立，则令i＝i+1并返回步骤S101；若成立，则IRS局部最优相移为
Figure FDA0003018791100000032
将其作为IRS相移矩阵
Figure FDA0003018791100000033
的主对角元素，重复以上步骤求取IRS相移矩阵
Figure FDA0003018791100000034
的其他对角元素；
其中，
Figure FDA0003018791100000035
为目标函数梯度为零时的相移；
Figure FDA0003018791100000036
为在点ζi的收缩；gradf(ζi+1)为黎曼梯度，|| ||2表示求矩阵的二范数；ε为迭代停止门限值。
6.根据权利要求5所述的一种IRS辅助下行多用户通信***的相移优化方法，其特征在于，IRS相移优化模型的流形优化表示为：
Figure FDA0003018791100000037
其中，f(ζ)为相移优化目标函数；ζ为优化对象，即IRS相移优化模型中相移矩阵
Figure FDA0003018791100000038
对角元素的收集向量，表示为
Figure FDA0003018791100000039
θn为第n个反射元件的相移；U＝B⊙CT为一个半定矩阵，B⊙CT表示矩阵B与矩阵CT的哈达玛乘积；α为一个正数用来控制收敛；IN为N*N元单位阵；v*为矩阵V对角元素的收集向量；Re{}表示取实部。
7.根据权利要求5所述的一种IRS辅助下行多用户通信***的相移优化方法，其特征在于，采用函数二阶泰勒展开构造信赖域子问题包括：
Figure FDA0003018791100000041
约束条件：||ηi||≤Δi，
其中，ζi表示第i次的迭代点；ηi为第i次迭代的搜索方向；
Figure FDA0003018791100000042
为在SN上点ζi的切空间，其中
Figure FDA0003018791100000043
是由N个单位模量约束条件所形成的乘积复圆流形；
Figure FDA0003018791100000044
为二次函数模型；f(ζi)为第i次迭代时的相移优化目标函数；gradf(ζi)为流形上的黎曼梯度；Hessf(ζi)[ηi]为f(ζi)的黎曼海塞矩阵；<a,b>表示a与b的内积，Δi为信赖域半径；||ηi||为第i次迭代搜索方向的二范数。
8.根据权利要求5所述的一种IRS辅助下行多用户通信***的相移优化方法，其特征在于，在第i次迭代时的实际下降量与根据二次模型函数的预测下降量的比值表示为：
Figure FDA0003018791100000045
其中，f(ζi)为点ζi的目标函数值；
Figure FDA0003018791100000046
为点ζi经过收缩后的目标函数值；
Figure FDA0003018791100000047
为0向量的二次模型函数值；
Figure FDA0003018791100000048
为搜索方向ηi的二次模型函数值。
9.根据权利要求5所述的一种IRS辅助下行多用户通信***的相移优化方法，其特征在于，根据获取的比值对信赖域进行调整包括：
Figure FDA0003018791100000049
其中，Δi+1为第i+1次迭代的信赖域半径；
Figure FDA00030187911000000410
为最大信赖域半径，设定为
Figure FDA00030187911000000411
Figure FDA00030187911000000412
根据搜索方向计算在第i次迭代时的相移优化目标函数实际下降量与根据二次模型函数的预测下降量的比值；||ηi||为第i次迭代搜索方向的二范数。
CN202110396539.3A 2021-04-13 2021-04-13 一种irs辅助下行多用户通信***的相移优化方法 Active CN113114317B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110396539.3A CN113114317B (zh) 2021-04-13 2021-04-13 一种irs辅助下行多用户通信***的相移优化方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110396539.3A CN113114317B (zh) 2021-04-13 2021-04-13 一种irs辅助下行多用户通信***的相移优化方法

Publications (2)

Publication Number Publication Date
CN113114317A true CN113114317A (zh) 2021-07-13
CN113114317B CN113114317B (zh) 2022-05-17

Family

ID=76716616

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110396539.3A Active CN113114317B (zh) 2021-04-13 2021-04-13 一种irs辅助下行多用户通信***的相移优化方法

Country Status (1)

Country Link
CN (1) CN113114317B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113709755A (zh) * 2021-08-25 2021-11-26 武汉大学 一种基于ris技术的异构网络公平共存方法
CN113852402A (zh) * 2021-09-06 2021-12-28 电子科技大学 一种irs辅助的noma-mimo大容量接入方法
CN114520980A (zh) * 2022-01-12 2022-05-20 山西大学 一种三维ris辅助的煤矿通感一体化网络优化方法及装置
CN115412141A (zh) * 2022-08-18 2022-11-29 南京邮电大学 一种irs辅助的空移键控调制***的相移优化方法
CN116319199A (zh) * 2023-03-02 2023-06-23 华南理工大学 无线功率通信网络最大吞吐量闭式解求解方法、装置及介质

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20130272445A1 (en) * 2010-12-07 2013-10-17 Huawei Technologies Co., Ltd. Precoding method and apparatus
CN105868162A (zh) * 2016-03-31 2016-08-17 北京工业大学 一种黎曼流形上的快速优化方法
CN111355520A (zh) * 2020-03-10 2020-06-30 电子科技大学 一种智能反射表面辅助的太赫兹安全通信***设计方法
CN111786708A (zh) * 2020-07-02 2020-10-16 电子科技大学 大规模mimo***的联合信道信息获取方法
CN111901802A (zh) * 2020-07-31 2020-11-06 东南大学 一种借助智能反射表面的miso***下行保密率优化方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20130272445A1 (en) * 2010-12-07 2013-10-17 Huawei Technologies Co., Ltd. Precoding method and apparatus
CN105868162A (zh) * 2016-03-31 2016-08-17 北京工业大学 一种黎曼流形上的快速优化方法
CN111355520A (zh) * 2020-03-10 2020-06-30 电子科技大学 一种智能反射表面辅助的太赫兹安全通信***设计方法
CN111786708A (zh) * 2020-07-02 2020-10-16 电子科技大学 大规模mimo***的联合信道信息获取方法
CN111901802A (zh) * 2020-07-31 2020-11-06 东南大学 一种借助智能反射表面的miso***下行保密率优化方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
LI YOU 等: "Energy Efficiency and Spectral Efficiency Tradeoff in RIS-Aided Multiuser MIMO Uplink Transmission", 《IEEE TRANSACTIONS ON SIGNAL PROCESSING》 *
ZHENGFENG LI等: "Weighted Sum-Rate Maximization for Multi-IRS Aided Cooperative Transmission", 《IEEE WIRELESS COMMUNICATIONS LETTERS》 *

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113709755A (zh) * 2021-08-25 2021-11-26 武汉大学 一种基于ris技术的异构网络公平共存方法
CN113709755B (zh) * 2021-08-25 2022-06-07 武汉大学 一种基于ris技术的异构网络公平共存方法
CN113852402A (zh) * 2021-09-06 2021-12-28 电子科技大学 一种irs辅助的noma-mimo大容量接入方法
CN113852402B (zh) * 2021-09-06 2023-04-18 电子科技大学 一种irs辅助的noma-mimo大容量接入方法
CN114520980A (zh) * 2022-01-12 2022-05-20 山西大学 一种三维ris辅助的煤矿通感一体化网络优化方法及装置
CN115412141A (zh) * 2022-08-18 2022-11-29 南京邮电大学 一种irs辅助的空移键控调制***的相移优化方法
CN115412141B (zh) * 2022-08-18 2023-07-25 南京邮电大学 一种irs辅助的空移键控调制***的相移优化方法
CN116319199A (zh) * 2023-03-02 2023-06-23 华南理工大学 无线功率通信网络最大吞吐量闭式解求解方法、装置及介质
CN116319199B (zh) * 2023-03-02 2024-05-31 华南理工大学 无线功率通信网络最大吞吐量闭式解求解方法、装置及介质

Also Published As

Publication number Publication date
CN113114317B (zh) 2022-05-17

Similar Documents

Publication Publication Date Title
CN113114317B (zh) 一种irs辅助下行多用户通信***的相移优化方法
CN112290995B (zh) 星地集成网络中基于安全能效的波束设计方法
Ni et al. Enhancing downlink transmission in MIMO HetNet with wireless backhaul
CN112383332A (zh) 一种基于智能反射表面的蜂窝基站通信***
CN111835406B (zh) 适用于多波束卫星通信的能效谱效权衡的鲁棒预编码方法
CN110011712B (zh) 一种面向毫米波大规模多输入多输出的混合预编码方法
CN113691295B (zh) 一种基于irs的异构网络中干扰抑制方法
CN114070365B (zh) 一种智能反射表面辅助的低射频复杂度多用户mimo上行频谱效率优化方法
CN113497649B (zh) 基于智能反射平面的太赫兹无线通信网络资源管控方法
CN104320169B (zh) 多用户3d‑mimo***中三维波束赋形设计方法
CN110138427A (zh) 基于部分连接的大规模多输入多输出混合波束赋形算法
CN111901812A (zh) 一种全双工蜂窝通信网基站和智慧反射面联合控制方法
CN104981004A (zh) 基于多用户双向中继***的收发器能效优化方法及装置
CN107241799B (zh) 异构网络中能量有效的资源分配和干扰对齐联合方法
CN113067610A (zh) 基于swipt的毫米波noma***安全能效最大化功率分配方法
Lu et al. Joint power control and passive beamforming in reconfigurable intelligent surface assisted user-centric networks
CN108337024B (zh) 一种基于能量采集的大规模mimo***能效优化方法
CN108521672B (zh) 一种分布式无线能量和信息传输***的资源分配方法
CN115379478B (zh) 一种基于ris辅助数能同传网络鲁棒能耗优化方法
CN115882911A (zh) 一种非理想硬件条件下的多用户通信***物理层控制方法
CN113890634B (zh) 一种智能反射面辅助的干扰抵消波束设计方法
CN113922849B (zh) 一种毫米波mimo-noma***下的用户分组和功率分配方法
CN116056210A (zh) 一种面向容量覆盖的irs辅助超密集网络资源分配方法
CN110048753A (zh) 基于mimo***加权能效最大化的分布式波束赋形优化方法
Mahesh et al. Design and Performance Analysis of Massive MIMO Modeling with Reflected Intelligent Surface to Enhance the Capacity of 6G Networks

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