## Quantized CSI-Based Tomlinson-Harashima Precoding in Multiuser MIMO Systems

Liang Sun, Member , IEEE and Ming Lei

## Abstract

This paper considers the implementation of Tomlinson-Harashima (TH) precoding for multiuser MIMO systems based on quantized channel state information (CSI) at the transmitter side. Compared with the results in [1], our scheme applies to more general system setting where the number of users in the system can be less than or equal to the number of transmit antennas. We also study the achievable average sum rate of the proposed quantized CSI-based TH precoding scheme. The expressions of the upper bounds on both the average sum rate of the systems with quantized CSI and the mean loss in average sum rate due to CSI quantization are derived. We also present some numerical results. The results show that the nonlinear TH precoding can achieve much better performance than that of linear zero-forcing precoding for both perfect CSI and quantized CSI cases. In addition, our derived upper bound on the mean rate loss for TH precoding converges to the true rate loss faster than that of zeroforcing precoding obtained in [2] as the number of feedback bits becomes large. Both the analytical and numerical results show that nonlinear precoding suffers from imperfect CSI more than linear precoding does.

## Index Terms

Tomlinson-Harashima precoding, QR decomposition, random vector quantization, zero-forcing, Givens transformation.

## I. INTRODUCTION

Since the pioneering work [3] and [4], multiple-input multiple-output (MIMO) communication systems have been extensively studied in both academic and industry communities and becomes the key technology of most emerging wireless standards. It is shown that significantly enhanced spectral efficiency and link reliability can be achieved compared with conventional single antenna

systems [3, 5]. In the downlink multiuser MIMO systems, multiple users can be simultaneously served by exploiting the spatial multiplexing capability of multiple transmit antennas, rather than trying to maximize the capacity of a single-user link.

The performance of a MIMO system with spatial multiplexing is severely impaired by the multi-stream interference due to the simultaneous transmission of parallel data streams. To reduce the interference between the parallel data streams, both the processing of the data streams at the transmitter (precoding) and the processing of the received signals (equalization) can be used. Precoding matches the transmission to the channel. Accordingly, linear precoding schemes with low complexity are based on zero-forcing (ZF) [6] or minimum mean-squareerror (MMSE) criteria [7] and their improved version of channel regularization [8]. In spite of very low complexity, the linear schemes suffer from capacity loss. Nonlinear processing at either the transmitter or the receiver provides an alternative approach that offers the potential for performance improvements over the linear approaches. This kind of approaches includes schemes employing linear precoding combined with decision feedback equalization (DFE) [5, 9], vector perturbation [10], Tomlinson-Harashima (TH) precoding [1, 11], and ideal dirty paper coding [12, 13] which is too complex to be implemented in practice. Vector perturbation has been proposed for multiuser MIMO channel model and can achieve rate near capacity [10]. It has superior performance to linear precoding techniques, such as zero-forcing beamforming and channel inversion, as well as TH precoding [10]. However, this method requires the joint selection of a vector perturbation of the signal to be transmitted to all the receivers, which is a multi-dimensional integer-lattice least-squares problem. The optimal solution with an exhaustive search over all possible integers in the lattice is complexity prohibited. Although some suboptimal solutions, such as sphere encoder [14], exist, the complexity is still much higher than TH precoding.

TH precoding can be viewed as a simplified version of vector perturbation by sequential generation of the integer offset vector instead of joint selection. This technique employs modulo arithmetic and has a complexity comparable to that of linear precoders. It was originally proposed to combat inter-symbol interference in highly dispersive channels [15] and can readily be extended to MIMO channels [1, 16]. Although it was shown in [10] that TH precoding does not perform nearly as well as vector perturbation for general SNR regime, it can achieve significantly

better performance than the linear pre-processing algorithm, since it limits the transmitted power increase while pre-eliminating the inter-stream interference [11]. Thus, it provides a good choice of tradeoff between performance and complexity and has recently received much attention [1, 11]. Note that TH precoding is strongly related to dirty paper coding. In fact, it is a suboptimal implementation of dirty paper coding proposed in [17].

As many precoding schemes, the major problem for systems with TH precoding is the availability of the channel state information (CSI) at the transmitter. In time division duplex systems, since the channel can be assumed to be reciprocal, the CSI can be easily obtained from the channel estimation during reception. In frequency division duplex (FDD) systems, the transmitter cannot estimate this information and the CSI has to be communicated from the receivers to the transmitter via a feedback channel. In this paper, we will focus on the implementation of TH precoding in FDD systems. In this context, for linear precoding, there have been extensive research results for MIMO systems with quantized CSI at the transmitter [2, 18, 19]. However, as far as we know, there has been very few works directed at the design of TH precoding based on the quantized CSI at the transmitter side. In this respect, the previous design in [20] is based on MMSE criteria. Since the MSE is a function of both statistics (moments) of the channels and the statistics of the channel quantization error, the computation of the MSE requires the exact distribution of the channels which can be very difficult to obtain in practical systems. In addition, even if the exact distribution function of channels could be obtained, the statistics of quantization error can be very difficult to obtain for more general channel fading other than uncorrelated Rayleigh fading even with simple random vector quantization (RVQ) codebook. Instead, we aim to design low complexity method which can be easily implemented in practical systems with arbitrary channel fading. Our scheme employs a more direct method which only depends on the quantized CDI of user channels.

In this paper, we design a multiuser spatial TH precoding based on quantized CSI and ZF criteria. As in [1], we focus on high spectral efficiency, in particular non-binary modulation alphabets and correspondingly we assume high signal-to-noise ratios (SNRs). In contrast to [1] where perfect CSI is at the transmitter side, we assume only quantized CDI is available at the transmitter. The feedforward filter as well as the feedback filter are computed at the transmitter only based on the available quantized CDI at the transmitter side. In addition, our scheme also

Fig. 1 . TH precoding for multiuser MIMO downlink with quantized CSI feedback.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0455v1-figure-930e6095-57c3-4bd9-bf40-cf2c6f0fe36d.jpg)

generalizes the results in [1] to more general system setting where the number of users K in system can be less than or equal to the number of transmit antennas n T . We also study the achievable average sum rate of the proposed quantized CSI-based TH precoding scheme by analytically characterizing the average sum rate and the rate loss due to quantized CSI as functions of the number of feedback bits per user. Our derived upper bound for TH precoding tracks the true rate loss quite closely and appears to converge faster than the upper bound for ZF precoding obtained in [2] as the number of feedback bits becomes large.

## II. SYSTEM MODEL

As shown in Fig. 1, we consider the multi-user downlink systems where TH precoding [15] is used at the transmitter for multi-user interference pre-subtraction. The transmitter is equipped with n T transmit antennas and K decentralized users each has a single antenna such that K ≤ n T . Let the vector s = [ s , 1 · · · , s K ] T ∈ C K × 1 represent the modulated signal vector for all users, where s k is the k -th modulated symbol stream for user k . Here we assume that an M -ary square constellation ( M is a square number) is employed in each of the parallel data streams and the constellation set is A = { s I + js Q | s , s I Q ∈ ± 1 √ 3 2( M -1) , ± 3 √ 3 2( M -1) , · · · , ± √ M √ 3 2( M -1) } .

In general, the average transmit symbol energy is normalized, i.e. E {| s k | 2 } = 1 . s is fed to the precoding unit, which consists of a backward square matrix B and a nonlinear operator MOD τ ( ) · which acts independently over the real and imaginary parts of its input as follows

$$\text{MOD} _ { \tau } ( x ) = x - \tau \left \lfloor \frac { x + \tau } { 2 \tau } \right \rfloor,$$

where τ = √ M √ 3 2( M -1) , /floorleft z /floorright is the largest integer not exceeding z . B must be strictly lower triangular to allow data precoding in a recursive fashion [1]. The construction of B will depend on the level of CSI of the supported users available at the transmitter side. If we temporarily neglect the nonlinear operator MOD τ ( ) · in Fig. 1, the channel signal vector x = [ x , 1 · · · , x K ] T can be generated as

$$x _ { 1 } = s _ { 1 },$$

$$x _ { k } = s _ { k } - \sum _ { l = 1 } ^ { k - 1 } [ \mathbf B ] _ { k, l } x _ { l }, \ k = 2, \cdots, K. \quad \quad \quad$$

In this way, if b = ( I + B ) -1 s become large in the presence of deep fading the transmit power can be greatly increased. TH precoding modulo (1) reduces the transmit symbols into the boundary square region of R = { x + jy x, y | ∈ ( -τ, τ ) } . With (1) and (2) the channel signals are equivalently given as

$$x _ { k } = s _ { k } + d _ { k } - \sum _ { l = 1 } ^ { k - 1 } [ \mathbf B ] _ { k, l } x _ { l }, \ k = 2, \cdots, K, \\ \$$

