A Low Complexity Detection Algorithm For Uplink Massive MIMO Systems Based on Alternating Minimization
Abstract
In this paper, we propose an algorithm based on the Alternating Minimization technique to solve the uplink massive MIMO detection problem. The proposed algorithm provides a lower complexity compared to the conventional MMSE detection technique, especially when the number of user equipment (UE) antennas is close to the number of base station (BS) antennas. This improvement is obtained without any matrix inversion. Moreover, the algorithm reformulates the maximum likelihood (ML) detection problem as a sum of convex functions based on decomposing the received vector into multiple vectors. Each vector represents the contribution of one of the transmitted symbols in the received vector. Alternating Minimization is used to solve the new formulated problem in an iterative manner with a closed form solution update in every iteration. Simulation results demonstrate the efficacy of the proposed algorithm in the uplink massive MIMO setting for both coded and uncoded cases.
I Introduction
Massive multipleinput multipleoutput (MIMO) is one of the most promising techniques for the 5th Generation (5G) networks due to its potential for enhancing throughput, spectra efficiency, and energy efficiency [1], [2]. Massive MIMO requires the BS to be equipped with arrays of hundreds of antennas to serve tens of user terminals with single or multi antennas.
Theoretical results of massive MIMO show that linear detectors such as zero forcing (ZF) and minimum mean square error (MMSE) can achieve optimum performance under the favorable propagation conditions [3]. The favorable propagation conditions mean that the number of BS antennas grows very large compared to the number of UE antennas, which leads to the columnvectors of the propagation matrix to be asymptotically orthogonal. As a result of this orthogonality, the ZF and MMSE detectors can be implemented with simple diagonal inversions [4].
The current practical number of the BS antennas in massive MIMO systems is in the order of tens to a hundred. This is far from the theoretical limit that leads to the orthogonality mentioned above [4]. Therefore, the linear detectors still need to perform a matrix inversion for the signal detection of the uplink massive MIMO system, which entails extensive computational complexity [5]. Neumann series expansion, Cholesky decomposition, and successive overrelaxation techniques are proposed in the literature to reduce the complexity of the matrix inversion process in the MMSE detector [6] [7]. These approaches require a lower computational complexity than the exact matrix inversion while delivering nearoptimal results only for theoretical massive MIMO configurations. In realistic massive MIMO scenarios the performance of implementing these approximations is far from the exact MMSE performance [6]. In this letter we present a novel formulation of a low complexity iterative algorithm based on Alternating Minimization, referred to as AltMin. This algorithm provides similar bit error rate (BER) performance to the exact matrix inverted MMSE technique, with one order less complexity.
The proposed algorithm reformulates ML detection problem as a sum of convex functions based on decomposing the received vector into multiple vectors. Each vector represents the contribution of one of the transmitted symbols in the received vector. Then, Alternating Minimization is used to solve the new formulated problem in an iterative manner with a closed form solution update in every iteration that does not require any matrix inversion.
Ii System Model and Problem Formulation
Iia System Model
Consider the uplink data detection in a multiuser (MU) massive MIMO system with BS antennas and UE antennas. The vector represents the complex transmitted signal, where is the transmitted symbol for user with . Each user transmits symbols over a flat fading channels and the signals are demodulated and sampled at the receiver. The vector represents the complex received signal, and the channel matrix can be represented as , where , and is the complex flat fading channel gain from transmit antenna to the receive antenna , with . The system can be modeled as; ,
where the is the complex additive white Gaussian noise (AWGN) vector whose elements are mutually independent and zero mean with variance . The corresponding realvalued system model is [8]. The equivalent ML detection problem of the real model can be written in the form , where , M is the QAM constellation size, and the normalization factor.
IiB Problem Formulation
First, we decompose the received vector y into a linear combination of vectors so that , where represents the contribution of the th transmitted symbol in the received vector. The element wise representation of the decomposed received vector is:
The th element of y can be represented as , . Let the column of the channel matrix H. Now, we relax the nonconvexity constraint on the feasible set , and reformulate the ML problem based on the above decomposition as follows:
(1) 
subject to
(2) 
(3) 
Where . The objective function in (1) is a sum of separable terms, each of which is a function of only one symbol and its contribution in the received vector. In the next section, we use AltMin to solve the proposed formulation.
Iii Proposed Algorithm
The optimization problem (1)(3) is strictly and jointly convex with respect to x and . Moreover, there is no common constraint that combines both x and . Therefore, in order to efficiently solve this problem, we first decompose it into the following two subproblems:

Given x, we obtain by solving
(4) 
Given , we obtain x by solving
(5)
Then, we propose, AltMin, an iterative algorithm that alternatively solves (4) for and (5) for x. Note, the respective KKT conditions of the above two subproblems form the complete set of the KKT conditions for the original problem.
Iiia Solving the subproblem (4)
The Lagrangian function of (4) can be written as:
(6) 
Therefore, by solving the above Lagrangian function, which is a function of and , we get the following closed form expression updates for every element of and :
(7) 
(8) 
Note, in the update of , we introduce the scaling factor, . We will show in Section IIIC that the proposed algorithm is optimal when . However, we empirically find that when , the number of iterations for convergence drops significantly, as depicted in Section IVA.
IiiB Solving the subproblem (5)
Since the objective function (5) is separable with respect to every element in the vector x and no constraint combines the elements of x. The update of the th element in the vector x () reduces to solving the following subproblem.
(9) 
The corresponding Lagrangian function of (9) is:
(10) 
Then, the following K.K.T. conditions, which are sufficient and necessary for the optimal solution to the convex optimization problem in (9) are:
(11) 
(12) 
To solve (11)(12) for every element in the vector x, among the following choices, we choose the one that minimizes (9):
(13) 
(14) 
(15) 
Note, we excluded the choice , and since cannot be equal to and at the same time.
To obtain the optimal solution to the proposed optimization problem (1)(3), AltMin solves (4) for and (5) for x. To perform the algorithm, we initially set x to 0, and solve (4) to obtain the initial ; with the updated , we then solve (5) to update x. In the next subsection, we will show that Algorithm 1 converges and the pairwise optimal x and can be obtained, which is also the optimal solution to the problem in (1)(3).
IiiC Optimality of AltMin Algorithm in solving (1)(3)
In order to show the optimality of the iterative algorithm (when ), we first give the following lemma.
Lemma 1
Given x, the optimal for the problem (4) is unique. Similarly, given the , the corresponding optimal x is unique.
Proof 1
Given a pair (, x), if is the optimal solution to (4) given x, and x is the optimal solution to (5) given , it can be said that and x are pairwise optimal for (1)(3). Furthermore, for each subproblem, the K.K.T. conditions form a subset of the K.K.T conditions of (1)(3) given the other primal variables, where the two subsets contain no common dual variables. Therefore, if the primal variables are pairwise optimal, the K.K.T. conditions in each corresponding subset are satisfied, all K.K.T. conditions of (1)(3) are satisfied, i.e., the pairwise optimal solution is also the optimal , x for (1)(3).
Lemma 2
Proof 2
Since , and x are updated by successively solving (4) and (5) in every iteration, the obtained objective value is nonincreasing over iterations [9]. At the convergence point, , and x are jointly optimal to the subproblems, since otherwise the objective value can be further decreased in the next iteration. Then, by Lemma 2, the obtained (, x) is the optimal solution to (1)(3).
IiiD Complexity Analysis of AltMin
The update of has a complexity of , while the update of x has a complexity of . The AltMin Algorithm performs iterations between these updates before converging to the optimal solution. Therefore, the overall complexity of AltMin in solving (1)(3) is in the order of , which is notably lower than of MMSE for large , as it is demonstrated in Section IV.
Iv Numerical Results
Simulation results for coded and uncoded uplink massive MIMO systems in a block flat fading channel is presented. We assume perfect knowledge of the channel state information at the receiver. QPSK modulations is considered for demonstration; however, the proposed algorithm can be extended for higher QAM modulations in a straightforward manner. The performance and computational complexity of the proposed algorithm based on AltMin Algorithm is compared with the linear MMSE detector.
Iva Number of AltMin Iterations
In this simulation experiment we examine the number of iterations required by the AltMin algorithm such that the BER performance of both the proposed algorithm and the MMSE technique are equal, for various massive MIMO configurations. We fix convergence tolerance at , while the maximum number of iterations is changing in a step size of 2. We also set the initial guess of x in the AltMin algorithm to zeros, and in the update of in (7).
Fig. 3 shows BER performance versus maximum number of iterations for , and configurations, at SNR=12 dB. It can be noticed that the larger the ratio is, the smaller the number of iterations is required by AltMin to reach MMSE performance. For example, 8 iterations is required for configuration, while 15 iterations is required for configuration. This results show that, on the average, to reach MMSE performance the number of iterations, , is much smaller than , which keeps the complexity order of the proposed algorithm around .
IvB BER performance comparison
In this subsection we present the BER performance of both the proposed algorithm and the MMSE technique with respect to SNR at various massive MIMO configuration. For each SNR value, we stop the AltMin algorithm at the iteration number at which its performance matches the MMSE performance based on the results from Fig. 3.
Fig. 3 shows that the BER of the proposed algorithm is upper bounded by that of the MMSE with the exact matrix inversion. For example, at higher ratio between and , such as or , the performance of the proposed algorithm and MMSE are the same. As the ratio becomes closer to 1, the BER of the proposed algorithm becomes slightly lower than that of MMSE.
IvC Turbo Coded BER Performance
The turbo coded BER performance of the proposed algorithm compared to MMSE is shown in Fig. 3 using coded QPSK modulation. In this simulation, all the above massive MIMO configurations are examined with rate1/2 turbo encoder and decoder of 10 iterations. output valued vector from all detectors is fed as an input to the BCJRbased turbo decoder. In Fig. 3, AltMin based algorithm performs similar to the MMSE detector for , and slightly better than MMSE for . As the number of uplink antennas increases, the coded AltMin based algorithm outperforms coded MMSE, for example the improvement for the case of at coded BER is about 1 dB compared to only 0.2 dB improvement in the case of .
IvD Computational Complexity Analysis
In this subsection we compare the computational complexity of the proposed algorithm with the MMSE technique in terms of the number of multiplication operations, as depicted in Table I. The comparison is based on the same SNR of dB and the same BER performance. More specifically, the number of iterations of AltMin is taken based on the results of Fig. 3, at which the BER of the two techniques coincides.
From Table I, it can be observed that at small , such as 16, MMSE outperforms the proposed algorithm by approximately a factor of 3. While for a large , such as 128, the proposed algorithm shows superior computational reduction by a factor of as compared to MMSE. Note, although our algorithm requires more computations than MMSE for small , it does not exhibit any matrix inversion, which is more advantageous in terms of hardware implementations [5].
=16  =32  =64  =128  
MMSE  0.057  0.311  2.195  16.97 
AltMin  0.204  0.409  1.409  2.818 
3.57  1.31  0.64  0.166 
V Conclusion
In this letter, we proposed an iterative low complexity algorithm based on Alternating Minimization approach. This algorithm is a good alternative for the MMSE technique in Massive MIMO applications where the number of BS antennas is not very large compared to the number of user equipment antennas. We show that the proposed algorithm avoids complicated matrix inversion by solving the reformulated ML problem in an iterative manner, in which each iterations performs a simple computations based on a closed form solution. The results reveal that the algorithm can provide lower computational complexity as compared to the exact matrix inverted MMSE with similar BER performance in both coded and uncoded cases.
References
 [1] F. Boccardi, R. W. Heath, A. Lozano, T. L. Marzetta, and P. Popovski, “Five disruptive technology directions for 5g,” IEEE Communications Magazine, vol. 52, no. 2, pp. 74–80, 2014.
 [2] T. Van Chien and E. Björnson, “Massive mimo communications,” in 5G Mobile Communications. Springer, 2017, pp. 77–116.
 [3] L. Lu, G. Y. Li, A. L. Swindlehurst, A. Ashikhmin, and R. Zhang, “An overview of massive mimo: Benefits and challenges,” IEEE journal of selected topics in signal processing, vol. 8, no. 5, pp. 742–758, 2014.
 [4] F. Rusek, D. Persson, B. K. Lau, E. G. Larsson, T. L. Marzetta, O. Edfors, and F. Tufvesson, “Scaling up mimo: Opportunities and challenges with very large arrays,” Signal Processing Magazine, IEEE, vol. 30, no. 1, pp. 40–60, 2013.
 [5] C. Zhang, X. Liang, Z. Wu, F. Wang, S. Zhang, Z. Zhang, and X. You, “On the lowcomplexity, hardwarefriendly tridiagonal matrix inversion for correlated massive mimo systems,” arXiv preprint arXiv:1802.10444, 2018.
 [6] X. Gao, L. Dai, Y. Hu, Z. Wang, and Z. Wang, “Matrix inversionless signal detection using sor method for uplink largescale mimo systems,” in Global Communications Conference (GLOBECOM), 2014 IEEE. IEEE, 2014, pp. 3291–3295.
 [7] B. Yin, M. Wu, C. Studer, J. R. Cavallaro, and C. Dick, “Implementation tradeoffs for linear detection in largescale mimo systems,” in IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP). IEEE, 2013, pp. 2679–2683.
 [8] A. Elgabli, A. Elghariani, A. O. AlAbbasi, and M. Bell, “Twostage lasso admm signal detection algorithm for large scale mimo,” in Signals, Systems, and Computers, 2017 51st Asilomar Conference on. IEEE, 2017, pp. 1660–1664.
 [9] Z. Wang, V. Aggarwal, and X. Wang, “Optimal energybandwidth allocation for energyharvesting networks in multiuser fading channels,” IEEE Journal on Selected Areas in Communications, vol. 34, no. 5, pp. 1565–1577, 2016.