mimo wireless communications via intelligent reflecting surface[j].ieee transactions on communications,2020,68(12):7851-7866.doi:10.1109/tcomm.2020.3024621.]针对多输入多输出(multiple-input multiple-output,mimo)通信系统，同样以保密速率最大化为目标，利用块坐标下降(block coordinate descent,bcd)和majorization-minimization算法，联合优化发射机预编码矩阵和irs相移矩阵的设计。文献[yu xianghao,xu dongfang,schober robert.enabling secure wireless communications via intelligent reflecting surfaces[c]//2019ieee global communications conference(globecom).waikoloa,hi,usa,2019:1-6.doi:10.1109/globecom38437.2019.9014322.]和[cui miao,zhang guangchi,zhang rui.secure wireless communication via intelligent reflecting surface[j].ieee wireless communications letters,2019,8(5):1410-1414.doi:10.1109/lwc.2019.2919685.]分别利用bcd和最小最大化算法求解保密速率最大化问题，后者同时还使用sdr算法和线性分式变换将优化问题转化为凸问题。文献[huang chong,chen gaojie,wang kaikit.multi-agent reinforcement learning-based buffer-aided relay selection in irs-assisted secure cooperative networks[j].ieee transactions on information forensics and security,2021,16:4101-4112.doi:10.1109/tifs.2021.3103062.]针对基站与用户和窃听者之间存在irs和多个缓存辅助中继的无线通信系统，提出了一种irs辅助的安全协作网络中缓存辅助中继选择方案，在irs反射信号的相移和幅度只能离散取值的场景下，在缓存辅助中继产生的延迟约束下，以最大化系统保密速率为目标，利用深度学习的方法求解缓存辅助中继选择和irs相移矩阵的联合优化问题。
[0003]
上述文献对波束成形、irs相移、人工噪声等的优化的前提是基站完美已知所有信道的状态信息(channel state information,csi)。是否能获得信道的csi对无线通信系统中传输方案的设计和系统性能有重要的影响。在实际的场景中，由于信道估计误差、信道时变等原因，获得的信道csi难免存在误差。文献[zheng beixiong,zhang rui.intelligent reflecting surface-enhanced ofdm:channel estimation and reflection optimization[j].ieee wireless communications letters,2020,9(4):518-522.doi:10.1109/lwc.2019.2961357.]提出了一种对irs增强的正交频分复用系统进行信道估计的传输协议，文献[zhao mingming,liu an,zhang rui.outage-constrained robust beamforming for intelligent reflecting surface aided wireless communication[j].ieee transactions on signal processing,2021,69:1301-1316.doi:10.1109/tsp.2021.3056899.]利用上述文献中的信道估计方法估计irs辅助的多用户下行通信系统中的csi误差，在用户端信干噪比小于阈值的概率和irs相位离散取值的约束下，利用连续凸逼近(successive convex approximation,sca)等算法联合优化基站发射信号波束成形矢量和irs相移矩阵，最小化基站发射功率。文献[zhao mingming,wu qingqing,zhao minjian,et al.exploiting amplitude control in intelligent reflecting surface aided wireless communication with imperfect csi[j].ieee transactions on communications,2021,69(6):4216-4231.doi:10.1109/tcomm.2021.3064959.]针对不准确csi下irs辅助的多用户通信系统，在irs反射单元的反射信号幅度不同的情况下，以多用户信息加权和速率最大化为目标，联合优化优化基站发射信号波束成形矢量和irs相移矩
阵。文献采用了迭代块连续上界最小化算法求解优化问题，并通过仿真的方式分析了反射信号幅度变化对多用户信息加权和速率的影响。
[0004]
在物理层安全的研究中，一般假设合法信道的csi准确已知，而窃听信道的csi则分为已知、部分已知或未知等几种情况，针对不同的情况，需要采用不同的安全方案。如前述的文献就是假设窃听信道的csi完美已知，上述文献则针对信道csi不完美的场景下irs辅助的物理层安全通信进行了研究。针对多用户、多窃听者、多irs的下行通信系统，文献[yu xianghao,xu dongfang,sun ying,et al.robust and secure wireless communications via intelligent reflecting surfaces[j].ieee journal on selected areas in communications,2020,38(11):2637-2652.doi:10.1109/jsac.2020.3007043.]在基站获得的irs与窃听者之间信道的csi不完美的情况下，对提高系统安全性能的优化问题进行了研究。文献假设基站与用户之间存在障碍物的遮挡，利用有界csi误差模型刻画csi误差，在窃听者的拦截速率小于给定值的约束下，联合优化基站发射波束成形矩阵、人工噪声协方差矩阵和irs相移矩阵，最大化系统和速率。文献[yu xianghao,xu dongfang,sun ying,et al.robust and secure wireless communications via intelligent reflecting surfaces[j].ieee journal on selected areas in communications,2020,38(11):2637-2652.doi:10.1109/jsac.2020.3007043.]针对基站与用户之间有障碍物遮挡，需要通过irs进行信号反射形成通信链路的miso系统，利用统计误差模型描述基站-irs-窃听者之间级联信道csi的误差，在合法用户信息传输速率和窃听者拦截中断概率约束下，最小化基站发射功率。文献[dong limeng,wang huiming and xiao haitao.secure cognitive radio communication via intelligent reflecting surface[j].ieee transactions on communications,2021,69(7):4678-4690.doi:10.1109/tcomm.2021.3073028.]针对irs辅助的认知无线电通信系统，在基站-窃听者的直连信道和基站-irs-窃听者级联信道的csi不完美情况下，对提高系统安全性能的优化问题进行了研究。文献使用有界csi误差模型刻画csi误差，利用辅助变量和sca等算法联合优化基站发射信号波束成形矢量和irs相移矩阵，最大化认知用户可获得的保密速率。
技术实现要素：[0005]
本发明的目的在于研究irs辅助的多用户下行系统中的物理层安全的优化问题。基站与用户之间缺少直传链路，利用irs反射形成传输链路。多个用户之间信息需要相互保密，每个时隙，非信息传输的目标用户视为窃听者。由于信道的时变性，基站拥有窃听信道的csi为过时信息，与真实的csi间存在误差。在此条件下，以最坏情况下的保密速率最大化为目标，联合优化基站发射信息信号和人工噪声的波束成形矢量以及irs相移矩阵。
[0006]
为了实现上述目的本发明采用如下技术方案：首先根据系统模型构建优化问题，然后将非凸的优化问题转化为两层优化问题，其中第1层优化问题的求解中包含第2层优化问题的求解，进一步利用一维搜索算法求解第1层优化问题，利用算法将第2层优化问题分解为两个交替迭代优化的非凸子问题分别求解，最后获得最大系统保密速率。具体包括以下步骤：
[0007]
(1)构建通信系统模型，根据构建的系统模型得到合法用户与窃听用户处的信息传输速率，进而得到系统保密速率；
[0008]
(2)以系统保密速率最大化为优化目标，构造在基站总发射功率和智能反射表面反射单元单位模约束下，对基站发射信息信号和人工噪声波束成形矢量，以及智能反射表面的相移矩阵进行联合优化的数学模型；
[0009]
(3)步骤(2)所述数学模型为非凸的优化问题，将其转化为求解两层优化问题，其中第1层优化问题的求解中包含第2层优化问题的求解；
[0010]
(4)第1层优化问题为单变量优化问题，利用一维搜索算法求解；
[0011]
(5)将第2层优化问题分解为两个交替迭代优化的非凸子问题分别求解。
[0012]
进一步，步骤(3)所述求解两层优化问题具体算法为：设置初始迭代次数i＝1，松弛变量β
(1)
＝1，β增加的步长δ，最优解的集合ω；将固定的β(i)代入到第2层优化问题中求解得到最优解，此最优解即为固定β(i)对应的将并入集合ω中；根据步长δ增加β(i)的大小后重复上述步骤，直到β(i)的数值超过其取值范围为止；从集合ω中获取一个最优值β
opt
使得的取值在集合内最大，即为系统最大保密速率。
[0013]
进一步，步骤(5)所述求解第2层优化问题交替迭代优化算法具体为：第一次迭代时设置初始迭代次数m＝0，φ
(0)
，误差容忍度ζ；固定φ＝φ
(m-1)
，求解第一个子问题，得到w
1(m)
与w
2(m)
；给定w1＝w
1(m)
与求解第二个子问题，得到φ
(m)
；将本次迭代得到的系统保密速率与上一次迭代得到的系统保密速率进行比较，判断迭代是否收敛，若未收敛则进行下一轮的迭代，否则结束迭代。
[0014]
与现有的不完美csi场景下irs辅助无线通信系统安全方案设计的相关研究相比，本发明具有以下有益技术效果：(1)当优化问题的目标为最大化系统保密速率时，目标函数和约束复杂且难以求解，所以文献多将系统的某个安全性能指标作为约束，而以合法用户速率最大化、基站发送功率最小化等作为优化目标。系统保密速率是评价物理层安全传输性能最直接、最重要的性能指标，以保密速率最大化为优化目标，理论和实际价值更高；(2)文献[dong limeng,wang huiming and xiao haitao.secure cognitive radio communication via intelligent reflecting surface[j].ieee transactions on communications,2021,69(7):4678-4690.doi:10.1109/tcomm.2021.3073028.]在窃听信道csi部分已知的情况下，没有使用人工噪声来增强保密性能，本发明在同样的csi条件下，采用了人工噪声，能进一步提升系统的保密速率，但优化问题的复杂度也更高。仿真实验表明，相较于基准方案，本发明所提出的优化算法能有效提高系统的保密速率。
附图说明
[0015]
图1为本发明的通信系统模型；
[0016]
图2为本发明的通信系统仿真模型；
[0017]
图3为基站发射功率对系统保密速率的影响；
[0018]
图4为基站发射天线数目对系统保密速率的影响；
[0019]
图5为irs反射单元数目对系统保密速率的影响。
具体实施方式
[0020]
本发明研究的下行多用户系统模型如图1所示。系统由基站(alice)、irs(rose)和k个用户组成，alice采用时分的方式轮流向这k个用户发送信息，发送给各用户的信息需要相互保密。不失一般性，将当前时刻信息的目标用户称为bob，其他用户则看作窃听者(eves)。alice配置m根天线，rose包含n个反射单元，用户则配置单根天线。alice与用户之间有障碍物遮挡，只能通过rose的反射形成传输链路。假设所有信道均为准静态平坦衰落信道。alice与rose、rose与bob、rose与eves之间的信道系数矩阵或矢量分别记为信道。alice与rose、rose与bob、rose与eves之间的信道系数矩阵或矢量分别记为其中rose的相移矩阵为其中an∈[0,1]、θn∈[0,2π]分别表示rose第n个反射单元的反射幅度与反射相位，n＝1,2,
…
,n。
[0021]
alice轮流向各用户传输信息，在信息传输开始前，alice会向用户发送导频序列，用户根据接收到的导频序列进行信道估计，并将结果发送回alice。因此，假设alice可以获得其与本次传输的目标用户，即bob间信道准确的csi。alice所拥有的rose与其他用户，也就是可能的窃听者间信道的csi，由于信道估计时间较早，而信道为时变衰落信道，与当前实际信道的csi会存在误差，并且估计时间越早，误差越大。因此假设在向bob传输信息的过程中，alice拥有alice-rose，rose-bob信道的完美csi，而拥有的rose-eves信道csi则存在误差。采用有界csi误差模型刻画信道csi误差，即真实的信道系数矢量为拥有的估计信道系数矢量与一个随机误差之和：
[0022][0023]
其中，h
re,k
表示irs与第k个窃听者间信道系数矢量，表示irs与第k个窃听者间信道系数矢量，分别表示alice获得的rose与第k个窃听者间信道过时的信道系数矢量，及其与当前真实信道系数矢量间的误差。εk表征alice对第k个窃听者的csi的不确定度，间隔上次通信时间越长，相应的不确定度越大。
[0024]
为了提高信息传输安全性能，alice在向bob发送信息的同时发送人工噪声干扰窃听者的窃听，alice的发送信号表示为
[0025]
x＝w1s+w2a
[0026]
其中，分别表示信息信号与人工噪声；分别为信息信号、人工噪声的波束成形矢量，满足功率约束tr(w1w
1h
)+tr(w2w
2h
)≤p，这里p是alice的最大发射功率。bob和第k个窃听者处的接收信号可以分别表示为
[0027][0028][0029]
其中，分别表示bob、第k个窃听者处的加性复高斯白噪声。将alice的发送信号代入上式并展开，可得
[0030]
[0031][0032]
根据上式，bob和第k个窃听者处的信息传输速率分别为
[0033][0034][0035]
分别表示bob以及第k个窃听者处的噪声功率。
[0036]
假设窃听者之间各自独立获取信息，系统保密速率为bob处的信息传输速率与所有窃听者的信息传输速率差值的最小值：
[0037][0038]
由于窃听信道的csi不完美已知，与真实的h
re,k
间有一个随机误差，其中基站和irs依据设置波束成形矢量和相移矩阵，相比较根据准确的csi进行设置，系统保密速率会有所下降，下降的程度与随机误差有关。以导致性能下降最严重的误差情况下的系统保密速率最大化为目标，联合优化alice信息信号波束成形矢量w1、人工噪声波束成形矢量w2与irs的相移矩阵φ，优化问题为
[0039][0040]
s.t.c1:tr(w1w
1h
+w2w
2h
)≤p,
[0041]
c2:|φ
n,n
|＝1,n＝1,2,
…
,n.
[0042]
问题中，表示在随机误差情况下，alice-rose-eves间的最大信息传输速率。其中，表示真实rose-eves信道系数取值范围；两个约束条件分别是alice的最大发射功率约束和irs反射幅度为1的约束。φ
n,n
表示矩阵φ中第n行第n列的值，tr(x)表示矩阵x的迹，上标h表示矩阵的共轭转置。
[0043]
优化问题中的目标函数与约束均为非凸函数，多个优化变量相互耦合，而且可能的h
re,k
有无限多个，该三个变量的优化问题直接求解十分困难，需要先对问题进行转化。首先引入松弛变量，使目标函数中取最大值的部分转化为等价的无穷多个不等式约束，然后将这无穷多个约束变换为有限个不等式约束，再将转化后的问题变换为两层优化问题。而第1层优化问题的求解中包含第2层优化问题的求解，其中第1层优化问题为单变量优化问题，利用一维搜索算法求解；第2层优化问题包含三个优化变量，利用算法将其分解为两个交替迭代优化的非凸子问题。两个非凸子问题通过利用charnes-cooper变换法和惩罚函数法转化为凸问题进行求解。
[0044]
优化问题的目标函数可写为
[0045][0046]
作代换变量和其中上式可化为
[0047][0048]
为了简化目标函数，引入松弛变量β，问题变换为等价问题
[0049][0050]
s.t.c1:tr(w1w
1h
+w2w
2h
)≤p,
[0051]
c2:|φ
n,n
|＝1,n＝1,2,...,n,
[0052][0053][0054][0055]
其中，log2β表示窃听信道信息速率在所有情况下的最大值。分别表示基站发送的信息信号功率和人工噪声功率均为正值。
[0056]
问题中的目标函数中的取最大值部分转化为了上式中的无穷多个不等式约束c3，需要将其转化为有限个不等式约束。将上述问题的c3改写为
[0057][0058]
将窃听信道误差代入到上式中，可得
[0059][0060][0061]
将上式展开，再合并同类项，可得
[0062][0063][0064]
式中仍然包含无穷多个二次不等式，下面利用文献中给出的方法把这无穷多个不等式转化为有限个不等式。
[0065]
定义f(x)＝xhax+xhb+bhx+c且下式等价成立
[0066][0067]
进行如下变量代换：进行如下变量代换：x＝δh
re,k
，in×n表示n
×
n维单位矩阵。再根据上式，就可将无穷多个约束转化为如下的有限个约束：
[0068][0069]
其中，tk(w1,w2,φ,β,tk)为)为tk为引入的辅助变量。
[0070]
将上式代入到优化问题中，可将其转化为
[0071]
[0072]
s.t.c1:tr(w1+w2)≤p,
[0073]
c2:|φ
n,n
|＝1,n＝1,2,...,n,
[0074][0075][0076]
虽然已将约束条件数量从无限转为有限，但优化问题中仍含有多个相互耦合的变量，直接求解该优化问题仍然十分困难。为了更加高效地解决优化问题，将β分离出来，使优化问题转换为一个等价的两层优化问题：第1层为在β的取值范围内寻找使得目标函数最大的β；第2层优化问题为在β给定的情况下求解使目标函数最大的w1、w2以及φ；在求解第1层问题时包含第2层优化问题的求解。下面先确定β的取值范围。为去掉优化优化问题目标函数的上标+，应保证合法信道速率不低于窃听信道速率，也即应满足
[0077][0078]
由于可进一步得到
[0079][0080]
又由于传输速率的非负性和tr(w1)≤p的约束，可将上式简化为
[0081][0082]
||x||f表示矩阵x的frobenius范数。
[0083]
在求得β取值范围后，优化问题等价的两层优化问题中的第1层优化问题为
[0084][0085][0086]
其中，为给定β下，的最大值，而如何获得该最大值，就是与优化问题等价的两层优化问题中的第2层优化问题：
[0087][0088]
s.t.c1:tr(w1+w2)≤p,
[0089]
c2:|φ
n,n
|＝1,n＝1,2,
…
,n,
[0090][0091]
[0092]
这是一个关于优化变量w1、w2以及φ为sdp问题。该问题的求解将在下文介绍。
[0093]
优化问题为关于β的单变量优化问题，β的取值范围有限，可以使用一维搜索算法进行求解。优化问题的两层求解算法总结如下表1所示。其中，δ为更新β增加的步长，ω为给定β后求解第2层优化问题得到最优解的集合，最后集合中最大值为原优化问题的解。
[0094]
算法1两层优化问题的求解算法
[0095]
(1)初始化参数：步长δ，i＝1，β
(1)
＝1.
[0096]
(2)循环：
[0097]
(3)给定β(i)，求解第2层优化问题，得到解
[0098]
(4)
[0099]
(5)更新β(i)＝β
(i-1)
+δ.
[0100]
(6)直到
[0101]
(7)输出ω中最大值及其对应的β
opt
.
[0102]
第2层优化问题中的优化变量w1、w2与φ仍然相互耦合，联合求解十分困难，因此采用交替迭代优化的方法：(1)固定irs相移矩阵φ，利用charnes-cooper变换等方法优化w1、w2；(2)固定w1、w2，利用惩罚函数和charnes-cooper变换等方法优化相移矩阵φ。两个优化交替迭代进行，直至收敛。
[0103]
在irs相移矩阵φ固定的情况下，第2层优化问题退化为
[0104][0105]
s.t.c1:tr(w1+w2)≤p,
[0106][0107][0108]
该优化问题中的目标函数仍为非凸函数。先进行如下的变量代换：该优化问题中的目标函数仍为非凸函数。先进行如下的变量代换：ξ＞0，其中ξ为引入的辅助变量，应用charnes-cooper变换，将问题变换为
[0109][0110][0111]
c2:tr(q+z)≤pξ,
[0112]
[0113][0114]
其中，定义其中，定义为z，q，λ
w,k
分别为w1，w2以及tk经过变量代换后的结果。
[0115]
上述问题已是关于q、z、λ
w,k
、ξ的凸问题，可以使用cvx工具箱进行求解。记该问题的最优解为q
opt
、z
opt
、ξ
opt
，则上述问题的解为，则上述问题的解为对w
1opt
与进行特征值分解，非零特征值对应的特征向量就是信息信号波束成形矢量与人工噪声矢量
[0116]
在w1、w2给定后，第2层优化问题退化为
[0117][0118]
s.t.c1:|φ
n,n
|＝1,n＝1,2,
…
,n,
[0119][0120]
该问题的目标函数与约束不是凸函数，需要进行一定的形式变换。首先，将相移矩阵φ的对角线元素组织为一个矢量其中，利用矢量v，可得到
[0121][0122]
另定义v＝vvh，对g作奇异值分解，得到
[0123][0124]
其中ai，bi分别表示对矩阵g进行奇异值分解后的列向量与行向量。为了将问题变换为关于v的sdp问题，利用v、v的定义和g的奇异值分解，可得
[0125][0126]
tk(w1,w2,φ,β,tk)定义中的用代换，上述问题的约束c2已经变换为关于v的形式。下面进一步将上述问题的目标函数也变换为关于v的函数。先利用和将上述问题的目标函数展开，将变量φ分离出来，得到
[0127]
再定义代入，则上述问题就可改写为
[0128][0129][0130]
c2:|v
n,n
|＝1,n＝1,2,
…
,n.
[0131]
其中，为
[0132][0133]
上述问题的目标函数仍为非凸函数，作变量代换ξ＞0，再应用charnes-cooper变换，上述问题变换为
[0134][0135][0136][0137]
c3:|ε
n,n
|＝ξ,n＝1,2,
…
,n,
[0138]
[0139]
其中，为
[0140][0141]
上述问题中的优化目标函数已为凸函数，但约束c3仍为非凸函数。将约束c3改写为如下等价约束
[0142][0143]
其中rank(ε)＝1为非凸函数，需要构造一个等价于rank(ε)＝1的凸约束。
[0144]
对于任意半正定矩阵a，下列不等式成立
[0145]
|i+a|≥1+tr(a)
[0146]
当且仅当rank(a)≤1,等号成立。
[0147]
应用上式，rank(ε)＝1可转化为
[0148][0149]
应用惩罚函数法，将约束作为惩罚加入到优化问题的目标函数中去，目标函数改写为惩罚函数，得到转换后的优化问题为
[0150][0151][0152][0153]
c3:vec(ε)＝ξn,
[0154][0155]
其中，κ为rank(ε)＝1的惩罚因子，当κ足够小时，两个优化问题的最优解相同。但是优化问题的目标函数中log2det(i+ε)关于ε为非凸函数，可将在上一轮迭代中得到的ε的解对其进行一阶泰勒近似而转化为
[0156]
log2det(i+ε)≤(log2e)tr{[(i+ε
(m)
)-1
]
*
(ε-ε
(m)
)}
[0157]
+(log2e)log2det(i+ε
(m)
)
[0158]
将上式带入到优化问题中得到
[0159][0160][0161][0162]
c3:vec(ε)＝ξn,
[0163]
[0164]
其中，ε
(m)
为第2层优化问题迭代求解过程中第m次迭代时优化问题中ε的解(这里假设当前是第m+1次迭代)。至此，将优化问题转换为了凸问题，可用cvx工具包求解得到最优解。记该问题的最优解为ε
opt
、ξ
opt
，则优化问题的解为对v
opt
进行特征值分解，非零特征值对应的特征向量就是v
opt
，再将v
opt
对角化可得相移矩阵φ
opt
。优化问题的求解算法归纳如表2所示。其中φ
(0)
为满足优化问题约束c2的任意对角矩阵，ζ表示误差容忍度，当两次迭代得到的保密速率相对差值的绝对值不大于误差容忍度时，迭代结束。
[0165]
算法2第2层优化问题的求解算法
[0166]
(1)初始化参数：迭代次数m＝0，φ
(0)
，ζ.
[0167]
(2)循环：
[0168]
(3)m＝m+1.
[0169]
(4)固定φ＝φ
(m-1)
，求解第一个子问题问题，得到w
1(m)
与w
2(m)
.
[0170]
(5)给定w1＝w
1(m)
与求解第二个子问题，得到φ
(m)
.
[0171]
(6)对w
1(m)
、进行特征值分解可得w
1(m)
、w
2(m)
.
[0172]
(7)直到
[0173]
(8)输出φ
opt
＝φ
(m)
.
[0174]
下面给出的优化方案的性能进行仿真。无特别说明时仿真中的用户数为k＝4，基站、irs和4个用户的位置如图2所示，单位为米(m)。各节点间信道为莱斯衰落信道，信道衰落包括路径损耗(大尺度衰落)和小尺度衰落。从alice到irs的信道矩阵的模型为
[0175][0176]
式中，为路径损耗，其中l0＝-30db表示参考距离为1米时的路径损耗，α
ar
为路径损耗指数，d
ar
为alice到rose之间的距离；为小尺度衰落部分，其中ρ
ar
为莱斯衰落因子，为非视距传输部分的信道系数矩阵，矩阵中的各元素为0均值、单位方差的复高斯随机变量，表示alice到rose之间的视距传输部分的信道系数矩阵：
[0177][0178]
其中，和表示发射天线的方位角和仰角，和表示irs的方位角和仰
角，与为
[0179][0180][0181]
式中(x
alice
,y
alice
,z
alice
)与(x
rose
,y
rose
,z
rose
)分别表示alice与rose的坐标，上式中为
[0182][0183]
与表示均匀线性阵列的导向矢量
[0184][0185][0186]
其中，d表示两根相邻天线的距离，λ为中心载波的波长，仿真中设定d/λ＝0.5。
[0187]
irs与用户间的信道模型与基站到irs的类似，信道系数矩阵为
[0188][0189][0190]
仿真中，信道模型的参数设置为：α
ar
＝3.5，α
rb
＝α
re,k
＝2.5，莱斯衰落因子为ρ
ar
＝1，ρ
rb
＝ρ
re,k
＝5；信道噪声功率为信道误差容忍度为ε1＝0.1、ε2＝0.15、ε3＝0.2；算法2中误差容忍度为γ＝10-3
，惩罚因子为κ＝5
×
10-6
。
[0191]
仿真中同时给出三种基准方案的结果进行性能对比。基准方案1——无人工噪声方案：即人工噪声矢量w2为0矢量，波束成形矢量和irs的相移矩阵采用与本发明类似的方法优化。基准方案2——irs随机相移方案：irs的相移矩阵随机取值，信号波束成形矢量w1以及人工噪声波束成形矢量w2采用与本发明求解优化问题类似的方法得到。基准方案3——最大比发射(maximum ratio transmission,mrt)：基站信息信号波束成形矢量为alice-rose-bob级联信道的共轭，同时在级联信道的正交方向发送人工噪声，信号波束成形矢量以及人工噪声协方差矩阵分别为
[0192]
[0193][0194]
其中，其中，与为alice给信息信号矢量和人工噪声分配的功率，表示合法信道的零空间。irs相移矩阵φ中的元素则在[0,2π]间随机选取。以保密速率为目标在总功率约束下优化与本节中各仿真图中给出的保密速率为1000组信道样本下的保密速率的平均值。
[0195]
图3给出了本发明方案与3种对比方案的系统保密速率随着alice发射功率变化的情况。仿真中，alice的天线数目m＝4，rose的反射单元数n＝8。从图3可以看出，本发明的设计方案优于其他基准方案。在alice发射功率较低的时候，无人工噪声方案与本发明方案性能比较接近，这是由于发射功率小的时候，alice为了保证与用户之间的通信，将大部分功率分配给信号，只将很少的功率分配给人工噪声。随着发射功率增加，无人工噪声方案的系统保密速率增速要明显低于其他几种方案，这是由于当alice有足够大的发射总功率时，可以给人工噪声分配更多的功率，干扰系统中的窃听者，降低窃听者的接收性能，从而增大系统的保密速率，说明加入人工噪声对提高系统保密速率有明显效果。irs随机相移方案的系统保密速率与本发明方案有类似的变化趋势，但始终低于本发明方案的系统保密速率。这是因为irs随机相移方案采用的是相位随机取值，没有优化irs相移矩阵，这也说明优化irs相移矩阵可以提高系统的保密速率。可以注意到，最大比发射方案的系统保密速率要明显低于本发明方案和irs随机相移方案，这是由于最大比发射方案中，基站的波束成形不是以系统保密速率最大化为目标进行优化，irs相移矩阵也是随机取值。
[0196]
图4给出系统保密速率随着alice处的发射天线数变化的情况。仿真中，alice发射功率p＝30dbm，rose的反射单元数n＝8。从图4可以看出，所有方案的系统保密速率均随着alice发射天线增多而增大，这是由于随着天线数目的增大，发射机有更大的空间自由度，能更精准地控制信号和噪声波束。可以注意到，随着发射天线的增加，最大比发射方案的保密速率的增长速度最小，这是由于最大比发射方案选用的mrt波束不是以系统保密速率最大化为目标，不能很好地利用alice端发射天线数增加带来的信道增益。
[0197]
图5给出系统保密速率随着rose处反射单元数变化的情况。仿真中，alice发射功率p＝30dbm，天线数目m＝4。可以看出，四种方案的系统保密速率均随着rose的反射单元数增大而增大。其中，无人工噪声方案的保密速率在反射单元数目增加到8以后的增长很小，这是由于相比于其他方案没有加入人工噪声，没有干扰系统中的窃听者，降低窃听者的接收性能，而且随着反射单元数目超过发射天线数目，irs不能有效利用反射单元增加所带来的空间自由度与信道增益，导致系统保密速率相比于其他三种方案提升很小。
技术特征：1.不准确信道状态信息条件下智能反射表面辅助的物理层安全优化方法，其特征在于，包括以下步骤：(1)构建通信系统模型，根据构建的系统模型得到合法用户与窃听用户处的信息传输速率，进而得到系统保密速率；(2)以系统保密速率最大化为优化目标，构造在基站总发射功率和智能反射表面反射单元单位模约束下，对基站发射信息信号和人工噪声波束成形矢量，以及智能反射表面的相移矩阵进行联合优化的数学模型；(3)步骤(2)所述数学模型为非凸的优化问题，将其转化为求解两层优化问题，其中第1层优化问题的求解中包含第2层优化问题的求解；(4)第1层优化问题为单变量优化问题，利用一维搜索算法求解；(5)将第2层优化问题分解为两个交替迭代优化的非凸子问题分别求解。2.根据权利要求1所述不准确信道状态信息条件下智能反射表面辅助的物理层安全优化方法，其特征在于：步骤(1)所述系统模型由基站、irs和k个用户组成，基站采用时分的方式轮流向这k个用户发送信息，发送给各用户的信息需要相互保密，基站配置m根天线，irs包含n个反射单元，用户配置单根天线，基站与用户之间有障碍物遮挡，只能通过irs的反射形成传输链路，所有信道均为准静态平坦衰落信道，将当前时刻基站发送信息的目标用户定义为合法用户，其他用户为窃听者。3.根据权利要求1所述不准确信道状态信息条件下智能反射表面辅助的物理层安全优化方法，其特征在于：基站在向合法用户传输信息的过程中，拥有基站-irs和irs-合法用户信道的完美csi，irs-窃听用户信道csi则存在误差；采用有界csi误差模型刻画信道csi误差，即真实的信道系数矢量为拥有的估计信道系数矢量与一个随机误差之和：其中，h
re,k
、δh
re,k
分别表示irs与第k个窃听者间信道系数矢量、基站获得的irs与第k个窃听者间信道过时的信道系数矢量，及其与当前真实信道系数矢量间的误差，ε
k
表征基站对第k个窃听者的csi的不确定度，间隔上次通信时间越长，相应的不确定度越大。4.根据权利要求1或3所述不准确信道状态信息条件下智能反射表面辅助的物理层安全优化方法，其特征在于：步骤(1)所述合法用户与窃听用户处的信息传输速率通过以下方法计算：合法用户和第k个窃听者处的接收信号分别表示为法计算：合法用户和第k个窃听者处的接收信号分别表示为式中，基站与irs、irs与合法用户、irs与窃听者之间的信道系数矩阵或矢量分别记为h
ar
、其中irs的相移矩阵为其中a
n
∈[0,1]、θ
n
∈[0,2π]分别表示irs第n个反射单元的反射幅度与反射相位，n＝1,2,
…
,n；n
b
、n
e,k
分别表示合法用户、第k个窃听者处的加
性复高斯白噪声；x为基站的发送信号，表示为x＝w1s+w2a其中，s,a分别表示信息信号与人工噪声；w1、w2分别为信息信号、人工噪声的波束成形矢量；合法用户和第k个窃听者处的信息传输速率分别为速率分别为所述系统保密速率为合法用户处的信息传输速率与所有窃听者的信息传输速率差值的最小值：5.根据权利要求4所述不准确信道状态信息条件下智能反射表面辅助的物理层安全优化方法，其特征在于：步骤(2)所述优化数学模型在基站总发射功率和智能反射表面反射单元单位模约束下，以系统保密速率最大化为优化目标，优化问题构造为s.t.c1:tr(w1w
1h
+w2w
2h
)≤p,c2:|φ
n,n
|＝1,n＝1,2,
…
,n.其中，表示在随机误差情况下，基站-irs-窃听者之间的最大信息传输速率，其中，表示真实irs-窃听者信道系数取值范围；两个约束条件分别是基站的最大发射功率约束和irs反射幅度为1的约束；φ
n,n
表示矩阵φ中第n行第n列的值，tr(x)表示矩阵x的迹，p表示基站的最大发射功率。6.根据权利要求5所述不准确信道状态信息条件下智能反射表面辅助的物理层安全优化方法，其特征在于：步骤(3)所述将非凸的优化问题转化为两层优化问题，具体包括：引入松弛变量，使目标函数中取最大值的部分转化为等价的无穷多个不等式约束，然后将这无穷多个约束变换为有限个不等式约束，再将转化后的问题变换为两层优化问题，两层优化问题中的第1层优化问题为：问题中的第1层优化问题为：其中，β为引入的松弛变量，为给定β下，
的最大值，而如何获得该最大值，就是两层优化问题中的第2层优化问题：s.t.c1:tr(w1+w2)≤p,c2:|φ
n,n
|＝1,n＝1,2,
…
,n,c3:c4:式中，和t
k
为引入的变量，t
k
(w1,w2,φ,β,t
k
)为i
n
×
n
表示n
×
n维单位矩阵。7.根据权利要求6所述不准确信道状态信息条件下智能反射表面辅助的物理层安全方案的设计方法，其特征在于：步骤(3)所述求解两层优化问题具体方法为：设置初始迭代次数i＝1，松弛变量β
(1)
＝1，β增加的步长δ，最优解的集合ω；将固定的β
(i)
代入到第2层优化问题中求解得到最优解，此最优解即为固定β
(i)
对应的将并入集合ω中；根据步长δ增加β
(i)
的大小后重复上述步骤，直到β
(i)
的数值超过其取值范围为止；从集合ω中获取一个最优值β
opt
使得的取值在集合内最大，即为系统最大保密速率。8.根据权利要求6所述不准确信道状态信息条件下智能反射表面辅助的物理层安全方案的设计方法，其特征在于：步骤(5)所述将第2层优化问题分解为两个交替迭代优化的非凸子问题，在irs相移矩阵φ固定的情况下，第一个子问题为s.t.c1:tr(w1+w2)≤p,c2:c3:在w1、w2给定后，第二个子问题为
9.根据权利要求8所述不准确信道状态信息条件下智能反射表面辅助的物理层安全方案的设计方法，其特征在于：步骤(5)所述求解第2层优化问题交替迭代优化算法具体为：第一次迭代时设置初始迭代次数m＝0，φ
(0)
，误差容忍度ζ；固定φ＝φ
(m-1)
，求解第一个子问题，得到w
1(m)
与w
2(m)
；给定w1＝w
1(m)
与求解第二个子问题，得到φ
(m)
；将本次迭代得到的系统保密速率与上一次迭代得到的系统保密速率进行比较，判断迭代是否收敛，若未收敛则进行下一轮的迭代，否则结束迭代。
技术总结本发明公开了一种不准确信道状态信息条件下智能反射表面辅助的物理层安全优化方法，利用智能反射表面对基站发送信号进行反射，从而实现信号从基站到用户的传输过程，再针对多用户系统物理层安全性能的优化，在基站拥有的窃听信道状态信息存在误差的条件下，以最坏情况下的系统保密速率最大化为目标，对基站发射信息信号和人工噪声的波束成形矢量、智能反射表面的相移矩阵进行联合优化。利用松弛变量、惩罚函数、Charnes-Cooper变换和交替迭代优化方法将原非凸半正定规划问题转化为凸优化问题，再采用标准的凸优化问题的求解方法进行求解。解。解。
技术研发人员：雷维嘉 翟泽旭
受保护的技术使用者：重庆邮电大学
技术研发日：2022.06.20
技术公布日：2022/11/1
