1

Enhancing massive MIMO: A new approach
for Uplink training based on heterogeneous
coherence times
arXiv:1707.03571v2 [cs.IT] 6 Oct 2017

Salah Eddine Hajri∗ , Mohamad Assaad∗ , Maialen Larrañaga∗
*TCL Chair on 5G, Laboratoire des Signaux et Systemes (L2S, CNRS),
CentraleSupelec
Gif-sur-Yvette, France
{Salaheddine.hajri, Mohamad.Assaad, Maialen.Larranaga}@centralesupelec.fr

Abstract
Massive multiple-input multiple-output (MIMO) is one of the key technologies in future generation
networks. Owing to their considerable spectral and energy efficiency gains, massive MIMO systems
provide the needed performance to cope with the ever increasing wireless capacity demand. Nevertheless,
the number of scheduled users stays limited in massive MIMO both in time division duplexing (TDD) and
frequency division duplexing (FDD) systems. This is due to the limited coherence time, in TDD systems,
and to limited feedback capacity, in FDD mode. In current systems, the time slot duration in TDD mode
is the same for all users. This is a suboptimal approach since users are subject to heterogeneous Doppler
spreads and, consequently, different coherence times. In this paper, we investigate a massive MIMO
system operating in TDD mode in which, the frequency of uplink training differs among users based
on their actual channel coherence times. We argue that optimizing uplink training by exploiting this
diversity can lead to considerable spectral efficiency gain. We then provide a user scheduling algorithm
that exploits a coherence interval based grouping in order to maximize the achievable weighted sum
rate.

2

I. I NTRODUCTION
Owing to the exponential increase in data traffic and to the emergence of data-hungry applications with a large number of connected devices, 5G networks need to be able to cope with huge
throughput demand. Several technologies have been considered in order to improve the wireless
networks performance with one noticeable example, namely, massive MIMO Marzetta [1]. A
typical massive MIMO base station (BS) is equipped with a large number of antennas that allow
the spatial multiplexing of a large number of user devices.
Massive MIMO systems have been intensively investigated and have shown to have an interesting potential in improving both the spectral and energy efficiencies of wireless networks.
Although promising, massive MIMO requires accurate channel state information (CSI) estimates
in order to improve the network performance. Typically, in a TDD system, CSI is estimated using
uplink training with orthogonal pilot sequences. Due to the limited coherence interval, these
sequences are reused which results in pilot contamination, see Ashikhmin et al. [2]. Another
reason for CSI inaccuracy, which has received less attention, is channel aging. This phenomenon
refers to variation of the channel between the instant when it is learned, via uplink training, and
the instant when it is applied for signal decoding in the uplink and precoding in the downlink.
The change of the wireless channel between the BS and the mobile user is mainly due to
the mobility of the latter. The impact of channel aging has been studied in a MIMO system
with coordinated multi-point transmission/reception (CoMP) in Thiele et al [3]. The authors
showed that there is hardly any difference between ideal and delayed feedback when utilizing
channel prediction filters with low mobility. The authors in Truong et al. [4] have analyzed the
achievable rate in uplink and downlink in the presence of channel aging in massive MIMO.
They have proposed a channel prediction scheme that overcomes the performance degradation
due to user mobility. In Papazafeiropoulos et al. [5], [6], the effect of channel aging combined
with channel prediction has been investigated in scenarios with regularized ZF precoders and
minimum-mean-square-error (MMSE) receivers, respectively. In Kong et al. [7], lower bounds
on the achievable sum rate with channel aging and prediction have been obtained for arbitrary
number of BS antennas and users. Moreover, the power scaling law has been derived for the
single cell downlink and multi cell uplink scenarios. The scenario we consider in this paper is
different to those proposed in the citations above, in that users in our model are assumed to
have different speeds. This means that the channel aging effect is not the same for all mobile

3

users. More related to our work is Vu et al. [8]. The authors in this paper have proposed efficient
multiuser models for massive MIMO systems that exploit the different user velocities in order
to optimize uplink training. The main idea is to allow users with faster changing channels to
send uplink training sequences with higher periodicity than users with low velocity.
In this paper, we aim at increasing the number of scheduled users per slot without incrementing
training resources. Users experience different Doppler spreads, it is therefore sensible to adapt
the frequency of their CSI estimation accordingly. By doing so we reduce the need to estimate
the same CSI at each time slot (for slow users). As a result, the network is able to serve more
users while increasing the achievable spectral efficiency. The training method we propose uses
an estimate of the Doppler spread or, equivalently, the autocorrelation of the channel in order
to acquire a channel estimate that can be reused in future slots. An estimate of the channel
autocorrelation can easily be obtained by exploiting the cyclic prefix in an OFDM symbol
[9]. In our scheme, before uplink training, copilot users are selected based on their channel
autocorrelation coefficients. With copilot users having comparable Doppler spreads, the impact
of channel aging on the achievable spectral efficiency is reduced. This enables the system to reuse
the same estimated CSI. While in Vu et al. [8], an equal average achievable rate for all users
is considered, we proceed differently by deriving a tight lower bound on the achievable rates
in a massive MIMO system when using the coherence-time based training briefly introduced
above. The derived expression provides more insights into the impact of an adaptive training
procedure where users have heterogeneous CSI delays. Using the derived lower bound, we are
able to provide a scheduling optimization algorithm which achieves a higher weighted sum-rate
by exploiting Doppler spreads diversity. We formulate a combinatorial optimization problem that
determines, for each copilot user group, whether it needs to perform uplink training or if it can
be processed using outdated CSI. We then propose a low complexity approximation algorithm
that achieves at least a (4 + α)-approximation of the optimal weighted sum rate.
This paper is organized as follows. We describe the considered system model in Section II.
We discuss the advantages of coherence-time based uplink training in Section III. We propose a
coherence time based scheduling algorithm in Section IV. Finally, in Section V, numerical results
are presented.

4

