Next Article in Journal
Approach for the Development of a Framework for the Identification of Activities of Daily Living Using Sensors in Mobile Devices
Next Article in Special Issue
A Custom Approach for a Flexible, Real-Time and Reliable Software Defined Utility
Previous Article in Journal
Sensor Compromise Detection in Multiple-Target Tracking Systems
Previous Article in Special Issue
Application of a Multimedia Service and Resource Management Architecture for Fault Diagnosis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Robust Rate Maximization for Heterogeneous Wireless Networks under Channel Uncertainties

1
School of Communication and Information Engineering, Chongqing University of Posts and Telecommunications, Chongqing 400065, China
2
State Key Laboratory of Integrated Services Networks, Xidian University, Xi’an 710071, China
*
Author to whom correspondence should be addressed.
Sensors 2018, 18(2), 639; https://doi.org/10.3390/s18020639
Submission received: 20 January 2018 / Revised: 9 February 2018 / Accepted: 17 February 2018 / Published: 21 February 2018
(This article belongs to the Special Issue Advances on Resources Management for Multi-Platform Infrastructures)

Abstract

:
Heterogeneous wireless networks are a promising technology in next generation wireless communication networks, which has been shown to efficiently reduce the blind area of mobile communication and improve network coverage compared with the traditional wireless communication networks. In this paper, a robust power allocation problem for a two-tier heterogeneous wireless networks is formulated based on orthogonal frequency-division multiplexing technology. Under the consideration of imperfect channel state information (CSI), the robust sum-rate maximization problem is built while avoiding sever cross-tier interference to macrocell user and maintaining the minimum rate requirement of each femtocell user. To be practical, both of channel estimation errors from the femtocells to the macrocell and link uncertainties of each femtocell user are simultaneously considered in terms of outage probabilities of users. The optimization problem is analyzed under no CSI feedback with some cumulative distribution function and partial CSI with Gaussian distribution of channel estimation error. The robust optimization problem is converted into the convex optimization problem which is solved by using Lagrange dual theory and subgradient algorithm. Simulation results demonstrate the effectiveness of the proposed algorithm by the impact of channel uncertainties on the system performance.

Graphical Abstract

1. Introduction

With the exponential growth of mobile data driven by various communication applications, such as WIFI and smartphones, the traditional wireless network via macrocell base stations (BSs) can not satisfy higher communication requirements (e.g., throughput and coverage). As a new candidate technology in the fifth generation (5G) wireless communication networks, heterogeneous network is proposed to improve network coverage and data rate [1,2,3]. In heterogeneous networks, many femtocells with low-cost and low-energy consumption are distributed around macro BSs randomly meanwhile femtocell users (FUs) share the same spectrum resource with macrocell users (MUs) to obtain better spectral efficiency (e.g., co-channel development) and supplement traditional single-tier cellular systems. However, the cross-tier interference from femtocell users to macrocell receivers should be controlled severely. Therefore, interference mitigation is very important for which power control-based resource allocation (i.e., power allocation) has been a feasible method for Heterogeneous wireless networks (HetNets) [4,5].
The design of resource allocation algorithms in HetNets has attracted much research interest owing to its importance. The main task in resource allocation of underlay femtocell networks is to reduce the interference power received at MUs and simultaneously obtain the expected performance of femtocells, such as maximize the throughput of femtocells. To achieve this objective, various resource allocation (i.e., power control) algorithms have been extensively studied for various scenarios in [6,7,8]. In [6], a resource allocation problem in both the uplink and the downlink for two-tier femtocell networks is considered for maximizing the capacity of sensitive FUs and delay-tolerant FUs under the cross-tier interference constraint of MU and the quality-of-service (QoS) constraint of delay-sensitive users. In [7], the authors propose an interference mitigation strategy to improve the uplink throughput by setting a fixed interference threshold and adjusting the maximum transmit power of femtocell user. To enhance the energy efficiency (EE) of HetNet, in [8], two EE resource allocation algorithms via game theory are proposed for downlink transmission in multichannel macro-femto networks. Note that all the literatures above are designed based on the assumption of perfect knowledge of channel state information (CSI) at the transmitters. Although they have generally been assumed that complete system information (i.e., perfect channel state information) are available to FUs, nevertheless, due to the random nature of wireless channels and inaccurate channel estimation, as well as channel delays, it is impossible for FUs to obtain exact values of system parameters, such as channel gains and interference power from other networks.
In fact, due to lack of cooperation between macrocell networks and femtocell networks, obtaining complete system information pertaining to MUs is difficult for FUs. For instance, uncertain channel gains between FUs and MUs’ receivers may cause the sum aggregated interference of FUs on the MUs’ receivers to exceed the tolerable threshold, which means that the QoS requirements of users in the macrocell can not be guaranteed. As a result, from the perspective of MUs, this would increase the outage probability of MUs. Additionally, uncertainties on channel gains between femtocell BS and FUs’ receivers may also reduce the minimum rate requirement or actual signal-to-interference-plus-noise-ratio (SINR) of each FU at its femtocell BS or receiver below the target threshold. Therefore, the robustness of resource allocation algorithms should be considered ahead of time.
Thanks to robust optimization theory [9,10], robust resource allocation algorithms under imperfect CSI have drawn considerable attentions to overcome the effect of uncertainties in HetNets. Generally, the robust resource allocation designs are developed by two types of CSI errors, i.e., worst-case model and stochastic model [11]. The deterministic model assumes that the instantaneous value of uncertain parameter is bounded by the worst-case upper boundary (e.g., ellipsoidal uncertainty sets), which aims to obtain the worst-case robustness and guarantee. On the contrary, the stochastic approach assumes that the statistical information of CSI or channel estimation errors is known at the transmitters and seeks to a suboptimal solution for maintaining a certain outage probability on average, such as probabilistic SINR constraint. In [12], with the consideration of ellipsoidal uncertainty sets, Vaezpour et al. propose a robust distributed resource allocation algorithm in orthogonal frequency division multiple access (OFDMA) based femtocell networks that maximizes the total rate of FUs under the constraints on the minimum data requirement of each FU and the co-tier as well as cross-tier interference. In [13], based on a robust Stackelberg game, a downlink power control scheme under column-wise uncertainty model is proposed for maximizing the capacities of each FU and each MU, where the impact of uncertainty and Nash equilibrium point are analyzed by using variational inequality (VI) and Stackelberg game theory respectively. In [14], taking the ellipsoidal channel uncertainties into account, Xiao et al. present a robust resource allocation algorithm in full-duplex based OFDMA femtocell HetNets considering the throughput maximization of the femtocell while avoiding harmful interference to the macrocell. In [15], based on game theory, a robust power control scheme for heterogeneous users is proposed by taking the bounded channel uncertainties into account. In [16], Xu et al. propose a robust resource allocation algorithm for multiuser HetNets with macrocells and microcells through geometrical programming method. However, in practice, the case of worst-case estimation error may not always exist so that the proposed algorithms sacrifice system performance in a large part.
Moreover, the upper bound of uncertainty can not be easily obtained in practical Heterogeneous systems. Nevertheless, the robust resource allocation problem with Bayesian approach (i.e., probability constraint [17]) is more appropriate for stochastic scenario because of the stochastic nature of measurements and estimation errors in wireless networks. In [18,19], based on the hierarchical game theory, Liu et al. propose a robust uplink power allocation algorithm to minimize the transmit power of FUs for two-tier femtocell networks sharing the same frequency with macrocells, where the uncertainties of channel gains are formulated as the probabilistic constraints. But these algorithms aim to minimize the power consumption, and thus cannot be used for the rate maximization in HetNets under channel uncertainties. In [20], Mokari et al. investigate the robust quantized ergodic resource allocation scheme to achieve the sum-rate maximization of femtocell network while guaranteeing the macro network interference requirements with any desired high probability. However, only channel gains of interference links subject to probability constraint are considered. They do not consider the uncertainty of signal links. In [21], for a two-tier HetNet, a robust MISO transmit power optimization problem is considered to design a minimum power downlink beamformer under the outage-based QoS constraints. However, all the aforementioned work only focuses on the protections of MUs in outage-based form (e.g., probabilistic cross-tier interference constraint of each MU). Under channel uncertainties, both the rate requirement of each FU and the performance protection of MUs should be considered simultaneously. Furthermore, the sensitivity analysis should be given to address the robustness of the scheme. The summary of the related work about resource allocation of HetNets is listed in Table 1.
In this paper, we consider the robust uplink resource allocation algorithm in a two-tier heterogeneous network under the stochastic model. We assume that there is a maco BS serving multiple MUs, and multiple femto BSs intending to serve their own FUs inside the macrocell network’s coverage. In our formulation, we maximize the sum data rate of all FUs under the outage probability constraints that the cross-tier interference of MU and the rate requirement of FU are considered. Our main contributions are summarized as follows:
  • A multiuser resource allocation model for a two-tier HetNet is built. In this paper, we consider that the perfect CSI in both interference links (i.e., from the FUs to MUs) and signal links (i.e., the signal links of femtocell network) are not obtained while the channel uncertainties of the cross-tier interference constraints of MUs and the minimum rate constraint of each FU are formulated into the probability constraints. The existing works only consider the uncertainties in interference links from FUs to MUs, which ignores the effect of channel uncertainties from inter-tier links among FUs. Due to channel uncertainties, the original problem is non-convex and computationally intractable.
  • Furthermore, the original non-convex optimization problem is transformed into two different problems under no CSI feedback and partial CSI. Unlike traditional Bernstein inequality approximation approach with high computational complexity, however, in this paper, both of outage-based interference and rate constraints are transformed into the closed-form approximations by some algebraic transformation. A sub-optimal algorithm is proposed to solve the robust mixed integer programming problem by using the Lagrange dual decomposition and the sub-gradient method.
  • Thirdly, the computational complexity and sensitivity analysis (i.e., the impacts of uncertain parameters) of the proposed algorithm are given in this paper. Simulation results demonstrate the effectiveness of the proposed algorithm in terms of converge performance by comparing with the existing robust power control algorithm.
