Enhanced Group Sparse Beamforming for Green CloudRAN: A Random Matrix Approach
Abstract
Group sparse beamforming is a general framework to minimize the network power consumption for cloud radio access networks (CloudRANs), which, however, suffers high computational complexity. In particular, a complex optimization problem needs to be solved to obtain the remote radio head (RRH) ordering criterion in each transmission block, which will help to determine the active RRHs and the associated fronthaul links. In this paper, we propose innovative approaches to reduce the complexity of this key step in group sparse beamforming. Specifically, we first develop a smoothed minimization approach with the iterative reweighted algorithm to return a KarushKuhnTucker (KKT) point solution, as well as enhancing the capability of inducing group sparsity in the beamforming vectors. By leveraging the Lagrangian duality theory, we obtain closedform solutions at each iteration to reduce the computational complexity. The wellstructured solutions provide the opportunities to apply the largedimensional random matrix theory to derive deterministic approximations for the RRH ordering criterion. Such an approach helps to guide the RRH selection only based on the statistical channel state information (CSI), which does not require frequent update, thereby significantly reducing the computation overhead. Simulation results shall demonstrate the performance gains of the proposed minimization approach, as well as the effectiveness of the large system analysis based framework for computing RRH ordering criterion.
I Introduction
Network densification [2, 3, 4] has been proposed as a promising way to provide ultrahigh data rates, achieve low latency, and support ubiquitous connectivity for the upcoming 5G networks [5]. However, to fully harness the benefits of dense wireless networks, formidable challenges arise, including interference management, radio resource allocation, mobility management, as well as high capital expenditure and operating expenditure. CloudRAN emerges as a disruptive technology to deploy costeffective dense wireless networks [6, 7]. It can significantly improve both energy and spectral efficiency, by leveraging recent advances in cloud computing and network function virtualization [8]. With shared computation resources at the cloud data center and distributed lowcost lowpower remote radio heads (RRHs), CloudRAN provides an ideal platform to achieve the benefits of network cooperation and coordination [9]. There is a unique characteristic in CloudRANs, namely, the highcapacity fronthaul links are required to connect the cloud center and RRHs [10, 11]. Such links will consume power comparable to that of each RRH, and thus brings new challenges to design green CloudRANs. To address this issue, a new performance metric, i.e., the network power consumption, which consists of both the fronthaul link power consumption and RRH transmit power consumption, has been identified in [7] for the green CloudRAN design.
Unfortunately, the network power minimization problem turns out to be a mixedinteger nonlinear programming problem [12], which is highly intractable. Specifically, the combinatorial composite objective function involves a discrete component, indicating which RRHs and the corresponding fronthaul links should be switched on, and a continuous component, i.e., beamformers to reduce RRH transmit power. To address this challenge, a unique group sparsity structure in the optimal beamforming vector has been identified in [7] to unify RRH selection and beamformer optimization. Accordingly, a novel threestage group sparse beamforming framework was proposed to promote group sparsity in the solution. Specifically, a mixed minimization approach was first proposed to induce group sparsity in the solution, thereby guiding the RRH selection via solving a sequence of convex feasibility problems, followed by coordinated beamforming to minimize the transmit power for the active RRHs.
Although the group sparse beamforming framework provides polynomialtime complexity algorithms via convex optimization, it has the limited capability to enhance group sparsity compared to the nonconvex approaches [13, 14]. In particular, the smoothed minimization supported by the iterative reweighted algorithm was developed in [15] to enhance sparsity for multicast group sparse beamforming. However, the computation burden of this iterative algorithm is still prohibitive in dense wireless networks with a large number of RRHs and mobile users. In [16], a generic twostage approach was proposed to solve largescale convex optimization problems in dense wireless cooperative networks via matrix stuffing and operator splitting, i.e., the alternating direction method of multipliers (ADMM) [17, 18]. This approach also enables parallel computing and infeasibility detection. However, as the proposed solutions need to be accomplished for each channel realization, and also due to the iterative procedures of the group sparse beamforming framework, it is still computationally expensive.
In this paper, we improve the performance of group sparse beamforming [7] via minimization, and a special emphasis is on reducing the computational complexity, with two ingredients: closedform solutions for each iteration, and the deterministic equivalents of the optimal Lagrange multipliers. Specifically, we first develop an iterative reweighted algorithm with closedform solutions at each iteration via leveraging the principles of the majorizationminimization (MM) algorithm [19] and the Lagrangian duality theory [20]. It turns out that the proposed iterative reweighted algorithm can find a KKT point for the smoothed nonconvex minimization problem. Furthermore, this reveals the explicit structures of the optimal solutions to each subproblem in the iterative reweighted algorithm, while the convex optimization approach in [21, 7, 15] fails to obtain the closedform solutions. Thereafter, the wellstructured closedform solutions provide opportunities to leverage the largedimensional random matrix theory [22, 23, 24] to perform the asymptotic analysis in the large system regimes [25, 26]. Specifically, the deterministic equivalents of the optimal Lagrange multipliers are derived based on the recent results of large random matrix analysis [27] to decouple the dependency of system parameters. These results are further used to perform an asymptotic analysis for computing the RRH ordering criterion, which only depends on longterm channel attenuation, thereby significantly reducing the computation overhead compared with previous algorithms that heavily depends on instantaneous CSI [28, 29, 30, 31, 7, 21, 15].
Based on the above proposal, we provide a threestage enhanced group sparse beamforming framework for green CloudRAN via random matrix theory. Specifically, in the first stage, we compute the enhanced RRH ordering criterion only based on the statistical CSI. This thus avoids frequent updates, thereby significantly reducing computational complexity, which is a sharp difference compared to the original proposal in [8]. This new algorithm is based on the principles of the MM algorithm and the Lagrangian duality theory, followed by the random matrix theory. With the obtained deterministic equivalents of RRH ordering criterion, a twostage largescale convex optimization framework [16] is adopted to solve a sequence of convex feasibility problems to determine the active RRHs in the second stage, as well as solving the transmit power minimization problem for the active RRHs in the final stage. Note that global instantaneous CSI is required for these two stages. Simulation results are provided to demonstrate the improvement of the enhanced group sparse beamforming framework. Moreover, the deterministic approximations turn out to be accurate even in the finitesized systems.
Ia Related Works
IA1 Sparse Optimization in Wireless Networks
The sparse optimization paradigm has recently been popular for complicated network optimization problems in wireless network design, e.g., the group sparse beamforming framework for green CloudRAN design [7, 30, 31, 32, 15], wireless caching networks [33, 34], user admission control [15, 35], as well as computation offloading [36]. In particular, the convex relaxation approach provides a principled way to induce sparsity via the minimization [35] for individual sparsity inducing and the mixed minimization for group sparsity inducing [7]. The reweighted algorithm [13] and the reweighted algorithm [15] were further developed to enhance sparsity. To enable parallel and distributed computing, the firstorder method ADMM algorithm was adopted to solve the group sparse beamforming problems [37]. A generic largescale convex optimization framework was further proposed to solve general largescale convex programs in dense wireless networks to enable scalability, parallel computing and infeasibility detection [16].
However, all the above algorithms need to be computed for each channel realization, which is computationally expensive. In this paper, we adopt large system analysis to compute the RRH ordering criterion for network adaptation only based on statistical CSI.
IA2 Large System Analysis via Random Matrix Theory
Random matrix theory [22] has been proven powerful for performance analysis, and the understanding and improving algorithms in wireless communications [38, 23], signal processing [39], and machine learning [40, 41], especially in large dimensional regimes for applications in the era of big data [42]. In dense wireless networks, random matrix theory provides a powerful way for performance analysis and algorithm design. In particular, the large system analysis was performed for simple precoding schemes, e.g., regularized zeroforcing in MISO broadcast channels with imperfect CSI [24]. A random matrix approach to the optimal coordinated multicell beamforming for massive MIMO was presented in [43] without closeform expressions for the optimal Lagrange multipliers. A simple channel model with two cells under different coordination levels was considered in [25].
However, all of the above results cannot be directly applied in the group sparse beamforming framework in dense CloudRAN due to the general channel models with heterogenous pathloss and the complicated beamformer structures with fully cooperative transmission. To determine the RRH ordering criterion based on statistical CSI, we adopt the technique in [27, 26] to compute the closedforms for the optimal Lagrange multipliers for each iteration in the procedure of reweighted minimization for sparsity inducing.
IB Organization
The remainder of the paper is organized as follows. Section II presents the system model and problem formulation, followed by performance analysis with the proposed threestage enhanced group sparse beamforming framework for green CloudRAN. In Section III, the iterative reweighted algorithm for group sparse beamforming is presented via the principle of MM algorithm and duality theory. The large system analysis for RRH ordering is performed in Section IV. Simulation results will be demonstrated in Section V. Finally, conclusions and discussions are presented in Section VI. To keep the main text clean and free of technical details, we divert most of the proofs to the Appendix.
Notations
Throughout this paper, is the norm. stands for either the size of a set or the absolute value of a scalar, depending on the context. is the spectral radius of the Hermitian matrix . Boldface lower case and upper case letters represent vectors and matrices, respectively. and denote the inverse, transpose, Hermitian and trace operators, respectively. We use to represent complex domain. denotes the expectation of a random variable. We denote and as a diagonal matrix of order and the identity matrix of order , respectively.
Ii System Model and Problem Formulation
Iia System Model
Consider a CloudRAN with RRHs and singleantenna mobile users (MUs), where each RRH is equipped with antennas. In CloudRANs, the BBU pool will perform the centralized signal processing and is connected to all the RRHs via highcapacity and lowlatency fronthaul links. In this paper, we will focus on the downlink signal processing. Specifically, let be the transmit beamforming vector from the th RRH to the th MU. The received signal at MU is given by
(1) 
where is the channel propagation between MU and RRH , with is the encoded transmission information symbol for MU , and is the additive Gaussian noise at MU .
We assume that ’s and ’s are mutually independent and all the MUs apply singleuser detection. The signaltointerferenceplusnoise ratio (SINR) for MU is given by
(2) 
where is the vector consisting of the channel coefficients from all the RRHs to MU , and is the aggregative beamforming vector for the MU from all the RRHs.
IiB Problem Formulation
In this paper, we aim at designing a green CloudRAN by minimizing the network power consumption, which consists of the fronthaul links power consumption, as well as the RRH transmit power consumption [7]. Specifically, let be the aggregative beamforming vector from all the RRHs to all the MUs. Define the support of the beamforming vector as , where is indexed by with . The relative fronthaul network power consumption is given by
(3) 
where is the relative fronthaul link power consumption [7] (i.e., the static power saving when both the fronthaul link and the corresponding RRH are switched off), , is a partition of , and is an indicator function that takes value 1 if and 0 otherwise. Note that is a nonconvex combinatorial function.
Furthermore, the total transmit power consumption is given by
(4) 
where is the drain inefficiency coefficient of the radio frequency power amplifier [7]. Therefore, the network power consumption is represented by the combinatorial composite function
(5) 
Given the QoS thresholds for all the MUs, in this paper, we aim at solving the following network power consumption minimization problem with the QoS constraints:
subject to  (6) 
Let form the beamforming coefficient group from RRH to all the MUs. Note that, when the RRH is switched off, all the beamforming coefficients in will be set to zeros simultaneously. Observing that there may be multiple RRHs being switched off to minimize the network power consumption, the optimal beamforming vector should have a groupsparsity structure. Therefore, problem is called a group sparse beamforming problem [7]. Note that, to simplify the presentation, we only impose the QoS constraints in problem . However, the proposed instantaneous CSI based iterative reweighted algorithm can be extended to the scenario with perRRH transmit power constraints, following the principles in [20]. More technical efforts are required to derive the asymptotic results using the random matrix method with more complicated structures of the optimal beamformers.
IiC Problem Analysis
Although the constraints in problem can be reformulated as convex secondorder cone constraints, the nonconvex objective function makes it highly intractable. To address this challenge, a weighted mixed norm minimization approach was proposed in [7] to convexify the objective and induce the group sparsity in beamforming vector , thereby guiding the RRH ordering to enable adaptively RRH selection. Specifically, we will first solve the norm minimization problem, and denote as the induced (approximated) group sparse beamforming vectors. Then the following RRH ordering criterion is adopted to determine which RRHs should be switched off [7]:
(7) 
where . In particular, the group sparsity structure information for beamforming vector is extracted from the squared norm of the beamforming vectors ’s, i.e., . The RRH with a smaller will have a higher priority to be switched off. Based on the RRH ordering result in (7), a bisection search approach can be used to find the optimal active RRHs [7] via solving a sequence of the following feasibility problems:
subject to  (8) 
where and . Problem turns out to be convex via reformulating the QoS constraints as secondorder cone constraints [7]. To further enhance the sparsity as well as to seek the quadratic forms of the beamforming vectors in the multicast transmission setting, a smoothed minimization approach was proposed in [15]. To scale to large problem sizes in dense CloudRANs, a twostage largescale parallel convex optimization framework was developed in [16] with the capability of infeasibility detection.
However, all of the above developed algorithms bear high computation overhead. In particular, while the feasibility problem (8) for RRH selection can be efficiently solved with the largescale optimization algorithm in [8], the RRH ordering criterion in (7) may be highly complicated to obtain, especially with the nonconvex formulation as in [15]. Moreover, the ordering criterion needs to be recomputed for each transmission slot, and depends on instantaneous CSI. Observing that the statistical CSI normally changes much slower than the instantaneous CSI, to reduce the computational burden, we propose to compute the RRH ordering criterion (7) only based on statistical CSI, i.e., the longterm channel attenuation. This is achieved by first developing a group sparsity penalty with quadratic forms in the aggregative beamforming vector , followed by an iterative reweighted algorithm with closedform solutions at each iteration via Lagrangian duality theory, as will be presented in Section III. Then asymptotic analysis is performed to obtain the RRH ordering criterion (7) based only on statistical CSI by leveraging the largedimensional random matrix theory [23, 24, 27], as will be presented in Section IV.
Overall, the proposed threestage enhanced group sparse beamforming framework is presented in Fig. 1. Specifically, in the first stage, the RRH ordering criterion is calculated only based on statistical CSI using Algorithm 2. In the second stage, the set of active RRHs is obtained based on instantaneous CSI via solving a sequence of convex feasibility problems using the largescale convex optimization framework in [16]. In the third stage, the transmit power is minimized by solving the convex program (IIB) with the fixed active RRHs using the largescale convex optimization algorithm in [16]. Overall, the proposed enhanced group sparse beamforming framework is scalable to large network sizes. This paper will focus on developing an effective algorithm for the first stage.
Remark 1
In this paper, we assume that, in the first stage, problem (IIB) is feasible for developing Algorithm 1 and Algorithm 2 to determine the RRH ordering criterion based on the instantaneous CSI and statistical CSI, respectively. We enable the capability of handing the infeasibility in the second stage based on the instantaneous CSI, using the largescale convex optimization algorithm [16]. Furthermore, when problem (IIC) is infeasible with all RRHs active, we adopt the user admission algorithm to find the maximum number of admitted users [15].
Iii Group Sparse Beamforming via Iterative Reweighted Algorithm
In this section, we develop a group sparse beamforming approach based on the smoothed minimization, supported by an iterative reweighted algorithm, thereby enhancing the group sparsity in the beamforming vectors. Instead of reformulating the QoS constraints in problem as secondorder cone constraints [7], we use the Lagrangian duality theory to reveal the structures of the optimal solutions at each iteration. The results will assist the large system analysis in Section IV.
Iiia Group Sparsity Inducing Optimization via Smoothed Minimization
To induce the group sparsity structure in the beamforming vector , thereby guiding the RRH ordering, the weighted mixed norm was proposed in [7]. However, the nonsmooth mixed norm fails to introduce the quadratic forms in the beamforming vector , in order to be compatible with the quadratic QoS constraints so that the Lagrangian duality theory can be applied [20]. To address this challenge, we adopt the following smoothed minimization () approach to induce group sparsity [14, 15]:
subject to  (9) 
where is some fixed regularizing parameter and is the weight for the beamforming coefficient group by encoding the prior information of system parameters [7]. Compared with the mixed norm minimization approach, the minimization approach can induce sparser solutions based on the fact . Unfortunately, problem (IIIA) is nonconvex due to the nonconvexity of both the objective function and the QoS constraints.
IiiB Iterative Reweighted Algorithm
We use the MM algorithm and the Lagrangian duality theory to solve problem (IIIA) with closedform solutions. Specifically, this approach generates the iterates by successively minimizing upper bounds of the objective function . We adopt the following upper bounds to approximate the smoothed norm by leveraging the results of the expectationmaximization (EM) algorithm [44].
Proposition 1
Given the value of at the th iteration, an upper bound for the objective function can be constructed as follows:
(10) 
where
(11) 
Proof:
The proof is mostly based on [15, Proposition 1]. \qed
Therefore, at the th iteration, we need to solve the following optimization problem:
subject to  (12) 
which is nonconvex due to the nonconvex QoS constraints. Although the QoS constraints can be reformulated as secondorder cone constraints as in [7], in this paper, we leverage the Lagrangian duality theory to obtain explicit structures of the optimal solution to problem , thereby reducing the computational complexity and further aiding the large system analysis in next section. The success of applying the Lagrangian duality approach is based on the fact that strong duality holds for problem [45], i.e., the gap between the primal optimal objective and the dual optimal objective is zero.
IiiB1 Simple Solution Structures
As the strong duality holds for problem [45], we solve it using duality theory. Specifically, let denote the Lagrange multipliers corresponding to the QoS constraints in problem . Define the Lagrangian function with as
(13)  
where is a block diagonal matrix with the scaled identity matrix as the th main diagonal block square matrix.
To find the optimal ’s, we take the gradient of the Lagrangian function with respect to and set it to zero, which implies
(14) 
By adding to both sides of (14), we have
(15) 
which implies . As is a scalar, the optimal must be parallel to the beamforming direction
(16) 
Therefore, the optimal beamforming vectors can be written as
(17) 
where denotes the optimal beamforming power. Observing that the beamforming powers ’s need to satisfy the SINR constraints with equality [20], i.e.,
(18) 
the optimal powers thus can be obtained via solving the following linear equation:
(19) 
where
Here, denotes the th element of the matrix .
To find the optimal parameter, by multiplying both sides by in (15), we obtain as
(20) 
which implies
(21) 
These fixedpoint equations can be computed using iterative function evaluation.
Based on (17), (19) and (21), we obtain the squared norm of the optimal solution to problem as follows:
(22) 
where is a block diagonal matrix with the identity matrix as the th main diagonal block square matrix and zeros elsewhere. Note that the optimal powers ’s (19) and Lagrangian multipliers ’s (21) should depend on the weights ’s (11) at each iteration .
The instantaneous CSI based iterative reweighted algorithm for group sparse beamforming is presented in Algorithm 1. We have the following result for its performance and convergence.
Theorem 1
Let be the sequence generated by the iterative reweighted algorithm. Then, every limit point of has the following properties:

