arxiv: v1 [cs.it] 26 Jan 2016

Transcription

1 SYMPTOTIC NLYSIS OF DOWNLIN MISO SYSTEMS OVER RICIN FDING CNNELS ugo Falconet,LucaSanguinetti,blaammoun +,MerouaneDebbah Mathematical and lgorithmic Sciences Lab, uawei France, Paris, France Large Networks and System Group LNES, CentraleSupélec, Université Paris-Saclay, Gif-sur-Yvette, France Dipartimento di Ingegneria dell Informazione, UniversityofPisa,Pisa,Italy + Electrical Engineering Department, ing bdullah University of Science and Technology, Thuwal, Saudi rabia arxiv: v [cs.it] 26 Jan 206 BSTRCT In this work, we focus on the ergodic sum rate in the downlink of asingle-celllarge-scalemulti-usermimosysteminwhichthe base station employs N antennas to communicate with single-antenna user equipments. regularized zero-forcing RZF scheme is used for precoding under the assumption that each link forms a spatially correlated MIMO Rician fading channel. The analysis is conducted assuming N and grow large with a non trivial ratio and perfect channel state information is available at the base station. Recent resultsfromrandommatrixtheoryand largesystemanalysis are used to compute an asymptotic expression of the signal-to-interferenceplus-noise ratio as a function of the system parameters, the spatial correlation matrix and the Rician factor. Numerical results are used to evaluate the performance gap in the finite system regime under different operating conditions.. INTRODUCTION Large-scale multiple-input multiple-output MIMO systems also known as massive MIMO systems are considered as one of the most promising technology for next generation wireless communication systems [ 3] because of their considerable spatial multiplexing gains. The use of large-scale MIMO systems is beneficial not only in terms of coverage and spectral efficiency but also in terms of energy-saving [4 6]. In this complex system model, a number of practical factors such as correlation effects and line-of-sight LOS components need to be included, which occur due to the space limitation of user equipments UEs and the densification of the antenna arrays resulting in a visible propagation path from the UEs, respectively. For typical systems of hundreds of antennas and tens of UEs, even computer simulations become challenging, which makes performance analysis of large-scale MIMO systems an important and a new subject of research. In this work, we consider the downlink of a single-cell largescale MIMO system in which the base station BS, equipped with N antennas, makes use of regularized zero-forcing RZF precoding to communicate with single-antenna UEs. In particular, we are interested in evaluating the ergodic sum rate of the system when a power constraint is imposed at the BS. The analysis is conducted assuming that N and grow large with a non trivial ratio under the assumption that perfect channel state information is available at the BS. Differently from most of the existing literature [7 ], we consider aspatiallycorrelatedmimoricianfadingmodel,whichismore This research has been supported by the ERC Starting Grant MORE and by the UWEI project RMTin5G. general and accurate to capture the fading variations when there is a LOS component. Compared to the Rayleigh fading channel, a Rician model makes the asymptotic analysis of large-scale MIMO systems much more involved. To overcome this issue, recent results from random matrix theory and large system analysis [ 3] are used to compute an asymptotic expression of the signal-to-interferenceplus-noise ratio SINR, which is eventually used to approximate the ergodic sum rate of the system. s shall be seen, the results are found to depend only on the system parameters, the spatial correlation matrix and the Rician factor. s a notable outcome of this work, the above analysis provides an analytical framework that can be used to evaluate the performance of the network under different settings without resorting to heavy Monte Carlo simulations and to eventually get insights on how the different parameters affect the performance. The main literature related to this work is represented by [7,0, 4 6]. Tools from random matrix theory are used in [7] to compute the ergodic sum rate in a single-cell setting with Rayleigh fading and different precoding schemes while the multicell case isanalyzed in [0]. In [4], the authors investigate a LOS-based conjugate beamforming transmission scheme and derive some expressions of the statistical SINR under the assumption that N grows large and maintain fixed. In [5], the authors study the fluctuations of the mutual information of a cooperative small cell network operating over a Rician fading channel under the form of a central limit theorem and provide explicit expression of the asymptotic variance. In [6], a deterministic equivalent of the ergodic sum rate and an algorithm for evaluating the capacity achieving input covariance matrices for the uplink of a large-scale MIMO are proposed for spatially correlated MIMO channel with LOS components. The analysis of the uplink rate with both zero-forcing and maximum ratio combining receivers is performed in [7]. The following notation is used throughout this work. The superscript stands for the conjugate transpose operation. The operator TrX denotes trace of matrix X whereas X [k] indicates that the kth column is removed from matrix X. TheN N identity matrix is denoted by I N whereas X =diag{x,...,x N } is used to denote a N N diagonal matrix of entries {x n}. randomvariablex is a standard complex Gaussian variable if CN0,. 2. SYSTEM ND CNNEL MODEL We consider the downlink of a network in which UEs are served by a single BS equipped with N antennas. The signal y k received

