Effect of AntennaMutual Coupling on MIMOChannel Estimation and Capacity

This paper describes investigations into the antenna mutual coupling (MC) effect on channel estimation and capacity of a multipleinput multiple-output (MIMO) wireless communication system. The presented investigations close the gap existing in the previous works which assessed the effect of mutual coupling on MIMO capacity under the assumption of availability of perfect channel state information (CSI) at the receiver. The new approach assumes that the perfect CSI is not available due to channel estimation errors. The investigations are carried out for different spacing between array antenna elements producing a varying effect of mutual coupling on the channel estimation and the resulting MIMO channel capacity.


Introduction
It has been shown in recent works that a wireless communication system can achieve an increased capacity by using multiple element antennas at transmitter and receiver by applying a suitable signal processing scheme [1].This unconventional wireless communication system is called a multiple-input multiple-output (MIMO) system.It offers an increased capacity without the need of increasing an operational bandwidth.
In order to realize the advantages of MIMO, two conditions have to be satisfied.One requires the presence of a rich scattering environment, and the other one entails accurate channel state information (CSI) to be available at the receiver [2].The rich scattering environment is necessary to support the formation of statistically independent virtual channels over which the parallel data transmission can take place.The lack of (spatial) correlation between the virtual channels leads to the increased MIMO capacity.The availability of accurate CSI is required to decode the received signal and to practically reach the MIMO capacity [1][2][3][4][5].In turn, an inaccurate CSI leads to an increased bit error rate (BER) that translates into a degraded capacity of the system [6][7][8].
It is known that the finite antenna spacing in array antennas introduces spatial correlation.This finite spacing is also responsible for mutual coupling which adversely affects signal transmission and reception due to the resulting antenna impedance mismatch.The mutual coupling effect is especially pronounced in tightly spaced arrays.Because of a considerable demand for compact size mobile station (MS) terminals, the effect of mutual coupling cannot be neglected and thus has to be taken into account while assessing the MIMO link performance.The problem of mutual coupling in MIMO systems for the case of peerto-peer communication has been addressed via simulations and measurements in [9][10][11][12][13][14][15].It has been shown that the mutual coupling may improve the channel capacity when the spacing of antenna elements in the transmit and receive array antennas is between 0.2 to 0.4λ (carrier wavelength).The reason is that in this case the mutual coupling decreases the spatial correlation and increases the channel effective degree of freedom (EDOF).However, the results obtained for MIMO capacity were based on the assumption that the perfect channel state information (CSI) was available at least at the receiver end.In practical cases, perfect CSI is not achievable due to channel estimation errors.To function International Journal of Antennas and Propagation properly, the MIMO system has to estimate properties of the channel before decoding the received data.
Obtaining an accurate CSI can be accomplished using suitable channel estimation methods.The methods based on the use of training sequences, known as the trainingbased channel estimation methods, are the most popular.In [16,17], several training-based methods including least square (LS) method, scaled least square (SLS) method and minimum mean square error (MMSE) method have been investigated.It has been shown that the accuracy of the investigated training-based channel estimation methods is influenced by the transmitted signal to noise ratio (SNR) in the training mode, and the number of antenna elements at the transmitter and receiver.In particular, it has been pointed out that when the transmitted SNR and the number of antenna elements are fixed, the SLS and MMSE methods offer better performance than the LS method.In [18][19][20], it has been shown that if the training-based SLS or MMSE channel estimation methods are used to estimate CSI, their performance can be improved when the channels exhibit a higher spatial correlation level.According to [9][10][11][12][13][14][15], when the antenna element spacing is within 0.2λ and 0.4λ, the mutual coupling can decrease the spatial correlation level.Therefore in such cases, the accuracy of training-based channel estimation methods is compromised.From these considerations one can see a tradeoff between the accuracy of channel estimation and the resulting capacity when the MIMO system operates under mutual coupling conditions.When the correlation increases, the channel estimation error becomes smaller.However, at the same time the increased correlation is responsible for decreasing the effective degree of freedom thus the MIMO channel capacity.
In this paper, the overall effect of antenna mutual coupling on performance of an MIMO wireless system is investigated.The channel capacity that provides the overall measure of MIMO system performance is determined assuming an imperfect CSI is available at the receiver.A closed-form mathematical expression for the MIMO channel capacity that includes channel estimation errors is derived.The effects of mutual coupling on channel estimation and resulting capacity are investigated via computer simulations using a suitable MIMO channel model.
The rest of the paper is organized as follows.Section 2 introduces MIMO channel and its capacity.Also described in this section are the training-based MIMO channel estimation methods, expressions for MIMO capacity including channel estimation errors, and the Kronecker model [21,22] of MIMO channel.Section 3 describes the mathematical model for the antenna mutual coupling.Section 4 describes computer simulation results for the MIMO system capacity under the assumption of perfect CSI, as well as when an imperfect CSI due to estimation errors is available at the receiver.Section 5 concludes the paper.

MIMO System
2.1.MIMO Channel Capacity.Generally, a wireless communication channel is described by the relationship between the input signal S and output signal Y .Assuming that S is a set of signals, p(s) is the probability of S given s, the entropy or a measure of uncertainty of S can be written as, in which I(s) is the self-information which is the entropy contribution of an individual message; E{•} stands for the expected value.
Similarly, the entropy of Y given y can be expressed as, The joint entropy of S and Y is given as, The conditional entropy (or conditional uncertainty) of S given Y is expressed as, ( The conditional entropy gives the uncertainty about the channel input after the channel output has been observed.In this case, the uncertainty of the channel input can be expressed as H(X) − H(X | Y ).This is commonly known as the mutual information of the channel which is given by, ( The mutual information that can be optimized in terms of p(s) is independent of the channel.Because the wireless communication channel capacity is defined as the maximum mutual information, then the channel capacity can be determined from, In a narrow band MIMO wireless communication system, the channel is described by a complex matrix H.The relationship between the input signal S and the output signal Y can be expressed as, where N is the Gaussian noise with the property that R N = E{NN H } = σ 2 n .σ 2 n is the noise power.The S is transmitted random i.i.d input signal having the property that Q = E{SS H } and trace(Q) is equal to the total transmit power P s .The transmitted signal to noise ratio (SNR) is The upper bond mutual information of MIMO can be written as, Assuming the number of transmitting and receiving antennas is the same (M t equal to M r ) and the transmitted power is equally allocated to the transmitting antennas, (8) turns into The upper bound ergodic capacity can be written as, An alternative expression for the capacity can be obtained by decomposing the MIMO channel into n = min(M r , M t ) equivalent single input single output (SISO) subchannels.
The gains of each subchannel can be represented by eigenvalues of channel correlation matrix R H = HH H as given by( 11) where ρ i and λ i are transmit SNR and gain on ith subchannel, respectively.In this expression the overall power constraint is used expressed as n i ρ i = ρ and λ i is the ith eigenvalue of the channel correlation matrix.
If power is equally allocated to each subchannel, which is easy to realize in practice, then (11) can be rewritten in the form( 12)

Training-Based MIMO Channel Estimation.
For training based channel estimation methods, the relationship between the received signal and the training sequence is assumed as given by: where P represents the M t × L complex training matrix and L is the length of the training sequence.
The goal is to estimate the complex channel matrix H from the knowledge of Y and P. The transmitted power in the training mode is constrained by the following expression where P p is a given power constant and • 2 F stands for the Frobenius norm.The training matrices are assumed to be orthogonal [16,17] and the transmitted signal to noise ratio (SNR) in the training mode is set to ρ equal to P p /σ 2 n .Using the LS method, the estimated channel matrix can be written in the form [23] where {•} † stands for pseudoinverse and The mean square error (MSE) of the estimated channel matrix in the LS method is given as According to [16,17], the optimal training sequence satisfying ( 17) is given by The minimum value of Mean Square Error (MSE) for the LS method is obtained as From ( 18) one can see that the optimal performance of the LS estimator is influenced by the number of antenna elements at the transmitter and the receiver.However, the channel matrix (and thus CSI) has no effect on the minimum value of MSE.
The SLS method further reduces the estimation error that is obtained in the LS method.The improvement is given by the scaling factor γ: The minimized value of MSE for the SLS method can be written as [16,17] By substituting (18) into (20), the minimum MSE is given by expression, The estimated channel matrix obtained using the SLS method is given as [16,17], Similar derivations apply for the MMSE channel estimation method.Using this method, the estimated channel matrix can be obtained from, By denoting the channel estimation error as, the MSE for the MMSE method can be expressed as In (25), Q is a unitary eigenvector matrix of R H and Λ is a diagonal matrix with eigenvalues of R H .They are obtained by applying the eigenvalue decomposition to R H :

MIMO Channel Capacity Including Channel Estimation
Errors.In practical cases, H has to be replaced by the estimated channel matrix, which carries an estimation error as described in the last section.By assuming that the channel estimation error is defined as e and the estimated channel matrix is H given by The received signal can accordingly be written as The variance of each element in the estimation error matrix is given as If the channel matrix describes an uncorrelated Rayleigh channel, the elements in the estimation error matrix e are all uncorrelated identical independent distribution (i.i.d) with zero mean.In this case, the variance of each element in the error matrix shares the same value as σ 2 e ; the error matrix correlation is diagonal and holds the property Using (25), one finds The channel capacity of MIMO system with an imperfectly known H at the receiver is defined as the maximum mutual information between Y and S: If the transmitter does not have any knowledge of the channel, the mutual information in (32) can be written as, Because adding any dependent term on Y does not change the entropy [24], then in which u is the MMSE estimator given as Combining the above with (28), one obtains If one assumes the special case of M t equal to M r and the transmitted signal power being equally allocated to the transmitting antennas, (36) becomes in which p = P s /M t is the power allocated to the signal transmitted through each transmit antenna.The conditioning decreases the entropy and thus, Next, the following can be applied: and thus, For the case of S | H and (S − uY | Y , H)having a Gaussian distribution, (40) can be expressed as [24][25][26], The lower bound of the ergodic channel capacity is then given as An uncorrelated Rayleigh channel is ideal.In practical cases, an MIMO channel is usually correlated due to the propagation environment.Consequently, the estimation error matrix is correlated rather than i.i.d.In this case, the variances of elements in the error matrix are not equal to each other and the error correlation matrix is not diagonal.Therefore, (30) and (31) are not true anymore.Thus, (37) needs to be modified as,  44) indicates that for a fixed value of ρ, the capacity is a function of the estimated channel matrix H and the channel estimation error correlation R e .As a result, the channel matrix and the quality of channel estimation influence the MIMO capacity.

MIMO Channel Model.
In the undertaken investigations, the Kronecker channel model [21,22] is postulated to construct the channel matrix H.In this model, the transmitter and receiver correlation matrices are assumed to be separable and the channel matrix is represented by: where G H is the matrix including identical independent distributed (i.i.d) Gaussian entries with the zero mean and the unit variance, and R R and R T are the spatial correlation matrices at the receiver and transmitter, respectively.Here, it is assumed that the transmitting and receiving sides of MIMO system are equipped with vertically polarized wire dipole antennas.The scattering environment is represented by circles of uniformly distributed scattering objects surrounding the transmitting and receiving nodes.In this case, the spatial correlation matrix elements can be obtained using the Clark's model as given by, where J 0 stands for the zero-order Bessel function, κ is a wave number and d i j is the distance between elements i and j of the uniform array antenna.The required correlation matrices R T and R R can then be generated using (46) Having determined R T and R R , the channel matrix can be calculated using (45).

Mutual Coupling Effects
The mutual coupling in an array of collinear side-by-side wire dipoles can be modeled using the theory described in [27].Assuming the array is formed by M rt wire dipoles, the mutual matrix can be calculated using the following relationship involving the impedance matrix where Z A is the element impedance in isolation, for example, when the wire dipole is λ/2, its value is Z A = 73 + j42.5[Ω]; Z T is impedance of the receiver at each element chosen as the complex conjugate of Z A to obtain the impedance match.
The mutual impedance matrix Z is given by Note that this expression provides the circuit representation for mutual coupling in array antennas.It is valid for single mode antennas.The assumed here wire dipoles fall into this category.For a side-by-side array configuration of wire dipoles having length l equal to 0.5λ, the expressions for {Z mn } can be adapted from [27,28] and are rewritten here as, where κ is the wave number equal to 2π/λ, d h is the horizontal distance between the two dipole antenna elements.C i (u) and S i (u) are the cosine and sine integrals, respectively.They are given as, In the presence of mutual coupling, the channel matrices H appearing in expression (45) have to be replaced by new channel matrices as given by  The modified transmit and receive correlation matrices (including mutual coupling effects) can be expressed by Using (54), the expression shown in (45) can be rewritten as (55) and (56), respectively, Expressions ( 55) and (56) are the fundamental expressions that can be used to assess the performance of the introduced training-based channel estimation methods and the resulting channel capacity for the MIMO system operating under the condition of mutual coupling being present at transmitter and receiver array antennas.CSI is available at the receiver.In the undertaken computer simulations, both transmitter and receiver are assumed to be equipped with a 4-element uniform array antenna.The elements are assumed to be wire dipoles having length of 0.5λ.Equations ( 45) and ( 55) are used to model the MIMO channel.Both cases when mutual coupling effects are taken into account or disregarded are considered.The transmit SNR is fixed at 20 dB.In the figure, the blue curves represent the case of taking into account mutual coupling effect while the blue curves represent the case when the mutual coupling effect is not taken into account.From the results presented in Figure 1, one can see that when the antenna spacing is smaller than 0.4λ, the capacity with mutual coupling effect is larger than the one without taking mutual coupling effect into account.When the spacing is larger than 0.4λ, the channel capacity when mutual coupling effect is taken into account is smaller than when the mutual coupling effect is disregarded.The cross-point occurs at the antenna spacing equal to 0.4λ.These results show that when the antenna element spacing is within 0.2λ and 0.4λ, the mutual coupling helps to increase the channel capacity when a perfect CSI is assumed to be available at the receiver.

Simulation Results
Figure 2 presents the MSE performance of the MMSE and SLS channel estimation methods versus the transmitter antenna element spacing.One can see that MSE for the MMSE method is smaller than for the SLS method.Therefore the MMSE shows a superior performance.This agrees with the findings in publications [16,17].On the other hand, for both MMSE and SLS methods, an apparent difference in MSE occurs when mutual coupling effects are taken into account or disregarded when the transmitter antenna element spacing varies from 0.1λ to 0.5λ.It is apparent that MSE when mutual coupling effects are taken into account is worse than the one when the mutual coupling is disregarded (not taken effect into account).When the antenna element spacing is larger than 0.5λ, the difference disappears.This finding agrees with the prediction that when the antenna element spacing is between 0.2λ to 0.4λ the mutual coupling compromises the accuracy of MIMO channel estimation.

Imperfect CSI due to Channel Estimation
Errors.This time, computer simulations for MIMO channel capacity include both mutual coupling effects as well as channel estimation errors.MMSE method is assumed to be used to estimate the channel matrix (CSI).Again, both the transmitter and receiver are assumed to be equipped with 4-element uniform array antennas.The antenna elements are assumed to be wire dipoles having length of 0.5λ.The antenna spacing at the receiver is set to 1λ while the one at transmitter is assumed to be 0.2λ.
Figure 3(a) shows the results for the 4 × 4 MIMO channel capacity assuming perfect and imperfect CSI at the receiver.One can see that when the channel estimation error is taken into account, the channel capacity is lower than the one assuming a perfect CSI at the receiver.Also, at the transmitter antenna element spacing equal to 0.2λ, the channel capacity taking into account the mutual coupling effect is always higher than when the mutual coupling is neglected.Figure 3(b) shows the results for MSE of MMSE method for the 4 × 4 MIMO system when mutual coupling effect is included or disregarded.One can see that for both cases, the MSE decreases as SNR becomes larger.One can see that MSE becomes worse when mutual coupling effect is included.The results presented in Figures 3(a) and 3(b) show that although the mutual coupling undermines the channel estimation accuracy, it does not adversely affect the channel capacity.For the presented cases, the mutual coupling decreasing the spatial correlation level improves the signal transmission quality.Therefore, the worsening channel estimation has only an insignificant impact on capacity.

Conclusion
In this paper, the effect of mutual coupling on MIMO channel estimation accuracy and the resulting channel capacity has been investigated.The mathematical analysis and simulation results have shown that when the antenna element spacing at either transmitter or receiver is within 0.2λ and 0.4λ, the mutual coupling decreases the spatial correlation level and undermines the estimation accuracy of the MIMO channel.However, the reduced accuracy of channel estimation does not have a significant impact on the channel capacity.The carried out computer simulations have shown that the decreased spatial correlation improves the capacity of MIMO system despite having an adverse effect on the channel estimation.
Capacity taking MC into account Capacity without taking MC into account (b) Capacity versus receiver antenna spacing

Figure 1 :
Figure 1: A 4 × 4 MIMO channel capacity considering mutual coupling effect versus transmitter or receiver antenna spacing at SNR equal to 20 dB.(a) Capacity versus transmitter antenna spacing, receiver antenna spacing is set to 1λ.(b) Capacity versus receiver antenna spacing, transmitter antenna spacing is set to 1λ.

Figure 2 :
Figure 2: MSE of training-based channel estimation for a 4 × 4 MIMO considering mutual coupling effect versus transmitter or receiver antenna spacing at SNR equal to 25 dB.(a) MMSE MSE versus transmitter antenna spacing, receiver antenna spacing is set to 1λ.(b) SLS MSE versus transmitter antenna spacing, receiver antenna spacing is set to 1λ.

Figure 3 :
Figure 3: (a) Channel capacity considering MMSE estimation error for a 4 × 4 MIMO under mutual coupling effect versus SNR.(b) MMSE MSE performance for a 4 × 4 MIMO under mutual coupling effect versus SNR; transmitter antenna spacing is set to 0. 1λ, receiver antenna spacing is set to 1λ.