1 Introduction
Cellfree Massive MIMO (multiple input multiple output) is a new promising wireless technology to deal with the mediocre celledge performance of cellular networks [1]. The system consists of many distributed access points (APs) that transmit coherently in the downlink and process the received signals coherently in the uplink [2, 3]
, leading to a higher signaltonoise ratio (SNR) without using more power. The distributed nature gives massive macrodiversity against pathloss and shadow fading. The spectral efficiency (SE) of Cellfree Massive MIMO has been characterized for Rayleigh fading channels and various implementations in
[4, 2, 3, 5]. The SE was optimized in [2, 3, 6, 7] (among others) to achieve maxmin fairness.In this paper, we instead consider that each user has an SE requirement that the system must satisfy to not cause service interruptions. Hence, the goal of the resource allocation is for the system to deliver the required SEs with as low power consumption as possible, leading to maximum energy efficiency. We stress that this approach to energy efficiency maximization is different from [6, 8]. In particular, we will consider the possibility to turn off APs that are not needed to serve the current set of users, bearing in mind that each user will mainly be served by its neighboring APs. This is an important feature since Cellfree Massive MIMO systems may have many APs [2, 3], where the large number is needed to provide consistent coverage but might not be needed at every time instant.
1.1 Relation to Prior Work
The problem of AP activation in distributed antenna systems has previously been analyzed for cloud radio access networks [9, 10, 11]. However, these prior works consider slowly fading channels and perfect channel state information (CSI). In contrast, we consider fast fading channels, imperfect CSI, and pilot contamination. Hence, the optimization problems considered in this paper are entirely different from those in previous work, and more useful for implementation.
1.2 Contributions
In this paper, we consider a Cellfree Massive MIMO system with multiantenna APs and maximum ratio transmission (MRT) in the downlink. We formulate a total power minimization problem where the active APs and transmit power allocation are the variables. This problem is nonconvex but can be solved as a mixedinteger secondorder cone (SOC) program. Since this approach is too complex for realtime applications, we develop a lowercomplexity solution by exploiting the inherent sparsity in the problem. The numerical results demonstrate that there are scenarios where only a subset of the APs are needed to satisfy the SE requirements and large power reductions can be achieved by turning off the remaining APs.
Notations:
We use boldface lowercase and uppercase letters to denote vectors and matrices, respectively. The transpose is denoted by the superscript
and the Hermitian transpose is denoted by . The expectation operator is anddenotes a circularly symmetric complex Gaussian distribution. The Euclidean norm,
norm, and norm of a vector is denoted as , , and , respectively. The cardinality of a set is denoted .2 System Model
We consider a Cellfree Massive MIMO system with distributed APs and singleantenna users. Each AP has antennas and is connected to a central processing unit (CPU) via unlimited fronthaul links. We consider the standard block fading model [12], where the channels are fixed and frequencyflat within a coherence interval of channel uses. The channel between AP and user is modeled by uncorrelated Rayleigh fading as
(1) 
where is the largescale fading coefficient.
The users are served by coherent joint transmission from the APs. We assume that each user has a required downlink SE value [b/s/Hz] that must be satisfied. The users and APs are arbitrarily distributed, thus it is likely that these SE requirements can be fulfilled without using all the APs. The main goal of this paper is to find a subset of active APs and the corresponding transmit powers that satisfy the SE requirements while minimizing the total power consumption.
We consider a timedivision duplex (TDD) protocol where each AP estimates the channels from itself to the users by using uplink pilot transmissions. As in [2, 3], orthonormal pilot signals are utilized and assigned to the users. We let denote the subset of users assigned to the same pilot as user . By following standard methods for channel estimation and downlink precoding [13, 14], we obtain the following result.
Lemma 1.
Suppose MRT is used for downlink transmission and AP allocates power to user , then an ergodic SE of user is
(2) 
where the effective SINR is
(3) 
the noise variance is
, is the uplink pilot power of user , and the meansquare of the channel estimates between AP and user is(4) 
The difference between Lemma 1 and the previous work [13, 14] is that only the subset of the APs are active. In (3), the numerator is proportional to which is the array gain achieved by having multiple antennas at each AP. The first term in the denominator is coherent interference from the pilotsharing users. The second term is conventional noncoherent interference and the third term is noise.
3 Total Power Minimization Problem
In this section, we formulate a total power minimization problem under the SE requirements of the users and a maximum transmit power per AP. We then rewrite it as a mixedinteger SOC program, which can be solved but with high complexity. Similar to [6], we model the total power consumption of the all active APs as
(5) 
where the first term is the transmit power multiplied with the inefficiency factor of the power amplifiers and the second term contains the fixed power per active AP.
The total power consumption minimization problem is
(6a)  
(6b)  
(6c) 
This is a complicated problem since all APs can transmit data to all users and also cause interference to all users. We first define and notice that the SE constraint in (6b) can be rewritten as
(7) 
To simplify the problem, we introduce the notation , , , , . These vectors and matrices are defined as
(8)  
(9)  
(10)  
(11)  
(12) 
where are the members of the set (i.e., the indices of the active APs). The th column of in (11) is denoted and the th row is denoted as . In (12), are the indices of the users belonging to the set , and is the cardinality of the set . The Hadamard product is denoted by . By using these notations and (7), we can obtain an equivalent epigraph representation of problem (6) as
(13a)  
(13b)  
(13c)  
(13d) 
The auxiliary variable moves the objective function of problem (6) to the first constraint in (13b). We observe that for a given , problem (6) reduces to an SOC program, as previously shown in [2, 3]. Hence, although (13) is nonconvex and NPhard, it can be solved by making an exhaustive search over all possible selections of . Since at least one AP needs to be active if there is users with nonzero SE requirements, there are different selections of the APs that need to be considered in an exhaustive search.
3.1 Globally Optimal Solution
Making an exhaustive search to solve (13) is very computationally costly even in a relatively small network. We will therefore further reformulate the problem to reduce the complexity, while guaranteeing to find the global optimum.
Let the binary optimization variable characterize the on/off activity of AP . We can then replace the maximum transmit power of AP by , which takes the original value when the AP is active and is zero when the AP is turned off. This feature is exploited to formulate a mixedinteger SOC program.
Lemma 2.
Consider the mixedinteger SOC program
(14a)  
(14b)  
(14c)  
(14d)  
(14e) 
where is the th row of matrix and . Moreover, the vectors and are defined as
(15)  
(16) 
with and .
Problems (13) and (14) are equivalent in the sense that they have the same optimal transmit powers. If we denote by
an optimal solution to the binary variables
in (14), the optimal set of active APs in problem (13) is(17) 
Proof.
The proof is omitted due to the space limitations. ∎
Problem (14) is a mixedinteger SOC program on standard form and can, thus, be solved by the generalpurpose toolbox CVX [15] using the MOSEK solver [16]. These softwares apply a branchandbound approach to deal with the binary variables.
The new binary variables provide the explicit link between the hardware and transmit power consumption, which is an important factor to obtain the global optimum to problem (14). A key reason that we can preserve the SOC structure despite adding the new binary variables is that the binary variables are not involved in the SINR constraints (14c). Instead there is an implicit connection via the zero maximum transmit power for inactive APs. This is different from previous approaches (e.g., [9]) where appears in the SINRs and therefore would break the SOC structure.
4 SparsityBased LowComplexity Algorithm
Motivated by the high computational complexity of solving the total power minimization problem using Lemma 2, we will now propose an algorithm that finds a good suboptimal solution with tolerable complexity. Since we are searching for a solution where many of the power variables are zero, we will reexpress (6) as a sparse reconstruction problem where we try to push many of the transmit power variables to become zero. We begin with the following result.
Lemma 3.
Proof.
The proof is omitted due to the space limitations. ∎
Lemma 3 shows that we do not need separate variables for optimizing the active APs set , but we can implicitly determine if AP is active or not by checking if or . The reformulated problem (18) has fewer variables but remains nonconvex due to the norm in the objective. A standard approach is to relax the norm to a convex norm, but we will take another approach that is known to utilize the sparsity more effectively [17]. We consider an norm relaxation of problem (18) for some :^{1}^{1}1From the range of the considered norms, the condition as in (21) leads to .
(21a)  
(21b)  
(21c) 
The objective function of problem (21) treats every vector as an entity in when seeking for a sparse solution. This groupsparse approach differs from the previous works that consider elementbased [18] or beamformingvectorbased [19] sparsity.
Even though problem (21) remains nonconvex after the norm relaxation, we can find a stationary point by adapting the iteratively reweighted least square approach [20], that was originally developed for componentwise sparsity. By dropping the exponent and in (21a), we achieve the following result.
Theorem 1.
Suppose problem (21) is feasible. Since the feasible set is convex, we can construct an iterative algorithm that starts with the given initial weight values and in iteration solves the following SOC program:
(22)  
to yield the solution , for which
(23) 
is the optimal transmit powers for AP at iteration . After that, the weight values are updated for the next iteration as
(24) 
where is a sufficient small positive damping constant with and . The proposed iterative process exhibits the properties below:

If an arbitrary AP has zero transmit power at the optimum of iteration , this AP will have zero transmit power in all the following iterations.
Proof.
The proof is omitted due to the space limitations. ∎
Theorem 1 guarantees a monotonically decreasing objective function and the main computational cost is to solve (22) in each iteration. The iterative process reaches a stationary point to problem (21) where APs are turned off along the iterations. The damping constant is introduced to cope with numerical issues that can appear when updating the weight values (24). The stopping criterion can be selected by comparing two consecutive iterations.
Due to the normrelaxation, we will not use the solution from Theorem 1 as the final solution but instead as an indication of which APs to turn off. More precisely, we compute the transmit power that the APs utilize at the solution from Theorem 1 and reorder the APs in increasing power order. We then make a bisectionlike search over how many APs should be turned on and always keep the ones that utilize the most power at the solution from Theorem 1. We thereby identify an AP set that minimizes the total power consumption.
5 Numerical Results
We have simulated a setup where APs and users are randomly distributed within a squared area of km
by a uniform distribution, given that the distance between two APs should be larger than
m. Each AP is equipped with antennas. The requested SE from each user is [b/s/Hz]. We apply the wraparound structure to get rid of edge effects and guarantee uniform simulation performance for the APs. Coherence intervals have symbols. There are orthogonal pilots and each is assigned to randomly selected users. The pilots are transmitted with uplink power W. We use the largescale fading model with correlated shadow fading from [5], which matches well with the 3GPP Urban Microcell model for a carrier frequency GHz. The maximum downlink power per AP is W and the noise variance is dBm. The power consumption is modeled similar to [6]: and W.We compare the following methods:
Fig. 1
shows the cumulative distribution function (CDF) of the total transmit power [W] achieved by the four different methods, where the randomness is due to the different AP and user locations. When all APs are active, the transmit power is
W on the average. The mixedinteger SOC program uses roughly more transmit power: W on the average. The proposed sparsitybased method utilizes about W, while the disjoint sparsity benchmark uses the highest transmit power level: about W.The proposed methods are not minimizing the transmit power but the total transmit power. Fig. 2 shows the CDF of the total power consumption [W]. Contrary to the previous figure, the benchmark where all APs are active has the the highest total power consumption: about W on average. By solving the proposed mixedinteger SOC program, we find the globally minimum total power consumption, which saves about compared with the baseline. The proposed sparsitybased method requires around more power than the global minimum. In contrast, the previous disjoint sparsity benchmark requires more power than the global minimum.
6 Conclusion
This paper has minimized the total downlink power consumption in cellfree Massive MIMO systems by jointly optimizing the number of active APs and their transmit powers, while satisfying the SEs requested by the users. A globally optimal solution can be found by formulating the considered problem as a mixedinteger SOC program. We observe a considerable reduction in the total power consumption (about ) compared with only minimizing the total transmit powers as in previous work. Since the mixedinteger SOC program has high computational complexity, we also developed a lowercomplexity algorithm that finds a good suboptimal solution with relatively low complexity by utilizing sparsity methods. In the simulation part, this method only requires 20% more power than the global minimum.
References
 [1] Giovanni Interdonato, Emil Björnson, Hien Quoc Ngo, Pål Frenger, and Erik G Larsson, “Ubiquitous cellfree massive MIMO communications,” EURASIP J. Wireless Commun..Netw., vol. 2019, no. 1, pp. 197, 2019.
 [2] Hien Quoc Ngo, Alexei Ashikhmin, Hong Yang, Erik G Larsson, and Thomas L Marzetta, “Cellfree massive MIMO versus small cells,” IEEE Trans. Wireless Commun., vol. 16, no. 3, pp. 1834–1850, Mar. 2017.
 [3] Elina Nayebi, Alexei Ashikhmin, Thomas L. Marzetta, Hong Yang, and Bhaskar D. Rao, “Precoding and power optimization in cellfree Massive MIMO systems,” IEEE Trans. Wireless Commun., vol. 16, no. 7, 2017.
 [4] Elina Nayebi, Alexei Ashikhmin, Thomas L Marzetta, and Bhaskar D Rao, “Performance of cellfree massive MIMO systems with MMSE and LSFD receivers,” in Asilomar Conference on Signals, Systems and Computers, 2016.
 [5] Emil Björnson and Luca Sanguinetti, “Making cellfree Massive MIMO competitive with MMSE processing and centralized implementation,” IEEE Trans. Wireless Commun., 2019, to appear. Available: https://arxiv.org/abs/1903.10611.
 [6] Hien Quoc Ngo, LeNam Tran, Trung Q Duong, Michail Matthaiou, and Erik G Larsson, “On the total energy efficiency of cellfree Massive MIMO,” IEEE Transactions on Green Communications and Networking, vol. 2, no. 1, pp. 25–39, 2018.
 [7] Manijeh Bashar, Kanapathippillai Cumanan, Alister G Burr, Mérouane Debbah, and Hien Quoc Ngo, “Enhanced maxmin SINR for uplink cellfree Massive MIMO systems,” arXiv preprint arXiv:1801.10188, 2018.
 [8] Long D Nguyen, Trung Q Duong, Hien Quoc Ngo, and Kamel Tourki, “Energy efficiency in cellfree Massive MIMO with zeroforcing precoding design,” IEEE Commun. Lett., vol. 21, no. 8, pp. 1871–1874, 2017.
 [9] Mingjie Feng, Shiwen Mao, and Tao Jiang, “BOOST: Base station onoff switching strategy for green massive MIMO HetNets,” IEEE Trans. Wireless Commun., vol. 16, no. 11, pp. 7319–7332, 2017.
 [10] QuangDoanh Vu, LeNam Tran, Markku Juntti, and EenKee Hong, “Energyefficient bandwidth and power allocation for multihoming networks,” IEEE Trans. Signal Process., vol. 63, no. 7, pp. 1684–1699, 2015.
 [11] Fengxia Han, Shengjie Zhao, Lu Zhang, and Jinsong Wu, “Survey of strategies for switching off base stations in heterogeneous networks for greener 5G systems,” IEEE Access, vol. 4, pp. 4959–4973, 2016.
 [12] Emil Björnson, Jakob Hoydis, and Luca Sanguinetti, “Massive MIMO networks: Spectral, energy, and hardware efficiency,” Foundations and Trends® in Signal Processing, vol. 11, no. 34, pp. 154–655, 2017.
 [13] Trinh Van Chien, Emil Björnson, and Erik G. Larsson, “Joint power allocation and user association optimization for Massive MIMO systems,” IEEE Trans. Wireless Commun., vol. 15, no. 9, pp. 6384 – 6399, 2016.
 [14] Giovanni Interdonato, Marcus Karlsson, Emil Björnson, and Erik G Larsson, “Downlink spectral efficiency of cellfree Massive MIMO with fullpilot zeroforcing,” in Proc. IEEE GlobalSIP, 2018.
 [15] CVX Research Inc., “CVX: Matlab software for disciplined convex programming, academic users,” http://cvxr.com/cvx, 2015.
 [16] Erling D. Andersen, Bo Jensen, Jens Jensen, Rune Sandvik, and Ulf Worsøe, “MOSEK version 6,” 2009, Technical Report TRC2009C3, MOSEK.
 [17] Zheng Zhang, Yong Xu, Jian Yang, Xuelong Li, and David Zhang, “A survey of sparse representation: algorithms and applications,” IEEE Access, vol. 3, pp. 490–530, 2015.
 [18] Emmanuel J Candes, Michael B Wakin, and Stephen P Boyd, “Enhancing sparsity by reweighted minimization,” Journal of Fourier analysis and applications, vol. 14, no. 56, pp. 877–905, 2008.
 [19] Shixin Luo, Rui Zhang, and Teng Joon Lim, “Downlink and uplink energy minimization through user association and beamforming in CRAN,” IEEE Trans. Wireless Commun., vol. 14, no. 1, pp. 494–508, 2014.
 [20] Demba Ba, Behtash Babadi, Patrick L Purdon, and Emery N Brown, “Convergence and stability of iteratively reweighted least squares algorithms,” IEEE Trans. Signal Process., vol. 62, no. 1, pp. 183–195, 2014.
 [21] Yuanming Shi, Jinkun Cheng, Jun Zhang, Bo Bai, Wei Chen, and Khaled B Letaief, “Smoothed minimization for green cloudRAN with user admission control,” IEEE J. Sel. Areas Commun., vol. 34, no. 4, pp. 1022–1036, 2016.
Comments
There are no comments yet.