SpringerOpen Newsletter

Receive periodic news and updates relating to SpringerOpen.

This article is part of the series Multicell Cooperation for Next Generation Communication Systems.

Open Access Open Badges Research

Interference mitigation techniques for clustered multicell joint decoding systems

Symeon Chatzinotas1* and Björn Ottersten12

Author Affiliations

1 Interdisciplinary Centre for Security, Reliability and Trust, University of Luxembourg, 6, rue Richard Coudenhove-Kalergi, 1359 Luxembourg, Luxembourg

2 Royal Institute of Technology (KTH), Osquldas v. 10, 100 44, Stockholm, Sweden

For all author emails, please log on.

EURASIP Journal on Wireless Communications and Networking 2011, 2011:132  doi:10.1186/1687-1499-2011-132

The electronic version of this article is the complete one and can be found online at: http://jwcn.eurasipjournals.com/content/2011/1/132

Received:25 March 2011
Accepted:14 October 2011
Published:14 October 2011

© 2011 Chatzinotas and Ottersten; licensee Springer.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


Multicell joint processing has originated from information-theoretic 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 per-cell ergodic sum-rate throughput as a figure of merit. In this process, the asymptotic eigenvalue distribution of the channel covariance matrices is analytically derived based on free-probabilistic arguments in order to quantify the sum-rate 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 interference-limited 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 LTE-advanced [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 R-transform 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 closed-form 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, <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M1">View MathML</a> 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 δ(·). In denotes an n × n identity matrix, <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M2">View MathML</a> an n × m matrix of ones, 0 a zero matrix and <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M3">View MathML</a> denotes n × m Gaussian matrix with entries drawn form a <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M4','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M4">View MathML</a> distribution. The figure of merit analyzed and compared throughout this paper is the ergodic per-cell sum-rate 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. "hyper-receiver") which is interconnected to all the BSs through a backhaul of wideband, delayless and error-free 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 path-loss models and user distributions were investigated in [6,7] providing closed-form 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 per-cell sum-rate 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 DS-CDMA [12], where chips act as multiple dimensions. Finally, linear MMSE filtering [13,14] followed by single-user decoding has been considered as an alternative to the optimal multiuser decoder which requires computationally-complex 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 sum-rate capacity performance is degraded by intercluster interference, especially affecting the individual rates of cluster-edge 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 cluster-edge cells [16]. This approach eliminates intercluster interference but at the same time limits the available degrees of freedom. In DS-CDMA 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 [18-20]. Its principle is based on aligning the interference on a signal subspace with respect to the non-intended 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 non-intended 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 non-intended 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 [21-23].

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 cluster-centre UTs. In a similar setting, the authors in [26] propose an IA-based resource allocation scheme which jointly optimizes the frequency-domain 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 limited-size cellular system without BS cooperation, showing that the interference-free dofs can be achieved as the number of UTs grows. Employing IA with unit multiplexing gain towards the non-intended BSs, they study the effect of multi-path channels and single-path 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 non-intended 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 cluster-edge 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 cluster-edge 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 per-cell sum-rate throughput.

thumbnailFigure 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 yi at any random BS can be expressed as follows:

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M5">View MathML</a>


where the n × 1 vector z denotes AWGN with <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M6','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M6">View MathML</a> and <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M7','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M7">View MathML</a>. The K n × 1 vector xi denotes the transmitted symbol vector of the ith UT group with <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M8','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M8">View MathML</a> where γ is the transmit Signal to Noise ratio per UT antenna. The n × Kn channel matrix <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M9','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M9">View MathML</a> 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 αGi, i+1(t)xi+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 systemc. Another intuitive description of the α factor is that it models the power imbalance between intra-cell and inter-cell signals.

Assuming a memoryless channel, the system channel model can be written in a vectorial form as follows:

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M10','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M10">View MathML</a>


where the aggregate channel matrix has dimensions Mn × (M + 1)Kn and can be modelled as:

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M11','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M11">View MathML</a>


with <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M12','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M12">View MathML</a> being a block-Toeplitz matrix and <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M13','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M13">View MathML</a>. In addition, <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M14','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M14">View MathML</a> is a M × M + 1 Toeplitz matrix structured as follows:

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M15','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M15">View MathML</a>


