CN114640379A - 一种基于智能反射面阵元分组的波束优化方法及*** - Google Patents

一种基于智能反射面阵元分组的波束优化方法及*** Download PDF

Info

Publication number
CN114640379A
CN114640379A CN202210143464.2A CN202210143464A CN114640379A CN 114640379 A CN114640379 A CN 114640379A CN 202210143464 A CN202210143464 A CN 202210143464A CN 114640379 A CN114640379 A CN 114640379A
Authority
CN
China
Prior art keywords
optimization problem
irs
optimization
array element
reflection
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
CN202210143464.2A
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.)
University of Science and Technology Beijing USTB
Shunde Graduate School of USTB
Original Assignee
University of Science and Technology Beijing USTB
Shunde Graduate School of USTB
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 University of Science and Technology Beijing USTB, Shunde Graduate School of USTB filed Critical University of Science and Technology Beijing USTB
Priority to CN202210143464.2A priority Critical patent/CN114640379A/zh
Publication of CN114640379A publication Critical patent/CN114640379A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • 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
    • 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/04013Intelligent reflective surfaces
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B17/00Monitoring; Testing
    • H04B17/30Monitoring; Testing of propagation channels
    • H04B17/391Modelling the propagation channel

Landscapes

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

Abstract

本发明公开了一种基于智能反射面阵元分组的波束优化方法及***,涉及无线通信技术领域。包括:构建基站与多个用户间通过智能反射面IRS进行信号传输的初始优化问题;采用克罗内克积方法对智能反射面IRS的反射阵元进行分组;将初始优化问题转化为分组后的IRS反射阵元的优化问题P1;利用分式规划理论将优化问题P1转化为复杂度降低的优化问题P2;引入辅助变量,并将优化问题P2分解为三个解耦的子问题,使用交替迭代优化方法求出每个子问题闭合形式的最优解;根据最优解,得到联合优化后的基站BS波束成形矢量和IRS被动波束成形矢量,最终得到智能反射面阵元分组的波束优化方法。本发明能够达到降低IRS被动波束成形矢量计算复杂度的目的。

Description

