1

Joint CSI Estimation, Beamforming and Scheduling
Design for Wideband Massive MIMO System

arXiv:1801.05349v1 [cs.IT] 16 Jan 2018

Tadilo Endeshaw Bogale+ , Long Bao Le+ , and Xianbin Wang++
Institute National de la Recherche Scientifique (INRS)
Université du Québec, Montréal, Canada+
University of Western Ontario, London, Canada++
Email: {tadilo.bogale, long.le}@emt.inrs.ca, xianbin.wang@uwo.ca
Abstract— This paper proposes a novel approach for designing
channel estimation, beamforming and scheduling jointly for wideband massive multiple input multiple output (MIMO) systems.
With the proposed approach, we first quantify the maximum
number of user equipments (UEs) that can send pilots which
may or may not be orthogonal. Specifically, when the channel
has a maximum of L multipath taps, and we allocate M̃ subcarriers for the channel state information (CSI) estimation, a
maximum of M̃ UEs CSI can be estimated (L times compared
to the conventional CSI estimation approach) in a massive
MIMO regime. Then, we propose to schedule a subset of these
UEs using greedy based scheduling to transmit their data on
each sub-carrier with the proposed joint beamforming and
scheduling design. We employ the well known maximum ratio
combiner (MRC) beamforming approach in the uplink channel
data transmission. All the analytical expressions are validated
via numerical results, and the superiority of the proposed design
over the conventional orthogonal frequency division multiplexing
(OFDM) transmission approach is demonstrated using extensive
numerical simulations in the long term evolution (LTE) channel
environment. The proposed channel estimation and beamforming
design is linear and simple to implement.
Index Terms— Beamforming, Channel estimation, Massive
MIMO, Multipath taps, OFDM, Pilot design, Scheduling

I. I NTRODUCTION
Massive multiple input multiple output (MIMO) technology is one of the promising means for achieving the energy and spectrum efficiency requirements of the future 5G
networks [1]–[3]. The beamforming gain envisaged by the
massive MIMO system depends on the availability of channel
state information (CSI). In general, the channel between the
transmitter and receiver is estimated from orthogonal pilot
sequences where the number of such sequences is limited
by the channel coherence time [4]–[7]. In a multiuser setup,
massive MIMO can be deployed both at the base stations (BSs)
and mobile stations (UEs). However, due to space, energy
consumption and cost constraints, deploying massive number
of antennas at the UEs is usually infeasible from practical
viewpoint particularly at microwave frequency bands. For this
reason, it is economical to deploy antenna arrays at the BSs
and single antennas at each UE for microwave frequency band
applications. This motivates us to consider a massive MIMO
system where each UE has single antenna as in [4].
In a massive MIMO system, channel estimation can be
performed either in frequency division duplex (FDD) or time
division duplex (TDD) approach. In both approaches, the

pilot sequences need to be orthogonal to maintain sufficient
channel estimation quality. In a traditional pilot transmission,
the number of pilots may need to be the same as the number
of transmitted antennas to maintain orthogonality. Thus, for
a massive MIMO system, one can estimate the maximum
number of channel coefficients using time division duplex
(TDD) approach by sending pilots from UEs. By doing so,
each BS can estimate and utilize the channel coefficients both
for the downlink and uplink transmissions. A number of approaches have been proposed to address channel estimation for
massive MIMO systems including a sub-space method which
exploits the sub-space information of the channel covariance
matrices [5], [7], [8], pilot optimization/scheduling/shifting
method which optimizes/schedules/shifts the pilot sequences
of each cell (UE) while utilizing appropriate signal dimensions
such as time and frequency [6], [9]–[11].
The problem of channel estimation for wideband massive
MIMO system is studied in the context of multicell system in
[12]. When the channel has L multipath taps, each cell has
M̃ sampling periods for pilot transmission and each BS serves
the number of UEs equal to K̃ = M̃
L which is the number
of UEs obtained with orthogonal pilots (see (4) of [4] and
Section II of [10]), it is shown in [12] that at least L cells can
reliably estimate the channels of their user equipments (UEs)
and perform beamforming without experiencing pilot contamination. This result has the same significance as increasing
the number of UEs by L times for the same time/frequency
resource. Although such an increment is useful, the design
approach of [12] has meaningful improvement compared to
the traditional orthogonal pilot in the massive MIMO regime.
The main reason the approach of [12] performs badly in a
small number of antennas regime is that in such a regime,
each UE will experience strong interference signal from the
other UEs which are generally large compared to the one
used in orthogonal pilots. This motivates us to schedule the
number of UEs adaptively while maximizing the total sum
spectrum efficiency (SE). This leads to the joint design of CSI
estimation, scheduling and beamforming for wideband massive
MIMO system which is the focus of the current paper. In
[12], pilot transmission is performed in time domain (i.e., nonorthogonal frequency division multiplexing (OFDM)). However, in the current paper, we apply pilot transmission in
frequency domain. As will be clear in the sequel, such a pilot

