[发明专利]一种D2D通信中基于布伦特方法的带有保障服务质量的能效优化发射功率控制方法

权利要求书

1.一种D2D通信中基于布伦特方法的带有保障服务质量的能效优化发射功率控制方法，其特征在于：所述控制方法包括以下步骤：

(1)在D2D通信中，设定最优化问题为：实现最大化D2D通信用户的能量效率同时保证包括D2D通信用户和蜂窝网用户的服务质量，采用控制D2D通信用户和蜂窝网用户的发射功率的方法，将该最优化问题描述为问题P1：

P1:maxp,qRd(p,q)q+h=1q+hlog2(1+qgddpgcd+n0)]]>

受限制于：

Rd(p,q)=log2(1+pgddqgcd+n0)≥Rdreq]]>

Pcmax≥P≥0]]>

Pdmax≥q≥0]]>

在问题P1中，各参数定义如下：

p：蜂窝网用户的发射功率；

q：D2D通信用户的发射功率；

h：D2D通信中发射机为其电路基本运作的功率消耗，发射功率除外；

蜂窝网用户的吞吐量要求；

D2D通信用户的吞吐量要求；

蜂窝网用户的最大发射功率；

D2D通信用户的最大发射功率；

g_dd：D2D通信中发射用户到D2D通信中接收用户之间的信道增益；

g_cB：蜂窝网用户到蜂窝网基站之间的信道增益；

g_dc：D2D通信中发射用户到蜂窝网基站之间的信道增益；

g_cd：蜂窝网用户到D2D通信中接收用户之间的信道增益；

R_c(p，q)：蜂窝网用户的上行链路的吞吐量；

R_d(p，q)：D2D通信用户的吞吐量；

n₀：背景噪声功率；

(2)利用公式将p用q代替，与此同时，引入q新的上下界变量与Q将问题P1等效地转化为一个对D2D通信用户发射功率的控制优化问题P2：

P2:max1q+hlog2(1+gddgcBqn0(gcB+αcgcd)+αcgdcgcdq)]]>

受限制于：

其中，

定义A＝g_cBg_dd，B＝α_cg_cdg_dc，C＝n₀(g_cB+α_cg_cd)，并定义函数将优化问题P2等效地转化为优化问题P3：

P3:maxQ‾≥q≥Q‾andz≥0z]]>

受限制于：F(z，q)≥0

在问题P2和问题P3中，各参数定义如下：

p：蜂窝网用户的发射功率；

q：D2D通信用户的发射功率；

h：D2D通信中发射机为其电路基本运作的功率消耗，发射功率除外；

蜂窝网用户的吞吐量要求；

D2D通信用户的吞吐量要求；

蜂窝网用户的最大发射功率；

D2D通信用户的最大发射功率；

g_dd：D2D通信中发射用户到D2D通信中接收用户之间的信道增益；

g_cB：蜂窝网用户到蜂窝网基站之间的信道增益；

g_dc：D2D通信中发射用户到蜂窝网基站之间的信道增益；

g_cd：蜂窝网用户到D2D通信中接收用户之间的信道增益；

n₀：背景噪声功率；

z：辅助待决定的控制变量；

(3)为求解问题P3，变形成如下子问题：

P3-A：q^*(z)＝arg max F(z，q)

受限制于：

然后利用q^*(z)，解决如下子问题得到{z^*，q^*}：

P3-B：max z

受限制于：F(z，q^*(z))≥0

z≥0

在子问题P3-A和子问题P3-B中，参数定义如下：

q^*：D2D通信中用户发射功率的最优值；

{z^*，q^*}：问题P3的最优解；

采用分层的思想对问题P3进行优化求解，步骤为：

步骤3.1：对于步骤(3)中的唯一最优值q^*(z)有：

q*(z)=Q‾,0<z<M‾Q‾,z>M‾-(A+2B)C+Δ2(A+B)B,M‾≤z≤M‾]]>

其中，A＝g_cBg_dd，B＝α_cg_cdg_dc，C＝n₀(g_cB+α_cg_cd)；

步骤3.2：问题P3-B中最优的情况在F(z，q^*(z))＝0下取得，另外定义D2D通信用户的发射机电路运作基本功耗h的两个上下临界值为：和

1)若h＜H，那么问题P3的最优解{z^*，q^*}为：q^*＝Q，从而获得问题P1的最优解为：D2D通信用户的最优发射功率q^*＝Q，蜂窝网用户的最优发射功率以及D2D通信用户的最大能量效率：

2)若那么问题P3的最优解{z^*，q^*}为：从而获得问题P1的最优解为：D2D通信用户的最优发射功率蜂窝网用户的最优发射功率以及D2D通信用户的最大能量效率：

3)若那么问题P3的最优解{z^*，q^*}通过如下过程求解：

步骤3.2.1：设置b为DU发送功率的上边界，a为下边界，z＝M，bx＝0.5(a+b)，ia＝min(a，b)，ib＝max(a，b)，v＝bx，w＝v，x＝v，e＝0，其中，q＝x，设置插值比例cgold，迭代精度epsilon，迭代次数IterTimes，为二分法计算误差的限度，初始化当前已迭代次数inter＝0；

步骤3.2.2：

步骤3.2.3：判断若inter≤IterTimes，则跳至步骤3.2.4，否则跳至步骤；

步骤3.2.4：设置xm＝0.5(ia+ib)，判断若abs(x-xm)＜epsilon·2-0.5(ib-ia)，则循环结束，跳至步骤3.2.12，否则跳至步骤3.2.5；

步骤3.2.5：判断若abs(e)≤epsilon，则跳至步骤3.2.7，否则，设置r＝(x-w)(fx-fv)，q＝(x-v)(fx-fw)，p＝(x-v)q-(x-w)r，q＝2(q-r)，其中，如果q大于零，则设置p＝-p，设置q＝abs(q)，etemp＝e，e＝d；

步骤3.2.6：判断若同时满足abs(p)＜abs(0.5·q·etemp)，p＞q(ia-x)，p＜q(ib-x)这三个条件，则设置d＝p/q，u＝x+d，其中，若u-ia＜epsilon·2并且ib-u＜epsilon·2，则设置d＝sign(xm-x)·epsilon，跳至步骤3.2.8，否则，跳至步骤3.2.7；

步骤3.2.7：当x≥xm时，设置e＝ia-x，反之，设置e＝ib-x，设置d＝cgold·e；

步骤3.2.8：判断若abs(d)≥epsilon，则设置u＝x+d，否则，设置u＝x+sign(d)·epsilon；

步骤3.2.9：设置q＝u，

步骤3.2.10：若fu≤fx，则判断u≥x，若是，则设置ia＝x，若否，则设置ib＝x，v＝w，fv＝fw，w＝x，fw＝fx，x＝u，fx＝fu，若fu＞fx，则判断u＜x，若是，则设置ia＝u，若否，则设置ib＝u，与此同时，如果fu≤fw并且w＝x，则设置v＝w，fv＝fw，w＝u，fw＝fu，如果fu≤fv并且v＝x以及v＝w，则设置v＝u，fv＝fu；

步骤3.2.11：inter＝inter+1，跳至步骤3.2.3；

步骤3.2.12：得到：当前最大值-fx，以及取到最大值时对应的x；

步骤3.2.13：如果则算法运行停止，输出：D2D通信用户的最优发射功率q^*＝x，蜂窝网用户的最优发射功率D2D通信用户的最大能量效率z^*＝z；否则执行步骤3.2.14；

步骤3.2.14：如果将z赋值给z，跳至步骤3.2.2；如果将z赋值给跳至步骤3.2.2。