16 Journal of Communications and Information Networks
section VI summarizes the main conclusions of this paper.
Notations: Throughout the paper, boldfaced letters are used
to represent vectors and matrices. The operator (·)
H
and (·)
∗
stand for the Hermitian and conjugate of a matrix, respec-
tively. We let I
N
denotes an N × N identity matrix, and
diag(A) denotes the diagonal matrix with the same diago-
nal elements as matrix A. The notation E
{
·
}
indicates the
expectation operator, whereas
k
·
k
stands for the Euclidean
norm. We use C N (·,·) to denote a multi-variate circularly-
symmetric complex Gaussian distribution.
II. SYSTEM MODEL
Consider the uplink of a single-cell multi-user MIMO sys-
tem formed by a base station (BS) equipped with an array of
M antennas and serving N single antenna user terminals in
the same time-frequency resource (see Fig. 1). The hexag-
onal cell is assumed to have radius of r
c
, and the users are
distributed uniformly and randomly over it with the exclusion
of a central disk of radius r
h
. We assume independent and
identically distributed (i.i.d.) Rayleigh block fading over T
symbols, which means that the channels are considered static
within each time-frequency coherent block and have indepen-
dent realizations in different blocks.
Figure 1 Illustration of a multi-user MIMO scenario: a BS is equipped with
an array of M antennas and simultaneously serves N randomly distributed
single-antenna user terminals in the same time-frequency resource
Single-carrier, narrowband transmission is considered, and
thus the channel matrix between BS and users are represented
by an M ×N matrix G = HD
1/2
, where H is the M ×N
matrix modeling small-scale fading coefficients between the
users and BS, and D is the N ×N diagonal matrix represent-
ing the large-scale fading (both path loss and shadow fad-
ing) coefficients. The (m, n) element of G can be modeled as
g
mn
= h
mn
p
β
n
, with h
mn
describing the instantaneous prop-
agation channel from the nth user to the mth antenna of BS
and β
n
being the large-scale fading component. We assume a
Rayleigh small-scale fading distribution such that h
mn
is zero
mean, circularly symmetric, complex Gaussian random vari-
able with variance 1. Besides, β
n
is assumed to be constant
across the antenna array
[13]
.
Linear processing is used at a BS to perform uplink data
detection. Under imperfect CSI assumption, the UEs send
orthogonal pilot sequences of length τ symbols to the BS
for channel estimation, which indicates that (T −τ) symbols
are left to enable uplink data transmission. The entire uplink
transmission process consists of the following three phases.
• First, the UEs send pilot sequences which are spatially
orthogonal and τ symbols long. As seen from Fig. 2, uplink
pilot transmission takes place first and occupies τ symbols.
• Second, based on the received pilot sequences at the BS,
channel estimation is performed for linear receiver processing.
In this paper, we ignore the processing time needed for this
phase.
• Third, all terminals send their uplink data to the BS using
the remaining (T −τ) symbols, which is shown in Fig. 2.
uplink pilot
uplink
transmission
τ symbols (T − τ) symbols
coherent time: T symbols
Figure 2 Illustration of the uplink transmission: uplink pilot signaling takes
place first and occupies τ symbols; the subsequent (T −τ) symbols are used
for uplink data transmission
A. Uplink Training and Data Transmission
In the uplink training phase, the UEs transmit a sequence
of pilot signals which can be grouped into a matrix Ψ =
[ψ
1
,··· ,ψ
i
,··· ,ψ
N
] ∈ C
τ×N
with each element ψ
n
= [ψ
n
] ∈
C
τ×1
being the pilot vector for the nth user and τ represent-
ing the number of training symbols. Assuming that each user
terminal transmits the pilot symbol with equal transmit power
p
t
, the signal of ith symbol at the receive antennas z
i
∈ C
M×1
is written as
z
i
=
√
τ p
t
N
∑
n=1
ψ
ni
g
n
+ w
i
, (1)
where g
n
denotes the nth column of the channel matrix G,
ψ
ni
is the ith element of ψ
n
, and w
i
represents the additive
white Gaussian noise (AWGN) vector with zero mean and unit
element-wise variance written as w
i
∼ C N (0, I
M
).
In practical systems, an automatic gain control (AGC) vari-
able gain amplifier is usually used before the ADC to adjust
the dynamic range of the received signal z
i
. We assume the
gain of AGC is set appropriately, and by adopting the additive
quantization noise model, the output of the quantizer can be
formulated as
ˆz
i
= αz
i
+ ˆw
i
= α
√
τ p
t
N
∑
n=1
ψ
ni
g
n
+ αw
i
+ ˆw
i
, (2)
where α = 1 −ρ with ρ being the inverse of the signal-to-
quantization-noise ratio, and ˆw
i
is the additive Gaussian quan-