The rest of this paper is organized as follows. Section 2 introduces the system model. Section 3 gives the robust resource allocation problem. Section 4 gives the transformation of outage probability constraint and shows the robust power allocation algorithm. And, the performance analysis and Simulation results are provided in Section 5, and conclusions are drawn in Section 6.

2. System Model

We consider a multiuser OFDM uplink HetNet comprising M FUs communicating with their femtocell BSs over K subcarriers. The FUs can utilize the spectrum resource of MUs by femtocell BSs opportunistically, as shown in Figure 1. Note that both M and K vary based on the number of active users and available vacant subcarriers dynamically, and indexed by m M { 1 , 2 , , M } and k K { 1 , 2 , , K } respectively. Suppose that both of BSs and user equipments (UEs) have single antenna. Further, we assume K M , and the bandwidth of each subcarrier is assumed to be B Hz which is much less than the coherent bandwidth of wireless channel. As a result, users undergo a flat fading. The notations are summarized in Table 2.
According to the definition of Table 1, based on Shannon theorem, the corresponding data rate of FU m over subcarrier k is formulated as
r m , k = B ρ m , k log 2 1 + p m , k h m , k σ m , k
where the ρ m , k denotes subcarrier assignment of FU m on subcarrier k, which can only be either 1 or 0, indicating whether the subcarrier k is used by the FU m or not.
Due to the limitation of battery capacity of the mth FU transmitter, the transmit power of each FU is not infinite, and the constraint is denoted by
k = 1 K ρ m , k p m , k p m m a x , m M
where p m m a x is the maximum transmit power of FU m.
Meanwhile, data rate should satisfy a minimum requirement to protect the QoS of FU m, given by
k = 1 K r m , k R m m i n , m M
where R m m i n denotes minimum rate requirement of FU m.
The total cross-tier interference constraint from femtocell networks to MU receiver can be expressed as
m = 1 M k = 1 K ρ m , k p m , k g m , k I t h
where I t h is the interference temperature level at MU-receiver.
The sum rate maximization based Power Allocation (PA) problem for HetNets can be written as
max ρ m , k , p m , k m = 1 M k = 1 K r m , k s . t . C 1 : k = 1 K ρ m , k = 1 , m M , C 2 : k = 1 K ρ m , k p m , k p m m a x , m M , C 3 : k = 1 K r m , k R m m i n , m M , C 4 : m = 1 M k = 1 K ρ m , k p m , k g m , k I t h , C 5 : ρ m , k { 0 , 1 } , m M , k K .
C 1 ensures that each subcarrier k is only assigned to one FU, where ρ m , k = 1 means that the kth subcarrier is used by FU m. Moreover, C 2 represents the transmission power constraint of FU m over subcarriers. C 3 can ensure the QoS of each FU. C 4 denotes the total interference power to the MU receiver. The problem with integer variable ρ m , k = 1 is a non-convex and mixed integer programming problem. Assume that the MU can provide the information of channel gains g m , k feedback to the FU, which means that channel gains can be accurately estimated by FUs. Most of current power allocation in HetNets focus on optimal power control under perfect CSI [22,23].
In practice, however, channel uncertainties are inevitable due to estimation errors and quantization errors which cause harmful interference to the MUs [24]. To reduce performance degradation of MUs, PA with channel uncertainties should be considered in advance. However, to the best of our knowledge, robust PA considering the QoS of both FUs and MUs has not been investigated in the existing works.

3. Robust Power Allocation Problem

