具体实施方式

下面结合附图对本发明作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案，而不能以此来限制本发明的保护范围。

实施例一：

请参阅图2，本发明系统模型包括：源节点(S)，目的节点(D)，全双工中继节点(FD-R)，和K个不可信能量采集节点(EHR_k)，其中S，D和EHR_k均配备单天线，FD-R配备M_r根接收天线和M_t根发送天线，且假设S和D之间不存在直达径。结合附图1-3对本发明公开的一种全双工无线携能中继通信系统高能效保密传输方法具体步骤描述如下：

构建不可信能量采集节点场景下的全双工无线携能中继节点信号传输模型；

以全双工无线携能中继节点信号传输模型的保密能量效率最大化为目标函数，建立基于人工噪声辅助的波束成形优化问题模型；

采用连续凸近似和惩罚函数联合方法实现对基于人工噪声辅助的波束成形优化问题模型的近似凸转换；

基于近似凸转换后的波束成形优化问题模型，设计一种快速收敛的迭代优化算法，获得使系统保密能量效率最大化的最佳中继节点波束成形矩阵和人工噪声协方差矩阵。

本专利提出的一种人工噪声辅助的全双工无线携能中继通信系统高能效保密传输方法，针对无线携能放大转发全双工中继通信系统，提出该系统中保密能效的定义。中继在采用放大转发协议下，则目的节点和第k个能量采集节点的可达速率分别为：

其中，P_s为源节点的发射功率，W为中继波束成形矩阵，Q为人工噪声协方差矩阵，h_sr是源节点到中继节点的信道矢量，h_rd和h_rk分别是中继节点到目的节点和能量采集节点的信道矢量，和分别为中继、目的节点和第k个能量采集节点的噪声功率。

定义该全双工无线携能中继节点信号保密能量效率η_SEE为：

其中，R_sec表示全双工无线携能中继通信系统的保密传输率，P_total是全双工无线携能中继通信系统的总功率，P_r是中继发射功率，P_c为硬件电路功耗，P_SIC为自干扰消除功耗，β为中继端的功率有效转换系数。

对于能量受限的中继系统，信号的总发送功率P_r应满足如下约束：

其中，为中继最大发射功率。同时为了满足能量采集节点的最小能量采集要求，第k个能量采集节点的采集功率P_k应满足以下约束：

其中，η_k表示第k个能量采集节点所需的最小采集能量，ξ表示能量转换效率。

基于此，该全双工无线携能中继通信系统的保密能量效率最大化问题可以建模为非凸问题(6)，包括：

P_k≥η_k，k∈K (6c)

H_rrW＝0，H_rrQ＝0 (6d)

其中，H_rr是中继自干扰信道，约束条件(6d)为中继节点处的自干扰迫零约束。

其次，假设H_rr的零空间为其中r＝Rank(H_rr)，并定义W＝PU，与且满足Q＝PVP^H，其中Q_U＝vec(U)vec(U)^H，U，V为新的辅助变量，然后通过辅助变量法进行等价转换，可将非凸问题(6)转换为优化问题(7)，包括：

s.t.ln(1+γ_d(Q_U，V))-ln(1+γ_k(Q_U，V))≥c²，k∈K (7b)

β(Tr(DQ_U)+Tr(V))+P_c+P_SIC≤b/ln 2 (7c)

c²/b≥a (7d)

ξ(Tr((A_k+B_k)Q_U)+Tr(C_kV))≥η_k，k∈K (7f)

Rank(Q_U)＝1 (7g)

其中，a，b，c为辅助变量，B_i，C_i，D为便于优化过程处理的信道变换矩阵，且 其中G_ri＝P^HH_riP， i∈{d，k}，I表示单位矩阵，表示克罗内克积，G_ri，H_sr，H_ri为信道矩阵，为H_rr的零空间，其中r为H_rr的秩，即r＝Rank(H_rr)；

优化问题(7)仍为非凸问题，对非凸约束条件(7g)，可等价表示为：

将公式(8)作为惩罚函数带入到目标函数(7a)并采用一阶泰勒展开对其进行处理，优化问题(7)可转换为以下问题，即：