2 by UE k takes the form [7] y k = h k g k s k + n k k {,...,} where h k C N is the random channel from the BS to user k, g k C N is the precoding vector associated to UE k, andn k CN0,σ 2 accounts for thermal noise. For analytic tractability, we assume that the BS is able to acquire perfect channel state information from the uplink pilots. The RZF scheme is used as precoding technology [7]. Therefore, the precoding matrix G =[g...g ] C N takes the form G = ξ + λi N 2 where =[h,...,h ] C N is the aggregate channel matrix, λ>0 is the so-called regularization parameter and ξ is chosen so as to satisfy the average total power constraint given by TrG G = where > 0 denotes the total available transmit power. Therefore, it follows that ξ 2 = Tr + λi N 2. 3 Under the assumption of Gaussian signaling, i.e., s k CN0, and single user detection with perfect channel state information at the receiver, the SINR γ k of user k takes the form γ k = h k g k 2 h k GG h k h k g k 2 + σ 2. 4 The rate r k of UE k is given by r k =log 2 +γk whereas the ergodic sum rate is defined as 3. MIN RESULT We aim to exploit the statistical distribution of the channel and the large dimensions of N and to compute the deterministic approximation of γ k,whichwillbeeventuallyusedtofindanapproximation of the ergodic sum rate. To begin with, we call D =diag{ β +ρ,..., β +ρ and compute the eigendecomposition of Θ /2 to obtain } 8 Θ /2 = U D /2 U 9 with U C N N being unitary. Then, we rewrite as follows = U Σ where Σ is defined as with X = UZ, = U Z and Σ = D /2 X D /2 + 0 ρ [ Z = ] β z,..., β z. +ρ lso, we assume the following grow rate of system dimensions: ssumption. The dimensions N and grow to infinity at the same pace, that is: lim inf N/ lim sup N/ <. 2 For technical reasons, the following reasonable assumptions are also imposed on the system settings [7, 2, 3]. ssumption 2. s N,,thematricesD and D have uniformly bounded spectral norm, i.e., sup D < sup D <. 3 r E = E [ ] log 2 +γk 5 Moreover, inf TrD > 0 inf Tr D > 0. 4 where the expectation is taken over the random channels {h k }.We assume that h k = β k w k 6 where β k accounts for the large scale channel fading of UE k and w k C N is the small scale fading channel. The latter is modelled as ρ w k = +ρ Θ/2 z k + +ρ z k 7 where z k C N is assumed to be Gaussian with zero mean and unit covariance, i.e., z k CN0 N, I N, and z k C N is a deterministic vector. The scalar ρ 0 is the Rician factor whereas the matrix Θ /2 is obtained from the Cholesky decomposition of Θ C N N, which accounts for the channel correlation matrix at the BS antennas. To make the problem analytically more tractable, we consider asystemwithcommonuechannelcorrelationmatrix[7,8,9]. lthough possible, the extension to the case in which UEs have different channel correlation matrices is mathematically muchmoreinvolved and it is left for future research. lso, we assume that the family of deterministic N matrices has bounded spectral norm: sup <. 5 Let us now introduce the fundamental equations that are needed to express a deterministic equivalent of γ k.westartwiththefollowing set of equations: δ = λ TrD I N + δd + I + δ D δ = Tr D λ I + δ D + I N + δd which admits a unique positive solution in the class of Stieltjes transforms of non-negative measures with support R + [2, 3]. The matrices T = λ I N + δd + I + δ D T = λ I + δ D + I N + δd