is a KKT point of problem (IIIA);

converges monotonically to for some KKT point .
Proof:
Please refer to Appendix A for details. \qed
Compared with the algorithm for the smoothed minimization in [15], the main novelty of the proposed iterative reweighted algorithm is revealing the explicit structures of the solutions at each iteration in Algorithm 1. Instead of using the interiorpoint algorithm to solve the convex subproblems at each iteration [46], Algorithm 1 with closedform solutions helps reduce the computational cost. The proposed iterative reweighted algorithm can only guarantee to converge to a KKT point, which may be a local minimum or the other stationary point (e.g., a saddle point and local maximum).
Remark 2
The main contributions of the developed Algorithm 1 include finding the closedform solutions in the iterations, in comparison with the work [15] using the interiorpoint algorithm, and proving the convergence of the lowcomplexity closedform iterative algorithm for the nonconvex group sparse inducing problem (IIIA), instead of the convex coordinated beamforming problem [45].
Unfortunately, computing the RRH ordering criterion ’s in (7) requires to run Algorithm 1 for each channel realization, which brings a heavy computation burden. In next section, we will find deterministic approximations for ’s to determine the RRH ordering only based on statistical CSI, which changes much more slowly than instantaneous channel states, and thus requires less frequent update.
Iv Group Sparse Beamforming via Large System Analysis
In this section, we present the large system analysis for the iterative reweighted algorithm in Algorithm 1, thereby enabling RRH ordering only based on statistical CSI. In this way, the ordering criterion will change only when the longterm channel attenuation is updated, and thus it can further reduce the computational complexity of group sparse beamforming. The main novelty of this section is providing the closedforms for the asymptotic analysis of the optimal Lagrangian multipliers based on the recent results in [27], thereby providing explicit expressions for the asymptotic RRH ordering results.
Iva Deterministic Equivalent of Optimal Parameters
In this subsection, we provide asymptotic analysis for the optimal beamforming parameters ’s (19) and ’s (21) when . In a CloudRAN with distributed RRHs, the channel can be modeled as , where is the smallscale fading and with as the pathloss from RRH to MU . With this channel model, we have the following result for the deterministic equivalent of the parameter.
Lemma 1 (Asymptotic Results for Parameter)
Assume
Let and satisfy and , respectively. We have almost surely, where is given by
(23) 
Here, is the unique solution of the following set of equations:
(24) 
Proof:
Please refer to Appendix B for details. \qed
Based on the above results, we further have the following asymptotic result for the optimal powers ’s.
Lemma 2 (Asymptotic Results for Optimal Powers)
Let be such that
(25) 
If and only if , then almost surely, where is given by
(26) 
Here , and are given as follows:
(27) 
and
(28) 
and
(29) 
respectively; and is given as
(30) 
where with
(31) 
Proof:
Please refer to Appendix C for details. \qed
IvB Statistical CSI Based Group Sparse Beamforming Algorithm
Based on the deterministic equivalents of the optimal beamforming parameters ’s and ’s in Lemma 1 and Lemma 2, respectively, we have the following theorem on the squared norm of the solution to problem , i.e., .
Theorem 2 (Asymptotic Results for the Beamformers)
At the th iteration, for the squared norm of the solution to problem , , we have almost surely with
(32) 
where .
Proof:
Please refer to Appendix D for details. \qed
We thus have the following deterministic equivalent of the weights ’s (11):
(33) 
Note that the asymptotic results of the powers ’s (26) and Lagrangian multipliers ’s (23) should depend on the weights ’s (33) at each iteration. Based on (7), we have the following asymptotic result for the RRH ordering criterion:
(34) 
where is the deterministic equivalent of the squared norm of the solution to problem using the iterative reweighted algorithm. Therefore, the RRH with a smaller will have a higher priority to be switched off. This ordering criterion will change only when the longterm channel attenuation is updated. Note that global instantaneous CSI is still needed in stage II to find the active RRHs, i.e., solving a sequence of convex feasibility problems (IIC).
The statistical CSI based iterative reweighted algorithm for group sparse beamforming is presented in Algorithm 2.
Remark 3
Although stage 2 and stage 3 still require instantaneous CSI, the RRH ordering criterion computed by Algorithm 2 in stage 1 is only based on statistical CSI and thus will be updated only when the longterm channel propagation is changed. Therefore, Algorithm 2 serves the purpose of further reducing the computation complexity of Algorithm 1 for RRH ordering. A promising future research direction is to select RRHs in stage 2 only based on statistical CSI. However, the main challenge is the infeasibility issue if only statistical CSI is available, as discussed in Remark 1.
V Simulation Results
In this section, we will simulate the proposed algorithms for network power minimization for CloudRANs. In all the realizations, we only account the channel realizations making the original problem feasible. If problem is infeasible, further works on user admission are required [15]. We set , and for all the simulations. The proposed iterative reweighted algorithms (Algorithm 1 and Algorithm 2) will terminate if either the number of iterations exceeds 30 or the difference between the objective values of consecutive iterations is less than .
Va Convergence of the Iteratively Reweighted Algorithm
Consider a network with 30antenna RRHs and 5 single antenna MUs uniformly and independently distributed in the square region meters. Fig. 2 shows the convergence of the iterative reweighted algorithm in different scenarios with different realizations of RRH and MU positions. For each scenario, the numerical results are obtained by averaging over 100 smallscale fading realizations. The largescale fading (i.e., statistical CSI) is randomly generated and fixed during simulations. In the two curves of each scenario, Algorithm 1 and Algorithm 2 are applied for Stage I of group sparse beamforming, respectively. It demonstrates that the large system analysis based Algorithm 2 provides accurate approximation even in a small system. Fast convergence is observed in the simulated setting.
VB Network Power Minimization