s.t.(7b)-(7f) (9b)

其中，和分别为矩阵的最大特征值与对应的特征矢量。

优化问题(9)仍为非凸问题，对非凸约束条件(7b)定义以下函数：

f_k(Q_U，V)＝ln(Tr((A_d+B_d)Q_U)+Tr(C_dV)+1)+ln(Tr((B_k+C_k)V)+1) (10a)

g_k(Q_U，V)＝ln(Tr(A_k+B_k)Q_U)+Tr(C_kV)+1)+ln(Tr((B_d+C_d)V)+1) (10b)

结合(10a)、(10b)并采用连续凸近似的方法对约束条件(7b)、(7d)进行处理，优化问题(9)可以转换成以下凸近似问题：

(7c)，(7e)-(7f) (11d)

其中，f_k(Q_U，V)，为优化过程中间函数变量，为优化问题(6)可行解，且：

最后，设计一种快速收敛的迭代优化算法，可获得系统保密能效最大化下的最佳中继节点波束成形矩阵和人工噪声协方差矩阵。所提优化算法如下：

步骤1：设置源节点的发送功率P_s，中继最大发送功率中继收发天线数M_r、M_t及各节点噪声功率等参数；

步骤2：设置收敛精度ε及迭代序号i＝0；

设置一个初始可行解

设置i＝i+1，基于迭代优化算法求解近似凸转换后的波束成形优化问题模型，得到全双工无线携能中继通信系统保密能效最大化下的最优解其中a,b,c为辅助变量；

当aⁱ-a^i-1>ε时，设置并重新设置i＝i+1，基于所述迭代优化算法求解近似凸转换后的波束成形优化问题模型，否则算法收敛；

基于收敛后的算法，计算最优解

基于最优解计算最佳中继节点波束成形矩阵W^*和人工噪声协方差矩阵Q^*：Q^*＝PV^*P^H，其中vec^-1(A)表示逆向量化运算。

图4是所提算法迭代性能与最优性能仿真图。从图4可以看出最优化过程可以在10以内的迭代次数得到最优性能，即达到收敛。验证所提出算法可在较小迭代次数下快速收敛。

图5是本发明所提出的保密能量效率最大化方案与半双工保密能量效率最大化(SEEM-HD)方案、全双工保密率最大化(SRM-FD)方案和全双工能效最大化(EEM-FD)三种方案的保密能效对比图。从图5可以看出，SEEM-FD方案的性能优于SEEM-HD，SRM-FD和EEM-FD方案。相较于SEEM-HD，SEEM-FD可获得更高的保密率，提升了保密能效性；对于SRM-FD，在系统保密能效性达到最佳值仍以全功率发送信号，导致系统能效性在后开始下降。而对于EEM-FD方案，因其未考虑系统安全性，无法有效应对能量采集节点的窃听，保密能效相较于SEEM-FD性能显著下降。

实施例二：

本发明实施例二提供了一种全双工无线携能中继通信系统高能效保密传输系统，包括：

模型构建模块，用于构建基于人工噪声辅助的波束成形优化问题模型，该模型以全双工无线携能中继节点处的保密能量效率最大化为目标；

凸转换模块，用于采用连续凸近似和惩罚函数，对波束成形优化问题模型进行近似凸转换；

模型计算模块，用于采用迭代优化算法，计算近似凸转换后的波束成形优化问题模型，获得全双工无线携能中继通信系统保密能量效率最大化下的最佳中继节点波束成形矩阵和人工噪声协方差矩阵。

实施例三：

本发明实施例三提供了一种全双工无线携能中继通信系统高能效保密传输装置，包括处理器及存储介质；

所述存储介质用于存储指令；

所述处理器用于根据所述指令进行操作以执行根据实施例一中所述方法的步骤。

实施例四：

本发明实施例四提供了计算机可读存储介质，其上存储有计算机程序，该程序被处理器执行时实现实施例一中所述方法的步骤。

本领域内的技术人员应明白，本申请的实施例可提供为方法、系统、或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。

本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。

这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。

这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。

以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明技术原理的前提下，还可以做出若干改进和变形，这些改进和变形也应视为本发明的保护范围。