CN102711262B - 一种异构网络下行链路资源分配方法 - Google Patents

Abstract

本发明公开了宏小区和毫微微小区混合网络下的通信资源管理技术领域中的一种异构网络下行链路资源分配方法。包括组建用户集合和可用资源块集合，计算并组成用户优先级列表；计算用户需求的资源块数量；从用户优先级列表中选择优先级最高的用户并从可用资源块集合中选择信道质量最好的资源块分配给该优先级最高的用户，直至分配给该优先级最高的用户的资源块等于该优先级最高的用户需求的资源块数量或者可用资源块集合为空；计算并更新宏小区/毫微微小区的基站使用可用资源块与用户通信时在该可用资源块上的传输功率。本发明降低了小区间的干扰，保证了每个用户每种业务的QoS需求，实现了系统吞吐量的最大化、用户之间以及业务之间的公平。

Description

一种异构网络下行链路资源分配方法

技术领域

本发明属于宏小区和毫微微小区混合网络下的通信资源管理技术领域，尤其涉及一种异构网络下行链路资源分配方法。

背景技术

移动用户对传输带宽和数据速率的需求呈指数增长趋势，IMT-Advanced（International Mobile Telecommunications-Advanced，高级国际移动通信）传统蜂窝网络架构已不能满足室内、中小企业等热点区域的无缝覆盖，因此，IMT-Advanced异构网络架构应运而生。

图1是异构网络系统构架图，如图1所示，在IMT-Advanced系统的Macro-Femto（Macrocell-Femtocell,宏小区-毫微微小区）异构网络中，毫微微小区通过在宏小区中插入大量的家庭基站拉近天线与用户之间的距离，实现小区分裂，大大降低发射功率，延长终端电池使用时间，并获得高的信噪比，从而提高了系统吞吐量，降低了中断概率。

家庭基站的密集部署、重叠覆盖会造成复杂的干扰结构，并且由于家庭基站根据用户的需求安装甚至移动，因此运营商无法进行家庭基站的网络规划与优化，这更加剧了干扰管理的难度。图2是异构网络下行链路干扰情境示意图，如图2所示，IMT-Advanced系统Macro-Femto异构网络架构面临的一个主要挑战就是：下行链路同信道同层干扰和跨层干扰。因此，更加有效的干扰管理技术就成为家庭基站标准化过程中首要解决的问题之一。

在异构网络中，由于各子层之间（即应用层、传输层、网络层、媒体接入控制层和物理层之间）为维护自身的独立性与完整性而存在的大量信息冗余，不仅过多消耗了带宽与功率等紧缺资源，而且直接制约了系统性能充分发挥。于是，高效跨层的、跨模块的干扰协调方法逐步成为业界关注的热点。

现有的Macro-Femto异构网络小区间干扰协调技术主要包含以下几类：

①智能功率控制，只考虑了物理层技术；

②Macro与Femto实行正交资源分配，频谱利用率低；

③多点协作传输与接收，MeNodeB（Mecro Node Base，宏基站）与HeNodeB（Home Node Base，家庭基站）以及HeNodeB与HeNodeB之间需要交换信息，信令开销大，实现复杂度高，对标准的改动较大，不具有很好的后向兼容性。

现有的跨层动态资源分配技术的缺点包括：

①大部分技术只考虑MAC层与物理层的信息交互，少部分技术考虑了网络层、MAC层、物理层之间的联合优化，这些技术都没有将TCP/IP五层联合设计，没有最大化系统性能提升；

②考虑的业务类型不全，用户发起的业务类型单一，而在实际通信系统中，用户可以同时发起多种不同QoS需求的业务；

③公平性指标只注重用户间吞吐量的平衡，大部分技术没有考虑时延、时延抖动、分组丢失率等QoS指标在用户之间的平衡。

发明内容

本发明的目的在于，提出一种异构网络下行链路资源分配方法，将动态跨层资源分配技术应用于Macro-Femto异构网络，达到降低小区间干扰、最大化系统吞吐量、保证每个用户每种业务的QoS需求、实现用户之间的公平以及业务之间的公平。

为实现上述目的，本发明提供的技术方案是，一种异构网络下行链路资源分配方法，将动态跨层资源分配技术应用于IMT-Advanced系统的Macro-Femto异构网络中，其特征是所述方法包括：

步骤1：将异构网络中宏小区/毫微微小区的用户放入集合中组成用户集合，将所有可用资源块放入集合中组成可用资源块集合，计算用户集合中每个用户的优先级，将用户按照优先级由高到低的顺序排序并组成用户优先级列表；

步骤2：在用户集合中选择优先级最高的用户；

步骤3：计算该优先级最高的用户需求的资源块数量；

步骤4：将该优先级最高的用户从用户集合中删除；

步骤5：判断用户集合是否为空，如果用户集合为空，则执行步骤6；否则，返回步骤2；

步骤6：从用户优先级列表中选择优先级最高的用户；

步骤7：从可用资源块集合中选择信道质量最好的资源块分配给该优先级最高的用户；

步骤8：将已经分配给该优先级最高的用户的信道质量最好的资源块从可用资源块集合中删除；

步骤9：判断可用资源块集合是否为空，如果可用资源块集合为空，则执行步骤13；否则，执行步骤10；

步骤10：判断分配给该优先级最高的用户的资源块的数量是否小于该优先级最高的用户需求的资源块数量，如果分配给该优先级最高的用户的资源块的数量小于该优先级最高的用户需求的资源块数量，则返回步骤7；否则，执行步骤11；

步骤11：将该优先级最高的用户从用户优先级列表中删除；

步骤12：判断用户优先级列表是否为空，如果用户优先级列表为空，则执行步骤13；否则，返回步骤6；

步骤13：计算并更新宏小区/毫微微小区的基站使用可用资源块与用户通信时在该可用资源块上的传输功率。

所述计算用户集合中每个用户的优先级采用公式W_i表示用户i的优先级，表示用户i的实时业务的权重，表示用户i的非实时业务的权重，表示用户i的尽力而为业务的权重，i∈M_u或i∈F_u，M_u表示所有宏小区的用户集合，F_u表示所有毫微微小区的用户集合，并且有

