Abstract
Multicell joint processing has originated from informationtheoretic principles as a means of reaching the fundamental capacity limits of cellular networks. However, global multicell joint decoding is highly complex and in practice clusters of cooperating Base Stations constitute a more realistic scenario. In this direction, the mitigation of intercluster interference rises as a critical factor towards achieving the promised throughput gains. In this paper, two intercluster interference mitigation techniques are investigated and compared, namely interference alignment and resource division multiple access. The cases of global multicell joint processing and cochannel interference allowance are also considered as an upper and lower bound to the interference alignment scheme, respectively. Each case is modelled and analyzed using the percell ergodic sumrate throughput as a figure of merit. In this process, the asymptotic eigenvalue distribution of the channel covariance matrices is analytically derived based on freeprobabilistic arguments in order to quantify the sumrate throughput. Using numerical results, it is established that resource division multiple access is preferable for dense cellular systems, while cochannel interference allowance is advantageous for highly sparse cellular systems. Interference alignment provides superior performance for average to sparse cellular systems on the expense of higher complexity.
1 Introduction
Currently cellular networks carry the main bulk of wireless traffic and as a result they risk being saturated considering the ever increasing traffic imposed by internet data services. In this context, the academic community in collaboration with industry and standardization bodies have been investigating innovative network architectures and communication techniques which can overcome the interferencelimited nature of cellular systems. The paradigm of multicell joint processing has risen as a promising way of overcoming those limitations and has since gained increasing momentum which lead from theoretical research to testbed implementations [1]. Furthermore, the recent inclusion of CoMP (Coordinated Multiple Point) techniques in LTEadvanced [2] serves as a reinforcement of the latter statement.
Multicell joint processing is based on the idea that signal processing does not take place at individual base stations (BSs), but at a central processor which can jointly serve the user terminals (UTs) of multiple cells through the spatially distributed BSs. It should be noted that the main concept of multicell joint processing is closely connected to the rationale behind Network MIMO and Distributed Antenna Systems (DAS) and those three terms are often utilized interchangeably in the literature. According to the global multicell joint processing, all the BSs of a large cellular system are assumed to be interconnected to a single central processor through an extended backhaul. However, the computational requirements of such a processor and the large investment needed for backhaul links have hindered its realization. On the other hand, clustered multicell joint processing utilizes multiple signal processors in order to form BS clusters of limited size, but this localized cooperation introduces intercluster interference into the system, which has to be mitigated in order to harvest the full potential of multicell joint processing. In this direction, reuse of time or frequency channel resources (resource division multiple access) could provide the necessary spatial separation amongst clusters, an approach which basically mimics the principles of the traditional cellular paradigm only on a cluster scale. Another alternative would be to simply tolerate intercluster signals as cochannel interference, but obviously this scheme becomes problematic in highly dense systems. Taking all this into account, the current paper considers the uplink of a clustered multicell joint decoding (MJD) system and proposes a new communication strategy for mitigating intercluster interference using interference alignment (IA). More specifically, the main contributions herein are:
1. the channel modelling of a clustered MJD system with IA as intercluster interference mitigation technique,
2. the analytical derivation of the ergodic throughput based on free probabilistic arguments in the Rtransform domain,
3. the analytical comparison with the upper bound of global MJD, the Resource Division Multiple Access (RDMA) scheme and the lower bound of clustered MJD with Cochannel Interference allowance (CI),
4. the comparison of the derived closedform expressions with Monte Carlo simulations and the performance evaluation using numerical results.
The remainder of this paper is structured as follows: Section 2 reviews in detail prior work in the areas of clustered MJD and IA. Section 3 describes the channel modelling, free probability derivations and throughput results for the following cases: (a) global MJD, (b) IA, (c) RDMA and (d) CI. Section 4 displays the accuracy of the analysis by comparing to Monte Carlo simulations and evaluates the effect of various system parameters in the throughput performance of clustered MJD. Section 6 concludes the paper.
1.1 Notation
Throughout the formulations of this paper, denotes expectation, (·)^{H }denotes the conjugate matrix transpose, (.)^{T }denotes the matrix transpose, ⊙ denotes the Hadamard product and ⊗ denotes the Kronecker product. The Frobenius norm of a matrix or vector is denoted by · and the delta function by δ(·). I_{n }denotes an n × n identity matrix, an n × m matrix of ones, 0 a zero matrix and denotes n × m Gaussian matrix with entries drawn form a distribution. The figure of merit analyzed and compared throughout this paper is the ergodic percell sumrate throughput.^{a}
2 Related work
2.1 Multicell joint decoding
This section reviews the literature on MJD systems by describing the evolution of global MJD models and subsequently focusing on clustered MJD approaches.
2.1.1 Global MJD
It was almost three decades ago when the paradigm of global MJD was initially proposed in two seminal papers [3,4], promising large capacity enhancements. The main idea behind global MJD is the existence of a central processor (a.k.a. "hyperreceiver") which is interconnected to all the BSs through a backhaul of wideband, delayless and errorfree links. The central processor is assumed to have perfect channel state information (CSI) about all the wireless links of the system. The optimal communication strategy is superposition coding at the UTs and successive interference cancellation at the central processor. As a result, the central processor is able to jointly decode all the UTs of the system, rendering the concept of intercell interference void.
Since then, the initial results were extended and modified by the research community for more practical propagation environments, transmission techniques and backhaul infrastructures in an attempt to more accurately quantify the performance gain. More specifically, it was demonstrated in [5] that Rayleigh fading promotes multiuser diversity which is beneficial for the ergodic capacity performance. Subsequently, realistic pathloss models and user distributions were investigated in [6,7] providing closedform ergodic capacity expressions based on the cell size, path loss exponent and geographical distribution of UTs. The beneficial effect of MIMO links was established in [8,9], where a linear scaling of the ergodic percell sumrate capacity with the number of BS antennas was shown. However, correlation between multiple antennas has an adverse effect as shown in [10], especially when correlation affects the BS side. Imperfect backhaul connectivity has also a negative effect on the capacity performance as quantified in [11]. MJD has been also considered in combination with DSCDMA [12], where chips act as multiple dimensions. Finally, linear MMSE filtering [13,14] followed by singleuser decoding has been considered as an alternative to the optimal multiuser decoder which requires computationallycomplex successive interference cancellation.
2.1.2 Clustered MJD
Clustered MJD is based on forming groups of M adjacent BSs (clusters) interconnected to a cluster processor. As a result, it can be seen as an intermediate state between traditional cellular systems (M = 1) and global MJD (M = ∞). The advantage of clustered MJD lies on the fact that both the size of the backhaul network and the number of UTs to be jointly processed decrease. The benefit is twofold; first, the extent of the backhaul network is reduced and second, the computational requirements of MJD (which depend on the number of UTs) are lower. The disadvantage is that the sumrate capacity performance is degraded by intercluster interference, especially affecting the individual rates of clusteredge UTs. This impairment can be tackled using a number of techniques as described here. The simplest approach is to just treat it as cochannel interference and evaluate its effect on the system capacity as in [15]. An alternative would be to use RDMA, namely to split the time or frequency resources into orthogonal parts dedicated to clusteredge cells [16]. This approach eliminates intercluster interference but at the same time limits the available degrees of freedom. In DSCDMA MJD systems, knowledge of the interfering codebooks has been also used to mitigate intercluster interference [12]. Finally, antenna selection schemes were investigated as a simple way of reducing the number of intercluster interferers [17].
2.2 Interference alignment
This section reviews the basic principles of IA and subsequently describes existing applications of IA on cellular networks.
2.2.1 IA preliminaries
IA has been shown to achieve the degrees of freedom (dofs) for a range of interference channels [1820]. Its principle is based on aligning the interference on a signal subspace with respect to the nonintended receiver, so that it can be easily filtered out by sacrificing some signal dimensions. The advantage is that this alignment does not affect the randomness of the signals and the available dimensions with respect to the intended receiver. The disadvantage is that the filtering at the nonintended receiver removes the signal energy in the interference subspace and reduces the achievable rate. The fundamental assumptions which render IA feasible are that there are multiple available dimensions (space, frequency, time or code) and that the transmitter is aware of the CSI towards the nonintended receiver. The exact number of needed dimensions and the precoding vectors to achieve IA are rather cumbersome to compute, but a number of approaches have been presented in the literature towards this end [2123].
2.2.2 IA and cellular networks
IA has been also investigated in the context of cellular networks, showing that it can effectively suppress cochannel interference [23,24]. More specifically, the downlink of an OFDMA cellular network with clustered BS cooperation is considered in [25], where IA is employed to suppress intracluster interference while intercluster interference has to be tolerated as noise. Using simulations, it is shown therein that even with unit multiplexing gain the throughput performance is increased compared to a frequency reuse scheme, especially for the clustercentre UTs. In a similar setting, the authors in [26] propose an IAbased resource allocation scheme which jointly optimizes the frequencydomain precoding, subcarrier user selection, and power allocation on the downlink of coordinated multicell OFDMA systems. In addition, authors in [24] consider the uplink of a limitedsize cellular system without BS cooperation, showing that the interferencefree dofs can be achieved as the number of UTs grows. Employing IA with unit multiplexing gain towards the nonintended BSs, they study the effect of multipath channels and singlepath channels with propagation delay. Furthermore, the concept of decomposable channel is employed to enable a modified scheme called subspace IA, which is able to simultaneously align interference towards multiple nonintended receivers over a multidimensional space. Finally, the effect of limited feedback on cellular IA schemes has been investigated and quantified in [25,27].
3 Channel model and throughput analysis
In this paper, the considered system comprises a modified version of Wyner's linear cellular array [4,12,28], which has been used extensively as a tractable model for studying MJD scenarios [29]. In the modified model studied herein, MJD is possible for clusters of M adjacent BSs while the focus is on the uplink. Unlike [23,24], IA is employed herein to mitigate intercluster interference between clusteredge cells. Let us assume that K UTs are positioned between each pair of neighboring BSs with path loss coefficients 1 and α, respectively (Figure 1). All BSs and UTs are equipped with n = K + 1 antennas ^{b }[10]. to enable IA over the multiple spatial dimensions for the clustered UTs. In this setting, four scenarios of intercluster interference are considered, namely global MJD, IA, RDMA and CI. It should be noted that only clusteredge UTs employ interference mitigation techniques, while UTs in the interior of the cluster use the optimal wideband transmission scheme with superposition coding as in [5]. Successive interference cancellation is employed in each cluster processor in order to recover the UT signals. Furthermore, each cluster processor has full CSI for all the wireless links in its coverage area. The following subsections explain the mode of operation for each approach and describe the analytical derivation of the percell sumrate throughput.
Figure 1. Graphical representation of the considered cellular system modelled as a modified version of Wyner's model. K UTs are positioned between each pair of neighboring BSs with path loss coefficients 1 and α respectively. All BSs and UTs are equipped with n = K + 1 antennas. The UTs positioned within the box shall be jointly processed. The red links denote intercluster interference.
3.1 Global multicell joint decoding
In global MJD, a central processor is able to jointly decode the signals received by neighboring clusters and, therefore, no intercluster interference takes place. In other words, the entire cellular system can be assumed to be comprised of a single extensive cluster. As it can be seen, this case serves as an upper bound to the IA case. The received n × 1 symbol vector y_{i }at any random BS can be expressed as follows:
where the n × 1 vector z denotes AWGN with and . The K n × 1 vector x_{i }denotes the transmitted symbol vector of the ith UT group with where γ is the transmit Signal to Noise ratio per UT antenna. The n × Kn channel matrix includes the flat fading coefficients of the ith UT group towards the ith BS modelled as independent identically distributed (i.i.d.) complex circularly symmetric (c.c.s.) random variables. Similarly, the term αG_{i, i+1}(t)x_{i+1}(t) represents the received signal at the ith BS originating from the UTs of the neighboring cell indexed i + 1. The scaling factor α < 1 models the amount of received intercell interference which depends on the path loss model and the density of the cellular system^{c}. Another intuitive description of the α factor is that it models the power imbalance between intracell and intercell signals.
Assuming a memoryless channel, the system channel model can be written in a vectorial form as follows:
where the aggregate channel matrix has dimensions Mn × (M + 1)Kn and can be modelled as:
with being a blockToeplitz matrix and . In addition, is a M × M + 1 Toeplitz matrix structured as follows:
Assuming no CSI at the UTs, the percell capacity is given by the MIMO multiple access (MAC) channel capacity:
Theorem 3.1. In the global MJD case, the percell capacity for asymptotically large n converges almost surely (a.s.) to the MarcenkoPastur (MP) law with appropriate scaling [6,10]:
Proof. For the sake of completeness and to facilitate latter derivations, an outline of the proof in [6,10] is provided here. The derivation of this expression is based on an asymptotic analysis in the number of antennas n → ∞:
where λ_{i }(X) and denote the eigenvalues and the asymptotic eigenvalue probability distribution function (a.e.p.d.f.) of matrix X respectively and denotes the Shannon transform of X with scalar parameter x. It should be noted that denotes the total UT transmit power normalized by the receiver noise power^{d}. The last step of the derivation is based on unit rank matrices decomposition and analysis on the Rtransform domain, as presented in [6,10]. The scaling factor
is the Frobenius norm of the Σ matrix normalized by the matrix dimensions and
where step (a) follows from [10, Eq.(34)]. □
3.2 Interference alignment
In order to evaluate the effect of IA as an intercluster interference mitigation technique, a simple precoding scheme is assumed for the clusteredge UT groups, inspired by [24]. Let us assume a n × 1 unit norm reference vector v with v^{2 }= n and
where y_{1 }and y_{M }represent the received signal vectors at the first and last BS of the cluster, respectively. The first UT group has to align its input x_{1 }towards the nonintended BS of the cluster on the left (see Figure A), while the Mth BS has to filter our the aligned interference coming from the M + 1th UT group which belongs to the cluster on the right. These two strategies are described in detail in the following subsections:
3.2.1 Aligned interference filtering
The objective is to suppress the term αG_{M, M+1}x_{M+1 }which represents intercluster interference. It should be noted that UTs of the M + 1th cell are assumed to have perfect CSI about the channel coefficients G_{M, M+1}. Let us also assume that and represent the transmitter vector and channel matrix of the jth UT in the ith group towards the BS. In this context, the following precoding scheme is employed to align interference:
where v_{j }= vv_{j }is a scaled version of v which satisfies the input power constraint . This precoding results in unit multiplexing gain and is by no means the optimal IA scheme^{e }[22] provide conditions for classifying a scenario as proper or improper, a property which is shown to be connected to feasibility., but it serves as a tractable way of evaluating the IA performance [23,24]. the feasibility of IA. Following this approach, the intercluster interference can be expressed as:
It can be easily seen that interference has been aligned across the reference vector and it can be removed using a K× n zeroforcing filter Q designed so that Q is a truncated unitary matrix [19] and Qv = 0. After filtering, the received signal at the Mth BS can be expressed as:
Assuming that the system operates in highSNR regime and is therefore interference limited, the effect of the AWGN noise colouring can be ignored, namely .
Lemma 3.1. The Shannon transform of the covariance matrix of QG_{M,M }is equivalent to that of a K × K Gaussian matrix G_{K×K}.
Proof. Using the property det(I + γAB) = det(I + γBA), it can be written that:
The K × n truncated unitary matrix Q has K unit singular values and therefore the matrix product Q^{H}Q has K unit eigenvalues and a zero eigenvalue. Applying eigenvalue decomposition on Q^{H}Q, the left and right eigenvectors can be absorbed by the isotropic Gaussian matrices and G_{M,M }respectively, while the zero eigenvalue removes one of the n dimensions. Using the definition of Shannon transform [30], Eq. (15) yields
□
Based on this lemma and for the purposes of the analysis, QG_{M,M }is replaced by G_{K × K }in the equivalent channel matrix.
3.2.2 Interference alignment
The Mth BS has filtered out incoming interference from the cluster on the right (Figure 1), but outgoing intercluster interference should be also aligned to complete the analysis. This affects the first UT group which should align its interference towards the Mth BS of the cluster on the left (Figure 1). Following the same precoding scheme and using Eq. (10)
where represents the fading coefficients of the jth UT of the first group towards the Mth BS of the neighboring cluster on the left. Since the exact eigenvalue distribution of the matrix product is not straightforward to derive, for the purposes of rate analysis it is approximated by a Gaussian vector with unit variance. This approximation implies that IA precoding does not affect the statistics of the equivalent channel towards the intended BS.
3.2.3 Equivalent channel matrix
To summarize, IA has the following effects on the channel matrix H used for the case of global MJD. The intercluster interference originating from the M + 1th UT group is filtered out and thus Kn vertical dimensions are lost. During this process, one horizontal dimension of the Mth BS is also filtered out, since it contains the aligned interference from the M + 1th UT group. Finally, the first UT group has to precode in order to align its interference towards the Mth BS of neighboring cluster and as a result only K out of Kn dimensions are preserved. The resulting channel matrix can be described as follows:
where _{}and
with being a n × (M  1)Kn + K matrix^{f}, being a (M  2)n × (M  1)Kn + K matrix and being a n  1 × (M  1)Kn + K matrix^{g}.
Since all intercluster interference has been filtered out and the effect of filter Q has been already incorporated in the structure of H_{IA}, the percell throughput in the IA case is still given by the MIMO MAC expression:
Theorem 3.2. In the IA case, the percell throughput can be derived from the Rtransform of the a.e.p.d.f. of matrix .
Proof. Following an asymptotic analysis where n→ ∞:
The a.e.p.d.f. of ^{}is obtained by determining the imaginary part of the Stieltjes transform for real arguments
considering that the Stieltjes transform is derived from the Rtransform [31] as follows
□
Theorem 3.3. The Rtransform of the a.e.p.d.f. of matrix is given by:
with k, β, q parameters given by:
Proof. Based on Eq.(19), the matrix _{}can be decomposed as the following sum:
where H_{1 }= Σ_{1 }⊙G_{n×(M1)Kn+K}, H_{2 }= Σ_{2 }⊙ G_{(M2)n×(M1)Kn+K }and H_{3 }= Σ_{3 }⊙ G_{n1×(M1)Kn+K}. Using the property of free additive convolution [30] and Theorem A.1 in Appendix A, Eq. (24) holds in the Rtransform domain. □
3.3 Resource division multiple access
RDMA entails that the available time or frequency resources are divided into two orthogonal parts assigned to clusteredge cells in order to eliminate intercluster interference [[16], Efficient isolation scheme]. More specifically, for the first part clusteredge UTs are inactive and the farright clusteredge _{BS }is active. For the second part, clusteredge UTs are active and the farright clusteredge _{BS }is inactive. While the available channel resources are cut in half for clusteredge UTs, double the power can be transmitted during the second part orthogonal part to ensure a fair comparison amongst various mitigation schemes. The channel modelling is similar to the one in global _{MJD }case (Eq. (1)), although in this case the throughput is analyzed separately for each orthogonal part and subsequently averaged. Assuming no _{CSI }at the UTs, the percell throughput in the RDMA case is given by:
where and denote the capacities for the first and second orthogonal part respectively. For the first part, the cluster processor receives signals from (M  1)K UTs through all M BSs and the resulting Mn × (M  1)Kn channel matrix is structured as follows:
Theorem 3.4. For the first part of the RDMA case, the percell throughput C_{RD1 }can be derived from the Rtransform of the a.e.p.d.f of matrix , where:
Proof. Following an asymptotic analysis where n→ ∞:
Using the matrix decomposition of Eq. (27) and free additive convolution [30]:
Eq. (28) follows from Eq. (42) with and theorem A.1. □
For the second part, the cluster processor receives signals from MK UTs through M  1 BSs and the resulting (M  1)n ×MKn channel matrix is structured as follows:
where the factor 2 is due to the doubling of the transmitted power.
Theorem 3.5. For the second part of the RDMA case, the percell throughput can be derived from the Rtransform of the a.e.p.d.f. of matrix , where:
Proof. Following an asymptotic analysis where n → ∞:
The rest of this proof follows the steps of Theorem 3.4. □
3.4 Cochannel interference allowance
CI is considered as a worst case scenario where no signal processing is performed in order to mitigate intercluster interference and thus interference is treated as additional noise [15]. As it can be seen, this case serves as a lower bound to the IA case. The channel modelling is identical with the one in global MJD case (Eq. (1)), although in this case the clusteredge UT group contribution αG_{M, M}_{+1}(t)x_{M}_{+1}(t) is considered as interference. As a result, the interference channel matrix can be expressed as:
Assuming no CSI at the UTs, the percell throughput in the CI case is given by [15,3234]:
where C_{I }denotes the throughput of the interfering UT group normalized by the cluster size:
Theorem 3.6. In the CI case, the percell throughput converges almost surely (a.s.) to a difference of two scaled versions of the the MP law:
Proof. Following an asymptotic analysis in the number of antennas n n →∞:
Eq. (37) follows from Eq. (35), (38) and Theorem 3.1. □
3.5 Degrees of freedom
This section focuses on comparing the degrees of freedom for each of the considered cases. The degrees of freedom determine the number of independent signal dimensions in the high SNR regime [35] and it is also known as prelog or multiplexing gain in the literature. It is a useful metric in cases where interference is the main impairment and AWGN can be considered unimportant.
Theorem 3.7. The degrees of freedom per BS antenna for the global MJD, I A, RDM A and CI cases are given by:
Proof. Eq. (39) can be derived straightforwardly by counting the receive dimensions of the equivalent channel matrices (Eq. (3) for global MJD, Eq. (18) for IA, Eqs. (27) and (31) for RDMA, Eq. (34) for CI) and normalizing by the number of BS antennas. □
Lemma 3.2. The following inequalities apply for the dofs of eq. (39):
Remark 3.1. It can be observed that d_{IA }= d_{RD }only for single UT per cell equipped with two antennas (K = 1, n = 2). For all other cases, d_{IA }> d_{RD}. Furthermore, it is worth noting that when the number of UTs K and antennas n grows to infinity, lim_{K,n }→ ∞ d_{IA }= d_{MJD }which entails a multiuser gain. However, in practice the number of served UTs is limited by the number of antennas (n = K+1) which can be supported at the BS and more importantly at UTside due to size limitations.
3.6 Complexity considerations
This paragraph discusses the complexity of each scheme in terms of decoding processing and required CSI. In general, the complexity of MJD is exponential with the number of users [36] and full CSI is required at the central processor for all users which are to be decoded. This implies that global MJD is highly complex since all system users have to be processed at a single point. On the other hand, clustering approaches reduce the number of jointlyprocessed users and as a result complexity. Furthermore, CI is the least complex since no action is taken to mitigate intercluster interference. RDMA has an equivalent receiver complexity with CI, but in addition it requires coordination between adjacent clusters in terms of splitting the resources. For example, time division would require intercluster synchronization, while frequency division could be even static. Finally, IA is the most complex since CSI towards the nonintended BS is also needed at the transmitter in order to align the interference. Subsequently, additional processing is needed at the receiver side to filter out the aligned interference.
4 Numerical results
This section presents a number of numerical results in order to illustrate the accuracy of the derived analytical expressions for finite dimensions and evaluate the performance of the aforementioned interference mitigation schemes. In the following figures, points represent values calculated through Monte Carlo simulations, while lines refer to curves evaluated based on the analytical expressions of section 3. More specifically, the simulations are performed by generating 10^{3 }instances of random Gaussian matrices, each one representing a single fading realization of the system. In addition, the variance profile matrices are constructed deterministically based on the considered α factors and used to shape the variance of the i.i.d. c.c.s. elements. Subsequently, the percell capacities are evaluated by averaging over the system realizations using: (a) Eq. (5) for global MJD, (b) Eq. (10)(14) and (20) for IA, (c) Eq. (26) for RDMA, (d) Eq. (35),(36) for CI. In parallel, the analytical curves are evaluated based on: (a) theorem 3.1 for global MJD, (b) theorems 3.2 and 3.3 for IA, (c) theorems 3.4 and 3.5 for RDMA, (d) theorems 3.1 and 3.6 for CI. Table 1 presents an overview of the parameter values and ranges used for producing the numerical results of the figures.
Table 1. Parameters for throughput results
Firstly, Figure 2 depicts the percell throughput versus the cluster size M for medium α factors. It should be noted that the α factor combines the effects of cell size and path loss exponent as explained in [37]. As expected the performance of global MJD does not depend on the cluster size, since it is supposed to be infinite. For all interference mitigation techniques, it can be seen that the penalty due to the clustering diminishes as the cluster size increases. Similar conclusions can be derived by plotting the degrees of freedom versus the cluster size M (Figure 3). In addition, it can be observed that the IA dofs approach the global MJD dofs as the number of UTs and antennas increases. Subsequently, Figure 4 depicts the percell throughput versus the α factor. For high α factors, RDMA performance converges to IA, whereas for low α factors RDMA performance degrades. It should be also noted that while the performance of global MJD and RDMA increase monotonically with α, the performances of IA and cochannel interference degrade for medium α factors. Finally, Figure 5 depicts the percell throughput versus the number of UTs per cell K. It should be noted that the number of antennas per UT n scale jointly with K. Based on this observation, a superlinear scaling of the performance can be observed, resulting primarily from the increase of spatial dimensions (more antennas) and secondarily from the increase of the system power (more UTs). As it can be seen, the slope of the linear scaling is affected by the selected interference mitigation technique.
Figure 2. Percell throughput scaling versus the cluster size M. The performance of global MJD does not depend on the cluster size, while for all interference mitigation techniques, the penalty due to the clustering diminishes as the cluster size increases. Parameters: α = 0.5, K = 5, n = 4, γ = 20 dB.
Figure 3. Degrees of freedom versus the cluster size M. The IA dofs approach the global MJD dofs as the number of UTs and antennas increases. Parameters: α = 0.5, K = 5, n = 4, γ = 20 dB.
Figure 4. Percell throughput scaling versus the α factor. For high α factors, RDMA performance converges to interference alignment, whereas for low α factors RDMA performance degrades even beyond the cochannel interference bound. While the performance of global MJD and RDMA increase monotonically with α, the performances of interference alignment and cochannel interference degrade for medium α factors. Parameters: M = 4, K = 5, n = 4, γ = 20 dB.
Figure 5. Percell throughput scaling versus the number of UTs per cell K. A superlinear scaling of the performance can be observed, resulting primarily from the increase of spatial dimensions (more antennas) and secondarily from the increase of the system power (more UTs). The slope of the linear scaling is affected by the selected interference mitigation technique. Parameters: M = 4, α = 0.5, γ = 20 dB.
6 Conclusion
In this paper, various techniques for mitigating intercluster interference in clustered MJD were investigated. The case of global MJD was initially considered as an upper bound, serving in evaluating the degradation due to intercluster interference. Subsequently, the IA scheme was analyzed by deriving the asymptotic eigenvalue distribution of the channel covariance matrix using freeprobabilistic arguments. In addition, the RDMA scheme was studied as a low complexity method for mitigating intercluster interference. Finally, the CI was considered as a worstcase scenario where no interference mitigation techniques is employed. Based on these investigations it was established that for dense cellular systems the RDMA scheme should be used as the best compromise between complexity and performance. For average to sparse cellular systems which is the usual regime in macrocell deployments, IA should be employed when the additional complexity and availability of CSI at transmitter side can be afforded. Alternatively, CI could be preferred especially for highly sparse cellular systems.
A Proof of theorem
Theorem A.1. Let A = [0 B 0] be the concatenation of the varianceprofiled Gaussian matrix B = C ⊙ G and a number of zero columns. Let also k be the ratio of nonzero to total columns of A, β be the ratio of horizontal to vertical dimensions of B and q the Frobenius norm of C normalized by the matrix dimensions. The Rtransform of A^{H }A is given by:
Proof. Let B = C ⊙ G_{n × m }be a varianceprofiled Gaussian matrix with β = m/n and q = C^{2}/nm. According to [10], the Rtransform of is given by:
Using eq. (23), the Stieltjes transform of can be expressed as:
Matrix has identical eigenvalues to plus a number of zero eigevalues with 0 < k < 1 defined as the ratio of nonzero eigenvalues over the total number of eigenvalues. As a result, the a.e.p.d.f. of ^{}can be written as:
Using the definition of the Stieltjes transform [30]:
and employing eq. (23), the proof is complete. □
5 Competing interests
The authors declare that they have no competing interests.
Notes
^{a}The term throughput is used instead of capacity since the described techniques are suboptimal in the informationtheoretic sense and lead to achievable sumrates except for MJD which leads to MIMO MAC capacity.
^{b}The multiple antennas are assumed to be uncorrelated although the analytical results can be extended in the correlated case based on the principles described in
^{c}For more details on the modelling of the α parameter, the reader is referred to [37].
^{d}For the purposes of the analysis the variable is kept finite as the number of antennas Mn grows large, so that the system power does not grow to infinity.
^{e}Depending on the signal dimensions and the channel coefficients, more than one degree of freedom per user could be achieved. The feasibility of higher multiplexing gain has been studied in [21,22]. More specifically, authors in [21] provide an algorithm which determines the achievable multiplexing gain by minimizing the interference leakage, while authors in
^{f}The structure of the first block of Σ_{1 }originates in the Gaussian approximation of .
^{g}The structure of the last block of Σ_{3 }is based on Lemma 3.1.
Acknowledgements
This work was partially supported by the National Research Fund, Luxembourg under the CORE project "CO^{2}SAT: Cooperative and Cognitive Architectures for Satellite Networks".
References

