CN116170102A - 一种可重构的多功能智能超表面及其应用 - Google Patents

一种可重构的多功能智能超表面及其应用 Download PDF

Info

Publication number
CN116170102A
CN116170102A CN202310187291.9A CN202310187291A CN116170102A CN 116170102 A CN116170102 A CN 116170102A CN 202310187291 A CN202310187291 A CN 202310187291A CN 116170102 A CN116170102 A CN 116170102A
Authority
CN
China
Prior art keywords
ris
sub
user
base station
penalty
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.)
Pending
Application number
CN202310187291.9A
Other languages
English (en)
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.)
Beijing University of Posts and Telecommunications
Original Assignee
Beijing University of Posts 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 Beijing University of Posts and Telecommunications filed Critical Beijing University of Posts and Telecommunications
Priority to CN202310187291.9A priority Critical patent/CN116170102A/zh
Publication of CN116170102A publication Critical patent/CN116170102A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B17/00Monitoring; Testing
    • H04B17/30Monitoring; Testing of propagation channels
    • H04B17/391Modelling the propagation channel
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B17/00Monitoring; Testing
    • H04B17/30Monitoring; Testing of propagation channels
    • H04B17/391Modelling the propagation channel
    • H04B17/3911Fading models or fading generators
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B7/00Radio transmission systems, i.e. using radiation field
    • H04B7/02Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
    • H04B7/04Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
    • H04B7/06Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station
    • H04B7/0613Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station using simultaneous transmission
    • H04B7/0615Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station using simultaneous transmission of weighted versions of same signal
    • H04B7/0617Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station using simultaneous transmission of weighted versions of same signal for beam forming
    • 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

Landscapes

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

Abstract

本发明公开了一种可重构的多功能智能超表面及其应用,多功能智能超表面可用于主动地改变无线信号传播环境,从而提升无线通信系统的信号传输性能。针对MF‑RIS辅助的非正交多址接入网络中的资源分配问题,本发明以最大化用户和速率为目标,构建了一个联合设计基站波束赋形与MF‑RIS系数的非凸优化问题;针对所构建的优化问题,首先采用惩罚函数法和连续凸逼近法将原始优化问题转化为多个半定规划子问题,然后使用交替优化方法获取基站波束赋形向量和MF‑RIS系数的次优解,从而最大化NOMA系统中的用户和速率。

Description

