CN111542121B - 一种应用于双向df中继系统的满足swipt的多维资源分配方法 - Google Patents

Abstract

本发明公开了一种应用于双向DF中继系统的满足SWIPT的多维资源分配方法。该方法旨在联合优化源端功率分配、中继信能功率分割比以及传输时间分配来提升系统的中断性能，其基本思想是首先依照系统模型建立最小化系统中断概率的优化模型，然后将优化中断概率的原优化问题等价转换为最大归一化信噪比问题，最后对优化问题进行降维处理并结合两步优化算法和放缩法求出各资源分配的闭式解。本发明降低了双向DF携能中继系统多维资源分配的复杂性，为中断性能的提升提供了一种有效的资源分配方案。

Description

一种应用于双向DF中继系统的满足SWIPT的多维资源分配 方法

技术领域

本发明涉及通信系统技术领域，具体涉及一种应用于双向DF中继系统的满足SWIPT的多维资源分配方法。

背景技术

随着物联网(Internet of Things,IoT)的发展，越来越多的小功率无线节点加入网络，大量的无线接入以及快速增长的通信需求导致能量迅速消耗，这严重制约着以能量受限传感器节点为基础的物联网的传输性能和使用寿命。而传统的电网充电和电池更换在一些特定传感器网络场景难以实现，如大规模传感器网络、嵌入式器件和偏远地区网络等，急需一种有效的供能机制为这类能量受限的传输网络提供持续的能量补给。这些年的研究表明利用无线射频能量收集的信息能量同传(Simultaneous Wireless Information andPower Transfer,SWIPT)技术是解决能量受限问题的重要技术，备受大量研究者的关注。

协作/中继网络是SWIPT技术的重要研究场景，通过利用RF能量收集技术，中继节点无需消耗自身能量便能提供可控、有效的能量传输和信息传输，这有利于吸引空闲节点担任系统中继，同时提升能量受限中继网络的传输持续性。目前，SWIPT技术在双向中继传输网络中的应用研究已有较坚实的基础，但是在联合资源优化部分仍然存在一些不完善的地方，大部分的研究集中在单一参数的优化改善上，例如通过折衷优化信息和能量传输提升系统的传输速率或者可靠性。通过优化功率分配提升系统的传输性能等。考虑到双向DF(译码转发协作方式，Decode and Forward)携能中继系统存在能量资源、时间资源、信息资源等多维度的资源，如何充分利用这些资源进行系统性能的改善，许多关键技术还需展开深入研究，因此重新设计合适的资源分配方案来提升系统的中断性能十分必要。

发明内容

针对现有技术存在的上述不足，本发明实际需要解决的问题是：如何通过优化传输设计进一步改善系统的中断性能。

为解决上述技术问题，本发明采用了如下的技术方案：

一种应用于双向DF中继系统的满足SWIPT的多维资源分配方法，所述双向DF中继系统包括两个源节点S₁和S₂，还包括一个配置PS能量收集接收机的中继R，其中，S₁和S₂的发送功率分别为P₁和P₂，P₁+P₂＝P_t，P_t为传输总功率，中继R为能量受限节点，不为传输提供自身能量，需要依靠无线能量收集获取能量再进行接下来的信息传输；所述应用于双向DF中继系统的满足SWIPT的多维资源分配方法包括：

S101、推导所述双向DF中继系统的中断概率公式，以最小化所述双向DF中继系统的中断概率为目标，构建传输总功率限制下联合优化功率分配、信能功率分割比分配和时间分配的优化问题OP0，

s.t.0≤α≤1，0≤ρ≤1，0≤t≤1

式中，α表示用户功率分配参数，ρ表示信能功率分割比值，t表示时间分配参数，α^o表示最优用户功率分配闭式解，ρ^o表示最优信能功率分割比值，t^o表示最优时间分配参数闭式解，P_out(α,ρ,t)表示由用户功率分配参数α、信能功率分割比值ρ和时间分配参数t决定的终端概率；

S102、基于信道状态信息，将OP0转换为最大归一化信噪比的联合优化问题OP1，

s.t.0≤α≤1，0≤ρ≤1，0≤t≤1

式中，Λ(α,ρ,t)为α、ρ及t对应的归一化信噪比；

S103、求出固定用户功率分配参数α和时间分配参数t的最优信能功率分割比值的表达式，然后将ρ^o代入Λ(α,ρ,t)，将OP1进行降维处理得到一个二维优化问题；