0  3  6  9  12  

GSBF [7]  27.671  28.528  29.924  32.759  43.855  

23.716  24.374  26.011  30.057  43.810  

23.789  24.463  26.092  30.500  43.783  

23.907  24.695  26.326  30.214  43.827  

23.905  24.689  26.315  30.205  43.821 
Consider a network with 10antenna RRHs and 6 single antenna MUs uniformly and independently distributed in the square region meters. The relative transport link power consumption are set to be . We average over 200 smallscale channel realizations. Fig. 3 demonstrates the average network power consumption using different algorithms. This figure shows that the proposed iterative reweighted algorithm outperforms the weighted mixed norm algorithm in [7]. Furthermore, it illustrates that the large system analysis provides accurate approximations for the iterative reweighted algorithm in finite systems with reduced computation overhead. In Table I, for the simulated scenarios, we see that the iterative reweighted algorithms achieve similar performance with different values of parameter and the large system analysis provides accurate approximations. Although the simulated results demonstrate that the performance is robust to different values of parameter , it is very interesting to theoretically identify the typical scenarios, where smaller values of parameter will yield much lower network power consumption.
Vi Conclusions and Discussions
In this paper, we developed an enhanced threestage group sparse beamforming framework with reduced computation overhead in CloudRANs. In particular, closedform solutions for the group sparse optimization problem were obtained by developing the iterative reweighted algorithm based on the MM algorithm and Lagrangian duality theory. This is the first effort to reduce the computation cost for RRH ordering in the first stage of group sparse beamforming. Based on the developed structured iterative algorithm, we further provided a large system analysis for the optimal Lagrangian multipliers at each iteration via random matrix theory, thereby computing the RRH ordering criterion only based on the statistical CSI. This is the second effort to enable computation scalability for RRH ordering in the first stage.
Several future directions of interest are listed as follows:

Although computing the RRH ordering criterion in the first stage with Algorithm 2 only needs statistical CSI, the second stage still requires global instantaneous CSI to find the active RRHs by solving a sequence of convex feasibility problems. It is thus particularly interesting to investigate efficient algorithms to select the active RRHs only based on the statistical CSI.

It is desirable to establish the optimality of the iterative reweighted algorithm for the network power minimization problem . However, considering the complicated problem structures of problem , this becomes challenging. It is also interesting to apply the developed statistical CSI based iterative reweighted algorithm for more complicated network optimization problems, e.g., the wireless caching problem [33, 34], the computation offloading problem [36], and beamforming problems with CSI uncertainty [32].
Appendix A Proof of Theorem 1
1) Observing that the phases of will not change the objective function and constraints of problem , the QoS constraints thus can be equivalently transformed to the secondorder cone constraints: , which are convex cones. The KKT points of problem should satisfy:
(35) 
where is the normal cone of the secondorder cone at point [47]. We shall show that any convergent subsequence of satisfies (35). Specifically, let be one such convergent subsequence with
(36) 
Based on the strong duality result for problem , the following KKT condition holds at :
(37) 
Based on (36) and (11), we have . Furthermore, based on [47, Proposition 6.6] and (36), we have . Therefore, by taking in (37), we have , which indicates that is a KKT point of problem . We thus complete the proof.
2) We first conclude the following fact:
(38)  
where the first inequality is based on the fact that function attains its maximum at [15, Proposition 1], and the second inequality follows from (37). Furthermore, as is continuous and is compact, the limit of the sequence is finite. Based on the convergence results in 1), we thus complete the proof.