3 s k =ū k [ λ t 2 kk a k T 2 a k F + δ d k + 2 α + ϑ V λ dk + F ϑα + V + λ2 λ t 2 kk d i λ t 2 i=,i k k,k a k TDTa k + δ d k ] a k Ta i 2 +δ d i +δ d + k 6 ψ = [ TrT λ TrT2 +2 F αv + ϑ ] V 2 + λ2 2 ϑα. 7 are approximations of the resolvent Q =ΣΣ zi N and the co-resolvent Q =Σ Σ zi.define ϑ = TrDTDT 8 ϑ = Tr D T D T 9 F = I TrTDT + δ D 2 D 20 = F λ 2 ϑ ϑ 2 α = TrDT2 22 V = D 2 TrT2 I + δ D. 23 lso, denote by a k the k-th column of matrix and call t kk, d k and d k the k-diagonal element of T, D and D, respectively. Then,the following theorem summarizes the main result of this work. Theorem. Let ssumptions 2 hold true. Then, we have that max k γ k γ k 0 almost surely with with ū k given by γ k = u 2 k s k + ψ + ū k σ 2 24 ū k = δ d k + a k Ta k λ t kk + δ d k whereas s k and ψ are given on the top of this page. 25 Proof. The proof is very much involved and relies on results in random matrix theory [2] as well as some recent ones on the deterministic equivalent of bilinear forms [3, Theorem ]. Due to the space limitations, it is omitted. complete proof is provided in [20]. Corollary. Let ssumptions 2 hold true. If Θ = I N,then s k and ψ k reduce to [ F a k T 2 a k s k =ū k λ t 2 kk +δ d + ϑ λ d k + k ] d i a k Ta i 2 + +δ d i +δ d 26 k and λ t 2 i=,i k kk ψ = TrT λ TrT2. 27 Proof. The result follows easily observing that if Θ = I N then F = V = TrT2 I + δ D 2 D and ϑ = α = TrT2. Let us now consider the case in which only the link of the user of interest, i.e., UE k, ischaracterizedbyaricianfadingchannel model whereas a Rayleigh fading channel model is considered for all the other UEs with indexes i k. ssumealsothat z k 2 = N such as for example when a uniform linear array is used. Then, the following result is obtained. Corollary 2. Let ssumptions 2 hold true. If Θ = I N, z k 2 = N and no LOS component is present i k,thenγ k reduces to γ k = λ N δ δβ i N i= +δβ i δβ k δβ k + δ + δβ k σ2 with δ being the unique positive solution of the following equation: δ = N λ + β i. 28 N + δβ i i= The above expression coincides with that in [7] obtained for the case in which all UEs experience Rayleigh fading. This indicates that the asymptotic expression of the SINR for a given user interest is independent from the underlying channel model if all the others are affected by Rayleigh fading and uniform linear array is used. We are ultimately interested in the individual rates {r k } and the ergodic sum rate r E.Sincethelogarithmisacontinuousfunction, by applying the continuous mapping theorem, from the almost sure convergence results of Theorem it follows that r k r k 0 almost surely with [7] r k =log 2 +γk. 29 n approximation of r E is obtained as follows [7] r E = log 2 +γk such that re re 0 holds true almost surely. 4. NUMERICL RESULTS 30 Monte-Carlo simulations are now used to validate the above asymptotic analysis for a network with finite size. We consider a cell To enable simple testing of other parameter values, the code is also available for download at the following address

