#### Abstract

A multiuser large-scale MIMO system with antenna correlation and mutual coupling is investigated in this paper. Based on the maximum signal-to-interference-plus-noise ratio (SINR) criteria, the optimal beamforming (BF) vector at the base station (BS) for each user is first obtained using statistical channel state information (CSI). Then, a closed-form expression for the achievable sum rate is derived in terms of a finite number of generalized Meijer-G functions, which is applicable to an arbitrary number of array elements and/or users, and provides an efficient means of evaluating the system performance. Finally, numerical results are provided to confirm the validity of the theoretical analysis and show the impact of various channel parameters on the system performance.

#### 1. Introduction

The pursuit of higher throughput and reliability has established multiple-input multiple-output (MIMO) technology an integral part of many advanced wireless communication systems. The application of MIMO in cellular architectures, including techniques such as multiuser MIMO (MU-MIMO), has gained significant research interest since it can offer spatial multiplexing gain even if each user only has a single antenna [1]. In such systems, a base station (BS) employed with multiple antennas simultaneously serves a number of cochannel users. In order to address the need for even greater capacity in future wireless networks, large-scale MIMO (LS-MIMO) or massive MIMO has been considered as a promising candidate to substantially improve spectral efficiency for next generation wireless systems. Unlike conventional small-scale MIMO systems (e.g., the LTE standard allows for only up to 8 antennas), massive MIMO suggests the use of hundreds of antennas at the BS, serving dozens of users at the same time [2–5].

In practical large-scale MIMO systems, however, a limited array aperture coupled with a large number of antennas may result in two major concerns: antenna correlation and mutual coupling (MC) effects. It is well-known that antenna correlation often exists in many practical environments due to insufficient antenna spacing and/or lack of scattering [6]. In addition, limited spacing between antenna elements in an array is also responsible for the mutual coupling, which can degrade signal transmission and reception due to antenna impedance mismatch [7]. The effects of antenna correlation and mutual coupling on traditional MIMO systems have been well-studied in [8–14], where the antenna arrays are commonly restrained less than 88 elements. Specifically, the authors in [8] applied random matrix theory to obtain analytical characterizations of the capacity of correlated multiantenna channels. In [9], the effects of fading correlations on MIMO systems were investigated by introducing the concept of sample-mean outage and information theoretic measures. The work in [10] analyzed the achievable rate of MIMO systems in the presence of mutual coupling and spatial correlation, which consisted of dipole antennas placed side-by-side in a linear pattern and in a very limited physical space. In [11], the impact of mutual coupling induced by two closely spaced minimum scattering antennas at the subscriber unit on MIMO systems was investigated, where both (de)correlation effects and variations of antenna gain resulting from coupling mechanisms were considered. The authors in [12] addressed the question of how coupling impacts bandwidth in the context of circular arrays. Moreover, the impact of mutual coupling on LTE-R MIMO capacity for antenna array configurations in high speed railway scenario was investigated in [13]. The authors in [14] provided a systematic review of the mutual coupling in MIMO systems, including the effects on performances of MIMO systems and various decoupling techniques.

By invoking a large number of antenna elements, massive MIMO or large-scale array has been viewed as a critical technique to overcome the combined effect of the mobile data traffic soaring trend and the favorable radio spectrum scarcity for the forthcoming communication systems which brings significant advantages including higher spectral efficiency, reduced radiated power consumption, and greater simplicity in signal processing. Considering both the antenna correlation and mutual coupling, the authors in [15] have investigated the performance of multiuser massive MIMO systems with uniform linear arrays (ULAs). The authors in [16] analyzed the performance of multiuser massive ULA MIMO Systems with mutual coupling calibration. Besides, the effect of mutual coupling in ULA large-scale antenna array with frugal calibration was studied in [17]. Compared with ULAs, uniform rectangular arrays (URAs) or uniform circular arrays (UCAs) are more space-efficient and offer numerous advantages [18]. Therefore, [19] has explored the impact of mutual coupling on uplink multiuser massive MIMO systems with a URA. Nevertheless, these aforementioned works assume perfect channel state information must be available at the BS, which is too restrictive in practice owing to fast fading, limited training sequences, and the heavy feedback burden.