S104、对所述二维优化问题，固定时间分配参数t，采用个例分析方法求出最优用户功率分配闭式解α^o；

S105、对所述二维优化问题，固定用户功率分配参数α，采用放缩法和近似法求解得到满足最大归一化信噪比的紧上界值的最优时间分配参数闭式解t^o；

S106、将求出的最优用户功率分配闭式解α^o及最优时间分配参数闭式解t^o代入最优信能功率分割比值的表达式得到最优信能功率分割比值ρ^o；

S107、基于最优用户功率分配闭式解α^o、最优时间分配参数闭式解t^o及最优信能功率分割比值ρ^o进行资源分配。

优选地，所述双向DF中继系统的传输过程包括MAC阶段和BC阶段，MAC阶段和BC阶段的持续时间分别为t和1-t，MAC阶段结束后中继R收集到的能量和需要进行解码传输的信息分别为E和Y_ID，

式中，η表示能量收集转化能量的效率，h₁表示节点S₁到R间的信道参数，h₂表示表示节点S₂到R间的信道参数，P₁表示节点S₁的发送功率，P₂表示节点S₂的发送功率，

表示信号经历宽带转换窄带过程中遭受的噪声，

表示信号经历宽带转换窄带过程中的噪声功率，

表示接收天线处生成的噪声，

表示天线处生成噪声功率，x₁表示节点S₁的传输信息，x₂表示节点S₂的传输信息；

在BC阶段，中继R利用网络编码理论解码y_ID并构造网络编码码字

然后，采用收集到的能量P_R＝E/(1-t)将x_R广播给节点S_i,i＝1,2，S_i收到的信息为

n_i表示第S_i处获得的噪声。

优选地，步骤S1中，所述双向DF中继系统的中断概率公式的推导过程如下：

定义P₁＝αP_t，P₂＝(1-α)P_t，求出端到端传输速率，

R_1R＝t·log₂(1+γ_1R)

R_R2＝(1-t)·log₂(1+γ_R2)

R_2R＝t·log₂(1+γ_2R)

R_R1＝(1-t)·log₂(1+γ_R1)

R_MA＝t·log₂(1+γ_MA)

式中，R_1R表示节点S₁到中继R的传输速率，γ_1R＝(1-ρ)αγ₁表示节点S₁到中继R的传输信噪比，R_2R表示节点S₂到中继R的传输速率，γ_2R＝(1-ρ)(1-α)γ₂表示节点S₂到中继R的传输信噪比，R_Ri,i＝1,2表示中继R到节点S_i,i＝1,2的传输速率，

表示中继R到节点S_i,i＝1,的传输信噪比，R_MA表示多址传输的可达传输速率，γ_MA＝(1-ρ)(αγ₁+(1-α)γ₂)表示多址信噪比，

当节点的速率阀值对(R_th1,R_th2)超出系统的可达传输速率域D则系统中断，则中断概率公式为：

式中，Pr{}表示概率求解公式，

优选地，步骤S103中：

当α和t固定时，最优问题OP1降解为一个仅由信能功率分割比值ρ决定的一维优化问题，该优化问题为：

Λ(ρ)＝min((1-ρ)g₁(α,t),ρg₂(α,t))，其中，

(1-ρ)g₁(α,t)和ρg₂(α,t)分别为关于ρ的单调递减函数和单调递增函数，最优信能功率分割比值在(1-ρ)g₁(α,t)＝ρg₂(α,t)时获得，

优选地，步骤S104中：

式中，

优选地，步骤S105中，以

及

代替g₁(α,t)及g₂(α,t)，

将

及

的改写为

式中，A_j和R_thj分别对应为：((A_j,R_thj))∈{(αγ₁,R_th1),((1-α)γ₂,R_th2),(αγ_-+γ₂,R_thΣ)}，j＝1,2，得到

式中，

表示将

代入Λ(α,ρ,t)后的归一化信噪比，且

存在一个紧上界值

将紧上界值作为优化目标进行求解得到最优时间分配参数闭式解t^o满足：

当|h₁|²＞|h₂|²，

当|h₁|²≤|h₂|²，

式中，