Due to lack of cooperation between femtocell network and macrocell network, it is a difficult task to precisely estimate the channel gain g m , k from FU-transmitter to MU-receiver. In this paper, our goal is to design a robust PA and subcarrier assignment scheme to ensure the communication quality of FUs and MUs under channel uncertainties. Hence, we reformulated the constraints C 3 and C 4 as probability form. The PA Problem (5) with outage probability constraints is formulated as
max ρ m , k , p m , k m = 1 M k = 1 K r m , k s . t . C 1 , C 2 , C 5 , C 6 : Pr k = 1 K r m , k < R m m i n ξ m , m M , C 7 : Pr m = 1 M k = 1 K ρ m , k p m , k g m , k > I t h ε .
where both of C 6 and C 7 ensure that the QoS of MU and FU are satisfied with a certain probability of outage event that is less than the outage probability threshold ξ m and ε , respectively. In general, Problem (6) is a challenging optimization problem for obtaining analytical solutions because probabilistic rate and interference constraints (i.e., C 6 and C 7 ) are intractable and are not convex. The key step in the development of a tractable robust PA algorithm is to transform the outage probability constraints into efficiently-computable representation (e.g., convex form).
In the following part, we will transform the outage probability constraints into the deterministic and convex ones under imperfect statistical distribution model of uncertain parameter and perfect probability distributions of channel estimation errors.

3.1. Transformation without Statistical Information

In this section, the probabilistic interference constraint and probabilistic rate constraint are converted into the deterministic ones.
Due to the feature of OFDM technology, there is no mutual interference among different subcarriers, so that the data of each FU is mutually independent for all subcarriers. Define the rate set
R k = r m , k R m m i n ,
and
R ¯ = k = 1 K r m , k R m m i n .
According to the above definition, the set R ¯ is a subset of the intersection of R k , i.e.,
R ¯ R = R 1 R 2 R K .
According to probability theory, we have the following relationship,
Pr R ¯ Pr R = k = 1 K Pr R k .
In other words, we have
Pr k = 1 K r m , k R m m i n k = 1 K Pr r m , k R m m i n .
If the upper bound of probabilistic rate constraint satisfies the outage probability requirement, based on the worst-case method, as we know, the stochastic constraint C 6 is satisfied. Therefore, we have
max Pr k = 1 K r m , k R m m i n = k = 1 K Pr r m , k R m m i n ξ m ,
From (12), we have the deterministic form of outage-based probability constraint C 6 , i.e.,
R m m i n B ρ m , k log 2 1 + p m , k σ m , k H h m , k 1 ξ m / K , m M .
If transmission power of SU satisfies the constraint (13), the outage probability of C 6 can be ensured. The proof of (13) is given in Appendix A.
Similarly, the probabilistic interference constraint C 7 can be transformed as
ρ m , k p m , k I t h K G g m , k 1 1 ε M K , m M , k K .
Therefore, the outage probability constraint C 7 becomes a deterministic one. To keep the presentation concise, the proof of (14) is given in Appendix B.
Based on (13) and (14), the PA problem without knowledge of statistical model is given as
max ρ m , k , p m , k m = 1 M k = 1 K r m , k s . t . C 1 , C 2 , C 5 , C 8 : B ρ m , k log 2 1 + p m , k σ m , k H h m , k 1 ξ m / K R m m i n , m M C 9 : K ρ m , k p m , k G g m , k 1 1 ε M K I t h .
Obviously, problem (6) becomes a deterministic one. If the inverse cumulative distribution functions (CDFs) of variables (e.g., g m , k and h m , k ) are well-known (i.e., G g m , k 1 ( . ) and H h m , k 1 ( . ) ), the problem (15) can be easily solved. For example, if the channel gains are follow Rayleigh fading models, the problem can be well solved. In real communication scenario, however, due to the mobility of UEs, it is impractical to assume same fading models during different environments. Therefore, it is necessary to find a general method to solve this problem.

3.2. Transformation with Perfect Statistical Model

In order to address the above problem, it is feasible to find the solution with knowledge of perturbation part of uncertain channel gains.
If there are some errors in the estimated channel gain from feedback quantization or channel estimation [25]. Channel uncertainties can be formulated as an additive model of uncertain parameter [11], i.e.,
g m , k = g ^ m , k + Δ g m , k , m M , k K , h m , k = h ^ m , k + Δ h m , k , m M , k K ,
where g ^ m , k is the estimated channel gain between FU-transmitter and MU-receiver. And h ^ m , k denotes the estimated channel gains from FUs to their BS. Those parameters are perfectly known for FUs. Δ g m , k and Δ h m , k are the corresponding perturbation terms (i.e., estimation error).
It is clear that FUs can obtain the CSI by estimating the channels between FUs and MUs so that the causes of CSI errors are from channel estimation, as a result, this type of estimation error can be formulated as independent Gaussian distribution model [25]. Therefore, uncertain parameter Δ g m , k is reasonably assumed to be distributed random with zero-mean and variance υ m , k 2 , i.e., Δ g m , k C N 0 , υ m , k 2 . However, the elements Δ h m , k in (16) should be assumed to follow the independent uniform distributed random variables, i.e., Δ h m , k [ δ m , k , δ m , k ] , where δ m , k denotes the upper bound of uncertain region. The reason is that the uplink channel gain from FU’s transmitter to the BS is obtained by a quantizer to quantize CSI and feed it back to its corresponding FU’s transmitter.
To protect the performance of MUs, each FU transmits less power for avoiding harmful interference to the MU over link k. Based on C 6 and Model (16), we have
Pr k = 1 K B ρ m , k log 2 1 + p m , k ( h ^ m , k + Δ h m , k ) σ m , k < R m m i n Pr C m B ρ m , k log 2 1 + p m , k ( h ^ m , k + Δ h m , k ) σ m , k R m m i n = Pr Δ h m , k D m , k ξ m ,
where C m denotes the subcarrier set of FU m, and C m represents the number of subcarrier that is chosen by FU m. And D m , k = σ m , k p m , k 2 R m m i n / C m B ρ m , k 1 h ^ m , k .
Due to uniform distribution of channel errors Δ h m , k , we have
2 R m m i n / C m B ρ m , k 1 + p m , k σ m , k h ^ m , k + 2 δ m , k ξ m ,
and
k C m B ρ m , k log 2 1 + p m , k σ m , k h ^ m , k + 2 δ m , k ξ m R m m i n .
As a result, the probabilistic rate constraint becomes a deterministic one as form of (19).
According to the relationship of (16) and constraint C 7 , we have
Pr m = 1 M k = 1 K ρ m , k p m , k Δ g m , k > I ¯ ε ,
where I ¯ = I t h m = 1 M k = 1 K ρ m , k p m , k g ^ m , k is a deterministic term for FU’s transmitter. The Constraint (20) can also be reformulated as
Pr m = 1 M k C m ρ m , k p m , k Δ g m , k I ¯ 1 ε .
To satisfy the outage probability requirement, the interference constraint of left side in (21) must hold under any channel estimation error.
Based on worst-case principle, we have
m = 1 M k C m ρ m , k p m , k Δ g m , k max Δ g m , k R g m = 1 M k C m ρ m , k p m , k Δ g m , k m = 1 M C m ρ m , k p m , k Δ g m , k ,
where k = arg max k g ^ m , k denotes the worst interference link between FU and MU. · denotes the number of element in set.
In other words, if the outage-based constraint can be ensured under the case of worst errors over all subcarriers, the outage performance of MU must be kept. Consequently, we have
Pr m = 1 M C m ρ m , k p m , k Δ g m , k I ¯ 1 ε .
Define B m = C m ρ m , k p m , k and B ˜ m = B m Δ g m , k . Since there is Δ g m , k C N ( 0 , υ m , k 2 ) , the variable B ˜ m follows the Gaussian distribution with zero mean and variance B m 2 υ m , k 2 , i.e., B ˜ m C N 0 , ( B m υ m , k ) 2 . As a result, we have
Pr m = 1 M C m ρ m , k p m , k Δ g m , k I ¯ = Pr m = 1 M B m Δ g m , k I ¯ = Pr m = 1 M B ˜ m I ¯ 1 ε .
Since the sum of Gaussian random variable B ˜ m is still obeying Gaussian distribution, i.e., B ^ = m = 1 M B ˜ m C N 0 , σ 2 and σ = m = 1 M B m υ m , k 2 , therefore, the expression of (24) can be reformulated as
Pr B ^ I ¯ = Q I ¯ 0 σ 1 ε ,
where Q ( x ) is a Gaussian Q-function with the expression of Q ( x ) = 1 π 0 π / 2 exp ( x 2 2 sin 2 θ ) d θ . As a result, we have
I ¯ Q 1 ( 1 ε ) m = 1 M B m υ m , k 2 ,
where Q 1 ( · ) denotes the inverse Gaussian Q-function. Combining with the inequality relationship i x i 2 i x i 2 , we can use the upper bound of right side in (26) to obtain the decomposable form, i.e.,
m = 1 M B m υ m , k 2 m = 1 M B m υ m , k 2
As a result, the outage-based probability constraint C6 becomes a deterministic form
m = 1 M k C m ρ m , k p m , k g ^ m , k + Q 1 ( 1 ε ) υ m , k I t h .
Define r ^ m , k = B ρ m , k log 2 1 + p m , k σ m , k h ^ m , k + 2 δ m , k ξ m , combining (6), (19) and (28), we obtain the considerable robust resource allocation problem as follows
max ρ m , k , p m , k m = 1 M k C m r ^ m , k s . t . C 1 , C 2 , C 5 , C 10 : m = 1 M k C m ρ m , k p m , k g ^ m , k + Q 1 ( 1 ε ) υ m , k I t h , C 11 : k C m r ^ m , k R m m i n , m M .
As a result, the probabilistic optimization problem (6) is transformed into a non-probabilistic one (29).

