[发明专利]基于效用的协作多中继认知系统的联合资源分配方法

摘要

基于效用的协作多中继认知系统的联合资源分配方法，属于无线通信技术领域。该发明方法建立了一种两相协作多中继通信系统模型，在主要用户无法达到其目标速率的情况下，次要用户主动作为AF中继协助主要用户达到其目标速率，作为回报次要用户可以利用未使用的OFDM子载波为自己传输数据，从而使协作通信系统资源利用率达到最大。该方法还考虑了各个次要用户的效用问题，以次要用户的α效用和为目标函数，主要用的传输速率和次要用户的总功率为限制条件，通过拉格朗日对偶原理求得最优的对偶因子，进而求得最优的子载波分配策略和功率控制因子，实现最优资源分配。

主权项

一种基于效用的协作多中继认知系统的联合资源分配方法，由以下系统来实现，该系统包括主要系统和次要系统，主要系统由主要用户发射端PT、主要用户接收端PR组成；次要系统由N个次要用户组成，每个次要用户包括一个次要用户发射端STn和一个次要用户接收端SRn，其中n∈U，表示第n个次要用户，集合U＝{1,2,3,...,N}，次要系统通过放大转发中继模式协助主要系统传输信号；本系统中共K个子载波，设子载波集合Ω＝{1,2,3,...,K}，设γk,1，γn,k,2，γn,k',3和γn,k',4分别为主要用户发射端对主要用户接收端、主要用户发射端对第n个次要用户发射端、第n个次要用户发射端对主要用户接收端和第n个次要用户发射端对第n个次要用户接收端链路的信道功率增益，其中k∈Ω，表示第k个子载波，n∈U；主要用户发送端在子载波k上发送信号的功率用次要用户发送端STn利用子载波k'给主要用户接收端和次要用户接收端发送信号的功率分别表示为和该分配方法的具体步骤如下：1)计算主要系统的平均速率第一传输阶段，主要用户发射端通过K个子载波发射信号至主要用户接收端，这时所有的次要用户发射端均能监听到信号，此时主要用户直接链路所能达到的瞬时传输速率为Rpd=ln(1+pkpγk,1)---(1)]]>第二传输阶段，次要用户发射端对接收的信号进行放大并转发，转发过程采用子载波配对策略，即接收子载波k上的信号并用子载波k'进行转发，子载波k、k'不一定相同；定义变量ρk,k′∈{0,1}为子载波配对变量，当ρk,k′＝1时表示第一传输阶段的子载波k和第二传输阶段中的子载波k'进行配对，并且ρk,k′必须满足约束符号表示任意的；由次要用户发送端STn转发主要用户获得的瞬时中继速率为：Rn,k,k′=ln(1+pkpγk,1+pkpγn,k,2pn,k′spγn,k′,31+pkpγn,k,2+pn,k′spγn,k′,3)---(2)]]>因此，在次要系统协作之下的主要系统的平均速率表示为：RP=E[12Σn=1NΣk∈GPΣk′=1Ktn,k,k′ρk,k′Rn,k,k′+12Σk∈G‾PRpd]---(3)]]>其中变量tn,k,k′∈{0,1}定义为子载波对分配变量，tn,k,k′＝1表示子载波对(k,k′)被分配给了次要用户发送端STn,并且tn,k,k′满足约束GP表示第一传输阶段被第二传输阶段中继的子载波集合，表示集合GP在集合Ω中的补集，符号E[·]表示对括号内表达式求数学期望；2)计算次要系统的平均速率与次要用户的平均功率每个次要用户的平均速率表示为：Rns=E[12Σk′∈G‾sηn,k′ln(1+pn,k′ssγn,k′,4)]---(4)]]>其中集合GS表示第二传输阶段用于中继主要用户信号的子载波集合，其中表示集合GS在集合Ω中的补集，变量ηn,k′∈{0,1}定义为子载波分配变量，表示第二传输阶段中未配对子载波的分配情况，ηn,k′＝1表示子载波k'被分配给了次要用户STn，并且ηn,k′满足约束次要系统的平均速率表示为：RS=12E[Σn=1NΣk′∈G‾sηn,k′ln(1+pn,k′ssγn,k′,4)]---(5)]]>次要用户的平均发射功率表示为：Rus=E[Σn=1N(Σk∈GPΣk′=1Ktn,k,k′ρk,k′E[pn,k′sp]+Σk′∈G‾Sηn,k′E[pn,k′ss])]---(6);]]>3)引入效用函数为了能使次要用户之间分配的资源更公平，我们引入了效用函数，采用的效用函数定义如下：Uα(rns)=(rns)1-α/(1-α)ifα≠1ln(rns)ifα=1---(7)]]>其中，指的是每个次要用户的可达到的速率，参数α≥0，随着α的增加，次要用户之间的资源分配会变得越来越公平；4)确定优化问题以次要用户平均功率的效用函数为目标函数，速率、功率限制条件为约束条件，构造如下优化问题：maxt,ρ,η,p,rΣn=1NUα(rns)---(8)]]>Subject to：RP>RT，Pus≤PS其中subject to是条件符号，表示受后面条件限制，RT表示主要用户需要达到的平均速率、PS表示次要系统的总功率；矢量t,ρ,η,p,r分别为符号表示定义，符号{}表示其中表达式的所有可能情况的集合；5)求解优化问题经验证，上述优化问题的目标函数是凸的，因此上述优化问题存在唯一的最优解，利用拉格朗日对偶理论，可以建立起原最小化问题即原问题与一个最大化问题即对偶问题之间的关联关系，我们研究的原问题具有强对偶性，因此可以通过求解对偶问题而得到原问题的最优值，原问题的对偶函数为：g(Ψ)=maxL(tn,k,k′,ρk,k′,ηn,k′,pn,k′sp,pn,k′ss,rns)---(10)]]>其中，max表示求最大值，是拉格朗日函数，其表达式如下：L(tn,k,k′,ρk,k′,ηn,k′,pn,k′sp,pn,k′ss,rns)=Σn=1N(U(rnS)-vnrnS)+((E[μ2Σn=1NΣk∈GPΣk′=1Ktn,k,k′ρk,k′ln(1+pkpγk,1+pkpγn,k,2pn,k′spγn,k′,3pkpγn,k,2+pn,k′spγn,k′,3)+μ2Σk∈G‾Pln(1+pkpγk,1)-λΣn=1NΣk∈GPΣk′=1Ktn,k,k′ρk,k′pn,k′sp])+(E[Σn=1NΣk′∈G‾svn2ηn,k′ln(1+pn,k′ssγn,k′,4)-λΣn=1NΣk′∈G‾Sηn,k′pn,k′ss]))+λPS-μRT---(11)]]>其中是对偶因子，符号表示定义，λ,ν＝{νn},μ为对偶因子，符号{}表示其中表达式的所有可能情况的集合，νn是对偶因子ν的集合中的某一可能情况；对偶函数对应的对偶问题如下：minΨ≥0g(Ψ)---(12)]]>其中min表示求最小值，即在对偶因子Ψ≥0的约束条件下，通过优化Ψ求解目标函数即对偶函数g(Ψ)的最小值，已知原问题具有强对偶性，所以通过对偶问题(12)式求得的最优值即为原问题的最优值，求解对偶问题最关键之处在于求解最优的对偶因子Ψ*，其求解过程具体如下：A)设置初始迭代次数t＝0，对偶因子初始值Ψ(0)为非负实数；B)当迭代次数为t时，用表示当前更新的对偶因子，其中v(t)＝{νn(t)}，基于当前对偶因子Ψ(t)求解对偶函数公式(11)，得到对应的最优次要用户发射功率第n个次要用户可达到的最大速率以及最优子载波对分配变量子载波分配变量ηn,k′，子载波配对变量ρk,k'；C)采用以下三式分别更新当前对偶因子λ(t)，μ(t)，νn(t)：λ(t+1)=[λ(t)-s_λ(t)(PS-Σn=1N(Σk∈GPΣk′=1Ktn,k,k′ρk,k′E[pn,k′sp]+Σk′∈G‾Sηn,k′E[pn,k′ss]))]+]]>μ(t+1)=[μ(t)-s_μ(t)(E[12Σn=1NΣk∈GPΣk′=1Ktn,k,k′ρk,k′Rn,k,k′+12Σk∈G‾PRpd]-RT)]+]]>vn(t+1)=[vn(t)+s_vn(t)(E[12Σk′∈G‾sηn,k′ln(1+pn,k′ssγn,k′,4)]-rns)]+]]>其中符号[]+表示[]中的部分取非负值，s_λ(t)、s_μ(t)、s_νn(t)表示相应对偶因子对应的迭代步长，t为迭代次数，λ(t+1)、μ(t+1)、νn(t+1)分别表示对三种对偶因子λ(t)、μ(t)、νn(t)的更新；D)令Ψ*＝Ψ(t+1)，若Ψ*满足预定义的数据精度，则输出最优对偶因子Ψ*，否则，令t＝t+1，跳转至步骤B)，继续迭代，直到满足预定义的数据精度；6)求得公平且最优的次要用户的速率设置效用函数参数α的值，并且将得到的最优对偶因子Ψ*一起带入对偶函数公式(11)得到公平且最优的次要系统效用。