综上所述，本发明公开的一种应用于双向DF中继系统的满足SWIPT的多维资源分配方法与现有技术相比，在进行传输模型建模时，充分考虑了影响该双向DF中继系统传输性能的多维度参数以及它们之间的相关性，构建了双向DF携能中继网络的多维资源优化问题；在性能优化上，以最小化系统中断概率为优化目标，联合优化了功率分配时间分配和信能功率分割比分配参数；在优化设计中，利用降维处理、放缩法和近似法以及两步优化算法，求出了参数值的闭式解，并证明了该求解方法得到的性能曲线是理论曲线的近似解，通过优化传输设计进一步改善了系统的中断性能。

附图说明

图1为本发明公开的一种应用于双向DF中继系统的满足SWIPT的多维资源分配方法的一个具体实施例的流程图；

图2为本发明中的双向DF中继系统的传输模型示意图；

图3为本发明中的双向DF中继系统的信号模型示意图；

图4为本发明中所提方案与数值仿真的比较图；

图5为本发明中所提优化方法以及其它对比优化方案与传输总功率的关系图；

图6为本发明中所提优化方法以及其它对比优化方案与传输距离的关系图。

具体实施方式

为了使本发明的目的、技术方案和优点更加清楚，下面将结合附图对本发明作进一步的详细描述说明。

术语解释：

SWIPT(Simultaneous Wireless Information and Power Transfer)：无线信息能量同传；

PS(Power Splitting)：功率分割；

DF(Decode and Forward)：解码转发；

MAC(Multiple Access)：多址；

BC(Broadcast)：广播

如图1所示，本发明公开了一种应用于双向DF中继系统的满足SWIPT的多维资源分配方法的一个具体实施例，如图2所示，为所述双向DF中继系统的传输和信号模型，所述双向DF中继系统包括两个源节点S₁和S₂，还包括一个配置PS能量收集接收机的中继R，其中，S₁和S₂的发送功率分别为P₁和P₂，P₁+P₂＝P_t，P_t为传输总功率，中继R为能量受限节点，不为传输提供自身能量，需要依靠无线能量收集获取能量再进行接下来的信息传输；所述应用于双向DF中继系统的满足SWIPT的多维资源分配方法包括：

S101、推导所述双向DF中继系统的中断概率公式，以最小化所述双向DF中继系统的中断概率为目标，构建传输总功率限制下联合优化功率分配、信能功率分割比分配和时间分配的优化问题OP0，

s.t.0≤α≤1，0≤ρ≤1，0≤t≤1

式中，α表示用户功率分配参数，ρ表示信能功率分割比值，t表示时间分配参数，α^o表示最优用户功率分配闭式解，ρ^o表示最优信能功率分割比值，t^o表示最优时间分配参数闭式解，P_out(α,ρ,t)表示由用户功率分配参数α、信能功率分割比值ρ和时间分配参数t决定的中断概率；

如图3所示，在步骤S101中，考虑整个传输阶段MAC和BC两个传输阶段，令MAC阶段的持续时间为t，BC阶段的持续时间为1-t。在MAC传输阶段中，中继节点R接收到的信息为

其中,P_i,i＝1,2和x_i分别为源节点S_i的发送功率和传输信息，且E{|x_i|²}＝1，h_i为节点S_i到R间的信道参数，

为接收天线处生成的噪声，

表示天线处生成噪声功率。

中继配备功率分割接收机，将接收到的信息y_R分为ρ:(1-ρ)两部分，

用于能量收集，

用于信息处理。能量收集转化能量的效率为η，则中继处收集到的能量为：

信息处理部分的信号可以表示为

其中，

n_r,b～CN(0,σ_b ²)为信号经历宽带转换窄带过程中遭受的噪声，

表示信号经历宽带转换窄带过程中的噪声功率。

在BC阶段，中继利用网络编码理论解码y_ID并构造出网络编码码字

然后采用收集到的能量P_R将x_R广播给源节点S_i,i＝1,2。源节点S_i在BC阶段收到的信息为

其中，

为S_i处获得的噪声。在接下来的计算中，不失一般性地忽略噪声

部分并假设各个节点收到的噪声影响相同，即令

源节点S_i,i＝1,2接收到信息y_Si后，首先解码信息包x_R，然后利用自身包含的信息x_i采用自干扰消除(Self-Cancellation)方法解码出需要的对方用户信息。例如，S₁解码S₂的信息：