II. S YSTEM M ODEL A ND P RELIMINARIES
We consider a two-dimensional hexagonal cellular network composed of C cells with one
macro base station (BS) and K users in each cell. Each BS is equipped with M antennas and
each user has a single omni-directional antenna. All users in the network have heterogeneous
movement directions and speeds. This means that they experience heterogeneous Doppler spreads
and, consequently, different wireless channel correlation across time slots. We assume that the
channel of each user can be decomposed as a product of large and small scale fading coefficients.
The wireless channel from the ith user in the bth cell to the j th BS, at time slot t, is given by:
q
[j] [j]
[j]
(1)
gib (t) = βib hib (t).
[j]

Under flat Rayleigh fading, hib ∈ CM ×1 follows a circular-symmetric complex Gaussian dis[j]

tribution hib ∼ CN (0, IM ), where IM is the identity matrix of size M × M . The large-scale
[j]

fading coefficient βib ∈ R+ depends on the pathloss and the shadowing. We suppose that the
large-scale fading coefficient stays constant for long coherence blocks and varies much slower
than small-scale fading. The time is slotted with a duration of Ts units for each time slot. Ts is a
system parameter that depends on the maximum Doppler spread supported by the network [11].
We consider a TDD mode where the entire frequency band is used for downlink and uplink
transmissions by all base stations and users. CSI estimation is performed using orthonormal
training sequences in the uplink. The same sets of pilot sequences are used in all cells. For τ
scheduled users for uplink training, in each cell, the orthonormal training sequences qi ∈ Cτ ×1
(qi† qj = δij ) are used. In practice, processing delays and user movement causes a variation
in the wireless channel between the instants when it is learned and used. This phenomenon
is referred to as channel aging. In order to take into consideration the autocorrelation of the
channel between two time instants, we adopt a time varying model for the wireless channel.
We consider a stationary ergodic Gauss-Markov block fading regular process (or auto regressive
model of order 1), where the channel remains constant for a slot duration and changes from slot
to slot according to the following equation
q


[j]
[j] [j]
[j]
[j]
gib (t) = βib ρib hib (t − 1) + εib (t) ,

(2)

[j]

where εib (t) denotes a temporally uncorrelated Gaussian noise with zero mean and variance
[j]2

[j]

(1 − ρib )IM . The temporal correlation parameter of the channel ρib is given by Jakes [10]:
[j]

[j]

ρib = J0 (2πfib Ts ),

(3)

5

[j]

where J0 (.) is the zeroth-order Bessel function of the first kind and fib represents the maximum
Doppler shift of the ith user in cell b with respect to the antennas of the j th BS. In previous
works on the channel aging issue, the same maximum Doppler shift has been considered for
all users. In our work, we adopt a more realistic setting in which, mobile users have different
frequency shifts since we consider different movement velocities and directions. For every user
i, b the maximum Doppler shift with respect to the antennas of the j th BS is given by:
[j]

fib =

νib fc
[j]
cos(θib ),
c

(4)

where νib is the velocity of user i, b in meters per seconds, c = 3 × 108 mps is the speed of light,
[j]

fc is the carrier frequency and θib represents the angle between the movement direction of the
mobile device and the direction of the incident wave. The channel autocorrelation is bounded
[j]

as, 0 ≤ ρib ≤ 1. The BS estimates the channels of its associated users thanks to uplink pilot
training. In this paper, MMSE channel estimation is used. The wireless channel of user i, b, at
slot t, can be written as follows:
[j]
gib (t)

=

[j]
ρib ĝib (t

− 1) +

[j]
ρib g̃ib (t

q
[j] [j]
− 1) + βib εib (t),

(5)

[j]

[j]

where
 ĝib (t − 21) is the MMSE channel estimate at time t − 1 which has a ĝib (t − 1) ∼
[j]
β
CN 0, 1 PibC [j] IM distribution with Pp the uplink pilot transmit power. The channel esti+ l,l6=b βil
Pp
 
 
[j]2
βib
[j]
[j]
[j]
mation error at time t−1 is denoted g̃ib (t−1) with g̃ib (t−1) ∼ CN 0, βib − 1 PC [j] IM .
Pp

Note that

[j]
ĝib (t

− 1) and

[j]
g̃ib (t

+

l,l6=b

βil

− 1) are mutually independent.

III. A NEW COHERENCE TIME BASED APPROACH FOR USER S CHEDULING I N M ASSIVE
MIMO S YSTEMS
The current massive MIMO models consider a fixed coherence interval Ts for all users in
the network. The duration of Ts is based on the maximum Doppler spread that the network can
handle [11]. Sidestepping the important fact that, users experience heterogeneous Doppler spreads
and, consequently, different coherence intervals, can only result in suboptimal exploitation of
the available resources. In fact, requiring users to perform uplink CSI estimation at the same
frequency without taking into consideration their actual channels autocorrelation, results in a
useless redundancy of information. If we take into consideration the channel aging effect between
time slots, we notice that, whenever the variation of channel is not important, we may reuse
outdated CSI estimates in order to reduce the need for further training. In practice, this means

6

that the network is not required to estimate the channels of all scheduled users at each time slot.
A more appropriate approach is to estimate the CSI of scheduled users whenever needed, as a
function of its actual coherence interval. The idea of reusing an already estimated information can
be quite tempting. It frees part of the training resources. These can be used for data transmission
or to schedule more users. In both cases, we can observe an increase in the achievable spectral
efficiency. Nevertheless, allowing the network to reuse the CSI from past slots can not provide
the needed capacity increase if no prior conditions on user selection is imposed. In this section,
we investigate the achievable rate with outdated CSI. We also derive an important condition
which ensures that a coherence-time based training scheme is able to provide a considerable
increase in the network performance.
A. An adaptive coherence time based user scheduling
We consider a massive MIMO system, where the network is not required to estimate the CSI
of all scheduled users at each time slot. The only requirement is that copilot users have the same
CSI delay. This guarantees a proportional aging effect for copilot users. The latter provides a
performance gain that will be explained later in this section. For all users in the same copilot
group to have the same CSI delay, we consider that all users in the same copilot group λg for
g = 1, ..., Ng are either scheduled for uplink training synchronously, and they have the same CSI
delay dg , g = 1, ..., Ng or are not scheduled at all. At each time slot, the network schedules all
covered users for data transmission, and a maximum of τ < Ng copilot groups are allowed to
perform uplink training. The rest will be using the last estimated version of their CSI. Instead of
requiring the scheduled users to perform uplink training with the same periodicity, we propose
a new Time Division Duplexing protocol where, the training frequency depends on the actual
coherence time. The proposed adaptive coherence-time based user scheduling procedure goes as
follows (see also Figure 1).
1) In the beginning of each large-scale coherence block, the network estimates the large scale
[j]