Assuming no CSI at the UTs, the per-cell capacity is given by the MIMO multiple access (MAC) channel capacity:

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M16">View MathML</a>


Theorem 3.1. In the global MJD case, the per-cell capacity for asymptotically large n converges almost surely (a.s.) to the Marcenko-Pastur (MP) law with appropriate scaling [6,10]:

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M17','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M17">View MathML</a>


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 → ∞:

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M18','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M18">View MathML</a>


where λi (X) and <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M19','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M19">View MathML</a> denote the eigenvalues and the asymptotic eigenvalue probability distribution function (a.e.p.d.f.) of matrix X respectively and <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M20','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M20">View MathML</a> denotes the Shannon transform of X with scalar parameter x. It should be noted that <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M21','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M21">View MathML</a> denotes the total UT transmit power normalized by the receiver noise powerd. The last step of the derivation is based on unit rank matrices decomposition and analysis on the R-transform domain, as presented in [6,10]. The scaling factor

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M22','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M22">View MathML</a>


is the Frobenius norm of the Σ matrix <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M23','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M23">View MathML</a> normalized by the matrix dimensions and

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M24','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M24">View MathML</a>


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 cluster-edge UT groups, inspired by [24]. Let us assume a n × 1 unit norm reference vector v with ||v||2 = n and

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M25','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M25">View MathML</a>


<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M26','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M26">View MathML</a>


where y1 and yM represent the received signal vectors at the first and last BS of the cluster, respectively. The first UT group has to align its input x1 towards the non-intended 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 αGM, M+1xM+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 GM, M+1. Let us also assume that <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M27','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M27">View MathML</a> and <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M28','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M28">View MathML</a> represent the transmitter vector and channel matrix of the jth UT in the ith group towards the <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M29','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M29">View MathML</a> BS. In this context, the following precoding scheme is employed to align interference:

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M30','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M30">View MathML</a>


where vj = vvj is a scaled version of v which satisfies the input power constraint <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M31','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M31">View MathML</a>. This precoding results in unit multiplexing gain and is by no means the optimal IA schemee [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:

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M32','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M32">View MathML</a>


It can be easily seen that interference has been aligned across the reference vector and it can be removed using a K× n zero-forcing 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:

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M33','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M33">View MathML</a>


Assuming that the system operates in high-SNR regime and is therefore interference limited, the effect of the AWGN noise colouring <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M34','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M34">View MathML</a> can be ignored, namely <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M35','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M35">View MathML</a>.

Lemma 3.1. The Shannon transform of the covariance matrix of QGM,M is equivalent to that of a K × K Gaussian matrix GK×K.

Proof. Using the property det(I + γAB) = det(I + γBA), it can be written that:

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M36','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M36">View MathML</a>


The K × n truncated unitary matrix Q has K unit singular values and therefore the matrix product QHQ has K unit eigenvalues and a zero eigenvalue. Applying eigenvalue decomposition on QHQ, the left and right eigenvectors can be absorbed by the isotropic Gaussian matrices <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M37','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M37">View MathML</a> and GM,M respectively, while the zero eigenvalue removes one of the n dimensions. Using the definition of Shannon transform [30], Eq. (15) yields

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M38','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M38">View MathML</a>


Based on this lemma and for the purposes of the analysis, QGM,M is replaced by GK × 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)

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M39','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M39">View MathML</a>


where <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M40','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M40">View MathML</a> 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 <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M41','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M41">View MathML</a> 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:

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M42','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M42">View MathML</a>


where <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M43','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M43">View MathML</a>and

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M44','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M44">View MathML</a>


with <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M45','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M45">View MathML</a> being a n × (M - 1)Kn + K matrixf, <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M46','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M46">View MathML</a> being a (M - 2)n × (M - 1)Kn + K matrix and <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M47','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M47">View MathML</a> being a n - 1 × (M - 1)Kn + K matrixg.

Since all intercluster interference has been filtered out and the effect of filter Q has been already incorporated in the structure of HIA, the per-cell throughput in the IA case is still given by the MIMO MAC expression:

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M48','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M48">View MathML</a>


Theorem 3.2. In the IA case, the per-cell throughput can be derived from the R-transform of the a.e.p.d.f. of matrix <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M49','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M49">View MathML</a>.

