1.一种大规模天线系统高能效资源优化方法,其特征在于该方法包括以下步骤:
1).设定初始发射功率设定传输速率约束罚因子σ1和干扰约束罚因子σ2,设定初始拉格朗日参数和其中,为用户i的初始发射功率,K表示用户数,为关于用户i的干扰约束的初始拉格朗日系数,为关于用户i的传输速率约束的初始拉格朗日系数;
2).计算在上述设定初始发射功率下对应的初始能效值q(0),
<mrow>
<msup>
<mi>q</mi>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
</msup>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>K</mi>
</munderover>
<msub>
<mi>log</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>N</mi>
<mi>t</mi>
</msub>
<msubsup>
<mi>p</mi>
<mi>i</mi>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
</msubsup>
</mrow>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
<mo>&NotEqual;</mo>
<mi>i</mi>
</mrow>
<mi>K</mi>
</munderover>
<msubsup>
<mi>p</mi>
<mi>j</mi>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
</msubsup>
<mo>+</mo>
<msup>
<mi>&sigma;</mi>
<mn>2</mn>
</msup>
<mo>/</mo>
<msub>
<mi>&beta;</mi>
<mi>i</mi>
</msub>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>K</mi>
</munderover>
<msubsup>
<mi>p</mi>
<mi>i</mi>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
</msubsup>
<mo>+</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>N</mi>
<mi>t</mi>
</msub>
</munderover>
<msub>
<mi>p</mi>
<mrow>
<mi>c</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
</mrow>
</mfrac>
</mrow>
其中,K为用户数,Nt为基站天线数,为用户i的初始发射功率,为用户j的初始发射功率,σ2为噪声功率谱密度,βi为用户i的信道大尺度衰落参数,pc,j表示天线j上的固定功率消耗;
3).针对用户i的发射功率进行迭代计算其中,k表示迭代次数,为用户i经过k+1次迭代后的发射功率,为用户i经过k次迭代后的发射功率,表示关于用户i的发射功率的表达式,迭代过程为:
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>p</mi>
<mi>i</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
<mo>=</mo>
<mi>T</mi>
<mrow>
<mo>(</mo>
<msubsup>
<mi>p</mi>
<mi>i</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mi>max</mi>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>,</mo>
<msubsup>
<mi>m</mi>
<mi>i</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>-</mo>
<msub>
<mi>&sigma;</mi>
<mn>2</mn>
</msub>
<mo>(</mo>
<mrow>
<msub>
<mover>
<mi>r</mi>
<mo>&OverBar;</mo>
</mover>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<msubsup>
<mi>p</mi>
<mi>i</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>M</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
<mfenced open = "{" close = "}">
<mtable>
<mtr>
<mtd>
<mrow>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
<mo>&NotEqual;</mo>
<mi>i</mi>
</mrow>
<mi>K</mi>
</munderover>
<mfrac>
<mn>1</mn>
<mrow>
<msub>
<mi>I</mi>
<mi>j</mi>
</msub>
<mrow>
<mo>(</mo>
<msubsup>
<mi>p</mi>
<mi>i</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msup>
<mi>&sigma;</mi>
<mn>2</mn>
</msup>
<mo>/</mo>
<msub>
<mi>&beta;</mi>
<mi>i</mi>
</msub>
</mrow>
</mfrac>
<mo>+</mo>
<msup>
<mi>q</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</msup>
<mi>ln</mi>
<mn>2</mn>
<mo>+</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>ln</mi>
<mn>2</mn>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
<mo>&NotEqual;</mo>
<mi>i</mi>
</mrow>
<mi>K</mi>
</munderover>
<mi>max</mi>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>,</mo>
<msubsup>
<mi>&omega;</mi>
<mi>j</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>-</mo>
<msub>
<mi>&sigma;</mi>
<mn>1</mn>
</msub>
<mo>(</mo>
<mrow>
<msub>
<mi>I</mi>
<mrow>
<mi>M</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>I</mi>
<mi>j</mi>
</msub>
<mrow>
<mo>(</mo>
<msubsup>
<mi>p</mi>
<mi>i</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mfrac>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中,为用户i经过k+1次迭代后的发射功率,分别为用户i和用户j第k次迭代后的拉格朗日参数,且σ1、σ2则分别为干扰约束及传输速率约束对应的罚因子,表示用户i经过k次迭代后的实际传输速率下限,其为关于用户i的发射功率的表达式,RM,i表示用户i的最低传输速率,K为用户数,表示用户j经过k次迭代后受到的实时干扰,其为关于用户i的发射功率的表达式,σ2为噪声功率谱密度,βi为用户i的信道大尺度衰落参数,q(k)为第k次迭代后的实际能效,IM,j为用户j可容忍的最大接收干扰;
4).对用户i的发射功率进行调整其中,为用户i经过k+1次迭代且进行调整后的发射功率,为用户i经过k+1次迭代后的发射功率,pT,i为用户i的最大发射功率;
5).用调整后的发射功率来计算能效值q(k+1)和拉格朗日参数和其中,q(k+1)为第k+1次迭代后的能效值,和分别为第k+1次迭代后的用户i的拉格朗日参数;
6).计算效用函数的值,对给定的量级在10-1之内的正常数ε,若|UEE|≤ε,则整个迭代计算则结束,否则,令k=k+1,返回步骤3)继续进行迭代计算,其中,UEE为效用函数,K为用户数,表示用户i的实时传输速率下限,q(k +1)为第k+1次迭代后的能效值,为用户i经过k+1次迭代且进行调整后的发射功率,Nt为基站天线数,pc,j表示天线j上的固定功率消耗,是用户i的下一次迭代的初始发射功率。