Motivated by the aforementioned observations, this paper investigates the performance of an uplink multiuser large-scale MIMO system in the presence of both antenna correlation and mutual coupling. Particularly, the contribution of this work can be summarized as follows.(i)Considering both the antenna correlation and mutual coupling, we establish a general framework of multiuser massive MIMO systems with URA. Compared to the existing works limited to the perfect CSI cases, our work employs statistical CSI to evaluate the system performance more practically.(ii)Based on the statistical CSI at the BS, we first formulate the maximum signal-to-interference-plus-noise (SINR) based optimization problem and obtain the optimal beamforming vector for each user by using the Rayleigh quotient. Using the generalized Meijer-G functions, a closed-form expression for the achievable sum rate of the large-scale MIMO system is derived.(iii)Theoretical analysis is confirmed through comparison with simulation results, and the effect of antenna correlation and mutual coupling on the LS-MIMO system is investigated. The findings in this work provide a useful guidance for the system design and performance evaluation of massive MIMO systems with joint antenna correlation and mutual coupling in the presence of imperfect CSI.

The rest of this paper is organized as follows. The system under consideration is described in Section 2. In Section 3, the channel statistical properties used to evaluate its performance are provided. In Section 4, the closed-form expression for the achievable sum rate is derived. Simulation and numerical results are presented in Section 5. Finally, conclusions are drawn in Section 6.

#### 2. System Model

We consider an uplink multiuser large-scale MIMO system as shown in Figure 1, where a BS equipped with an -element URA receives signals from single antenna users. The users are assumed to be randomly located within a single cell according to a uniform distribution, and the size of the URA is assumed to be with being wavelength of the transmitted signal, and and are the length and width of the URA in wavelengths, respectively. The number of elements in the azimuth and elevation dimensions are and , respectively. Therefore, we have and and .

**(a) Uplink multiuser large-scale MIMO system**

**(b) Diagram of URA**

##### 2.1. Uplink Multiuser Large-Scale MIMO System

In this paper, we focus on the uplink of multiuser large-scale MIMO systems, where the received signal at the BS can be expressed aswhere denotes the channel matrix with being the channel vector between the -th user and BS, and with representing the transmit power of the -th user. The vector represents zero mean additive white Gaussian noise (AWGN) satisfying , and is the vector of the individual symbols simultaneously transmitted by the users.

Throughout this paper, all the wireless channels are assumed to experience cluster-based correlated fading, which is a widely adopted channel model for the outdoor propagation environment [6]. In this model, the channel vector between the -th user and BS can be written aswhere is the total number of multipath signals, and represents the fading coefficient of the -th path channel. The vector denotes the 2D array steering vector for a - element URA, which can be rewritten aswhere the -th element of the can be expressed as [20]where is the wave number, and , respectively, represent the azimuth and elevation angle of arrival of the -th path signal, with and indicating the mean of angle-of-arrival (AOA), and and are the corresponding angular spreads.

##### 2.2. Channel Modeling

With the above array model and assuming mutual coupling, the correlated channel vector of the -th user can be expressed as [6, 15, 19]where(i) denotes the channel correlation matrix (CCM) of each link, which is given by with the eigendecomposition of the CCM being given by where is a diagonal matrix with the eigenvalues arranged in a nonincreasing order, and are the corresponding eigenvectors.(ii) has independent identically distributed (i.i.d) entries following .

The mutual coupling matrix can be calculated by using the following relationship involving the impedance matrix as [15]where and denote the antenna impedance and load impedance, respectively. According to the URA in Figure 1, the mutual impedance matrix can be described byWhere, for a side-by-side array configuration of wire dipoles, can be determined using the classical induced electromotive force (EMF) method in [21].

##### 2.3. Beamforming Scheme with Statistical CSI

The received signals are separated into individual data streams by multiplying with the linear detector aswhere denotes the beamforming matrix. For the purpose of our performance analysis, we write the output signal of the -th user after beamforming aswhere and are the -th columns of the beamforming matrix and channel matrix , respectively. By considering both of the channel correlation and mutual coupling effects, the average output SINR can be expressed aswhere we have applied the fact thatand obtain

Proposition 1. *Suppose and are the two sized Hermitian matrices, and is positive definite; the following inequality holds [22]:The inequality is satisfied with equality at the optimum only whenwhere can be any nonzero constant. In the above equations, and denote the largest eigenvalue and the corresponding eigenvector of matrix , respectively.*

