[发明专利]一种协作认知网络中公平且能效高的联合资源分配方法

权利要求书

1.一种协作认知网络中公平且能效高的联合资源分配算法，由以下协作认知无线电系统来实现：该系统包括主要系统和次要系统两部分，主要系统由主要用户发射端PT、主要用户接收端PR组成，次要系统由N对次要用户发射端ST_n和次要用户接收端SR_n组成，其中

n∈U，表示第n个次要用户，集合U＝{1,2,3,...,N}，次要系统在主要系统运作过程中作为其中继，协助其传输信号，中继模式为译码转发；本系统中共N个次要用户，K个子载波，设子载波集合S＝{1,2,3,...,K}，设γ_k,0、γ_n,k,1、γ_n,k,2和γ_n,k,3分别为主要用户发射端对主要用户接收端、主要用户发射端对第n个次要用户发射端、第n个次要用户发射端对主要用户接收端和第n个次要用户发射端对第n个次要用户接收端链路的信道功率增益，其中k∈S，表示第k个子载波，n∈U，各链路对应的信号发射功率分别为p_k,0，p_k,0，p_n,k,2和p_n,k,3；该方法的具体步骤如下：

1)计算主要系统的平均速率

第一传输阶段，主要用户发射端通过K个子载波发射信号至主要用户接收端，这时所有的次要用户发射端均能监听到信号，每个次要用户发射端接收子载波集合表示为满足其中符号U表示对集合求并集，因此，主要用户发射端至各次要用户发射端的平均速率可表示为：

R 1 = E [ 1 2 Σ n = 1 N Σ k ∈ Ω n 1 l n ( 1 + γ n , k , 1 p k , 0 ) ] - - - ( 1 ) ]]>

其中符号E[]是对括号中的部分求数学期望，符号∑表示对其上下标所限制的范围内进行求和；

第二传输阶段，次要用户发射端对接收的信号进行重新编码并重传，因此之前分配的子载波也被打乱重新分配，重新分配后的子载波集合表示为满足其中

表示第n个次要用户发射端用来给主要用户接收端传送信号使用的子载波集合，而表示剩余的用来给第n个次要用户接收端传送信号的子载波集合，设集合因此，主要用户接收端处的平均速率可以表示为：

R 2 = E [ 1 2 Σ n = 1 N Σ k ∈ Ω n p l n ( 1 + γ n , k , 2 p n , k , 2 + γ k , 0 p k , 0 ) + 1 2 Σ k ∈ Π ln ( 1 + γ k , 0 p k , 0 ) ] - - - ( 2 ) ]]>

因此，在次要系统协作之下的主要系统的平均速率可以表示为：

R_P＝min{R₁,R₂} (3)

其中min{}是对括号中部分取最小值；

2)计算次要系统的平均速率与次要用户的平均功率

次要系统的平均速率可以表示为：

R S = E [ 1 2 Σ n = 1 N Σ k ∈ Ω n S l n ( 1 + γ n , k , 3 p n , k , 3 ) ] - - - ( 4 ) ]]>

次要用户的平均发射功率可以表示为：

P n a v e r = E [ Σ k ∈ Ω n P p n , k , 2 + Σ k ∈ Ω n S p n , k , 3 ] - - - ( 5 ) ]]>

3)引入效用函数

为了能使次要用户之间分配的资源更公平，我们引入了效用函数，采用的效用函数定义如下：

V β ( P ‾ n ) = ( P ‾ n ) 1 + β 1 + β - - - ( 6 ) ]]>

其中，指的是每个次要用户的平均功率，参数β≥0，随着β的增加，次要用户之间的资源分配会变得越来越公平；

4)确定优化问题

以次要用户平均功率的效用函数为目标函数，速率、功率限制条件为约束条件，构造如下优化问题：

min i m i z e : Σ n = 1 N V β ( P ‾ n ) s u b j e c t t o : E [ 1 2 Σ n = 1 N Σ k = 1 K α n , k , 1 ln ( 1 + γ n , k , 1 p k , 0 ) ] ≥ R T 1 E [ 1 2 Σ n = 1 N Σ k = 1 K α n , k , 2 ln ( 1 + γ n , k , 2 p n , k , 2 + γ k , 0 p k , 0 ) + 1 2 Σ n = 1 N Σ k = 1 K α n , k , 3 ln ( 1 + γ k , 0 p k , 0 ) ] ≥ R T 1 E [ 1 2 Σ n = 1 N Σ k = 1 K α n , k , 3 ln ( 1 + γ n , k , 3 p n , k , 3 ) ] ≥ R T 2 P ‾ n ≥ E [ Σ k = 1 K α n , k , 2 p n , k , 2 + Σ k = 1 K α n , k , 3 p n , k , 3 ] , ∀ n - - - ( 7 ) ]]>

