搜索 分析 新世界 法规 图书 网址导航 更多
高级用户登录 | 登录 | |

多用户OFDM中继系统能量效率优化算法
审中-实审

申请号:201610073946.X 申请日:2016-01-29
摘要:本发明公开了一种多用户OFDM中继系统能量效率优化算法,在满足系统最小传输速率的前提下,采用源在TS2也发送数据的策略以提供分集增益,构建了包含子载波配对以及功率分配的系统能效优化模型。该模型是一个典型的非线性混合整数规划问题,这类问题通常求解非常复杂。我们通过定义等效信道增益,将复杂的最优化问题简化为准凸规划问题,并通过Dinkelbach方法、匈牙利算法以及次梯度算法来求得最优解,大大减小了计算复杂度。仿真结果表明,我们提出的算法能在保证满足系统最小传输速率要求的情况下,提高系统的能量效率,并且具有很好的全局收敛性。
申请人: 安徽工程大学
地址: 241000 安徽省芜湖市北京中路8号
发明(设计)人: 梁广俊
主分类号: H04W52/26(2009.01)I
分类号: H04W52/26(2009.01)I H04W52/34(2009.01)I H04W52/46(2009.01)I H04L5/00(2006.01)I
  • 法律状态
2016-07-27  实质审查的生效IPC(主分类):H04W 52/26申请日:20160129
2016-06-29  公开
注:本法律状态信息仅供参考,即时准确的法律状态信息须到国家知识产权局办理专利登记簿副本。
  • 其他信息