Proof. Following an asymptotic analysis where n→ ∞:

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M50','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M50">View MathML</a>


The a.e.p.d.f. of <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M49','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M49">View MathML</a>is obtained by determining the imaginary part of the Stieltjes transform <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M51','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M51">View MathML</a> for real arguments

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M52','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M52">View MathML</a>


considering that the Stieltjes transform is derived from the R-transform [31] as follows

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M53','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M53">View MathML</a>


Theorem 3.3. The R-transform of the a.e.p.d.f. of matrix <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M49','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M49">View MathML</a>is given by:

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M54','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M54">View MathML</a>


with k, β, q parameters given by:

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M55','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M55">View MathML</a>

and <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M56','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M56">View MathML</a>given by theorem A.1.

Proof. Based on Eq.(19), the matrix <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M57','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M57">View MathML</a>can be decomposed as the following sum:

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M58','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M58">View MathML</a>


where H1 = Σ1 Gn×(M-1)Kn+K, H2 = Σ2 G(M-2)n×(M-1)Kn+K and H3 = Σ3 Gn-1×(M-1)Kn+K. Using the property of free additive convolution [30] and Theorem A.1 in Appendix A, Eq. (24) holds in the R-transform domain.   □

3.3 Resource division multiple access

RDMA entails that the available time or frequency resources are divided into two orthogonal parts assigned to cluster-edge cells in order to eliminate intercluster interference [[16], Efficient isolation scheme]. More specifically, for the first part cluster-edge UTs are inactive and the far-right cluster-edge BS is active. For the second part, cluster-edge UTs are active and the far-right cluster-edge BS is inactive. While the available channel resources are cut in half for cluster-edge 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 per-cell throughput in the RDMA case is given by:

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M59','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M59">View MathML</a>


where <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M60','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M60">View MathML</a> and <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M61','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M61">View MathML</a> 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:

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M62','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M62">View MathML</a>


Theorem 3.4. For the first part of the RDMA case, the per-cell throughput CRD1 can be derived from the R-transform of the a.e.p.d.f of matrix <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M63','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M63">View MathML</a>, where:

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M64','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M64">View MathML</a>


Proof. Following an asymptotic analysis where n→ ∞:

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M65','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M65">View MathML</a>


Using the matrix decomposition of Eq. (27) and free additive convolution [30]:

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M66','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M66">View MathML</a>


Eq. (28) follows from Eq. (42) with <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M67','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M67">View MathML</a> 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:

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M68','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M68">View MathML</a>


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 per-cell throughput <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M61','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M61">View MathML</a>can be derived from the R-transform of the a.e.p.d.f. of matrix <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M69','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M69">View MathML</a>, where:

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M70','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M70">View MathML</a>


Proof. Following an asymptotic analysis where n → ∞:

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M71','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M71">View MathML</a>


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 cluster-edge UT group contribution αGM, M+1(t)xM+1(t) is considered as interference. As a result, the interference channel matrix can be expressed as:

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M72','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M72">View MathML</a>


Assuming no CSI at the UTs, the per-cell throughput in the CI case is given by [15,32-34]:

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M73','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M73">View MathML</a>


where CI denotes the throughput of the interfering UT group normalized by the cluster size:

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M74','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M74">View MathML</a>


Theorem 3.6. In the CI case, the per-cell throughput converges almost surely (a.s.) to a difference of two scaled versions of the the MP law:

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M75','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M75">View MathML</a>


Proof. Following an asymptotic analysis in the number of antennas n n →∞:

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M76','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M76">View MathML</a>


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:

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M77','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M77">View MathML</a>


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):

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M78','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M78">View MathML</a>


Remark 3.1. It can be observed that dIA = dRD only for single UT per cell equipped with two antennas (K = 1, n = 2). For all other cases, dIA > dRD. Furthermore, it is worth noting that when the number of UTs K and antennas n grows to infinity, limK,n → ∞ dIA = dMJD 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 UT-side 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 jointly-processed 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 inter-cluster synchronization, while frequency division could be even static. Finally, IA is the most complex since CSI towards the non-intended 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 103 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 per-cell 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 per-cell 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 per-cell 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 per-cell 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.