where d k ∈ { 2 τ ( p I + jp Q ) | p , p I Q ∈ Z } is properly selected to ensure the real and imaginary parts of x k are constrained into R [1]. The constellation of the modified data symbols v k = s k + d k is simply the periodic extension of the original constellation along the real and imaginary axes. Equivalently, the effective data symbols v k ( k = 2 , · · · , K ) are passed into B , which is implemented by the feedback structure. Thus, we have

$$\mathbf v = \mathbf C \mathbf x,$$

where v = [ v , v 1 2 , · · · , v K ] T and

$$C = B + I.$$

We will make the standard observation that the elements of x are almost uncorrelated and uniformly distributed over the Voronoi region of the constellation R , and that such a model becomes more precise as n T increases [11, Theorem 3.1]. With E { ss H } = I , the covariance of x can be accurately approximated as R x = M M -1 I [11]. Moreover, the induced shaping loss by the non-Gaussian signaling leads to the fact that the achievable rate can be up to 1 53 . dB from the channel capacity [21]. However, as indicated in [1], the so-called shaping loss can be bridged by higher-dimensional precoding lattices. A scheme named 'inflated lattice' precoding has been proved to be capacity-achieving in [22]. Thus, following [1], we will ignore the shaping gap in this work.

A spatial channel pre-equalization is performed at the transmitter side using a feedforward precoding matrix F ∈ C n T × K . Throughout this work, we assume equal power allocation to all supported users. Then the received signal can be written as

$$r = \sqrt { \frac { P } { \kappa } } \text{HF} x + \text{n},$$

where H = [ h T 1 , · · · , h T K ] T is the compact flat fading channel matrix consisting of all users's channel vectors and h k ∈ C n T × 1 is the channel from the transmitter to user k 1 . κ is used for transmit power normalization. F ∈ C n T × K satisfies transmit power constraint P κ Tr { FR F x H } = P κ M M -1 Tr { FF H } = P . As for B F , is also designed based on the level of CSI available to the transmitter. We assume n is the white additive noise at all the receivers with the covariance R n = I without loss of generality. Each receiver compensates for the channel gain by dividing by a factor g k prior to the modulo operation as follows:

$$y = G \left ( \sqrt { \frac { P } { \kappa } } \text{HF} x + \text{n} \right ),$$

where G = diag ( g 1 1 , , · · · , g K,K ) .

Throughout this work, we assume each receiver can obtain perfect CSI of his own through channel estimation and feeds back this information to the transmitter via a zero-delay feedback link with possible rate-constraints. In this part, we assume this feedback information is perfect

1 The ordering of the users' channel vectors in H will affect the precoding order of the users' information signals and further affect the performance of each user. However, at this stage we assume the user channel vectors are randomly ordered. Thus the TH precoding order of the users is 1 2 , , · · · , K .

at the transmitter. The design of TH precoding for MU MIMO systems with perfect CSI has been studied in [1] where, for simplicity, the authors restrict the work to the systems with equal number of transmit antennas and users. In this part, we review and extend that construction for systems with arbitrary number of users no more that of the transmit antennas.

With perfect CSI at the transmitter, let the QR decomposition of the compact channels be H = RQ , where R = [ r i,j ] ∈ C K × K is a lower left triangular matrix and Q ∈ C K × n T is a semiunitary matrix with orthonormal rows which satisfies QQ H = I . Then the precoding matrix F is given as F = Q H , the scaling matrix G is given as G = √ κ P ∆ with ∆ = diag ( r -1 1 1 , , · · · , r -1 K,K ) and the feedback matrix reads B = ∆HF -I = ∆R -I . According to the transmit power constraint P κ M M -1 K = P , we have κ = M M -1 K . With the processing, the effective received data symbols y corrupted by additive noise can be written as [1]

$$y = v + G n.$$

At the receivers, each symbol in y is firstly modulo reduced into the boundary region of the signal constellation A . A quantizer of the original constellation will follow the modulo operation to detect the received signals. The SNR ξ k for receiver k can be written as

$$\xi _ { k } = \frac { P } { \kappa } | r _ { k, k } | ^ { 2 }.$$

In the following part, we will describe how to implement the precoding with quantized CSI obtained at the transmitter.

## III. SYSTEM WITH QUANTIZED TRANSMIT CSI

In practical systems, perfect CSI is never available at the transmitter. For example, in a FDD system, the transmitter obtains CSI for the downlink through the limited feedback of B bits by each receiver. Following the studies of quantized CSI feedback in [2, 19], channel direction vector is quantized at each receiver, and the corresponding index is fed back to the transmitter via an error and delay-free feedback channel. Given the quantization codebook W = { w 1 , · · · , w n } ( w i ∈ C 1 × n T ), which is known to both the transmitter and all the receivers, the k -th receiver

selects the quantized channel direction vector of its own channel as follows:

$$\hat { h } _ { k } = \arg \max _ { w _ { i } \in \mathbb { W } } \{ | \bar { h } _ { k } w _ { i } | ^ { 2 } \},$$

where ¯ h k = h k ‖ h k ‖ is the channel direction vector of user k .

In this work, we use RVQ codebook, in which the n quantization vectors are independently and isotropically distributed on the n T -dimensional complex unit sphere. Although RVQ is suboptimal for a finite-size system, it is very amenable to analysis and also its performance is close to the optimal quantization [2]. Using the result in [2], for user k we have

$$\bar { h } _ { k } = \hat { h } _ { k } \cos \theta _ { k } + \tilde { h } _ { k } \sin \theta _ { k },$$

where cos 2 θ k = | ¯ h h k ˆ H k | 2 , ˜ h k ∈ C 1 × n T is a unit norm vector isotropically distributed in the orthogonal complement subspace of ˆ h k and independent of sin θ k . Then H can be written as

$$H = \Gamma \left ( \Phi \hat { H } + \Omega \tilde { H } \right ), \quad \quad \quad ( 1 2 ) \\. \quad. \quad.. \. \. \. \. \. \. \. \. \. \.$$

/negationslash

With the quantized CDI at the transmitter side, the transmitter obtains the feedforward precoding matrix F and feedback matrix B through the QR decomposition of compact channel matrix ˆ H in the same way as the QR decomposition of matrix H , i.e. ˆ H = ˆ ˆ RQ , where the matrices ˆ R and ˆ Q have the same structure as the matrices R and Q respectively. Then we have F = ˆ Q H and B = ( diag { ˆ R }) -1 ˆ R -I . In addition, the scaling matrix at the receivers now becomes where Γ = diag ( ρ , 1 · · · , ρ K ) with ρ k = ‖ h k ‖ , Φ = diag (cos θ , 1 · · · , cos θ K ) and Ω = diag ( sin θ , 1 · · · , sin θ K ) , ˆ H = [ ˆ h T 1 , · · · , ˆ h T K ] T and ˜ H = [ ˜ h T 1 , · · · , ˜ h T K ] T . For simplicity of analysis, in this work we consider the quantization cell approximation used in [19, 23], where each quantization cell is assumed to be a Voronoi region of a spherical cap with surface area approximately equal to 1 n of the total surface area of the n T -dimensional unit sphere. For a given codebook W , the actual quantization cell for vector w i , R i = { ¯ : h | ¯ hw i | 2 ≥ | ¯ hw j | 2 , ∀ i = j } , is approximated as ˜ R ≈ i { ¯ : h | ¯ hw i | ≥ 1 -δ } , where δ = 2 -B n T -1 .

$$G = \sqrt { \frac { \kappa } { P } } \left ( \Gamma \Phi \text{ diag} \left \{ \hat { R } \right \} \right ) ^ { - 1 }. \\.. \. \. \. \. \. \. \. \ \. \ \.$$

Using the same operation at the receiver side as that in perfect CSI case to detect the received

signals, the detected signal vector ˆ y can be further written as

$$\text{signs, use usereceipts} \, \text{$y$ can we runen when as} \\ \hat { y } & = G \left ( \sqrt { \frac { P } { \kappa } } \text{H} F x + n \right ) \\ & = G \sqrt { \frac { P } { \kappa } } \Gamma \left ( \hat { \Phi } \hat { H } + \Omega \tilde { H } \right ) F x + G n \\ & = v + \left ( \Phi \, \text{diag} \left \{ \hat { R } \right \} \right ) ^ { - 1 } \Omega \tilde { H } \hat { Q } ^ { H } x + \sqrt { \frac { \kappa } { P } } \left ( \Gamma \Phi \, \text{diag} \left \{ \hat { R } \right \} \right ) ^ { - 1 } n, \\ \text{where we have used the relationship} \, v = \left ( \text{diag} \left \{ \hat { R } \right \} \right ) ^ { - 1 } \hat { R } x. \, \text{In} \, ( 14 ), \, \text{the first term is the useful}$$

where we have used the relationship v = ( diag { ˆ R }) -1 ˆ Rx . In (14), the first term is the useful signal vector for all the users and the second term is interference signal caused by the quantized CSI.

According to (14), the output signal-to-interference-plus-noise ratio (SINR) γ k for receiver k can be written as

$$\gamma _ { k } & = \frac { 1 } { \frac { \sin ^ { 2 } \theta _ { k } } { | \hat { r } _ { k, k } | ^ { 2 } \cos ^ { 2 } \theta _ { k } } } \| \tilde { \mathbf h } _ { k } \hat { Q } ^ { H } \| ^ { 2 } + \frac { \kappa } { P } \frac { 1 } { \rho _ { k } ^ { 2 } | \hat { r } _ { k, k } | ^ { 2 } \cos ^ { 2 } \theta _ { k } } } \\ & = \frac { \frac { P } { \kappa } \rho _ { k } ^ { 2 } | \hat { r } _ { k, k } | ^ { 2 } \cos ^ { 2 } \theta _ { k } } { \frac { P } { \kappa } \rho _ { k } ^ { 2 } \| \tilde { \mathbf h } _ { k } \hat { Q } ^ { H } \| ^ { 2 } \sin ^ { 2 } \theta _ { k } + 1 }. \\ \text{AVERAGE SUM RATE ANALYSIS UNDER QUANTIZED CSI FEEDBACK}$$

## IV. AVERAGE SUM RATE ANALYSIS UNDER QUANTIZED CSI FEEDBACK

In this section we will study the achievable average sum rate of the proposed quantized CSI feedback TH precoding scheme. Although the exact distribution of each term in the expression of the output SINR γ k in (15) can be obtained (see for the detailed information), these terms are located at both the numerator and the denominator in (15) . Thus, to obtain the exact closed-form expression of the distribution of output SINR γ k can be very difficult if not impossible, not to mention the exact closed-form expression of the average sum rate. Thus, to simplify analysis, we have appealed to studying some bounds of the average sum rate and the average sum rate loss instead of exact results. For tractability, throughout this section we assume each user's channel is Rayleigh-faded. In the following subsection, we will first study the statistical distribution of the power of interference signal at each user caused by quantized CSI.

## A. Interference Part

In this subsection, assuming Rayleigh fading channel and RVQ for quantized CSI feedback, we will derive the statistical distribution of interference part P κ ρ 2 k ‖ ˜ h Q k ˆ H ‖ 2 sin 2 θ k in (15). It is

well known that ρ 2 k has a χ 2 2 n T distribution and the distribution of sin 2 θ k is given in [2, 18]. However, since ˜ h k ⊥ ˆ h k ( k = 1 , · · · , K ) and ˆ Q is determined by ˆ h k ( k = 1 , · · · , K ), ˜ h k for k = 1 , · · · , K are not independent of ˆ Q . The distribution of the term ‖ ˜ h Q k ˆ H ‖ 2 is still unknown and to obtain the exact result is not trivial. The following lemma presents the exact distribution of this interference term. It is one of the key contributions of this paper.

Lemma 1: For 1 &lt; K &lt; n T , the random variables ε k = ‖ ˜ h Q k ˆ H ‖ 2 for k = 1 , · · · , K follow the same beta distribution with shape ( K -1) and ( n T -K ) which is denoted as ε k ∼ Beta( K -1 , n T -K ) . In addition, the probability density function (p.d.f.) of ε k is given as

$$f _ { \varepsilon _ { k } } ( x ) & = \frac { 1 } { \beta ( K - 1, n _ { T } - K ) } x ^ { K - 2 } ( 1 - x ) ^ { n _ { T } - K - 1 }. \\.$$

where β a, b ( ) = ∫ 1 0 t a -1 t b -1 d t is beta function [24]. Specially, when K = 1 there is no interference term. When K = n T , ε k = ‖ ˜ h Q k ˆ H ‖ 2 is equal to 1 which is a constant.

## Proof: See Appendix A.

Lemma 1 implies a very interesting result that, with randomly ordered user channel vectors, the signal of the user which is precoded ahead suffers from the same interference signal power as the signals of the users which are precoded afterwards. In the following we will only focus on the general situation that 1 &lt; K &lt; n T . However, it is easy to check that all the obtained results also apply to the special cases of K = 1 and K = n T .

The expectation of the logarithm of the interference term ε k , which is shown to be useful in the following theorems, is obtained in the following lemma.

Lemma 2: The expectation of the logarithm of the interference term ε k is given by

$$\mathbb { E } _ { \mathbf h, \mathbf w } & \left [ - \log _ { 2 } \left ( \| \tilde { \mathbf h } _ { k } \hat { Q } ^ { \mathbf \mu } \| ^ { 2 } \right ) \right ] = \log _ { 2 } e \sum _ { m = K - 1 } ^ { n _ { T } - 2 } \sum _ { l = 0 } ^ { n _ { T } - m - 2 } \frac { ( n _ { T } - 2 )! } { m! \ l! \ ( n _ { T } - m - 2 - l )! } ( - 1 ) ^ { l } \frac { 1 } { m + l } \quad ( 1 7 ) \\ \mathfrak { s } \quad \ e \ \circ \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ a c e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad 0 \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ u n e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad \ e \quad$$

Proof:

$$\gamma \colon \text{See Appendix B}.$$

## B. Upper Bounds on the Average Sum Rate Loss and Sum Rate

The instantaneous achievable rates for user k with perfect CSI and quantized CSI feedback are given as

$$R _ { P, k } = \log _ { 2 } \left ( 1 + \xi _ { k } \right )$$

and

$$R _ { Q, k } = \log _ { 2 } ( 1 + \gamma _ { k } ),$$

respectively. The following theorem quantifies the average sum rate performance degradation as a function of the feedback rate.

Theorem 1: With B feedback bits per user, the average sum rate loss of user k due to quantized CSI feedback can be upper bounded by 2

$$\Delta R _ { k } & = \mathbb { E } _ { \mathbf H, W } \{ R _ { P, k } - R _ { Q, k } \} \\ & \leq \Delta R = \log _ { 2 } \left ( 1 + c P \ 2 ^ { - \frac { B } { n _ { T } - 1 } } \right ) + \frac { \log _ { 2 } ( e ) } { n _ { T } - 1 } \sum _ { \colon, \, \cdot } ^ { n _ { T } - 1 } \beta \left ( n, \frac { i } { n _ { T } - 1 } \right ),$$

$$\leq \Delta R = \log _ { 2 } \left ( 1 + c P \ 2 ^ { - \frac { B } { n _ { T } - 1 } } \right ) + \frac { \log _ { 2 } ( e ) } { n _ { T } - 1 } \sum _ { i = 1 } ^ { n _ { T } - 1 } \beta \left ( n, \frac { i } { n _ { T } - 1 } \right ), \\ \left ( \tau _ { - 1 } \aleph _ { - } \right ) \,. \quad. \quad. \quad. \quad. \quad. \quad.$$

where c = ( K -1) n T κ n ( T -1) and n = 2 B is the size of codebook.

Proof:

$$\gamma \colon \text{See Appendix C.}$$

According to the results in [2, Theorem 1 ], the average sum rate loss due to quantized feeback for ZF precoding is upper bounded by ∆ R zf &lt; log 2 ( 1 + P 2 -B n T -1 ) . We find the first term at the right hand side (RHS) of (20) can be approximated as log 2 ( 1 + P 2 -B n T -1 ) with high order constellation, large number of transmit antennas and large number of supported users. Thus, the second term at the RHS of (20) can be seen as the sum rate degradation of nonlinear precoding compared with that of linear precoding when only quantized CSI is available at the transmitter side. In addition, similar to the results for the linear ZF beamforming in [2], the rate loss for nonlinear precoding is also an increasing function of the system SNR ( P ). Thus, the system with fixed feedback rate is interference-limited at high SNR regime, which is shown in the following theorem.

Theorem 2: The average sum rate of user k achieved by quantized CSI-based TH precoding with B feedback bits per user is bounded as

$$R _ { Q, k } & \leq \log _ { 2 } e \left ( \sum _ { m = K - 1 } ^ { n _ { T } - 2 } \sum _ { l = 0 } ^ { n _ { T } - m - 2 } \frac { ( n _ { T } - 2 )! } { m! \ l! \ ( n _ { T } - m - 2 - l )! } ( - 1 ) ^ { l } \frac { 1 } { m + l } + \frac { 1 } { n _ { T } - 1 } \sum _ { l = 1 } ^ { n } \frac { 1 } { l } \right ), \quad ( 2 1 ) \\. & \quad \. \. \. \. \. \. \.$$

where n = 2 B is the size of codebook.

2 Note that, in contrast to ZF precoding, for TH precoding different users have different average sum rate loss. Interestingly, simulation results show that, for finite SNR, the users precoded earlier will suffer from greater sum rate loss. However, in this work will adopt the average sum rate loss over all supported users.

## Proof: See Appendix D.

We can see from this theorem that, with fixed feedback bits per user, as the interference and signal power both increase linearly with P , the system becomes interference-limited and the average sum rate converges to an upper bound. These can also be observed from the simulation results in Fig. 3.

In the context of linear ZF precoding in [2], the author showed the interference-limited scenario can be avoided by scaling the feedback rate linearly with the SNR P dB (in decibels). Particularly, it is showed in [2, Theorem 3] that in order to maintain a constant average sum rate loss no greater than log 2 b bits per user between the system with perfect CSI and the system with finite-rate feedback, it is sufficient to scale the number of feedback bits per user according to

$$B = ( n _ { T } - 1 ) \frac { \log _ { 2 } 1 0 } { 1 0 } P _ { d \mathrm B } - ( n _ { T } - 1 ) \log _ { 2 } ( b - 1 ).$$

However, for nonlinear TH precoding, the explicit relationship between the feedback rate and the SNR to maintain a constant average sum rate loss cannot be easily obtained. This is mainly due to the fact that the expression of average sum rate loss in (20) is a much more complex function of n ( B ) than the corresponding expression for linear ZF precoding given in [2, Theorem 1]. In the following we will derive a corresponding relationship for the system employing TH precoding. First, in Appendix E we show that the second term at the RHS of (20) can be bounded by a decreasing function of n for a fixed n T . In addition, this upper bound approaches zero as n → ∞ . Thus, as n scales linearly with SNR (in decibels), for an arbitrary given constant 0 &lt; ε &lt; 1 , we can always find a positive integer N ε ( ) such that, whenever n ≥ N ε ( ) ,

$$\Delta R \leq \log _ { 2 } \left ( 1 + c P \ 2 ^ { - \frac { B } { n _ { T } - 1 } } \right ) + \varepsilon. \\ \dots \dots \dots \dots \dots \dots$$

To characterize a sufficient condition of the scaling of feedback rate, we set the RHS of (23) to be the maximum allowable gap of log 2 b . After some simple manipulations, we get

$$B & = ( n _ { T } - 1 ) \log _ { 2 } P - \log _ { 2 } ( b - 2 ^ { \varepsilon } - 1 ) + \log _ { 2 } c \\ & = ( n _ { T } - 1 ) \frac { \log _ { 2 } 1 0 } { 1 0 } P _ { \text{dB} } - \log _ { 2 } ( b - 2 ^ { \varepsilon } - 1 ) + \log _ { 2 } c.$$

Fig. 2 . 4 × 4 system with increasing number of feedback bits.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0455v1-figure-a4be74ca-60d8-4ae4-83ab-61e66247f891.jpg)

In Fig. 2, the average sum rate curves are shown for a system with n T = 4 and K = 4 . The feedback rate is assumed to scale according to the relationship given in (24). Notice that, since ε can be set to be a small number when B is large enough, in the simulation we set ε = 0 to get a stronger condition than (24). Quantized CSI-based TH precoding is seen to perform within around 4 dB and 5 5 . dB of TH precoding with perfect CSI for b = 3 and b = 4 respectively.

## V. NUMERICAL RESULTS

In this section we present some numerical results. We assume n T = K = 4 . Here the SNR of the systems is defined to be equal to P .

Fig. 3 shows the average sum rate performance of TH precoding and linear ZF precoding with both perfect CSI and quantized CSI, and 4 8 , and 15 feedback bits per user. We can see TH precoding performs better than linear precoding in both perfect CSI and quantized CSI cases. When the SNR is small and moderate, the average sum rate achieved by quantized CSI-based TH precoding can even be better than that of perfect CSI-based linear ZF precoding.

Fig. 3 . The average sum rate performance of TH precoding and ZF precoding for both perfect CSI and quantized CSI.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0455v1-figure-a6c52847-5e27-4731-a112-6f8ecd948c51.jpg)

Fig. 4 plots the average sum rate loss per user as a function of the number of feedback bits for both ZF precoding and TH precoding in a system at an SNR of 25 dB. We also plot the upper bound from Theorem 1 in this paper and the upper bound from Theorem 1 in [2]. From the figure we can see that nonlinear precoding suffers from imperfect CSI more than linear precoding does. However, the performance of nonlinear precoding can still be better than linear precoding when SNR is not large or the feedback quantization resolution is high enough. In addition, we notice that the upper bound for TH precoding tracks the true rate loss quite closely, and appears to converge faster than the upper bound for linear precoding obtained in [2] as B increases.

## VI. CONCLUSION

In this paper, we have investigated the implementation of TH precoding in the downlink multiuser MIMO systems with quantized CSI at the transmitter side. In particular, our scheme generalized the results in [1] to more general system setting where the number of users K in the systems can be less than or equal to the number of transmit antennas n T . In addition,

Fig. 4 . The average sum rate loss per user and corresponding upper bounds against the number of feedback bits. P = 25 dB.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0455v1-figure-bae8078a-fcd7-4d4b-aa20-8a5f77391e7c.jpg)

we studied the achievable average sum rate of the proposed scheme by deriving expressions of upper bounds on both the average sum rate and the mean loss in sum rate due to CSI quantization. Our numerical results showed that the nonlinear TH precoding could achieve much better performance than that of linear zero-forcing precoding for both perfect CSI and quantized CSI cases. In addition, our derived upper bound for TH precoding converged to the true rate loss faster than the upper bound for zero-forcing precoding obtained in [2] as the number of feedback bits increased.

## APPENDIX A

## PROOF OF Lemma 1

The results for the special cases that K = 1 and n T are trivial. In the following we will consider the cases that 1 &lt; K &lt; n T . Since the user channel vectors in H are unordered, so are the quantized channel direction vectors in ˆ H = [ ˆ h T 1 , · · · , ˆ h T K ] T . According to the QR

decomposition of ˆ H we have

$$\hat { h } _ { k } & = \sum _ { l = 1 } ^ { k } \hat { r } _ { k, l } \hat { q } _ { l }, \\ + \quad \nu \quad \alpha \colon \alpha \dots \dots \colon \dots \dots \dots \dots$$

If we require ˆ r i,i &gt; 0 for i = 1 , · · · , K , this decomposition is unique . Particularly, we have ˆ r 1 1 , = 1 and ˆ q 1 = ˆ h 1 . In addition, ˜ h k is isotropically distributed in the null space of ˆ h k [2]. Thus, for k = 1 we have ˜ h 1 ⊥ ˆ q 1 or equivalently ˜ h 1 is an isotropically distributed unit vector in the null space of ˆ q 1 .

With the assumption of RVQ, the quantized channel direction vectors ˆ h k ( k = 1 , · · · , K ) are independently and isotropically distributed on the n T -dimensional complex unit sphere due to the assumption of i.i.d. Rayleigh fading. Thus we can conclude that the orthonormal basis ˆ q 1 , · · · , ˆ q K of the subspace spanned by quantized channel vectors ˆ h k ( k = 1 , · · · , K ) have no preference of direction, i.e., [ ˆ q T 1 , · · · , ˆ q T K ] T is isotropically distributed in the K × n T semiunitary space. Thus, to derive the distribution of ε 1 = ‖ ˜ h Q 1 ˆ H ‖ 2 , we can assume ˆ q i = e i for i = 1 , · · · , K without loss of generality, where e i is the i -th row of the identity matrix I n T . Recall that ˜ h 1 ⊥ ˆ q 1 , thus the random vector ˜ h 1 can be written in the form of ˜ h 1 = [0 , v ] , where the vector v = [ v , v 1 2 , · · · , v n T -1 ] is isotropically distributed on the ( n T -1) -dimensional complex unit sphere. Then ε 1 = ‖ ˜ h Q 1 ˆ H ‖ 2 = ∑ K -1 l =1 | v l | 2 . Let t l = | v l | 2 . It has been obtained in [25] that the joint p.d.f. of t 1 , · · · , t K -1 is

$$f ( t _ { 1 }, \dots, t _ { K - 1 } ) = \begin{cases} \begin{array} { c c } \frac { \Gamma ( n _ { T - 1 } ) } { \Gamma ( n _ { T } - K ) } \left ( 1 - \sum _ { i = 1 } ^ { K - 1 } t _ { i } \right ) ^ { n _ { T } - K - 1 }, & t _ { i } \geq 0 \text{ for } i = 1, \cdots, K - 1, \sum _ { i = 1 } ^ { K - 1 } t _ { i } = 1 \\ 0, & \text{ otherwise } \end{array}. \end{cases}$$

 Now we want to obtain the distribution of u 1 = ∑ K -1 l =1 t l . We define the following transformation of variables

$$u _ { 1 } = \sum _ { l = 1 } ^ { K - 1 } t _ { l }, \ \ u _ { i } = t _ { i } \quad \text{for $i=2,\cdots,K-1$.} \\ \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \ \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \cdots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \ dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \ Dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \Dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \frown \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \ { \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \$$

It is easy to obtain the corresponding Jacobian is J = 1 . Thus the joint p.d.f. of u , 1 · · · , u K -1 is

$$f _ { u _ { 1 }, \dots, u _ { K - 1 } } \left ( x _ { 1 }, \cdots, x _ { K - 1 } \right ) & = \frac { \Gamma ( n _ { T } - 1 ) } { \Gamma ( n _ { T } - K ) } \left ( 1 - x _ { 1 } \right ) ^ { n _ { T } - K - 1 }. \\ & \quad \,. \.$$

Since 0 ≤ t i ≤ 1 , we have 0 ≤ t 1 = u 1 -∑ K -1 i =2 u i ≤ 1 . The region of the random variables after transformation can be obtained as D = ( u , 1 · · · , u K -1 ) | 0 ≤ K -1 l =2 u l ≤ u 1 ≤ 1 0 , ≤

.

u i ≤ 1 for i = 2 , · · · , K -1 } . Then the marginal distribution of u 1 can be obtained as

$$- \quad \cdot \quad \cdot \quad, \quad \Lambda \quad \cdot \\ f _ { u _ { 1 } } ( x ) \quad = \int \dots \int _ { \mathcal { D } } f ( x, x _ { 2 } \dots, x _ { K - 1 } ) \, \mathrm d x _ { 2 } \dots \mathrm d x _ { K - 1 } \\ = \int \dots \int _ { \mathcal { D } } \frac { \Gamma ( n _ { T } - 1 ) } { \Gamma ( n _ { T } - K ) } \left ( 1 - x _ { 1 } \right ) ^ { n _ { T } - K - 1 } \, \mathrm d x _ { 2 } \dots \mathrm d x _ { K - 1 } \\ \stackrel { ( a ) } { = } \frac { \Gamma ( n _ { T } - 1 ) } { \Gamma ( n _ { T } - K ) } ( 1 - x ) ^ { n _ { T } - K - 1 } \frac { x ^ { K - 2 } } { ( K - 2 )! } \\ \text{ich is given by (16), where in (a) we have used the identity } \, \int \, \dots \int \, \mathrm d t _ { 1 } \cdots \, \mathrm d t _ { n } =$$

which is given by (16), where in (a) we have used the identity ∫ ∫ · · · ∫ ∑ n i =1 t i ≤ h t 1 ≥ 0 , ··· ,t n ≥ 0 d t 1 · · · d t n = h n n ! [24]. We find that ε 1 = u 1 follows beta distribution with shape ( K -1) and ( n T -K ) . In the following we will prove ε k s have the same distribution.

Let π be an arbitrary and channel-independent permutation of (1 , 2 , · · · , K ) . P π = [ 1 π (1) , · · · , 1 π ( K ) ] T is the permutation matrix corresponding to π and 1 π ( ) i is the π ( ) i -th column of identity matrix. We denote ˆ H π = P H π ˆ = [ h T π (1) , · · · , h T π ( K ) ] T the matrix obtained by permutating the row vector of matrix ˆ H according to the permutation π . Then the QR decomposition of ˆ H π can be written as ˆ H π = P RQ π ˆ ˆ = ˆ R Q π ˆ π . With the assumption that ˆ R π has positive diagonal elements, the above QR decomposition of ˆ H π is unique. Using Givens transformation, there is a series of Givens matrices G 1 , · · · , G K -1 ∈ C K × K which satisfy P RG π ˆ 1 · · · G K -1 = ¯ R π [26], where ¯ R π ∈ C K × K is a lower triangular matrix with positive diagonal elements. Since Givens matrix is unitary, we have G 1 · · · G K -1 G H K -1 · · · G H 1 = I . So ˆ H π can be written as ˆ H π = ¯ R G π H K -1 · · · G Q H 1 ˆ . Let ¯ Q π = G H K -1 · · · G Q H 1 ˆ . Then we have ˆ H π = ¯ R Q π ¯ π where ¯ Q π is unitary. Thus ˆ H π = ¯ R Q π ¯ π is also a QR decomposition of ˆ H π . Using the uniqueness of QR decomposition, we conclude that ¯ Q π = ˆ Q π and ¯ R π = ˆ R π . Thus we have

$$\varepsilon _ { k } & = \| \tilde { h } _ { k } \hat { Q } ^ { H } \| ^ { 2 } \\ & = \| \tilde { h } _ { k } \hat { Q } ^ { H } _ { \pi } G ^ { H } _ { K - 1 } \cdots G ^ { H } _ { 1 } \| ^ { 2 } \\ & \stackrel { ( b ) } { = } \| \tilde { h } _ { k } \bar { Q } ^ { H } _ { \pi } \| ^ { 2 },$$

where ( b ) is due to the fact that the matrix G i is unitary for i = 1 , · · · , K -1 . If we let π (1) = k , ˆ h k will be the first row of ˆ H π . According to the previous derivation in the proof, we know ε k for k = 2 , · · · , K -1 have the same distribution as ε 1 whose p.d.f. is give by (16).

## APPENDIX B

## PROOF OF Lemma 2

Let Y = ‖ ˜ h Q k ˆ H ‖ 2 . As shown in Lemma 1, Y follows the beta distribution with shape ( K -1) and ( n T -K ) and the cumulative distribution function (c.d.f.) is given by Pr ( Y ≤ y ) = I x ( K -1 , n T -K ) , where I x ( · , · ) is the regularized incomplete beta function. Using the facts that

$$I _ { x } \left ( a, b \right ) & = \sum _ { m = a } ^ { a + b - 1 } \frac { ( a + b - 1 )! } { m! ( a + b - m - 1 )! } x ^ { m } ( 1 - x ) ^ { a + b - m - 1 }, \\ \intertext { r \infty } \varepsilon \quad \varepsilon \quad \varepsilon \quad \cdot \quad \varepsilon \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdots$$

E [ X ] = ∫ ∞ 0 Pr ( X ≥ x ) d x for nonnegative random variables and binomial expansion, we have

$$\mathbb { E } [ X ] = \int _ { 0 } ^ { \infty } \Pr \left ( X \geq x \right ) d x & \text{ for nonnegative random variables and binomial expansion, we have } \\ \mathbb { E } \left [ - \ln Y \right ] & = \int _ { 0 } ^ { \infty } \Pr \left ( Y \leq e ^ { - x } \right ) d x \\ & = \int _ { 0 } ^ { \infty } I _ { e - x } ( K - 1, n _ { T } - K ) d x \\ & = \int _ { 0 } ^ { \infty } \left \{ \sum _ { m = K - 1 } ^ { n _ { T } - 2 } \frac { ( n _ { T } - 2 )! } { m! ( n _ { T } - 2 - m )! } e ^ { - m x } ( 1 - e ^ { - x } ) ^ { n _ { T } - 2 - m } \right \} d x \\ & = \int _ { 0 } ^ { \infty } \left \{ \sum _ { m = K - 1 } ^ { n _ { T } - 2 } \frac { ( n _ { T } - 2 )! } { m! ( n _ { T } - 2 - m )! } e ^ { - m x } \sum _ { l = 0 } ^ { n _ { T } - 2 } \binom { n _ { T } - m - 2 } { l } ( - 1 ) ^ { l } e ^ { - l x } \right \} d x \\ & = \sum _ { m = K - 1 } ^ { n _ { T } - 2 } \sum _ { l = 0 } ^ { n _ { T } - m - 2 } \frac { ( n _ { T } - 2 )! } { m! ( n _ { T } - 2 - m )! } \binom { n _ { T } - m - 2 } { l } ( - 1 ) ^ { l } \frac { 1 } { m + l }. \\ & = \sum _ { m = K - 1 } ^ { n _ { T } - 2 } \sum _ { l = 0 } ^ { n _ { T } - m - 2 } \frac { ( n _ { T } - 2 )! } { m! \, l! \, ( n _ { T } - m - 2 - l )! } ( - 1 ) ^ { l } \frac { 1 } { m + l }.$$

Thus (17) is proved.

## APPENDIX C

## PROOF OF Theorem 1

First we will prove the fact that | r k,k | 2 and ρ 2 k | ˆ r k,k | 2 have the same distribution. Let the QR decomposition of matrix ˇ H = ΦH ˆ be ˇ ˇ QR . It is easy to see | ˇ r k,k | 2 = ρ 2 k | ˆ r k,k | 2 , where ˇ r k,k is the k -th diagonal element of ˇ R . Since we assume using RVQ, ˇ H has the same distribution as H . Thus ρ 2 k | ˆ r k,k | 2 has the same distribution as | r k,k | 2 .

Using (9), (15) (18) and (19), we can write

$$\text{Using (9), (15) (18) and (19), we can write } \\ \Delta R _ { k } & = \mathbb { E } _ { \mathbb { H }, w } \left \{ \log _ { 2 } \left ( 1 + \frac { P } { \kappa } | _ { r _ { k }, k } | ^ { 2 } \right ) - \log _ { 2 } \left ( 1 + \frac { \frac { P } { \kappa } \rho _ { k } ^ { 2 } | \tilde { r } _ { k, k } | ^ { 2 } \cos ^ { 2 } \theta _ { k } } { \frac { \kappa } { \kappa } \rho _ { k } ^ { 2 } \| \hbar { k } } \right ) \right \} \\ & = \mathbb { E } _ { \mathbb { H }, w } \left \{ \log _ { 2 } \left ( 1 + \frac { P } { \kappa } | _ { r _ { k }, k } | ^ { 2 } \right ) \right \} - \mathbb { E } _ { \mathbb { H }, w } \left \{ \log _ { 2 } \left ( \frac { P } { \kappa } \rho _ { k } ^ { 2 } \left ( | \tilde { r } _ { k, k } | ^ { 2 } \cos ^ { 2 } \theta _ { k } + \| \tilde { h } _ { k } \tilde { Q } ^ { \theta } \| ^ { 2 } \sin ^ { 2 } \theta _ { k } \right ) + 1 \right ) \right \} \\ & \quad + \mathbb { E } _ { \mathbb { H }, w } \left \{ \log _ { 2 } \left ( \frac { P } { \kappa } \rho _ { k } ^ { 2 } \| \hbar { k } \tilde { Q } ^ { \theta } \| ^ { 2 } \sin ^ { 2 } \theta _ { k } \right ) \right \} \\ & \quad \leq \mathbb { E } _ { \mathbb { H }, w } \left \{ \log _ { 2 } \left ( 1 + \frac { P } { \kappa } | _ { r _ { k }, k } | ^ { 2 } \right ) \right \} - \mathbb { E } _ { \mathbb { H }, w } \left \{ \log _ { 2 } \left ( 1 + \frac { P } { \kappa } \rho _ { k } ^ { 2 } | \tilde { r } _ { k, k } | ^ { 2 } \cos ^ { 2 } \theta _ { k } \right ) \right \} \\ & \quad + \mathbb { E } _ { \mathbb { H }, w } \left \{ \log _ { 2 } \left ( 1 + \frac { P } { \kappa } \rho _ { k } ^ { 2 } \| \tilde { h } _ { k } \tilde { Q } ^ { \theta } \| ^ { 2 } \sin ^ { 2 } \theta _ { k } \right ) \right \} \\ & \quad \approx - \mathbb { E } _ { \mathbb { H }, w } \left \{ \log _ { 2 } \left ( \cos ^ { 2 } \theta _ { k } \right ) \right \} + \mathbb { E } _ { \mathbb { H }, w } \left \{ \log _ { 2 } \left ( 1 + \frac { P } { \kappa } \rho _ { k } ^ { 2 } \| \tilde { h } _ { k } \tilde { Q } ^ { \theta } \| ^ { 2 } \sin ^ { 2 } \theta _ { k } \right ) \right \}$$

$$\leq \log _ { 2 } \left ( 1 + \frac { P } { \kappa } \mathbb { E } _ { \mathbf H, W } \left \{ \rho _ { k } ^ { 2 } \| \tilde { \mathbf h } _ { k } \hat { Q } ^ { H } \| ^ { 2 } \sin ^ { 2 } \theta _ { k } \right \} \right ) - \mathbb { E } _ { \mathbf H, W } \left \{ \log _ { 2 } \left ( \cos ^ { 2 } \theta _ { k } \right ) \right \}. \\ \intertext { H a r a t i n d e w a l i m i n t a n d t h a n o n a n o t i v a t i m c i l h. } \| \tilde { \mathbf h }. \hat { \mathbf A } ^ { H \| ^ { 2 } \sin ^ { 2 } A. \ i n t h a c o n d t. }$$

$$& - \stackrel { - \mathcal { M }, \mathcal { N } } { \Phi } \left \{ \stackrel { - \mathcal { M } } { \mathcal { M } } ^ { \mathcal { M } } \stackrel { \mathcal { M }, \mathcal { K }, \mathcal { K } } { \Phi } \right \} \right \} \stackrel { - \mathcal { M }, \mathcal { N } } { \Phi } \stackrel { \mathcal { M }, \mathcal { K } } { \Phi } \stackrel { \mathcal { M } } { \Phi } \right \} \right \} \\ & \quad + \mathbb { E } _ { \mathbf H, W } \left \{ \log _ { 2 } \left ( 1 + \frac { P } { \kappa } \rho _ { k } ^ { 2 } \| \tilde { \mathbf h } _ { k } \hat { Q } ^ { H } \| ^ { 2 } \sin ^ { 2 } \theta _ { k } \right ) \right \} \\ & \quad \approx - \mathbb { E } _ { \mathbf H, W } \left \{ \log _ { 2 } \left ( \cos ^ { 2 } \theta _ { k } \right ) \right \} + \mathbb { E } _ { \mathbf H, W } \left \{ \log _ { 2 } \left ( 1 + \frac { P } { \kappa } \rho _ { k } ^ { 2 } \| \tilde { \mathbf h } _ { k } \hat { Q } ^ { H } \| ^ { 2 } \sin ^ { 2 } \theta _ { k } \right ) \right \} \\ & \quad \leq \log _ { 2 } \left ( 1 + \frac { P } { \kappa } \mathbb { E } _ { \mathbf H, W } \left \{ \rho _ { k } ^ { 2 } \| \tilde { \mathbf h } _ { k } \hat { Q } ^ { H } \| ^ { 2 } \sin ^ { 2 } \theta _ { k } \right \} \right ) - \mathbb { E } _ { \mathbf H, W } \left \{ \log _ { 2 } \left ( \cos ^ { 2 } \theta _ { k } \right ) \right \}. \\ \quad \longrightarrow \quad \dots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots.$$

Here (31) holds by eliminating the non-negative terms ‖ ˜ h Q k ˆ H ‖ 2 sin 2 θ k in the second term of (30). (32) follows by using high SNR approximation and the fact that | r k,k | 2 and ρ 2 k | ˆ r k,k | 2 have the same distribution which has been proved above. (33) follows by applying Jensen's inequality.

Since the norm of the channel vector ρ k and the direction of channel vector ¯ h k are independent and sin 2 θ k and ˜ h k ( ˆ h k ) are also independent with each other [2], we have

$$\mathbb { E } _ { \mathbf H, W } \left \{ \rho _ { k } ^ { 2 } \| \tilde { \mathbf h } _ { k } \hat { Q } ^ { \mathcal { H } } \| ^ { 2 } \sin ^ { 2 } \theta _ { k } \right \} = \mathbb { E } _ { \mathbf H, W } \left \{ \rho _ { k } ^ { 2 } \right \} \mathbb { E } _ { \mathbf H, W } \left \{ \| \tilde { \mathbf h } _ { k } \hat { Q } ^ { \mathcal { H } } \| ^ { 2 } \right \} \mathbb { E } _ { \mathbf H, W } \left \{ \sin ^ { 2 } \theta _ { k } \right \}. \quad \text{(34)} \\ \text{Each the main of this and this of $2\partial\cos_{bosh_{used_recommend{cases}$\theta$} \, \text{$a$} \, \text{$a$} \, \text{$a$} \, \text{$a$} \, \text{$a$} \, \text{$a$} \, \text{$a$} \, \text{$a$} \, \text{$a$} \, \text{$a$} \, \text{$a$} \, \text{$a$} \, \text{$a$} \, \text{$a$} \, \text{$a$} \, \text{$a$} \, \text{$a$} \, \text{$a$} \, \text{$a$} \, \text{$a$} \, \text{$a$} \, \text{$a$} \, \text{$a$ } \, \text{$a$} \, \text{$a$} \, \text{$a$} \, \text{$a$} \, \text{$a$} \, \text{$a$} \, \text{$a$} \, \text{$a$} \, \text{$a$} \, \text{$a$} \, \text{$a$} \, \text{$a$} \, \text{$a$} \, \text{$a$} \, \text{$a$} \, \text{$a$} \, \text{$a$} \, \text{$a$} \, \text{$a$} \, \text{$a$} \, \text{$a$} \, \text{$a$} \. \text{$a$} \, \text{$a$} \, \text{$a$} \, \text{$a$} \, \text{$a$} \, \text{$a$} \, \text{$a$} \, \text{$a$} \, \text{$a$} \, \text{$a$} \, \text{$a$} \, \text{$a$} \, \text{$a$} \, \text{$a$} \, \text{$a$} \, \text{$a$} \, \text{$a$} \, \text{$a$} \, \text{$a$} \, \text{$a$} \, \text{$a$} \, \text{$a$} \, \end{$$

Each term of right hand side of (34) can be obtained respectively as follows.

$$\mathbb { E } _ { \mathbf H, W } \left \{ \rho _ { k } ^ { 2 } \right \} & = \mathbb { E } _ { \mathbf H, W } \left \{ \chi _ { 2 n _ { T } } ^ { 2 } \right \} = n _ { T }, & \\ \int _ { \left \| \tilde { K } _ { \mathbf H } \right \| 2 } \right \| _ { \mathbb { P } } \quad \int _ {$$

$$\mathbb { E } _ { \mathrm H, W } \left \{ \| \tilde { \mathbf h } _ { k } \hat { \mathbf Q } ^ { H } \| ^ { 2 } \right \} & = \mathbb { E } _ { \mathrm H, W } \left ( \text{Beta} \left ( K - 1, n _ { T } - 1 \right ) \right ) = \frac { K - 1 } { n _ { T } - 1 }, \\ & \quad \quad \quad$$

$$\mathbb { E } _ { \mathrm H, W } \left \{ \sin ^ { 2 } \theta _ { k } \right \} & \leq 2 ^ { - \frac { B } { M - 1 } }, \\ \intertext { a b t i n d k i n s o n d f o m i n d f o m i n d f o m i n d f o m i n d f o m i n d f o m i n d f o m i n d f o m i n d f o m i n d f$$

where (36) can be easily obtained by using p.d.f. result in (16) and (37) is given in [18] and [2, Lemma 1] respectively. In [27], the second term in (33) was obtained as

$$\mathbb { E } _ { \mathbf H, W } \left \{ \log _ { 2 } \left ( \cos ^ { 2 } \theta _ { k } \right ) \right \} = \log _ { 2 } ( e ) \sum _ { i = 1 } ^ { n } \binom { n } { i } ( - 1 ) ^ { i } \sum _ { l = 1 } ^ { i ( n _ { T } - 1 ) } \frac { 1 } { l },$$

and (38) is rewritten in [28] as

$$\mathbb { E } _ { \mathbf H, W } \left \{ \log _ { 2 } \left ( \cos ^ { 2 } \theta _ { k } \right ) \right \} = - \frac { \log _ { 2 } ( e ) } { n _ { T } - 1 } \sum _ { i = 1 } ^ { n _ { T } - 1 } \beta \left ( n, \frac { i } { n _ { T } - 1 } \right ). \\ \intertext { n a t r e c t f o$$

The final result follows by combing (34)-(39) .

## APPENDIX D

PROOF OF Theorem 2

The average sum rate for user k can be upper bounded as

$$\text{$ne average sum rate for user $\kappa$ can be upper nounae as} \\ \\ \mathbb { E } _ { \text{H}, w } \left \{ R _ { Q, k } \right \} & = \mathbb { E } _ { \text{H}, w } \left \{ \log _ { 2 } \left ( 1 + \frac { \underline { P } _ { \kappa } \rho _ { k } ^ { 2 } | \hat { r } _ { k, k } | ^ { 2 } \cos ^ { 2 } \theta _ { k } } { \underline { P } _ { \kappa } ^ { 2 } \| \tilde { h } _ { k } \tilde { Q } ^ { H } \| ^ { 2 } \sin ^ { 2 } \theta _ { k } + 1 } \right ) \right \} \\ & \leq \mathbb { E } _ { \text{H}, w } \left \{ \log _ { 2 } \left ( 1 + \frac { | \hat { r } _ { k, k } | ^ { 2 } \cos ^ { 2 } \theta _ { k } } { \| \tilde { h } _ { k } \tilde { Q } ^ { H } \| ^ { 2 } \sin ^ { 2 } \theta _ { k } } \right ) \right \} \\ & = \mathbb { E } _ { \text{H}, w } \left \{ \log _ { 2 } \left ( \frac { | \hat { r } _ { k, k } | ^ { 2 } \cos ^ { 2 } \theta _ { k } + \| \tilde { h } _ { k } \tilde { Q } ^ { H } \| ^ { 2 } \sin ^ { 2 } \theta _ { k } } { \| \tilde { h } _ { k } \tilde { Q } ^ { H } \| ^ { 2 } \sin ^ { 2 } \theta _ { k } } \right ) \right \} \\ & \leq - \mathbb { E } _ { \text{H}, w } \left \{ \log _ { 2 } \left ( \| \tilde { h } _ { k } \tilde { Q } ^ { H } \| ^ { 2 } \sin ^ { 2 } \theta _ { k } \right ) \right \}, \\ \\ \text{where (40) is obtained using the facts that $|\hat { r } _ { k, k } | ^ { 2 } \leq 1$ and $|\tilde { h } _ { k } \tilde { Q } ^ { H } \| ^ { 2 } \leq 1.$ By using [2, Lemma$} \\ \mathbb { A } \stackrel { \sim } { \sim } \mathbb { P } _ { \text{H} } \quad \text{$\mathbb { A } \stackrel { \sim } { \sim } \mathbb { A } \stackrel { \sim } { \sim } \mathbb { A } \stackrel { \sim } { \sim } \mathbb { A } \stackrel { \sim } { \sim } \mathbb { A } \stackrel { \sim } { \sim }$$

where (40) is obtained using the facts that | ˆ r k,k | 2 ≤ 1 and ‖ ˜ h Q k ˆ H ‖ 2 ≤ 1 . By using [2, Lemma 3] and (36), E H , W { R Q,k } can be upper bounded as

$$\mathbb { E } _ { \mathrm H, W } \left \{ R _ { Q, k } \right \} \leq - \mathbb { E } _ { \mathrm H, W } \left \{ \log _ { 2 } \left ( \left \| \tilde { \mathbf h } _ { k } \hat { Q } ^ { H } \right \| ^ { 2 } \right ) \right \} - \mathbb { E } _ { \mathrm H, W } \left \{ \log _ { 2 } \left ( \sin ^ { 2 } \theta _ { k } \right ) \right \}. \quad \, \quad \, ( 4 1 ) \\ \intertext { T h a n d } \mathbb { m } ( \mathcal { U } \colon \, \text{all some $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machines $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ Machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$Machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ machine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ Engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$Engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engines $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ machine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engines $t$ engines $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engines $t$Engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engines $t$ Engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engines $t$engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engine $t$ engines $t$ legging _ \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \ \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \$$

Then (21) follows by using Lemma 2 and [2, Lemma 3].

## APPENDIX E

THE PROOF THAT THE RHS OF (20) CAN BE BOUNDED BY A DECREASING FUNCTION OF n FOR A FIXED n T

$$\text{Let } \mathcal { J } \coloneqq \frac { \log _ { 2 } ( e ) } { n _ { T } - 1 } \sum _ { i = 1 } ^ { n _ { T } - 1 } \beta \left ( n, \frac { i } { n _ { T } - 1 } \right ). & \beta \left ( n, \frac { i } { n _ { T } - 1 } \right ) \text{ can be written as} \\ \quad \_ \int \quad i \quad \bigvee \quad \Gamma ( n ) \Gamma ( \frac { i } { n _ { T } - 1 } )$$

$$\beta \left ( n, \frac { i } { n _ { T } - 1 } \right ) = \frac { \Gamma ( n ) \Gamma ( \frac { i } { n _ { T } - 1 } ) } { \Gamma \left ( n + \frac { i } { n _ { T } - 1 } \right ) }.$$

By applying Kershaw's inequality for the gamma function [29],

$$\begin{array} { c c c } & & \frac { \Gamma ( x + s ) } { \Gamma \left ( x + 1 \right ) } < \left ( x + \frac { s } { 2 } \right ) ^ { s - 1 }, \forall x > 0, 0 \leq s \leq 1. \\ & & \\. & \colon &. &. \end{array} \quad.$$

With x = n -1 + i n T -1 and s = 1 -i n T -1 , we have

$$\frac { \Gamma ( n ) } { \Gamma \left ( n + \frac { i } { n _ { T } - 1 } \right ) } \leq \left ( n - \frac { 1 } { 2 } + \frac { i } { 2 ( n _ { T } - 1 ) } \right ) ^ { - \frac { i } { n _ { T } - 1 } } \leq \left ( n - \frac { 1 } { 2 } \right ) ^ { - \frac { i } { n _ { T } - 1 } }. \\ \text{Thus, $\mathcal{J}$ can be upper bounded as}$$

Thus, J can be upper bounded as

$$\mathcal { J } \leq \sum _ { i = 1 } ^ { n _ { T } - 1 } \Gamma ( \frac { i } { n _ { T } - 1 } ) \left ( n - \frac { 1 } { 2 } \right ) ^ { - \frac { i } { n _ { T } - 1 } }, \\ \intertext { w a d e d i n d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d o d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d t e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e  d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d  e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d ed e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e s t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t i n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n s t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n s s t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n e s s t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e s s t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t i n t t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t u n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e s s s t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e$$

where the RHS is a decreasing function of n for a fixed n T .

## REFERENCES

- [1] C. Windpassinger, R. F. H. Fischer, T. Vencel, and J. B. Huber, 'Precoding in multiantenna and multiuser communications,' IEEE Trans. Wireless Commun. , vol. 3, no. 4, pp. 1305-1316, Jul. 2004.
- [2] N. Jindal, 'MIMO broadcast channels with finite-rate feedback,' IEEE Trans. Inform. Theory , vol. 52, pp. 5045-5060, Nov. 2006.
- [3] ˙ . I E. Telatar, 'Capacity of multi-antenna Gaussian channels,' Europ. Trans. Commun. , pp. 585-595, Nov.-Dec. 1999.
- [4] G. J. Foschini and J. E. Hall, 'Layered space-time architecture for wireless communication in a fading environment when using multi-element antennas,' Bell Labs Tech. J. , pp. 41-59, 1996.
- [5] P. W. Wolniansky, G. J. Foschini, G. D. Golden, and R. A. Valenzuela, 'V-BLAST: An architecture for realizing very high data rates over the rich-scattering wireless channel,' in Proc. URSI Int. Symposium on Signals, Systems, and Electronics , pp. 295-300, Pisa, Italy 1998.
- [6] T. Haustein, C. von Helmolt, E. Jorswieck, V. Jungnickel, and V. Pohl, 'Performance of MIMO systems with channel inversion,' in Proc. 55th IEEE Veh. Technol. Conf. , pp. 35-39, Birmingham, AL, May 2002.
- [7] M. Joham, K.Kusume, M. H. Gzara, and W. Utschick, 'Transmit Wiener filter for the downlink of TDD DS-CDMA systems,' in Proc. IEEE 7th Symp. Spread-Spectrum Technol., Applicat. , pp. 9-13, Prague, Czech Republic, Sep. 2002.
- [8] C. B. Peel, B. M. Hochwald, and A. L. Swindlehurst, 'A vector-perturbation technique for near-capacity multiantenna multi-user communication -Part I: Channel inversion and regularization,' IEEE Trans. Commun. , vol. 53, no. 1, pp. 195-202, Jan. 2005.
- [9] J. Yang and S. Roy, 'Joint transmitter-receiver optimization for multi-input multi-output systems with decision feedback,' IEEE Trans. Inform. Theory , vol. 40, no. 5, pp. 1334-1347, Sep. 1994.
- [10] B. M. Hochwald, C. B. Peel, and A. L. Swindlehurst, 'A vector-perturbation technique for near-capacity multiantenna multiuser communication - Part II: Perturbation,' IEEE Trans. Commun. , vol. 53, no. 3, pp. 537-544, Mar. 2005.
- [11] R. F. H. Fischer, Precoding and Signal Shaping for Digital Transmission , 1st ed. USA: New York: Wiley, 2002.
- [12] M. Costa, 'Writing on dirty paper,' IEEE Trans. Inform. Theory , vol. 29, no. 1, pp. 439-441, May 1983.

- [13] H. Weingarten, Y. Steinberg, and S. Shamai (Shitz), 'The capacity region of the Gaussian multiple-input multiple-output broadcast channel,' IEEE Trans. Inform. Theory , vol. 52, no. 9, pp. 3936-3964, Sep. 2006.
- [14] M. O. Damen, A. Chkeif, and J.-C. Belfiore, 'Lattice code decoder for space-time codes,' IEEE Commun. Lett. , vol. 4, pp. 161-163, May. 2000.
- [15] H. Harashima and H. Miyakawa, 'Matched-transmission technique for channels with intersymbol interference,' IEEE Trans. Commun. , vol. 20, pp. 774-780, Aug. 1972.
- [16] A. A. D'Amico, 'Tomlinson-Harashima precoding in MIMO systems: A unified approach to transceiver optimization based on multiplicative schur-convexity,' IEEE Trans. Signal Process. , vol. 56, no. 8, pp. 3662-3677, Aug. 2008.
- [17] G. Caire and S. Shamai (Shitz), 'On the achievable throughput of a multi-antenna Gaussian broadcast channel,' IEEE Trans. Inform. Theory , vol. 49, no. 7, pp. 1691-1706, Jul. 2003.
- [18] C. K. Au-Yeung and D. J. Love, 'On the performance of random vector quantization limited feedback beamforming in a MISO system,' IEEE Trans. Wireless Commun. , vol. 6, pp. 458-462, Feb. 2007.
- [19] T. Yoo, N. Jindal, and A. Goldsmith, 'Multi-antenna downlink channels with limited feedback and user selection,' IEEE J. Sel. Areas Commun. , vol. 25, no. 7, pp. 1478-1491, Sep. 2007.
- [20] I. Slim, A. Mezghani, and J. A. Nossek, 'Quantized CDI based Tomlinson Harashima precoding for broadcast channels,' in Proc. IEEE Int. Conf. on Commun. (ICC) , Jun. 2011, pp. 1-5.
- [21] R. D. Wesel and J. M. Cioffi, 'Achievable rates for Tomlinson-Harashima precoding,' IEEE Trans. Inform. Theory , vol. 44, no. 2, pp. 824-831, Mar. 1998.
- [22] U. Erez and R. Zamir, 'Achieving 1 2 / log (1 + SNR ) on the AWGN channel with lattice encoding and decoding,' IEEE Trans. Inform. Theory , vol. 50, no. 10, pp. 2293-2314, Oct. 2004.
- [23] K. K. Mukkavilli, A. Sabharwal, E. Erkip, and B. Aazhang, 'On beamforming with finite rate feedback in multiple antenna systems,' IEEE Trans. Inform. Theory , vol. 50, no. 10, pp. 2562-2579, Oct. 2003.
- [24] I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Series, and Products , 6th ed. New York: Academic, 2000.
- [25] L. Sun and M. R. McKay, 'Eigen-based transceivers for the MIMO broadcast channel with semi-orthogonal user selection,' IEEE Trans. Signal Process. , vol. 58, no. 10, pp. 5246-5261, Oct. 2010.
- [26] G. H. Golub and C. F. V. Loan, Matrix Computations , 3rd ed. Baltimore: Johns Hopkins Univ. Press, 1996.
- [27] R. Bhagavatula and J. R. W. Heath, 'Adaptive limited feedback for sum-rate maximizing beamforming in cooperative multicell systems,' IEEE Trans. Signal Process. , vol. 59, no. 2, pp. 800-811, Feb. 2011.
- [28] R. Bhagavatula and R. W. Heath, 'Adaptive bit partitioning for multicell intercell interference nulling with delayed limited feedback,' IEEE Trans. Signal Process. , vol. 59, no. 8, pp. 3824-3836, Aug. 2011.
- [29] D. Kershaw, 'Some extensions of W. Gautschi's inequalities for the gamma function,' Math. Comput. , vol. 41, no. 164, pp. 607-611, Oct. 1983.
Liang Sun Ming Lei
Quantized CSI-Based Tomlinson-Harashima Precoding in Multiuser MIMO Systems
This paper considers the implementation of Tomlinson-Harashima (TH) precoding
for multiuser MIMO systems based on quantized channel state information (CSI)
at the transmitter side. Compared with the results in [1], our scheme applies
to more general system setting where the number of users in the system can be
less than or equal to the number of transmit antennas. We also study the
achievable average sum rate of the proposed quantized CSI-based TH precoding
scheme. The expressions of the upper bounds on both the average sum rate of the
systems with quantized CSI and the mean loss in average sum rate due to CSI
quantization are derived. We also present some numerical results. The results
show that the nonlinear TH precoding can achieve much better performance than
that of linear zero-forcing precoding for both perfect CSI and quantized CSI
cases. In addition, our derived upper bound on the mean rate loss for TH
precoding converges to the true rate loss faster than that of zeroforcing
precoding obtained in [2] as the number of feedback bits becomes large. Both
the analytical and numerical results show that nonlinear precoding suffers from
imperfect CSI more than linear precoding does.