4. Robust Resource Allocation Algorithm

Obviously, Problem (29) is still not a convex problem due to the integer variable ρ m , k . Since both the real variable p m , k and integer variable ρ m , k are in the optimization problem, it conducts a mixed integer programming problem. Thus, we relax the subcarrier assignment factor into continuous one and introduce a variable s m , k = ρ m , k p m , k for any FU and subcarrier, where ρ m , k [ 0 , 1 ] indicates a time-division multiple access (TDMA) strategy for different FUs. Therefore, the Problem (29) becomes a convex optimization problem as follows,
max ρ m , k , s m , k m = 1 M k C m r ˜ m , k s . t . C 1 , C 12 : k C m s m , k p m m a x , m M , C 13 : k C m r ˜ m , k R m m i n , m M , C 14 : m = 1 M k C m s m , k g ^ m , k + Q 1 ( 1 ε ) υ m , k I t h ,
where r ˜ m , k = B ρ m , k log 2 1 + s m , k ( h ^ m , k + 2 δ m , k ξ m ) σ m , k ρ m , k denotes the effective data rate due to unknown channel gain h m , k (i.e., true physical channel gain). The objective function in (30) is concave since the corresponding hessian matrix can be proved to be negative semi-definiteness by using the Lemma 1 in [26]. In addition, the constraints are linear so that the robust optimization Problem (30) is a convex problem.
We can obtain the analytical solutions by using Lagrangian methods [27]. To deal with the optimization problem (30), we first define a Lagrangian function as
L { s m , k } , { ρ m , k } , { λ m } , { μ m } , { β m } , υ = m = 1 M k C m r ˜ m , k m = 1 M λ m k C m s m , k p m m a x m = 1 M β m k C m ρ m , k 1 υ m = 1 M k C m s m , k g ^ m , k + Q 1 ( 1 ε ) υ m , k I t h m = 1 M μ m R m m i n k C m r ˜ m , k ,
where { λ m } , { μ m } , { β m } and υ are the nonnegative Lagrangian multipliers for the corresponding constraints in (30), respectively. The Lagrangian dual function can be defined as
g { λ m } , { β m } , { μ m } , υ = max { s m , k } , { ρ m , k } L { s m , k } , { ρ m , k } , { λ m } , { β m } , { μ m } , υ .
And the dual problem is formulated as
min { λ m , β m , μ m , υ } g { λ m } , { β m } , { μ m } , υ s . t . λ m 0 , β m 0 , μ m 0 , υ 0 .
The Lagrangian dual Problems (32) and (33) can be decomposed into a master problem and M × K subproblems [27]. As a result, the Lagrangian function can be rewritten as
L { s m , k } , { ρ m , k } , { λ m } , { β m } , { μ m } , υ = m = 1 M k C m L m , k { s m , k } , { ρ m , k } , { λ m } , { β m } , { μ m } , υ + m = 1 M λ m p m m a x m = 1 M R m m i n μ m + m = 1 M β m + I t h υ ,
where
L m , k { s m , k } , { ρ m , k } , { λ m } , { β m } , { μ m } , υ = 1 + μ m r ˜ m , k λ m s m , k β m ρ m , k υ s m , k g ^ m , k + Q 1 ( 1 ε ) υ m , k .
The Karush-Kuhn-Tucker (KKT) conditions of variables s m , k and ρ m , k can be obtained by calculating the derivatives of (35) with the optimal variables, i.e.,
0 s m , k L m , k ( ) s m , k 0 ,
0 ρ m , k L m , k ( ) ρ m , k 0 ,
where ⊥ denotes the orthogonality relationship of the corresponding variables. Additionally, the derivatives L m , k ( ) s m , k and L m , k ( ) ρ m , k are given as
L m , k ( ) s m , k = 1 + μ m r ˜ m , k ( s m , k , ρ m , k ) s m , k λ m υ g ^ m , k + Q 1 ( 1 ε ) υ m , k ,
L m , k ( ) ρ m , k = 1 + μ m r ˜ m , k ( s m , k , ρ m , k ) ρ m , k β m .
As a result, according to KKT conditions, we obtain the optimal transmit power
p m , k = s m , k ρ m , k = 1 ln 2 B ( 1 + μ m ) λ m + υ g ^ + Q 1 ( 1 ε ) υ m , k σ m , k h ^ m , k + 2 δ m , k ξ m + ,
where [ x ] + = max ( 0 , x ) . Moreover, the subcarrier k is assigned to the optimal user m , i.e.,
ρ m , k = 1 m = max m β ^ m , k ,
where
β ^ m , k = 1 + μ m r ˜ m , k ( s m , k , ρ m , k ) ρ m , k = B 1 + μ m log 2 ( 1 + p m , k h ^ m , k σ m , k ) p m , k h ^ m , k ln 2 σ m , k + p m , k h ^ m , k .
The Lagrangian multipliers can be updated as follows
λ m t + 1 = λ m t d 1 t p m m a x k C m p m , k t + ,
μ m t + 1 = μ m t d 2 t k C m r ˜ m , k t R m m i n + ,
υ t + 1 = υ t d 3 t I t h m = 1 M k C m p m , k t g ^ m , k + Q 1 ( 1 ε ) υ m , k + .
where t denotes the iteration number. d 1 , d 2 and d 3 are the corresponding small step sizes. When the step sizes are sufficiently small, Lagrange multipliers can converge to equilibrium points [28]. The implementation procedure of the algorithm is given in Algorithm 1 and the algorithm flow chart is shown in Figure 2.
Algorithm 1: The implementation procedure.
1:
Initialize maximum iteration number T m a x ; Set: iterations t = 0 , M > 0 , K > 0 ; Lagrangian multipliers λ m ( 0 ) > 0 , μ m ( 0 ) > 0 , and υ ( 0 ) > 0 ; Outage probability thresholds ε ( 0 , 1 ) , ξ m ( 0 , 1 ) ; Upper bound of channel estimation error in FU’s link is δ m , k [ 0 , 1 ] ; Variance of estimation error in FU-to-MU link is υ m , k [ 0 , 1 ] ;
2:
Set maximum transmit power p m m a x > 0 and initialize power p m , k > 0 with same initialization values among all subcarriers. Define interference I t h , randomly generate g ^ m , k and h ^ m , k .
3:
Initialize ρ m , k with subcarrier assignment method in [29], group subcarrier set C m and calculate C m .
4:
repeat
5:
for t = 1 to T m a x do
6:
for m = 1 to M do
7:
  for k = 1 to K do