定义P₁＝αP_t，P₂＝(1-α)P_t，其中0≤α≤1，定义R_ij,i,j∈[1,2,R]为节点i到节点j的端到端(End-to-End)可达传输速率，R_MA为多址传输的可达传输速率。根据上述模型可求出R_ij,i,j∈[1,2,R]公式如下：

根据中断概率的定义：当源节点的速率阀值对(R_th1,R_th2)超出系统的可达传输速率域D则系统中断。写出此模型下的中断概率公式为：

其中，P rx{为x的概率求解公式,

根据上述公式可以看出当传输速率阀值和信道条件固定时，中断概率P_out由用户功率分配参数α、信能功率分割比值ρ和时间分配参数t决定，以最小化系统中断概率为目标，构建优化问题如下：

S102、基于信道状态信息，将OP0转换为最大归一化信噪比的联合优化问题OP1，

s.t.0≤α≤1，0≤ρ≤1，0≤t≤1

式中，Λ(α,ρ,t)为α、ρ及t对应的归一化信噪比；

在步骤S102中，当每个节点均知道所有的瞬时信道状态信息(InstantaneousChannel State Information)时，则公式(4)的优化问题可以等价于最大归一化信噪比Λ的问题，等价优化问题如下：

其中，

是包含多个变量(α、ρ和t)的复杂函数。

S103、求出固定用户功率分配参数α和时间分配参数t的最优信能功率分割比值的表达式，然后将ρ^o代入Λ(α,ρ,t)，将OP1进行降维处理得到一个二维优化问题；

当固定参数α和t时，最优问题OP1降解为一个仅由信能功率分割比值ρ决定的一维优化问题，该优化问题为：

由于Λ(ρ)可写为Λ(ρ)＝min((1-ρ)g₁(α,t),ρg₂(α,t))，其中，g₁(α,t)，g₂(α,t)为

由公式(7)可以看出，当给定变量α和t时，(1-ρ)g₁(α,t)和ρg₂(α,t)分别为关于ρ的单调递减函数和单调递增函数。因此，最优信能功率分割比值在(1-ρ)g₁(α,t)＝ρg₂(α,t)时获得，即：

S104、对所述二维优化问题，固定时间分配参数t，采用个例分析方法求出最优用户功率分配闭式解α^o；

由公式(8)可以看出，所求出的ρ^o是变量α和t的函数。将所求出的ρ^o代入Λ可以将原优化目标函数改写为一个关于α和t的函数：

进一步地采用分步优化求解所得到的的优化目标函数Λ(α,t)中的最优α和t。

当固定时间分配参数t时，Λ(α,t)降解为一个仅由功率分配参数α决定的一维函数，令此时的目标函数为

并将优化问题等价于最小化目标函数Λ(α)的倒数

则得到：

分析可知，通过判断α₁₃和α₂₃之间的大小关系，可以将

分为两种情形，分别为：

情形1：当α₁₃≥α₂₃

情形2：当α₁₃＜α₂₃

其中，

α₁₂，α₁₃和α₂₃分别为q₁与q₂，q₁与q₃和q₂与q₃之间的交点，

由公式(11)至(12)可知

是一个满足凸函数性质的分段函数，因此它的最优值受边界值α＝0、α＝1，分段交点值α₁₂、α₁₃、α₂₃，以及令

得到的值影响，因此，当给定任意时间分配值t，得到满足该系统在传输总功率限制下的最优功率分配值为：

其中，

并且

以下是对上述最优用户功率分配闭式解α^o的证明：

结合上述得到的

的分段数学模型，可以计算得到

的导数

如下所示：

情形1：α₁₃≥α₂₃

情形2：α₁₃＜α₂₃

由公式(17)和(18)可知，

的单调性受γ₁和γ₂的大小关系影响。因此，

的单调性存在如下三种子情况：

子情况1:γ₁＞γ₂

当γ₁＞γ₂，情形1中的

在区间α∈(0，α₁₂]内

单调递减并在α＝α₁₂处取得最小值；在区间α∈[α₁₂,1)内，

的值受α的变化影响，使得在区间

单调递减，在区间

单调递增，因此，在区间，α∈[α₁₂,1)内，

在

处取得最小值。综上所述，在整个区间α∈(0，1)，

在

处取得最小值：

其中

令

的

第二段函数等于零求出α^*如公式(15)所示。

当γ₁＞γ₂，情形2中的

在区间α∈(0，α₁₃]和区间α∈[α₁₃,α₂₃]内均满足

在这两个区间单调递减，由于

是连续函数，因此，在区间α∈(0，α₂₃]内，

在α＝α₂₃处取得最小值；在区间α∈[α₂₃,1)内，

的值受α的变化影响，随着α的增加先减小再增加(或单调递增)，并在

处取得最小值，

其中，

子情况2:γ₁＜γ₂

当γ₁＜γ₂，情形1中的

在区间α∈[α₁₂,1)内

单调递增并在α＝α₁₂处取得最小值；在区间α∈(0，α₁₂]内，

的值受α的变化影响，使得在区间

单调递减，在区间

单调递增，因此，在区间，α∈(0，α₁₂]内，

在

处取得最小值。综上所述，在整个区间

在

处取得最小值：

其中

令

的第一段函数为零求出α⁺如公式(16)所示。

当γ₁＜γ₂，情形2中的

在区间α∈[α₁₃,α₂₃]和区间α∈[α₂₃,1)内均满足

在这两个区间单调递增，由于

是连续函数，因此，在区间α∈[α₁₃,1)内，

在α＝α₁₃处取得最小值；在区间α∈(0，α₁₃]内，

的值受α的变化影响，随着α的增加先减小再增加(或单调递减)，并在

处取得最小值，

其中，

子情况3:γ₁＝γ₂

当γ₁＝γ₂，g₂(α)为一个常数，

的单调性仅由g₁(α)决定。情形1中的

在区间α∈[α₁₂,1)内

单调递增，在区间α∈(0，α₁₂]内

单调递减，因此，最小值在α＝α₁₂处获得，α^o＝α₁₂。情形2中，在区间α∈(0，α₁₃]内

单调递减，在区间α∈[α₁₃,α₂₃]内，

在区间α∈[α₂₃,1)内

单调递减，因此，最小值在区间α∈[α₁₃,α₂₃]内获得，

由于α₁₂∈[α₁₃,α₂₃]，可令α^o＝α₁₂。

S105、对所述二维优化问题，固定用户功率分配参数α，采用放缩法和近似法求解得到满足最大归一化信噪比的紧上界值的最优时间分配参数闭式解t^o；

由于g₁(α,t)、g₂(α,t)中包含时间的部分是一个指数函数，考虑到当x远大于1时，满足2^x-1～2^x，因此，为方便求出闭式解，将g₁(α,t)、g₂(α,t)用以下近似值替代：

对于任一

可以写为

任一

可以写为

其中，A_j、R_thj依次对应于((A_j,R_thj))∈{(αγ₁,R_th1),((1-α)γ₂,R_th2),(αγ_-+γ₂,R_thΣ)}，j＝1,2。结合上述得到的

可以得到

公式(20)是一个凸函数，可以看出只要分母取得最小值，则

便能取得最大值。考虑到分母中的第二项中包含信道参数|h_3-i|²，当存在大尺度衰落时，信道参数的值十分的小，此时上式存在一个紧上界值：

以紧上界值作为

的值，即令

分析可得，当给定任意功率分配值α，可求得最优时间分配参数的近似优化闭式解为：

当|h₁|²＞|h₂|²时，则

当|h₁|²≤|h₂|²时，则

其中，

以下是对上述最优时间分配参数闭式解to的证明：由紧上界值可知，当固定α时，

是一个关于t的一维凸函数，由于变量t仅存在于分母中，最大化

等价于最小化

已知Q＝e^lnQ，e^lnQ的单调性与lnQ的单调性一致，因此，可通过求解lnQ的最小值，获得Q的最小值。将lnQ展开，可得到

对其求导，得到

分析上式可知，lnQ是一个凹函数并在

时取得最小值，进而得到

在

取得最大值。

另一方面，由于

中的Q是由

提供，并且Q与

呈反比。由公式(19)可知，

是一个分段函数，因此Q也是一个分段函数，并且满足

当|h₁|²＞|h₂|²时，则

当|h₁|²≤|h₂|²时，则

其中，

分别讨论上述两种情形下Q的最小值可知，当|h₁|²＞|h₂|²时，Q的第一段分段函数是一个凹函数并在

处取得最小值，Q的第二段分段函数是一个凹函数并在

处取得最小值。比较t₀、t₁、t₂可知，若t₁＜t₀，则t₂＜t₀；若t₂＞t₀，则t₁＞t₀。因此，可以得出结论，当t₁＜t₀，Q的最小值在t₁处获得；当t₂＞t₀，Q的最小值在t₂处获得；当t₂＜t₀＜t₁，Q的最小值在t₀处获得。将上述推导总结得到最优时间分配值t^o如公式(22)所示。

同理，采用相同方法推导可知，当|h₁|²≤|h₂|²时，令Q取得最小值的最优时间分配值t^o如公式(23)所示。

S106、将求出的最优用户功率分配闭式解α^o及最优时间分配参数闭式解t^o代入最优信能功率分割比值的表达式得到最优信能功率分割比值ρ^o；

S107、基于最优用户功率分配闭式解α^o、最优时间分配参数闭式解t^o及最优信能功率分割比值ρ^o进行资源分配。

如图4所示，为本发明公开的应用于双向DF中继系统的满足SWIPT的多维资源分配方法与数值仿真方法得到的中断概率结果对比图。由图可以看出，本方法得到的性能曲线是数值仿真得到的理论曲线的近似解。如图5及图6所示，为本发明公开的应用于双向DF中继系统的满足SWIPT的多维资源分配方法与其他方法的中断概率结果对比图。其中，图4给出了本发明中所提优化方法以及其它对比优化方案与传输总功率P_t的关系图，仿真过程中选取用户传输速率R_th1＝0.5bit/s/Hz，R_th2＝2bit/s/Hz,图中d₁和d₂分别表示S₁到R的距离，以及S₂到R的距离。图5给出了本发明中所提优化方法以及其它对比优化方案与传输距离d₁的关系图，仿真过程中选取传输总功率为P_t＝20dBm，图中R_th1和R_th2分别表示源端用户S₁的传输速率，以及源端用户S₂的传输速率。由图4和图5可知，随着传输总功率的增加，系统的中断概率降低，并且，中继R的位置影响系统的中断概率。此外，在整个仿真过程中可以看出，本发明的方法(图4及图5中的多维联合优化)与其他方法相比，中断概率更低。

综上所述，本发明公开的一种应用于双向DF中继系统的满足SWIPT的多维资源分配方法与现有技术相比，在进行传输模型建模时，充分考虑了影响该双向DF中继系统传输性能的多维度参数以及它们之间的相关性，构建了双向DF携能中继网络的多维资源优化问题；在性能优化上，以最小化系统中断概率为优化目标，联合优化了功率分配时间分配和信能功率分割比分配参数；在优化设计中，利用降维处理、放缩法和近似法以及两步优化算法，求出了参数值的闭式解，并证明了该求解方法得到的性能曲线是理论曲线的近似解，通过优化传输设计进一步改善了系统的中断性能。

最后说明的是，以上实施例仅用以说明本发明的技术方案而非限制，尽管通过参照本发明的优选实施例已经对本发明进行了描述，但本领域的普通技术人员应当理解，可以在形式上和细节上对其作出各种各样的改变，而不偏离所附权利要求书所限定的本发明的精神和范围。

Claims (1)

1.一种应用于双向DF中继系统的满足SWIPT的多维资源分配方法，其特征在于，所述双向DF中继系统包括两个源节点S₁和S₂，还包括一个配置PS能量收集接收机的中继R，其中，S₁和S₂的发送功率分别为P₁和P₂，P₁+P₂＝P_t，P_t为传输总功率，中继R为能量受限节点，不为传输提供自身能量，需要依靠无线能量收集获取能量再进行接下来的信息传输；所述应用于双向DF中继系统的满足SWIPT的多维资源分配方法包括：

S101、推导所述双向DF中继系统的中断概率公式，以最小化所述双向DF中继系统的中断概率为目标，构建传输总功率限制下联合优化功率分配、信能功率分割比分配和时间分配的优化问题OP0，

s.t.0≤α≤1，0≤ρ≤1，0≤t≤1

式中，α表示用户功率分配参数，ρ表示信能功率分割比值，t表示时间分配参数，α^o表示最优用户功率分配闭式解，ρ^o表示最优信能功率分割比值，t^o表示最优时间分配参数闭式解，P_out(α,ρ,t)表示由用户功率分配参数α、信能功率分割比值ρ和时间分配参数t决定的中断概率；

S102、基于信道状态信息，将OP0转换为最大归一化信噪比的联合优化问题OP1，

s.t.0≤α≤1，0≤ρ≤1，0≤t≤1

式中，Λ(α,ρ,t)为α、ρ及t对应的归一化信噪比；

S103、求出固定用户功率分配参数α和时间分配参数t的最优信能功率分割比值的表达式ρ^o，然后将ρ^o代入Λ(α,ρ,t)，将OP1进行降维处理得到一个二维优化问题；

S104、对所述二维优化问题，固定时间分配参数t，采用个例分析方法求出最优用户功率分配闭式解α^o；

S105、对所述二维优化问题，固定用户功率分配参数α，采用放缩法和近似法求解得到满足最大归一化信噪比的紧上界值的最优时间分配参数闭式解t^o；

S106、将求出的最优用户功率分配闭式解α^o及最优时间分配参数闭式解t^o代入最优信能功率分割比值的表达式得到最优信能功率分割比值ρ^o；

S107、基于最优用户功率分配闭式解α^o、最优时间分配参数闭式解t^o及最优信能功率分割比值ρ^o进行资源分配；

所述双向DF中继系统的传输过程包括MAC阶段和BC阶段，MAC阶段和BC阶段的持续时间分别为t和1-t，MAC阶段结束后中继R收集到的能量和需要进行解码传输的信息分别为E和Y_ID，

式中，η表示能量收集转化能量的效率，h₁表示节点S₁到R间的信道参数，h₂表示节点S₂到R间的信道参数，P₁表示节点S₁的发送功率，P₂表示节点S₂的发送功率，y_R表示R接收到的信息，

表示信号经历宽带转换窄带过程中遭受的噪声，

表示信号经历宽带转换窄带过程中的噪声功率，

表示接收天线处生成的噪声，

表示天线处生成噪声功率，x₁表示节点S₁的传输信息，x₂表示节点S₂的传输信息；

在BC阶段，中继R利用网络编码理论解码y_ID并构造网络编码码字

然后，采用收集到的能量P_R＝E/(1-t)将x_R广播给节点S_i,i＝1,2，S_i收到的信息为

表示第S_i处获得的噪声；

步骤S1中，所述双向DF中继系统的中断概率公式的推导过程如下：

定义P₁＝αP_t，P₂＝(1-α)P_t，求出端到端传输速率，

R_1R＝t·log₂(1+γ_1R)

R_R2＝(1-t)·log₂(1+γ_R2)

R_2R＝t·log₂(1+γ_2R)

R_R1＝(1-t)·log₂(1+γ_R1)

R_MA＝t·log₂(1+γ_MA)

式中，R_1R表示节点S₁到中继R的传输速率，γ_1R＝(1-ρ)αγ₁表示节点S₁到中继R的传输信噪比，R_2R表示节点S₂到中继R的传输速率，γ_2R＝(1-ρ)(1-α)γ₂表示节点S₂到中继R的传输信噪比，R_Ri,i＝1,2表示中继R到节点S_i,i＝1,2的传输速率，

表示中继R到节点S_i,i＝1,2的传输信噪比，R_MA表示多址传输的可达传输速率，γ_MA＝(1-ρ)(αγ₁+(1-α)γ₂)表示多址信噪比，

当节点的速率阀值对(R_th1,R_th2)超出系统的可达传输速率域D则系统中断，则中断概率公式为：

式中，Pr{}表示概率求解公式，

步骤S103中：

当α和t固定时，最优问题OP1降解为一个仅由信能功率分割比值ρ决定的一维优化问题，该优化问题为：

Λ(ρ)＝min((1-ρ)g₁(α,t),ρg₂(α,t))，其中，

(1-ρ)g₁(α,t)和ρg₂(α,t)分别为关于ρ的单调递减函数和单调递增函数，最优信能功率分割比值在(1-ρ)g₁(α,t)＝ρg₂(α,t)时获得，

步骤S104中：

式中，

步骤S105中，以

及

代替g₁(α,t)及g₂(α,t)，

将

及

的改写为

式中，A_j和R_thj分别对应为：((A_j,R_thj))∈{(αγ₁,R_th1),((1-α)γ₂,R_th2),(αγ_-+γ₂,R_thΣ)}，j＝1,2，得到

式中，

表示将

代入Λ(α,ρ,t)后的归一化信噪比，且

存在一个紧上界值

将紧上界值作为优化目标进行求解得到最优时间分配参数闭式解t^o满足：

当|h₁|²＞|h₂|²，

当|h₁|²≤|h₂|²，

式中，

Images