P Marsch, G Fettweis, Coordinated MultiPoint in Mobile Communications From Theory to Practice (Cambridge University Press, 2011)

3GPP:, Simulation Assumptions for Evaluation of UL CoMP.. R1110803 (2011) [http://www.3gpp.org/ftp/tsg_ran/wg1_rl1/TSGR1_64/Docs/R1110803.zip] webcite

SV Hanly, PA Whiting, Informationtheoretic capacity of multireceiver networks. Telecommun Syst 1, 1–42 (1993). Publisher Full Text

A Wyner, Shannontheoretic approach to a Gaussian cellular multipleaccess channel. IEEE Trans Inf Theory 40(6), 1713–1727 (1994). Publisher Full Text

O Somekh, S Shamai, Shannontheoretic approach to a Gaussian cellular multipleaccess channel with fading. IEEE Trans Inf Theory 46(4), 1401–1425 (2000). Publisher Full Text

S Chatzinotas, M Imran, C Tzaras, On the capacity of variable density cellular systems under multicell decoding. IEEE Commun Lett 12(7), 496–498 (2008)

S Chatzinotas, MA Imran, C Tzaras, Optimal information theoretic capacity of the planar cellular uplink channel. IEEE 9th Workshop on Signal Processing Advances in Wireless Communications (SPAWC'08) (Pernambuco, Brazil, 2008), pp. 196–200

D Aktas, M Bacha, J Evans, S Hanly, Scaling results on the sum capacity of cellular networks with MIMO links. IEEE Trans Inf Theory 52(7), 3264–3274 (2006)

S Chatzinotas, MA Imran, C Tzaras, Uplink capacity of MIMO cellular systems with multicell processing. IEEE International Symposium on Wireless Communication Systems (ISWCS'08) (Reykjavik, Iceland, 2008), pp. 453–457

S Chatzinotas, M Imran, R Hoshyar, On the multicell processing capacity of the cellular MIMO uplink channel in correlated rayleigh fading environment. IEEE Trans Wireless Commun 8(7), 3704–3715 (2009)

O Simeone, O Somekh, H Poor, S Shamai, Distributed MIMO in multicell wireless systems via finitecapacity links. Communications, Control and Signal Processing, 2008. ISCCSP 2008. 3rd International Symposium on, 203–206 (2008)

O Somekh, B Zaidel, S Shamai, Spectral efficiency of joint multiple cellsite processors for Randomly Spread DSCDMA Systems. IEEE Trans Inf Theory 53(7), 2625–2637 (2007)

BL Ng, J Evans, S Hanly, On the Capacity of cellular networks with global LMMSE receiver. IEEE International Conference on Communications, 2007. ICC '07, 870–876 (2007)

S Chatzinotas, MA Imran, R Hoshyar, B Ottersten, multicell LMMSE filtering capacity under correlated Multiple BS Antennas. IEEE 72nd Vehicular Technology Conference (VTC'10Fall), Ottawa, Canada (2010)

S Chatzinotas, B Ottersten, Free probability based capacity calculation of multiantenna Gaussian fading channels with cochannel interference. Physical Communication 4(3), 206–217 (2011). Publisher Full Text

E Katranaras, M Imran, R Hoshyar, Sum rate of linear cellular systems with clustered joint processing. Vehicular Technology Conference, 2009. VTC Spring 2009. IEEE 69th, 1–5 (2009)

W Choi, J Andrews, Downlink performance and capacity of distributed antenna systems in a multicell environment. IEEE Trans Wireless Commun 6, 69–73 (2007)

V Cadambe, S Jafar, Interference alignment and degrees of freedom of the K User interference channel. IEEE Trans Inf Theory 54(8), 3425–3441 (2008)

V Cadambe, S Jafar, Interference alignment and the degrees of freedom of wireless X networks. IEEE Trans Inf Theory 55(9), 3893–3908 (2009)

S Jafar, S Shamai, Degrees of freedom region of the MIMO X channel. IEEE Trans Inf Theory 54, 151–170 (2008)

K Gomadam, V Cadambe, S Jafar, Approaching the capacity of wireless networks through distributed interference alignment. Global Telecommunications Conference, 2008. IEEE GLOBECOM 2008. IEEE, 1–6 (2008)

CM Yetis, T Gou, SA Jafar, AH Kayran, On feasibility of interference alignment in MIMO interference networks. IEEE Trans Signal Process 58(9), 4771–4782 (2010)

R Tresch, M Guillaud, E Riegler, On the achievability of interference alignment in the Kuser constant MIMO interference channel. Statistical Signal Processing, 2009. SSP '09. IEEE/SP 15th Workshop on, 277–280 (2009)

C Suh, D Tse, Interference alignment for cellular networks. Communication, Control, and Computing, 2008 46th Annual Allerton Conference on, 1037–1044 (2008)

R Tresch, M Guillaud, Cellular interference alignment with imperfect channel knowledge. Communications Workshops, 2009. ICC Workshops 2009. IEEE International Conference on, 1–5 (2009)

B Da, R Zhang, Exploiting interference alignment in multicell cooperative OFDMA resource allocation.. Submitted to GlobeComm 2011 (2011) [http://arxiv.org/abs/1103.3093] webcite

N Lee, W Shin, B Clerckx, Interference alignment with limited feedback on twocell interfering TwoUser MIMOMAC.. ICC 2011 (2010) [http://arxiv.org/abs/1010.0933] webcite

O Somekh, O Simeone, Y BarNess, AM Haimovich, U Spagnolini, SS Shitz, An information theoretic view of distributed antenna processing in cellular systems. Distributed Antenna Systems: Open Architecture for Future Wireless Communications, ed. by Zhang Y, Hu H, Luo J (Auerbach Publications, CRC Press, New York, 2007)

J Xu, J Zhang, J Andrews, When does the wyner model accurately describe an uplink cellular network? GLOBECOM 2010, 2010 IEEE Global Telecommunications Conference, 1–5 (2010)

AM Tulino, S Verdú, Random matrix theory and wireless communications. Commun Inf Theory 1, 1–182 (2004)

N Raj Rao, A Edelman, The polynomial method for random matrices. Found Comput Math 8(6), 649–702 (2007)

GJ Foschini, MJ Gans, On limits of wireless communications in a fading environment when Using multiple antennas. Wirel Pers Commun 6(3), 311–335 (1998). Publisher Full Text

R Blum, MIMO capacity with interference. IEEE J Select Areas Commun 21(5), 793–801 (2003). Publisher Full Text

A Lozano, A Tulino, Capacity of multipletransmit multiplereceive antenna architectures. Information Theory, IEEE Transactions on 48(12), 3117–3128 (2002). Publisher Full Text

R Etkin, D Tse, H Wang, Gaussian Interference Channel Capacity to Within One Bit. IEEE Trans Inf Theory 54(12), 5534–5562 (2008)

S Verdú, Multiuser Detection (Cambridge University Press, UK, 1998)

S Chatzinotas, M Imran, C Tzaras, Capacity Limits in Cooperative Cellular Systems. Cooperative Wireless Communications, ed. by Zhang Y, Chen HH, Guizani M (Auerbach Publications, Taylor & Francis Group, Boca Raton, 2009)