AU et al.: ANALYTICAL PERFORMANCE OF MIMO-SVD SYSTEMS IN RICEAN FADING CHANNELS 1317
Letting σ
2
s
|ρ(d)|
2
K+1
σ
2
e
+
1
|ρ(d)|
2
− 1
, t he post-processing
SINR for a particular realization is given by
SINR
i
=
P
i
μ
2
(l
i
)
σ
2
s
m
q=1,q=i
P
q
μ
2
(l
q
)+σ
2
w
(λ
i
+ σ
2
s
).(11)
With equal power allocation, namely P
i
= P
T
/m for all i,
such a realization is simplified as
SINR
i
=
SINR
i
˜
λ
i
(12)
where
˜
λ
i
λ
i
+ σ
2
s
and
SINR
i
=
σ
2
s
m
q=1,q=i
μ
2
(l
q
)
μ
2
(l
i
)
+
mσ
2
w
P
T
μ
2
(l
i
)
−1
. (13)
A. Statistical Characterization of the Post-Processing SINR
To analyze the BER and outage probability performance
of each eigen-subchannel, we require the PDF of SINR
i
for all i, or equivalently, the marginal PDF of a (translated)
ordered eigenvalue
˜
λ
i
of a complex noncentral Wishart ma-
trix. Although there are many results available for complex
central Wishart matrices, there are very few results related
to the noncentral counterpart. Recently, [11] derived new
exact closed-form expressions for marginal CDFs and new
asymptotic first-order expansions for marginal PDFs
1
of the
ordered eigenvalues. In particular, given the first-order ex-
pansion f or the marginal PDF of λ
i
[11, eq. (27)], namely
f
λ
i
(λ
i
), the marginal PDF of
˜
λ
i
can be derived by f(
˜
λ
i
)=
f
λ
i
(
˜
λ
i
− σ
2
s
)
dλ
i
/d
˜
λ
i
. Hence, we generalize their results to
our problem as follows.
Lemma 1: Let ϕ
1
,...,ϕ
r
o
be the non-zero eigenvalues of
K/(K +1)H
u,t
H
H
u,t
with rank r
0
≤ m, and let λ
i
be the
i-th largest eigenvalue of
ˆ
H
t
ˆ
H
H
t
. Then the marginal density
function of
˜
λ
i
= λ
i
+ σ
2
s
is given by
f(
˜
λ
i
)=a
i
˜
λ
i
− σ
2
s
d
i
. (14)
Here, d
i
=(m −i + 1)(n −i +1)−1, and a
i
is given in [11].
In addition to the performance of individual eigen-
subchannels, the average system performance can be derived
by using statistical unordered eigenvalue distributions [7], [8].
In contrast to the first-order expansions for marginal PDFs
of ordered eigenvalues, exact closed-form expressions for
marginal PDFs of unordered eigenvalues are derived in [12].
Thus, given the exact expression for the marginal PDF of
λ
a
[12, Lemma 2], namely, f
λ
a
(λ
a
), the marginal PDF of
˜
λ
a
can be derived by f(
˜
λ
a
)=f
λ
a
(
˜
λ
a
− σ
2
s
)
dλ
a
/d
˜
λ
a
.
Lemma 2: Let φ
m−r
u
+1
,...,φ
m
be the non-zero eigen-
values of KH
u,t
H
H
u,t
with rank r
u
≤ m, and let
λ
a
be an arbitrary (unordered) non-zero eigenvalue of
ˆ
H
t
ˆ
H
H
t
, with a being randomly chosen from {1,...,m}.
Then the marginal PDF of
˜
λ
a
= λ
a
+ σ
2
s
is given
by (15) where κ = e
−
m
l=m−r
u
+1
φ
l
m(τ!)
m−1
m−r
u
l=1
(l −
1)!
m
l=m−r
u
+1
φ
m−r
u
l
m−r
u
+1≤k<l≤m
(φ
l
− φ
k
)
−1
, τ =
1
It is worth noticing that exact marginal PDFs can be trivially obtained
from the exact marginal CDFs by differentiation [9]. However, the resulting
expressions, which are in terms of Nuttall Q-functions [13], are very cum-
bersome to work with.
n − m,
˜
D
i,j
is the (i, j)-th cofactor of an m × m matrix
˜
A
shown in (16) with [s]
t
=(s + t − 1)!/(s − 1)! being the
Pochhammer symbol and
p
F
q
(s
1
,...,s
p
; t
1
,...,t
q
; x) being
the generalized hypergeometric function.
III. B
IT ERROR RAT E ANA LYS I S
For an M-ary constellation, the conditional BER of the i-
th eigen-subchannel is given in the following generic form.
P
E
i; μ(l
i
)
= E
˜
λ
i
C
1
(M)
v
2
v=v
1
c
2
c=c
1
C
2
(v, c, M )
Q
C
3
(c, M)
SINR
i
˜
λ
i
(17)
where v
1
, v
2
, c
1
, c
2
, C
1
(M), C
2
(v, c, M ), C
3
(c, M) are
constellation-dependent parameters and Q(·) is the Gaussian
Q-function. The main advantage of using (17) is its ease
of application to various kinds of modulations such as M-
ary quadrature amplitude modulation (MQAM) and
√
M-ary
pulse amplitude modulation (PAM) [14]. However, for M -
ary phase shift keying (MPSK), where M>4, the exact
expression in (17) reduces to a tight closed-form approxima-
tion [15]
2
. Given (14), we now analyze the BER performance
of an individual eigen-subchannel as follows
3
.
Theorem 1: The BER of the i-th largest eigen-subchannel
of MIMO-SVD systems in Ricean fading, with channel es-
timation error and feedback delay, is given by the following
closed-form high-SINR approximation
P
(∞)
E
(i) ≈
1
M
M
l
i
=1
1
M
m−1
M
l
q
=1,
q=1,...,m,q=i
P
(∞)
E
i; μ(l
i
)
(18)
with P
(∞)
E
i; μ(l
i
)
being the high-SINR conditional BER
approximation
P
(∞)
E
i; μ(l
i
)
≈
a
i
2
d
i
Γ(d
i
+
3
2
)
√
π(d
i
+1)
C
1
(M)
v
2
v=v
1
c
2
c=c
1
C
2
(v, c, M )
C
3
(c, M)
SINR
i
K +1
−d
i
−1
(19)
where Γ(z)=
∞
0
x
z−1
e
−x
dx is the Gamma function [16].
Proof: Given the first-order expansion for the marginal PDF
of an ordered eigenvalue in Lemma 1, the high-SINR con-
ditional BER approximation (19) is obtained by using [20,
Proposition 1] along with [11, Section I V.A] and Appendix III.
Then, (18) can be obtained by un-conditioning (19) over the
amplitudes of constellation symbols of all m − 1 cochannel
interferences followed by averaging over μ(l
i
) of all M
constellation symbols of the i-th eigen-subchannel.
In addition to the individual eigen-subchannels, it is also
important to analyze the average system performance so as
to gain a more thorough insight. It has been shown in [17]
that a significant performance improvement can be obtained
2
For M =2and 4, i.e., BPSK and QPSK, exact expressions can be
obtained by considering them as binary PAM and 4QAM, respectively.
3
It is noteworthy that unequal power allocation can be straight-forwardly
employed in Theorem 1. However, in order to have a fair comparison with the
average BER derived by using marginal unordered eigenvalue distributions in
Theorem 2, we restrict the analysis to systems with uniform power allocation.