本发明提出一种RIS辅助MIMO‑NOMA通信系统中的发送预编码和相移矩阵设计方法，在提出的通信系统中，首先设置RIS的初始相移矩阵，根据基站与两用户间的信道增益大小将用户分为强用户和弱用户；然后分别设计两个用户的发送预编码矩阵，并在此基础上通过二分法得到此时最优的基站发送功率分配方案；接下来，利用已有的结果设计RIS的相移矩阵并根据优化后的相移矩阵重新对发送预编码矩阵及功率分配方案进行迭代优化；最后当系统的和速率收敛时得到RIS相移矩阵、发送预编码矩阵及功率分配方案的全局最优解。本发明在通信环境受限的场景下，利用RIS和NOMA的优点，在保证用户通信公平性的同时，有效地提高系统的总传输速率。
1.一种RIS辅助MIMO‑NOMA通信系统中的发送预编码和相移矩阵设计方法，其特征在于，包括以下步骤：S1、构建基于统计CSI的RIS辅助MIMO‑NOMA下行无线传输系统，所述系统包括一个具有N根天线的基站、一块配备L个被动反射单元的RIS以及两个多天线用户；
S2、在系统中，设置RIS的初始相移矩阵，根据基站与两用户间的信道增益大小将用户分为强用户和弱用户，分别设计强用户和弱用户的发送预编码矩阵，并在此基础上通过二分法得到此时最优的基站发送功率分配方案；设定RIS相移矩阵和强用户的发送预编码矩阵的初始值为单位矩阵，然后设计两个用户的发送预编码矩阵和功率分配方案，具体操作如下：S201、设Q1和Q2是大小为N×N的矩阵，Q1、Q2分别表示强用户和弱用户的发送预编码矩阵，RIS的初始相移矩阵为Φ＝IL，其中，IL是L×L的单位矩阵；
S202、当强用户的发送预编码矩阵Q1已知时，设计弱用户的发送预编码矩阵Q2：其中，B2、Λ2均是辅助变量，具体表达式如下：‑1
式中， 分别是对矩阵B2进行奇异值分解得到的特征向量矩阵与特征值矩阵，IN为N×N的单位矩阵,μ2是令Q2满足基站发送功率限制的归一化参数，A21为辅助变量；
S203、基于步骤S202中得到的Q2，设计强用户的发送预编码矩阵Q1：GSVD
其中，Λ 是表示最佳功率分配的对角矩阵， 是对矩阵进行广义奇异值分解得到的特征向量矩阵，具体的广义奇异值分解形式为：
其中， 和 均为广义奇异值分解得到的特征向量矩阵， 和 均为广义奇异值分解得到的特征值矩阵，矩阵D、A1、B1均为信道近似相关矩阵，其表达式分别如下：‑1
式中， 为优化前的Q1，其初值为单位矩阵，A22为辅助变量，μ1、 均为强用户的等效信道参数，T1、T2分别是N×N、L×L的确定性非负矩阵， 分别为L×N、n×L的确定性矩阵，Γs、Ξs均为计算相关的辅助变量；
S204、基于步骤S202和S203中得到的最优的发送预编码矩阵Q1、Q2，设计强弱两用户间最优的功率分配方案，具体操作如下：S2041、假设基站的发送总功率为P，强用户的发送功率为P1，则弱用户的发送功率为P‑P1；
S2042、令P1＝P ,其中 P1,min＝0,P1,max＝P,根据P计算弱用户的传输速率R2：式中，Im为m×m的单位矩阵， 为m×L的确定性矩阵，R1、R3分别是L×L、m×m的确定性非负矩阵，T3是L×L的确定性非负矩阵， 是弱用户的等效信道参数，2
是强用户的等效信道参数，Θw、Θu均是计算相关的辅助变量，σ是噪声项；
如果R2≤R0，令P1,max＝P；反之，则令P1,min＝P ；其中R0是设定的用户正常通信所需的最小传输速率；
S2043、重复步骤S2042，直至P1,max‑P1,min≤ε，得到最终收敛的P ，那么强用户最优的发*送功率为P1＝P，弱用户最优的发送功率为P‑P1；其中，ε是表示收敛阈值的参数；
S3、利用已有的结果设计RIS的相移矩阵并根据优化后的相移矩阵重新对发送预编码矩阵及功率分配方案进行迭代优化；
S4、当系统的和速率收敛时得到RIS相移矩阵、发送预编码矩阵及功率分配方案的全局最优解。
2.根据权利要求1所述RIS辅助MIMO‑NOMA通信系统中的发送预编码和相移矩阵设计方法，其特征在于，所述步骤S2中，基站根据与用户间的信道增益大小，将用户划分为强用户和弱用户，其中强用户有n根天线，弱用户有m根天线；基站与用户之间通过RIS将基站发送的信号反射给用户，因此，基站与强用户之间的信道Hs和与弱用户之间的信道Hw分别为：Hs＝H2ΦH1,
式中，H1是L×N的矩阵，H2是n×L的矩阵，H3是m×L的矩阵，H1、H2、H3分别表示基站与RIS之间的信道、基站与强用户之间的信道以及基站与弱用户之间的信道；Φ是L×L的对角矩阵，表示RIS的相移矩阵；
为了进一步表示每条信道的散射分量和视距分离，每条信道可作如下划分：式中， 表示每条信道的视距分量； 表示每条信道的散射分量且
3.根据权利要求2所述RIS辅助MIMO‑NOMA通信系统中的发送预编码和相移矩阵设计方法，其特征在于，所述步骤S3中，设计RIS相移矩阵，具体操作如下：S301、给定发送预编码矩阵Q1、Q2，计算系统和速率Rsum关于RIS相移矩阵中对角元素φl的导数，得到梯度方向向量其中，φl为RIS相移矩阵中第l个对角元素，且l＝1,2,…,L；对于任意的l＝1,2,…,L，有如下表达式：其中， 为Rsum关于φl的偏导数, 是强用户的等效信道参数，{Γs,Ξs,Θs,Ψs,Πs}、{Γw,Ξw,Θw,Ψw,Πw}、{Γu,Ξu,Θu,Ψu,Πu}均是计算相关的中间变；
其中，j为虚数单位，α为每次迭代的步长，θ是由Φ的对角元素组成的列向量且是第t次迭代时RIS相移矩阵中第l个对角元素值；
S303、将步骤S302求得的RIS相移矩阵Φ代入步骤S202、S203和S204中，基于新的RIS相移矩阵重新设计发送预编码矩阵Q1、Q2以及功率分配方案，并利用所得结果再次优化RIS的相移矩阵，经过数次交替优化，计算系统的和速率；
其中，R1为强用户的传输速率，其计算表达式为：其中，In为n×n的单位矩阵；R2为弱用户的传输速率，可由步骤S2042中的表达式计算得到；
当系统的和速率Rsum收敛时，即可得到RIS相移矩阵、发送预编码矩阵Q1、Q2以及功率分配方案的全局最优解。
4.根据权利要求3所述RIS辅助MIMO‑NOMA通信系统中的发送预编码和相移矩阵设计方法，其特征在于，所述步骤S202中，A21的表达式为：式中， 分别为L×N、m×L的确定性矩阵， 均为若用户的等效信道参数，Γw、Ξw均为与计算相关的辅助变量，T1、T3分别是N×N、L×L的确定性非负矩阵。
5.根据权利要求4所述RIS辅助MIMO‑NOMA通信系统中的发送预编码和相移矩阵设计方法，其特征在于，弱用户的等效信道参数 其表达式分别如下：计算相关的辅助变量{Γw,Ξw,Θw,Ψw,Πw}，其表达式分别如下：2
式中，R1、R3分别是L×L、m×m的确定性非负矩阵，σ是噪声项，tr(·)表示矩阵的迹。
6.根据权利要求3所述RIS辅助MIMO‑NOMA通信系统中的发送预编码和相移矩阵设计方法，其特征在于，所述步骤S203中，A22的表达式为：式中， 为m×L的确定性矩阵，Γu、Ξu均为计算相关的辅助变量， 均为强用户的等效信道参数，T3是L×L的确定性非负矩阵。
7.根据权利要求6所述RIS辅助MIMO‑NOMA通信系统中的发送预编码和相移矩阵设计方法，其特征在于，强用户的等效信道参数 其表达式分别如下：计算相关的辅助变量{Γs,Ξs,Θs,Ψs,Πs}和{Γu,Ξu,Θu,Ψu,Πu}，其表达式分别如下：2
式中，R1、R2、R3分别是L×L、n×n、m×m的确定性非负矩阵，σ是噪声项，tr(·)表示矩阵的迹。
8.根据权利要求1所述RIS辅助MIMO‑NOMA通信系统中的发送预编码和相移矩阵设计方GSVD法，其特征在于，Λ 的对角元素为：
其中， 分别为对角矩阵 中的对角元素，μ1是令Q1满足基站发送功率限制的归一化参数，νi是矩阵 中的对角元素。
9.根据权利要求1所述RIS辅助MIMO‑NOMA通信系统中的发送预编码和相移矩阵设计方法，其特征在于，所述步骤S2041中，{Ωs,Ωw,Ωu}的表达式分别为其中，Εll为第l行l列的元素为1，其余元素全为0的L×L矩阵。
RIS辅助MIMO‑NOMA通信系统中的发送预编码和相移矩阵设计
方法
技术领域
[0001] 本发明属于物联网(Internet of Things，IoT)技术领域，具体涉及一种可重构智能反射表面(Reconfigurable Intelligent Surface，RIS)辅助多输入多输出非正交多址接入(Multiple Input Multiple Output‑Non‑Orthogonal Multiple Access，MIMO‑NOMA)通信系统中的发送预编码和相移矩阵设计方法。
背景技术
[0002] 随着IoT技术的空前发展，大量具有传感、计算和无线通信能力的IoT设备不断涌现，但由于尺寸和成本的限制，许多IoT设备携带的电池容量非常有限，在更换设备或电池困难且成本高昂的情况下，IoT设备的寿命成为其是否能长期使用的关键。因此，有效利用频谱和能量资源来延长IoT的寿命至关重要。非正交多址接入(Non‑Orthogonal Multiple Access，NOMA)作为IoT中提出的关键技术之一，因其能够显著提高IoT系统的频谱效率而受到广泛的研究。NOMA的基本概念是在发送端将用户信号叠加在同一资源块(如频率、时间)上，在接收端通过连续干扰消除技术过滤不需要的信号，实现正确解调。因此，与传统的正交多址方案相比，NOMA提供了更高的频谱效率以及更好的用户公平性。
[0003] 然而，在有些场景中，比如基站与用户之间因为存在障碍物阻挡，无法建立直接传输链路，仅凭NOMA无法满足正常传输的要求。而RIS可以通过控制信号的幅值或相位反射信号，且工作在全双工状态下，没有自干扰，能够在保证正常通信的同时，有效地提高频谱效率和系统的能量效率。
[0004] 本发明结合NOMA和RIS的优点提出一种RIS辅助MIMO‑NOMA通信系统中的发送预编码和相移矩阵设计方法，其难点在于发送预编码矩阵、RIS相移矩阵、用户间功率方案的设计。
发明内容
[0005] 本发明所要解决的技术问题是，克服现有技术的不足而提供一种RIS辅助的MIMO‑NOMA通信系统中的发送预编码和相移矩阵设计方法，在通信环境受限的场景下，利用RIS和NOMA的优点，在保证用户通信公平性的同时，有效地提高系统的通信性能。
[0006] 本发明提供一种RIS辅助MIMO‑NOMA通信系统中的发送预编码和相移矩阵设计方法，包括以下步骤：
[0007] S1、构建基于统计CSI的RIS辅助MIMO‑NOMA下行无线传输系统，所述系统包括一个具有N根天线的基站、一块配备L个被动反射单元的RIS以及两个多天线用户；
[0008] S2、在系统中，设置RIS的初始相移矩阵，根据基站与两用户间的信道增益大小将用户分为强用户和弱用户；
[0009] S3、分别设计强用户和弱用户的发送预编码矩阵，并在此基础上通过二分法得到此时最优的基站发送功率分配方案；
[0010] S4、利用已有的结果设计RIS的相移矩阵并根据优化后的相移矩阵重新对发送预编码矩阵及功率分配方案进行迭代优化；
[0011] S5、当系统的和速率收敛时得到RIS相移矩阵、发送预编码矩阵及功率分配方案的全局最优解。
[0012] 本发明在提出的通信系统中，首先设置RIS的初始相移矩阵，根据基站与两用户间的信道增益大小将用户分为强用户和弱用户；然后分别设计两个用户的发送预编码矩阵，并在此基础上通过二分法得到此时最优的基站发送功率分配方案；接下来，利用已有的结果设计RIS的相移矩阵，并根据优化后的相移矩阵重新对发送预编码矩阵以及功率分配方案进行迭代优化；最后，当系统的和速率收敛时，得到RIS相移矩阵、发送预编码矩阵以及功率分配方案的全局最优解。本发明可以在通信环境受限的场景下，利用RIS和NOMA的优点，在保证用户通信公平性的同时，有效地提高系统的总传输速率。
[0013] 本发明进一步优化的技术方案如下：
[0014] 所述步骤S2中，基站根据与用户间的信道增益大小，将用户划分为强用户和弱用户，其中强用户有n根天线，弱用户有m根天线。基站与用户之间因障碍物阻碍没有直接链路，需要通过RIS将基站发送的信号反射给用户，因此，基站与强用户之间的信道Hs和与弱用户之间的信道Hw分别为：
[0015] Hs＝H2ΦH1,
[0016] Hw＝H3ΦH1,
[0017] 式中，H1是L×N的矩阵，H2是n×L的矩阵，H3是m×L的矩阵，H1、H2、H3分别表示基站与RIS之间的信道、基站与强用户之间的信道以及基站与弱用户之间的信道；Φ是L×L的对角矩阵，表示RIS的相移矩阵。为了进一步表示每条信道的散射分量和视距分离，每条信道可作如下划分：
[0018]
[0019] 式中， 表示每条信道的视距分量， 分别为L×N、n×L、m×L的确定性矩阵； 表示每条信道的散射分量且 其中R1、R2、R3分别是L×L、n×n、m×m的确定性非负矩阵，表示接收天线相关矩阵；T1、T2、T3分别是N×N、L×L、L×L的确定性非负矩阵，表示发射天线相关矩阵；X1、X2、X3表示信道的随机分量部分，且X1服从均值为0方差为的复高斯分布，X2和X3服从均值为0方差为 的复高斯分布； 表示矩阵的平方根运算。
[0020] 所述步骤S2中，设定RIS相移矩阵和强用户的发送预编码矩阵的初始值为单位矩阵，然后设计两个用户的发送预编码矩阵和功率分配方案，具体操作如下：
[0021] S201、设Q1和Q2是大小为N×N的矩阵，Q1、Q2分别表示强用户和弱用户的发送预编码矩阵，两个用户的功率分配方案为等功率分配，RIS的初始相移矩阵为Φ＝IL，其中，IL是L×L的单位矩阵；
[0022] S202、当强用户的发送预编码矩阵Q1已知时，设计弱用户的发送预编码矩阵Q2：
[0023]
[0024] 其中，B2、Λ2均是辅助变量，具体表达式如下：
[0025] B2＝(IN+A21Q1)‑1A21,
[0026]
[0027] 式中， 分别是对矩阵B2进行奇异值分解 得到的特征向量矩阵与特征值矩阵，IN为N×N的单位矩阵,μ2是令Q2满足基站发送功率限制的归一化参数，H ‑1 +
(·) 表示矩阵的共轭转置，(·) 表示对矩阵求逆，(·) 表示取括号中数据与0相比的最大值，A21为辅助变量；
[0028] S203、基于步骤S202中得到的Q2，设计强用户的发送预编码矩阵Q1：
[0029]GSVD
[0030] 其中，Λ 是表示最佳功率分配的对角矩阵， 是对矩阵进行广义奇异值分解得到的特征向量矩阵，具体的广义奇异
值分解形式为：
[0031]
[0032]
[0033] 其中， 和 均为广义奇异值分解得到的特征向量矩阵， 和 均为广义奇异值分解得到的特征值矩阵，矩阵D、A1、B1均为信道近似相关矩阵，其表达式分别如下：
[0034]
[0035]
[0036] B1＝(IN+A21Q2)‑1A21,
[0037] 式中， 为优化前的Q1，其初值为单位矩阵；A22为辅助变量，μ1、 均为强用户的等效信道参数，T1、T2分别是N×N、L×L的确定性非负矩阵， 分别为L×N、n×L的确定性矩阵，Γs、Ξs均为计算相关的辅助变量；
[0038] S204、基于步骤S202和S203中得到的最优的发送预编码矩阵Q1、Q2，设计强弱两用户间最优的功率分配方案，具体操作如下：
[0039] S2041、假设基站的发送总功率为P，强用户的发送功率为P1，则弱用户的发送功率为P‑P1；
[0040] S2042、令P1＝P*,其中 P1,min＝0,P1,max＝P,根据P*计算弱用户的传输速率R2：
[0041]
[0042] 式中，Im为m×m的单位矩阵， 为m×L的确定性矩阵，R1、R3分别是L×L、m×m的确定性非负矩阵，T3是L×L的确定性非负矩阵， 是弱用户的等效信道参数，2
是强用户的等效信道参数，Θw、Θu均是计算相关的辅助变量，σ是噪声项；
[0043] 如果R2≤R0，令P1,max＝P*；反之，则令P1,min＝P*；其中R0是设定的用户正常通信所需的最小传输速率；
[0044] S2043、重复步骤S2042，直至P1,max‑P1,min≤ε，得到最终收敛的P*，那么强用户最优*的发送功率为P1＝P，弱用户最优的发送功率为P‑P1；其中，ε是表示收敛阈值的参数。
[0045] 所述步骤S4中，设计RIS相移矩阵，具体操作如下：
[0046] S301、给定发送预编码矩阵Q1、Q2，计算系统和速率Rsum关于RIS相移矩阵中对角元素φl的导数，得到梯度方向向量
[0047]
[0048] 其中，φl为RIS相移矩阵中第l个对角元素，且l＝1,2,…,L；对于任意的l＝1,2,…,L，有如下表达式：
[0049]
[0050] 其中， 为Rsum关于φl的偏导数， 是强用户的等效信道参数，{Γs,Ξs,Θs,Ψs,Πs}、{Γw,Ξw,Θw,Ψw,Πw}、{Γu,Ξu,Θu,Ψu,Πu}均是计算相关的中间变量；
[0051] S302、计算RIS的相移矩阵Φ：
[0052]
[0053] 其中，j为虚数单位，α为每次迭代的步长，θ是由Φ的对角元素组成的列向量且是第t次迭代时RIS相移矩阵中第l个对角元素值，diag(·)表示括号中元素为对角线元素的对角矩阵，exp(·)表示以自然常数e为底的指数函数，arg(·)表示括号中元素的辐角；
[0054] S303、将步骤S302求得的RIS相移矩阵Φ代入步骤S202、S203和S204中，基于新的RIS相移矩阵重新设计发送预编码矩阵Q1、Q2以及功率分配方案，并利用所得结果再次优化RIS的相移矩阵，经过数次交替优化，计算系统的和速率
[0055] 所述步骤S4中，计算系统的和速率：
[0056] Rsum＝R1+R2,
[0057] 其中，R1为强用户的传输速率，其计算表达式为：
[0058]
[0059] 其中，In为n×n的单位矩阵；R2为弱用户的传输速率，可由步骤S2042中的表达式计算得到；
[0060] 当系统的和速率Rsum收敛时，即可得到RIS相移矩阵、发送预编码矩阵Q1、Q2以及功率分配方案的全局最优解。
[0061] 所述步骤S202中，A21的表达式为：
[0062]
[0063] 式中， 分别为L×N、m×L的确定性矩阵， 均为若用户的等效信道参数，Γw、Ξw均为与计算相关的辅助变量，T1、T3分别是N×N、L×L的确定性非负矩阵。
[0064] 进一步的，弱用户的等效信道参数 其表达式分别如下：
[0065]
[0066]
[0067]
[0068]
[0069] 计算相关的辅助变量{Γw,Ξw,Θw,Ψw,Πw}，其表达式分别如下：
[0070]
[0071] Ξw＝σ2IN+w2R3,
[0072] Θw＝IL+w1ΨwR1,
[0073]
[0074] 式中，R1、R3分别是L×L、m×m的确定性非负矩阵，σ2是噪声项,tr(·)表示矩阵的迹。
[0075] 所述步骤S203中，A22的表达式为：
[0076]
[0077] 式中， 为m×L的确定性矩阵，Γu、Ξu均为计算相关的辅助变量， 均为强用户的等效信道参数，T3是L×L的确定性非负矩阵。
[0078] 进一步的，强用户的等效信道参数 其表达式分别如下：
[0079]
[0080]
[0081]
[0082]
[0083]
[0084]
[0085]
[0086]
[0087] 计算相关的辅助变量{Γs,Ξs,Θs,Ψs,Πs}和{Γu,Ξu,Θu,Ψu,Πu}，其表达式分别如下：
[0088]
[0089] Ξs＝σ2IN+s2R2,
[0090] Θs＝IL+s1ΨsR1,
[0091]
[0092]
[0093]
[0094] Ξu＝σ2IN+u2R3,
[0095] Θu＝IL+u1ΨuR1,
[0096]
[0097]
[0098] 式中，R1、R2、R3分别是L×L、n×n、m×m的确定性非负矩阵，σ2是噪声项,tr(·)表示矩阵的迹。
[0099] 进一步的，ΛGSVD的对角元素为：
[0100]
[0101] 其中， 分别为对角矩阵 中的对角元素，μ1是令Q1满足基站发送功率限制的归一化参数，νi是矩阵 中的对角元素。
[0102] 所述步骤S2041中，{Ωs,Ωw,Ωu}的表达式分别为
[0103]
[0104]
[0105]
[0106] 其中，Εll为第l行l列的元素为1，其余元素全为0的L×L矩阵。
[0107] 本发明采用以上技术方案与现有技术相比，具有以下技术效果：
[0108] (1)本发明结合了RIS和NOMA的优点，在通信受限的场景下，利用RIS能够反射信号的优点，建立了更贴合实际的MIMO‑NOMA通信系统；
[0109] (2)本发明设计发送预编码和相移矩阵时，考虑了发送天线与接收天线的空间相关性，通过迭代优化获得了得到RIS相移矩阵、发送预编码矩阵以及功率分配方案的全局最优解，使得系统的总传输速率能够得到最大程度的提升。
附图说明
[0110] 图1为本发明的流程图。
[0111] 图2为本发明中RIS辅助MIMO‑NOMA通信系统的示意图。
具体实施方式
[0112] 下面结合附图对本发明的技术方案做进一步的详细说明：本实施例在以本发明技术方案为前提下进行实施，给出了详细的实施方式和具体的操作过程，但本发明的保护权限不限于下述的实施例。
[0113] 本实施例提出了一种RIS辅助MIMO‑NOMA通信系统中的发送预编码和相移矩阵设计方法，如图1所示，包括以下步骤：
[0114] 步骤1：构建基于统计CSI的RIS辅助MIMO‑NOMA下行无线传输系统，设置RIS的初始相移矩阵，根据基站与两用户间的信道增益大小将用户分为强用户和弱用户，然后分别设计两个用户的发送预编码矩阵。
[0115] 步骤1.1：构建基于统计CSI的RIS辅助MIMO‑NOMA下行无线传输系统，如图2所示所述系统包括一个具有N根天线的基站、一块配备L个被动反射单元的RIS以及两个多天线用户，基站根据与用户间的信道增益大小，将用户划分为强用户和弱用户，其中强用户有n根天线，弱用户有m根天线。基站与用户之间因障碍物阻碍没有直接链路，需要通过RIS将基站发送的信号反射给用户，因此，基站与强用户之间的信道Hs和与弱用户之间的信道Hw分别为：
[0116] Hs＝H2ΦH1,
[0117] Hw＝H3ΦH1,
[0118] 式中，H1是L×N的矩阵，H2是n×L的矩阵，H3是m×L的矩阵，H1、H2、H3分别表示基站与RIS之间的信道、基站与强用户之间的信道以及基站与弱用户之间的信道；Φ是L×L的对角矩阵，表示RIS的相移矩阵。为了进一步表示每条信道的散射分量和视距分离，每条信道可作如下划分：
[0119]
[0120] 式中， 表示每条信道的视距分量， 分别为L×N、n×L、m×L的确定性矩阵； 表示每条信道的散射分量且 其中R1、R2、R3分别是L×L、n×n、m×
m的确定性非负矩阵，表示接收天线相关矩阵；T1、T2、T3分别是N×N、L×L、L×L的确定性非负矩阵，表示发射天线相关矩阵；X1、X2、X3表示信道的随机分量部分，且X1服从均值为0方差为 的复高斯分布，X2和X3服从均值为0方差为 的复高斯分布； 表示矩阵的平方根运算。
[0121] 步骤1.2：设定RIS相移矩阵和强用户的发送预编码矩阵的初始值为单位矩阵，然后设计两个用户的发送预编码矩阵和功率分配方案，具体操作如下：
[0122] 步骤1.21：设Q1和Q2是大小为N×N的矩阵，Q1、Q2分别表示强用户和弱用户的发送预编码矩阵，两个用户的功率分配方案为等功率分配，RIS的初始相移矩阵为Φ＝IL，其中，IL是L×L的单位矩阵；
[0123] 步骤1.22：当强用户的发送预编码矩阵Q1已知时，设计弱用户的发送预编码矩阵Q2：
[0124]
[0125] 其中，B2、Λ2均是辅助变量，具体表达式如下：
[0126] B2＝(IN+A21Q1)‑1A21,
[0127]
[0128] 式中， 分别是对矩阵B2进行奇异值分解 得到的特征向量H
矩阵与特征值矩阵，μ2是令Q2满足基站发送功率限制的归一化参数，(·) 表示矩阵的共轭‑1 +
转置，(·) 表示对矩阵求逆，(·)表示取括号中数据与0相比的最大值，A21为辅助变量，其表达式为：
[0129]
[0130] 其中， 是弱用户的等效信道参数，表达式分别如下：
[0131]
[0132]
[0133]
[0134]
[0135] 其中，{Γw,Ξw,Θw,Ψw,Πw}均为计算相关的辅助变量，表达式分别如下：
[0136]
[0137] Ξw＝σ2IN+w2R3,
[0138] Θw＝IL+w1ΨwR1,
[0139]
[0140] 式中，IN为N×N的单位矩阵,σ2是噪声项,tr(·)表示矩阵的迹。
[0141] 步骤2：强用户的发送预编码矩阵和两用户间功率分配方案的设计,包括以下步骤：
[0142] 步骤2.1：基于步骤1.22中得到的Q2，设计强用户的发送预编码矩阵Q1：
[0143]
[0144] 其中， 是对矩阵 进行广义奇异值分解得到的特征向量矩阵，具体的广义奇异值分解形式为：
[0145]
[0146]
[0147] 其中， 和 均为广义奇异值分解得到的特征向量矩阵， 和 均为广义奇异值分解得到的特征值矩阵，矩阵D、A1、B1均为信道近似相关矩阵，其表达式分别如下：
[0148]
[0149]
[0150] B1＝(IN+A21Q2)‑1A21,
[0151] 式中， 为优化前的Q1，其初值为单位矩阵；A22为辅助变量，其表达式为：
[0152]
[0153] 其中， 是强用户的等效信道参数，其表达式分别如下：
[0154]
[0155]
[0156]
[0157]
[0158]
[0159]
[0160]
[0161]
[0162] 其中，{Γs,Ξs,Θs,Ψs,Πs}和{Γu,Ξu,Θu,Ψu,Πu}是计算相关的辅助变量，其表达式分别如下：
[0163]
[0164] Ξs＝σ2IN+s2R2,
[0165] Θs＝IL+s1ΨsR1,
[0166]
[0167]
[0168]
[0169] Ξu＝σ2IN+u2R3,
[0170] Θu＝IL+u1ΨuR1,
[0171]
[0172]
[0173] 其中，ΛGSVD是表示最佳功率分配的对角矩阵，其对角元素为：
[0174]
[0175] 其中， 分别为对角矩阵 中的对角元素，μ1是令Q1满足基站发送功率限制的归一化参数，νi是矩阵 中的对角元素。
[0176] 步骤2.2：基于步骤:2.1和步骤1.22中得到的最优的发送预编码矩阵Q1、Q2，设计强弱两用户间最优的功率分配方案，包括以下步骤：
[0177] 步骤2.21：假设基站的发送总功率为P，强用户的发送功率为P1，则弱用户的发送功率为P‑P1；
[0178] 步骤2.22：令P1＝P*,其中 P1,min＝0,P1,max＝P,根据P*计算弱用户的传输速率R2：
[0179]
[0180] 如果R2≤R0，令P1,max＝P*；反之，则令P1,min＝P*；其中，Im为m×m的单位矩阵，R0是设定的用户正常通信所需的最小传输速率；
[0181] 步骤2023：重复步骤2.22，直至P1,max‑P1,min≤ε，得到最终收敛的P*，那么强用户最*优的发送功率为P1＝P，弱用户最优的发送功率为P‑P1；其中，ε是表示收敛阈值的参数。
[0182] 步骤3：RIS相移矩阵的设计，包括以下步骤：
[0183] 步骤3.1：给定发送预编码矩阵Q1、Q2，计算系统和速率Rsum关于RIS相移矩阵中对角元素φl的导数，得到梯度方向向量
[0184]
[0185] 其中，φl为RIS相移矩阵中第l个对角元素，且l＝1,2,…,L；对于任意的l＝1,2,…,L，有如下表达式：
[0186]
[0187] 其中， 为Rsum关于φl的偏导数，{Ωs,Ωw,Ωu}的表达式分别为
[0188]
[0189]
[0190]
[0191] 其中，Εll为第l行l列的元素为1，其余元素全为0的L×L矩阵。
[0192] 步骤3.2：计算RIS的相移矩阵Φ：
[0193]
[0194] 其中，j为虚数单位，α为每次迭代的步长，θ是由Φ的对角元素组成的列向量且是第t次迭代时RIS相移矩阵中第l个对角元素值，diag(·)表示括号中元素为对角线元素的对角矩阵，exp(·)表示以自然常数e为底的指数函数，arg(·)表示括号中元素的辐角；
[0195] 步骤3.3：将求得的RIS相移矩阵Φ代入步骤1.2中，基于新的RIS相移矩阵重新设计发送预编码矩阵Q1、Q2以及功率分配方案，并利用所得结果再次优化RIS的相移矩阵，经过数次交替优化，计算系统的和速率：
[0196] Rsum＝R1+R2,
[0197] 其中，R1为强用户的传输速率，其计算表达式为：
[0198]
[0199] 其中，In为n×n的单位矩阵；R2为弱用户的传输速率，可由步骤2.22中的表达式计算得到。当系统的和速率Rsum收敛时，即可得到RIS相移矩阵、发送预编码矩阵Q1、Q2以及功率分配方案的全局最优解。
[0200] 以上所述，仅为本发明中的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉该技术的人在本发明所揭露的技术范围内，可理解想到的变换或替换，都应涵盖在本发明的包含范围之内，因此，本发明的保护范围应该以权利要求书的保护范围为准。