[j]

fading and channel autocorrelation coefficients, i.e., βib and ρib for all i = 1, . . . , K, and
b, j = 1, . . . , C.
2) Next, the network constructs Nc user clusters, based on the autocorrelation coefficients,
using the K-mean algorithm with K = Nc , see Young et al [12]. Each cluster will be
characterized by an average autocorrelation coefficient or, equivalently, an average Doppler

7

spread. Defining the number of clusters Nc is of paramount importance. In this work, we
choose to define Nc as
Nc = d

Tmax
e,
Ts

(6)

where Tmax represents the maximum coherence time experienced by the users. Defining
Nc as in (6) is motivated by the need to have an average coherence slot per cluster that
approach a multiple of Ts . Doing so guarantees a proper definition of the periodicity of
CSI estimation as a function of Ts .
3) Next, all users in the network (K per cell) will be allocated to Ng copilot user groups.
Each group contains at maximum C users from the same channel autocorrelation cluster
and from different cells.
4) The network then schedules τ copilot users for uplink training synchronously. Selecting the
users to be scheduled for uplink training can be done using various types of scheduling
algorithms depending on the desired performance criteria. The algorithm we propose is to
be found in Section IV.
Copilot groups from cluster 1
Copilot groups from cluster 2
Copilot groups from cluster 3

Figure 1: Example of the adaptive coherence-time based user scheduling procedure with C =
4, K = 2, Ng = 6, Nc = 3. Left: illustration of Step 2). Right: illustration of Step 3).

B. Spectral efficiency with outdated CSI
We now investigate the achievable spectral efficiency when the aforementioned procedure is
applied. For the sake of analytical traceability, we consider that all copilot groups contain exactly
C users. We, henceforth, refer to each user by its copilot group and serving BS indexes g, l.
During uplink data transmission, the j th BS receives the following data signal at time slot t:
Yu[j] (t)

Ng
C X
X
p
[j]
=
Pu ggl (t)Sgl + Wu (t),
l=1 g=1

(7)

8

where Pu denotes the reverse link transmit power, Wu (t) ∼ CN (0, IM ) is the additive noise
and Sgl denotes the uplink signal of the user from copilot group g in cell l. In what follows we
derive lower bounds on the achievable sum rate with a matched filter receiver when outdated
CSI is used. At the reception, the j th BS applies a matched filter receiver based on the latest
version of estimated CSI for the scheduled users. The resulting average achievable sum rate in
the system is given in the following theorem. The proof can be found in Appendix A in [13].
Theorem 1. The network serves Ng copilot groups, τ of which are scheduled for uplink training.
Using a matched filter receiver based on the latest available CSI estimates for each group, the
average achievable sum rate Ru in the uplink is lower bounded by:
!

Ng 
[l]2 [l]2dg
C X
X
(M − 1)βgl ρgl
τ
,
log 1 +
Ru ≥
1−
p
n
T
(M
−
1)
×
I
+
I
s
gl
gl
l=1 g=1

(8)

where dg , g = 1...Ng represents the CSI delays of users using the same pilot sequence. Iglp and
Igln are given by:
Iglp =

C
X

2dg

ρ[l]
gc

2

[l]
βgc
,

(9)

c6=l

Igln

C X
C
X
X
2dg
[l]
[l]
=(
βkc +
(βgc
− ρ[l]
gc
c=1 k6=g

+

1
1
)×( +
Pu
Pp

c=1
C
X

[l]

βgb ).

[l]2

1
Pp

+

βgc
PC

[l]

)

b=1 βgb

(10)

b=1

Equation (8) provides further insights into the evolution of the achievable sum rate as a function
of the CSI time offset. The achievable sum rate for a given copilot user group decreases as the
time offset increases since the correlation between the estimated CSI and the actual channel fades
over time. In order to investigate the potential gain that the proposed approach can provide, we
compare it with a reference model in which all of the Ng copilot groups take part in the uplink
training at each time slot, as in classical TDD protocols. This comparison is done in the next
theorem where we consider the asymptotic regime (as M grows large). The proof can be found
in Appendix B in [13].
Theorem 2. In the asymptotic regime, the proposed training approach enables to improve the
achievable spectral efficiency of each user if the following condition is satisfied ∀ λg , g =

9

1, ..., Ng :
[min]2

ρ̄g

[max]2



!dg

1+

≥

[∞]
SINRg,l

g
 TTs −N
−τ
s

[∞]

ρ̄g

−1
,

(11)

SINRg,l

with
[l]2

[∞]
SINRg,l
[min]

and ρ̄g

[max]

and ρ̄g

=P

βgl

[l]2

,

(12)

b6=l βgb

the minimum and maximum channel autocorrelation coefficients in group

g, respectively.
Condition (11) ensures that the achievable spectral efficiency of all users is improved when
outdated CSI is used. We can see that the gain in spectral efficiency is maintained as long as
the SINR degradation over time is compensated by the spared resources from uplink training.
It also implies that the ratio between the maximum and minimum autocorrelation should be
bounded. The bound becomes tighter as the allowed CSI delay increases. This condition is quite
intuitive. It means that copilot users are required to experience comparable channel aging effects
in order to be able to improve the achievable spectral efficiency using less uplink training. In
order for Condition (11) to be satisfied we need users from the same copilot group to have
equivalent coherence times and, consequently, they proceed to uplink training with the same
periodicity. This is the reason why in Step 2) of the scheme proposed in Section III.A we group
the mobile users according to channel autocorrelation coefficients. Next, we provide an efficient
user scheduling method for uplink training. It exploits the autocorrelation based grouping in
order to improve the achievable weighted sum rate. We consider that all copilot groups are
active during data transmission. The network is required to select which groups will refresh
their CSI by scheduling them for training.
IV. I MPROVING THE ACHIEVABLE WEIGHTED SUM RATE
Assuming a constant coherence interval for all users results in a suboptimal resources exploitation since in practice, users experience heterogeneous Doppler frequencies. We showed,
in Theorem 2, that allowing the network to use outdated CSI estimates can actually improve
the achievable spectral efficiency. In this section, we study the impact of the proposed adaptive
training scheme on the achievable weighted sum rate. In this case, scheduling users for uplink
training does not depend solely on the estimated CSI but also on the users traffic pattern. In

10

order to take into consideration the QoS, we consider optimizing the weighted sum of rates. The
weights enable to introduce fairness in the scheduling procedure [14]. We denote the weights
associated with the users by wgl , ∀g = 1, ..., Ng , l = 1, ..., C.
The aim here is to schedule copilot groups for uplink training in order to maximize the
achievable weighted sum rate. As in previous sections we assume that all users are scheduled
for data transmission and τ < Ng copilot groups are selected to update their CSI. Grouping
users based on their autocorrelation coefficients simplifies the scheduling problem. Instead of
deciding, at each time slot, to which user each pilot sequence is going to be allocated, decisions
are made on predefined collections of copilot users. We define the vector Y = (y1 , . . . , yNg ),
with yg , ∀g = 1, ..., Ng given by


1 if group g is scheduled for training,
yg =

0 otherwise.

(13)

The copilot group scheduling problem when outdated CSI use is permitted, is formulated as
follow
PNg
max(1 −
Y

g=1

Ts

Ng
C X
yg X
~ Y ),
wgl Rgl (d,
)

(14)

l=1 g=1

Ng

subject to

X

yg ≤ τ,

(14a)

g=1

~ Y ) is the lower bound on the achievable sum rate for the user in copilot group g
where Rgl (d,
and cell l (see Eq. (8)), and d~ is the delay vector. Here, (14a) captures the fact that the number
of scheduled users for uplink training in each cell, is at maximum equal to τ . In the following
theorem, we show that problem (14) is equivalent to maximizing a submodular set function
subject to a cardinality constraint. This formulation enables to apply a low complexity algorithm
that solves the problem. The proof can be found in Appendix C in [13].
Theorem 3. Problem (14) is equivalent to maximizing a submodular set function with a cardinality constraint.
In order to solve problem (14), we use the Submod-Max-Cardinality algorithm in Gupta et
al. [15]. This algorithm starts by running the commonly used greedy approach, see Nemhauser
et al. [16]. This provides a first candidate solution S1 . The approximation algorithm of Feige et
al. [17] is then used on S1 in order to obtain a second candidate solution S2 . The same steps

11

are repeated on the remaining set of copilot groups that were not selected in S1 . The algorithm
returns then the best found solution. When used to solve a non-monotone submodular function
maximization subject to a cardinality constraint, Submod-Max-Cardinality yields a (4 + α)approximation of the optimal solution [15]. The detailed algorithm can be found in Section IV
[13]
V. N UMERICAL R ESULTS
In this section we provide some numerical results demonstrating the performance of the
proposed training/copilot group scheduling scheme. Our training procedure is compared with
a massive MIMO system operating according to the classical TDD protocol, considered as a
reference model. We consider an hexagonal cell network with C = 7 cells. Each cell has a
1.5 Km radius. The mobile users are located uniformly at random in each cell and we assume
that no user is closer than r0 = 10 m to its serving BS. We consider heterogeneous movement
speeds and directions. User speeds are randomly generated in the interval [20 Km/h , 80 Km/h].
User movement direction, which is defined by the angle between the movement direction of
the mobile device and the direction of the incident wave, is also randomly generated over the
interval [0 , 2π]. We consider a path-loss exponent σ = 3.5 and a relatively short coherence
block of Ts = 200 symbols is considered. We suppose a 200 M hz bandwidth [18]. The spectral
efficiency is measured in bits/s/hz. The CSI delays for each copilot group dg , g = 1, ..., Ng are
selected randomly over the interval [0, dmax ].
In Figure 2, we illustrate a comparison of the CDFs of the achievable spectral efficiency
between the reference model and the proposed outdated CSI based training scheme for different
numbers of antennas at the BS. Figure 2 shows that, for 50 receive antennas at the BS, the
proposed training scheme achieves 5%-outage rate around 71 bit/s/hz. This performance represents a gain of 5 bit/s/hz or, equivalently, 1000 Mbit/s for a system bandwidth of 200 M hz
[18], compared with the reference model that achieves 5% outage rate around 66

bit/s/hz.

For 120 antennas at the BS, the gain in the 5%-outage rate when, outdated CSI is used,
attains 1400 Mbit/s. In fact, while the proposed training scheme achieves 5%-outage rate around
87 bit/s/hz, the reference model achieves the same outage rate around 80 bit/s/hz. This gain
comes from the spared resources during uplink training. In fact, using outdated CSI enables to
increase the number of channel uses that are allocated to data transmission.

12

Figure 2: Comparison of the CDFs of spectral efficiency (Ng = 30, τ = 15, dmax = 2 )

Figure 3: Spectral efficiency for varying values of M (Ng = 30, τ = 15, dmax = 2)

In Figure 3, we illustrate the relation between the achievable spectral efficiency and the
number of BS antennas. We compare the lower bound in Theorem 1 with a simulated curve.

13

Figure 3 shows that the derived is tight for a wide range of M . This figure also shows that,
when compared with the reference model, the proposed training scheme increases the achievable
spectral efficiency by 6.8 bit/s/hz for M = 50. This is equivalent to a rate gain of 1360 Mbit/s for
a system bandwidth of 200 M hz [18]. This gain attains 9.9 bit/s/hz or, equivalently, 1980 Mbit/s
for M = 150.

Figure 4: Weighted Spectral efficiency for varying value of τ (Ng = 20, M = 100, dmax = 3 )

In Figure 4, we illustrate the performance of the proposed training/copilot group scheduling
scheme. The weights wgl , ∀g = 1, ..., Ng , l = 1, ..., C are randomly generated in the interval
[0, 1]. Figure 4 shows that the optimal weighted sum rate saturate at τ = 8. This means
that the optimal number of copilot groups scheduled for uplink training is 8 in this setting.
Consequently, scheduling more users for training will not result in a weighted sum rate gain.
This actually concords with the suggestion that not all users who are active for data transmission
are required to participate in uplink training. Consequently, an adaptive training scheme that takes
into consideration CSI aging is justified.

14

VI. C ONCLUSION
In this paper, we studied the achievable sum rate in a massive MIMO system with a matched
filter receiver taking into account the channel aging effect. We proposed a training scheme
in which, we deliberately utilize outdated CSI estimates in order to optimize uplink training.
The derived results show that, although outdated CSI degrades the achievable channel gain,
we are able to achieve better spectral efficiency when training frequency is adapted to the
channel’s autocorrelation. We show that, grouping the users in autocorrelation based clusters
and optimizing their scheduling accordingly, provides a substantial increase in the achievable
weighted sum rate. Future works will include investigating more sophisticated algorithms that
enable to further leverage heterogeneous coherence times.
VII. A PPENDIX
A. Proof of Theorem 1
The network serves Ng copilot groups, τ of which are scheduled for uplink training. At
the reception, each BS uses a matched filter receiver that is based on the latest available CSI
estimates. BS l detects the signal of user g in cell l by applying the following filter
[l]

ugl (t) =

ĝgl (t − dg )
[l]

kĝgl (t − dg )k

, t ≥ dg ,

(17)

[l]

where ĝgl (t − dg ) denotes the latest available CSI estimate for user g in cell l. Consequently,
the detected signal of user g in cell l is given by the following
[l]
Yu (t)
†
ugl (t) √
Pu

=

Ng C
X
X

u†gl (t)(

k=1

=

Wu (t)
[l]
gkc (t)Skc + √ )
Pu
c=1

[l]
u†gl (t)(ggl (t)Sgl

+

C
X

[l]
ggc
(t)Sgc

+

c6=l

= u†gl (t)((ρgl )dg ĝgl (t − dg )Sgl +
[l]

[l]

(18)

Ng C
X
X
k6=g

C
X
c6=l

Wu (t)
[l]
gkc (t)Skc + √ )
Pu
c=1

dg [l]
(ρ[l]
gc ) ĝgc (t − dg )Sgc +

C
X

dg [l]
(ρ[l]
gc ) g̃gc (t − dg )Sgc

c=1

q
Ng C
g −1
C dX
X
X
X [l]
Wu (t)
[l] [l]
[l] j
+
(ρgc ) βgc εgc (t − j)Sgc +
gkc (t)Skc + √ )
Pu
c=1 j=0
k6=g c=1
= u†il (t)(I1 (t) + I2 (t) + I3 (t)),

15

with
[l]

[l]

I1 (t) = (ρgl )dg ĝgl (t − dg )Sgl ,
I2 (t) =

C
X

dg [l]
(ρ[l]
gc ) ĝgc (t − dg )Sgc ,

c6=l

I3 (t) =

C
X

dg [l]
(ρ[l]
gc ) g̃gc (t

c=1

q
Ng C
g −1
C dX
X
X
X [l]
Wu (t)
[l] [l]
[l] j
− dg )Sgc +
(ρgc ) βgc εgc (t − j)Sgc +
gkc (t)Skc + √
Pu
c=1 j=0
k6=g c=1

The third equality in Equation (18) follows from the fact that
[l] [l]

[l]

[l]

[l]

[l]
gic (t)

q
[l] [l]
[l]
= βic hic (t), hic (t) =

[l]

ρic hic (t − 1) + εic (t) for all t and gic (t) = ĝic (t) + g̃ic (t) for all t.
We note that I1 (·) refers to the useful signal, I2 (·) represents the impact of pilot contamination
and I3 (·) regroups the impact of the white noise, channel estimation error, non correlated
interference due to users with different pilot sequences and the impact of channel aging. The
instant spectral efficiency attained by user g in cell l is:
!


|u†gl (t)I1 (t)|2
τ
.
Rg,l = 1 −
log 1 + †
Ts
|ugl (t)I2 (t)|2 + |u†gl (t)I3 (t)|2

(19)

We now define Rg,l to be the average achievable sum rate of user g in cell l, namely,
!!


|u†gl (t)I1 (t)|2
τ
Rg,l = E
1−
log 1 + †
(20)
Ts
|ugl (t)I2 (t)|2 + |u†gl (t)I3 (t)|2
!
!!


|u†gl (t)I1 (t)|2
τ
[l]
ĝgl (t − dg )
,
log 1 + †
=E E
1−
Ts
|ugl (t)I2 (t)|2 + |u†gl (t)I3 (t)|2
0

the last equality follows from the law of total expectation. Let us define Rg,l such that
!
!


|u†gl (t)I1 (t)|2
τ
0
[l]
1−
ĝgl (t − dg ) ,
Rg,l = E
log 1 + †
†
2
2
Ts
|ugl (t)I2 (t)| + |ugl (t)I3 (t)|

(21)

therefore,
0

Rg,l = E(Rg,l ).

(22)

1
Based on the convexity of log(1 + x+a
), and Jensen’s inequality the following inequality can be

obtained
0
Rg,l

!


[l]
[l]
|u†gl (t)(ρgl )dg ĝgl (t − dg )|2
τ
≥ 1−
log 1 +
,
[l]
[l]
Ts
E(|u†gl (t)I2 (t)|2 |ĝgl (t − dg )) + E(|u†gl (t)I3 (t)|2 |ĝgl (t − dg ))
(23)

16

since
E(|u†gl (t)I1 (t)|2 |ĝgl (t − dg )) = |u†gl (t)(ρgl )dg ĝgl (t − dg )|2 ,
[l]

[l]

[l]

(24)

by the property E(f (Z)|Z) = f (Z) for a random variable Z. We now aim at computing
E(|u†gl (t)Ij (t)|2 |ĝgl (t − dg )) for j = 2, 3. In order to do so, we are first going to obtain an
[l]

alternative expression for I2 (t), that is,
I2 (t) =

C
X

[l]
ĝgc
(t

− dg )Sgc =

[l]
ĝgl (t

− dg )

c6=l

[l]

[l]

βgl

[l]

c6=l

[l]

βgc

[l]

since ĝgc (t − dg ) = ĝgl (t − dg )

C
[l]
X
βgc

βgl

Sgc ,

(25)

[l]

. Therefore, I2 (t) and ĝgl (t − dg ) are correlated. Consequently,

we obtain
C
[l]2
h
i
2X
[l]
†
†
2 [l]
[l] 2dg βgc
.
E |ugl (t)I2 (t)| |ĝgl (t − dg ) = ugl (t)ĝgl (t − dg )
(ρgc )
[l]2
β
c6=l
gl

(26)

We will now compute E(|u†gl (t)I3 (t)|2 |ĝgl (t − dg )). First note that, I3 (t) is independent of ĝgl (t −
[l]

[l]

dg ) and since u†gl (t) has unit norm, we have that E(|u†gl (t)I3 (t)|2 |ĝgl (t − dg )) = E(|I3 (t)|2 ),
[l]

therefore we obtain
 X
g −1 q
C
C dX
X
[l]
j [l]
[l] dg [l]
βgc (ρ[l]
E(|I3 (t)| ) = E |
(ρgc ) g̃gc (t − dg )Sgc +
gc ) εgc (t − j)Sgc
2

c=1

+

Ng C
X
X

(27)

c=1 j=0

[l]
gkc (t)Skc

k6=g c=1

Wu (t)
+ √ |2
Pu



X
g −1 q
C
C dX
X
[l]
j [l]
2
[l] dg [l]
2
=E
|(ρgc ) g̃gc (t − dg )| +
| βgc (ρ[l]
gc ) εgc (t − j)|
c=1

+

Ng C
X
X

c=1 j=0
[l]
|gkc (t)|2

k6=g c=1


Wu (t) 2
+| √ | ,
Pu

where the last equality follows from noting the following four properties; (i) Skc · Sic0 = 0 for all
k 6= i and all c, c0 ∈ {0, . . . , C}, (ii) E(ZWu (t)) = E(Z)E(Wu (t)) = 0 for all random variables
Z that are independent of Wu (t) (zero mean complex Gaussian noise), (iii) similar to the previous
[l]

[l]

[l]

property, E(Zεic (t)) = E(Z)E(εic (t)) = 0 for all Z independent of εic (t) (zero mean complex
[l]

[l]

white Gaussian noise) and finally (iv) gkc and g̃k0 c0 are independent for all (k, c) 6= (k 0 , c0 ). We
now compute the four terms in Equation (27). The last term, i.e.,
p
1
E(|Wu (t)/ Pu |2 ) =
.
Pu

(28)

17

Next we compute the second term in Equation (27), namely,

C
C
X
X
dg [l]
2
2dg  [l]
βgc −
E(
|(ρ[l]
(ρ[l]
gc ) g̃gc (t − dg )| ) =
gc )
c=1

c=1

[l]
(βgc )2
P
1
+ C
b=1
Pp


[l]

βgb

.

(29)

[l]

(βgc )2
P
[l]
+ C
b=1 βgb

[l]

the latter is satisfied due to the fact that the variance of g̃gc (t−dg ) being given by βgc −

[l]

1
Pp

for all g and c.
We now compute the third term in Equation (27), that is,
!
g −1
g −1 q
C dX
C dX
X
X
[l] [l] j [l]
2
[l] [l] 2j
2
βgc
(ρgc ) (1 − (ρ[l]
E
| βgc (ρgc ) εgc (t − j)| =
gc ) )

(30)

c=1 j=0

c=1 j=0

=

C
X

[l] 1
βgc

1−

c=1

=

C
X

[l]

− (ρgc )2dg
[l]
(ρgc )2

2
(1 − (ρ[l]
gc ) )

[l]
2dg
βgc
(1 − (ρ[l]
),
gc )

c=1

[l]

for the second equality we have used the expression of finite geometric sums since (ρgc )2 < 1
for all g and c. We are left with the first term in Equation (27), that is,
Ng C
X
X

[l]
E(|gkc (t)|2 )

=

Ng C
X
X

Ng C
q
X
X [l]
[l] [l]
2
E(| βkc hkc (t)| ) =
βkc ,

k6=g c=1

k6=g c=1

(31)

k6=g c=1

Combining all four terms, that is, Equations (28), (29), (30) and (31), we obtain


Ng C
C
C
[l] 2
h
i X
X
X
X
(βgc )
[l]
[l]
[l] 2dg
2dg  [l]

βgc
(1 − (ρgc
) )+
(ρ[l]
βgc −
E |u†gl (t)I3 (t)|2 =
βkn +
PC
gc )
[l]
1
+
β
c=1
c=1
k6=g c=1
b=1 gb
Pp
(32)
=

Ng C
X
X
k6=g c=1

[l]
βkc

+

C
X
c=1

[l]2

[l]
(βgc

−

2dg
ρ[l]
gc
1
Pp

+

βgc
PC

[l]
b=1 βgb

)+

1
.
Pu

Substituting the results in Equations (24), (26) and (32) in Equation (21), we obtain
!


[l]
[l]
(ρgl )2dg |u†gl (t)ĝgl (t − dg )|2
τ
0
Rg,l ≥ 1 −
log 1 +
,
Ts
F

(33)

18

with
F

[l]
=|u†gl (t)ĝgl (t

2

− dg )|

C
X

[l]2

2dg
(ρ[l]
gc )

c6=l

βgc

[l]2
βgl

+

Ng C
X
X
k6=g c=1

[l]
βkc

+

C
X
c=1

[l]2

[l]
(βgc

−

[l]2dg
ρgc
1
Pp

+

βgc
PC

[l]
b=1 βgb

)+

1
.
Pu
(34)

From Equation (20) and Equation (33) we obtain
!!


[l]
[l]
(ρil )2dg |u†il (t)ĝil (t − dg )|2
τ
0
log 1 +
Rg,l = E(Rg,l ) = E
1−
Ts
F
!!


[l]
(ρgl )2dg
τ
=E
1−
log 1 +
,
Ts
G

(35)

where
G=

F
|u†gl (t)ĝgl (t − dg )|2
[l]

.

In order to compute the final expression of the bound on the average rate, it now suffices to
compute the explicit expression of the right hand side (RHS) of Equation (35). In order to do
so, we apply Jensen’s inequality to the RHS in Eq. (35), that is,
!


[l]
(ρg,l )2dg
τ
Rg,l ≥ 1 −
log 1 +
,
Ts
E(G)

(36)

with
!
C
[l]2
X
β
1
gc
2dg
(ρ[l]
E(G) =
+E
gc )
[l]
[l]2
|u†gl (t)ĝgl (t − dg )|2
β
c6=l
gl


Ng C
C
[l]2
X
X [l] X
1
βgc
[l]
[l]2dg
)+
·
βkc +
(βgc
.
− ρgc
PC
[l]
1
Pu
β
+
c=1
gb
k6=g c=1
b=1
Pp
Note that

[l]
u†gl (t)ĝgl (t

sequently, the mean value of
1

to
(M −1)×

[l]2
β
gl
P
[l]
1 + C
β
Pp
b=1 gb

[l]2

2

− dg )

(37)

has a Gamma distribution with parameters (M,
1
2
[l]
u†gl (t)ĝgl (t−dg )

1
+
Pp

βgl
PC

b=1

[l]

βgb

). Con-

(that has an inverse Gamma distribution) is equal

. Combining this together with the results in Equations (37) and (35) we

obtain the desired lower bound on the average achievable spectral efficiency of user g, l, that is,
!


[l]
[l]
(M − 1)(βgl )2 (ρgl )2dg
τ
Rg,l ≥ 1 −
log 1 +
,
(38)
Ts
(M − 1)Iglp + Igln

19

where Iglp and Igln are given by
Iglp

=

C
X

2dg
[l] 2
(ρ[l]
(βgc
) , and
gc )

c6=l

Igln = (

1
+
Pp

C
X


[l]

βgb ) · 

Ng C
X
X

[l]

βkc +

k6=g c=1

b=1

C
X

2dg

[l]
(βgc
− ρ[l]
gc

c=1

[l]2
βgc
PC
1
+
b=1
Pp


[l]
βgb

)+

1
.
Pu

Summing the achievable spectral efficiency of all grouped users concludes the proof.
B. Proof of Theorem 2
In order to prove theorem 2, we consider the assymptotic regime where the number of BS
antennas M grows very large. In this case the lower bound on achievable spectral efficiency of
each user g, l converges to the following limit:


1−

τ
Ts


log 1 +

[l]2 [l]2dg
(M − 1)βgl ρgl
(M − 1) × Iglp + Igln

!





[l]2 [l]2dg
βgl ρgl
τ
 . (39)
log 1 + P
−→ 1 −
[l]2dg [l]2
C
Ts
β
ρ
b6=l

gb

gb

The proposed framework is compared with a reference massive MIMO system where, all scheduled users are required to perform uplink training. In the asymptotic regime, the lower bound
on the achievable spectral efficiency of each user g, l in the reference system converges to the
following limit:


Ng
1−
Ts

[l]2

!

βgl



log 1 + PC

[l]2

.

(40)

b6=l βgb

The aim here, is to improve the achievable spectral efficiency of each scheduled users. Consequently, the spectral efficiency of each user in the two considered systems should verify,
∀ g = 1...Ng , l = 1...C:


!




[l]2 [l]2dg
[l]2
β
ρ
βgl
τ
N
g
gl
gl
≥ 1−
1−
log 1 + P
log 1 + PC [l]2 .
[l]2dg [l]2
C
Ts
Ts
β
ρ
β
b6=l

gb

b6=l

gb

(41)

gb

(27) is equivalent to the following condition:
[l]2 [l]2dg
βgl ρgl

PC

[l]2dg

b6=l ρgb

[l]2

≥

1+

βgb

[l]2
βgl
PC [l]2
b6=l βgb

[min]2

We consider the extreme case where ρ2gl = ρ̄g
[max]
ρ̄g

g
! TTss−N
−τ

− 1.

(42)
[max]2

and ∀b 6= l, ρ2gb = ρ̄g

[min]

. Here ρ̄g

and

denote respectively the minimum and maximum channel autocorrelation coefficients in

20

group g. This means that we assume the worst case scenario for each user where, its coherence
time is always lower than its fellow copilot users. Finally, by considering

[∞]
SIN Rg,l

[l]2

=

βgl
PC

b6=l

[l]2

βgb

,

we obtain (11) which finishes the proof.
C. Proof of Theorem 3
In order to prove theorem 3, we start by demonstrating that the objective function f in (14),
is submodular. We consider two copilot groups allocations, X and Z such that X ⊆ Z. We

define the following ground set G = x[g] , g = 1...Ng , where x[g] , g = 1...Ng is an abstract
element denoting the scheduling of copilot group g for uplink training. We need to prove that
the marginal value of scheduling a new copilot group x[n] in X and Z verifies:




f X ∪ x[n] − f (X) ≥ f Z ∪ x[n] − f (Z) .

(43)

We start by computing the marginal values of adding an element x[n] to X and Z which are,
respectively, given by:
f X ∪ x[n]





− f (X) =

!
 X
C X
[l]2
|X| |X| + 1
(M − 1)βil
−
{
wil log 1 +
Ts
Ts
(M − 1) × Iilp + Iiln
l=1 [i]
x ∈X

(44)
[l]2 [l]2di

[l]2

!

(M − 1)βnl
(M − 1)βil ρil
|X| + 1
)wnl log 1 +
+
wil log 1 +
} + (1 −
p
p
n
n
(M − 1) × Iil + Iil
Ts
(M − 1) × Inl
+ Inl
x[i] ∈G\[X∪x[n] ]
!
[l]2 [l]2dn
|X|
(M − 1)βnl ρnl
− (1 −
)wnl log 1 +
,
p
n
Ts
(M − 1) × Inl
+ Inl

!

X

and,
f Z ∪ x[n]





− f (Z) =

|Z| |Z| + 1
−
Ts
Ts

 X
C X
{
wil log 1 +
l=1 x[i] ∈Z

[l]2

(M − 1)βil
(M − 1) × Iilp + Iiln

!

(45)
[l]2 [l]2di

!

[l]2

|Z| + 1
(M − 1)βnl
(M − 1)βil ρil
+
wil log 1 +
} + (1 −
)wnl log 1 +
p
p
n
n
+ Inl
(M − 1) × Iil + Iil
Ts
(M − 1) × Inl
x[i] ∈G\[Z∪x[n] ]
!
[l]2 [l]2dn
|Z|
(M − 1)βnl ρnl
− (1 −
)wnl log 1 +
.
p
n
Ts
(M − 1) × Inl
+ Inl
X

!

21

The difference between the two marginal values is given by:




(46)
f X ∪ x[n] − f (X) − (f Z ∪ x[n] − f (Z)) =
!
!
[l]2
[l]2 [l]2dn
|Z| − |X|
(M − 1)βnl
|X| − |Z|
(M − 1)βnl ρnl
[
]wnl log 1 +
+[
]wnl log 1 +
p
p
n
n
Ts
(M − 1) × Inl
+ Inl
Ts
(M − 1) × Inl
+ Inl
!
!!
C
[l]2 [l]2di
[l]2
X
X wil
(M − 1)βil ρil
(M − 1)βil
− log 1 +
.
+
[ ] log 1 +
p
p
n
n
+
I
+
I
T
(M
−
1)
×
I
(M
−
1)
×
I
s
il
il
il
il
l=1 [i]
x ∈Z\X

The difference in marginal values is positive. Consequently, the objective function f is submodular. In addition, it is clear that (14a) is a cardinality constraint. Problem (14) is then
equivalent to maximizing a submodular function with a cardinality constraint.
R EFERENCES
[1] T. L. Marzetta, Noncooperative cellular wireless with unlimited numbers of base station antennas, IEEE Trans. Wireless.
Commun., vol. 9, no. 11, pp. 3590-3600, Nov. 2010.
[2] Alexei Ashikhmin, Thomas L. Marzetta and Liangbin LiInterference Reduction in Multi-Cell Massive MIMO Systems I:
Large-Scale Fading Precoding and Decoding, November 2014.
[3] L. Thiele, M. Olbrich, M. Kurras, and B. Matthiesen, Channel aging effects in CoMP transmission: Gains from linear
channel prediction, in Proc. Asilomar Conf. Signals Syst. Comput., Nov. 2011, pp. 1924-1928.
[4] K. Truong and R. Heath, Effects of channel aging in massive MIMO systems, J. Commun. Netw., vol. 15, no. 4, pp.
338-351, Sep. 2013.
[5] A. K. Papazafeiropoulos and T. Ratnarajah,Linear precoding for downlink massive MIMO with delayed CSIT and channel
prediction, in Proc. IEEE Wireless Commun. Netw. Conf. (WCNC), Apr. 2014, pp. 809-914.
[6] A. K. Papazafeiropoulos and T. Ratnarajah, Uplink performance of massive MIMO subject to delayed CSIT and anticipated
channel prediction, in Proc. IEEE Int. Conf. Acoust. Speech Signal Process. (ICASSP), May 2014, pp. 3162-3165.
[7] C. Kong, C. Zhong, A. K. Papazafeiropoulos, M. Matthaiou, and Z. Zhang, Sum-rate and power scaling of massive MIMO
systems with channel aging, IEEE Trans. Commun., vol. 63, no. 12, pp. 4879-4893, Dec. 2015.
[8] Thang X. Vu, Trinh Anh Vu, Symeon Chatzinotas and Bjorn Ottersten, Spectral-Efficient Model for Multiuser Massive
MIMO: Exploiting User Velocity, IEEE ICC’17, Mai 2017
[9] J. Cai, W. Song, and Z. Li, Doppler spread estimation for mobile OFDM systems in rayleigh fading channel, IEEE Trans.
Consumer Electron., vol. 49, pp. 973.977, Nov. 2003.
[10] W. C. Jakes, Ed., Microwave Mobile Communications. New York: Wiley, 1974.
[11] S. Sesia, I. Toufik, and M. Baker, Eds., LTE: The UMTS Long Term Evolution. John Wiley and Sons, 2009.
[12] T. Y. Young and T. W. Calvert, Classification, estimation and pattern recognition, American Elsevier Publishing Company,
1974.
[13] S. Hajri and M. Assaad, M. Larrañaga, Enhancing massive MIMO: A new approach for Uplink training based on
heterogeneous coherence times, http://www.l2s.centralesupelec.fr/perso/salah.hajri/publications
[14] L. Tassiulas and A. Ephremides, Stability properties of constrained queueing systems and scheduling policies for maximum
throughput in multihop radio networks, IEEE Trans. Automat. Contr., vol. 4, pp. 1936-1948, December 1992.

22

[15] A. Gupta, A. Roth, G. Schoenebeck and K. Talwar. Constra ined non-monotone submodular maximiza- tion: Offline and
secretary algorithms. Proc. of WINE , 2010. A more detailed version available at http://arxiv.org/abs/1003.1517v2.
[16] G. L. Nemhauser, L. A. Wolsey, and M. L. Fisher, An analysis of approximations for maximizing submodular set functions
i, Mathematical Programming, 14:265-294, 1978.
[17] U. Feige, V. Mirrokni, and J. Vondrak. Maximizing non-monotone submodular functions. In Proceedings of 48th Annual
IEEE Symposium on Foundations of Computer Science (FOCS) , 2007.
[18] Geoff Varrall, 5G Spectrum and Standards, Artech House, 31 mai 2016

