CN115378652A - 一种ris辅助的wpcn系统物理层安全通信方法 - Google Patents

一种ris辅助的wpcn系统物理层安全通信方法 Download PDF

Info

Publication number
CN115378652A
CN115378652A CN202210868899.3A CN202210868899A CN115378652A CN 115378652 A CN115378652 A CN 115378652A CN 202210868899 A CN202210868899 A CN 202210868899A CN 115378652 A CN115378652 A CN 115378652A
Authority
CN
China
Prior art keywords
ris
base station
user
factor
representing
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
CN202210868899.3A
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.)
South China University of Technology SCUT
Original Assignee
South China University of Technology SCUT
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 South China University of Technology SCUT filed Critical South China University of Technology SCUT
Priority to CN202210868899.3A priority Critical patent/CN115378652A/zh
Publication of CN115378652A publication Critical patent/CN115378652A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L63/00—Network architectures or network communication protocols for network security
    • H04L63/20—Network architectures or network communication protocols for network security for managing network security; network security policies in general
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L69/00—Network arrangements, protocols or services independent of the application payload and not provided for in the other groups of this subclass
    • H04L69/30—Definitions, standards or architectural aspects of layered protocol stacks
    • H04L69/32—Architecture of open systems interconnection [OSI] 7-layer type protocol stacks, e.g. the interfaces between the data link level and the physical level
    • H04L69/322—Intralayer communication protocols among peer entities or protocol data unit [PDU] definitions
    • H04L69/323—Intralayer communication protocols among peer entities or protocol data unit [PDU] definitions in the physical layer [OSI layer 1]
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04W—WIRELESS COMMUNICATION NETWORKS
    • H04W28/00—Network traffic management; Network resource management
    • H04W28/02—Traffic management, e.g. flow control or congestion control
    • H04W28/08—Load balancing or load distribution
    • H04W28/09—Management thereof
    • H04W28/0917—Management thereof based on the energy state of entities

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Security & Cryptography (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Computer Hardware Design (AREA)
  • Computing Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Mobile Radio Communication Systems (AREA)

Abstract

本发明公开了一种RIS辅助的WPCN系统物理层安全通信方法，该方法通过采用RIS辅助和WPCN无线供能技术，可实现WPCN系统的安全通信和终端的无限能量补充。该通信方法包括以下步骤：S1:构建RIS辅助的WPCN系统网络模型；S2:确定模型优化的原问题，在充分考虑非线性能量收集模型的情形下，建立以最大化系统安全传输速率为目标的优化问题；S3:简化所述原问题，通过采用AO、SCA和SDR方法将原问题转化为若干凸子问题；S4:建立迭代算法，通过所述迭代算法对所述若干子问题进行交替迭代求解从而得到最优基站波束赋形因子、RIS能量与信息反射系数、以及时间分配因子，实现系统安全传输速率最大化。

Description

一种RIS辅助的WPCN系统物理层安全通信方法
技术领域
本发明涉及无线通信技术领域，具体涉及一种适用于无线供能通信网络(Wireless Powered Communication Networks,WPCN)系统的物理层安全通信方法。
背景技术
物理层安全通信利用无线信道的随机性来实现保密通信，与传统加密方式相比，具有更高的安全性。WPCN利用收集技术提取接收信号中的能量，可以为无线终端节点的能量受限问题提供便利的解决方法。可重构智能表面(Reconfigurable IntelligentSurface,RIS)由大量低成本的无源元件组成的一个无源反射表面阵列，可对无线电波进行相位、幅度和频率等的调整，使网络达到性能最佳化、损耗最低化的目的。将RIS应用于物理层安全，可望进一步提高无线网络的安全性。
目前已有一些学者对WPCN物理层安全进行一定的研究。胡凯等申请的专利《一种基于智能反射面WPCN系统的物理层安全通信方法》提出了RIS辅助的WPCN系统物理层安全通信方法，但其未考虑能量收集电路本身的能量消耗，并且考虑线性能量收集模型，与实际的能量收集电路不符。Zhong Y在《IEEE Communications Letters》发表的论文“Cooperative Jamming-Aided Secure Wireless Powered Communication Networks:AGame Theoretical Formulation”中，提出利用干扰节点协作WPCN系统实现安全传输的方法，其不足之处在于该方案同样考虑线性能量收集模型，并且在上行链路用户向基站传输信息时，干扰源生成的人工噪声对窃听端和用户端都会造成干扰，从而影响系统的安全传输速率。
发明内容
本发明的目的是为了解决现有技术中的上述缺陷，提供一种RIS辅助的WPCN系统物理层安全通信方法。在考虑电路硬件损耗和非线性能量收集模型的情形下，以系统的安全传输速率最大化为目标，对基站波束赋形因子、RIS能量与信息反射系数以及时间分配因子进行联合优化，以此达到在为终端节点提供无线充电功能的同时，提高合法终端的信息传输速率、降低窃听端的窃听能力，从而实现安全通信的目的。
本发明的目的可以通过采取如下技术方案达到：
一种RIS辅助的WPCN系统物理层安全通信方法，该WPCN系统包括1个配备N根发射天线构成均匀线形阵列的基站、1个配备M个反射单元构成均匀矩形阵列的控制器RIS、1个配备单天线的合法用户User和1个配备单个天线窃听端Eve，所述物理层安全通信方法的实现步骤如下：
S1、在传输开始前，基站向包括用户User和控制器RIS在内的所有接收端发送广播指令，所有接收端收到广播指令后向基站发送训练信号；
S2、基站根据接收的训练符号对WPCN系统上行合法信道进行估计，得到基站与控制器RIS的信道信息
Figure BDA0003760489260000021
基站与用户User的信道信息
Figure BDA0003760489260000022
其中
Figure BDA0003760489260000023
表示基站的第i根天线与控制器RIS之间的信道信息，hi表示基站的第i根天线与User之间的信道信息，
Figure BDA0003760489260000024
表示复数域，通过信道信息的对称性获取基站与用户User的下行合法信道；
S3、定义最大化系统安全传输速率的优化问题如下
Figure BDA0003760489260000031
s.t.C1:0≤τk≤1,k＝0,1
C2:τ0+τ1＝1
C3:||f||2≤Pmax
Figure BDA0003760489260000032
Figure BDA0003760489260000033
式中f＝[f1,f2,...,fi...,fN]表示基站波束赋形因子，fi表示基站第i根天线的波束赋形因子，τ＝[τ0,τ1]T表示时间分配因子，τ0表示基站传输能量信号的时间因子，τ1表示用户User传输信息信号的时间因子，e＝[e1,e2,...,em,...,eM]T表示控制器RIS能量反射系数向量，
Figure BDA0003760489260000034
q＝[q1,q2,...,qm,...,qM]T表示控制器RIS信号反射系数向量，
Figure BDA0003760489260000035
R＝[τ1(RU-RE)]+＝maxτ1(RU-RE),0表示WPCN系统的安全传输速率，RU＝log2(1+γU)表示基站接收到用户User的信息速率，
Figure BDA0003760489260000036
表示信噪比，
Figure BDA0003760489260000037
表示上行链路用户User向基站发送信息信号的总信道增益，
Figure BDA0003760489260000038
表示用户User的噪声功率，
Figure BDA0003760489260000039
表示下行链路基站向用户User发送能量信号的总信道增益，μr|hUf|2表示收集电路本身的功率消耗，μr表示电路损耗因子；RE＝log2(1+γE)表示用户User被窃听端Eve窃听到的信息速率，
Figure BDA00037604892600000310
表示窃听端信噪比，
Figure BDA00037604892600000311
表示窃听端Eve窃听用户User发出信号的信道增益，
Figure BDA00037604892600000312
表示窃听端Eve的噪声功率；Pmax表示基站最大发送功率；P＝min|hUf|2,Psat表示用户User接收功率，其中Psat表示用户User最大接收功率，|hUf|2表示基站给用户传输的功率；约束条件C1和C2表示时间分配因子的可能取值；约束条件C3表示基站波束赋形因子f可能的取值；约束条件C4和C5分别表示em和qm可能的取值；
S4、基站生成信号s，该信号s服从均值为0方差为1的高斯分布，在传输时隙t0＝τ0T中，其中T表示一个传输时隙的时间，利用基站波束赋形因子f，对信号s进行波束成形，基站向用户User传输的能量信号为yU＝(f·s)。
进一步地，首先基站发送训练信号进行信道估计继而获取控制器RIS与用户User之间的信道信息
Figure BDA0003760489260000041
控制器RIS与窃听端Eve之间的信道信息
Figure BDA0003760489260000042
通过准静态平稳衰落特性和信道的统计信息计算得到窃听端Eve和用户User之间的信道信息
Figure BDA0003760489260000043
进一步地，所述物理层安全通信方法中，首先利用AO、SCA和SDR方法对步骤S3中所述的非凸优化问题进行分解，将其分解为若干凸子问题；
进一步地，利用凸优化工具包CVX对步骤S3中所述的优化问题进行迭代求解；
最终，获取最优基站波束赋形因子
Figure BDA0003760489260000044
其中fi *表示基站第i根天线的最优波束赋形因子、控制器RIS能量反射系数向量
Figure BDA0003760489260000045
其中
Figure BDA0003760489260000046
表示控制器RIS第m个最优的能量反射因子、RIS信号反射系数向量
Figure BDA0003760489260000047
其中，
Figure BDA0003760489260000048
表示控制器RIS第m个最优的信息反射因子、时间分配因子
Figure BDA0003760489260000049
其中，
Figure BDA00037604892600000410
表示基站传输能量信号的最优时间因子，
Figure BDA00037604892600000411
表示用户User传输信息信号的最优时间因子,过程如下：
S3.1、初始化收敛容忍度ε、迭代变量k、波束赋形因子f＝[f1,f2,...,fi...,fN]、控制器RIS能量反射系数向量e＝[e1,e2,...,em,...,eM]T、控制器RIS信号反射系数向量q＝[q1,q2,...,qm,...,qM]T和时间分配因子τ＝[τ0,τ1]T，获取k＝1时的初始值
Figure BDA0003760489260000051
其中fi (1)表示基站第i根天线的k＝1次迭代值，
Figure BDA0003760489260000052
Figure BDA0003760489260000053
表示控制器RIS第m个能量反射元件的k＝1次迭代值，
Figure BDA0003760489260000054
Figure BDA0003760489260000055
表示控制器RIS第m个信息反射元件的k＝1次迭代值，
Figure BDA0003760489260000056
其中
Figure BDA0003760489260000057
表示基的时间因子站传输能量信号的时间因子k＝1次迭代值，
Figure BDA0003760489260000058
表示用户User传输信息信号k＝1次迭代值，计算出k＝1时用户收集到的功率
Figure BDA0003760489260000059
其中
Figure BDA00037604892600000510
表示k＝1基站向用户User传输的功率，Psat为用户User最大收集功率，从而计算出系统最大安全传输速率R(1)；
S3.2、固定在第k＝1时RIS能量反射系数向量e(k)、RIS信息反射系数向量q(k)和时间分配因子τ(k)，令
Figure BDA00037604892600000511
使用CVX工具包求解步骤S3中最大化系统安全传输速率的优化问题得到在第k+1次波束赋形因子的最优解f(k+1)；
S3.3、利用上述求得的第k+1次波束赋形因子f(k+1)，固定f(k+1),q(k),τ(k)，令
Figure BDA00037604892600000512
使用凸优化工具包CVX求解步骤S3中最大化系统安全传输速率的优化问题得到在第k+1次RIS能量反射系数向量的最优解e(k+1)；
S3.4、利用上述求得的第k+1次波束赋形因子f(k+1)和RIS能量反射系数向量e(k+1)，固定f(k+1),e(k+1),τ(k)，令
Figure BDA00037604892600000513
使用凸优化工具包CVX求解步骤S3中优化问题得到在第k+1次RIS信息反射系数向量的最优解q(k+1)；
S3.5、利用上述求得的第k+1次波束赋形因子f(k+1),RIS能量反射系数向量e(k+1)和RIS信息反射系数向量q(k+1),令
Figure BDA00037604892600000514
使用凸优化工具包CVX求解步骤S3中优化问题得到在第k+1次时间分配因子的最优解τ(k+1)，计算出k+1时系统最大安全传输速率R(k+1)；
S3.6、判断R(k+1)-R(k)≤ε是否成立，若R(k+1)-R(k)≤ε成立，则停止迭代，执行步骤S3.7；否则令τ(k)＝τ(k+1),e(k)＝e(k+1),q(k)＝q(k+1),f(k)＝f(k+1),k＝k+1返回步骤S3.2；
S3.7、通过CVX求解，最终得到WPCN通信系统达到最大化安全速率时对应最优的f*、e*、q*和τ*。
本发明相对于现有技术具有如下的优点及效果：
1、本发明在考虑电路硬件损耗和非线性能量收集模型的情形下，提出一种RIS辅助的WPCN物理层安全通信方法。通过对基站波束赋形因子f、RIS能量反射单元e与信息反射单元q和时间分配因子τ进行联合优化使得系统安全传输速率最大化。
2、本发明利用WPCN技术解决传统网络能量受限问题节点以非线性能量收集模型进行能量收集，利用RIS技术在物理层安全层面提高通信网络的安全性，在此基础上构建了一个非凸优化问题。通过将非凸多变量优化问题拆分成多个含有单一优化变量的凸子问题进行求解，实现简单，复杂性较低。
3、本发明提出的方法，在考虑用户User能量收集受限和非线性能量收集情况的情况，通过对比基站传输功率和用户User最大接收功率的大小从而确定用户User在实际应用中真实收集功率，继而计算出更加符合实际的WPCN系统安全传输速率。
附图说明
此处所说明的附图用来提供对本发明的进一步理解，构成本申请的一部分，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。在附图中：
图1是本发明中RIS辅助的WPCN通信系统模型示意图；
图2是本发明中公开的RIS辅助的WPCN系统物理层安全通信方法的实现流程图；
图3是本发明获取系统安全传输速率的实现流程图；
图4是本发明中公开的RIS辅助的WPCN非线性与线性能量收集模型对比图；
图5是本发明固定基站传输速率为Pmax＝30dBm时，在不同实施例下系统安全传输速率对比图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
实施例1
本实施例公开一种RIS辅助的WPCN系统物理层安全通信方法，该WPCN系统包括1个配备N根发射天线构成均匀线形阵列的基站、1个配备M个反射单元构成均匀矩形阵列的控制器RIS、1个配备单天线的合法用户User和1个配备单个天线窃听端Eve。
本实施例中，具体参数设置如下：
设定WPCN系统在三维坐标系中，基站位于(0m,0m,10m)处、控制器RIS位于(5m,0m,10m)处、用户User位于(5m,2m,10m)和窃听端Eve位于在以(2m,5m,10m)；基站配有N＝6根天线，最大传输功率Pmax＝30dBm，传输时隙T＝1，用户最大接收功率Psat＝10dBm；控制器RIS配备M＝32个反射元，收敛容忍度ε＝10-4、噪声功率
Figure BDA0003760489260000081
本实施例中的所有信道均为莱斯信道，信道模型为
Figure BDA0003760489260000082
其中β表示莱斯因子，β110；GLoS表示视距无线传输部分，GNLoS表示非视距无线传输部分；PL为路径损失，其计算公式为PL＝-30-10αlog10(d)，其中α为路径损失参数，d为链路的欧式距离，与控制器RIS相关的路径损失参数α为2.2，基站到用户User和基站到窃听端Eve的路径损失参数α为3.6。
该物理层安全通信方法实现步骤如下：
S1、在传输开始前，基站向包括用户User和控制器RIS在内的所有接收端发送广播指令，所有接收端收到广播指令后向基站发送训练信号；
S2、基站根据接收的训练符号对WPCN系统上行合法信道进行估计，得到基站与控制器RIS的信道信息
Figure BDA0003760489260000083
基站与用户User的信道信息
Figure BDA0003760489260000084
其中
Figure BDA0003760489260000085
表示基站的第i根天线与控制器RIS之间的信道信息，hi表示基站的第i根天线与User之间的信道信息，
Figure BDA0003760489260000086
表示复数域，通过信道信息的对称性获取基站与用户User的下行合法信道；
S3、定义最大化系统安全传输速率的优化问题如下
Figure BDA0003760489260000087
s.t.C1:0≤τk≤1,k＝0,1
C2:τ0+τ1＝1
C3:||f||2≤Pmax
Figure BDA0003760489260000088
Figure BDA0003760489260000089
其中f＝[f1,f2,...,fn...,f6]是基站波束赋形因子，其中fi表示基站第i根天线的波束赋形因子，τ＝[τ0,τ1]T表示时间分配因子，其中τ0表示基站传输能量信号的时间因子，τ1表示用户User传输信息信号的时间因子，e＝[e1,e2,...,em,...,e32]T是控制器RIS能量反射系数向量，其中
Figure BDA0003760489260000091
q＝[q1,q2,...,qm,...,q32]T是RIS信号反射系数向量，其中
Figure BDA0003760489260000092
R＝[τ1(RU-RE)]+＝maxτ1(RU-RE),0表示WPCN系统的安全传输速率，其中RU＝log2(1+γU)表示基站接收到用户User的信息速率，
Figure BDA0003760489260000093
表示信噪比，
Figure BDA0003760489260000094
表示上行链路用户User向基站发送信息信号的总信道增益，
Figure BDA0003760489260000095
表示用户User的噪声功率，
Figure BDA0003760489260000096
表示下行链路基站向用户User发送能量信号的总信道增益，μr|hUf|2表示收集电路本身的功率消耗，μr表示电路损耗因子；RE＝log2(1+γE)表示用户User被窃听端Eve窃听到的信息速率，
Figure BDA0003760489260000097
表示窃听端信噪比，
Figure BDA0003760489260000098
表示窃听端Eve窃听用户User发出信号的信道增益，
Figure BDA0003760489260000099
表示窃听端Eve的噪声功率；Pmax表示基站最大发送功率；P＝min|hUf|2,Psat表示用户User接收功率，其中Psat表示用户User最大接收功率，|hUf|2表示基站给用户传输的功率；约束条件C1和C2表示时间分配因子的可能取值；约束条件C3表示基站波束赋形因子f可能的取值；约束条件C4和C5分别表示em和qm可能的取值；
S3.1、初始化收敛容忍度ε、迭代变量k、波束赋形因子f＝[f1,f2,...,fi...,f6]、控制器RIS能量反射系数向量e＝[e1,e2,...,em,...,e32]T、控制器RIS信号反射系数向量q＝[q1,q2,...,qm,...,q32]T和时间分配因子τ＝[τ0,τ1]T，获取k＝1时的初始值
Figure BDA0003760489260000101
其中fi (1)表示基站第i根天线的k＝1次迭代值，
Figure BDA0003760489260000102
其中
Figure BDA0003760489260000103
表示控制器RIS第m个能量反射元件的k＝1次迭代值，
Figure BDA0003760489260000104
其中
Figure BDA0003760489260000105
表示控制器RIS第m个信息反射元件的k＝1次迭代值，
Figure BDA0003760489260000106
其中
Figure BDA0003760489260000107
表示基站传输能量信号的时间因子k＝1次迭代值，
Figure BDA0003760489260000108
表示用户User传输信息信号的时间因子k＝1次迭代值，计算出k＝1时用户收集到的功率
Figure BDA0003760489260000109
其中
Figure BDA00037604892600001010
表示k＝1基站向用户User传输的功率，Psat为用户User最大收集功率，从而计算出系统最大安全传输速率R(1)；
S3.2、固定在第k＝1时RIS能量反射系数向量e(k)、RIS信息反射系数向量q(k)和时间分配因子τ(k)，令
Figure BDA00037604892600001011
使用CVX工具包求解步骤S3中最大化系统安全传输速率的优化问题得到在第k+1次波束赋形因子的最优解f(k+1)；
S3.3、利用上述求得的第k+1次波束赋形因子f(k+1)，固定f(k+1),q(k),τ(k)，令
Figure BDA00037604892600001012
使用凸优化工具包CVX求解步骤S3中最大化系统安全传输速率的优化问题得到在第k+1次RIS能量反射系数向量的最优解e(k+1)；
S3.4、利用上述求得的第k+1次波束赋形因子f(k+1)和RIS能量反射系数向量e(k+1)，固定f(k+1),e(k+1),τ(k)，令
Figure BDA00037604892600001013
使用凸优化工具包CVX求解步骤S3中优化问题得到在第k+1次RIS信息反射系数向量的最优解q(k+1)；
S3.5、利用上述求得的第k+1次波束赋形因子f(k+1),RIS能量反射系数向量e(k+1)和RIS信息反射系数向量q(k+1),令
Figure BDA00037604892600001014
使用凸优化工具包CVX求解步骤S3中优化问题得到在第k+1次时间分配因子的最优解τ(k+1)，计算出k+1时系统最大安全传输速率R(k+1)；
S3.6、判断R(k+1)-R(k)≤ε是否成立，若R(k+1)-R(k)≤ε成立，则停止迭代，执行步骤S3.7；否则令τ(k)＝τ(k+1),e(k)＝e(k+1),q(k)＝q(k+1),f(k)＝f(k+1),k＝k+1返回步骤S3.2；
S3.7、通过CVX求解，最终得到WPCN通信系统达到最大化安全速率时对应最优的f*、e*、q*和τ*。
图4是RIS辅助的WPCN非线性与线性能量收集模型对比图；从图中可以看出，在实施例1所表示的曲线中，随着基站的发功率率逐渐增大，WPCN系统的安全传输速率也逐渐增大最后趋于平缓，说明基站的发射功率对WPCN系统的安全传输速率有影响，但是当基站的发射功率超过某一界限WPCN系统的安全传输速率将不受基站功率影响。
实施例2
本实施例公开一种RIS辅助的WPCN系统物理层安全通信方法，该WPCN系统包括1个配备N根发射天线构成均匀线形阵列的基站、1个配备M个反射单元构成均匀矩形阵列的控制器RIS、1个配备单天线的合法用户User和1个配备单个天线窃听端Eve。
本实施例中，具体参数设置如下：
设定WPCN系统在三维坐标系中，基站位于(0m,0m,10m)处、控制器RIS位于(5m,0m,10m)处、用户User位于(5m,2m,10m)和窃听端Eve位于在以(2m,5m,10m)；基站配有N＝6根天线，最大传输功率Pmax＝30dBm，传输时隙T＝1，用户最大接收功率Psat＝10dBm；控制器RIS配备M＝16个反射元，收敛容忍度ε＝10-4、噪声功率
Figure BDA0003760489260000111
本实施例中的所有信道均为莱斯信道，信道模型为
Figure BDA0003760489260000121
其中β表示莱斯因子，β110；GLoS表示视距无线传输部分，GNLoS表示非视距无线传输部分；PL为路径损失，其计算公式为PL＝-30-10αlog10(d)，其中α为路径损失参数，d为链路的欧式距离，与控制器RIS相关的路径损失参数α为2.2，基站到用户User和基站到窃听端Eve的路径损失参数α为3.6。
该物理层安全通信方法实现步骤如下：
S1、在传输开始前，基站向包括用户User和控制器RIS在内的所有接收端发送广播指令，所有接收端收到广播指令后向基站发送训练信号；
S2、基站根据接收的训练符号对WPCN系统上行合法信道进行估计，得到基站与控制器RIS的信道信息
Figure BDA0003760489260000122
基站与用户User的信道信息
Figure BDA0003760489260000123
其中
Figure BDA0003760489260000124
表示基站的第i根天线与控制器RIS之间的信道信息，hi表示基站的第i根天线与User之间的信道信息，
Figure BDA0003760489260000125
表示复数域，通过信道信息的对称性获取基站与用户User的下行合法信道；
S3、定义最大化系统安全传输速率的优化问题如下
Figure BDA0003760489260000126
s.t.C1:0≤τk≤1,k＝0,1
C2:τ0+τ1＝1
C3:||f2≤Pmax
Figure BDA0003760489260000127
Figure BDA0003760489260000128
其中f＝[f1,f2,...,fn...,f6]是基站波束赋形因子，其中fi表示基站第i根天线的波束赋形因子，τ＝[τ0,τ1]T表示时间分配因子，其中τ0表示基站传输能量信号的时间因子，τ1表示用户User传输信息信号的时间因子，e＝[e1,e2,...,em,...,e16]T是控制器RIS能量反射系数向量，其中
Figure BDA0003760489260000131
q＝[q1,q2,...,qm,...,q16]T是RIS信号反射系数向量，其中
Figure BDA0003760489260000132
R＝[τ1(RU-RE)]+＝maxτ1(RU-RE),0表示WPCN系统的安全传输速率，其中RU＝log2(1+γU)表示基站接收到用户User的信息速率，
Figure BDA0003760489260000133
表示信噪比，
Figure BDA0003760489260000134
表示上行链路用户User向基站发送信息信号的总信道增益，
Figure BDA0003760489260000135
表示用户User的噪声功率，
Figure BDA0003760489260000136
表示下行链路基站向用户User发送能量信号的总信道增益，μr|hUf|2表示收集电路本身的功率消耗，μr表示电路损耗因子；RE＝log2(1+γE)表示用户User被窃听端Eve窃听到的信息速率，
Figure BDA0003760489260000137
表示窃听端信噪比，
Figure BDA0003760489260000138
表示窃听端Eve窃听用户User发出信号的信道增益，
Figure BDA0003760489260000139
表示窃听端Eve的噪声功率；Pmax表示基站最大发送功率；P＝min|hUf|2,Psat表示用户User接收功率，其中Psat表示用户User最大接收功率，|hUf|2表示基站给用户传输的功率；约束条件C1和C2表示时间分配因子的可能取值；约束条件C3表示基站波束赋形因子f可能的取值；约束条件C4和C5表示em和qm可能的取值；
S3.1、初始化收敛容忍度ε、迭代变量k、波束赋形因子f＝[f1,f2,...,fi...,f6]、控制器RIS能量反射系数向量e＝[e1,e2,...,em,...,e16]T、控制器RIS信号反射系数向量q＝[q1,q2,...,qm,...,q16]T和时间分配因子τ＝[τ0,τ1]T，获取k＝1时的初始值
Figure BDA00037604892600001310
其中fi (1)表示基站第i根天线的k＝1次迭代值，
Figure BDA00037604892600001311
其中
Figure BDA00037604892600001312
表示控制器RIS第m个能量反射元件的k＝1次迭代值，
Figure BDA00037604892600001313
其中
Figure BDA00037604892600001314
表示控制器RIS第m个信息反射元件的k＝1次迭代值，
Figure BDA0003760489260000141
其中
Figure BDA0003760489260000142
表示基站传输能量信号的时间因子k＝1次迭代值，
Figure BDA0003760489260000143
表示用户User传输信息信号的时间因子k＝1次迭代值，计算出k＝1时用户收集到的功率
Figure BDA0003760489260000144
其中
Figure BDA0003760489260000145
表示k＝1基站向用户User传输的功率，Psat为用户User最大收集功率，从而计算出系统最大安全传输速率R(1)；
S3.2、固定在第k＝1时RIS能量反射系数向量e(k)、RIS信息反射系数向量q(k)和时间分配因子τ(k)，令
Figure BDA0003760489260000146
使用CVX工具包求解步骤S3中最大化系统安全传输速率的优化问题得到在第k+1次波束赋形因子的最优解f(k+1)；
S3.3、利用上述求得的第k+1次波束赋形因子f(k+1)，固定f(k+1),q(k),τ(k)，令
Figure BDA0003760489260000147
使用凸优化工具包CVX求解步骤S3中最大化系统安全传输速率的优化问题得到在第k+1次RIS能量反射系数向量的最优解e(k+1)；
S3.4、利用上述求得的第k+1次波束赋形因子f(k+1)和RIS能量反射系数向量e(k+1)，固定f(k+1),e(k+1),τ(k)，令
Figure BDA0003760489260000148
使用凸优化工具包CVX求解步骤S3中优化问题得到在第k+1次RIS信息反射系数向量的最优解q(k+1)；
S3.5、利用上述求得的第k+1次波束赋形因子f(k+1),RIS能量反射系数向量e(k+1)和RIS信息反射系数向量q(k+1),令
Figure BDA0003760489260000149
使用凸优化工具包CVX求解步骤S3中优化问题得到在第k+1次时间分配因子的最优解τ(k+1)，计算出k+1时系统最大安全传输速率R(k+1)；
S3.6、判断R(k+1)-R(k)≤ε是否成立，若R(k+1)-R(k)≤ε成立，则停止迭代，执行步骤S3.7；否则令τ(k)＝τ(k+1),e(k)＝e(k+1),q(k)＝q(k+1),f(k)＝f(k+1),k＝k+1返回步骤S3.2；
S3.7、通过CVX求解，最终得到WPCN通信系统达到最大化安全速率时对应最优的f*,e*,q*和τ*。
图5是本发明固定基站传输速率为Pmax＝30dBm时，在不同实施例下系统安全传输速率对比图；从图中可以看出当基站的最大发射功率固定时，WPCN系统的安全传输速率随着迭代次数增加而逐渐增加最终趋于平缓，说明本发明所提出方法是可收敛的，从图中可以看出在迭代6次左右，WPCN系统的安全传输速率就已经平缓，说明本发明提出的方法收敛速度是较快的。
上述实施例为本发明较佳的实施方式，但本发明的实施方式并不受上述实施例的限制，其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化，均应为等效的置换方式，都包含在本发明的保护范围之内。

Claims (4)

1.一种RIS辅助的WPCN系统物理层安全通信方法，其特征在于，该WPCN系统包括1个配备N根发射天线构成均匀线形阵列的基站、1个配备M个反射单元构成均匀矩形阵列的控制器RIS、1个配备单天线的合法用户User和1个配备单个天线窃听端Eve，所述物理层安全通信方法的实现步骤如下：
S1、在传输开始前，基站向包括用户User和控制器RIS在内的所有接收端发送广播指令，所有接收端收到广播指令后向基站发送训练信号；
S2、基站根据接收的训练符号对WPCN系统上行合法信道进行估计，得到基站与控制器RIS的信道信息
Figure FDA0003760489250000011
基站与用户User的信道信息
Figure FDA0003760489250000012
其中
Figure FDA0003760489250000013
表示基站的第i根天线与控制器RIS之间的信道信息，hi表示基站的第i根天线与User之间的信道信息，
Figure FDA0003760489250000014
表示复数域，通过信道信息的对称性获取基站与用户User的下行合法信道；
S3、定义最大化系统安全传输速率的优化问题如下
P1:
Figure FDA0003760489250000015
s.t.C1:0≤τk≤1,k＝0,1
C2:τ0+τ1＝1
C3:||f||2≤Pmax
C4:|em|2＝1,
Figure FDA0003760489250000016
C5:|qm|2＝1,
Figure FDA0003760489250000017
式中f＝[f1,f2,...,fi...,fN]表示基站波束赋形因子，fi表示基站第i根天线的波束赋形因子，τ＝[τ0,τ1]T表示时间分配因子，τ0表示基站传输能量信号的时间因子，τ1表示用户User传输信息信号的时间因子，e＝[e1,e2,...,em,...,eM]T表示控制器RIS能量反射系数向量，|em|2＝1,
Figure FDA0003760489250000021
q＝[q1,q2,...,qm,...,qM]T表示控制器RIS信号反射系数向量，|qm|2＝1,
Figure FDA0003760489250000022
R＝[τ1(RU-RE)]+＝maxτ1(RU-RE),0表示WPCN系统的安全传输速率，RU＝log2(1+γU)表示基站接收到用户User的信息速率，
Figure FDA0003760489250000023
表示信噪比，
Figure FDA0003760489250000024
表示上行链路用户User向基站发送信息信号的总信道增益，
Figure FDA0003760489250000025
表示用户User的噪声功率，
Figure FDA0003760489250000026
表示下行链路基站向用户User发送能量信号的总信道增益，μr|hUf|2表示收集电路本身的功率消耗，μr表示电路损耗因子；RE＝log2(1+γE)表示用户User被窃听端Eve窃听到的信息速率，
Figure FDA0003760489250000027
表示窃听端信噪比，
Figure FDA0003760489250000028
Figure FDA0003760489250000029
表示窃听端Eve窃听用户User发出信号的信道增益，
Figure FDA00037604892500000210
表示窃听端Eve的噪声功率；Pmax表示基站最大发送功率；P＝min|hUf|2,Psat表示用户User接收功率，其中Psat表示用户User最大接收功率，|hUf|2表示基站给用户传输的功率；约束条件C1和C2表示时间分配因子的可能取值；约束条件C3表示基站波束赋形因子f可能的取值；约束条件C4和C5分别表示em和qm可能的取值；
S4、基站生成信号s，该信号s服从均值为0方差为1的高斯分布，在传输时隙t0＝τ0T中，其中T表示一个传输时隙的时间，利用基站波束赋形因子f，对信号s进行波束成形，基站向用户User传输的能量信号为yU＝(f·s)。
2.根据权利要求1所述的一种RIS辅助的WPCN系统物理层安全通信方法，其特征在于，所述物理层安全通信方法中，通过发送训练信号进行信道估计的方法继续获取控制器RIS与用户User之间的信道信息
Figure FDA0003760489250000031
控制器RIS与窃听端Eve之间的信道信息
Figure FDA0003760489250000032
通过准静态平稳衰落特性和信道的统计信息计算得到窃听端Eve和用户User之间的信道信息
Figure FDA0003760489250000033
3.根据权利要求1所述的一种RIS辅助的WPCN系统物理层安全通信方法，其特征在于，所述物理层安全通信方法中，利用凸优化工具包CVX对步骤S3中最大化系统安全传输速率的优化问题进行求解，获取最优基站波束赋形因子
Figure FDA0003760489250000034
其中fi *表示基站第i根天线的最优波束赋形因子、控制器RIS能量反射系数向量
Figure FDA0003760489250000035
其中
Figure FDA0003760489250000036
表示控制器RIS第m个最优的能量反射因子、RIS信号反射系数向量
Figure FDA0003760489250000037
其中，
Figure FDA0003760489250000038
表示控制器RIS第m个最优的信息反射因子、时间分配因子
Figure FDA0003760489250000039
其中，
Figure FDA00037604892500000310
表示基站传输能量信号的最优时间因子，
Figure FDA00037604892500000311
表示用户User传输信息信号的最优时间因子。
4.根据权利要求1所述的一种RIS辅助的WPCN系统物理层安全通信方法，其特征在于，所述步骤S3过程如下：
S3.1、初始化收敛容忍度ε、迭代变量k、波束赋形因子f＝[f1,f2,...,fi...,fN]、控制器RIS能量反射系数向量e＝[e1,e2,...,em,...,eM]T、控制器RIS信号反射系数向量q＝[q1,q2,...,qm,...,qM]T和时间分配因子τ＝[τ0,τ1]T，获取k＝1时的初始值
Figure FDA00037604892500000312
其中fi (1)表示基站第i根天线的k＝1次迭代值，
Figure FDA00037604892500000313
Figure FDA00037604892500000314
表示控制器RIS第m个能量反射元件的k＝1次迭代值，
Figure FDA00037604892500000315
Figure FDA00037604892500000316
表示控制器RIS第m个信息反射元件的k＝1次迭代值，
Figure FDA00037604892500000317
其中
Figure FDA00037604892500000318
表示基的时间因子站传输能量信号的时间因子k＝1次迭代值，
Figure FDA00037604892500000319
表示用户User传输信息信号k＝1次迭代值，计算出k＝1时用户收集到的功率
Figure FDA0003760489250000041
其中
Figure FDA0003760489250000042
表示k＝1基站向用户User传输的功率，Psat为用户User最大收集功率，从而计算出系统最大安全传输速率R(1)；
S3.2、固定在第k＝1时RIS能量反射系数向量e(k)、RIS信息反射系数向量q(k)和时间分配因子τ(k)，令
Figure FDA0003760489250000043
使用CVX工具包求解步骤S3中最大化系统安全传输速率的优化问题得到在第k+1次波束赋形因子的最优解f(k+1)；
S3.3、利用上述求得的第k+1次波束赋形因子f(k+1)，固定f(k+1),q(k),τ(k)，令
Figure FDA0003760489250000044
使用凸优化工具包CVX求解步骤S3中最大化系统安全传输速率的优化问题得到在第k+1次RIS能量反射系数向量的最优解e(k+1)；
S3.4、利用上述求得的第k+1次波束赋形因子f(k+1)和RIS能量反射系数向量e(k+1)，固定f(k+1),e(k+1),τ(k)，令
Figure FDA0003760489250000045
使用凸优化工具包CVX求解所述优化问题得到在第k+1次RIS信息反射系数向量的最优解q(k+1)；
S3.5、利用上述求得的第k+1次波束赋形因子f(k+1)，RIS能量反射系数向量e(k+1)和RIS信息反射系数向量q(k+1)，令
Figure FDA0003760489250000046
使用凸优化工具包CVX求解所述优化问题得到在第k+1次时间分配因子的最优解τ(k+1)，计算出k+1时系统最大安全传输速率R(k+1)；
S3.6、判断R(k+1)-R(k)≤ε是否成立，若R(k+1)-R(k)≤ε成立，则停止迭代，执行步骤S3.7；否则令τ(k)＝τ(k+1),e(k)＝e(k+1),q(k)＝q(k+1),f(k)＝f(k+1)，k＝k+1返回步骤S3.2；
S3.7、通过CVX求解，最终得到WPCN通信系统达到最大化安全速率时对应最优的f*，e*，q*和τ*。
CN202210868899.3A 2022-07-22 2022-07-22 一种ris辅助的wpcn系统物理层安全通信方法 Pending CN115378652A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210868899.3A CN115378652A (zh) 2022-07-22 2022-07-22 一种ris辅助的wpcn系统物理层安全通信方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210868899.3A CN115378652A (zh) 2022-07-22 2022-07-22 一种ris辅助的wpcn系统物理层安全通信方法

Publications (1)

Publication Number Publication Date
CN115378652A true CN115378652A (zh) 2022-11-22

Family

ID=84062101

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210868899.3A Pending CN115378652A (zh) 2022-07-22 2022-07-22 一种ris辅助的wpcn系统物理层安全通信方法

Country Status (1)

Country Link
CN (1) CN115378652A (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116489654A (zh) * 2023-05-18 2023-07-25 北京航空航天大学 智能超表面辅助的对抗导频欺骗攻击的物理层安全通信方法
CN117499962A (zh) * 2023-12-29 2024-02-02 南京邮电大学 一种有源可重构智能表面辅助通信的安全能效优化方法
CN118075852A (zh) * 2024-04-18 2024-05-24 南京邮电大学 一种双有源ris辅助的pls传输系统功耗优化方法及系统

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116489654A (zh) * 2023-05-18 2023-07-25 北京航空航天大学 智能超表面辅助的对抗导频欺骗攻击的物理层安全通信方法
CN117499962A (zh) * 2023-12-29 2024-02-02 南京邮电大学 一种有源可重构智能表面辅助通信的安全能效优化方法
CN117499962B (zh) * 2023-12-29 2024-03-15 南京邮电大学 一种有源可重构智能表面辅助通信的安全能效优化方法
CN118075852A (zh) * 2024-04-18 2024-05-24 南京邮电大学 一种双有源ris辅助的pls传输系统功耗优化方法及系统

Similar Documents

Publication Publication Date Title
CN110830097B (zh) 一种基于反射面的主被动互惠共生传输通信系统
CN115378652A (zh) 一种ris辅助的wpcn系统物理层安全通信方法
Mirza et al. Channel estimation method and phase shift design for reconfigurable intelligent surface assisted MIMO networks
CN111818533B (zh) 一种基于智能反射面的无线通信系统设计方法
CN109981153B (zh) 一种人工噪声辅助的大规模mimo安全统计预编码方法
CN112565121B (zh) 用于智能反射面辅助的大规模mimo系统的半盲信道估计方法
CN111787558A (zh) 一种mimome系统及其安全速率优化方法
CN113556164B (zh) Irs辅助的swipt系统中基于能效优先的波束成型优化方法
CN113825159B (zh) 基于智能反射面的无线携能通信系统鲁棒资源分配方法
CN112994770B (zh) 基于部分csi的ris辅助多用户下行鲁棒无线传输方法
CN116015374B (zh) 智能反射面辅助的非正交多址接入携能网络安全波束赋形方法
CN114039706B (zh) 一种基于新型可重构智能表面的空时编码发送方法
CN116436502A (zh) 一种智能反射面选择与相位矩阵调整方法
US20230308200A1 (en) Wireless telecommunications network
CN112888059B (zh) 一种智能反射面辅助通信的快速数据传输协议
CN114626229A (zh) 智能反射面辅助的无线通信与能量收集系统性能分析方法
CN110690913A (zh) 一种基于不完全信道信息的协作空间调制系统中功率分配方法
CN115022146B (zh) 一种可重构智能表面辅助的通信与定位一体化全双工系统
CN115334524B (zh) 一种基于全向智能超表面的通信和雷达目标检测方法
CN116367192A (zh) 一种基于irs辅助的多天线通信系统鲁棒能效优化方法
CN114157392B (zh) 一种分布式irs辅助通信系统安全传输的优化方法
CN113766492B (zh) 一种安全Massive MIMO网络资源配置方法
CN113726396B (zh) 一种全双工无线携能中继通信系统高能效保密传输方法
CN115696437A (zh) 基于irs的无线能量传输网络的用户总速率最大化方法
CN114124185A (zh) 一种优化irs辅助通信系统中相移矩阵的低复杂度方法

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