transmission approach naturally leads the scheduling design to
select only a subset of the UEs for transmission.
In particular, the current paper proposes a novel approach
for designing channel estimation, beamforming and scheduling jointly for wideband massive MIMO systems. With the
proposed approach, we first quantify the maximum number
of UEs that can send pilots which may or may not be
orthogonal. Specifically, when the channel has a maximum
of L multipath taps, and allocate M̃ sub-carriers for the CSI
estimation, a maximum of M̃ UEs can use these sub-carriers
for their CSI estimation (i.e., L times more UEs compared
to the traditional minimum mean square error (MMSE) CSI
estimation approach) in a massive MIMO regime. Then, we
select only a subset of these UEs to transmit their data for
each sub-carrier with the proposed joint beamforming and
scheduling design. We employ the well known maximum
ratio combiner (MRC) beamforming approach in the uplink
channel data transmission. All the analytical expressions are
validated via numerical results, and the superiority of the
proposed design over the conventional transmission approach
is demonstrated using extensive numerical simulations for
LTE channel models. The proposed channel estimation and
beamforming design is linear and simple to implement.
This paper is organized as follows. Section II discusses the
system and channel models. In Section III, a brief summary
of conventional OFDM based pilot and data transmission is
discussed. In Sections IV and V, the proposed joint channel
estimation, beamforming and scheduling design is provided. In
Section VI, extensive simulation results are presented. Finally,
Section VII draws conclusions.
II. S YSTEM AND C HANNEL M ODEL
We consider a multiuser system with a transmission scheme
having a sampling period Ts and a channel having a maximum
delay spread Td . It is assumed that the BS and each UE are
equipped with N antennas and 1 antenna, respectively. For
this setting, the number of multipath channel impulse response
(CIR) taps L is approximated as L = TTds [13]. The multipath
coefficients between the kth UE and nth BS antenna is denoted
as [12]
h̄kn = [h̄kn1 , h̄kn2 , · · · , h̄knL ]T .

(1)

The analytical model for H̄k = [h̄k1 , h̄k2 , · · · , h̄kN ] may
vary from one standard to another. In this paper, we assume
that h̄kn are correlated temporarily over multipath taps and
spatially over the BS antenna arrays. We further assume that
the temporal correlation matrices corresponding to all antennas
are the same, and the spatial correlation matrix of all multipath
components are the same. Under such assumptions, H̄k is
modeled as [14]
√
√
H̄k = Pk H̃k Rk
(2)
where each entry of H̃k is an i.i.d ZMCSCG random variable
with unit variance, Pk ∈ C L×L and Rk ∈ C N ×N are
positive semidefinite temporal and spatial correlation matrices,
respectively. Furthermore, the correlation between an antenna
and itself (i.e., diagonal elements of Rk ) is assumed to be 1.

Fig. 1: Existing pilot allocation.
In most wireless standards (e.g., Long term evolution
(LTE)), OFDM based transmission is commonly adopted. For
such a transmission, the channel coefficient of each sub-carrier
has practical importance which can be obtained by linearly
combining the multipath tap coefficients (1). To this end, the
channel between the kth UE and nth antenna BS in sub-carrier
s can be expressed as [13], [15]
hkns = f̃sH h̄kn
2π

2π

(3)
2π

where f̃sH = [1, e−j M̃ s , e−j M̃ 2s , · · · , e−j M̃ (L−1)s ] with M̃
being the fast Fourier transform (FFT) size of the OFDM
during pilot transmission. As can be seen from this expression,
for fixed L, h̄kn and sampling period Ts , the CSI coefficient
at sub-carrier s is different for different M̃ .
III. C ONVENTIONAL PILOT AND DATA TRANSMISSION
APPROACH

Here we discuss the conventional pilot and data transmission
approaches. Such an approach first utilizes orthogonal pilots
for each sub-carrier during channel estimation as shown in
Fig. 1 for the settings L = 4 and M̃ = 16. As can be seen
from this figure, each UE may need to spend at least L = 4
pilots for this setup to maintain orthogonality (For example,
the first UE utilizes sub-carriers 0, 4, 8 and 12). This shows
that with the conventional approach, the BS can estimate the
CSI information for a maximum of K̃ = M̃
L = 4 UEs.
In fact, such an estimator yields the anticipated theoretical
beamforming performance during the data transmission stage,
and maintains a non-decreasing signal to interference plus
noise ratio (SINR) as N increases and leads to unbounded
SINR when N → ∞. However, if K̃ > M̃
L (i.e., a reusing of
pilots in one or more sub-carriers), the beamforming performance suffers from an SINR floor even if N → ∞. The reason
of this phenomena, the reuse of pilot sequences in the same
time (frequency) resource, is termed as pilot contamination
[4].
We would like to emphasize here that pilot transmission
can also be performed without OFDM. With this pilot transmission, one also gets K̃ = M̃
L = 4 (see the discussion after
(6) of [12], (4) of [4] and Section II of [10]). Thus, with
the conventional pilot transmission, transmitting pilots either
with OFDM or without OFDM approach does not bring any
difference on the maximum number of supported UEs. In
addition, the conventional approach estimates each UE’s subcarrier CSI without taking into account the data transmission
strategy. Furthermore, with the conventional approach, one
can estimate the CSI of each sub-carrier independently for

As we can see from this expression, ĥks depends on the
linear combination vector tks and the pilot matrix Xk which
are designed based on the beamforming strategy during data
transmission as discussed below.
B. Data Transmission
As discussed above, we consider the uplink channel for data
transmission using the OFDM approach. For this reason, the
received signal of each sub-carrier can be obtained independently at each antenna of all BSs. To this end, the nth BS
antenna receives the following signal in sub-carrier s (yins )

Fig. 2: Proposed pilot allocation.
each antenna. In the following, we discuss the proposed joint
channel estimation and data transmission approach.
IV. P ROPOSED J OINT CSI ESTIMATION AND
B EAMFORMING
This section discusses the proposed channel estimation and
beamforming approach. We employ the commonly adopted
maximum ratio combiner (MRC) beamforming1 .

The main idea of the proposed pilot transmission is to let
every UE to use all the sub-carriers to transmit their pilot
symbols as shown in Fig. 2, and then to compute each UE’s
sub-carrier CSI as a linear combination of the received signals
of all sub-carriers, where the linear combination coefficients
and pilot sequences are selected carefully to ensure that each
UE’s sub-carrier SINR increases as N increases (and it grows
unbounded as N → ∞)
If we assume the total number of supported UEs as K, the
received signal by the nth antenna at sub-carrier s becomes
rns =

K
X

fsH h̄kn dks + w̃ns

k=1

where w̃ns is the noise sample at the nth BS antenna and subcarrier s during data transmission and dks is the data symbol
transmitted in sub-carrier s of the kth UE. The overall received
signal at the sth sub-carrier becomes
ys =

K
X

HTi fs∗ dis + w̃s = HTk fs∗ dks +

i=1

A. Proposed Pilot Transmission

K
X

yns =

f̃sH hin xis + wns

R = [r1 , r2 , · · · , rN ] =

K
X

where
ys
=
[y1s , y2s , · · · , yN s ]T ,
fsH
=
2π
2π
s
−j
2s
−j
(L−1)s
−j 2π
[1, e M , e M , · · · , e M
]
and
w̃s
=
[w̃1s , w̃2s , · · · , w̃N s ]T . It is assumed that each entry of
w̃s is an i.i.d ZMCSCG random variable with variance σ̃ 2 .
Now, let us assume that we are interested in estimating dks
using a beamforming vector aks ∈ C N ×1 as
H
dˆks = aH
ks ys = aks (

T

ĥks =R tks .

(4)

1 The extension of the proposed technique for other beamforming techniques
is discussed in [16].

HTi fs∗ dis + w̃s ).

(5)

As can be seen from this expression, there are two main indeT ∗
H
pendent components: aH
ks Hi fs and aks w̃s , and the average
powers of these terms are given as (see Appendix A of [16])
,θkis

z
}|
{
T ∗ 2
H
∗ T
∗ 2
E{|aH
ks Hi fs | } = |tks Xi F̃ Cii fs |
+

z
K
X

,µkis

}|

z

2
2
E{|aH
ks w̃s | } =σ

+ tH
ks

}|



{

X∗u F̃T Cimus F̃∗ XTu tks

m=1 u6=m

i=1

z

,µ̄kis

X
K
K X

z
K
X

{

∗
T
∗ T
tH
ks Xm F̃ Cimms F̃ Xm tks

m6=i

i=1

where xi = [xi1 , xi2 , · · · , xiM̃ ]T , rn = [rn1 , rn2 , · · · , rnM̃ ]T ,
wn = [wn1 , wn2 , · · · , wnM̃ ], Xi = diag(xi ), W =
[w1 , w2 , · · · , wN ], Hi = [hi1 , hi2 , · · · , hiN ] and F̃ =
[f̃1 , f̃2 , · · · , f̃M̃ ].
We express the kth UE sth sub-carrier CSI estimate as

K
X
i6=k

+ tH
ks

Xi F̃H Hi + W

HTi fs∗ dis + w̃s

i6=k

i=1

where xis is the pilot symbol sent by the ith UE in subcarrier s and wns the additive noise at the nth BS antenna
and sth sub-carrier which is assumed to be independent and
identically distributed (i.i.d) zero mean circularly symmetric
complex Gaussian (ZMCSCG) random variable with variance
σ2 .
The received signals of all antennas and sub-carriers can be
expressed as

K
X

,ξks

}|

{

∗ T
∗ T
tH
ks Xi F̃ Cii F̃ Xi tks
,ξ̄ks

X
K X
K

i=1 m6=i

}|



{

X∗i F̃T Cim F̃∗ XTm tks

(6)

where Cki
=
E{H∗k HTi } and Ckims
=
∗ T ∗ T
∗ T
E{Hi Hk fs fs Hk Hm }. With the proposed approach,
one can achieve the following average SINR with dˆks
T ∗ 2
E{|aH
ks Hk fs | }

γ̄ks = PK

i6=k

= PK

i6=k

T ∗ 2
H
2
E{|aH
ks Hi fs | } + E{|aks w̃s | }
θkks + µkks + µ̄kks
.
(θkis + µkis + µ̄kis ) + ξks + ξ¯ks

(7)

V. P ROPOSED TWO PHASE SCHEDULING DESIGN
Once the joint CSI estimation and beamforming procedures are determined, the next phase will be to perform UE
scheduling. In general, a scheduler is designed to optimize
some performance criteria. In the current paper, we consider
the maximization of the total sum SE achieved by all subcarriers with a per UE power constraint since it is an uplink
transmission. This problem can be mathematically formulated
as
max

K,K,tks

K
M X
X

s=1 k=1

log(1 + γ̄ks ), Pks ≤ Pmax , k ∈ Kmax (8)

where Pmax is the maximum transmitted power at each subcarrier of each UE and K is the scheduled set of UEs out
of Kmax sets. This problem is a mixed integer programming
non-convex problem whose global optimal solution is difficult
to obtain. In the following, we discuss the proposed two phase
approach to solve this problem sub-optimally.
For the appropriately selected K and xk , one can exploit the
following behaviors in (6): θkis scales with N 2 , µkis and ξks
scale with N L, and µ̄kis and ξ¯ks decreases as N increases. For
this reason, we propose a two stage scheduling approach. In
the first stage, we employ θkis to carefully select the maximum
number of UEs Kmax and xk while ensuring that γ̄ks with
N and becomes unbounded as N → ∞. This is due to the
fact that N ≫ L in a massive MIMO regime for a fixed L.
In the second phase, we tune 1 ≤ K ≤ Kmax and select the
corresponding UEs using greedy based scheduling by taking
into account θkis , µkis and ξks as the contributions of these
terms scale with LN .
In the first phase, we select Kmax to maintain a nondecreasing γ̄ks with N and becomes unbounded as N → ∞.
In fact one can express γ̄ks
1

γ̄ks = PKN

θkks + β0

1
i6=k N θkis

(9)

+ β1
PK

+ µ̄kks ) and β1 =
where β0 =
i6=k (µkis +
µ̄kis ) + ξks + ξ¯ks ) which are the terms independent of N .
It follows that to maintain a non-decreasing γ̄ks with N and
becomes unbounded as N → ∞, θkks should increase with N 2
while maintaining θkis = 0, i 6= k. One straightforward way
of achieving this is by designing xk and tks while ensuring
1
N (µkks

max

K,xk ,tks

1
N(

∗
|tH
ks F̄ks xk |

∗
s.t tH
ks F̄is xi = 0, ∀i 6= k

(10)

where F̄ks = diag(F̃T Ckk fs∗ ). For arbitrary Xk and Ckk ,
a linear combination vector tks ensuring the equality constraints does not exist if K > M̃ . This can be verified by
considering i.i.d channel Ckk = I. For this reason, we allow
a maximum of Kmax = M̃ UEs to transmit their pilots
at the pilot transmission phase2 . In fact, one can maintain
an increasing average SINR with increasing N if and only
∗
H
∗
if |tH
ks F̄ks xk | > 0, tks F̄is xi = 0, ∀i 6= k is maintained.
This is possible when we select xk which ensures Qs ,
[F̄1s x∗1 , F̄2s x∗2 , · · · , F̄Ks x∗K ] is a full rank matrix. For a
general setup, it is not trivial to select xk to ensure a full
rank Qs for any Ckk and M̃ = M (if it exists). However, we
have observed that Qs , ∀s achieves higher rank when xk , ∀k
are selected from the orthogonal Zadoff-Chu sequence. This
motivates us to select xk , ∀k as the K orthogonal Zadoff-Chu
sequences.
In the second phase, we select 1 ≤ K ≤ Kmax by
taking into account θkis , µkis and ξks as the contributions
of these terms scale with LN . From the first phase, although
K = Kmax UEs can be scheduled for data transmission such
a choice of K may not yield the maximum sum SE (9)3 . Thus,
for a given N , SNR and channel environment, we search the
appropriate K exhaustively that yields the maximum total sum
SE over all sub-carriers. One way of selecting K UEs is by
considering all possible combinations of the UEs. However,
such an approach has an exponential complexity. For this
reason, we incrementally set K and select the corresponding
UEs by applying the well known greedy based scheduling
approach to significantly reduce the design complexity [18].
We would like to mention here that one can also further
improve the total sum SE by greedily selecting the UEs per
sub-carrier from the K UEs. More details in this direction can
be found in [16].
We would like to mention here that with the existing
approach, we have obtained K̃. One can also apply greedy
based approach to schedule only a subset of the UEs while
maximizing the total sum rate. Such an optimization, however,
may need to be performed for each channel realization which
makes this approach computationally expensive compared to
that of the proposed approach which uses channel statistics to
schedule the UEs. This is due to the fact that the proposed
scheduling is performed when the channel covariance matrix
is modified which occurs in a much slower rate compared to
the CSI variation.
VI. S IMULATION R ESULTS
This section provides simulation results. In general, different multipath channel propagation models are suggested for
different standards, where each model characterizes both the
small scale and large scale fadings. In a typical environment,
the large scale fading captures shadowing, power delay profile
(PDP), and distance dependent path loss, spatial correlation
over multiple antennas and temporal correlation over multipath
2 Note that in some scenario, it is possible to allow K > M̃ (for example
when Cii Ckk = 0, ∃i as discussed in [5] and [17]). However, such a special
case is not considered here.
3 This will be demonstrated in the simulation section.

L̃
X
i=1

¯ g((l − Q)T − τ ), l = 1, · · · , L
q̄i h̄
kni
s
i

(11)

where L̃ is the size of PDP coefficients, Q is a shift parameter
selected to maintain h̄knl = 0, l ≤ 0, g(t) is the overall
transceiver pulse shaping filter, τi is the delay coefficient
¯
corresponding to the ith path, and h̄
kni is the Rayleigh fading
coefficient corresponding to the ith path, kth UE and nth
BS antenna. One can observe from (11) that L depends on
several factors including g(t), Ts and τi . In fact, for a given
Td = maxi,j |τi − τj |, different sets of τi , q̄i may lead to a
different statistical properties of h̄knl .
All results are obtained by averaging 10000 channel real2
i| }
. We set K̃ = 4
izations. The SNR is defined as γ = E{|d
σ2
and Rk = I, the temporal correlation matrices of multipath
channels are computed from the PDP, and g(t) is set as a
raised cosine filter (RCF) with roll off factor α = 0.254 .
A. Effect of PDP and bandwidth on rank(Qs )
The performance of the scheduler depends on rank(Qs ),
and L generally increases with B for a fixed delay spread Td
and g(t) which ultimately help increase Kmax . In this subsection, we examine rank(Qs ) for dense and sparse channel
paths. In general, one may experience Ckk be closer to identity
in the dense channel paths which yields rank(Qs ) ≈ Kmax .
On the other hand, for the case of sparse channel paths,
rank(Qs ) will likely be much less than Kmax .
To demonstrate this, we plot the CDF of the SVD of Qs
(i.e., the percentage of SVD(Qs ) > T h) averaged over all subcarriers and UEs for different threshold values T h when the
temporal covariance matrices of all UEs are computed from
the PDP parameters of the LTE and custom designed PDPs
as shown in Fig. 3. One can observe from this figure that the
the percentage of SVD(Qs ) > T h increases as the multipath
density increases. This can be confirmed by comparing the
CDF(SVD(Qs ) > T h) of uniformly spaced delays with
∆τ = 10 and ∆τ = 50ns, and the LTE environments LTEEPA, LTE-ETA, and LTE-ETU. Next we demonstrate the size
of SVD(Qs ) > T h averaged over all sub-carriers and UEs for
different values of B as shown in Fig. 4. As can be seen from
this figure, increasing B helps increase L and the increment
has a significant performance gain when L̃ is large which is
expected.
The figures 3 and 4 assume that the total number of UEs
sending their pilot is set to K = Kmax = M̃ = K̃L. For
this setup, in fact, CDF(SVD(Qs ) > T h) decreases as T h
increases in all environments which may not be desirable for a
4 Note that one can obtain g(t) by employing a square root RCF with
α = 0.25 at the transmitter and receiver.

0.8

CDF(SVD(Q s)>Th)

h̄knl =

1

0.6

0.4
LTE-EPA
LTE-ETA
LTE-ETU
Uniformly spaced delay ( =50ns)
Uniform spaced delay ( =10ns)

0.2

0
-6

-3

0

3

6

9

12

15

Threshold (Th) in dB

Fig. 3: The CDF of the SVD(Qs ) > T h averaged over all
sub-carriers and UEs when Ts = 7.68µs (i.e., B = 5MHz,
M = 256) and different environments. For the uniformly
spaced delay setup, we use q̄i = 0 and [τ1 , τ2 , · · · ] = [0 :
∆τ : 5000].
450
LTE-ETU (Th=6dB)
LTE-ETU (Th=12dB)
Uniformly spaced delay (
Uniformly spaced delay (

400
350

|SVD(Q s)>Th| 0

channel taps [19]. In this simulation, we ignore the effect of
shadowing, Doppler effect, and assume that all UEs experience
the same distance dependent path loss and PDP. Further, we
utilize the PDP information to get the channel coefficients by
assuming that the small scale fading corresponding to each
delay tap is modeled by a Rayleigh fading. With these settings,
h̄knl , l = 1, · · · , L becomes

=10ns, Th=6dB)
=10ns, Th=12dB)

300
250
200
150
100
50
0
0

2.5

5

7.5

10

12.5

15

17.5

20

Bandwidth B (MHz)

Fig. 4: The |SVD(Qs ) > T h|0 averaged over all subcarriers and UEs for different values of B with LTE-ETU and
uniformly spaced delay parameter setups.
high SNR transmission applications. For arbitrary T h, one way
of increasing the CDF of SVD(Qs ) > T h is by decreasing
K as desired. Fig. 5 shows the CDF of SVD(Qs ) > T h for
different values of T h and Ts = 7.68µs (i.e., B = 5MHz,
M = 256) for LTE-ETU and uniformly spaced delay setup
with ∆τ = 10ns. One can observe from this figure and Fig.
3 that decreasing K definitely helps increase the number of
UEs experiencing higher SNR.
B. Comparison of OFDM and proposed design
In this simulation, we compare the SE performance of the
proposed channel estimation, beamforming and two phase
scheduling design with that of the conventional OFDM approach. In this regard, we consider the PDP’s used in LTEETU environment with B = 5MHz, SN R = 0dB and
different values of N and K as shown in Fig. 6. One can
observe from this figure that different choices of K yields

K for data transmission using the greedy based scheduling
for each sub-carrier with the proposed joint beamforming and
scheduling design. All the analytical expressions have been
demonstrated via numerical results, and the superiority of the
proposed design over the conventional transmission approach
is demonstrated using extensive numerical simulations for LTE
channel models. In particular, for larger K, the proposed
approach requires larger N to achieve better total SE compared
to the conventional OFDM approach. Consequently, the benefit
of the proposed design can be exploited by tuning K according
to the channel environment and number of BS antennas.

1.2

CDF(SVD(Q s)>Th)

1

0.8

0.6
LTE-ETU (K=0.25Kmax )

0.4

LTE-ETU (K=0.5Kmax )

0.2

0
-6

Uniformly spaced delay (

=10ns, K=0.25Kmax )

Uniformly spaced delay (

=10ns, K=0.5Kmax )

-3

0

3

6

9

12

15

Threshold (Th) in dB

Fig. 5: The CDF of the SVD(Qs ) > T h averaged over all
sub-carriers and UEs when Ts = 7.68µs (i.e., B = 5MHz,
M = 256) for different values of K. For the uniformly spaced
delay setup, we use q̄i = 0 and [τ1 , τ2 , · · · ] = [0 : ∆τ : 5000].
35
Conventional OFDM
Proposed (K=0.1Kmax )

Total SE (bps/Hz)

30

Proposed (K=0.3Kmax )

25

Proposed (K=0.5Kmax )

20
15
10
5
0
0

10

20

30

40

50

60

70

80

Number of BS antennas (N)

Fig. 6: Comparison of OFDM and proposed design when Ts =
7.68µs (i.e., B = 5MHz, M = 256) for different values of K
in the LTE-ETU environment.
different average spectrum efficiency (SE). Furthermore, when
K increases the proposed approach requires large N to achieve
better total SE compared to the conventional OFDM approach.
This validates that the benefit of the proposed design can be
exploited by tuning K according to the channel environment
and number of BS antennas.
VII. C ONCLUSIONS
This paper proposes a novel approach for designing channel
estimation, beamforming and scheduling jointly for wideband
massive MIMO systems. With the proposed approach, we first
quantify the maximum number of UEs that can send pilots
which may or may not be orthogonal. Specifically, when the
channel has a maximum of L multipath taps, and allocate M̃
sub-carriers for the CSI estimation, a maximum of M̃ UEs’
CSI can be estimated (i.e., L times larger compared to the traditional MMSE based CSI estimation approach) in a massive
MIMO regime. Then, we select only a subset of these UEs

R EFERENCES
[1] F. Boccardi, R. W. Heath, A. Lozano, T. L. Marzetta, and P. Popovski,
“Five disruptive technology directions for 5G,” IEEE Commun. Mag.,
Feb. 2014.
[2] L. B. Le, V. Lau, E. Jorswieck, N.-D. Dao, A. Haghighat, D.-I. Kim, and
T. Le-Ngoc, “Enabling 5G mobile wireless technologies,” EURASIP J.
Wireless Commun. Networking., vol. 218, pp. 1 – 14, 2015.
[3] T. E. Bogale and L. B. Le, “Massive MIMO and millimeter wave for 5G
wireless HetNet: Potentials benefits and challenges,” IEEE Veh. Techno.
Mag., 2016.
[4] 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.
[5] H. Yin, D. Gesbert, M. C. Filippou, and Y. Liu, “A coordinated approach
to channel estimation in large-scale multiple-antenna systems,” IEEE J.
Select. Areas in Commun., vol. 31, no. 2, pp. 264 – 273, Feb. 2013.
[6] T. E. Bogale and L. B. Le, “Pilot optimization and channel estimation
for multiuser massive MIMO systems,” in Proc. IEEE Conference on
Information Sciences and Systems (CISS), Mar. 2014, pp. 1 – 5.
[7] Q. N. Hien and E. G. Larsson, “EVD-based channel estimation in
multicell multiuser MIMO systems with very large antenna arrays,” in
ICASSP, 2012, Kyoto, Japan, 2012, pp. 3249 – 3252.
[8] R. R. Muller, L. Cottatellucci, and M. Vehkapera, “Blind pilot decontamination,” IEEE J. Select. Topics in Signal Process., vol. 8, no. 5, pp.
773 – 786, Oct. 2014.
[9] T. X. Vu, T. A. Vu, and T. S.Q Quek, “Successive pilot contamination elimination in multi-antenna multi-cell networks,” IEEE Wireless
Commun. Letters, pp. 617 – 620, Nov. 2014.
[10] F. Fernandes, A. Ashikhmin, and T. L. Marzetta, “Inter-cell interference
in noncooperative TDD large scale antenna systems,” IEEE J. Select.
Areas in Commun., vol. 31, no. 2, pp. 192 – 201, Feb. 2013.
[11] W. A. W. M. Mahyiddin, P. A. Martin, and P. J. Smith, “Pilot contamination reduction using time-shifted pilots in finite massive MIMO
systems,” in Proc. IEEE Vehicular Technology Conference (VTC Fall),
Sep. 2014, pp. 1 – 5.
[12] T. E. Bogale, L. B. Le, X. Wang, and L. Vandendorpe, “Pilot
contamination in wideband massive MIMO system: Number of cells
vs multipath,” in Proc. IEEE Global Telecommun. Conf., 2015.
[13] D. Tse and P. Viswanath, Fundamentals of Wireless Communication,
Cambridge University Press, 2005.
[14] F. Talebi and T. Pratt, “Channel sounding and parameter estimation for a
wideband correlation-based MIMO model,” IEEE Trans. Veh. Technol.,
pp. 1 – 9, 2015.
[15] H. Zarrinkoub, Understanding LTE with MATLAB: From Mathematical
Modeling to Simulation and Prototyping, John W&S, Ltd, 2014.
[16] T. E. Bogale, L. B. Le, and X. Wang, “Joint CSI estimation, beamforming and scheduling design for wideband massive MIMO system,”
2017,
https://www.dropbox.com/s/511rzcqfu0uwnyq/
JIEEEDP1.pdf?dl=0 (Available).
[17] H. Yin, D. Gesbert, M. C. Filippou, and Y. Liu, “Decontaminating
pilots in massive MIMO systems,” in Proc. IEEE Conference on
Commnications (ICC), Jun. 2013, pp. 3170 – 3175.
[18] T. E. Bogale, X. Wang, and L. B. Le, “Adaptive channel prediction,
beamforming and scheduling design for 5G V2I network,” in VTC-Fall,
Sep. 2017.
[19] 3GPP TS 36.104, “Base station (BS) radio transmission and reception,”
2012, 3rd Generation Partnership Project; Technical Specification Group
Radio Access Network; Evolved Universal Terrestrial Radio Access (EUTRA).