$W_{\mathrm{RT}}^i=\left\{\begin{array}{cc}{\mathrm{\&beta;}}_{\mathrm{RT}}(1+\frac{1}{T_{\mathrm{RT}\_\max }^i-\frac{L_{\mathrm{RT}\_\mathrm{HOL}}^i}{r_{\mathrm{RT}}^i}+\mathrm{\&epsiv;}}),& T_{\mathrm{RT}\_\max }^i\&GreaterEqual;\frac{L_{\mathrm{RT}\_\mathrm{HOL}}^i}{r_{\mathrm{RT}}^i},\mathrm{\&epsiv;}={10}^{-6}\\ 0,& T_{\mathrm{RT}\_\max }^i<\frac{L_{\mathrm{RT}\_\mathrm{HOL}}^i}{r_{\mathrm{RT}}^i}\end{array}\right.,$

$W_{\mathrm{NRT}}^i=\left\{\begin{array}{cc}{\mathrm{\&beta;}}_{\mathrm{NRT}}(1+\frac{R_{\mathrm{NRT}\_\min }^i}{r_{\mathrm{NRT}}^i-R_{\mathrm{NRT}\_\min }^i+\mathrm{\&epsiv;}}),& r_{\mathrm{NRT}}^i\&GreaterEqual;R_{\mathrm{NRT}\_\min }^i,\mathrm{\&epsiv;}={10}^{-6}\\ 0,& r_{\mathrm{NRT}}^i<R_{\mathrm{NRT}\_\min }^i\end{array}\right.,$

$W_{\mathrm{BE}}^i=\left\{\begin{array}{cc}{\mathrm{\&beta;}}_{\mathrm{BE}}& N\_\mathrm{RB}\_\mathrm{remain}>0\\ 0& N\_\mathrm{RB}\_\mathrm{remain}=0\end{array}\right.,$

β_RT表示实时业务的优先级，表示用户i的实时业务的队列头分组的大小，ε为设定值，表示实时业务所能容忍的最大传输时延，表示用户i的实时业务当前时刻的数据传输速率，β_NRT表示非实时业务的优先级，表示非实时业务需求的最小数据传输速率，表示用户i的非实时业务当前时刻的数据传输速率，β_BE表示尽力而为业务的优先级，N_RB_remain表示满足服务宏小区或毫微微小区内所有用户实时业务和非实时业务的QoS需求之后剩余的资源块数量。

所述计算该优先级最高的用户需求的资源块数量采用公式

$N_{\mathrm{req}}^i=\frac{r_{\mathrm{req}}^i}{\underset{i={\mathrm{\&Sigma;}}_{m^{\&prime;}<m}\left|M_u^{m^{\&prime;}}\right|+1}{\overset{{\mathrm{\&Sigma;}}_{m^{\&prime;}<m}\left|M_u^{m^{\&prime;}}\right|+\left|M_u^m\right|}{\mathrm{\&Sigma;}}}{r}_{\mathrm{req}}^i}\&times;\left|J\right|,i\&Element;M_u^m,\&ForAll;m\&Element;M,$

其中，为宏小区用户i需求的资源块数量，M为所有宏小区组成的集合，为第m个宏小区的用户集合，为第m个宏小区的用户集合的势，为第m个宏小区的第1个用户在所有宏小区用户中的编号，为第m个宏小区的最后1个用户在所有宏小区用户中的编号，m′=1，2,...,m-1，为所有宏小区用户中的第i个用户请求的数据速率，为第m个宏小区的所有用户请求的数据速率之和，J为系统可用资源块集合，|J|为系统可用资源块集合的势。

所述计算该优先级最高的用户需求的资源块数量采用公式

$N_{\mathrm{req}}^{i^{\&prime;}}=\frac{r_{\mathrm{req}}^{i^{\&prime;}}}{\underset{i^{\&prime;}={\mathrm{\&Sigma;}}_{f^{\&prime;}<f^m}\left|F_u^{f^{\&prime;}}\right|+1}{\overset{{\mathrm{\&Sigma;}}_{f^{\&prime;}<f^m}\left|F_u^{f^{\&prime;}}\right|+\left|F_u^{f^m}\right|}{\mathrm{\&Sigma;}}}{r}_{\mathrm{req}}^{i^{\&prime;}}}\&times;\left|J\right|,i^{\&prime;}\&Element;F_u^{f^m},\&ForAll;m\&Element;M,$

其中，为毫微微小区用户i'需求的资源块数量，M为所有宏小区组成的集合，f^m为第m个宏小区中的毫微微小区且f^m∈F^m，F^m为第m个宏小区的毫微微小区集合，为第m个宏小区的第f^m个毫微微小区的用户集合，为集合的势，为第m个宏小区的第f^m个毫微微小区的第1个用户在所有毫微微小区用户中的编号，为第m个宏小区的第f^m个毫微微小区的最后1个用户在所有毫微微小区用户中的编号，f′=1,2,...,f^m-1，为所有毫微微小区用户中的第i'个用户请求的数据速率，为第m个宏小区的第f^m个毫微微小区的所有用户请求的数据速率之和，J为系统可用资源块集合，|J|为系统可用资源块集合的势。

所述计算分配给每个用户的每个资源块上的传输功率采用公式

$p_{m,i}^j={(\frac{B}{{\mathrm{J\&lambda;}}_i}-\frac{N_{0}B/J+I}{{\mathrm{\&Gamma;}}_i{\left|h_{m,i}^j\right|}^2})}^{+}={(u_i-\frac{N_{0}B/J+I}{{\mathrm{\&Gamma;}}_i{\left|h_{m,i}^j\right|}^2})}^{+},$

其中，为第m个宏小区基站使用第j个可用资源块与宏小区用户i通信时在所述第j个可用资源块上的传输功率，B为整个系统的可用带宽，N₀为高斯白噪声功率谱密度，J为可用资源块的数量，I为相邻宏小区和毫微微小区基站发射的下行链路信号对宏小区用户产生的同信道干扰，Γ_i为由误码率需求决定的系统参数，为第m个宏小区基站使用第j个可用资源块与宏小区用户i通信时的信道增益，λ_i为拉格朗日乘子并且有 为第m个宏小区基站使用第j个可用资源块与宏小区用户i通信的指示因子并且或 表示第m个宏小区基站不使用第j个可用资源块与宏小区用户i通信，表示第m个宏小区基站使用第j个可用资源块与宏小区用户i通信，为宏小区基站的最大允许传输功率，u_i为注水水平，函数(x)⁺等价于max(0,x)。

所述计算分配给每个用户的每个资源块上的传输功率采用公式

$p_{f,i^{\&prime;}}^j={(\frac{B}{{\mathrm{J\&lambda;}}_{i^{\&prime;}}}-\frac{N_{0}B/J+I^{\&prime;}}{{\mathrm{\&Gamma;}}_{i^{\&prime;}}{\left|h_{f,i^{\&prime;}}^j\right|}^2})}^{+}={(u_{i^{\&prime;}}-\frac{N_{0}B/J+I^{\&prime;}}{{\mathrm{\&Gamma;}}_{i^{\&prime;}}{\left|h_{f,i^{\&prime;}}^j\right|}^2})}^{+},$

其中，为第f个毫微微小区基站使用第j个可用资源块与毫微微小区用户i′通信时在所述第j个可用资源块上的传输功率，B为整个系统的可用带宽，N₀为高斯白噪声功率谱密度，J为可用资源块的数量，I′为相邻宏小区和毫微微小区基站发射的下行链路信号对毫微微小区用户产生的同信道干扰，Γ_i′为由误码率需求决定的系统参数，为第f个毫微微小区基站使用第j个可用资源块与毫微微小区用户i'通信时的信道增益；λ_i′为拉格朗日乘子并且有 是第f个毫微微小区基站使用第j个可用资源块与毫微微小区用户i'通信的指示因子并且或 表示第f个毫微微小区基站不使用第j个可用资源块与毫微微小区用户i′通信，表示第f个毫微微小区基站使用第j个可用资源块与毫微微小区用户i′通信，为毫微微小区基站的最大允许传输功率，u_i'为注水水平，函数(x)⁺等价于max(0,x)。

本发明降低了小区间的干扰，保证了每个用户每种业务的QoS需求，实现了系统吞吐量的最大化、用户之间的公平以及业务之间的公平。

附图说明

图1是异构网络系统构架图；

图2是异构网络下行链路干扰情境示意图；

图3是动态跨层资源分配示意图；

图4是异构网络下行链路资源分配方法流程图；

图5是主要参数传递示意图。

具体实施方式

下面结合附图，对优选实施例作详细说明。应该强调的是，下述说明仅仅是示例性的，而不是为了限制本发明的范围及其应用。

本发明针对IMT-Advanced系统Macro-Femto异构融合网络多小区、多业务、多用户下行链路情境，设计一种动态跨层资源分配方法，旨在实现三个保证：保证每个用户不同业务QoS需求；保证不同用户之间的公平性；保证同一用户不同业务之间的公平性，同时最大化整体网络吞吐量。图3是动态跨层资源分配示意图，如图3所示，跨层设计的基本思想在于：通过层与层之间信息的交互，实现整体系统性能优化。图3中，设计TCP/IP模型五层之间的信息交互以及采取的相应措施。①应用层：应用层向其他层传递其QoS需求，例如：时延范围、时延抖动、吞吐量及分组丢失率。应用层根据从底层获得的信道状态信息调整其工作方式，降低对带宽的需求。②传输层：传输层利用从网络层、媒体接入控制层以及物理层传递过来的信息降低拥塞误判概率，提高传输层吞吐量。例如：传输层利用媒体接入控制层的重传机制调整重传定时器的取值。传输层利用从应用层获得的QoS信息调整其滑动窗口的大小。传输层将数据分组丢失率及吞吐量信息提供给应用层，辅助应用层调整发送速率。③网络层：网络层根据应用层对QoS的不同需求，将数据分组路由到不同的物理网络接口上。网络层向传输层提供移动IP的切换信息来控制其定时器。网络层根据媒体接入控制层的无线信道信号强度减少移动IP切换带来的时延。④媒体接入控制层：媒体接入控制层根据应用层的QoS信息对数据帧进行不同的处理，赋予低时延需求的数据帧更高的优先级，赋予可靠性需求高的数据帧更强的纠错编码及更多的重传次数。利用传输层提供的往返时间(RTT)及重传定时器(RTO)来控制媒体接入控制层的重传机制。⑤物理层：物理层向网络层提供信道状态信息辅助网络层进行路由选择。物理层根据媒体接入控制层提供的控制指令调整相应的参数。物理层向应用层提供信道状态信息辅助应用层调整业务编码方式及传输速率。

首先，针对IMT-Advanced系统Macro-Femto异构网络，建立动态跨层资源分配模型。

在IMT-Advanced系统中，将宏小区集合用Μ={1,..,M}表示，第m（m∈Μ）个宏小区内的毫微微小区集合用F^m={1,..,F^m}表示。所有宏小区的所有毫微微小区的集合用来表示。第m个宏小区内的毫微微小区f^m(f^m∈F^m)的用户集合用来表示。其中，表示毫微微小区f^m的第l个用户，表示集合的势(一个集合的势等于这个集合的元素个数)。第m个宏小区内的所有毫微微小区用户用方法进行编号：第m个宏小区内毫微微小区1的用户编号集合用来表示，第m个宏小区内的毫微微小区2的用户编号集合用 $F_u^{f^m=2}=\left\{\right|F_u^{f^m=1}|+1,...,|F_u^{f^m=1}|+|F_u^{f^m=2}\left|\right\}$ 来表示，以此类推。表示第m个宏小区内所有毫微微小区用户编号构成的集合。所有宏小区的所有毫微微小区内的用户编号构成的集合用来表示。宏小区m(m∈M)的用户集合用来表示，其中，m_l表示宏小区m的第l个用户，表示集合的势。所有宏小区的用户用方法进行编号：宏小区1的用户编号集合用来表示，宏小区2的用户编号集合用来表示，以此类推。表示所有宏小区用户编号构成的集合。J={1,...,J}表示基站可用资源块构成的集合。表示宏小区第i个用户在第j个可用资源块上的SINR（Signal to interference andnoise ratio，信干噪比），相应的向量表示为 表示毫微微小区第i'个用户在第j个可用资源块上的SINR，相应的向量表示为 ${\mathrm{\&gamma;}}_{i^{\&prime;}}={[{\mathrm{\&gamma;}}_{i^{\&prime;}}^1,...,{\mathrm{\&gamma;}}_{i^{\&prime;}}^J]}^T,$ ${\&ForAll;i}^{\&prime;}\&Element;F_u.$

最优化问题建模为：

$\left(\max (\underset{i=1}{\overset{\left|M_u\right|}{\mathrm{\&Sigma;}}}{W}_i{U}_i+\underset{m=1}{\overset{M}{\mathrm{\&Sigma;}}}\underset{i^{\&prime;}=1}{\overset{\left|F_u^m\right|}{\mathrm{\&Sigma;}}}{W}_{i^{\&prime;}}{U}_{i^{\&prime;}})\right)---\left(1\right)$

上式的物理意义是最大化网络中所有用户加权效用之和。其中，

U_i=ln R_i (2)

U_i'＝ln R_i' (3)

$R_i=\underset{j=1}{\overset{J}{\mathrm{\&Sigma;}}}{b}_{m,i}^j\frac{B}{J}{\log }_2(1+{\mathrm{\&Gamma;}}_i{\mathrm{\&gamma;}}_i^j)=\underset{j=1}{\overset{J}{\mathrm{\&Sigma;}}}{b}_{m,i}^i\frac{B}{J}{\log }_2(1+{\mathrm{\&Gamma;}}_i\frac{p_{m,i}^j{\left|h_{m,i}^j\right|}^2}{N_{0}B/J+I})---\left(4\right)$

Γ_i=-1.5/ln(5BER) (5)

$I=\underset{k\&Element;M,k\&NotEqual;m}{\mathrm{\&Sigma;}}{b}_k^j{p}_k^j{\left|h_{k,i}^j\right|}^2+\underset{k^{\&prime;}\&Element;F}{\mathrm{\&Sigma;}}{b}_{k^{\&prime;}}^j{p}_{k^{\&prime;}}^j{\left|h_{k^{\&prime;},i}^j\right|}^2---\left(6\right)$

$R_{i^{\&prime;}}=\underset{j=1}{\overset{J}{\mathrm{\&Sigma;}}}{b}_{f,i^{\&prime;}}^j\frac{B}{J}{\log }_2(1+{\mathrm{\&Gamma;}}_{i^{\&prime;}}{\mathrm{\&gamma;}}_{i^{\&prime;}}^j)=\underset{j=1}{\overset{J}{\mathrm{\&Sigma;}}}{b}_{f,i^{\&prime;}}^i\frac{B}{J}{\log }_2(1+{\mathrm{\&Gamma;}}_{i^{\&prime;}}\frac{p_{f,i^{\&prime;}}^j{\left|h_{f,i^{\&prime;}}^j\right|}^2}{N_{0}B/J+I^{\&prime;}})---\left(7\right)$

$I^{\&prime;}=\underset{k\&Element;M}{\mathrm{\&Sigma;}}{b}_k^j{p}_k^j{\left|h_{k,i^{\&prime;}}^j\right|}^2+\underset{k^{\&prime;}\&Element;F{k}^{\&prime;}\&NotEqual;f}{\mathrm{\&Sigma;}}{b}_{k^{\&prime;}}^j{p}_{k^{\&prime;}}^j{\left|h_{k^{\&prime;},i^{\&prime;}}^j\right|}^2---\left(8\right)$

物理层约束条件为：

$p_m^j\&GreaterEqual;0,j\&Element;J,\&ForAll;m\&Element;M---\left(9\right)$

$p_f^j\&GreaterEqual;0,j\&Element;J,\&ForAll;f\&Element;F---\left(10\right)$

$\underset{j\&Element;J}{\mathrm{\&Sigma;}}{p}_m^j\&le;P_m^{\max },\&ForAll;m\&Element;M---\left(11\right)$

$\underset{j\&Element;J}{\mathrm{\&Sigma;}}{p}_f^j\&le;P_f^{\max },\&ForAll;f\&Element;F---\left(12\right)$

$b_m^j\&Element;\left\{\mathrm{0,1}\right\},j=\mathrm{1,2},...,J,m\&Element;M---\left(13\right)$

$b_{m,i}^j\&Element;\left\{\mathrm{0,1}\right\},j=\mathrm{1,2},...,J,m\&Element;M,i\&Element;M_u^m---\left(14\right)$

$\underset{i={\mathrm{\&Sigma;}}_{m^{\&prime;}<m}\left|M_u^{m^{\&prime;}}\right|+1}{\overset{{\mathrm{\&Sigma;}}_{m^{\&prime;}<m}\left|M_u^{m^{\&prime;}}\right|+\left|M_u^m\right|}{\mathrm{\&Sigma;}}}{b}_{m,i}^j\&le;1,j\&Element;J,m\&Element;M---\left(15\right)$

$b_f^j\&Element;\left\{\mathrm{0,1}\right\},j=\mathrm{1,2},...,J,\&ForAll;f\&Element;F---\left(16\right)$

$b_{f,i^{\&prime;}}^j\&Element;\left\{\mathrm{0,1}\right\},j\&Element;J,\&ForAll;f\&Element;F---\left(17\right)$

$\underset{i^{\&prime;}={\mathrm{\&Sigma;}}_{f^{\&prime;}<f^m}\left|F_u^{f^{\&prime;}}\right|+1}{\overset{{\mathrm{\&Sigma;}}_{f^{\&prime;}<f^m}\left|F_u^{f^{\&prime;}}\right|+\left|F_u^{f^m}\right|}{\mathrm{\&Sigma;}}}{b}_{f,i^{\&prime;}}^j\&le;1,j\&Element;J,\&ForAll;f^m\&Element;F^m---\left(18\right)$

MAC层约束条件为：

i∈M_u或i∈F_u (19)

i∈M_u或i∈F_u (20)

网络层约束条件：

i∈M_u或i∈F_u (21)

传输层约束条件：

Pr_pd≤PER (23)

应用层约束条件为：用户发起的各种业务的QoS需求。

上述公式中，W_i表示宏小区用户i的权重；U_i表示宏小区用户i的效用；W_i′表示毫微微小区用户i'的权重；U_i'表示毫微微小区用户i′的效用；R_i表示宏小区用户i的吞吐量；R_i'表示毫微微小区用户i′的吞吐量；是第m个宏小区基站使用第j个可用资源块与宏小区用户i通信的指示因子，表示第m个宏小区基站不使用第j个可用资源块与宏小区用户i通信，表示第m个宏小区基站使用第j个可用资源块与宏小区用户i通信；B表示整个系统的可用带宽；Γ_i和Γ_i'是由误码率需求决定的系统参数；是第m个宏小区基站使用第j个可用资源块与宏小区用户i通信时在第j个可用资源块上的传输功率；表示第m个宏小区基站使用第j个可用资源块与宏小区用户i通信时的信道增益；N₀表示高斯白噪声功率谱密度；I表示相邻宏小区和毫微微小区基站发射的下行链路信号对宏小区用户产生的同信道干扰；是第m个宏小区基站使用第j个可用资源块的指示因子，表示第m个宏小区基站不使用第j个可用资源块，表示第m个宏小区基站使用第j个可用资源块；表示第m个宏小区基站在第j个可用资源块上的传输功率；是第f个毫微微小区基站使用第j个可用资源块的指示因子，表示第f个毫微微小区基站不使用第j个可用资源块，表示第f个毫微微小区基站使用第j个可用资源块；表示第f个毫微微小区基站在第j个可用资源块上的传输功率；表示第f个毫微微小区基站使用第j个可用资源块与毫微微小区用户i′通信时的信道增益；是第f个毫微微小区基站使用第j个可用资源块与毫微微小区用户i′通信的指示因子，表示第f个毫微微小区基站不使用第j个可用资源块与毫微微小区用户i'通信，表示第f个毫微微小区基站使用第j个可用资源块与毫微微小区用户i'通信；是第f个毫微微小区基站使用第j个可用资源块与毫微微小区用户i'通信时在第j个可用资源块上的传输功率；I′表示相邻宏小区和毫微微小区基站发射的下行链路信号对毫微微小区用户产生的同信道干扰。

上述公式（9）指明第m个宏小区基站分配给第j个可用资源块的功率应该为一个正实数。

上述公式（10）指明第f个毫微微小区基站分配给第j个可用资源块的功率应该为一个正实数。

上述公式（11）指明第m个宏小区基站在所有资源块上的传输功率之和不应超过基站所能够传输的最大功率，表示宏小区基站的最大允许传输功率。

上述公式（12）指明第f个毫微微小区基站在所有资源块上的传输功率之和不应超过基站所能够传输的最大功率，表示毫微微小区基站的最大允许传输功率。

上述公式（13）指明第m个宏小区基站使用第j个可用资源块的指示因子只能取0或1。

上述公式（14）指明第m个宏小区基站使用第j个可用资源块与宏小区用户i通信的指示因子只能取0或1。

上述公式（15）指明第m个宏小区基站只能将第j个可用资源块至多分配给本小区内的一个宏小区用户。

上述公式（16）指明第f个毫微微小区基站使用第j个可用资源块的指示因子只能取0或1。

上述公式（17）指明第f个毫微微小区基站使用第j个可用资源块与毫微微小区用户i'通信的指示因子只能取0或1。

上述公式（18）指明第f^m个毫微微小区基站只能将第j个可用资源块至多分配给本小区内的一个毫微微小区用户。

上述公式（19）指明用户i实时业务的HOL(head of line,队列的头分组)的传输时延不应超过实时业务分组所能容忍的最大传输时延

上述公式（20）指明用户i非实时业务的数据传输速率不应小于非实时业务的数据传输速率最小需求

上述公式（21）指明用户i实时业务、非实时业务和尽力而为业务的数据传输速率之和不应超过网络分配给用户i的可用资源块所能够承载的吞吐量，和分别表示用户i实时业务、非实时业务和尽力而为业务当前时刻的数据传输速率。尽力而为业务是指没有服务质量要求的业务，这类业务的优先级低于实时业务和非实时业务，系统优先将资源分配给实时业务和非实时业务。如果有资源剩余，尽力而为业务才会得到服务，比如ftp(filetransmission protocol,文件传送协议)业务。

上述公式（22）指明传输层的重传时间T_重传不应超过TCP协议定时器所规定的时间

上述公式（23）指明分组丢失概率Pr_pd不应超过误包率要求PER。

实时业务（RT，real time）、非实时业务(NRT，non real time)和尽力而为业务(BE，best effort)的优先级分别为β_RT、β_NRT和β_BE(β_RT>β_NRT>β_BE)。

用户的权重如下：

$W_i=W_{\mathrm{RT}}^i+W_{\mathrm{NRT}}^i+W_{\mathrm{BE}}^i,$ i∈M_u或i∈F_u(24)

$W_{\mathrm{RT}}^i=\left\{\begin{array}{cc}{\mathrm{\&beta;}}_{\mathrm{RT}}(1+\frac{1}{T_{\mathrm{RT}\_\max }^i-\frac{L_{\mathrm{RT}\_\mathrm{HOL}}^i}{r_{\mathrm{RT}}^i}+\mathrm{\&epsiv;}}),& T_{\mathrm{RT}\_\max }^i\&GreaterEqual;\frac{L_{\mathrm{RT}\_\mathrm{HOL}}^i}{r_{\mathrm{RT}}^i},\mathrm{\&epsiv;}={10}^{-6}\\ 0,& T_{\mathrm{RT}\_\max }^i<\frac{L_{\mathrm{RT}\_\mathrm{HOL}}^i}{r_{\mathrm{RT}}^i}\end{array}\right.,$ i∈M_u或i∈F_u(25)

$W_{\mathrm{NRT}}^i=\left\{\begin{array}{cc}{\mathrm{\&beta;}}_{\mathrm{NRT}}(1+\frac{R_{\mathrm{NRT}\_\min }^i}{r_{\mathrm{NRT}}^i-R_{\mathrm{NRT}\_\min }^i+\mathrm{\&epsiv;}}),& r_{\mathrm{NRT}}^i\&GreaterEqual;R_{\mathrm{NRT}\_\min }^i,\mathrm{\&epsiv;}={10}^{-6}\\ 0,& r_{\mathrm{NRT}}^i<R_{\mathrm{NRT}\_\min }^i\end{array}\right.,$ i∈M_u或i∈F_u(26)

$W_{\mathrm{BE}}^i=\left\{\begin{array}{cc}{\mathrm{\&beta;}}_{\mathrm{BE}}& N\_\mathrm{RB}\_\mathrm{remain}>0\\ 0& N\_\mathrm{RB}\_\mathrm{remain}=0\end{array}\right.,$ i∈M_u或i∈F_u(27)

其中，W_i表示用户i的权重，等于用户i的所有业务的权重之和；表示用户i的实时业务的权重；表示用户i的非实时业务的权重；表示用户i的尽力而为业务的权重；表示用户i实时业务HOL(head of line,队列的头分组)的大小；ε是一个极小正数，用来防止分母为0的情况出现；N_RB_remain表示满足服务小区内所有用户实时业务和非实时业务的QoS需求之后剩余的资源块数目。

传统的公平性观点：保证系统中的每个用户的可获速率不低于门限值。IMT-Advanced系统业务种类多样化，每种业务的QoS要求不同。因此，传统的公平性观点不再适用。新型公平性观点：保证每个用户每种业务的QoS满足用户的需求。传统的资源分配算法可以实现系统效率与用户速率公平性之间的折中。但是这些算法在设计时忽略了一个事实，那就是对于不同的应用，即使获得的数据速率是相同的，但是用户的效用值(也就是用户的满意度)是不同的。此时，为不同的用户分配相同的速率并不能保证用户之间的真正公平。例如，如果所有用户都为“尽力而为”(Best Effort,BE)业务，保证速率公平就是保证了用户之间的公平性，如果其中一些用户为视频业务，当网络负载较重时，为视频业务用户和BE业务用户分配相同的数据速率，对视频业务用户是非常不公平的。由此可得，当用户采用不同应用时，效用公平才是系统优化时应该考虑的优化目标。在传统宏蜂窝网络中，使用Jain’s公平指数衡量算法公平性，其数学定义为：

$x=\frac{{\left(\underset{i=1}{\overset{K}{\mathrm{\&Sigma;}}}{r}_i\right)}^2}{K\underset{i=1}{\overset{K}{\mathrm{\&Sigma;}}}{r}_i^2}$

其中，r_i为用户i分配到的资源，如速率。对于一个包含K个类型相同的用户的网络，Jain’s公平指数的取值范围为[1/K,1],其值越大说明越公平。Jain’s公平指数不再适应于多业务的异构融合网络，应当对Jain’s公平指数进行修正。修正后的Jain’s公平指数如下：

修正后的Jain’s公平指数则能够更精确地反映未来多业务网络环境中各种资源分配算法的公平性。

基于上述动态跨层资源分配模型，本发明提供的一种异构网络下行链路资源分配方法如图4所示，包括：

步骤1：将异构网络中宏小区或者毫微微小区的用户放入集合中组成用户集合，将所有可用资源块放入集合中组成可用资源块集合，计算用户集合中每个用户的优先级，将用户按照优先级由高到低的顺序排序并组成用户优先级列表。

用户集合中每个用户的优先级通过用户权重体现，用户权重越高，其优先级越高。因此，计算用户集合中每个用户的优先级可以采用上述公式（24）-（27），通过计算用户权重获得用户的优先级。其中，β_RT,β_NRT,β_BE由基站应用层传递给基站MAC层，由移动台应用层传递给基站MAC层，由基站传输层传递给基站MAC层，由基站网络层传递给基站MAC层，N_RB_remain由基站物理层传递给基站MAC层。

步骤2：在用户集合中选择优先级最高的用户。

步骤3：计算该优先级最高的用户需求的资源块数量。

对于宏小区用户，优先级最高的用户需求的资源块数量采用公式

$N_{\mathrm{req}}^i=\frac{r_{\mathrm{req}}^i}{\underset{i={\mathrm{\&Sigma;}}_{m^{\&prime;}<m}\left|M_u^{m^{\&prime;}}\right|+1}{\overset{{\mathrm{\&Sigma;}}_{m^{\&prime;}<m}\left|M_u^{m^{\&prime;}}\right|+\left|M_u^m\right|}{\mathrm{\&Sigma;}}}{r}_{\mathrm{req}}^i}\&times;\left|J\right|,i\&Element;M_u^m,\&ForAll;m\&Element;M---\left(28\right)$

进行计算。其中，为宏小区用户i需求的资源块数量，M为所有宏小区组成的集合，为第m个宏小区的用户集合，为第m个宏小区的用户集合的势，为第m个宏小区的第1个用户在所有宏小区用户中的编号，为第m个宏小区的最后1个用户在所有宏小区用户中的编号，m′=1,2,...,m-1，为宏小区用户i请求的数据速率，为第m个宏小区的所有用户请求的数据速率之和，J为系统可用资源块集合，|J|为系统可用资源块集合的势。

对于毫微微小区用户，优先级最高的用户需求的资源块数量采用公式

$N_{\mathrm{req}}^{i^{\&prime;}}=\frac{r_{\mathrm{req}}^{i^{\&prime;}}}{\underset{i^{\&prime;}={\mathrm{\&Sigma;}}_{f^{\&prime;}<f^m}\left|F_u^{f^{\&prime;}}\right|+1}{\overset{{\mathrm{\&Sigma;}}_{f^{\&prime;}<f^m}\left|F_u^{f^{\&prime;}}\right|+\left|F_u^{f^m}\right|}{\mathrm{\&Sigma;}}}{r}_{\mathrm{req}}^{i^{\&prime;}}}\&times;\left|J\right|,i^{\&prime;}\&Element;F_u^{f^m},\&ForAll;m\&Element;M---\left(29\right)$

进行计算。其中，为毫微微小区用户i'需求的资源块数量，M为所有宏小区组成的集合，f^m为第m个宏小区中的毫微微小区且f^m∈F^m，F^m为第m个宏小区的毫微微小区集合，为第m个宏小区的第f^m个毫微微小区的用户集合，为集合的势，为第m个宏小区的第f^m个毫微微小区的第1个用户在所有毫微微小区用户中的编号，为第m个宏小区的第f^m个毫微微小区的最后1个用户在所有毫微微小区用户中的编号，f′=1,2,...,f^m-1，为毫微微小区用户i′请求的数据速率，为第m个宏小区的第f^m个毫微微小区的所有用户请求的数据速率之和，J为系统可用资源块集合，|J|为系统可用资源块集合的势。

计算得到的或由移动台应用层传递给移动台物理层，之后，通过空中接口反馈给基站，基站物理层成功译码此信息后，将参数传递给基站MAC层。

步骤4：将该优先级最高的用户从用户集合中删除。

步骤5：判断用户集合是否为空，如果用户集合为空，则执行步骤6；否则，返回步骤2。

步骤6：从用户优先级列表中选择优先级最高的用户。

步骤7：从可用资源块集合中选择信道质量最好的资源块分配给该优先级最高的用户。

在资源块分配算法中，假设功率分配算法为等功率分配，即 信道质量最好的资源块由基站测量用户发送的导频信息，从而估计出信道质量最好的资源块。

步骤8：将已经分配给该优先级最高的用户的信道质量最好的资源块从可用资源块集合中删除。

步骤9：判断可用资源块集合是否为空，如果可用资源块集合为空，则执行步骤13；否则，执行步骤10。

步骤10：判断分配给该优先级最高的用户的资源块的数量是否小于该优先级最高的用户需求的资源块数量，如果分配给该优先级最高的用户的资源块的数量小于该优先级最高的用户需求的资源块数量，则返回步骤7；否则，执行步骤11。

步骤11：将该优先级最高的用户从用户优先级列表中删除。

步骤12：判断用户优先级列表是否为空，如果用户优先级列表为空，则执行步骤13；否则，返回步骤6。

步骤13：计算并更新宏小区或毫微微小区基站使用可用资源块与小区用户通信时在该可用资源块上的传输功率。

本发明的目的是在基站下行链路传输功率限制条件下，获得最大的吞吐量。当用户的资源块分配确定后，此最优化问题用公式(30)-(31)或(32)-(33)表示,公式(30)-(31)针对宏小区，公式(32)-(33)针对毫微微小区。

$\underset{p_{m,i}}{\max }{R}_i---\left(30\right)$

其中， $p_{m,i}=(b_{m,i}^1{p}_{m,i}^1,...,b_{m,i}^j{p}_{m,i}^j,...,b_{m,i}^J{p}_{m,i}^J),j\&Element;J.$

约束条件：

$\underset{j=1}{\overset{J}{\mathrm{\&Sigma;}}}{b}_{m,i}^j{p}_{m,i}^j\&le;\frac{\underset{j=1}{\overset{J}{\mathrm{\&Sigma;}}}{b}_{m,i}^j}{\left|J\right|}{P}_m^{\max },m\&Element;M,i\&Element;M_u---\left(31\right)$

$\underset{p_{f,i^{\&prime;}}}{\max }{R}_{i^{\&prime;}}---\left(32\right)$

其中， $p_{f,i^{\&prime;}}=(b_{f,i^{\&prime;}}^1{p}_{f,i^{\&prime;}}^1,...,b_{f,i^{\&prime;}}^j{p}_{f,i^{\&prime;}}^j,...,b_{f,i^{\&prime;}}^J{p}_{f,i^{\&prime;}}^J),j\&Element;J.$

约束条件：

$\underset{j=1}{\overset{J}{\mathrm{\&Sigma;}}}{b}_{f,i^{\&prime;}}^j{p}_{f,i^{\&prime;}}^j\&le;\frac{\underset{j=1}{\overset{J}{\mathrm{\&Sigma;}}}{b}_{f,i^{\&prime;}}^j}{\left|J\right|}{P}_f^{\max },i^{\&prime;}\&Element;F_u,f\&Element;F---\left(33\right)$

上述公式(30)中的目标函数是凹函数，此最优化问题可以通过拉格朗日方法求解。通过求解KKT条件，可得和λ_i如下所示：

$p_{m,i}^j={(\frac{B}{J{\mathrm{\&lambda;}}_i}-\frac{N_{0}B/J+I}{{\mathrm{\&Gamma;}}_i{\left|h_{m,i}^j\right|}^2})}^{+}={(u_i-\frac{N_{0}B/J+I}{{\mathrm{\&Gamma;}}_i{\left|h_{m,i}^j\right|}^2})}^{+}---\left(34\right)$

${\mathrm{\&lambda;}}_i=\frac{B\underset{j=1}{\overset{J}{\mathrm{\&Sigma;}}}{b}_{m,i}^j}{\underset{j=1}{\overset{J}{\mathrm{\&Sigma;}}}{b}_{m,i}^j{P}_m^{\max }+J\&times;\underset{j=1}{\overset{J}{\mathrm{\&Sigma;}}}\frac{N_{0}B/J+I}{{\mathrm{\&Gamma;}}_i{\left|h_{m,i}^j\right|}^2}}---\left(35\right)$

其中，u_i表示注水水平，

同理可得和λ_i′如下所示：

$p_{f,i^{\&prime;}}^j={(\frac{B}{J{\mathrm{\&lambda;}}_{i^{\&prime;}}}-\frac{N_{0}B/J+I^{\&prime;}}{{\mathrm{\&Gamma;}}_{i^{\&prime;}}{\left|h_{f,i^{\&prime;}}^j\right|}^2})}^{+}={(u_{i^{\&prime;}}-\frac{N_{0}B/J+I^{\&prime;}}{{\mathrm{\&Gamma;}}_{i^{\&prime;}}{\left|h_{f,i^{\&prime;}}^j\right|}^2})}^{+}---\left(36\right)$

${\mathrm{\&lambda;}}_{i^{\&prime;}}=\frac{B\underset{j=1}{\overset{J}{\mathrm{\&Sigma;}}}{b}_{f,i^{\&prime;}}^j}{\underset{j=1}{\overset{J}{\mathrm{\&Sigma;}}}{b}_{f,i^{\&prime;}}^j{P}_f^{\max }+J\&times;\underset{j=1}{\overset{J}{\mathrm{\&Sigma;}}}\frac{N_{0}B/J+I^{\&prime;}}{{\mathrm{\&Gamma;}}_{i^{\&prime;}}{\left|h_{f,i^{\&prime;}}^j\right|}^2}}---\left(37\right)$

其中，u_i′0为注水水平，

公式（30）-（37）中各个变量的含义在建立模型时已经说明，这里不再赘述。

对于宏小区，通过计算可以得到第m个宏小区基站使用第j个可用资源块与宏小区用户i通信时在第j个可用资源块上的传输功率，通过宏小区基站更新值，即可使系统获得最大的吞吐量。

对于毫微微小区，通过计算可以得到第f个毫微微小区基站使用第j个可用资源块与毫微微小区用户i'通信时在第j个可用资源块上的传输功率，通过毫微微小区基站（即家庭基站）更新值，即可使系统获得最大的吞吐量。

动态跨层资源分配模型中，各主要参数的传递过程如图5所示。

本发明的技术效果在于：

（1）联合优化TCP/IP五层参数与相应机制。现有的最好技术考虑了网络层、MAC层、物理层之间的联合优化，没有将TCP/IP五层联合设计，没有最大化系统性能，并且网络层、MAC层的约束并不全面。本发明的最优化数学模型的目的是最大化整个网络加权吞吐量，同时保证用户间的公平和业务间的公平。本发明的MAC层约束条件不仅考虑了实时业务的时延约束，而且考虑了非实时业务最小传输速率需求。本发明的传输层约束条件指明数据包的的重传时间不应超过定时器规定的时间，并且分组丢失概率不能大于业务的误包率要求。

（2）所设计的实现算法，分三步执行，易于工程实现，解决了最优化模型在多项式时间内不可解的问题。

（3）本发明建立的网络环境为多个macro小区，多个femto小区，每个用户可以同时执行多种类型的业务，物理层传输技术采用MIMO-OFDM。这种网络环境更符合实际网络情况，因此，所建立的模型和算法能更好地解决实际问题。

（4）在所建立的模型中，用户可以同时发起QoS不同的业务。而在现有最好技术中，每个用户只能发起一种类型的业务。因此，本发明采用的用户流量模型更符合真实通信情景。

（5）公平性指标更加全面，综合考虑了各种不同业务的QoS要求，包括了吞吐量公平、延时公平、丢包率公平等。而现有最好技术只局限于吞吐量公平和延时公平。修正后的Jain’s公平指数与传统的Jain’s公平指数相比，更适合于评估多媒体通信的公平性。在本修正公式提出之前，尚未发现本领域技术人员提出新的公平性评价指标。

（6）本发明更好地阐述了新型公平性观点：以业务为粒度，保障每种业务的QoS要求，实现业务之间的公平。改进了传统的以用户速率为粒度的公平性观点，更符合未来多业务的网络通信环境。现有最好的技术虽然意识到用户执行不同QoS需求的业务，但是对公平性阐述模糊，仍然采用传统的公平性评价指标。

以上所述，仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应该以权利要求的保护范围为准。

Claims (5)

1.一种异构网络下行链路资源分配方法，将动态跨层资源分配技术应用于IMT-Advanced系统的Macro-Femto异构网络中，其特征是所述方法包括：

步骤1：将异构网络中宏小区或毫微微小区的用户放入集合中组成用户集合，将所有可用资源块放入集合中组成可用资源块集合，计算用户集合中每个用户的优先级，将用户按照优先级由高到低的顺序排序并组成用户优先级列表；

其中，计算用户集合中每个用户的优先级采用公式，W_i表示用户i的优先级，表示用户i的实时业务的权重，表示用户i的非实时业务的权重，表示用户i的尽力而为业务的权重，i∈M_u或i∈F_u，M_u表示所有宏小区的用户集合，F_u表示所有毫微微小区的用户集合，并且有

β_RT表示实时业务的优先级，表示用户i的实时业务的队列头分组的大小，ε为设定值，表示实时业务所能容忍的最大传输时延，表示用户i的实时业务当前时刻的数据传输速率，β_NRT表示非实时业务的优先级，表示 非实时业务需求的最小数据传输速率，表示用户i的非实时业务当前时刻的数据传输速率，β_BE表示尽力而为业务的优先级，N_RB_remain表示满足服务宏小区或毫微微小区内所有用户实时业务和非实时业务的QoS需求之后剩余的资源块数量；

步骤2：在用户集合中选择优先级最高的用户；

步骤3：计算该优先级最高的用户需求的资源块数量；

步骤4：将该优先级最高的用户从用户集合中删除；

步骤5：判断用户集合是否为空，如果用户集合为空，则执行步骤6；否则，返回步骤2；

步骤6：从用户优先级列表中选择优先级最高的用户；

步骤7：从可用资源块集合中选择信道质量最好的资源块分配给该优先级最高的用户；

步骤8：将已经分配给该优先级最高的用户的信道质量最好的资源块从可用资源块集合中删除；

步骤9：判断可用资源块集合是否为空，如果可用资源块集合为空，则执行步骤13；否则，执行步骤10；

步骤10：判断分配给该优先级最高的用户的资源块的数量是否小于该优先级最高的用户需求的资源块数量，如果分配给该优先级最高的用户的资源块的数量小于该优先级最高的用户需求的资源块数量，则返回步骤7；否则，执行步骤11；

步骤11：将该优先级最高的用户从用户优先级列表中删除；

步骤12：判断用户优先级列表是否为空，如果用户优先级列表为空，则执行步骤13；否则，返回步骤6；

步骤13：计算并更新宏小区或毫微微小区的基站使用可用资源块与用户通信时在该可用资源块上的传输功率。

2.根据权利要求1所述的方法，其特征是所述计算该优先级最高的用户需求的资源块数量采用公式

其中，为宏小区用户i需求的资源块数量，M为所有宏小区组成的集合，为第m个宏小区的用户集合，为第m个宏小区的用户集合的势，为第m个宏小区的第1个用户在所有宏小区用户中的编号，为第m个宏小区的最后1个用户在所有宏小区用户中的编号，m'＝1,2,...,m-1，为所有宏小区用户中的第i个用户请求的数据速率，为第m个宏小区的所有用户请求的数据速率之和，J为系统可用资源块集合，|J|为系统可用资源块集合的势。

3.根据权利要求2所述的方法，其特征是所述计算宏小区的基站使用可用资源块与用户通信时在该可用资源块上的传输功率采用公式

其中，为第m个宏小区基站使用第j个可用资源块与宏小区用户i通信时在所述第j个可用资源块上的传输功率，B为整个系统的可用带宽，N₀为高斯白噪声功率谱密度，J为可用资源块的数量，I为相邻宏小区和毫微微小区基站发射的下行链路信号对宏小区用户产生的同信道干扰，Γ_i为由误码率需求决定的系统参数，为第m个宏小区基站使用第j个可用资源块与宏小区用户i通信时的信道 增益，λi为拉格朗日乘子并且有，为第m个宏小区基站使用第j个可用资源块与宏小区用户i通信的指示因子并且或 表示第m个宏小区基站不使用第j个可用资源块与宏小区用户i通信，表示第m个宏小区基站使用第j个可用资源块与宏小区用户i通信，为宏小区基站的最大允许传输功率，u_i为注水水平，函数(x)⁺等价于max(0,x)。

4.根据权利要求1所述的方法，其特征是所述计算该优先级最高的用户需求的资源块数量采用公式

其中，为毫微微小区用户i'需求的资源块数量，M为所有宏小区组成的集合，f^m为第m个宏小区中的毫微微小区且f^m∈F^m，F^m为第m个宏小区的毫微微小区集合，为第m个宏小区的第f^m个毫微微小区的用户集合，为集合的势， 为第m个宏小区的第f^m个毫微微小区的第1个用户在所有毫微微小区用户中的编号，为第m个宏小区的第f^m个毫微微小区的最后1个用户在所有毫微微小区用户中的编号，f'＝1,2,...,f^m-1，为所有毫微微小区用户中的第i'个用户请求的数据速率，为第m个宏小区的第f^m个毫微微小区的所有用户请求的数据速率之和，J为系统可用资源块集合，|J|为系统可用资源块集合的势。

5.根据权利要求4所述的方法，其特征是所述计算毫微微小区的基站使用可用资源块与用户通信时在该可用资源块上的传输功率采用公式

其中，为第f个毫微微小区基站使用第j个可用资源块与毫微微小区用户i'通信时在所述第j个可用资源块上的传输功率，B为整个系统的可用带宽，N₀为高斯白噪声功率谱密度，J为可用资源块的数量，I'为相邻宏小区和毫微微小区基站发射的下行链路信号对毫微微小区用户产生的同信道干扰，Γ_i'为由误码率需求决定的系统参数，为第f个毫微微小区基站使用第j个可用资源块与毫微微小区用户i'通信时的信道增益；λ_i'为拉格朗日乘子并且有 ，是第f个毫微微小区基站使用第j个可用资源块与毫微微小区用户i'通信的指示因子并且或表示第f个毫微微小区基站不使用第j个可用资源块与毫微微小区用户i'通信，表示第f个毫微微小区基站使用第j个可用资源块与毫微微小区用户i'通信，为毫微微小区基站的最大允许传输功率，u_i'为注水水平，函数(x)⁺等价于max(0,x)。