thumbnailFigure 2. Per-cell 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.

thumbnailFigure 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.

thumbnailFigure 4. Per-cell 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.

thumbnailFigure 5. Per-cell 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 free-probabilistic arguments. In addition, the RDMA scheme was studied as a low complexity method for mitigating intercluster interference. Finally, the CI was considered as a worst-case 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 variance-profiled Gaussian matrix B = C G and a number of zero columns. Let also k be the ratio of non-zero 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 R-transform of AH A is given by:

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M79','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M79">View MathML</a>


Proof. Let B = C Gn × m be a variance-profiled Gaussian matrix with β = m/n and q = ||C||2/nm. According to [10], the R-transform of <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M80','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M80">View MathML</a> is given by:

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M81','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M81">View MathML</a>


Using eq. (23), the Stieltjes transform of <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M80','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M80">View MathML</a> can be expressed as:

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M82','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M82">View MathML</a>


Matrix <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M83','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M83">View MathML</a> has identical eigenvalues to <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M80','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M80">View MathML</a> plus a number of zero eigevalues with 0 < k < 1 defined as the ratio of non-zero eigenvalues over the total number of eigenvalues. As a result, the a.e.p.d.f. of <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M83','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M83">View MathML</a>can be written as:

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M84','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M84">View MathML</a>


Using the definition of the Stieltjes transform [30]:

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M85','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M85">View MathML</a>


and employing eq. (23), the proof is complete.   □

5 Competing interests

The authors declare that they have no competing interests.


aThe term throughput is used instead of capacity since the described techniques are suboptimal in the information-theoretic sense and lead to achievable sum-rates except for MJD which leads to MIMO MAC capacity.

bThe multiple antennas are assumed to be uncorrelated although the analytical results can be extended in the correlated case based on the principles described in

cFor more details on the modelling of the α parameter, the reader is referred to [37].

dFor the purposes of the analysis the variable <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M86','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M86">View MathML</a> is kept finite as the number of antennas Mn grows large, so that the system power does not grow to infinity.

eDepending 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

fThe structure of the first block of Σ1 originates in the Gaussian approximation of <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M87','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2011/1/132/mathml/M87">View MathML</a>.

gThe structure of the last block of Σ3 is based on Lemma 3.1.


This work was partially supported by the National Research Fund, Luxembourg under the CORE project "CO2SAT: Cooperative and Cognitive Architectures for Satellite Networks".


  1. P Marsch, G Fettweis, Coordinated Multi-Point in Mobile Communications From Theory to Practice (Cambridge University Press, 2011)

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

  3. SV Hanly, PA Whiting, Information-theoretic capacity of multi-receiver networks. Telecommun Syst 1, 1–42 (1993). Publisher Full Text OpenURL

  4. A Wyner, Shannon-theoretic approach to a Gaussian cellular multiple-access channel. IEEE Trans Inf Theory 40(6), 1713–1727 (1994). Publisher Full Text OpenURL

  5. O Somekh, S Shamai, Shannon-theoretic approach to a Gaussian cellular multiple-access channel with fading. IEEE Trans Inf Theory 46(4), 1401–1425 (2000). Publisher Full Text OpenURL

  6. 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)

  7. 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

  8. 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)

  9. 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

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

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

  12. O Somekh, B Zaidel, S Shamai, Spectral efficiency of joint multiple cell-site processors for Randomly Spread DS-CDMA Systems. IEEE Trans Inf Theory 53(7), 2625–2637 (2007)

  13. 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)

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

  15. 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 OpenURL

  16. 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)

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

  18. 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)

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

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

  21. 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)

  22. 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)

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

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

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

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

  27. N Lee, W Shin, B Clerckx, Interference alignment with limited feedback on two-cell interfering Two-User MIMO-MAC.. ICC 2011 (2010) [http://arxiv.org/abs/1010.0933] webcite OpenURL

  28. O Somekh, O Simeone, Y Bar-Ness, 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)

  29. 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)

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

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

  32. 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 OpenURL

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

  34. A Lozano, A Tulino, Capacity of multiple-transmit multiple-receive antenna architectures. Information Theory, IEEE Transactions on 48(12), 3117–3128 (2002). Publisher Full Text OpenURL

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

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

  37. 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)