-
作为无线通信网络管理的关键问题之一,网络节能或绿色通信受到了广泛的关注[1]。最常用的节能方法是功率分配,它在一定的服务质量(quality of service,QoS)限制条件下,通过优化各链路的传输功率使网络的功耗最小[2-3]。随着无线网络的快速发展,网络管理变得更加复杂和更具挑战,单纯依赖功率分配来实现网络节能越来越困难。事实上,人们普遍认为应该将功率分配和其他的网络管理手段(如波束形成[4-6]、用户调度[7]和协作处理[8-11]等)结合起来,共同实现网络节能。
除了上述方法外,接入控制也是一种非常有效的网络管理手段[12-16]。它通过有选择地接入某些用户来为网络优化提供更多的灵活性和更大的性能提升空间。由于问题本身的复杂性,接入控制往往作为独立的网络管理问题被单独研究,如分别以负载均衡[12]、系统容量[13]或阻塞概率[14]为目标制定相应的接入控制方案等。近年来,研究人员开始考虑将接入控制和其他网络管理方法联合起来进行网络优化。如文献[15]利用接入控制和波束形成等方法进行联合优化,提升网络频谱效率。文献[16]将接入控制和功率分配结合起来降低网络的功耗。
本文研究了与文献[16]相似的问题,即通过功率分配和接入控制的联合优化来实现网络节能,但采用了不同的接入控制策略和求解算法。文献[16]要求接入尽可能多的用户,并为此设计了基于交替紧缩搜索和功率分配的求解算法。具体而言,它首先固定发射功率,逐个排除无法满足QoS条件的用户;然后,基于固定的接入用户求解传统功率分配问题,更新发射功率;以此类推,直到交错优化过程收敛。与文献[16]不同,本文从稀疏优化的角度来描述功率分配和接入控制的联合优化问题,并在接入控制中采用一种更灵活的软策略。与之相对,文献[16]强制要求接入所有QoS可能被满足的用户,这可以看作是一种硬策略。软策略的灵活性主要体现在两方面:1) 当网络无法满足所有用户的QoS条件时,接入控制促使网络转而去服务尽可能多的用户,而非简单地拒绝所有用户。这提高了网络服务时间和用户接入概率。2) 在此基础上可以通过权重因子控制接入用户的数量。如果某个用户处于较差的信道环境,网络需要大大增加对应链路的发射功率以满足该用户的QoS条件。而这又对其他链路造成干扰,相应地其他链路也会增大发射功率,最终显著增加整个网络的功耗。软策略可以有选择地挑选用户接入,以相对较低的功耗代价满足接入用户的QoS条件,提升了功效。显然,与绝对功耗相比,功效是更为合理的节能指标。
算法方面,本文采用分布式求解策略。随着无线网络规模的增加,优化问题维数越来越大,分布式求解策略显然更具吸引力,即依赖于局部信道信息,或进行简单数据交换后,各基站可以自行决定是否接入对应用户。基于这一思路,本文对原始优化问题进行转换,使用交错方向乘子法(alternating direction method of multipliers,ADMM[5, 17])将其分解成几个简单子问题进行分布式迭代求解。求解算法具有低运算量,每一步都有闭合解。
本文的主要工作为:1) 在传统功率分配问题中引入接入控制,进行联合优化实现网络节能,并通过接入控制软策略和稀疏优化完成用户分类和挑选,提高网络的传输功效。2) 基于ADMM框架设计分布式求解算法,在进行少量数据交换后,各链路自行决定是否建立连接,且分布式算法复杂度低,每一步都有闭合解。
-
考虑一个由多个单天线基站和用户组成的无线通信网络。简便起见,假设基站和用户的数量均为M,且用户m指定接入到基站m,m = 1,2,…,M。考虑下行链路,pm表示基站m的实际发射功率,Pm表示该基站发射功率上限。定义p = [p1,p2,…,pM]T为网络的实际发射功率矢量。
在此基础上,假设hn,m为基站m与用户n之间的传输增益,m,n = 1,2,…,M。以用户端的信号干扰噪声比(signal to interference plus noise ratio,SINR)来衡量网络服务质量QoS,假设用户m的QoS门限为gm,则QoS条件表示为:
$$\text{SIN}{{\text{R}}_{m}}=\frac{{{p}_{m}}{{h}_{m,m}}}{{{\sigma }^{2}}+\sum\limits_{n\ne m}{{{p}_{n}}{{h}_{m,n}}}}\ge {{\gamma }_{m}}\quad \forall \ m$$ (1) 式中,σ2表示用户的接收噪声功率。
-
传统的网络节能问题可以描述为:
$$\begin{matrix} ({{\text{P}}_{1}})\quad \underset{\mathbf{p}}{\mathop{\min }}\,\ \sum\limits_{m=1}^{M}{{{p}_{m}}} \\ \text{s}\text{.t}\text{.}\ \ \text{ SIN}{{\text{R}}_{m}}\ge {{\gamma }_{m}},\ \ 0\le {{p}_{m}}\le {{P}_{m}},\ \ \forall \ m \\ \end{matrix}$$ 在实际应用中该策略有两点不足。首先,由于发射功率存在上限,网络可能无法满足所有用户的QoS条件,导致(P1)没有可行解,整个网络停止工作。这种由于少量用户QoS无法满足,导致所有用户都不能接入的处理方式并不合理。更具实际意义的策略是,在网络无法满足所有用户的QoS条件时,应该尝试服务尽可能多的用户。其次,与绝对功耗相比,功效更能体现网络的节能程度,它描述了消耗单位功率实现的传输速率。如果某个用户位于恶劣的信道环境中,网络需要显著增加传输功率以满足该用户的QoS条件,并对抗由此给其他用户造成的干扰,会降低整个网络的功效。因此,需要在接入用户数和功效之间折衷,有选择地挑选某些用户优先接入。
基于上述考虑,在(P1)中融合接入控制,并以能否满足QoS条件作为用户接入依据。式等价为:
$${{\gamma }_{m}}\left( {{\sigma }^{2}}+\sum\limits_{n\ne m}{{{p}_{n}}{{h}_{m,n}}} \right)-{{p}_{m}}{{h}_{m,m}}\le 0\quad \forall \ m$$ (2) 如果用户m的QoS条件满足,则式成立;反之式不成立,它的左边应为正数。为了覆盖这两种情况,引入变量s = [s1,s2,…,sM]T,sm≥0,m =1,2,…,M,得到如下功率分配和接入控制联合优化问题,有:
$$\begin{matrix} ({{\text{P}}_{2}})\quad \underset{\mathbf{p},\ \mathbf{s}}{\mathop{\min }}\,\ \sum\limits_{m=1}^{M}{{{p}_{m}}}+\lambda ||\mathbf{s}|{{|}_{0}} \\ \,\text{s}\text{.t}\text{.}\ \text{ }\ {{\gamma }_{m}}\left( {{\sigma }^{2}}+\sum\limits_{n\ne m}{{{p}_{n}}{{h}_{m,n}}} \right)-{{p}_{m}}{{h}_{m,m}}={{s}_{m}} \\ 0\le {{p}_{m}}\le {{P}_{m}},\ \ {{s}_{m}}\ge 0,\ \ \forall \ m \\ \end{matrix}$$ 对于优化后的结果,如果sm = 0,有SINRm = γm;如果sm > 0,则SINRm < γm,用户m的QoS无法满足,不能接入。矢量s的稀疏度||s||0指示了接入的用户数量;λ > 0为权重因子,用于平衡接入用户数和传输功率,体现了接入控制软策略。
稀疏惩罚项中含有非凸的l0范数,常用方法是将它松弛成l1范数,即:
$$\begin{matrix} ({{\text{P}}_{3}})\quad \underset{\mathbf{p},\ \mathbf{s}}{\mathop{\min }}\,\ \sum\limits_{m=1}^{M}{{{p}_{m}}}+\sum\limits_{m=1}^{M}{{{\lambda }_{m}}}\left| {{s}_{m}} \right| \\ \text{s}\text{.t}\text{.}\ \text{ }\ {{\gamma }_{m}}\left( {{\sigma }^{2}}+\sum\limits_{n\ne m}{{{p}_{n}}{{h}_{m,n}}} \right)-{{p}_{m}}{{h}_{m,m}}={{s}_{m}} \\ 0\le {{p}_{m}}\le {{P}_{m}},\ \ {{s}_{m}}\ge 0,\ \ \forall \ m \\ \end{matrix}$$ 在(P3)中为各个|sm|分配了不同的权重因子lm > 0,使问题的处理更加灵活。
-
(P3)是一个凸问题,可以利用软件包(如CVX[18]等)直接进行集中式求解。另一方面,随着用户和基站的增加,网络规模越来越大,(P3)的维数变高。集中式地求解高维优化问题增大了主处理器的负担。在这种情况下,分布式求解算法更具吸引力,其目的是依赖局部信息,在简单数据交换后,基站能自行决定是否接入对应用户。但是,由于功率变量pm和接入变量sm在QoS条件中耦合,直接对(P3)分布式求解十分困难。
-
为了书写简便,将(P3)写成矩阵和矢量的形式,则有:
$$\begin{matrix} ({{\text{P}}_{4}})\quad \underset{\mathbf{p},\ \mathbf{s}}{\mathop{\min }}\,\ \sum\limits_{m=1}^{M}{{{p}_{m}}}+\sum\limits_{m=1}^{M}{{{\lambda }_{m}}\left| {{s}_{m}} \right|} \\ \text{s}\text{.t}\text{.}\ \ \text{ }\mathbf{Bp}+\mathbf{d}=\mathbf{s} \\ 0\le {{p}_{m}}\le {{P}_{m}},\ \ \forall \ m \\ \end{matrix}$$ 式中,矩阵B RM x M和矢量d RM x 1定义如下:
$$\begin{matrix} \mathbf{B}=[{{\mathbf{b}}_{1}},\ \ {{\mathbf{b}}_{2}},\ \ \cdots ,\ \ {{\mathbf{b}}_{M}}]={{[\mathbf{\bar{b}}_{1}^{\text{T}},\ \ \mathbf{\bar{b}}_{2}^{\text{T}},\ \ \cdots ,\ \ \mathbf{\bar{b}}_{M}^{\text{T}}]}^{\text{T}}}= \\ \left[ \begin{matrix} -{{h}_{1,1}} & {{\gamma }_{1}}{{h}_{1,2}} & \cdots & {{\gamma }_{1}}{{h}_{1,M}} \\ {{\gamma }_{2}}{{h}_{2,1}} & -{{h}_{2,2}} & \cdots & {{\gamma }_{2}}{{h}_{2,M}} \\ \vdots & \vdots & \ddots & \vdots \\ {{\gamma }_{M}}{{h}_{M,1}} & {{\gamma }_{M}}{{h}_{2,M}} & \cdots & -{{h}_{M,M}} \\ \end{matrix} \right] \\ \mathbf{d}={{[{{\gamma }_{1}}{{\sigma }^{2}},\ \ {{\gamma }_{2}}{{\sigma }^{2}},\ \ \cdots ,\ \ {{\gamma }_{M}}{{\sigma }^{2}}]}^{\text{T}}} \\ \end{matrix}$$ (3) 式中,bm和bm分别为B的第m个列矢量和行矢量。
与(P3)相比,(P4)去掉了sm的非负限制。事实上,去掉限制sm≥0并不影响(P3)的求解,或者说(P4)的最优解自动满足非负性。利用反证法,假设 $({{\mathbf{p}}^{*}},\ {{\mathbf{s}}^{*}})$ 是(P4)的最优解,且存在sk*<0,1≤k≤M,即:
$${{\gamma }_{k}}\left( {{\sigma }^{2}}+\sum\limits_{n\ne k}{p_{n}^{*}{{h}_{k,n}}} \right)-p_{k}^{*}{{h}_{k,k}}=s_{k}^{*}<0$$ (4) 那么,总可以找到 $0\le \bar{p}_{k}^{*}<p_{k}^{*}$ 和 $\bar{s}_{k}^{*}=0$ ,使:
$${{\gamma }_{k}}\left( {{\sigma }^{2}}+\sum\limits_{n\ne k}{p_{n}^{*}{{h}_{k,n}}} \right)-\bar{p}_{k}^{*}{{h}_{k,k}}=\bar{s}_{k}^{*}=0$$ (5) 对于用户m,m ≠ k,用户k的传输功率由pk*降为pk*,相当于减少了干扰,因此其QoS条件仍然满足。即将(p*,s*)中的(pk*,sk*)换成(pk*,sk*),可以得到一个新的矢量对,表示为(p*,s*),它仍然位于(P4)的可行集内。显然,(p*,s*)对应的目标函数值更小,因此,(p*,s*)不是(P4)的最优解,推出矛盾。
-
使用增广拉格朗日(augmented lagrangian)方法[17],(P4)的解可以通过求解(P5)获得,有:
$$\begin{align} & ({{\text{P}}_{5}})\quad \underset{\mathbf{p},\ \mathbf{s},\ \mathbf{\mu }}{\mathop{\min }}\,\ {{L}_{c}}(\mathbf{p},\ \mathbf{s},\ \mathbf{\mu }) \\ & \text{s}\text{.t}\text{.}\ \ \text{ }0\le {{p}_{m}}\le {{P}_{m}},\ \ \forall \ m \\ \end{align}$$ 式中,c > 0为惩罚因子;μ = [μ1,μ2,…,μM]T为限制条件Bp + d = s对应的拉格朗日乘子;Lc(p,s,μ)为增广拉格朗日函数,定义为:
$${{\mathbf{\mu }}^{\text{T}}}(\mathbf{Bp}+\mathbf{d}-\mathbf{s})+\frac{c}{2}\left\| \mathbf{Bp}+\mathbf{d}-\mathbf{s} \right\|_{2}^{2}$$ (6) 于是,可以将变量{p,s}分成{p}和{s},在此基础上利用如下的ADMM框架[17]迭代求解(P5),有:
$$\left\{ \begin{align} & \mathbf{p}(l+1)=\arg \underset{\mathbf{p}}{\mathop{\min }}\,\ {{L}_{c}}(\mathbf{p},\ \ \mathbf{s}(l),\ \ \mathbf{\mu }(l)) \\ & \text{s}\text{.t}\text{.}\ \,\text{ }0\le {{p}_{m}}\le {{P}_{m}},\ \ \forall \ m \\ & \mathbf{s}(l+1)=\arg \underset{\mathbf{s}}{\mathop{\min }}\,\ {{L}_{c}}(\mathbf{p}(l+1),\ \ \mathbf{s},\ \ \mathbf{\mu }(l)) \\ & \mathbf{\mu }(l+1)=\mathbf{\mu }(l)+c[\mathbf{Bp}(l+1)+\mathbf{d}-\mathbf{s}(l+1)] \\ \end{align} \right.$$ (7) 式中,l为迭代序号。由此变量p和变量s可以分开求解。进一步地,更新p、s和μ的过程还可以通过分布式的方式完成。
-
更新p等效于求解如下问题,有:
$$\begin{matrix} \underset{\mathbf{p}}{\mathop{\min }}\,\ {{L}_{c}}(\mathbf{p},\ \ \mathbf{s}(l),\ \ \mathbf{\mu }(l)) \\ \ \text{s}\text{.t}\text{.}\ \text{ }\,0\le {{p}_{m}}\le {{P}_{m}},\ \ \forall \ m \\ \end{matrix}$$ 由于hn,m为独立分布的随机信道,矩阵B通常是满秩的,因此Lc(p,s(l),μ(l))为关于p的严格凸函数。此外,各个pm具有独立的可行集。以上条件确保了坐标下降(coordinate descend,CD)法的收敛性[19]。因此,可以使用CD法求解该问题,即依次求解pm,m= 1,2,…,M,直到p收敛。在求解pm时,将其他的pn,n ≠ m,作为已知参数。
具体而言,以k为CD法的迭代序号,pm(l,k)为第l轮ADMM迭代过程中,第k次CD迭代后基站m的传输功率,第k+1次CD迭代过程如下:
1) 引入临时变量g = Bp(l,k) + d - s(l);
2) 开始迭代,令m =1,2,…,M,重复如下操作:
① 令g = g - pm(l,k)bm;
② 求解如下优化问题,获得pm(l,k+1),有:
$$\begin{matrix} \underset{{{p}_{m}}}{\mathop{\min }}\,\ \ {{p}_{m}}+{{\mathbf{\mu }}^{\text{T}}}(l)\left( {{p}_{m}}{{\mathbf{b}}_{m}}+\mathbf{g} \right)\,+\frac{c}{2}\left\| {{p}_{m}}{{\mathbf{b}}_{m}}+\mathbf{g} \right\|_{2}^{2} \\ \,\text{s}\text{.t}\text{.}\ \ \ 0\le {{p}_{m}}\le {{P}_{m}} \\ \end{matrix}$$ 利用一阶最优条件,求得pm(l,k+1)的最优解为:
$${{p}_{m}}(l,\ k+1)=\text{Proj}{{\,}_{[0,\ {{P}_{m}}]}}\,\left( \frac{-1-\mathbf{b}_{m}^{\text{T}}[\mathbf{\mu }(l)+c\mathbf{g}]}{c\mathbf{b}_{m}^{\text{T}}{{\mathbf{b}}_{m}}} \right)$$ (8) 式中, $\text{Proj}{{\,}_{[0,\ {{P}_{m}}]}}$ (∙)表示在[0,Pm]上的投影运算。
③ 更新g = g +(l,k + 1)bm;
3) 如果收敛,输出p(l + 1) = p(l,k);否则,k = k+1,回到步骤2);
因此,更新p的工作可以由各个基站依次分布式完成,在当前基站完成更新后,只需要将g传输给下一个基站即可。
-
更新s等效于求解如下的无限制优化问题,有:
$$\underset{\mathbf{s}}{\mathop{\min }}\,\ {{L}_{c}}(\mathbf{p}(l+1),\ \ \mathbf{s},\ \ \mathbf{\mu }(l))$$ 它可以分解成M个关于sm的独立子问题,有:
$$\underset{{{s}_{m}}}{\mathop{\min }}\,\ \left\{ \begin{matrix} {{\lambda }_{m}}\left| {{s}_{m}} \right|+{{\mu }_{m}}(l)[\ {{{\mathbf{\bar{b}}}}_{m}}\mathbf{p}(l+1)+{{d}_{m}}-{{s}_{m}}]+ \\ \frac{c}{2}{{[\ {{{\mathbf{\bar{b}}}}_{m}}\mathbf{p}(l+1)+{{d}_{m}}-{{s}_{m}}]}^{2}} \\ \end{matrix} \right\}$$ 同样,利用一阶最优条件,最优解sm(l + 1)满足:
$${{\mu }_{m}}(l)-c[{{s}_{m}}(l+1)-{{\mathbf{\bar{b}}}_{m}}\mathbf{p}(l+1)-{{d}_{m}}]\in {{\lambda }_{m}}\partial \left| {{s}_{m}}(l+1) \right|$$ (9) 式中,∂|sm(l + 1)|表示非连续的绝对值函数|sm(l + 1)|的次梯度(sub-gradient)[20],表示为:
$$\partial \left| {{s}_{m}}(l+1) \right|=\left\{ \begin{align} & \frac{{{s}_{m}}(l+1)}{\left| {{s}_{m}}(l+1) \right|}\quad \ \,\,\,{{s}_{m}}(l+1)\ne 0 \\ & \left\{ x\left| \ \right.|x|\,\,\le 1 \right\}\ \ \ {{s}_{m}}(l+1)=0 \\ \end{align} \right.$$ (10) 为了表述简洁,令:
$${{y}_{m}}={{\mu }_{m}}(l)+c{{\mathbf{\bar{b}}}_{m}}\mathbf{p}(l+1)+c{{d}_{m}}$$ (11) 将式和式代入式(9),可得:
$${{s}_{m}}(l+1)=\left\{ \begin{align} & 0\text{ }\quad \ \left| {{y}_{m}} \right|\le {{\lambda }_{m}} \\ & \frac{(\left| {{y}_{m}} \right|-{{\lambda }_{m}}){{y}_{m}}}{c\left| {{y}_{m}} \right|}\text{ }\quad \ \\ \end{align} \right.$$ (12) 因此,可以通过式更新sm,m =1,2,…,M。上述过程可以分布式进行,即各基站基于式(11)并行地计算ym,并根据|ym|的大小更新sm,由此自行决定是否接入用户m。另外,式表明了s的稀疏性,并且lm越大,表明稀疏惩罚项的权重越大,sm越趋向为0,用户m接入的概率也越大。
-
利用式更新μ,它同样可以分布式进行,即独立地计算μm(l + 1),m = 1,2,…,M,有:
$${{\mu }_{m}}(l+1)={{\mu }_{m}}(l)+c\left[ {{{\mathbf{\bar{b}}}}_{m}}\mathbf{p}(l+1)+{{d}_{m}}-{{s}_{m}}(l+1) \right]$$ (13) -
基于ADMM的分布式功率分配和接入控制联合算法如下:
1) 令l = 0;
2) 更新p
① k = 0;
② g = Bp(l,k) + d - s(l);
③ for m = 1: M
g = g - pm(l,k)bm;
$${{p}_{m}}(l,\ k+1)=\text{Proj}{{\,}_{[0,\ {{P}_{m}}]}}\left( \frac{-1-\mathbf{b}_{m}^{\text{T}}[\mathbf{\mu }(l)+c\mathbf{g}]}{c\mathbf{b}_{m}^{T}{{\mathbf{b}}_{m}}} \right)$$ g = g + pm(l,k + 1)bm;
end
④ 如果p收敛,p(l + 1) = p(l,k + 1);
否则,k = k + 1,转到②;
3) 更新s
for m = 1: M
$${{y}_{m}}={{\mu }_{m}}(l)+c{{\mathbf{\bar{b}}}_{m}}\mathbf{p}(l+1)+c{{d}_{m}}$$ $${{s}_{m}}(l+1)=\left\{ \begin{align} & 0\quad \ \left| {{y}_{m}} \right|\le {{\lambda }_{m}} \\ & \frac{(\left| {{y}_{m}} \right|-{{\lambda }_{m}}){{y}_{m}}}{c\left| {{y}_{m}} \right|}\quad \text{否则}\ \\ \end{align} \right.$$ end
4) 更新μ
for m = 1: M
$${{\mu }_{m}}(l+1)={{\mu }_{m}}(l)+c[{{\mathbf{\bar{b}}}_{m}}\mathbf{p}(l+1)+{{d}_{m}}$$ $$-{{s}_{m}}(l+1)]$$ end
5) 如果p、s和μ收敛,输出结果,停止;
否则,l = l + 1,转到步骤2)。
该算法的优点在于其主要步骤,更新p、s和μ等,都可以交给M个基站分布式进行,同时仅需要很少的数据即可完成。其中,更新p需要在各个基站间依次进行,而更新s和μ则可以完全并行地分布式进行。该算法的每一步都通过闭合式直接计算,因此复杂度低。在每轮迭代中,每个基站的运算量仅为O(M),是一种高效的分布式求解算法。
另一方面,与并行更新接入矢量s不同,在更新功率矢量p的每一轮迭代中,需要依次更新pm,m =1,2,…,M。尽管更新pm的运算非常简单,但当M较大时,每一轮迭代等待时间会较长。因此,本文的后续工作是进一步研究并行更新pm的算法,其主要难点是如何确保并行更新算法的收敛性。
Distributed Method for Joint Power Allocation and Admission Control Based on ADMM Framework
-
摘要: 传统网络节能问题在给定的服务质量(QoS)限制下实现传输功耗最小化。在此基础上引入接入控制,通过联合优化接入用户和发射功率进一步实现网络节能。当网络无法满足所有用户的QoS要求时,接入控制使网络服务尽可能多的用户。它还能对用户进行分类和挑选,以较低功耗代价满足接入用户的QoS条件,提高传输功效。在将原问题转换成一个近似的凸稀疏优化问题后,利用交错乘子法(ADMM)对其进行分布式迭代求解。该算法的每一步都具有闭合解,因此运算量很低。计算机仿真验证了该算法的正确性和有效性。Abstract: The traditional green networks are usually achieved by minimizing the transmit power under some quality of service (QoS) constraints. In this paper, admission control is integrated with power allocation to pursue further performance improvement in green communication. This technique enables the network to serve as many users as possible when the network cannot guarantee all the users' QoS requirements. Moreover, it classifies the users and picks out a subset of admissible users whose QoS constraints can be easily satisfied with relatively low power consumption; thus the power efficiency is improved. After approximating the original joint power and admission control problem by a convex sparse optimization problem, an efficient distributed algorithm is developed for the approximated problem by fitting it into the alternating direction method of multipliers (ADMM) framework. Specifically, each iteration of ADMM can be computed in closed form, thus giving it very low complexity. The effectiveness of the proposed algorithm is validated by a series of numerical simulations.
-
[1] 郭秉义. 绿色通信网络的节能算法研究[D]. 广州:华南理工大学, 2014. GUO Bing-yi. Energy saving research for green communications networks[D]. Guangzhou:South China University of Technology, 2014. [2] YATES R D. A framework for uplink power control in cellular radio systems[J]. IEEE Journal on Selected Areas in Communications, 1995, 13(7):1341-1347. doi: 10.1109/49.414651 [3] FOSCHINI G J, MILJANIC Z. A simple distributed autonomous power control algorithm and its convergence[J]. IEEE Trans Vehicular Technology, 1993, 42(4):641-646. doi: 10.1109/25.260747 [4] SCHUBERT M, BOCHE H. Solution of multiuser downlink beamforming problem with individual SINR constraints[J]. IEEE Trans Vehicular Technology, 2004, 53(1):18-28. doi: 10.1109/TVT.2003.819629 [5] LIN J R, LI Y B, PENG Q C. Joint power allocation, base station assignment and beamformer design for an uplink SIMO heterogeneous network[C]//IEEE Int'l Conf on Acoustics Speech, and Signal Processing (ICASSP). Florence:IEEE, 2014:434-438. [6] SIDIROPOULOS N D, DAVIDSON T N, LUO Z Q. Transmit beamforming for physical-layer multicasting[J]. IEEE Trans Signal Processing, 2006, 54(6):2239-2251. doi: 10.1109/TSP.2006.872578 [7] 黄晓燕, 毛玉明, 吴凡, 等. 中继增强的无线蜂窝多小区系统的联合调度与功率控制算法[J]. 电子与信息学报, 2012, 37(7):1665-1671. http://www.cnki.com.cn/Article/CJFDTOTAL-DZYX201207025.htm HUANG Xiao-yan, MAO Yu-ming, WU Fan, et al. Coordinated scheduling and power control algorithm for relay-enhanced cellular network with multicells[J]. Journal of Electronics & Information Technology, 2012, 37(7):1665-1671. http://www.cnki.com.cn/Article/CJFDTOTAL-DZYX201207025.htm [8] SARAYDA C U, MANDAYAM N B, GOODMAN D J. Pricing and power control in a multicell wireless date network[J]. IEEE Journal on Selected Areas in Communications, 2001, 19(10):1883-1892. doi: 10.1109/49.957304 [9] 冯文江, 蒋卫恒, 邓艺娜, 等. 基于非信任中继协作的保密通信联合功率控制[J]. 通信学报, 2014, 35(11):59-68. http://www.cnki.com.cn/Article/CJFDTOTAL-TXXB201411007.htm FENG Wen-jiang, JIANG Wei-heng, DENG Yi-na, et al. Joint power control for untrusted relay cooperation-based confidential communication[J]. Journal on Communications, 2014, 35(11):59-68. http://www.cnki.com.cn/Article/CJFDTOTAL-TXXB201411007.htm [10] GESBERT D, HANLY S, HUANG H, et al. Multi-cell MIMO cooperative networks:a new look at interference[J]. IEEE Journal on Selected Areas in Communications, 2010, 28(9):1380-1408. doi: 10.1109/JSAC.2010.101202 [11] DAHROUJ H, YU W. Coordinated beamforming for the multi-cell multi-antenna wireless system[J]. IEEE Trans Wireless Communications, 2010, 13(7):1332-1340. http://cn.bing.com/academic/profile?id=2151896803&encoded=0&v=paper_preview&mkt=zh-cn [12] 何蓉. 无线Mesh网络高效公平接入控制技术研究[D]. 成都:西南交通大学, 2011. http://cdmd.cnki.com.cn/article/cdmd-10613-1012389674.htm HE Rong. Research on efficient and fair access control technologies in wireless mesh network[D]. Chengdu:Southwest Jiaotong University, 2011. http://cdmd.cnki.com.cn/article/cdmd-10613-1012389674.htm [13] 苏伟, 刘琪, 袁坚, 等. 基于终端移动与业务达到认知的自适应无线接入控制机制[J]. 电子学报, 2011, 39(9):2148-2153. http://www.cnki.com.cn/Article/CJFDTOTAL-DZXU201109033.htm SU Wei, LIU Qi, YUAN Jian, et al. Adaptive radio access control scheme based on terminal mobility and service arriving cognition[J]. Acta Electronica Sinica, 2011, 39(9):2148-2153. http://www.cnki.com.cn/Article/CJFDTOTAL-DZXU201109033.htm [14] 王庆敏. 绿色节能通信网络中基站与用户的联合最优接入控制[D]. 大连:大连理工大学, 2013. http://cdmd.cnki.com.cn/article/cdmd-10141-1013198988.htm WANG Qing-min. Optimal joint base station and user equipment admission control for green wireless cellular networks[D]. Dalian:Dalian University of Technology, 2013. http://cdmd.cnki.com.cn/article/cdmd-10141-1013198988.htm [15] KUANG Q, SPEIDEL J, DROSTE H. Joint base-station association, channel assignment, beamforming and power control in heterogeneous networks[C]//IEEE the 75th Vehicular Technology Conference. Yokuhama:IEEE, 2012:1-5. [16] LIU Y F, DAI Y H, LUO Z Q. Joint power and admission control via linear programming deflation[J]. IEEE Trans Signal Processing, 2013, 61(6):1327-1338. doi: 10.1109/TSP.2012.2236319 [17] FUKUSHIMA M. Application of the alternating direction method of multipliers to separable convex programming problems[J]. Computational Optimization and Applications, 1992, 1(1):93-111. doi: 10.1007/BF00247655 [18] GRANT M C, BOYD S P. The CVX users' guide[[J/OL].[2014-10-24]. http://cvxr.com/cvx/. [19] BERTSEKAS D P. Nonlinear programming[M]. Second ed. Belmont, MA, USA:Athena Scientific, 1999. [20] KIWIEL K. Methods of descent for nondifferentiable optimization[M]. Berlin:Springer, 1985.