一种可重构的多功能智能超表面及其应用
技术领域
本发明涉及无线资源分配技术领域,尤其涉及一种可重构的多功能智能超表面(Multi-Functional Reconfigurable Intelligence Surface,MF-RIS)及其应用。
背景技术
与正交多址(Orthgonal Multiple Access,OMA)相比,非正交多址接入(Non-Orthgonal Multiple Access,NOMA)能够实现更高的频谱效率并支持大规模连接。先前有研究表明,用户信道条件之间的差异可以用于增强NOMA系统性能。然而,大规模网络中的用户可能具有较差或相似的信道条件,这阻碍了串行干扰消除(Successive InterferenceCancellation,SIC)的应用和NOMA的有效部署。因此,在实际网络中,调整信道条件和增强信道差异能够有效地发挥NOMA的潜力。
最近,凭借重构无线传播环境的能力,智能超表面(ReconfigurableIntelligence Surface,RIS)已成为提高NOMA网络性能的关键技术。通过适当地设计反射系数,RIS能够灵活地调整组合信道条件以增强用户之间的差异,从而提高NOMA在大规模网络中的性能。NOMA和RIS的结合在相关文献中已经进行了大量的研究,并且验证了RIS辅助NOMA网络在增强信号强度、保密通信和降低发射功率方面的优势,也对比了RIS辅助的NOMA和OMA网络中的功耗问题。然而,现有的关于RIS辅助NOMA网络的技术大多使用只支持信号反射的单功能RIS(Single Functional RIS,SF-RIS)。这意味着位于另一侧的用户不能被SF-RIS有效地服务。
为了克服这一限制,提出了双功能RIS(Dual-Functional RIS,DF-RIS)的概念。与SF-RIS不同,DF-RIS是指可以同时进行信号反射和折射的可重构双功能表面,如同时折射和反射RIS(Simultaneous Transmitting and Reflecting RIS,STAR-RIS)和智能全向表面(Intelligent Omni-Surface,IOS)。由于其支持全空间覆盖的有利特性,DF-RIS吸引了很多关注。然而由DF-RIS中继的信号需要通过级联信道,这导致严重的双重衰落问题。这个问题可能严重限制RIS辅助NOMA网络的性能,因为级联信道的质量下降,使得信道之间的差异不像预期的那样明显。因此,有必要设计新的RIS架构来改善无线信道的信道质量,同时减轻现有无源RIS所面临的双重衰落问题。
发明内容
本发明针对现有级联信道的质量下降、双重衰落等问题,提出一种可重构的多功能智能超表面及其在NOMA网络资源分配中的应用,构建得到的多功能智能超表面,不仅可以基于场等效原理将入射信号分为折射和反射两部分,还可以借助有源负载放大输入信号。因此,MF-RIS能够实现全空间覆盖并克服双重衰落问题。
为了实现上述目的,本发明提供如下技术方案:
一方面,本发明提供了一种可重构的多功能智能超表面,其基于场等效原理,利用表面等效电阻抗和磁阻抗实现入射信号的反射和折射,利用负电阻元件实现信号放大。
另一方面,本发明还提供了上述的可重构的多功能智能超表面在NOMA网络资源分配中的应用,包括以下步骤:
S1、设计包含基站有源波束赋形向量和MF-RIS的系数,以最大化所有用户的可实现和速率为目标,构建优化问题和约束条件;
S2、根据步骤S1构建的约束条件确定初始用户和速率表达式,对于其中高度耦合的MF-RIS系数和基站有源波束赋形向量,利用连续凸逼近和惩罚函数法转换为半定规划子问题;
S3、对步骤S2转换后的半定规划子问题进行交替优化,当和速率变化量满足收敛条件,得到有源波束赋形向量和MF-RIS的系数,从而实现最优性能。
进一步地,步骤S1优化问题为:
Figure BDA0004104335590000021
约束条件为:
Figure BDA0004104335590000031
Figure BDA0004104335590000032
Figure BDA0004104335590000033
Figure BDA0004104335590000034
其中M表示MF-RIS的元件个数,N表示基站天线的个数,K表示用户个数,
Figure BDA0004104335590000035
分别表示基站-用户、基站-RIS和RIS-用户的信道,
Figure BDA0004104335590000036
Figure BDA0004104335590000037
分别表示第m个元件的幅度和相移响应,βmax表示最大放大系数,Pmax表示基站的最大发射功率,Po表示MF-RIS的最大放大功率,Rk表示用户k的数据速率,
Figure BDA0004104335590000038
表示用户k的最小速率要求,wk表示用户k的波束向量,Θk表示用户k的MF-RIS的系数矩阵,
Figure BDA0004104335590000039
表示MF-RIS处引入的加性高斯白噪声,hk为基站-MF-RIS-用户的等效组合信道。
进一步地,用户k的波束向量wk和用户k的MF-RIS的系数矩阵Θk在优化问题中是高度耦合MF-RIS系数,对wk和Θk分别优化。
进一步地,对于wk子问题中的目标函数,步骤S2半定规划子问题转化方法为:利用辅助变量Ak和Bk代替wk和hk,等价的表示有源波束赋形问题的表达式,引入
Figure BDA00041043355900000310
Dk=(HHΘk)(HHΘk)H,Wk=wkwk H,对表达式进行进一步的表示,其中Wk≥0且存在非凸秩一约束rank(Wk)=1,此时接着对于转化后的表达式进行一阶泰勒展开,分离Ak和Bk;对于Wk=wkwk H的非凸秩一约束,引入‖Wk*-‖Wk2=0将其转化为目标函数中的惩罚项,其中
Figure BDA00041043355900000311
和‖Wk2=ε1(Wk)分别表示Wk的核范数和谱范数,接着利用一阶泰勒展开得到惩罚项的凸上界:
Figure BDA00041043355900000312
其中
Figure BDA00041043355900000313
Figure BDA00041043355900000314
是第τ1次迭代中
Figure BDA00041043355900000315
的最大特征值对应的特征向量,将其引入目标函数作为惩罚函数,得到一个关于wk子问题的标准的半定规划子问题。
进一步地,关于wk子问题的半定规划子问题的求解过程为:首先初始化惩罚因子η,η=μη,μ<1,然后随着逐渐减小η,以获得一个总体次优解,当惩罚项满足以下标准时,该过程终止:
Figure BDA00041043355900000316
其中∈1为预定义的最大违反量。
进一步地,对于MF-RIS的系数矩阵Θk子问题,步骤S2半定规划子问题转化方法为:定义vk和Vk,其中vk=[uk;1],
Figure BDA0004104335590000041
rank(Vk)=1,定义gk和Gk,gk=[gk,1,gk,2,...,gk,M]H,Gk=Hwk,Qk=diag([|gk,1|2,|gk,2|2,...,|gk,M|2]),
Figure BDA0004104335590000042
得到
Figure BDA0004104335590000043
Figure BDA0004104335590000044
定义
Figure BDA0004104335590000045
Figure BDA0004104335590000046
代替关于P0的约束,对于秩一约束,基于MF-RIS系数表达式,采用以下形式进行替换:
Figure BDA0004104335590000047
其中‖Vk*和‖Vk2分别表示矩阵Vk的核范数和谱范数,
Figure BDA0004104335590000048
并且
Figure BDA0004104335590000049
是对应于第τ2次迭代中
Figure BDA00041043355900000410
的最大特征值的特征向量,将其引入目标函数作为惩罚函数,得到一个关于MF-RIS的系数矩阵子问题的标准的半定规划子问题。
进一步地,关于Θk子问题的半定规划子问题的求解过程为:首先初始化惩罚因子ξ,ξ>0是确保Vk秩为1的惩罚因子,然后随着逐渐减小ξ,以获得一个总体次优解,当惩罚项满足以下标准时,该过程终止:
Figure BDA00041043355900000411
其中∈2为预定义的最大违反量。
进一步地,步骤S3中MF-RIS系数和基站有源波束赋形子问题交替优化的方法为:初始化
Figure BDA00041043355900000412
误差容限Δ,最大迭代次数T0,max,惩罚因子η和ξ,以及预定义的阈值∈;设置迭代指数τ0=0,给定
Figure BDA00041043355900000413
通过求解基站有源波束赋形子问题更新
Figure BDA00041043355900000414
的值,接着给定
Figure BDA00041043355900000415
通过求解MF-RIS系数子问题更新
Figure BDA00041043355900000416
的值,更新τ0=τ0+1,重复上述过程直到
Figure BDA00041043355900000417
或是τ0>T0,max,接着利用
Figure BDA00041043355900000418
更新
Figure BDA00041043355900000419
直到MF-RIS系数和有源波束赋形向量的约束满足预定义的阈值∈,则输出收敛解
Figure BDA00041043355900000420
与现有技术相比,本发明的有益效果为:
本发明提出的可重构的多功能智能超表面,不仅可以基于场等效原理将入射信号分为折射和反射两部分,还可以借助有源负载放大输出信号,主动地改变无线信号传播环境,从而提升无线通信系统的信号传输性能,实现全空间覆盖并克服双重衰落问题。此外,针对非正交多址接入(Non-Orthgonal Multiple Access,NOMA)网络中的资源分配问题,将本发明的MF-RIS应用于NOMA网络,以最大化用户和速率为目标,构建了一个联合设计基站波束赋形与MF-RIS系数的非凸优化问题;针对所构建的优化问题,首先采用惩罚函数法和连续凸逼近法(Successive Convex Approximation,SCA)将原始优化问题转化为多个半定规划子问题,然后使用交替优化方法获取基站波束赋形向量和MF-RIS系数的次优解,可以调整用户之间的信道条件,从而扩大不同用户信道增益之间的差异,最终提高所有NOMA用户的可达和速率。仿真结果表明,所提出的MF-RIS辅助NOMA网络能够提供比传统RIS高约57%的和速率,并且MF-RIS更倾向于部署在用户侧以获得更好的性能。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明中记载的一些实施例,对于本领域普通技术人员来讲,还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的MF-RIS辅助NOMA网络原理示意图。
图2为本发明实施例提供的MF-RIS应用于模拟场景中的仿真示意图。
图3为本发明实施例提供的MF-RIS在模拟场景中仿真得到的功率结果示意图。
具体实施方式
为了更好地理解本技术方案,下面结合附图对本发明的方法做详细的说明。
请参阅图1,本发明首先提出了一种可重构的多功能智能超表面(MF-RIS),基于场等效原理,利用表面等效电阻抗和磁阻抗实现入射信号的反射和折射,利用负电阻元件实现信号放大,从而实现全空间覆盖并克服双重衰落问题。
然后,我们研究了下行链路多用户NOMA网络中的资源分配问题,其中MF-RIS被部署来协助基站和用户之间的通信,扩大不同用户信道增益之间的差异,从而提高所有NOMA用户的可达和速率。具体来说,通过利用SCA和惩罚函数的方法,联合优化有源波束赋形向量和MF-RIS系数以最大化可实现的和速率。
MF-RIS应用在下行链路多用户NOMA网络中的资源分配问题的步骤如下:
S1设计包含基站有源波束赋形向量和MF-RIS的系数,以最大化所有用户的可实现和速率为目标,构建优化问题和约束条件;
我们考虑MF-RIS辅助的NOMA下行链路网络,其中部署由M个元件组成的MF-RIS来辅助从N天线BS到K个单天线用户的通信,请参阅图1的右侧。元件和用户的集合分别由
Figure BDA0004104335590000061
Figure BDA0004104335590000062
表示。基站-用户、基站-RIS和RIS-用户的信道分别由
Figure BDA0004104335590000063
表示。此外,我们定义
Figure BDA0004104335590000064
作为折射(p=t)或反射(p=r)波束赋形向量,其中p∈{t,r}表示到达特定用户的信号被MF-RIS折射或反射,
Figure BDA0004104335590000065
Figure BDA0004104335590000066
分别表示第m个元件的幅度和相移响应,最大放大系数满足βmax≥1。由于能量守恒定律,放大功能受到提供给放大器的能量的限制,即,
Figure BDA0004104335590000067
如果用户k位于反射空间,则用户k的MF-RIS的对角矩阵由Θk=diag(ur)给出,否则Θk=diag(ut)。
我们假设所有信道的完美信道状态信息在基站处是可获取的。那么在用户k处接收的信号被表示为:
Figure BDA0004104335590000068
其中x=∑kwksk,x表示发射信号,wk
Figure BDA0004104335590000069
分别表示用户k的发射预编码波束和信息符号。
Figure BDA00041043355900000610
表示MF-RIS处的噪声,每个元件的噪声功率为
Figure BDA00041043355900000611
Figure BDA00041043355900000612
表示用户k处的加性高斯白噪声,功率为
Figure BDA00041043355900000613
通过使用SIC,信道质量高的用户可以消除来自信道质量较差用户的干扰,以提高信干噪比。因此,我们假设用户的索引相对于他们的信道增益以递增的顺序排列,即:
‖h12≤‖h22≤…≤‖hK2
其中,
Figure BDA0004104335590000071
为等效组合信道。对于固定的SIC解码顺序,用户k对应的可达和速率由Rk=log2(1+γk)给出,其中γk可由下式获得:
Figure BDA0004104335590000072
我们的目标是通过联合优化基站的有源波束赋形向量和MF-RIS的系数来最大化所有用户的和速率。在BS发射功率和MF-RIS放大功率以及用户的最小服务质量(Qualityof Service,QoS)约束下,所考虑的优化问题可以表述为:
Figure BDA0004104335590000073
Figure BDA0004104335590000074
Figure BDA0004104335590000075
Figure BDA0004104335590000076
Figure BDA0004104335590000077
其中Pmax和P0分别表示基站和MF-RIS的最大发射和放大功率。
Figure BDA0004104335590000078
表示用户k的最小速率要求。具体地,在第2-5个公式中分别给出了发射功率、放大功率、QoS要求和解码顺序的约束。
可以观察到,由于非凸的目标函数和约束,优化问题是难以处理的。此外,有源波束赋形向量和MF-RIS系数高度耦合,难以直接求解。因此,我们的目标是将优化问题转化为一些易处理的凸子问题,并通过迭代进行交替优化。
S2、根据步骤S1构建的约束条件确定初始用户和速率表达式,对于其中高度耦合的MF-RIS系数和基站有源波束功率约束引入辅助变量集进行转化,应用一阶泰勒展开获得转化后的下界;
给定MF-RIS系数,有源波束赋形优化问题仍然是非凸的。为了便于解决,我们首先引入一个辅助变量集
Figure BDA0004104335590000079
其中Ak和Bk定义为:
Figure BDA00041043355900000710
Figure BDA00041043355900000711
因此,数据速率可以改写为:
Rk=log2(1+(AkBk)-1)
那么,前文优化问题中的有源波束赋形优化问题可以等价地表示为:
Figure BDA0004104335590000081
Figure BDA0004104335590000082
Figure BDA0004104335590000083
Figure BDA0004104335590000084
Figure BDA0004104335590000085
我们进一步定义:
Figure BDA0004104335590000086
Dk=(HHΘk)(HHΘk)H
Figure BDA0004104335590000087
Wk≥0,rank(Wk)=1
为了处理非凸约束,我们采用一阶泰勒展开,然后我们得到如下下界:
Figure BDA0004104335590000088
其中
Figure BDA0004104335590000089
Figure BDA00041043355900000810
分别是第τ1次迭代中Ak和Bk的可行点。
S3、以步骤S2转换后的子问题为基础引入惩罚函数。将该子问题转化为半正定规划子问题。
对于S4中Wk的非凸秩一约束,我们假设将其转化为目标函数中的惩罚项。因此,我们首先引入一个等式:
Figure BDA00041043355900000811
其中
Figure BDA00041043355900000812
和‖Wk2=ε1(Wk)分别表示Wk的核范数和谱范数。εi(Wk)是矩阵Wk的第i个最大奇异值。因此,当矩阵Wk秩为1时,等式成立。
接下来,我们将该式添加到目标函数上作为惩罚项。由于该等式作为惩罚项使得目标函数非凸,我们应用一阶泰勒展开来获得凸上界,如下:
Figure BDA0004104335590000091
其中
Figure BDA0004104335590000092
Figure BDA0004104335590000093
是第τ1次迭代中
Figure BDA0004104335590000094
的最大特征值对应的特征向量。
通过将凸上界引入目标函数,我们得到以下问题:
Figure BDA0004104335590000095
Figure BDA0004104335590000096
Figure BDA0004104335590000097
Figure BDA0004104335590000098
Figure BDA0004104335590000099
Figure BDA00041043355900000910
Figure BDA00041043355900000911
Figure BDA00041043355900000912
其中η>0是惩罚因子,如果Wk秩不为1,则惩罚目标函数。可以验证,当η→0时,上述问题的解{Wk}总是满足等式
Figure BDA00041043355900000913
转化后的问题是一个标准的SDP问题,它可以通过CVX有效地解决。为了获得高质量的解,我们首先初始化一个大的η,以找到一个可行的初始点,然后随着η=μη,μ<1逐渐减小η,直到一个足够小的值,以获得一个总体次优解。当惩罚项满足以下标准时,该过程终止:
Figure BDA00041043355900000914
其中∈1表示等式
Figure BDA00041043355900000915
预定义的最大违反值。
对于MF-RIS的系数设计,我们定义如下:
vk=[uk;1]
Figure BDA00041043355900000916
rank(Vk)=1
设gk=[gk,1,gk,2,...,gk,M]H,Gk=Hwk
则有Qk=diag([|gk,1|2,|gk,2|2,...,|gk,M|2])
Figure BDA00041043355900000917
给定
Figure BDA0004104335590000101
可以得到
Figure BDA0004104335590000102
Figure BDA0004104335590000103
因此,根据上述转化,非凸约束可以转化为凸约束。
为了处理非凸约束‖h12≤‖h22≤...≤‖hK2,和
Figure BDA0004104335590000104
我们定义
Figure BDA0004104335590000105
则有:
Figure BDA0004104335590000106
Figure BDA0004104335590000107
根据上述变换,我们可以得到:
Figure BDA0004104335590000108
Figure BDA0004104335590000109
基于上式,可以将‖h12≤‖h22≤...≤‖hK2中的解码顺序重新表述为:
Figure BDA00041043355900001010
然后,给定有源波束赋形向量,MF-RIS系数设计的子问题可以由下式给出:
Figure BDA00041043355900001011
Figure BDA00041043355900001012
Figure BDA00041043355900001013
Figure BDA00041043355900001014
Figure BDA00041043355900001015
Figure BDA00041043355900001016
Figure BDA00041043355900001017
Figure BDA00041043355900001018
Figure BDA00041043355900001019
Figure BDA0004104335590000111
其中Fi表示Fk中的wk被wi代替时的值。
类似于基站有源波束功率约束,我们用以下形式替换上文中的秩一约束:
Figure BDA0004104335590000112
其中‖Vk*和‖Vk2分别表示矩阵Vk的核范数和谱范数。且
Figure BDA0004104335590000113
并且
Figure BDA0004104335590000114
是对应于第τ2次迭代中
Figure BDA0004104335590000115
的最大特征值的特征向量。
通过将
Figure BDA0004104335590000116
引入
Figure BDA0004104335590000117
本发明的上述问题可以被重新表述为:
Figure BDA0004104335590000118
Figure BDA0004104335590000119
Figure BDA00041043355900001110
Figure BDA00041043355900001111
Figure BDA00041043355900001112
Figure BDA00041043355900001113
Figure BDA00041043355900001114
Figure BDA00041043355900001115
Figure BDA00041043355900001116
其中ξ>0是确保Vk秩为1的惩罚因子。
上述变换后的优化问题是标准的SDP问题。CVX可以有效地解决这个问题。迭代结束标准由下式给出:
Figure BDA00041043355900001117
其中∈2表示预定义的最大违反量。
基于上述推导,我们提出了一种基于惩罚的迭代算法来有效地解决问题。细节在下述算法中给出。
具体的,在算法中基于惩罚函数法和逐次凸逼近法对上述MF-RIS系数和基站有源波束赋形的子问题进行交替优化,初始化
Figure BDA0004104335590000121
误差容限Δ,最大迭代次数T0,max,惩罚因子η和ξ,以及预定义的阈值∈。设置迭代次数τ0=0,给定
Figure BDA0004104335590000122
通过求解基站有源波束赋形的子问题更新
Figure BDA0004104335590000123
的值,接着给定
Figure BDA0004104335590000124
通过求解MF-RIS系数的子问题更新
Figure BDA0004104335590000125
的值,更新τ0=τ0+1,重复上述过程直到
Figure BDA0004104335590000126
或是τ0>T0,max,接着更新
Figure BDA0004104335590000127
Figure BDA0004104335590000128
直到MF-RIS系数和有源波束赋形向量的约束满足预定义的阈值∈,则输出收敛解
Figure BDA0004104335590000129
由于两个优化问题的目标在迭代中是非递减的,并且系统吞吐量有上限,所提出的算法可以保证收敛。此外,算法的复杂度为
Figure BDA00041043355900001210
其中Iin和Iout分别表示收敛所需的内部和外部迭代的次数。
请参阅图2,在运用场景中,基站和MF-RIS分别位于(0,0,0)米和(0,50,20)米处。此外,用户被平均分成两部分,分别分布在以(0,45,0)米和(0,55,0)米为圆心、半径r=3米的圆上。我们假设所有信道为莱斯衰落,并设置K=6,N=16,M=100,Pmax=20dBm。接着我们设置三种对比方案,分别部署Reflecting-only RIS,Active RIS(一种能够同时支持信号反射和信号放大的RIS元件)以及STAR-RIS来协助下行NOMA系统通信。
由此进行仿真的结果请参阅图3,具体来说,当Pmax=10dBm时,MF-RIS辅助的NOMA方案提供比仅反射RIS的方案高57%的和速率。这个结果可以解释如下。一方面,信号反射和折射功能使得MF-RIS能够服务于整个空间中的所有用户,从而获得更高的系统和速率。另一方面,通过提供额外的功率来放大入射信号,MF-RIS能够有效地提升级联信道质量。此外,MF-RIS方案优于STAR-RIS和active RIS方案的原因是,后两者仅能部分解决仅反射RIS所面临的问题(即,半空间覆盖和双重衰落问题)。
数值结果验证了MF-RIS的有效性和MF-RIS相对于传统RIS的优越性。仿真结果表明,所提出的MF-RIS辅助NOMA网络能够提供比传统RIS高约57%的和速率,并且MF-RIS更倾向于部署在用户侧以获得更好的性能。
以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换,但这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。

Claims (9)

1.一种可重构的多功能智能超表面,其特征在于,其基于场等效原理,利用表面等效电阻抗和磁阻抗实现入射信号的反射和折射,利用负电阻元件实现信号放大。
2.根据权利要求1所述的可重构的多功能智能超表面在NOMA网络资源分配中的应用,其特征在于,包括以下步骤:
S1、设计包含基站有源波束赋形向量和MF-RIS的系数,以最大化所有用户的可实现和速率为目标,构建优化问题和约束条件;
S2、根据步骤S1构建的约束条件确定初始用户和速率表达式,对于其中高度耦合的MF-RIS系数和基站有源波束赋形向量,利用连续凸逼近和惩罚函数法转换为半定规划子问题;
S3、对步骤S2转换后的半定规划子问题进行交替优化,当和速率变化量满足收敛条件,得到有源波束赋形向量和MF-RIS的系数,从而实现最优性能。
3.根据权利要求2所述的可重构的多功能智能超表面在NOMA网络资源分配中的应用,其特征在于,步骤S1优化问题为:
Figure FDA0004104335570000011
约束条件为:
Figure FDA0004104335570000012
Figure FDA0004104335570000013
Figure FDA0004104335570000014
Figure FDA0004104335570000015
其中M表示MF-RIS的元件个数,N表示基站天线的个数,K表示用户个数,
Figure FDA0004104335570000016
分别表示基站-用户、基站-RIS和RIS-用户的信道,
Figure FDA0004104335570000017
Figure FDA0004104335570000018
分别表示第m个元件的幅度和相移响应,βmax表示最大放大系数,Pmax表示基站的最大发射功率,Po表示MF-RIS的最大放大功率,Rk表示用户k的数据速率,
Figure FDA0004104335570000019
表示用户k的最小速率要求,wk表示用户k的波束向量,Θk表示用户k的MF-RIS的系数矩阵,
Figure FDA00041043355700000110
表示MF-RIS处引入的加性高斯白噪声,hk为基站-MF-RIS-用户的等效组合信道。
4.根据权利要求3所述的可重构的多功能智能超表面在NOMA网络资源分配中的应用,其特征在于,用户k的波束向量wk和用户k的MF-RIS的系数矩阵Θk在优化问题中是高度耦合MF-RIS系数,对wk和Θk分别优化。
5.根据权利要求4所述的可重构的多功能智能超表面在NOMA网络资源分配中的应用,其特征在于,对于wk子问题中的目标函数,步骤S2半定规划子问题转化方法为:利用辅助变量Ak和Bk代替wk和hk,等价的表示有源波束赋形问题的表达式,引入
Figure FDA0004104335570000021
Dk=(HHΘk)(HHΘk)H,Wk=wkwk H,对表达式进行进一步的表示,其中Wk≥0且存在非凸秩一约束rank(Wk)=1,此时接着对于转化后的表达式进行一阶泰勒展开,分离Ak和Bk;对于Wk=wkwk H的非凸秩一约束,引入‖Wk*-‖Wk2=0将其转化为目标函数中的惩罚项,其中
Figure FDA0004104335570000022
和‖Wk2=ε1(Wk)分别表示Wk的核范数和谱范数,接着利用一阶泰勒展开得到惩罚项的凸上界:
Figure FDA0004104335570000023
其中
Figure FDA0004104335570000024
Figure FDA0004104335570000025
是第τ1次迭代中
Figure FDA0004104335570000026
的最大特征值对应的特征向量,将其引入目标函数作为惩罚函数,得到一个关于wk子问题的标准的半定规划子问题。
6.根据权利要求5所述的可重构的多功能智能超表面在NOMA网络资源分配中的应用,其特征在于,关于wk子问题的半定规划子问题的求解过程为:首先初始化惩罚因子η,η=μη,μ<1,然后随着逐渐减小η,以获得一个总体次优解,当惩罚项满足以下标准时,该过程终止:
Figure FDA0004104335570000027
其中∈1为预定义的最大违反量。
7.根据权利要求4所述的可重构的多功能智能超表面在NOMA网络资源分配中的应用,其特征在于,对于MF-RIS的系数矩阵Θk子问题,步骤S2半定规划子问题转化方法为:
定义vk和Vk,其中vk=[uk;1],
Figure FDA0004104335570000028
Vk≥0,rank(Vk)=1,定义gk和Gk,gk=[gk,1,gk,2,...,gk,M]H,Gk=Hwk,Qk=diag([|gk,1|2,|gk,2|2,...,|gk,M|2]),
Figure FDA0004104335570000029
得到
Figure FDA00041043355700000210
Figure FDA00041043355700000211
定义
Figure FDA00041043355700000212
Figure FDA0004104335570000031
Figure FDA0004104335570000032
代替关于P0的约束,对于秩一约束,基于MF-RIS系数表达式,采用以下形式进行替换:
Figure FDA0004104335570000033
其中‖Vk*和‖Vk2分别表示矩阵Vk的核范数和谱范数,
Figure FDA0004104335570000034
并且
Figure FDA0004104335570000035
是对应于第τ2次迭代中
Figure FDA0004104335570000036
的最大特征值的特征向量,将其引入目标函数作为惩罚函数,得到一个关于MF-RIS的系数矩阵子问题的标准的半定规划子问题。
8.根据权利要求7所述的可重构的多功能智能超表面在NOMA网络资源分配中的应用,其特征在于,关于Θk子问题的半定规划子问题的求解过程为:首先初始化惩罚因子ξ,ξ>0是确保Vk秩为1的惩罚因子,然后随着逐渐减小ξ,以获得一个总体次优解,当惩罚项满足以下标准时,该过程终止:
Figure FDA0004104335570000037
其中∈2为预定义的最大违反量。
9.根据权利要求4所述的可重构的多功能智能超表面在NOMA网络资源分配中的应用,其特征在于,步骤S3中MF-RIS系数和基站有源波束赋形子问题交替优化的方法为:初始化
Figure FDA0004104335570000038
误差容限Δ,最大迭代次数T0,max,惩罚因子η和ξ,以及预定义的阈值∈;设置迭代指数τ0=0,给定
Figure FDA0004104335570000039
通过求解基站有源波束赋形子问题更新
Figure FDA00041043355700000310
的值,接着给定
Figure FDA00041043355700000311
通过求解MF-RIS系数子问题更新
Figure FDA00041043355700000312
的值,更新τ0=τ0+1,重复上述过程直到
Figure FDA00041043355700000313
或是τ0>T0,max,接着利用
Figure FDA00041043355700000314
更新
Figure FDA00041043355700000315
直到MF-RIS系数和有源波束赋形向量的约束满足预定义的阈值∈,则输出收敛解
Figure FDA00041043355700000316
CN202310187291.9A 2023-02-22 2023-02-22 一种可重构的多功能智能超表面及其应用 Pending CN116170102A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310187291.9A CN116170102A (zh) 2023-02-22 2023-02-22 一种可重构的多功能智能超表面及其应用

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310187291.9A CN116170102A (zh) 2023-02-22 2023-02-22 一种可重构的多功能智能超表面及其应用

Publications (1)

Publication Number Publication Date
CN116170102A true CN116170102A (zh) 2023-05-26

Family

ID=86421775

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310187291.9A Pending CN116170102A (zh) 2023-02-22 2023-02-22 一种可重构的多功能智能超表面及其应用

Country Status (1)

Country Link
CN (1) CN116170102A (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116600311A (zh) * 2023-07-17 2023-08-15 南京邮电大学 一种双ris辅助上行noma系统的和速率最大化方法及系统
CN117439673A (zh) * 2023-12-21 2024-01-23 厦门大学 基于非正交多址的智能超表面辅助水声数据传输方法
CN117499961A (zh) * 2023-12-29 2024-02-02 南京邮电大学 一种双ris辅助noma上行传输能效优化方法及系统
CN117997402A (zh) * 2024-01-19 2024-05-07 安徽大学 多层智能超表面辅助的多用户上行通信方法和装置

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116600311A (zh) * 2023-07-17 2023-08-15 南京邮电大学 一种双ris辅助上行noma系统的和速率最大化方法及系统
CN116600311B (zh) * 2023-07-17 2023-11-10 南京邮电大学 一种双ris辅助上行noma系统的和速率最大化方法及系统
CN117439673A (zh) * 2023-12-21 2024-01-23 厦门大学 基于非正交多址的智能超表面辅助水声数据传输方法
CN117439673B (zh) * 2023-12-21 2024-03-12 厦门大学 基于非正交多址的智能超表面辅助水声数据传输方法
CN117499961A (zh) * 2023-12-29 2024-02-02 南京邮电大学 一种双ris辅助noma上行传输能效优化方法及系统
CN117499961B (zh) * 2023-12-29 2024-03-15 南京邮电大学 一种双ris辅助noma上行传输能效优化方法及系统
CN117997402A (zh) * 2024-01-19 2024-05-07 安徽大学 多层智能超表面辅助的多用户上行通信方法和装置
CN117997402B (zh) * 2024-01-19 2024-08-06 安徽大学 多层智能超表面辅助的多用户上行通信方法和装置

Similar Documents

Publication Publication Date Title
CN116170102A (zh) 一种可重构的多功能智能超表面及其应用
CN112672375B (zh) 一种智能反射面辅助的非正交多址接入网络中安全通信方法
CN113193894B (zh) 可重构智能表面辅助的多用户miso下行无线通信谱效能效联合优化方法
CN106253956B (zh) 基于码本的模数混合预编码方法
CN114828258B (zh) 一种智能反射面辅助认知无线电系统资源分配优化方法
CN114221724B (zh) 智能反射面辅助的认知无线携能通信网络性能优化方法
CN109714091B (zh) 一种在毫米波mimo系统中基于分层设计的迭代混合预编码方法
Ren et al. Ris-assisted cooperative NOMA with SWIPT
CN112822703B (zh) 智能反射面辅助的非正交多址接入系统性能增益优化方法
CN114070365B (zh) 一种智能反射表面辅助的低射频复杂度多用户mimo上行频谱效率优化方法
CN105743559B (zh) 一种Massive MIMO混合波束形成和空时编码多用户下行传输方法
CN108809383B (zh) 一种用于massive MIMO上行系统信号的联合检测方法
Zhang et al. RIS-assisted self-interference mitigation for in-band full-duplex transceivers
CN113852402B (zh) 一种irs辅助的noma-mimo大容量接入方法
CN110191476B (zh) 一种基于可重构天线阵列的非正交多址接入方法
Han et al. Secure transmission for STAR-RIS aided NOMA against internal eavesdropping
CN114785388A (zh) 智能全向面辅助的多用户大规模simo上行m阶调制加权和速率优化方法
CN114765785B (zh) 一种基于最大信噪比的多智能反射面选择方法
Huang et al. Joint microstrip selection and beamforming design for mmwave systems with dynamic metasurface antennas
Jalali et al. Power-efficient joint resource allocation and decoding error probability for multiuser downlink MISO with finite block length codes
Guo et al. Hybrid MRT and ZF Learning for Energy-Efficient Transmission in Multi-RIS-Assisted Networks
CN116033461B (zh) 一种基于star-ris辅助的共生无线电传输方法
CN116208971B (zh) 有源ris辅助的非正交多址接入系统上行传输方法
CN115733533B (zh) 基于流形的irs辅助去蜂窝大规模mimo系统波束赋形方法
CN111740766A (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