其中α_n,k,1,α_n,k,2,α_n,k,3表示的是子载波分配，之前在式(1)、(2)、(4)、(5)中，用表示分配的子载波集合，为了方便，我们定义了符号α_n,k,1,α_n,k,2,α_n,k,3∈[0,1]，当时，α_n,k,1＝1，当时，α_n,k,1＝0；当时，α_n,k,2＝1，当时，α_n,k,2＝0；当时，α_n,k,3＝1，当时，α_n,k,3＝0；符号表示对于所有的n都要满足条件；表示主要系统的目标速率，表示次要系统的目标速率，(7)式中的subject to符号及其后面的式子表示为约束式，subject to表示为约束符号，符号minimize表示求最小值符号，(7)式表示在约束式中对主要系统速率、次要系统速率及次要用户功率进行限制的条件下，求解目标函数即符号minimize后的部分的最小值，该最小化问题在下面的描述中也称为原问题；

5)求解优化问题

经验证，上述优化问题的目标函数是凸的，因此上述优化问题存在唯一的最优解，利用拉格朗日对偶理论，可以建立起原最小化问题即原问题与一个最大化问题即对偶问题之间的关联关系，我们研究的原问题具有强对偶性，因此可以通过求解对偶问题而得到原问题的最优值，原问题的对偶函数为：

D ( Λ ) = min i m i z e : Σ n = 1 N V β ( P ‾ n ) + λ ( R T 1 - E [ 1 2 Σ n = 1 N Σ k = 1 K α n , k , 1 ln ( 1 + γ n , k , 1 p k , 0 ) ] ) + ϵ ( R T 1 - E [ 1 2 Σ n = 1 N Σ k = 1 K α n , k , 2 ln ( 1 + γ n , k , 2 p n , k , 2 + γ k , 0 p k , 0 ) + 1 2 Σ n = 1 N Σ k = 1 K α n , k , 3 ln ( 1 + γ k , 0 p k , 0 ) ] ) + μ ( R T 2 - E [ 1 2 Σ n = 1 N Σ k = 1 K α n , k , 3 ln ( 1 + γ n , k , 3 p n , k , 3 ) ] ) + Σ n = 1 N δ n ( E [ Σ k = 1 K ( α n , k , 2 p n , k , 2 + α n , k , 3 p n , k , 3 ) ] - P ‾ n ) - - - ( 8 ) ]]>

其中Λ:＝{λ,ε,μ,δ}是对偶因子，其中符号:＝表示定义，λ,ε,μ,δ表示(7)中四个限制条件的对偶因子，δ_n是对偶因子δ中的第n个元素，对偶函数对应的对偶问题如下：

m a x i m i z e : D ( Λ ) s u b j e c t t o : Λ ≥ 0 - - - ( 9 ) ]]>

即在对偶因子Λ≥0的约束条件下，通过优化Λ求解目标函数即对偶函数D(Λ)的最大值，已知原问题具有强对偶性，所以通过对偶问题(9)式求得的最优值即为原问题的最优值，求解对偶问题最关键之处在于求解最优的对偶因子Λ^*，Λ^*的求解过程具体如下：

A)设置初始迭代次数t＝0，对偶因子初始值Λ(0)为非负实数；

B)当迭代次数为t时，用Λ(t)表示当前更新的对偶因子，基于当前对偶因子Λ(t)求解对偶函数公式(8)，得到对应的最优次要用户发射功率以及最优子载波分配

C)采用以下4式分别更新4种对偶因子：

λ ( t + 1 ) = [ λ ( t ) + s _ λ ( t ) ( R T 1 - E [ 1 2 Σ n = 1 N Σ k = 1 K α n , k , 1 ln ( 1 + γ n , k , 1 p k , 0 ) ] ) ] + ϵ ( t + 1 ) = [ ϵ ( t ) + s _ ϵ ( t ) ( R T 1 - E [ 1 2 Σ n = 1 N Σ k = 1 K α n , k , 2 ln ( 1 + γ n , k , 2 p n , k , 2 + γ k , 0 p k , 0 ) + 1 2 Σ n = 1 N Σ k = 1 K α n , k , 3 ln ( 1 + γ k , 0 p k , 0 ) ] ) ] + μ ( t + 1 ) = [ μ ( t ) + s _ μ ( t ) ( R T 2 - E [ 1 2 Σ n = 1 N Σ k = 1 K α n , k , 3 ln ( 1 + γ n , k , 3 p n , k , 3 ) ] ) ] + δ n ( t + 1 ) = [ δ n ( t ) + s _ δ n ( t ) ( E [ Σ k = 1 K ( α n , k , 2 p n , k , 2 + α n , k , 3 p n , k , 3 ) ] - P ‾ n ) ] + - - - ( 10 ) ]]>

其中符号[]⁺表示[]中的部分取非负值，s_λ(t)、s_ε(t)、s_μ(t)、s_δ_n(t)表示相应对偶因子对应的迭代步长，t为迭代次数，δ_n(t)为更新后的对偶因子δ(t)中的第n个元素；

D)令Λ^*＝Λ(t+1)，若Λ^*满足预定义的数据精度，则输出最优对偶因子Λ^*，否则，令t＝t+1，跳转至步骤B)，继续迭代，直到满足预定义的数据精度；

6)求得公平且最优的次要用户平均功率

设置效用函数参数β的值，并且将得到的最优对偶因子Λ^*一起带入对偶函数公式(8)中得到公平且最优的次要用户平均功率