8:
   Calculate transmit power p m , k according to (40)
9:
   Calculate β ^ m , k by (42)
10:
   Calculate ρ m , k by using (41)
11:
   Calculate Lagrange multipliers λ m , μ m , and υ from (4345)
12:
  end for
13:
end for
14:
t = t + 1
15:
end for
16:
until t = T m a x or transmit power convergence

5. Performance Analysis

In this section, we will analyze the complexity, sensitivity and simulation results of the proposed algorithm.

5.1. Complexity Analysis

The the computational complexity of the proposed algorithm is analyzed in this subsection. Obviously, the calculation of (38) for each femtocell user on each subcarrier in each femtocell k implies M K operations, so that the worst-case complexity of search (37) (i.e., find the optimal ρ m , k ) needs M K operations in each iteration. If the subgradient methods used in (39)–(41) require t m iterations to converage, then the updates of Lagrange multipliers λ m and μ m for m , k need O M K operations [30], and the computation of υ needs O 1 , so that the t m is a polynomial function of M 2 K 2 . Thus the total complexity of the algorithm is O M 2 K 2 t m [27]. Additionally, t m can be small enough when the initial values of Lagrange multipliers λ m ( 0 ) , μ m ( 0 ) and υ ( 0 ) are well chosen [6,22].

5.2. Sensitivity Analysis

In this section, we will give the exact perturbation version of sum rate gap under uncertainties Δ g m , k and Δ h m , k . The reduction in the sum transmission rate of FUs can be formulated as
U g a p = m = 1 M k C m 2 μ m δ m , k ξ m + υ Q 1 ( 1 ε ) υ m , k ,
where μ m and υ denotes optimal Lagrangian multipliers. The proof is given in Appendix C.

5.3. Simulation Results

In this subsection, we use computer simulations to demonstrate the performance of the proposed algorithm in different scenarios. The simulation parameters are given in Table 3.
Figure 3 shows the convergence performance of the proposed algorithm in terms of transmit power and interference power versus the number of iterations t. We can see that the proposed algorithm takes only small iterations to converge, which indicates that it has a good real-time performance for practical applications. The transmit power of each FU is constrained by the maximum transmit power level. In addition, our proposed algorithm can well protect the performance of MU because the corresponding interference power does not exceed the interference temperature threshold.
In Figure 4, the transmission rate of FUs is depicted as a function of ε under different channel estimation errors Δ g m , k (i.e., the corresponding variance υ m , k ). The outage probability of FU is defined as ξ m = 0.1 . The variance of channel estimation error Δ h m , k is δ m , k = 0.01 . As can be seen, the total data rate of FUs increases as the maximum transmit power threshold p m m a x increases. This can be explained, as for lower values of p m m a x , the total transmit power of FUs p m , k , k is limited. Increasing the transmit power threshold p m m a x can increase the feasible region of transmit power p m , k . Hence, it enables the proposed algorithm to improve the sum data of FUs. In addition, the total data rate of FUs under the higher outage threshold of MU ε is higher than that of lower outage probability, e.g., rate ε = 0.2 , υ m , k = 0.5 > rate ε = 0.1 , υ m , k = 0.5 . The reason is that, where MUs allow the bigger outage probability, the communication quality of MUs means that it is not easy to be interrupted. Therefore, it allows that FUs transmit more power to improve their performance. Moreover, in the case of same outage probability, we can find that the sum rate of FUs with small variance of estimation error is better than that of bigger variance, e.g., υ m , k = 0.2 . Because the bigger error variance means that the estimated channel is not exact, the channel estimation value g ^ m , k deviates from the true channel gain g m , k seriously.
Figure 5 depicts the effect of FU’s outage probability ξ m and the variance of channel estimation error ( Δ h m , k ) on the FU performance. The outage probability and the variance of channel uncertainty of MU are defined as ε = 0.1 and υ m , k = 0.01 . The sum data rate of FUs increases with the increasing transmit power level p m m a x . In addition, as the value of δ m , k increases, i.e., the estimation error increases, the transmission rate of FUs increases accordingly, since it needs more transmit power to overcome the effect of channel uncertainty so that the basic rate requirement of each FU can be satisfied. Additionally, the data rate of FUs increases with the bigger value of ε m .
Figure 6 shows the interference power to the MU under different channel uncertainties. The outage probability of MU is ε = 0.1 . The outage probability and variance of FU are defined as ξ m = 0.1 and δ m , k = 0.1 separately. From Figure 6, the interference power to the MU of the proposed algorithm and the non-robust algorithm in [6] increases with the bigger transmit power threshold p m m a x . Because the bigger transmit power can provide more wider feasible region. Additionally, the interference power received by MU in [6] exceeds the threshold I t h . But the interference to the MU of our method is below the interference power threshold. Additionally, the bigger upper bound of uncertainty υ m , k is, more estimation errors are in the communication system. Therefore, the interference power to MUs under υ m , k = 0.1 is bigger than the interference power under υ m , k = 0.01 .
In Figure 7, the achievable data rate of FUs versus the total power level for the above mentioned schemes is given under different variances of channel estimation errors over FU-to-MU links. From Figure 7, the sum data rate of FUs under the proposed algorithm and the non-robust algorithm increases with the bigger maximum transmit power level. Additionally, as expected, the sum data of FUs under the non-robust algorithm in [6] is higher than our proposed algorithm for a given transmit power upper bound if perfect CSI is assumed, because the allowable transmit power is higher. Moreover, the sum data of FUs converges to an equilibrium point (i.e., saturation state). The reason is that the optimal power is limited by the interference power threshold under the higher transmit power region. In order to protect the QoS of MUs and the limitation of maximum transmit power, it does not allow the FUs to increase their transmission power endlessly.