一种基于智能反射面阵元分组的波束优化方法及***
技术领域
本发明涉及无线通信技术领域,特别是指一种基于智能反射面阵元分组的波束优化方法及***。
背景技术
IRS(Intelligent Reflecting Surface,智能反射面)是一项通过软件控制反射重新配置无线传播环境的新技术,能够显著提高无线通信网络的性能。IRS是一个由大量低成本被动反射阵元组成的平面,每个阵元都能独立调整入射信号的振幅和相位,进而可以调控BS(Base Station,基站)和用户间信号传输方向,达到主动控制无线信道的目的。
然而,IRS通常包含大量反射阵元,随着IRS反射阵元数量的增加,会给通信***带来许多技术挑战。最重要的是以下两点。一是大量的IRS反射阵元将大大增加信道估计时间,这将给通信***带来很大时延,降低用户的通信质量。另一点是,随着IRS反射阵元数量增加,计算IRS被动波束成形矢量的复杂性将增加。对IRS相邻反射阵元进行分组,并根据相关优化算法得出分组后的IRS被动波束成形矢量,在保证用户通信质量的情况下,可以大大降低IRS被动波束成形矢量的计算复杂度。
IRS作为一种新的无线通信辅助传输技术,在无线通信的覆盖扩展领域已有一定的研究。现有研究主要采用基于交替迭代优化算法来优化BS端波束成形矢量和IRS被动波束成形矢量。
现有研究表明,IRS辅助下的无线通信***比传统的通信***获得了更高的用户速率,特别是在用户和基站间存在障碍物造成严重阻塞时,提升效果更为明显。但在实际应用中,IRS需要根据CSI(Channel State Information,信道状态信息)的变化实时的调整反射阵元的反射系数,由于IRS通常具有多个反射阵元,这将给被动波束成形矢量的优化带来很高的计算复杂度,造成一定的时延。
发明内容
本发明针对现有技术在IRS辅助的多用户通信场景下,IRS被动波束成形矢量计算复杂度大的问题,提出了本发明。
为解决上述技术问题,本发明提供如下技术方案:
一方面,本发明提供了一种基于智能反射面阵元分组的波束优化方法,该方法由基于智能反射面阵元分组的波束优化***实现,该***包括基站、智能反射面RIS以及多个用户,该方法包括:
S1、构建基站与多个用户间通过智能反射面IRS进行信号传输的初始优化问题;初始优化问题为优化多个用户的总速率。
S2、采用克罗内克积方法对智能反射面IRS的反射阵元进行分组;将初始优化问题转化为分组后的IRS反射阵元的优化问题P1。
S3、利用分式规划理论将优化问题P1转化为复杂度降低的优化问题P2。
S4、引入辅助变量,并将优化问题P2分解为三个解耦的子问题,使用交替迭代优化方法求出每个子问题闭合形式的最优解。
S5、根据最优解,得到联合优化后的基站BS波束成形矢量和IRS被动波束成形矢量,最终得到智能反射面阵元分组的波束优化方法。
可选地,初始优化问题中,多个用户中的第k个用户的接收信号,如下式(1)所示:
Figure BDA0003507654130000021
其中,
Figure BDA0003507654130000022
为IRS和第k个用户之间的信道,M是一个IRS具有的反射阵元个数,hM,k是第M个反射阵元和第k个用户之间的信道,T为矩阵的转置,H为矩阵的转置共轭;Φ=ηdiag([θ1,...,θM]T)为IRS反射系数表示的矩阵,η∈[0,1]表示IRS中第M个反射阵元的开或关状态,
Figure BDA0003507654130000023
Figure BDA0003507654130000024
表示第M个阵元的反射系数;
Figure BDA0003507654130000025
是BS和IRS之间的信道,N是基站配备有天线根数,
Figure BDA0003507654130000026
是BS和第M个反射阵元之间的信道,设定所有信道都是Saleh-Valenzuela毫米波信道;P=[P1,...PK]表示波束成形矩阵,其中
Figure BDA0003507654130000031
表示发送给第k个用户信号的波束成形矢量,K是用户个数;sk表示BS发送给第k个用户的信号,满足
Figure BDA0003507654130000032
k∈{1,2,...,K};
Figure BDA0003507654130000033
表示发送给第j个用户信号的波束成形矢量;sj表示BS发送给第j个用户的信号,满足
Figure BDA0003507654130000034
j∈{1,2,...,K};μk表示均值为0和方差为
Figure BDA0003507654130000035
的加性高斯白噪声。
可选地,分组后的IRS反射阵元的优化问题中,多个用户中的第k个用户的接收信号,如下式(2)所示:
Figure BDA0003507654130000036
其中,
Figure BDA0003507654130000037
为IRS分组后的相移因子表示成的矩阵,
Figure BDA0003507654130000038
表示第J组相移因子,J为IRS阵元分组个数;
Figure BDA0003507654130000039
为IRS与第k个用户之间的信道;
Figure BDA00035076541300000310
为BS与IRS之间的信道。
可选地,分组后的IRS反射阵元的优化问题中,多个用户中的第k个用户对应的信干噪比,如下式(3)所示:
Figure BDA00035076541300000311
可选地,分组后的IRS反射阵元的优化问题P1,如下式(4)-(6)所示:
Figure BDA00035076541300000312
s.t.|φj|=1,n=1,2,...,J (5)
tr(PPH)≤Pmax (6)
其中,|φj|=1表示将IRS每一组反射阵元的反射幅度设为1;Pmax表示BS最大发射功率。
可选地,复杂度降低的优化问题P2,如下式(7)-(10)所示:
Figure BDA00035076541300000313
s.t.|φj|=1,n=1,2,...,J (8)
tr(PPH)≤Pmax (9)
αk≥0,k=1,2,...,K (10)
其中,α=[α1,...,αK]T为引入的辅助变量。
可选地,αk的最优解的求解方法包括:
令变量P和
Figure BDA0003507654130000041
固定,对上述公式(6)中的变量αk求偏导,得到αk的最优解;
将αk的最优解代入上式(6)中,将优化问题P2转换为优化问题P3。
可选地,P的最优解的求解方法包括:
Figure BDA0003507654130000042
为等效信道,将优化问题P3重新表示为优化问题P4。
给定αk
Figure BDA0003507654130000043
引入辅助变量β=[β1,...,βK]T,对优化问题P4进行二次转换,得到优化问题P5。
给定P,通过对优化问题P5求偏导,得到βk的最优解。
固定βk,通过拉格朗日乘子法得到pk的最优解。
可选地,
Figure BDA0003507654130000044
的最优解的求解方法包括:
给定α和P,引入辅助变量ρ=[ρ1,...,ρK]T,对优化问题P5进行二次转换,得到优化问题P6。
给定ρ,将优化问题P6重新表示为优化问题P7。
通过约束松弛法将优化问题P7重新表示为优化问题P8。
根据拉格朗日对偶分解法,将优化问题P8重新表示为优化问题P9。
对优化问题P9求解,最终得到
Figure BDA0003507654130000045
的最优解。
另一方面,本发明提供了一种基于智能反射面阵元分组的波束优化***,该***应用于实现基于智能反射面阵元分组的波束优化方法,该***包括基站、智能反射面RIS以及多个用户;其中:
基站,用于向多个用户传输信号。
智能反射面RIS,用于采用克罗内克积方法对智能反射面IRS的反射阵元进行分组;将初始优化问题转化为分组后的IRS反射阵元的优化问题P1;利用分式规划理论将优化问题P1转化为复杂度降低的优化问题P2;引入辅助变量,并将优化问题P2分解为三个解耦的子问题,使用交替迭代优化方法求出每个子问题闭合形式的最优解;根据最优解,得到联合优化后的BS波束成形矢量和IRS被动波束成形矢量,最终得到智能反射面阵元分组的波束优化方法。
多个用户,用于接收基站传输的信号。
可选地,初始优化问题中,多个用户中的第k个用户的接收信号,如下式(1)所示:
Figure BDA0003507654130000051
其中,
Figure BDA0003507654130000052
为IRS和第k个用户之间的信道,M是一个IRS具有的反射阵元个数,hM,k是第M个反射阵元和第k个用户之间的信道,T为矩阵的转置,H为矩阵的转置共轭;Φ=ηdiag([θ1,...,θM]T)为IRS反射系数表示的矩阵,η∈[0,1]表示IRS中第M个反射阵元的开或关状态,
Figure BDA0003507654130000053
Figure BDA0003507654130000054
表示第M个阵元的反射系数;
Figure BDA0003507654130000055
是BS和IRS之间的信道,N是基站配备有天线根数,
Figure BDA0003507654130000056
是BS和第M个反射阵元之间的信道,设定所有信道都是Saleh-Valenzuela毫米波信道;P=[P1,...PK]表示波束成形矩阵,其中
Figure BDA0003507654130000057
表示发送给第k个用户信号的波束成形矢量,K是用户个数;sk表示BS发送给第k个用户的信号,满足
Figure BDA0003507654130000058
k∈{1,2,...,K};
Figure BDA0003507654130000059
表示发送给第j个用户信号的波束成形矢量;sj表示BS发送给第j个用户的信号,满足
Figure BDA00035076541300000510
j∈{1,2,...,K};μk表示均值为0和方差为
Figure BDA00035076541300000511
的加性高斯白噪声。
可选地,分组后的IRS反射阵元的优化问题中,多个用户中的第k个用户的接收信号,如下式(2)所示:
Figure BDA00035076541300000512
其中,
Figure BDA00035076541300000513
为IRS分组后的相移因子表示成的矩阵,
Figure BDA00035076541300000514
表示第J组相移因子,J为IRS阵元分组个数;
Figure BDA00035076541300000515
为IRS与第k个用户之间的信道;
Figure BDA00035076541300000516
为BS与IRS之间的信道。
可选地,分组后的IRS反射阵元的优化问题中,多个用户中的第k个用户对应的信干噪比,如下式(3)所示:
Figure BDA00035076541300000517
可选地,分组后的IRS反射阵元的优化问题P1,如下式(4)-(6)所示:
Figure BDA00035076541300000518
s.t.|φj|=1,n=1,2,...,J (5)
tr(PPH)≤Pmax (6)
其中,|φj|=1表示将IRS每一组反射阵元的反射幅度设为1;Pmax表示BS最大发射功率。
可选地,复杂度降低的优化问题P2,如下式(7)-(10)所示:
Figure BDA0003507654130000061
s.t.|φj|=1,n=1,2,...,J (8)
tr(PPH)≤Pmax (9)
αk≥0,k=1,2,...,K (10)
其中,α=[α1,...,αK]T为引入的辅助变量。
可选地,智能反射面,进一步用于:
令变量P和
Figure BDA0003507654130000062
固定,对上述公式(6)中的变量αk求偏导,得到αk的最优解;
将αk的最优解代入上式(6)中,将优化问题P2转换为优化问题P3。
可选地,智能反射面,进一步用于:
Figure BDA0003507654130000063
为等效信道,将优化问题P3重新表示为优化问题P4。
给定αk
Figure BDA0003507654130000064
引入辅助变量β=[β1,...,βK]T,对优化问题P4进行二次转换,得到优化问题P5。
给定P,通过对优化问题P5求偏导,得到βk的最优解。
固定βk,通过拉格朗日乘子法得到pk的最优解。
可选地,智能反射面,进一步用于:
给定α和P,引入辅助变量ρ=[ρ1,...,ρK]T,对优化问题P5进行二次转换,得到优化问题P6。
给定ρ,将优化问题P6重新表示为优化问题P7。
通过约束松弛法将优化问题P7重新表示为优化问题P8。
根据拉格朗日对偶分解法,将优化问题P8重新表示为优化问题P9。
对优化问题P9求解,最终得到
Figure BDA0003507654130000065
的最优解。
本发明实施例提供的技术方案带来的有益效果至少包括:
上述方案中,提出的用于无线通信的智能反射面分组优化模型,考虑了在IRS不同阵元分组数量和不同用户数量时对***下行速率的影响,并将IRS阵元不同分组数量与IRS不分组两种情况在平均每次优化时长方面作出了对比。可以在保证下行用户总速率的前提下,大大降低模型优化时间。
附图说明
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本发明实施例提供的基于智能反射面阵元分组的波束优化方法流程示意图;
图2是本发明实施例提供的IRS反射阵元分组示意图;
图3是本发明实施例提供的IRS阵元不同分组数量与IRS不分组平均每次优化时长对比图;
图4是本发明实施例提供的IRS阵元不同分组数量不同用户时可达到的下行用户总速率示意图;
图5是本发明实施例提供的基于智能反射面阵元分组的波束优化***框图。
具体实施方式
为使本发明要解决的技术问题、技术方案和优点更加清楚,下面将结合附图及具体实施例进行详细描述。
如图1所示,本发明实施例提供了一种基于智能反射面阵元分组的波束优化方法,该方法可以由基于智能反射面阵元分组的波束优化***实现。如图1所示的基于智能反射面阵元分组的波束优化方法流程图,该方法的处理流程可以包括如下的步骤:
S1、构建基站与多个用户间通过智能反射面进行信号传输的初始优化问题;初始优化问题为优化多个用户的总速率。
一种可行的实施方式中,考虑下行多用户无线通信***,用户和BS之间的直接链路被障碍物阻挡,其中基站配备有N根天线,K个用户且每个用户配有一根天线,一个IRS具有M个反射阵元,多个用户中的第k个用户的接收信号如下式(1)所示:
Figure BDA0003507654130000081
其中,
Figure BDA0003507654130000082
为IRS和第k个用户之间的信道,M是一个IRS具有的反射阵元个数,hM,k是第M个反射阵元和第k个用户之间的信道,T为矩阵的转置,H为矩阵的转置共轭;Φ=ηdiag([θ1,...,θM]T)为IRS反射系数表示的矩阵,η∈[0,1]表示IRS中第M个反射阵元的开或关状态,
Figure BDA0003507654130000083
Figure BDA0003507654130000084
表示第M个阵元的反射系数;
Figure BDA0003507654130000085
是BS和IRS之间的信道,N是基站配备有天线根数,
Figure BDA0003507654130000086
是BS和第M个反射阵元之间的信道,设定所有信道都是Saleh-Valenzuela毫米波信道;P=[P1,...PK]表示波束成形矩阵,其中
Figure BDA0003507654130000087
表示发送给第k个用户信号的波束成形矢量,K是用户个数;sk表示BS发送给第k个用户的信号,满足
Figure BDA0003507654130000088
k∈{1,2,...,K};
Figure BDA0003507654130000089
表示发送给第j个用户信号的波束成形矢量;sj表示BS发送给第j个用户的信号,满足
Figure BDA00035076541300000810
j∈{1,2,...,K};μk表示均值为0和方差为
Figure BDA00035076541300000811
的加性高斯白噪声。
如图2-4所示,对于IRS阵元分组的情况,假设所有阵元被分成J组,每组包含相同的反射阵元数量B,
Figure BDA00035076541300000812
则(j-1)B+1到jB之间的反射阵元属于第j组,假定一组中所有阵元都具有相同的相移因子,所以IRS相移矩阵可以重新表示为下式(2):
Figure BDA00035076541300000813
其中,
Figure BDA00035076541300000814
代表IRS分组后的相移因子表示成的矩阵,
Figure BDA00035076541300000815
表示第J组相移因子。分组后,IRS与第k个用户之间的信道可以表示为
Figure BDA00035076541300000816
BS与IRS之间的信道可以表示为
Figure BDA00035076541300000817
所以分组后的IRS反射阵元的优化问题中,多个用户中的第k个用户的接收信号被重写为下式(3):
Figure BDA00035076541300000818
Figure BDA00035076541300000819
道;
Figure BDA0003507654130000091
为BS与IRS之间的信道。
可选地,分组后的IRS反射阵元的优化问题中,多个用户中的第k个用户对应的SINR(Signal-to-interference-plus-noise Ratio,信干噪比)可由下式(4)求得:
Figure BDA0003507654130000092
S2、采用克罗内克积方法对智能反射面的反射阵元进行分组;将初始优化问题转化为分组后的IRS反射阵元的优化问题P1。
一种可行的实施方式中,在数学上,相应的优化问题可以表示为下式(5)-(7):
Figure BDA0003507654130000093
s.t.|φj|=1,n=1,2,...,J (6)
tr(PPH)≤Pmax (7)
其中,|φj|=1表示将IRS每一组反射阵元的反射幅度设为1,Pmax表示BS最大发射功率,上述问题是一个非凸问题。
S3、利用分式规划理论将优化问题P1转化为复杂度降低的优化问题P2。
S4、引入辅助变量,并将优化问题P2分解为三个解耦的子问题,使用交替迭代优化方法求出每个子问题闭合形式的最优解。
一种可行的实施方式中,为了解决上述这个非凸优化问题,可以利用分式规划理论将其转化为一个复杂度低的多的问题,如下式(8)-(11)所示:
Figure BDA0003507654130000094
s.t.|φj|=1,n=1,2,...,J (9)
tr(PPH)≤Pmax (10)
αk≥0,k=1,2,...,K (11)
其中α=[α1,...,αK]T是问题(P1)转化为(P2)时引入的辅助变量。本申请将该问题分解为三个解耦的子问题,并使用交替迭代优化方法求出每个子问题闭合形式的最优解。
可选地,αk的最优解的求解方法包括:
首先,当变量P和
Figure BDA0003507654130000101
固定时,将公式(8)对变量αk求偏导,令求导结果为0,可以得到αk的最优解为
Figure BDA0003507654130000102
得到αk的最优解后,将其带入公式(8)可以将问题(P2)转化为下式(12)-(14):
Figure BDA0003507654130000103
s.t.|φj|=1,n=1,2,...,J (13)
tr(PPH)≤Pmax (14)
因此,可以通过交替优化P,
Figure BDA0003507654130000104
来求解问题(P3)。
可选地,P的最优解的求解方法包括:
为了便于表示,将
Figure BDA0003507654130000105
表示为等效信道,所以,问题(P3)可以被重新表示为下式(15)-(16):
Figure BDA0003507654130000106
s.t.tr(ppH)≤Pmax (16)
上述优化问题在给定αk
Figure BDA0003507654130000107
时可以通过二次转换表示为下式(17)-(18):
Figure BDA0003507654130000108
s.t.tr(ppH)≤Pmax (18)
其中β=[β1,...,βK]T是利用二次转换时引入的辅助变量,当P给定时,公式(17)对βk求偏导,并令求导结果等于0,可以求得βk的最优解如下式(19)所示:
Figure BDA0003507654130000109
接下来,固定βk,pk的最优解可以通过拉格朗日乘子法得到,如下式(20)所示:
Figure BDA00035076541300001010
其中,μ≥0是功率约束(18)的拉格朗日乘子,且μ的最优解如下式(21)所示:
μop={μ≥0:tr(PPH)=Pmax} (21)
可选地,
Figure BDA00035076541300001011
的最优解的求解方法包括:
为了便于表示,将
Figure BDA0003507654130000111
重写下式(22):
Figure BDA0003507654130000112
其中,θ=[θ1,...,θJ]T。定义
Figure BDA0003507654130000113
当α和P给定时,通过二次规划可以将原问题写为下式(23):
Figure BDA0003507654130000114
其中
Figure BDA0003507654130000115
上述问题通过约束松弛法写为下式(24)-(25):
Figure BDA0003507654130000116
Figure BDA0003507654130000117
上述问题可以用拉格朗日对偶分解法得到,如下式(26)-(27)所示:
Figure BDA0003507654130000118
Figure BDA0003507654130000119
将公式(26)对θ求导,令求导结果为0,可以得到θ的最优解为下式(28):
Figure BDA00035076541300001110
其中,
Figure BDA00035076541300001114
将式(28)带入式(26)可以得到下式(29):
Figure BDA00035076541300001111
上述问题可以采用现有技术中的舒尔补码改写为半定编程问题并用CVX工具包求解,得到下式(30)-(31):
Figure BDA00035076541300001112
Figure BDA00035076541300001113
S5、根据最优解,得到联合优化后的基站BS波束成形矢量和IRS被动波束成形矢量,最终得到智能反射面阵元分组的波束优化方法。
本发明实施例中,提出的用于无线通信的智能反射面分组优化模型,考虑了在IRS不同阵元分组数量和不同用户数量时对***下行速率的影响,并将IRS阵元不同分组数量与IRS不分组两种情况在平均每次优化时长方面作出了对比。可以在保证下行用户总速率的前提下,大大降低模型优化时间。
如图5所示,本发明实施例提供了一种基于智能反射面阵元分组的波束优化***,该***应用于实现基于智能反射面阵元分组的波束优化方法,该***包括基站、智能反射面RIS以及多个用户。如图5所示的基于智能反射面阵元分组的波束优化***框图,其中:
基站,用于向多个用户传输信号。
智能反射面RIS,用于采用克罗内克积方法对智能反射面IRS的反射阵元进行分组;将初始优化问题转化为分组后的IRS反射阵元的优化问题P1;利用分式规划理论将优化问题P1转化为复杂度降低的优化问题P2;引入辅助变量,并将优化问题P2分解为三个解耦的子问题,使用交替迭代优化方法求出每个子问题闭合形式的最优解;根据最优解,得到联合优化后的BS波束成形矢量和IRS被动波束成形矢量,最终得到智能反射面阵元分组的波束优化方法。
多个用户,用于接收基站传输的信号。
可选地,初始优化问题中,多个用户中的第k个用户的接收信号,如下式(1)所示:
Figure BDA0003507654130000121
其中,
Figure BDA0003507654130000122
为IRS和第k个用户之间的信道,M是一个IRS具有的反射阵元个数,hM,k是第M个反射阵元和第k个用户之间的信道,T为矩阵的转置,H为矩阵的转置共轭;Φ=ηdiag([θ1,...,θM]T)为IRS反射系数表示的矩阵,η∈[0,1]表示IRS中第M个反射阵元的开或关状态,
Figure BDA0003507654130000123
Figure BDA0003507654130000124
表示第M个阵元的反射系数;
Figure BDA0003507654130000125
是BS和IRS之间的信道,N是基站配备有天线根数,
Figure BDA0003507654130000126
是BS和第M个反射阵元之间的信道,设定所有信道都是Saleh-Valenzuela毫米波信道;P=[P1,...PK]表示波束成形矩阵,其中
Figure BDA0003507654130000127
表示发送给第k个用户信号的波束成形矢量,K是用户个数;sk表示BS发送给第k个用户的信号,满足
Figure BDA0003507654130000128
k∈{1,2,...,K};
Figure BDA0003507654130000129
表示发送给第j个用户信号的波束成形矢量;sj表示BS发送给第j个用户的信号,满足
Figure BDA00035076541300001210
j∈{1,2,...,K};μk表示均值为0和方差为
Figure BDA00035076541300001211
的加性高斯白噪声。
可选地,分组后的IRS反射阵元的优化问题中,多个用户中的第k个用户的接收信号,如下式(2)所示:
Figure BDA0003507654130000131
其中,
Figure BDA0003507654130000132
为IRS分组后的相移因子表示成的矩阵,
Figure BDA0003507654130000133
表示第J组相移因子,J为IRS阵元分组个数;
Figure BDA0003507654130000134
为IRS与第k个用户之间的信道;
Figure BDA0003507654130000135
为BS与IRS之间的信道。
可选地,分组后的IRS反射阵元的优化问题中,多个用户中的第k个用户对应的信干噪比,如下式(3)所示:
Figure BDA0003507654130000136
可选地,分组后的IRS反射阵元的优化问题P1,如下式(4)-(6)所示:
Figure BDA0003507654130000137
s.t.|φj|=1,n=1,2,...,J (5)
tr(PPH)≤Pmax (6)
其中,|φj|=1表示将IRS每一组反射阵元的反射幅度设为1;Pmax表示BS最大发射功率。
可选地,复杂度降低的优化问题P2,如下式(7)-(10)所示:
Figure BDA0003507654130000138
s.t.|φj|=1,n=1,2,...,J (8)
tr(PPH)≤Pmax (9)
αk≥0,k=1,2,...,K (10)
其中,α=[α1,...,αK]T为引入的辅助变量。
可选地,智能反射面,进一步用于:
令变量P和
Figure BDA0003507654130000139
固定,对上述公式(6)中的变量αk求偏导,得到αk的最优解;
将αk的最优解代入上式(6)中,将优化问题P2转换为优化问题P3。
可选地,智能反射面,进一步用于:
Figure BDA00035076541300001310
为等效信道,将优化问题P3重新表示为优化问题P4。
给定αk
Figure BDA0003507654130000141
引入辅助变量β=[β1,...,βK]T,对优化问题P4进行二次转换,得到优化问题P5。
给定P,通过对优化问题P5求偏导,得到βk的最优解。
固定βk,通过拉格朗日乘子法得到pk的最优解。
可选地,智能反射面,进一步用于:
给定α和P,引入辅助变量ρ=[ρ1,...,ρK]T,对优化问题P5进行二次转换,得到优化问题P6。
给定ρ,将优化问题P6重新表示为优化问题P7。
通过约束松弛法将优化问题P7重新表示为优化问题P8。
根据拉格朗日对偶分解法,将优化问题P8重新表示为优化问题P9。
对优化问题P9求解,最终得到
Figure BDA0003507654130000142
的最优解。
本发明实施例中,提出的用于无线通信的智能反射面分组优化模型,考虑了在IRS不同阵元分组数量和不同用户数量时对***下行速率的影响,并将IRS阵元不同分组数量与IRS不分组两种情况在平均每次优化时长方面作出了对比。可以在保证下行用户总速率的前提下,大大降低模型优化时间。
本领域普通技术人员可以理解实现上述实施例的全部或部分步骤可以通过硬件来完成,也可以通过程序来指令相关的硬件完成,所述的程序可以存储于一种计算机可读存储介质中,上述提到的存储介质可以是只读存储器,磁盘或光盘等。
以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (10)

1.一种基于智能反射面阵元分组的波束优化方法,其特征在于,所述方法由基于智能反射面阵元分组的波束优化***实现,所述***包括基站、智能反射面RIS以及多个用户,所述方法包括:
S1、构建所述基站与多个用户间通过所述智能反射面IRS进行信号传输的初始优化问题;所述初始优化问题为优化多个用户的总速率;
S2、采用克罗内克积方法对智能反射面IRS的反射阵元进行分组;将所述初始优化问题转化为分组后的IRS反射阵元的优化问题P1;
S3、利用分式规划理论将所述优化问题P1转化为复杂度降低的优化问题P2;
S4、引入辅助变量,并将所述优化问题P2分解为三个解耦的子问题,使用交替迭代优化方法求出每个子问题闭合形式的最优解;
S5、根据所述最优解,得到联合优化后的基站BS波束成形矢量和IRS被动波束成形矢量,最终得到智能反射面阵元分组的波束优化方法。
2.根据权利要求1所述的方法,其特征在于,所述初始优化问题中,多个用户中的第k个用户的接收信号,如下式(1)所示:
Figure FDA0003507654120000011
其中,hk=[h1,k,...,hM,k]T∈£M×1为IRS和第k个用户之间的信道,M是IRS具有的反射阵元个数,hM,k是第M个反射阵元和第k个用户之间的信道,T为矩阵的转置,H为矩阵的转置共轭;Φ=ηdiag([θ1,...,θM]T)为IRS反射系数表示的矩阵,η∈[0,1]表示IRS中第M个反射阵元的开或关状态,
Figure FDA0003507654120000012
Figure FDA0003507654120000013
表示第M个阵元的反射系数;W=[w1,...,wM]T∈£M×N是BS和IRS之间的信道,N是基站配备有天线根数,wM∈£1×N是BS和第M个反射阵元之间的信道,设定所有信道都是Saleh-Valenzuela毫米波信道;P=[P1,...PK]表示波束成形矩阵,其中pk∈£N×1表示发送给第k个用户信号的波束成形矢量,K是用户个数;sk表示BS发送给第k个用户的信号,满足E{|sk|2}=1,k∈{1,2,K,K};pj∈£N×1表示发送给第j个用户信号的波束成形矢量;sj表示BS发送给第j个用户的信号,满足E{|sj|2}=1,j∈{1,2,K,K};μk表示均值为0和方差为
Figure FDA0003507654120000014
的加性高斯白噪声。
3.根据权利要求2所述的方法,其特征在于,分组后的IRS反射阵元的优化问题中,多个用户中的第k个用户的接收信号,如下式(2)所示:
Figure FDA0003507654120000021
其中,
Figure FDA0003507654120000022
为IRS分组后的相移因子表示成的矩阵,
Figure FDA0003507654120000023
表示第J组相移因子,J为IRS阵元分组个数;
Figure FDA0003507654120000024
为IRS与第k个用户之间的信道;
Figure FDA0003507654120000025
为BS与IRS之间的信道。
4.根据权利要求3所述的方法,其特征在于,分组后的IRS反射阵元的优化问题中,多个用户中的第k个用户对应的信干噪比,如下式(3)所示:
Figure FDA0003507654120000026
5.根据权利要求4所述的方法,其特征在于,分组后的IRS反射阵元的优化问题P1,如下式(4)-(6)所示:
Figure FDA0003507654120000027
s.t.|φj|=1,n=1,2,...,J (5)
tr(PPH)≤Pmax (6)
其中,|φj|=1表示将IRS每一组反射阵元的反射幅度设为1;Pmax表示BS最大发射功率。
6.根据权利要求5所述的方法,其特征在于,所述复杂度降低的优化问题P2,如下式(7)-(10)所示:
Figure FDA0003507654120000028
s.t.|φj|=1,n=1,2,...,J (8)
tr(PPH)≤Pmax (9)
αk≥0,k=1,2,...,K (10)
其中,α=[α1,...,αK]T为引入的辅助变量。
7.根据权利要求6所述的方法,其特征在于,所述αk的最优解的求解方法包括:
令变量P和
Figure FDA0003507654120000029
固定,对上述公式(6)中的变量αk求偏导,得到αk的最优解;
将所述αk的最优解代入上式(6)中,将优化问题P2转换为优化问题P3。
8.根据权利要求7所述的方法,其特征在于,所述P的最优解的求解方法包括:
Figure FDA0003507654120000031
为等效信道,将优化问题P3重新表示为优化问题P4;
给定αk
Figure FDA0003507654120000032
引入辅助变量β=[β1,...,βK]T,对优化问题P4进行二次转换,得到优化问题P5;
给定P,通过对所述优化问题P5求偏导,得到βk的最优解;
固定βk,通过拉格朗日乘子法得到pk的最优解。
9.根据权利要求8所述的方法,其特征在于,所述
Figure FDA0003507654120000033
的最优解的求解方法包括:
给定α和P,引入辅助变量ρ=[ρ1,...,ρK]T,对优化问题P5进行二次转换,得到优化问题P6;
给定ρ,将优化问题P6重新表示为优化问题P7;
通过约束松弛法将优化问题P7重新表示为优化问题P8;
根据拉格朗日对偶分解法,将优化问题P8重新表示为优化问题P9;
对所述优化问题P9求解,最终得到
Figure FDA0003507654120000034
的最优解。
10.一种基于智能反射面阵元分组的波束优化***,所述***用于实现基于智能反射面阵元分组的波束优化方法,所述基于智能反射面阵元分组的波束优化***包括基站、智能反射面RIS以及多个用户;其中:
所述基站,用于向所述多个用户传输信号;
所述智能反射面RIS,用于采用克罗内克积方法对智能反射面IRS的反射阵元进行分组;将所述初始优化问题转化为分组后的IRS反射阵元的优化问题P1;利用分式规划理论将所述优化问题P1转化为复杂度降低的优化问题P2;引入辅助变量,并将所述优化问题P2分解为三个解耦的子问题,使用交替迭代优化方法求出每个子问题闭合形式的最优解;根据所述最优解,得到联合优化后的BS波束成形矢量和IRS被动波束成形矢量,最终得到智能反射面阵元分组的波束优化方法;
所述多个用户,用于接收所述基站传输的信号。
CN202210143464.2A 2022-02-16 2022-02-16 一种基于智能反射面阵元分组的波束优化方法及*** Pending CN114640379A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210143464.2A CN114640379A (zh) 2022-02-16 2022-02-16 一种基于智能反射面阵元分组的波束优化方法及***

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210143464.2A CN114640379A (zh) 2022-02-16 2022-02-16 一种基于智能反射面阵元分组的波束优化方法及***

Publications (1)

Publication Number Publication Date
CN114640379A true CN114640379A (zh) 2022-06-17

Family

ID=81945812

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210143464.2A Pending CN114640379A (zh) 2022-02-16 2022-02-16 一种基于智能反射面阵元分组的波束优化方法及***

Country Status (1)

Country Link
CN (1) CN114640379A (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115361043A (zh) * 2022-07-14 2022-11-18 北京交通大学 高铁毫米波通信***的通信控制方法及控制***
WO2024032614A1 (zh) * 2022-08-12 2024-02-15 中兴通讯股份有限公司 波束管理方法、波束配置方法、接入节点、中继设备

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115361043A (zh) * 2022-07-14 2022-11-18 北京交通大学 高铁毫米波通信***的通信控制方法及控制***
CN115361043B (zh) * 2022-07-14 2023-09-12 北京交通大学 高铁毫米波通信***的通信控制方法及控制***
WO2024032614A1 (zh) * 2022-08-12 2024-02-15 中兴通讯股份有限公司 波束管理方法、波束配置方法、接入节点、中继设备

Similar Documents

Publication Publication Date Title
Qi et al. Integrating sensing, computing, and communication in 6G wireless networks: Design and optimization
CN114640379A (zh) 一种基于智能反射面阵元分组的波束优化方法及***
CN113726383A (zh) 一种智能反射面辅助的无线通信***
CN114070365B (zh) 一种智能反射表面辅助的低射频复杂度多用户mimo上行频谱效率优化方法
CN114285445B (zh) 电磁辐射约束的ris-dma辅助多用户mimo上行谱效优化方法
CN113904742A (zh) 非正交多址接入网络中基于智能反射面的资源分配方法
CN113556164A (zh) Irs辅助的swipt***中基于能效优先的波束成型优化方法
CN113765553A (zh) 一种基于机器学习的多波束卫星通信***鲁棒预编码方法
CN111901812A (zh) 一种全双工蜂窝通信网基站和智慧反射面联合控制方法
CN115173901A (zh) 基于irs辅助的miso无线携能通信***的能效最大化方法
CN114785388A (zh) 智能全向面辅助的多用户大规模simo上行m阶调制加权和速率优化方法
CN116033461B (zh) 一种基于star-ris辅助的共生无线电传输方法
CN115664471B (zh) 基于宽学习的毫米波mimo基站协作波束选择方法
CN114844538B (zh) 一种基于宽学习的毫米波mimo用户增量协作波束选择方法
CN114337902B (zh) 一种irs辅助的毫米波多小区间干扰的抑制方法
CN116419245A (zh) 一种基于智能反射面辅助的速率分割多址接入的多小区通信***能效优化方法
CN114285443B (zh) 动态超表面天线辅助的近场宽带上行mimo传输方法
CN115334524B (zh) 一种基于全向智能超表面的通信和雷达目标检测方法
CN114785387B (zh) 智能全向面辅助的多用户miso下行加权和速率优化方法
CN113992254B (zh) 智能反射面辅助无线通信***波束成形方法
CN115802466A (zh) 一种基于分布式ris辅助多用户***的联合功率分配和相移设计方法
CN115733533A (zh) 基于流形的irs辅助去蜂窝大规模mimo***波束赋形方法
CN114844537A (zh) 深度学习辅助鲁棒大规模mimo收发联合方法
CN114710186A (zh) 一种低成本广覆盖的毫米波预编码方法
Papazafeiropoulos et al. Achievable Rate Optimization for Large Stacked Intelligent Metasurfaces Based on Statistical CSI

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