The classical detection techniques for multiple-input multiple-output (MIMO) systems are usually designed with the assumption that the additive complex Gaussian noise is uncorrelated. However, for closely spaced antennas, the additive noise is correlated due to the mutual antenna coupling. This letter analyzes an improved zero-forcing (ZF) technique for MIMO channels in colored environments. The additive noise is assumed to be correlated and the Rayleigh MIMO channel is considered doubly correlated. The improved ZF detector, based on the generalized least squares estimator (GLS), takes into account the noise covariance matrix and provides an unbiased estimator of the transmitted symbol vectors. We introduce some novel bounds on the achievable sum rate, on the normalized mean square error at the receiver output, and on the outage probability. The derived expressions are compared to Monte Carlo simulations.
1. Introduction
Multiple-input multiple-output (MIMO) systems have great potential for increasing spectral efficiency [1]. To achieve the promised capacity, efficient and reliable reception techniques are required. The linear zero-forcing (ZF) technique is one of the most popular receivers [2–6]. In fact, despite its low-complexity, the ZF detector yields suboptimal performances [2]. The analysis of the MIMO-ZF detector has attracted considerable attention in the literature [2–6]. The effect of channel estimation error of the MIMO-ZF detector performance has been investigated in [2–4]. The SNR loss due to the spatial correlation and its effect on bit error rate (BER) of the MIMO-ZF detector are discussed in [4–6]. Although previous works present important contributions concerning the MIMO-ZF analysis, they assume that the channel additive complex Gaussian noise is uncorrelated. The noise vector entries are then modeled as independent and identically distributed (i. i. d.) random variables. This assumption is valid when the antennas are sufficiently spaced and thus coupling between them is negligible. However, for MIMO systems with closely spaced antennas, there is a mutual coupling between antenna elements [7]. This mutual coupling leads to a correlated thermal noise at the receiver front-end.
Motivated by the discussion above, this letter proposes and analyzes a ZF detection technique for MIMO systems in correlated environments. We assume that the additive noise is correlated and the Rayleigh MIMO channel has both transmit and receive correlations. We also assume that the fading channel matrix is perfectly known at the receiver. When the noise covariance matrix is known, the ZF detector corresponds to the generalized least squares estimator (GLS) of the transmitted symbols. When the noise covariance is unknown, we propose to replace the true matrix by the estimated one. The ZF detector is then called the feasible generalized least squares estimator (FGLS). This letter presents some bounds for the achievable sum rate (ASR), the normalized mean square error (NMSE) at the GLS receiver output, and the global outage probability (GOP). It is noteworthy that the scenario with doubly correlated channel has been rarely studied [2]. In fact, for this case, it is relatively difficult to describe the statistical properties of the postprocessing SNR [2]. The few relevant studies are given in [8] for 2×2 MIMO systems and several fading models in [2]. Some results on the eigenstatistics for doubly correlated MIMO channels are reported in [9].
2. System Model and ZF Receiver2.1. System Model
Consider an uncoded MIMO system with Nt transmit and Nr receive antennas with Nr>Nt. The received signal vector yt∈CNr×1 at time t is
(1)yt=Hxt+nt,t=1,2,…,T,
where T is the number of received samples in the observation window or the frame length, xt∈CNt×1 is the transmitted vector, H∈CNr×Nt is the channel matrix, and nt∈CNr×1 is a circularly symmetric complex Gaussian noise vector with zero mean and covariance matrix Σ=E[ntntH]=σn2Ω. We assume that the noise covariance matrix is given by the widely used noise model (see [10] and the references therein):
(2)Σ=EntntH=σn2Ω,
where Ω is a positive definite matrix, with the (p,q)th entry
(3)Ωp,q=0.9p-qexpjπ2p-q.
The transmitted symbols are uncorrelated with covariance Σx=E[xtxtH]=EsINt. We consider a doubly spatially correlated Rayleigh flat-fading channel. With this assumption, spatial correlation occurs on both transmit and receive sides and the channel matrix H is modeled as
(4)H=Φr1/2HwΦt1/2,
where Hw is an Nr×Nt matrix with i. i. d. zero-mean complex Gaussian entries with unit variance. The spatial correlation matrices Φt and Φr are defined as [Φt]i,j=ρt|i-j| and [Φr]i,j=ρr|i-j|.
When the channel matrix H is available at the receiver, the classical ZF detector estimates the transmitted vector as
(5)x^to=H†yt=xt+H†nt,
where H† is the pseudoinverse of the channel matrix. The estimator x^to is an unbiased estimator with covariance matrix [11]
(6)varx^to=HHH-1HHΣHHHH-1.
The estimator x^to is called the ordinary least squares (OLS) estimator [12]. If the noise is uncorrelated, that is, Σ=σn2INr, the OLS estimator corresponds to the best linear unbiased estimator (BLUE) [12]. However, when the noise is correlated there is no guarantee that the OLS estimator is the BLUE.
2.2. Generalized Least Squares Estimator with Known Ω
Since the OLS estimator is no longer BLUE, the ZF receiver degrades the system BER. Using the Gauss-Markov theorem and assuming that H is the full rank matrix, we can prove that the BLUE estimator of xt is given by [12]
(7)x^tg=HHΣ-1H-1HHΣ-1yt.
The GLS estimator x^tg is unbiased with covariance matrix [12]
(8)varx^tg=Ext-x^tgxt-x^tgH=HHΣ-1H-1.
We can easily show that (the matrix inequality A≤B means that B-A is positive semidefinite)(9)varx^tg≤varx^to.
As can be seen in (7), the GLS estimator needs the knowledge of the noise covariance matrix Σ. By combining (7) and (2), we can write the GLS estimator as
(10)x^tg=HHΩ-1H-1HHΩ-1yt,
which does not require the knowledge of σn2.
3. Performance Analysis3.1. Derivation of the Sum Rate
The analysis of the ergodic sum rate of the ZF receiver has been widely discussed [13, 14]. However, previous works have considered systems with uncorrelated noise. In this section, we introduce a novel upper and lower bounds on the ASR of the improved ZF detector. Assuming independent decoding at the receiver, the ASR is
(11)R=∑k=1NtElog21+γk,
where γk is the SNR of the kth stream. Since the covariance matrix of (8) is the MSE of the data estimation, we have
(12)γk=γ0Z-1kk,
where γ0=Es/σn2 and Z=HHΩ-1H. When Hw is a complex Gaussian matrix, it is well known that HwHHw follows the central complex Wishart distribution. However, it is generally hard to derive the distribution of HHΩ-1H. In this case, we will use the decomposition det(HHΩ-1H)=det(Φt)det(HwHΦrH/2Ω-1Φr1/2Hw) and the first moment of the generalized variance given in [13]
(13)EdetHwHAHw=Nt!detW∏i<jNrλj-λi,
where A is an arbitrary Hermitian positive definite matrix, λi are the eigenvalues of A, and W is defined by
(14)Wi,j=λij-1,j=1,…,Nr-Nt,λij,j=Nr-Nt+1,…,Nr.
Let us define A=ΦrH/2Ω-1Φr1/2.
Proposition 1.
The ASR of the GLS receiver in a correlated environment is upper bounded by
(15)Ru=∑m=1NtNtlog2Nt-1!detW~+Ntγ0/σmdetW∏i<jNrλj-λi-Ntln2∑k=1Nt-1ψk+∑k=Nr-Nt+2NrdetVk∏i<jNrλj-λi,
where ψ(·) is Euler’s digamma function, σm is the mth diagonal entry of Φt-1, W~ is defined as
(16)W~i,j=λij-1,j=1,…,Nr-Nt+1,λij,j=Nr-Nt+2,…,Nr,
and V is an Nr×Nr matrix with elements
(17)Vki,j=λij-1,j≠k,λij-1lnλi,j=k.
Furthermore, the ASR is lower bounded by
(18)Rl=∑m=1Ntlog21+γ0expψNt-lnσmdetVNr-Nt+1∏i<jNrλj-λiWWWWWWW000WW+detVNr-Nt+1∏i<jNrλj-λi.
Proof.
The proof is based on a similar approach as that in [13, Appendix III], more specifically, by assuming A=ΦrH/2Ω-1Φr1/2 instead of A=Φr.
3.2. Normalized Mean Square Error
Let us define the conditioned normalized MSE (NMSE) at the ZF receiver output, given the channel matrix, as
(19)NMSEH=trExt-x^tgxt-x^tgHNt=trZ-1Nt.
The NMSE is obtained by averaging (19) with respect to H.
Proposition 2.
The NMSE at the GLS receiver output in a correlated environment is lower bounded by
(20)NMSE≥1Nt∑k=1NtexpψNt-lnσm+detYNr-Nt+1∏i<jNrγj-γi.
Proof.
Since trace and expectation commute, we can write
(21)NMSE=1Nt∑k=1NtEexplnZ-1kk(22)=1Nt∑k=1NtEexplndetHmHHmdetHHH(23)≥1Nt∑k=1NtexpElndetHmHHmdetHHH,
where Hm is the (m,m)th minor of H and (22) is obtained from matrix properties, while (23) is obtained thanks to Jensen’s inequality since exp(·) is a convex function. The proof follows by invoking a lemma form [15, Lemma 4] and a result from [13, equations (15), (46)] which give
(24)ElndetHHH-lndetHmHHm=ψNt-lnσm+detYNr-Nt+1∏i<jNrγj-γi.
Substituting (24) into (23), we obtain (20).
3.3. Global Outage Probability
There are several definitions of the GOP which depend on the application of interest [16]. In this letter, we assume that the GOP is defined as the probability that all instantaneous SNRs γk fall below a fixed threshold γ¯. This means that the system is considered available if at least one of the substreams achieves the desired performance. This definition, called all-outage probability in [16], can be expressed as
(25)Poγ¯=PrD≤γ¯INt,
where D=Diagγ1,…,γNt.
Proposition 3.
The GOP of the GLS receiver in a correlated environment is upper bounded by
(26)Poγ¯≤1lnγ¯∑k=1NtψNt-lnσmγ0+detYNr-Nt+1∏i<jNrγj-γi.
Proof.
We start by writing Po(γ¯)=PrlnD≤lnγ¯INt since ln(·) is a strictly monotonic increasing function and using the Markov’s inequality [17](27)PrgD≤αI≤1αtrEgD.
Then, we can write the following:
(28)Poγ¯≤1lnγ¯trElnD=1lnγ¯EtrlnD(29)≤1lnγ¯∑k=1NtElnγ0Z-1kk(30)≤Ntlnγ0lnγ¯-1lnγ¯∑k=1NtlnEdetHmHHmdetHHH,
where (29) results from Jensen’s inequality since ln(·) is concave function. The bound in (26) is obtained using (24).
4. Improved ZF Receiver with Unknown Ω
When Ω is unknown at the receiver, we estimate the covariance matrix Σ and use the estimator Σ^T (the subscript T refers to the frame length) to evaluate the feasible generalized least squares estimator (FGLS) [12]
(31)x^tf=HHΣ^T-1H-1HHΣ^T-1yt.
To obtain the FGLS estimator, we propose to use the following iterative process (with ni iterations).
Perform the OLS estimation x^t(1)=H†yt.
Compute the residual vector et(1) for (t=1,…,T)(32)et1=yt-Hx¯t1,
where x¯t(1) is the hard estimate of xt obtained by mapping x^t(1) to the closest constellation point using the classical decision rule for the considered modulation.
Estimate Σ by Σ^T(1)=1/T∑t=1Tet(1)(et(1))H.
For any iteration k=2,…,ni, compute the FGLS
(33)x^tk=HHΣ^Tk-1H-1HHΣ^Tk-1yt.
Update the residual error vector as et(k)=yt-Hx¯t(k).
Update the covariance matrix estimation
(34)Σ^Tk=1T∑t=1TetketkH.
Iterate Steps (3), (4), (5), and (6).
When Σ is replaced by Σ^T, the Gauss-Markov theorem no longer holds [12]. This means that the estimator x^tT is not necessarily BLUE. In general, it is difficult to evaluate the statistical properties of x^tT. Monte Carlo simulations can be used to check the accuracy of the FGLS estimator.
5. Numerical Results
In this section, Monte Carlo simulations are employed to compare the OLS, GLS, and FGLS estimators and support our derivation of the bounds. We consider a 2×4 MIMO system with frame size T=104. All simulation curves are obtained over 104 independent MIMO channel realizations. The transmitted symbols are 16-QAM modulated. The correlation coefficients are fixed to ρt=ρr=0.6.
In Figure 1, the simulated ASR for the three detection techniques are plotted. For the FGLS case, the plotted curve is obtained after the fifth iteration. For the GLS method, the ASR obtained by the simulation is also compared against the bounds of (15) and (18). It is shown that the GLS and the FGLS outperform the OLS receiver. For higher SNRs, the ASR obtained when Σ is unknown is nearly the same with that when Σ is known. Furthermore, we can see that the analytical bounds exhibit a negligible difference with the simulation curves. However, the lower bound is tighter across all SNRs.
Sum rate of OLS, GLS (simulated and bounds), and FGLS receivers.
Figure 2 compares the NMSE performance of the OLS, GLS, and FGLS approaches. It shows that the lower bound of (20) works well for γ0=6 dB. For the FGLS method, we have plotted the NMSE after each iteration. It can be seen that, from the first iteration to the fourth one, there is an important improvement in the NMSE. However, after the fourth iteration, the NMSE improvement is not noticeable. It is thus possible to shut down the iteration process after the fourth iteration.
NMSE performance of OLS, GLS, and FGLS receivers.
The GOPs of the different receivers are plotted in Figure 3. The upper bound of (26) is also provided. The threshold γ¯ is fixed in order to assure a target BER of 10-4. We can see that the upper bound of (26) is fairly tight and can well predict the GOP. At a GOP of 10-4, the GLS technique allows a significant gain of about 3 dB with respect to OLS estimator. The GOP of the FGLS receiver after the fifth iteration is close to that of the GLS.
GOP versus SNR for target BER of 10-4.
6. Conclusion
A ZF detector is analyzed in the presence of correlated noise and doubly correlated channel. When the noise covariance matrix is known, the GLS algorithm allows a BLUE estimation of the transmitted sequence. In this case, we provide upper and lower bounds of the ASR, a lower bound of the NMSE, and an upper bound of the GOP. When the noise covariance matrix is unknown, we propose to use the iterative FGLS algorithm. Computer simulations have demonstrated the efficiency of the proposed FGLS detector which leads to a substantial improvement in the system performance. After four iterations, the FGLS performances approach those of the GLS. We have validated the derived analytic bounds by computer simulations.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
FoschiniG. J.GansM. J.On limits of wireless communications in a fading environment when using multiple antennas19986331133510.1023/A:10088892227842-s2.0-0032022195WangC.AuE. K. S.MurchR. D.MowW. H.ChengR. S.LauV.On the performance of the MIMO zero-forcing receiver in the presence of channel estimation error20076380581010.1109/TWC.2007.053842-s2.0-33947576428GoreD. A.HeathR. W.Jr.PaulrajA. J.Transmit selection in spatial multiplexing systems200261149149310.1109/LCOMM.2002.8055172-s2.0-0036864743SiriteanuC.MiyanagaY.BlosteinS. D.KurikiS.ShiX.MIMO zero-forcing detection analysis for correlated and estimated rician fading20126173087309910.1109/TVT.2012.22021332-s2.0-84866508676SalehH. A. A.HamoudaW.Performance of zero-forcing detectors over MIMO flat-correlated Ricean fading channels200931101610.1049/iet-com:200801102-s2.0-58149481426McKayM. R.CollingsI. B.Error performance of MIMO-BICM with zero-forcing receivers in spatially-correlated Rayleigh channels20076378779210.1109/TWC.2006.052592-s2.0-33947587685KrusevacS. M.KennedyR. A.RapajicP. B.Effect of signal and noise mutual coupling on MIMO channel capacity200740331732810.1007/s11277-006-9193-82-s2.0-33846570303KiesslingM.SpeidelJ.Analytical performance of MIMO zeroforcing receivers in correlated Rayleigh fading environmentsProceedings of the 4th IEEE Workshop Signal Processing Advances in Wireless Communications (SPAWC '03)2003383387McKayM. R.GrantA. J.CollingsI. B.Performance analysis of MIMO-MRC in double-correlated Rayleigh environments200755349750710.1109/TCOMM.2007.8924502-s2.0-34147190205VibergM.Maximum likelihood array processing in spatially correlated noise fields using parameterized signals1997454996100410.1109/78.5641872-s2.0-0031123069KayS. M.19931stPrentice HallTimmN. H.2002Berlin, GermanySpringerMatthaiouM.ZhongC.RatnarajahT.Novel generic bounds on the sum rate of MIMO ZF receivers20115994341435310.1109/TSP.2011.2157148MR28659882-s2.0-80051777051MatthaiouM.ZhongC.McKayM. R.RatnarajahT.Sum rate analysis of ZF receivers in distributed MIMO systems201331218019110.1109/JSAC.2013.1302072-s2.0-84872942680JinS.McKayM.ZhongC.WongK.-K.Ergodic capacity analysis of amplify-and-forward MIMO dual-hop systems20105652204222410.1109/TIT.2010.2043765MR2729781OrdonezL. G.PalomarD. P.FonollosaOrdered eigenvalues of a general class of Hermitian random matrices with application to the performance analysis of MIMO systems200957267268910.1109/TSP.2008.2007602MR26033932-s2.0-60649087881AhlswedeR.WinterA.Strong converse for identification via quantum channels200248356957910.1109/18.985947MR18899692-s2.0-0036495139