6. Conclusions

In this paper, to well protect the QoS of macrocell cellular networks and improve system capacity, we have formulated a sum rate maximization problem for a two-tier heterogeneous network with one macrocell and multiple femtocells. Due to the effect of channel uncertainties in inter-tier interference channels and the forward channel of femtocell user, we have imposed the robust interference constraint and the robust rate constraint under outage-based uncertainty models to protect the transmission performance of the macrocell and the femtocells. With the scenarios of no CSI feedback and partial CSI feedback, the outage-based probability constraints are transformed into closed-from expressions by approximate methods. With the Lagrange dual method, the robust power allocation problem has been decomposed into a dual problem. The robust power allocation algorithm is proposed by using the sub-gradient method to solve the primal problem. Simulation results have demonstrated that the proposed algorithm can fast converge to the optimal value. The efficient resource allocation with the consideration of channel uncertainties in heterogeneous networks will be investigated in our future work. Green communication emphasizes on incorporating energy awareness in communication systems. Therefore, energy-efficient RA has attracted more attention in both industry and academia recently. In the future, we will focus on the issue of energy efficiency improvement.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (Grant No. 61601071, 61301124), Open Foundation of ISN State Key Lab of Xidian University (Grant No. ISN17-01), Scientific and Technological Research Program of Chongqing Municipal Education Commission (Grant No. KJ1600412), Municipal Natural Science Foundation of Chongqing (Grant No. CSTC2016JCYJA2197), and Dr. Startup Founds of Chongqing University of Posts and Telecommunications (Grant No. A2016-12).

Author Contributions

Yongjun Xu designed the algorithm, performed the simulations, and prepared the manuscript. Guoquan Li is responsible for coordinating and writing. Yuan Hu is revising the paper. Both authors discussed the results and approved the publication.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

According to (1) and (8), we have
k = 1 K Pr B ρ m , k log 2 1 + p m , k h m , k σ m , k R m m i n ξ m ,
and
k = 1 K Pr h m , k σ m , k p m , k 2 R m m i n / ( B ρ m , k ) 1 ξ m .
If we do not know the distribution function of channel gain h m , k , similarly, we have the following result
k = 1 K H h m , k Z m , k ξ m ,
where Z m , k = σ m , k p m , k 2 R m m i n / ( B ρ m , k ) 1 . And H h m , k · denotes the CDF of channel h m , k . Since we can not directly obtain the relationship between transmit power p m , k and outage probability ξ m , therefore, we reform the constraint into the following approximation type, i.e.,
k = 1 K H h m , k Z m , k K H h m , k ¯ Z m , k ¯ ,
where k ¯ = arg max H h m , k Z m , k is a channel selection factor, which is equivalent to
k ¯ = arg min k h m , k σ m , k .
In practice, considering that the outage probability requirement of each subcarrier can be satisfied under the worst-case channel environment, the sum of outage set in (A2) can limit the outage probability under the threshold ξ m . Therefore, from (A3) to (A5), we obtain the following expression
K H h m , k ¯ Z m , k ¯ = Δ K H h m , k σ m , k p m , k 2 R m m i n / ( B ρ m , k ) 1 ξ m ,
σ m , k 2 R m m i n / ( B ρ m , k ) 1 H h m , k 1 ξ m / K p m , k .
Therefore, if transmit power of SU over each subcarrier satisfies the constraint (A7), the outage constraint in C 7 can be ensured. For the sake of simplification, the constraint (A7) also can be expressed as
R m m i n B ρ m , k log 2 1 + p m , k σ m , k H h m , k 1 ξ m / K .
The proof is completed.

Appendix B

To prove the inequality (10), the upper bound of outage-based interference constraint satisfies
Pr k = 1 K m = 1 M ρ m , k p m , k g m , k I t h Pr max m , k ρ m , k p m , k g m , k I t h K = 1 Pr max m , k ρ m , k p m , k g m , k I t h K = 1 k = 1 K m = 1 M Pr ρ m , k p m , k g m , k I t h K = 1 Pr ρ m , k p m , k g m , k I t h K M K ε .
Since the outage probability level ε satisfies the interval of [0,1], we have
1 ε M K Pr ρ m , k p m , k g m , k I t h K .
Without the knowledge of PDF of channel gain g m , k , we have
1 ε M K G g m , k I t h ρ m , k p m , k K ρ m , k p m , k I t h K G g m , k 1 1 ε M K ,
where G g m , k · and G g m , k 1 · denotes the cumulative distribution function (CDF) and the corresponding inverse function of variable g m , k , respectively. The proof is completed.

Appendix C

Based on the formula of Taylor series of the two element function, we have
u h ^ m , k + Δ h m , k , g ^ m , k + Δ g m , k = u h ^ m , k , g ^ m , k + m = 1 M k C m u h ^ m , k , g ^ m , k + Δ g m , k Δ h m , k Δ h m , k + o ( Δ h m , k ) n + m = 1 M k C m u h ^ m , k + Δ h m , k , g ^ m , k Δ g m , k Δ g m , k + o ( Δ g m , k ) n , ( Δ h m , k 0 , Δ g m , k 0 ) ,
where o ( Δ g m , k ) n and o ( Δ h m , k ) n denote the corresponding high order infinitesimal small quantities. In addition, u h ^ m , k , g ^ m , k is the optimal utility function without estimation errors (assuming that the estimated channel gains are equal to the true channel gains).
Ignoring the effect of high order small variables, the sum rate gap between our robust power allocation algorithm and the non-robust algorithm can be written as
U g a p = u r o b u n o n = m = 1 M k C m u h ^ m , k , g ^ m , k + Δ g m , k Δ h m , k Δ h m , k + m = 1 M k C m u h ^ m , k + Δ h m , k , g ^ m , k Δ g m , k Δ g m , k ,
where u r o b = u h ^ m , k + Δ h m , k , g ^ m , k + Δ g m , k is the total data rate with imperfect CSI and u n o n = u h ^ m , k , g ^ m , k is the total rate with perfect CSI. According to the sensitivity analysis in [32], under the consideration of small uncertainties, we have
u h ^ m , k , g ^ m , k + Δ g m , k Δ h m , k μ m u h ^ m , k + Δ h m , k , g ^ m , k Δ g m , k υ ,
According to the relationship (A13) and (A14), we have the following expresssion
U g a p = u r o b u n o n = m = 1 M k C m 2 μ m δ m , k ξ m + υ Q 1 ( 1 ε ) υ m , k ,
where Δ h m , k = 2 δ m , k ξ m and Δ g m , k = Q 1 ( 1 ε ) υ m , k . The proof is completed.

References

  1. Wang, F.; Chen, W.; Wu, Q.Q.; Tang, H.Y. Robust MISO beamforming under the deterministic model in two-tier heterogeneous networks. IEEE Access 2017, 5, 14616–14625. [Google Scholar] [CrossRef]
  2. Andrews, J.G.; Claussen, H.; Dohler, M.; Rangan, S.; Reed, M.C. Femtocells: Past, Present, and Future. IEEE J. Sel. Area. Commun. 2012, 30, 497–508. [Google Scholar] [CrossRef]
  3. Lopez-Perez, D.; Guvenc, I.; De la Roche, G.; Kountouris, M.; Quek, T.Q.S.; Zhang, J. Enhanced intercell interference coordination challenges in heterogeneous networks. IEEE Wirel. Commun. 2011, 18, 22–30. [Google Scholar] [CrossRef]
  4. Wang, H.N.; Wang, J.H.; Ding, Z. Distributed power control in a two-tier heterogeneous network. IEEE Trans. Wirel. Commun. 2015, 14, 6509–6523. [Google Scholar] [CrossRef]
  5. Romero, J.; Sanchez-Gonzalez, J.; Agusti, R.; Lorenzo, B.; Glisic, S. Power-efficient resource allocation in a heterogeneous network with cellular and D2D capabilities. IEEE Trans. Veh. Technol. 2016, 65, 9272–9286. [Google Scholar] [CrossRef]
  6. Zhang, H.J.; Jiang, C.; Beaulieu, N.C.; Chu, X.; Wen, X.; Tao, M. Resource allocation in spectrum-sharing OFDMA femtocells with heterogeneous services. IEEE Trans. Commun. 2014, 62, 2366–2377. [Google Scholar] [CrossRef]
  7. Jo, H.S.; Mun, C.; Moon, J.; Yook, J.G. Interference mitigation using uplink power control for two-tier femtocell networks. IEEE Trans. Wirel. Commun. 2009, 8, 4906–4910. [Google Scholar]
  8. Mao, T.; Feng, G.; Liang, L.; Qin, S.; Wu, B. Distributed energy-efficient power control for macrofemto networks. IEEE Trans. Veh. Technol. 2016, 65, 718–731. [Google Scholar] [CrossRef]
  9. Gershman, A.B.; Sidiropoulos, N.D. Space-Time Processing for MIMO Communications; John Wiley and Sons: Hoboken, NJ, USA, 2005. [Google Scholar]
  10. Ben-Tal, A.; Nemirovski, A. Selected Topics in Robust Convex Optimization. Math. Program. 2007, 112, 125–158. [Google Scholar] [CrossRef]
  11. Xu, Y.J.; Zhao, X.H.; Liang, Y.C. Robust power control and beamforming in cognitive radio networks: A survey. IEEE Commun. Surv. Tutor. 2015, 17, 1834–1857. [Google Scholar] [CrossRef]
  12. Vaezpour, E.; Dehghan, M.; Yousefizadeh, H. Robust distributed resource allocation in OFDMA femtocell networks. Comput. Commun. 2017, 109, 1–12. [Google Scholar] [CrossRef]
  13. Zhu, K.; Hossain, E.; Anpalagan, F.A. Downlink power control in two-tier cellular OFDMA networks under uncertainties: A robust stackelberg game. IEEE Trans. Commun. 2015, 63, 1123–1129. [Google Scholar] [CrossRef]
  14. Xiao, S.; Zhou, X.G.; Wu, Y.Y.; Li, G.Y.; Guo, W. Robust resource allocation in full-duplex-enabled OFDMA femtocell networks. IEEE Trans. Wirel. Commun. 2017, 16, 235–242. [Google Scholar] [CrossRef]
  15. Parsaeefard, S.; Van der Schaar, M.; Sharafat, A.R. Robust power control for heterogeneous users in shared unlicensed bands. IEEE Trans. Wirel. Commun. 2014, 13, 3167–3182. [Google Scholar] [CrossRef]
  16. Xu, Y.J.; Hu, Y.; Chen, Q.; Song, T.; Lai, R. Robust resource allocation for multi-tier cognitive heterogeneous networks. In Proceedings of the IEEE International Conference on Communications (ICC), Paris, France, 21–25 May 2017; pp. 1–6. [Google Scholar]
  17. Parsaeefard, S.; Sharafat, A.R. Robust distributed power control in cognitive radio networks. IEEE Trans. Mob. Comput. 2013, 12, 213–220. [Google Scholar] [CrossRef]
  18. Liu, Z.X.; Li, S.Y.; Ka, M.; Guan, X.P.; Li, X.B. Robust power allocation based on hierarchical game with consideration of different user requirements in two-tier femtocell networks. Comput. Netw. 2017, 122, 179–190. [Google Scholar] [CrossRef]
  19. Liu, Z.X.; Zhang, P.; Guan, X.P.; Yang, H.J. Robust power control for femtocell networks under outage-based QoS constraints. Comput. Netw. 2016, 102, 145–156. [Google Scholar] [CrossRef]
  20. Mokari, N.; Saeedi, H.; Azmi, P. Quantized ergodic radio resource allocation in cognitive femto networks with controlled collision and power outage probabilities. IEEE J. Sel. Areas Commun. 2014, 32, 2090–2104. [Google Scholar] [CrossRef]
  21. Wang, K.Y.; Jacklin, N.; Ding, Z.; Chi, C.Y. Robust MISO transmit optimization under outage-based QoS constraints in two-tier heterogeneous networks. IEEE Trans. Wirel. Commun. 2013, 12, 1883–1897. [Google Scholar] [CrossRef]
  22. Zhang, H.J.; Liu, H.; Cheng, J.; Leung, V.C.M. Downlink energy efficiency of power allocation and wireless backhaul bandwidth allocation in heterogeneous small cell networks. IEEE Trans. Commun. 2017, 95, 1–13. [Google Scholar]
  23. Arani, A.H.; Mehbodniya, A.; Omidi, M.J.; Adachi, F.; Saad, W.; Guven, I. Distributed learning for energy efficient resource management in selforganizing heterogeneous networks. IEEE Trans. Veh. Technol. 2017, 66, 9287–9303. [Google Scholar] [CrossRef]
  24. Zhang, Y.; Wang, S. Resource allocation for cognitive radio-enabled femtocell networks with imperfect spectrum sensing and channel uncertainty. IEEE Trans. Veh. Technol. 2016, 65, 7719–7728. [Google Scholar] [CrossRef]
  25. Shenouda, M.B.; Davidson, T.N. Probabilistically-constrained approaches to the design of the multiple antenna downlink. In Proceedings of the IEEE 42nd Asilomar Conference on Signals, Systems and Computers, Pacific Grove, CA, USA, 2–5 October 2008; pp. 1120–1124. [Google Scholar]
  26. Zhang, H.J.; Jiang, C.; Mao, X.; Chen, H.H. Interference-limited resource optimization in cognitive femtocells with fairness and imperfect spectrum sensing. IEEE Trans. Veh. Technol. 2016, 65, 1761–1771. [Google Scholar] [CrossRef]
  27. Boyd, S.; Vandenberghe, L. Convex Optimization; Cambridge University Press: Cambridge, UK, 2004. [Google Scholar]
  28. Horst, R.; Pardalos, P.M.; Thoai, N.V. Introduction to Global Optimization, 2nd ed.; Kluwer Academic Publishers: Boston, MA, USA, 2000. [Google Scholar]
  29. Shen, Z.; Andrews, J.G.; Evans, B.L. Adaptive resource allocation in multiuser OFDM systems with proportional rate constraints. IEEE Trans. Wirel. Commun. 2005, 4, 2726–2737. [Google Scholar] [CrossRef]
  30. Yu, W.; Lui, R. Dual methods for nonconvex spectrum optimization of multicarrier systems. IEEE Trans. Commun. 2006, 54, 1310–1322. [Google Scholar] [CrossRef]
  31. Setoodeh, P.; Haykin, S. Robust transmit power control for cognitive radio. Proc. IEEE 2009, 97, 915–939. [Google Scholar] [CrossRef]
  32. Cacuci, D.G. Sensitivity and Uncertainty Analysis, Volume II: Applications to Large-scale Systems; Chapman and Hall/CRC: Boca Raton, FL, USA, 2003. [Google Scholar]
Figure 1. An example of the two-tier HetNet consisting of one macrocell and multiple femtocells.
Figure 1. An example of the two-tier HetNet consisting of one macrocell and multiple femtocells.
Sensors 18 00639 g001
Figure 2. Algorithm flow chart.
Figure 2. Algorithm flow chart.
Sensors 18 00639 g002
Figure 3. Convergence performance of the proposed algorithm.
Figure 3. Convergence performance of the proposed algorithm.
Sensors 18 00639 g003
Figure 4. Sum data rate of FUs versus the maximum transmit power under different MU outage probabilities.
Figure 4. Sum data rate of FUs versus the maximum transmit power under different MU outage probabilities.
Sensors 18 00639 g004
Figure 5. Sum data rate of FUs versus the maximum transmit power under different FU outage probabilities.
Figure 5. Sum data rate of FUs versus the maximum transmit power under different FU outage probabilities.
Sensors 18 00639 g005
Figure 6. Comparison of interference power to the MU versus maximum transmit power.
Figure 6. Comparison of interference power to the MU versus maximum transmit power.
Sensors 18 00639 g006
Figure 7. The effect of channel uncertainty on the data rate of FUs.
Figure 7. The effect of channel uncertainty on the data rate of FUs.
Sensors 18 00639 g007
Table 1. A summary of the related work about RA algorithm of HetNets.
Table 1. A summary of the related work about RA algorithm of HetNets.
CSISystem ModelProblem FocusCommon Solution Approach
PerfectUplink[6,7]Rate maximization[6,7]Dual decomposition theory[6]
Interference mitigation strategy[7]
Downlink[6,8]EE maximization[8]
Game strategy[8]
ImperfectWorst-case[12,13,14,15,16]Rate maximization[12,13,14,16,20]MINLP[12]
Stackelberg[13,15,18,19]
Power minimization[15]Lagrangian game theory[14]
Stochastic model[18,19,20,21]Geometrical programming method[16]
Increase social utility[18,19,21]Iterative approach[20]
Semidefinite programming[21]
Problem analysisFrom the aspect of multiuser network scenario, the robust SINR of FU and robust interference constraint of MU are not considered simultaneously, which can not be ignored.The problem considered in this paper is more comprehensive and complex, more satisfying the heterogeneous network robust transmission scenario.
Table 2. Notation of Symbol.
Table 2. Notation of Symbol.
SymbolExplanation
ρ m , k subcarrier assignment of FU m on subcarrier k
σ m , k background noise of FU m on subcarrier k
p m , k transmit power of the mth FU over subcarrier k
h m , k direct channel gain of FU m over subcarrier k
r m , k data rate of FU m over subcarrier k
I t h interference temperature level at MU-receiver
g m , k channel gain of FU-to-MU link on subcarrier k
p m m a x maximum transmit power of FU m
R m m i n minimum rate requirement of FU m
Bbandwidth of each subcarrier
Table 3. Simulation parameters.
Table 3. Simulation parameters.
System ParametersValues
Number of FUs M4
Number of subcarriers K128
Bandwidth of each subcarrier B10 KHz
The background noise σ m , k 1 × 10 8 mW
Allowable interference level I t h 1 × 10 5 mW
Estimated channel gains g ^ m , k , h ^ m , k [0,1] [31]
Outage probability threshold of MU ε [0,1]
Outage probability threshold of FU ξ m [0,1]
Minimum data requirement R m m i n 2 × 10 5 bit/s
Maximum transmit power p m m a x 5 mW

Share and Cite

MDPI and ACS Style

Xu, Y.; Hu, Y.; Li, G. Robust Rate Maximization for Heterogeneous Wireless Networks under Channel Uncertainties. Sensors 2018, 18, 639. https://doi.org/10.3390/s18020639

AMA Style

Xu Y, Hu Y, Li G. Robust Rate Maximization for Heterogeneous Wireless Networks under Channel Uncertainties. Sensors. 2018; 18(2):639. https://doi.org/10.3390/s18020639

Chicago/Turabian Style

Xu, Yongjun, Yuan Hu, and Guoquan Li. 2018. "Robust Rate Maximization for Heterogeneous Wireless Networks under Channel Uncertainties" Sensors 18, no. 2: 639. https://doi.org/10.3390/s18020639

APA Style

Xu, Y., Hu, Y., & Li, G. (2018). Robust Rate Maximization for Heterogeneous Wireless Networks under Channel Uncertainties. Sensors, 18(2), 639. https://doi.org/10.3390/s18020639

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop