HW 3
HW 3
Question 1) We know that the capacity of a MIMO system with Gaussian iid channel
entries under fast fading and CSIR is,
C = E log(det(I + P
1
HH )),
Nt
(1)
where the expectation is over the fading matrix H. Let us assume that Nr Nt . Our
discussion applies to all kind of matrices with iid entries.
(a) Assume that H has arbitrary distribution, but fast-fading with CSIR is the underlying
model. What is the achievable rate if we still employ independent Gaussian distributions
at the transmit antennas.
(b) From (1),
Nr
C = E log(1 + P i ),
i=1
where the expectation is over the eigen-values (random) of N1t HH . We will now show that
C
Nr is easily evaluated in the limiting sense, i.e. by assuming Nr and Nt to be large. To this
end, denote
1 Nr
FH (x) =
1{i x} .
Nr i=1
Compute Fh ().
(c) Recall the idea of Riemann integral in calculus. Assume that the maximal Eigen value
of N1t HH , say max , is bounded almost surely (mild technical assumption), thus 0 i .
into n intervals, and call the lth interval Il . Argue that for any
Cut the interval [0, ]
continuous function g(),
Nr
n Nr
g(i ) g( l)1{i Il } ,
n
i=1
l=1 i=1
and
n
1 Nr
which is called the limiting cdf of the eigen values. This will tell us the asymptotic fraction
of Eigen values which lie below a given number. Argue that
max
1 Nr
g(x)dF (x)
g(i )
Nr i=1
0
(e) In the special case where Nr = Nt , the limiting Eigen-value distribution of N1t HH has a
very appealing form, a result shown by Marcenko and Pastur, known as the quarter-circle
law. In particular
x
F (x) = f (u)du,
0
where
1
4 u2 , 0 u 2.
(f ) When Nr Nt , we can still have a closed form characterization of f (u), given by the
Marcenko-Pastur law (we do not state it here, amble references are available online.) In
general, show that
max
C
log(1 + P )f ()d.
Nr
min
Thus we can compute the rate per receive antenna by simply evaluating the line-integral.
Compare this with the complexity of integrating over the fading matrix.
C
for m m Gaussian
(g) Is the above formula good? Let us get a surprise by evaluating m
random matrices when m = 1, 2, . We expect the value to converge when m gets to a few
hundreds or thousands. Here is the plot with m, where P = 100 and the channel entries
are normalized Gaussians.
C
m
3.5
2.5
2
0
20
40
60
80
m
100
120
C
m.
140
160
Question 2) You have to submit Exercise 8.6 of the book by Tse and Viswanath, which
(2)
x
h + z = x + z.
h
This is known as channel inversion. For the above model, we expect the same error
probability behaviour as an AWGN. If this surprises you, then you have some thinking to
do for the remaining part of the question.
Find the flaw in the argument which guarantees an exponentially decaying error probability by channel inversion.
Hint: Watch your transmissions, assuming that you are allowed to spend an average
transmit-power of P
(b) Consider the model in the previous part and assume that the fading distribution is that
of max h1 , h2 , where hi , i = 1, 2 is normalized complex Gaussian. Can we now achieve an
exponentially decaying error probability with SNR by channel inversion.
Question 6) Consider a 3 3 MIMO flat fading system, The communication model for
every channel use is y = Hx + z, with the fixed fading matrix given by
1
0
1 0
61
H=
4 2
5
6+1
0 4 2
6+1
4 2
61
4 2
Let the noise process be distributed as Nc (0, Q), where Q is the noise covariance matrix.
Consider an average transmit power constraint of
Ex2 100.
Find the capacity when
1 0 0
Q = 0 1 41
1
0
4 1