By converting the average SINR expression into the form of Rayleigh quotient in Proposition 1, we can obtain the optimal beamforming vector in terms of the maximum SINR criterion, namely,where and . Here, and represent the maximum eigenvalue and corresponding eigenvector of the matrix.

By substituting (19) into (11), after some algebraic manipulations, we havewhere , , with , representing the average SNR of each link.

#### 3. Channel Statistical Property

This section studies the statistical properties of each link, which will be used in the subsequent derivations. First, the probability density function (PDF) of is given by [23]To obtain the PDF of , we first haveThen, following methods similar to those in [24], can be derived asBased on (20), the PDF of is given byHence, by substituting (21) and (23) into (24), we haveBy using [25, eq. ] along with the identity [26, eq. ], the PDF of can be expressed aswhere denotes the Meijer-G function of a single variable [26].

#### 4. Achievable Sum Rate

The ergodic achievable rate of the -th user is given by [24, 27]Utilizing (26) in (27) and applying the Meijer-G function representation of [26, eq. ],we obtain the closed-form expression of aswhere represents the Meijer-G functions of two variables [28]. In deriving (29), we have applied [29, eq. 2.6.2)]. Hence, the total achievable sum rate for all the users can be obtained by where the exact results can be accordingly obtained from (29) with algebraic manipulations.

#### 5. Numerical Results

In this section, we investigate the validity of the performance analysis and the effect of mutual coupling and antenna correlation on the LS-MIMO system through numerical simulations. Here, without loss of generality, the simulation parameters are chosen as follows: the ratio between the length and width of the URA is , element distance , and , [13], and the transmit SNR is 10dB. In addition, the users are uniformly distributed in both the azimuth and elevation directions, and the number of multipath components per user is , the number of users is , and the Monte Carlo simulations are obtained by averaging over channel realizations.

Figure 2 illustrates the achievable sum rate of the uplink multiuser large-scale MIMO system versus the number of antennas for different AOAs and angular spreads. Here, the AOAs are assumed to satisfy or , and the angular spreads are set as or , which denote the low and high antenna correlation, respectively. As can be observed, the theoretical results obtained in (30) agree well with the Monte Carlo simulations, which validates our derived sum rate expressions. Meanwhile, it can also be found that the achievable sum rate decreases with the increase in angular spread, implying that the performance of LS-MIMO can be improved with more correlated channels, when statistical CSI is used to fulfill optimal precoding.

By considering a fixed array with size of , Figure 3 plots the achievable sum rate versus interelement spacing for different angular spreads with . The numbers of antennas in each row and column of the array are calculated as and . The fluctuation of curves is due to the up-and-down of mutual impendence when antenna distance increases [30]. It is observed that when only the antenna correlation is taken into account (Correlation Only), an increase in the element spacing, which corresponds to a decrease in the number of antennas and the antenna correlation, leads to a decrease in the achievable sum rate. On the other hand, for the cases where both antenna correlation and mutual coupling are considered (Correlation and Coupling), the maximal achievable sum rate is obtained when . Although more antenna elements can be deployed with associated higher antenna correlation in the presence of smaller interelement spacing, the mutual coupling becomes dominant factor and reduces the sum rate of the system. In addition, the number of antennas and channel correlation are reduced with an increase in the interelement spacing with the array size fixed. With these two competing factors, is the best choice in terms of the achievable sum rate.

#### 6. Conclusions

In this paper, we have investigated the achievable sum rate of a multiuser large-scale MIMO system taking into account mutual coupling and antenna correlation. Specifically, with optimal SINR-based beamforming at the BS, an exact closed-form expression for the system sum rate is derived in terms of Meijer-G functions of one or two variables, which provides an efficient approach to indicate the impact of various key parameters on the system performance. Simulations validate the theoretical results and also show the favorable conditions for achievable sum rate when both antenna correlation and mutual coupling are considered.

#### Data Availability

The MATLAB code data used to support the findings of this study were supplied by NSFC 61471392 under license and so cannot be made freely available. Requests for access to these data should be made to Prof. Yusheng Li.

#### Conflicts of Interest

The authors declare that they have no conflicts of interest.

#### Acknowledgments

This work is supported by the Research Project of NUDT under grant ZK18-02-11 and the National Natural Science Foundation of China under Grant 61471392.