主权项  一种多用户OFDM中继系统能量效率优化算法,其特征在于,所述算法包括如下步骤:步骤1:建立系统优化模型:步骤1.1:系统参数设定,一个两跳OFDM单中继网络,基站通过单个中继发送信息给K个用户,通信模式采用TDD模式,系统共划分为N个子载波,每一个子载波的都占用相同的带宽并经历相互独立的瑞利衰落,假设每一个中继节点可以获得自己前向和后向的瞬时信道信息,并且所有的信道瞬时信息在源节点S和K个用户节点处都是可知的,定义分别表示第一跳被占用的子载波i在基站和中继之间的信道增益,第二跳被占用的子载波i在中继和用户k之间的信道增益,第一跳直接链路下占用子载波i基站和用户k之间的信道增益以及第二跳直接链路下占用的子载波i基站和用户k之间的信道增益,中继和用户k在每个子载波上的高斯白噪声(AWGN)的功率分别为步骤1.2:优化指标设定,包括a)信道分配如果第一跳基站通过子载波m将信号发送给中继,中继通过子载波n将信号发送给用户k,则表示路径P(m,n,k)存在,那么代表子载波m与n配对、子载波对(m,n)?用户分配的三维分配因子可以表示为:任意用户可以分配到多个子载波对,而为防止自干扰,对于任意子载波对只能被一个用户占用,可以表示为<mrow><munderover><mi>&Sigma;</mi><mrow><mi>n</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>K</mi></munderover><msub><mi>&phi;</mi><mrow><mi>m</mi><mo>,</mo><mi>n</mi><mo>,</mo><mi>k</mi></mrow></msub><mo>=</mo><mn>1</mn><mo>,</mo><mo>&ForAll;</mo><mi>m</mi></mrow>m与n;b)功率分配对于任意路径P(m,n,k),基站第一时隙、第二时隙的发射功率和中继的发射功率分别表示为任意路径P(m,n,k)的最优功率分配可以描述为整个系统的功率限制可以表示为:其中Ptotal表示基站和中继两者总的功率之和;c)中继策略中继采用DF转发方式,路径P(m,n,k)上最大接收速率可以表示为:<mrow><mi>R</mi><mrow><mo>(</mo><mi>m</mi><mo>,</mo><mi>n</mi><mo>,</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mi>m</mi><mi>i</mi><mi>n</mi><mo>{</mo><msub><mi>log</mi><mn>2</mn></msub><mrow><mo>(</mo><mn>1</mn><mo>+</mo><msub><mi>a</mi><mi>m</mi></msub><msubsup><mi>P</mi><mrow><mi>m</mi><mi>n</mi><mi>k</mi></mrow><mrow><mi>s</mi><mn>1</mn></mrow></msubsup><mo>)</mo></mrow><mo>,</mo><msub><mi>log</mi><mn>2</mn></msub><mrow><mo>(</mo><mn>1</mn><mo>+</mo><msub><mi>b</mi><mrow><mi>n</mi><mi>k</mi></mrow></msub><msubsup><mi>P</mi><mrow><mi>m</mi><mi>n</mi><mi>k</mi></mrow><mi>r</mi></msubsup><mo>+</mo><msub><mi>c</mi><mrow><mi>m</mi><mi>k</mi></mrow></msub><msubsup><mi>P</mi><mrow><mi>m</mi><mi>n</mi><mi>k</mi></mrow><mrow><mi>s</mi><mn>1</mn></mrow></msubsup><mo>+</mo><msub><mi>a</mi><mrow><mi>n</mi><mi>k</mi></mrow></msub><msubsup><mi>P</mi><mrow><mi>m</mi><mi>n</mi><mi>k</mi></mrow><mrow><mi>s</mi><mn>2</mn></mrow></msubsup><mo>)</mo></mrow><mo>}</mo><mo>,</mo></mrow>其中<mrow><msub><mi>a</mi><mi>m</mi></msub><mo>=</mo><mfrac><msup><mrow><mo>|</mo><msubsup><mi>h</mi><mi>m</mi><mrow><mi>s</mi><mi>r</mi></mrow></msubsup><mo>|</mo></mrow><mn>2</mn></msup><msubsup><mi>&sigma;</mi><mi>r</mi><mn>2</mn></msubsup></mfrac><mo>,</mo><msub><mi>b</mi><mrow><mi>n</mi><mi>k</mi></mrow></msub><mo>=</mo><mfrac><msup><mrow><mo>|</mo><msubsup><mi>h</mi><mi>n</mi><mrow><mi>r</mi><mi>k</mi></mrow></msubsup><mo>|</mo></mrow><mn>2</mn></msup><msubsup><mi>&sigma;</mi><mi>k</mi><mn>2</mn></msubsup></mfrac><mo>,</mo><msub><mi>c</mi><mrow><mi>m</mi><mi>k</mi></mrow></msub><mo>=</mo><mfrac><msup><mrow><mo>|</mo><msubsup><mi>h</mi><mi>m</mi><mrow><mi>s</mi><mi>k</mi><mn>1</mn></mrow></msubsup><mo>|</mo></mrow><mn>2</mn></msup><msubsup><mi>&sigma;</mi><mi>k</mi><mn>2</mn></msubsup></mfrac><mo>,</mo><msub><mi>a</mi><mrow><mi>n</mi><mi>k</mi></mrow></msub><mo>=</mo><mfrac><msup><mrow><mo>|</mo><msubsup><mi>h</mi><mi>n</mi><mrow><mi>s</mi><mi>k</mi><mn>2</mn></mrow></msubsup><mo>|</mo></mrow><mn>2</mn></msup><msubsup><mi>&sigma;</mi><mi>k</mi><mn>2</mn></msubsup></mfrac><mo>,</mo></mrow>1/2的意义是1个信息的传输需要2个时隙;步骤1.3:优化模型设定,当给定路径P(m,n,k)上总功率为Pmnk时,最大接收速率可以转化为:其中amnk表示路径P(m,n,k)上的等效增益,表示为整个系统的速率可以表示为:<mrow><msub><mi>R</mi><mrow><mi>t</mi><mi>o</mi><mi>t</mi><mi>a</mi><mi>l</mi></mrow></msub><mrow><mo>(</mo><mi>&Phi;</mi><mo>,</mo><mi>P</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>K</mi></munderover><munderover><mi>&Sigma;</mi><mrow><mi>m</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><munderover><mi>&Sigma;</mi><mrow><mi>n</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>&phi;</mi><mrow><mi>m</mi><mi>n</mi><mi>k</mi></mrow></msub><msub><mi>log</mi><mn>2</mn></msub><mrow><mo>(</mo><mn>1</mn><mo>+</mo><msub><mi>a</mi><mrow><mi>m</mi><mi>n</mi><mi>k</mi></mrow></msub><msub><mi>P</mi><mrow><mi>m</mi><mi>n</mi><mi>k</mi></mrow></msub><mo>)</mo></mrow><mo>,</mo></mrow>其中P=[Pmnk]N×N×K,Φ=[φmnk]N×N×K且功率Pmnk和因子φmnk分别满足先前的约束,相应地,基站和中继所消耗的功率和以及系统由于传输所消耗的总的功率可以分别表示为:<mrow><msub><mi>P</mi><mrow><mi>u</mi><mi>s</mi><mi>e</mi><mi>r</mi><mi>s</mi></mrow></msub><mrow><mo>(</mo><mi>&Phi;</mi><mo>,</mo><mi>P</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>m</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><munderover><mi>&Sigma;</mi><mrow><mi>n</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>K</mi></munderover><msub><mi>&phi;</mi><mrow><mi>m</mi><mi>n</mi><mi>k</mi></mrow></msub><msub><mi>P</mi><mrow><mi>m</mi><mi>n</mi><mi>k</mi></mrow></msub></mrow>其中PC是固定为常数的环路电流,表示功率转换效率的倒数,最大化系统平均能量效率的最优化模型可以表示为:<mfenced open="" close=""><mtable><mtr><mtd><munder><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow><mrow><mi>&Phi;</mi><mo>,</mo><mi>P</mi></mrow></munder></mtd><mtd><mrow><msub><mi>&eta;</mi><mi>E</mi></msub><mrow><mo>(</mo><mi>&Phi;</mi><mo>,</mo><mi>P</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>R</mi><mrow><mi>t</mi><mi>o</mi><mi>t</mi><mi>a</mi><mi>l</mi></mrow></msub><mrow><mo>(</mo><mi>&Phi;</mi><mo>,</mo><mi>P</mi><mo>)</mo></mrow></mrow><mrow><msub><mi>P</mi><mrow><mi>o</mi><mi>t</mi><mi>a</mi><mi>l</mi></mrow></msub><mrow><mo>(</mo><mi>&Phi;</mi><mo>,</mo><mi>P</mi><mo>)</mo></mrow></mrow></mfrac></mrow></mtd></mtr></mtable></mfenced>s.t.?(3),(4)and(5)<mrow><msub><mi>&phi;</mi><mrow><mi>m</mi><mi>n</mi><mi>k</mi></mrow></msub><mo>=</mo><mo>{</mo><mn>0</mn><mo>,</mo><mn>1</mn><mo>}</mo><mo>,</mo><msub><mi>P</mi><mrow><mi>m</mi><mi>n</mi><mi>k</mi></mrow></msub><mo>&GreaterEqual;</mo><mn>0</mn><mo>,</mo><mo>&ForAll;</mo><mi>m</mi><mo>,</mo><mi>n</mi><mo>,</mo><mi>k</mi></mrow>Rtotal(Φ,S)≥Rreq;步骤2:获取瞬时信道信息,目的节点通过训练序列获得各信道的瞬时信道信息;根据式计算子载波用户对(m,n,k)的等效信道增益amnk,并使用黄金分割法在区间(0,1]求得使得等效信道增益amnk取得最大值;步骤3:初始化外循环参数,包括:初始化外循环迭代终止常数εouter以及外循环最大允许迭代次数q0=0,n=0;步骤4:内循环算法主体部分;步骤4.1:选取合适的拉格朗日因子λ的初始值λ(0);步骤4.2:定义决策矩阵Φ={φmnk},φmnk=1表示第m个子载波和第n个子载波对进行配对并指派给用户k,φmnk=0表示其它情况,计算子载波配对决策因子A(i,j),然后通过子载波用户对三维分配算法得到决策矩阵Φ;步骤4.3:根据拉格朗日因子算法计算(m,n,k)计算Smnk以及步骤4.4:根据次梯度算法更新拉格朗日因子;步骤4.5:重复步骤4.1至步骤4.4直到相邻两次的拉格朗日因子的差值的绝对值小于εinner或内循环迭代次数大于同时相应的决策矩阵Φ也被确定;步骤5:外循环算法收敛判断;步骤5.1:计算<mrow><msub><mi>R</mi><mrow><mi>t</mi><mi>o</mi><mi>t</mi><mi>a</mi><mi>l</mi></mrow></msub><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>K</mi></munderover><munderover><mi>&Sigma;</mi><mrow><mi>m</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><munderover><mi>&Sigma;</mi><mrow><mi>n</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>&omega;</mi><mi>k</mi></msub><msub><mi>&phi;</mi><mrow><mi>m</mi><mi>n</mi><mi>k</mi></mrow></msub><msub><mi>log</mi><mn>2</mn></msub><mrow><mo>(</mo><mn>1</mn><mo>+</mo><msub><mi>a</mi><mrow><mi>m</mi><mi>n</mi><mi>k</mi></mrow></msub><msub><mi>P</mi><mrow><mi>m</mi><mi>n</mi><mi>k</mi></mrow></msub><mo>)</mo></mrow></mrow>以及然后根据公式来计算此时对应的能量效率;步骤5.2:如果相邻两次计算得到的能量效率的差值的绝对值小于εout或外循环迭代次数大于转入步骤6,否则重复步骤3至步骤5;步骤6:输出算法计算结果,针对满足φmnk=1的所有子载波用户对(m,n,k)计算以及以获得功率分配信息。
公开号  105722206A
公开日  2016-06-29
专利代理机构  合肥鼎途知识产权代理事务所(普通合伙) 34122
代理人  叶丹
颁证日  
优先权  
国际申请  
国际公布  
进入国家日期