4 verage rate per UE [bps/z] ρ = ρ =/2 ρ =0 verage rate per UE [bps/z] ν =0.9 =8 ν =0 =8 ν =0 =6 ν =0.9 =6 2 E log 2 + γ k log 2 + γ k E log 2 + γ k log 2 + γ k Number of transmit antennas, N Number of transmit antennas, N Fig.. verage rate per UE vs N when =6and the Rician factor ρ is 0, /2 and Fig. 2. veragerateperuevsn when =8, 6 and ρ =for different values of the correlation factor ν among BS antennas. of radius R = 250 m. We use the standard correlation model [Θ] i,j = ν i j and assume that the large scale coefficients β k are obtained as β k =2L x + x κ k/ x κ.theparameterκ>2is the path loss exponent, x k denotes the distance of UE k from the BS, x >0 is some cut-off parameter and L x is a constant that regulates the attenuation at distance x. We assume κ =3.5, L x = 86.5 db and x =25m. The results are obtained for 000 different channel realizations assuming that the UEs are randomly distributed in the cell. We assume that a uniform linear array is adopted at the BS such that nth entry of z k is [ z k ] n =exp in π sin θ k where θ k is the arrival angle of UE k taking values in the interval [ π,π. Thetransmitpower is fixed to 0 Watt and the regularization parameter λ is computed as λ = σ 2 E [ ] β k /PT,wherethe expectation is taken over the distribution of the users positions. Fig. illustrates the average per UE rate when N grows large and is kept fixed to 6. Thecorrelationfactorissettoν =0.9 and the Rician factor is ρ =0, /2 and. sseen,thetheoretical asymptotic results match very well with Monte Carlo simulations over the entire of N, alsowhenn takes relatively small values. s expected, increasing the Rician factor ρ improves the system performance. The impact of the correlation matrix is analyzed in Fig. 2 for ρ =and =8, 6. TheUEsaredistributeduniformlyonacircle of radius 2/3R around the BS. s expected, reducing the correlation factor ν improves the system performance as it increases the spatial multiplexing capabilities of the channel model. 5. CONCLUSION In this work, we analyzed the ergodic sum rate in the downlink of a single-cell large-scale MIMO system operating over a Rician fading channel. regularized zero-forcing precoding scheme under the assumption of perfect channel state information is used. Recent results from large-scale random matrix theory allowed us to give concise approximations of the SINRs. Such approximations turned out to depend only on the long-term channel statistics, the Rician factor and the deterministic component. Numerical results indicated that these approximations are very accurate. pplied to practical net- works, such results may lead to important insights on how the different parameters affect the performance and allow to simulate the network behavior without the need of extensive Monte Carlo simulations. More details and insights on these aspects will be given in the extended version [20]. 6. REFERENCES [] T. Marzetta, Noncooperative cellular wireless with unlimited numbers of base station antennas, IEEE Transactions on Wireless Communications, vol.9,no.,pp ,Nov [2] F. Rusek, D. Persson, B.. Lau, E. Larsson, T. Marzetta, O. Edfors, and F. Tufvesson, Scaling up MIMO: Opportunities and challenges with very large arrays, IEEE Signal Processing Magazine, vol.30,no.,pp.40 60,Jan203. [3] E. Larsson, O. Edfors, F. Tufvesson, and T. Marzetta, Massive MIMO for next generation wireless systems, IEEE Communications Magazine,vol.52,no.2,pp.86 95,February204. [4]. Ngo, E. Larsson, and T. Marzetta, Energy and spectral efficiency of very large multiuser MIMO systems, IEEE Trans. Commun., vol.6,no.4,pp ,203. [5] E. Björnson, L. Sanguinetti, J. oydis, and M. Debbah, Optimal design of energy-efficient multi-user MIMO systems: Is massive MIMO the answer? IEEE Trans. Wireless Commun., vol. 4, no. 6, pp , 205. [6] E. Björnson, L. Sanguinetti, and M. ountouris, Deploying dense networks for maximal energy efficiency: Small cells meet massive MIMO, 206, to appear. [Online]. vailable: [7] S. Wagner, R. Couillet, M. Debbah, and D. T. M. Slock, Large system analysis of linear precoding in correlated MISO broadcast channels under limited feedback, IEEE Trans. Inf. Theory, vol.58,no.7,pp ,July202.

5 [8] J. oydis, S. ten Brink, and M. Debbah, Massive MIMO in the UL/DL of cellular networks: ow many antennas do we need? IEEE J. Sel. reas Commun., vol. 3, no. 2, pp. 60 7, Feb [9] L. Sanguinetti, R. Couillet, and M. Debbah, Large system analysis of base station cooperation in the downlink, IEEE Trans. Wireless Comm., 205,underreview.[Online]. vailable: [0] J. Zhang, C.-. Wen, S. Jin, X. Gao, and.-. Wong, Large system analysis of cooperative multi-cell downlink transmission via regularized channel inversion with imperfect CSIT, IEEE Trans. Wireless Commun.,vol.2,no.0,October203. [] R. Couillet and M. Debbah, Random matrix methods for wireless communications. Cambridge, New York: Cambridge University Press, 20. [2] W. achem, M. harouf, J. Najim, and J. W. Silverstein, CLT for information-theoretic statistics of non-centered Gram random matrices, Random Matrices: Theory and pplications, vol., no. 2, pp. 50, 202. [Online]. vailable: [3] W. achem, P. Loubaton, J. Najim, and P. Vallet, On bilinear forms based on the resolvent of large random matrices, nnales de l Institut enri Poincar, Probabilities et Statistiques, vol. 49, no., pp , 203. [4] D. Yue and G. Y. Li, LOS-based conjugate beamforming and power-scaling law in massive-mimo systems, CoRR, vol.abs/ ,204.[online].vailable: [5] J. oydis,. ammoun, J. Najim, and M. Debbah, Outage performance of cooperative small-cell systems under Rician fading channels, in IEEE 2th International Workshop on Signal Processing dvances in Wireless Communications SPWC,June20,pp [6] J. Zhang, C.-. Wen, S. Jin, X. Gao, and.-. Wong, On capacity of large-scale MIMO multiple access channels with distributed sets of correlated antennas, IEEE J. Sel. reas Commun.,vol.3,no.2,pp.33 48,February203. [7] Q. Zhang, S. Jin,.-. Wong,. Zhu, and M. Matthaiou, Power scaling of uplink massive MIMO systems with arbitrary-rank channel means, IEEE J. Sel. Topics Signal Process.,vol.8,no.5,pp ,Oct204. [8]. dhikary, J. Nam, J.-Y. hn, and G. Caire, Joint spatial division and multiplexing the large-scale array regime, IEEE Trans. Inf. Theory, vol.59,no.0,pp ,Oct203. [9]. ammoun,. Muller, E. Bjornson, and M. Debbah, Linear precoding based on polynomial expansion: Large-scale multicell mimo systems, IEEE J. Sel. Topics Signal Process., vol.8, no. 5, pp , Oct 204. [20]. ammoun, L. Sanguinetti, and M. Debbah, symptotic analysis of the downlink of large-scale MIMO systems over Rician fading channels, IEEE Trans. Wireless Commun., Jan. 206, under review. [Online]. vailable: