1

Multipath Multiplexing for Capacity
Enhancement in SIMO Wireless Systems
(Draft with more Results and Details)
arXiv:1707.07067v1 [cs.IT] 21 Jul 2017

Tadilo Endeshaw Bogale, Member, IEEE, Long Bao Le, Senior Member, IEEE,
Xianbin Wang Fellow, IEEE and Luc Vandendorpe Fellow, IEEE

Abstract
This paper proposes a novel and simple orthogonal faster than Nyquist (OFTN) data transmission
and detection approach for a single input multiple output (SIMO) system. It is assumed that the signal
having a bandwidth B is transmitted through a wireless channel with L multipath components. Under
this assumption, the current paper provides a novel and simple OFTN transmission and symbol-bysymbol detection approach that exploits the multiplexing gain obtained by the multipath characteristic
of wideband wireless channels. It is shown that the proposed design can achieve a higher transmission
rate than the existing one (i.e., orthogonal frequency division multiplexing (OFDM)). Furthermore, the
achievable rate gap between the proposed approach and that of the OFDM increases as the number of
receiver antennas increases for fixed L. This implies that the performance gain of the proposed approach
can be very significant for a large-scale multi-antenna wireless system. The superiority of the proposed
approach has been shown theoretically and confirmed via numerical simulations. Specifically, we have
found upper-bound average rates of 15 bps/Hz and 28 bps/Hz with the OFDM and proposed approaches,
respectively, in a Rayleigh fading environment with 32 receive antennas and signal to noise ratio (SNR)
of 15.3 dB. The extension of the proposed approach for different system setups and associated research
problems is also discussed.
Index Terms
Part of this work has been presented in IEEE GLOBECOM 2016 conference. Tadilo Endeshaw Bogale and Long Bao Le are
with the Institute National de la Recherche Scientifique (INRS), Université du Québec, Montréal, Canada and Xianbin Wang is
with the University of Western Ontario, London, Ontario, and Luc Vandendorpe is with the University Catholique de Louvain,
Belgium. Email: {tadilo.bogale, long.le}@emt.inrs.ca, xianbin.wang@uwo.ca and luc.vandendorpe@uclouvain.be

Faster than Nyquist, FTN, Orthogonal faster than Nyquist, OFTN, Massive MIMO, Multipath fading,
Multipath multiplexing, OFDM, Pulse shaping

I. I NTRODUCTION
The study in Wireless World Research Forum (WWRF) has predicted that 7 trillion wireless
devices will be supported by wireless networks, both for human and machine-type communications, in 2017. Furthermore, the mobile data traffic has increased dramatically over the years,
which is mainly driven by the massive demand of data-hungry devices such as smart phones,
tablets and broadband wireless applications such as multimedia, 3D video games, e-Health, Car2X
communications [1], [2]. It is expected that the future 5G wireless network will deliver about
1000 times more capacity than that of the current 4G system. To meet this demand different
enabling techniques have been under development including massive multiple input multiple
output (MIMO), millimeter wave, full duplex communication and faster than Nyquist (FTN)
transmission [3]–[6].
The first work on FTN transmission was proposed by Mazo in 1975 where it is claimed
that it is possible to achieve a transmission rate beyond what can be achieved with the Nyquist
criterion without requiring additional bandwidth [7]. In general, the symbol error probability plays
a fundamental role in understanding a communication strategy which can be directly related to
the minimum Euclidean distance dmin between two distinct symbols. The key idea of [7] is that
for the scenario with Nyquist symbol duration TN y and sinc pulse shaping filter, dmin remains the
same even when we pack a symbol in τ TN y with τ ∈ [0.802, 1]. This result implies that one can
still maintain the same error probability as that of the Nyquist transmission by selecting τ ≈ 0.82

which in turn helps to improve the spectrum efficiency by

1
0.802

≈ 25%. It is also shown that

different pulse shaping filters result in different Mazo limits (i.e., minimum τ ). For instance, the

minimum τ with a root raised cosine filter with excess bandwidth β = 0.3 yields a 42% increase
in bandwidth efficiency (i.e., τmin = 0.703) [8].
These results develop new perspectives about the achievable rate and have recently attracted
researchers (see for example [3], [4], [9]). Having said this, however, there is one major difference
between Nyquist and FTN transmissions. The data symbols transmitted with the Nyquist approach
can be decoded independently on a symbol-by-symbol basis with appropriate matched filtering
which consequently facilitates low complexity receiver. However, as the FTN transmission does

not enable orthogonal symbol-by-symbol detection, it requires more complicated symbol detections [3], [4], [10]. In this regard, low-complexity receiver algorithms are suggested by several
researchers such as a successive interference cancellation approach in [11] (see also [12], [13] and
their references). A detailed survey of FTN transmission for different system parameters and its
implementation aspects can be found in [4] where it is stated that such a transmission approach
can provide up to twice the spectrum efficiency compared to the existing Nyquist transmission
(i.e., 2 times the rate achieved by the single input single output (SISO) system).
The FTN transmission can be realized both for single-antenna and multi-antenna systems [14].
It is shown that for the single input multiple output (SIMO) and MIMO channels, the FTN
transmission also improves the performance of the system [15]. In fact, the FTN allows the
transmitter to send data symbols at a rate higher than the Nyquist transmission which can also
be interpreted as a novel approach to exploit the additional degrees of freedom of the channel
(i.e., ”multiplexing gain”). This recalls us the two fundamental aspects of a wireless channel (i.e.,
diversity and multiplexing) which were studied in the existing MIMO communication literature
[16], [17]. The multi-antenna system can be designed as a SIMO, MIMO and MISO (multiple
input single output). It is well known that when either the transmitter or receiver has single
antenna (i.e., MISO and SIMO), the communication system can achieve only a diversity (power)
gain [16]–[18].
For MIMO systems, one can achieve both diversity and multiplexing gains under favorable
channel conditions (i.e., it is possible to transmit and reliably detect on a symbol-by-symbol basis
more than one symbol per Nyquist period). On the other hand, it has been shown in [16] that there
is a fundamental tradeoff between the diversity and multiplexing gains in a MIMO system (i.e.,
one can increase diversity (multiplexing) gain while decreasing multiplexing (diversity) gain).
A detailed analysis of the MIMO system performance for different parameter setups and signal
to interference plus noise ratio (SINR) levels for flat fading channels can be found in various
wireless communication books (see Chapters 5-7 of [16] for extensive study).
The design of efficient transmission strategies for flat fading channels can be extended almost
straightforwardly to the frequency selective channels. This is due to the fact that any signal
transmitted through a frequency selective channel can be splitted into orthogonal chunks of signals
(i.e., sub-carriers) each experiencing a flat fading channel by using the orthogonal frequency
division multiplexing (OFDM) transmission approach. By employing the OFDM scheme, the

achievable rate per Hertz of a SIMO system can be expressed as [16]
REx

Ns
1 X
i
log2 (1 + SN REx
(N )) bps/Hz
=
Ns i=1

(1)

where Ns is the number of sub-carriers which increases with the transmission bandwidth B Hz for
i
fixed sub-carrier spacing, SN REx
(N ) is the signal to noise ratio (SNR) achieved at sub-carrier

i of the OFDM transmission approach, which is a function of the number of receive antennas
N . The SNR achieved at each sub-carrier depends on different factors ranging from the number
of antennas to the channel condition. In a multi-antenna system, REx is usually quantified by
assuming that the channel between any transmit receive pair experiences independent fading [16].
Also, the deployment of large antenna arrays at the base station (BS) has been prototyped and
studied where recent measurement results [19]–[24] suggest that despite the statistical difference
between the measured channels and the independent and identically distributed (i.i.d) channels,
most of the theoretical conclusions made under the independence assumption are still valid for
the realistic massive MIMO channels1 . Now when the channel between the transmitter and each
receiver has P unit-variance independent multipath components, the existing OFDM transmission
can have the following upper-bound average rate [17]2


NP
REx = log2 1 + 2
σ

(2)

where σ 2 is the noise variance. As we can see from this expression, increasing the number of
receiver antennas N enables to enhance the channel achievable rate (2) logarithmically only [17],

[26].
Motivation and Objective: The FTN transmission can be realized for different systems such
as massive MIMO, mmWave and full-duplex wireless systems. Thus, such a transmission strategy
will undoubtedly improve the overall capacity of the future 5G network. However, the existing
FTN transmission can improve the overall point-to-point achievable rate at most two times only3 .
For instance, if we consider a SIMO system, the capacity of FTN transmission will be at most
1
2

This is also stated as a ”favorable propagation condition” which is commonly adopted in MIMO channels [16], [17].


2
As can be seen in (1), the accurate rate expression of the ith sub-carrier is Ri = log2 1 + |hσoi2| , and the true average

value of Ri is given as Riav = E{Ri } where hoi is the channel between the transmitter and receiver of sub-carrier i. Since the
exact analytical expression of Riav is difficult to derive for most practically relevant statistical models of hoi (for example, i.i.d
case), we use the well known Jensen’s inequality to get the upper-bound expression for Riav by replacing |hoi |2 with its average

value [25].
3
This can be translated to a maximum of 2 times system capacity compared to the Nyquist transmission.

2REx . Furthermore, the receiver may need to use a non-linear (joint) decoding approach which is
often not desirable in practice. This motivates us to examine the following problem: Is it possible
to employ a FTN transmission approach while ensuring that the point-to-point link achievable
rate is improved more than 2 times with a simple symbol-by-symbol decoding? If yes, in what
scenario and how can we enable the symbol-by-symbol decoding?
In the current paper, we propose a novel FTN transmission and symbol-by-symbol detection
approach for SIMO systems over wideband channels (i.e., frequency selective fading channels).
The proposed transmission scheme exploits the inherent multipath characteristics of wideband
channels. We show that when the channel has P independent multipath components (i.e., a
favorable propagation condition) which is commonly adopted in different MIMO channels [16],
[17] and the receiver has multiple antennas N ≥ P , one can achieve the following upper-bound
rate4

RP r




P2
3P − 2 N − P − 1
=
log2 1 +
.
3P − 2
P2
σ2

(3)

For fixed P and σ 2 , one can notice that RGap = RP r −REx increases as N increases which renders

the proposed approach appealing for large receiver antenna array systems. The complexity of the
proposed approach is almost the same as that of the OFDM transmission approach and it is
simple to implement. As will be clear in Section IV, the proposed design allows the estimated
symbols to be decoded independently without experiencing inter-symbol-interference (ISI). For
this reason, one can consider the proposed transmission strategy as orthogonal faster than Nyquist
(OFTN). The rate expressions RP r and REx have been given for correlated multipath channel
components as well. The proposed approach also does not suffer from the peak to average power
ratio (PAPR) problem which provides additional advantage compared to OFDM. On the other
hand, one can straightforwardly extend the proposed transmission approach to multiuser (massive)
MIMO systems.
We would like to emphasize here that the proposed approach relies solely on the number and
characteristics of multipath components, and number of antennas at the receiver side. Moreover,
the proposed OFTN transmission exploiting the multiplexing gain which is incorporated in multipath channels has never been revealed in the existing works (to the best of our knowledge).
4

As will be clear in the simulation section, RP r and REx are almost the same as their true average values for practically

relevant N and P .

This motivates us to use the term multipath multiplexed transmission. Indeed, different possible
bandwidths are utilized in the current long term evolution (LTE) network (see [27], [28] for
more details). For instance, if we consider a B = 2.5 MHz bandwidth, the number of multipath
components is around 18 which is significant [29]. Increasing B would increase the number of
multipath components in a given transmission environment.
This paper is organized as follows. Sections II and III discuss the considered signal and channel
models, and a brief summary of the OFDM and problem formulation of the proposed design,
respectively. The detailed descriptions of the proposed OFTN transmission and detection approach,
and its performance analysis are presented in Sections IV and V, respectively. The effect of power
adaptation on the performance of OFDM and OFTN designs, and some of the practical issues
of the OFTN scheme including its similarities with that of FTN is explained in Sections VI VII. Simulation results are provided in Section VIII. Finally, some open research problems and
conclusions are given in Sections IX and X, respectively.
Notations: In this paper, upper (lower-case) boldface letters denote matrices (column) vectors.
X(i,j) , XT , XH and E(X) denote the (i, j)th element, transpose, conjugate transpose and expected
value of X, respectively. The convmtx(.), 1N , I and CN ×M (<N ×M ) denote convolution matrix,
an N sized vector of ones, identity matrix and N × M complex (real) entries, respectively.

Notations: In this paper, upper (lower-case) boldface letters denote matrices (column) vectors.

X(i,j) , XT , XH and E(X) denote the (i, j)th element, transpose, conjugate transpose and expected
value of X, respectively. The convmtx(.), 1N , I and CN ×M (<N ×M ) denote convolution matrix,
an N sized vector of ones, appropriate size identity matrix and N × M complex (real) entries,

respectively.

II. S IGNAL AND C HANNEL M ODEL
We consider a single user system where the transmitter has one antenna whereas, the receiver
has N  1 antennas. The transmitted signal has a transmission bandwidth B which is sufficiently

large making the channel to be frequency selective (i.e., wideband signal). The channel impulse
response between the transmitter and the nth receive antenna is denoted by hn (t) which is assumed
to be linear time invariant. In the following, we provide the discrete time baseband equivalent of
hn (t). For better exposition, we first consider a SISO channel (i.e., ignoring the subscript n) and
the result is extended further to a SIMO channel.

TABLE I: Frequently used variables, vectors and matrices
Var.

Definition

Q

Spatial channel correlation matrix

R

Temporal channel correlation matrix

N

Number of receive antenna

S

Number of transmitted symbols in 3P-2 sampling period

P

Number of multipath components of the channel impulse response

L

Number of multipath taps in OFDM transmission (L = P when g(t) = rect(t))

g(t)

Transmitter and receiver pulse shaping filter in existing OFDM (i.e., g(t) = rect(t) of width Ts

g̃(t)

Transmitter and receiver pulse shaping filter of proposed design (i.e., g̃(t) = rect(t) of width LT̃s

Q

Width of the transmitter and receiver pulse shaping filters (Q = L = P when g̃(t) = rect(t))

B

Signal bandwidth (B =

Td

Channel delay spread

T̃s

Temporal spacing of equally spaced independent channel coefficients (T̃s =

Tb

Block duration of OFDM

Ns

Number of sub-carriers of OFDM (it is assumed that all sub-carriers are used)

s[k]

Transmitted data signal in kth symbol period of OFDM

di , ∀i

Transmitted symbols of proposed approach which are temporally separated by T̃s

1
)
Ts

Td
)
P

yl

Received sample at time slot l obtained at each antenna (i.e., (24))

H̄I
¯
H̄

Equivalent random channel between transmitter and all receivers when yl of all antennas are stacked
¯ = Q1/2 H̄ R1/2
Overall channel between transmitter and all receivers when yl of all antennas are stacked H̄
I

hi , ∀i

Discrete time channel impulse response coefficients with temporal spacig T̃s

If we transmit an arbitrary baseband waveform x(t) having bandwidth B and carrier frequency
fc , the received signal in baseband form can be expressed as [17], [30]
y(t) = h(t) ∗ x(t) + w(t)

(4)

where h(t) is the baseband equivalent channel impulse response, w(t) is the additive noise and
(∗) is the convolution operator. In general, one can model h(t) as (see (3.10) of [30], (2.19) of
[17] and [31])
h(τ ) =

P
−1
X
i=0

αi ejφi (τ ) δ(τ − τi )

(5)

where P is the total number of multipath components, and αi , φi and τi are the attenuation
(due to pathloss and shadowing which is varied quite slowly), phase and propagation delay of
the channel from the transmitter to the receiver in path i, respectively. Here we assume, without
loss of generality, that each multipath component corresponds to a single reflector/scatterer and

TABLE II: LTE Extended Typical Urban model (ETU)
Delay (ns)

0

50

120

200

230

500

1600

2300

5000

Relative PDP (dB)

-1

-1

-1

0

0

0

-3

-5

-7

φi (τ ) = 2πfc τi + φfDi (τ ) with φfDi (τ ) = 2πτ fDi , fDi =

v
λ

cos θi (i.e., Doppler shift) and θi is

the angle of arrival (AOA) of path i. As can be seen from the above expression, a small change
in τi can lead to a significant variation in the phase φi . Thus, the rate of change of φi can be in
the order of very few milliseconds [17], [30].
In general, the channel impulse response (5) is a random process due to the random amplitude,
phase and delay components. Thus, h(t) can be characterized only statistically or from experimental results. In a practical setup, the channel response corresponding to τi is uncorrelated with
that of τj , j 6= i. This is mainly due to the fact that these two channel responses are the results of

different scatterers (i.e., uncorrelated scattering) [30] (see page 181). Furthermore, as fc τi =

di
λ

with di  λ, where di is the propagation distance of the ith path, φi can be reasonably modeled

as a uniformly distributed random variable U[0, 2π] [17], [30]. The power delay profile (PDP)
of a multipath channel can therefore be given as

Ah (τi , τj ) = E{h(τi )h(τj )H } = E{αi αj }δ|i − j| , Ah (τi ) = αi2 .

(6)

In a wireless channel, the delay spread Td is defined as the time index τi where Ah (τi ) ≈ 05 .

Different standards suggest different delay parameters and PDP’s for different environments. For
instance, the Extended Typical Urban model (ETU) of the LTE network uses the specification
given in Table II [32]. As can be seen from this table, the LTE ETU channel model contains
P = 9 independent paths where each of them has a different PDP (see also [33] for other channel
models). Furthermore, one can observe that the time delays may not be spaced uniformly.
Thus, without loss of generality, one can represent the discrete time baseband equivalent form
of h(t) as
h = [h1 , h2 , · · · , hP ]T = R1/2 h̄

(7)

where h̄ ∈ C P ×1 is a vector containing independent channel coefficients (i.e., as in Table II),

R ∈ C P ×P is a temporal correlation matrix and h1 , · · · , hP are equally spaced channel coefficients
5

The multipath delay spread Td can also be defined as the difference in the propagation delay between the longest and shortest

path having a significant energy, i.e., Td = maxi,j |τi − τj |.

with temporal spacing given by T̃s ,

Td
P

(for example Td = 5µs and P = 9 if we apply the

ETU channel model given in Table II). From these explanations, one can understand that Td is
determined by the effects of reflector and scatterers which are mainly related to the propagation
environments. Specifically, one may expect to experience very large Td in a densely populated
urban areas, and very small value of Td in an indoor environment [17], [27], [28], [30], [34]. The
coherence bandwidth which explains how fast the channel changes in frequency is given by
Wc =

1
.
Td

(8)

The channel coherence bandwidth Wc can also be treated as the bandwidth in which the channel
is considered as ’almost constant’. Now when B  Wc , the wireless channel is termed as a

frequency selective channel (or wideband channel) which is the focus of the current paper. When
the receiver has N antennas, one can model the channel coefficient between the transmitter and
each receiver with (7).
III. E XISTING OFDM T RANSMISSION AND P ROBLEM S TATEMENT
This section provides a brief summary of the well known OFDM transmission and formally
states the considered problem.
A. Summary of OFDM Transmission
When the channel is frequency selective, OFDM is known to be an efficient transmission
strategy to maximize the achievable rate [17]. In fact, OFDM can be used to harness the effect
of ISI that the transmitted signal experiences when B  Wc . For brevity, this section briefly

summarizes the main idea of OFDM transmission and reception. In the existing transmission
approach, the baseband transmitted signal x(t) can be expressed as
x(t) =

X
k

(9)

s[k]g(Bt − k)

where s[k], ∀k are the transmitted symbols with period Ts =

1
B

and g(t) is the pulse shaping

filter. In most wireless systems, g(t) is designed as a rectangular, raised cosine or sinc function.
In this section, we also start with a SISO system (i.e., by dropping the antenna subscript n) and
then extend the analysis for SIMO system. When the receiver has single antenna, the received

signal y(t) can be expressed as
y(t) = h(t) ∗ x(t) + w(t) =

X
k

s[k]

P
−1
X
i=0

(10)

hi g(B(t − iT̃s ) − k) + w(t).

The recovered signal at the mth symbol period y[m] = y(mTs ) can be rewritten as (see (2.35)
of [17])
y[m] =

L−1
X
l=0

fl s[m − l] + w[m] =

L−1
X
l=0

(11)

h̃L−1−l s[m − l] + w[m]

where w[m] is the noise sample and


P
−1
X
T̃s
fl =
hi g l − i , [h̃0 , h̃1 , · · · , h̃L−1 ] = [fL−1 , fL−2 , · · · , f0 ].
Ts
i=0

(12)

As can be seen from (11), y[m] contains not only s[m] but also additional interference terms

s[m − l], ∀l, where l depends on the channel and the support width of g(t). If we consider
g(t) = rect(t) (i.e., g(t) = 1, t ∈ [0, Ts ] and g(t) = 0, t > Ts ), we will have L =

T̃s P
Ts

=

Td
Ts

=P

(i.e., the number of multipath taps) [35]. Now if the OFDM signal has a block length Tb = Ns Ts ,
the received samples y[1], y[2], y[3], · · · , y[Ns ] can be stacked in vector form y which can be
expressed as

(13)

y = H̃s + w

where s = [s[1], s[2], · · · , s[Ns ]], w = [w[1], w[2], · · · , w[Ns ]] and H̃ is a Toeplitz matrix formed
by using the channel coefficients [h̃0 , h̃1 , · · · , h̃L−1 ] [17], [18], i.e.,

0
0
0
0
0
0
0
0
0
 h̃L−1

 h̃L−2 h̃L−1 0
0
0
0
0
0
0
0

 .
.
.
..
..
 ..
0
0
0
0
0
0
0


 h̃0
h̃1 h̃3 · · · h̃L−2 h̃L−1
0
0
0
0


H̃ = 
h̃0 h̃1 h̃3 · · · h̃L−2 h̃L−1
0
0
0
 0

 0
0
h̃0 h̃1 h̃3
· · · h̃L−2 h̃L−1
0
0


 0
0
0 h̃0 h̃1
h̃3
· · · h̃L−2 h̃L−1
0

 .
.
..
 .
..
.
 .

0
0
··· 0
h̃0
h̃1
h̃3
· · · h̃L−2 h̃L−1














.












(14)

Since H̃ is a full-rank matrix, the receiver can reliably estimate s by utilizing appropriate

equalization techniques (for example with a zero forcing (ZF) approach which is a simple

channel inversion technique). However, such an approach creates undesirable effects such as
noise enhancement.
The existing OFDM transmission and reception approach aims to circumvent this drawback,
and if we append the last L symbols of s in the beginning of the transmission (i.e., the transmitted
signal becomes [s[Ns − L + 1], s[Ns − L + 2], · · · , s[Ns ], s]) and if the receiver discards the first

L received samples (i.e., cyclic prefix removal), the input-output relation can be expressed as

(15)

ỹ = H̃c s + w

where H̃c becomes a circulant matrix which has a number of important characteristics. Specifically, H̃c can be equivalently expressed as H̃c = FH GF, where F is a discrete Fourier transform
(DFT) matrix of size Ns and G is a diagonal matrix where its elements are computed from h̃l , ∀l.

Due to this reason, one can precode s as s̄ = FH s before transmission (i.e., inverse fast Fourier

transform (IFFT) operation) and then perform post processing the received signal with F which
consequently yields [17], [18]
ỹ = FH GFFH s + w ⇒ Fỹ = Gs + w̃ ⇒ ŝ[m] = s[m] +

g[m]∗ w̃[m]
, m = 1, 2, · · · , Ns (16)
|g[m]|2

where w̃ = Fw = [w̃[1], w̃[2], · · · , w̃[Ns ]] and each entry of w̃ is still an i.i.d zero mean circularly

symmetric complex Gaussian (ZMCSCG) random variable with variance σ 2 , g[m] is the mth

diagonal element of G (i.e., the channel gain corresponding to the mth symbol) and ŝ[m] is the
estimate of s[m]. As we can see Gs is a parallel channel where each element of s does not
experience any ISI signal. Since F is a Fourier matrix, one can interpret s̄ as the time domain
version of the frequency domain signal s. For this reason, the precoding and post processing
operations of the OFDM transmission are usually stated as a signal transformation from frequency
domain to time domain and vice verse, respectively.
Now when the receiver has N antennas (SIMO system), each antenna will transform its received
signal into a parallel channel as in (16) to coherently combine ŝ[m] with the maximum ratio
combining (MRC) approach. By doing so, one can express the estimate of s[m] as
ŝ[m] = s[m] +

˜
g̃[m]H w̃[m]
, m = 1, 2, · · · , Ns
|g̃[m]|2

(17)

˜
where g̃[m] = [g1 [m], g2 [m], · · · , gN [m]] and w̃[m]
= [w̃1 [m], w̃2 [m], · · · , w̃N [m]]. As we can

see from the above expression, the existing OFDM transmission scheme is able to transmit Ns

symbols in Tb , (Ns + L)Ts seconds. And the SNR of each symbol depends on N ; increasing
(decreasing) N increases (decreases) the SNR of the symbol s[m].
B. Problem Statement
To better present the objective of the work, we consider the time-domain input-output expression
given in (13). When N ≥ L and hn = [hn0 , hn1 , · · · , hn(L−1) ] (i.e., the elements of the channel
vector defined in (7) corresponding to the nth receive antenna) are independent random variables,
the received samples obtained at time LTs in all antennas can be expressed as
˜
yL = H̃s
(1:L) + wL

(18)

˜ = [h̃ ; h̃ ; · · · ; h̃ ], w = [w [L], w [L], · · · , w [L]] and s
T
where H̃
1
2
N
L
1
2
N
(1:L) = [s[1], s[2], · · · , s[L]]

with h̃n = [h̃n0 , h̃n1 , h̃n2 , · · · , h̃n(L−2) , h̃n(L−1) ] is the channel vector defined in (12) corresponding
to the nth receive antenna. Similarly, the received samples at 2LTs can be expressed as
˜
y2L = H̃s
(L+1:2L) + w2L .

(19)

The expressions (18) and (19) have the same mathematical structure as that of the MIMO channel.
Furthermore, one can decode s(1:L) and s(L+1:2L) from (18) and (19) efficiently using different
existing MIMO detection schemes such as ZF, MRC etc. This demonstrates that the received
samples yi , i 6= kL are utilized just to improve the SNR of the transmitted symbols yi , i = kL.

In fact, the relation between upper-bound rate, SNR, Tb and B is




Lγ0
Lγ0
Ns 1
R = log2 1 + 2 ≈
log2 1 + 2
σ
Tb B
σ

(20)

where γ0 is the equivalent SNR of the transmitted symbols si , ∀i when they are decoded by

employing yi , i = nL only. Furthermore, for fixed Tb and B, as each sub-carrier transmits only
one symbol, the total transmitted symbols in an OFDM frame is Ns (i.e., number of transmitted
symbols is the same as that of sub-carriers). Thus, a linear increase in Ns (number of transmitted
symbols) increases the rate ”linearly” whereas, a linear increase in SNR increases the rate ”only
logarithmically”. Therefore, for a given Tb and B, increasing Ns is more advantageous than
increasing Lγ0 .
The above discussions show that the existing transmission and decoding approach employs
yi , i 6= kL just to improve SNR which is not the best strategy from the achievable rate point of

view (i.e., the benefits of the time slots yi , i 6= kL are not exploited efficiently). This motivates

us to set the following objective: For the given symbol constellation, B and Tb , is it possible to
transmit more than Ns symbols and reliably decode all of the symbols in a symbol-by-symbol
basis like in the OFDM approach while achieving better achievable rate than that of the existing
OFDM transmission approach? In other words, can we take advantage of the inefficiently used
symbol time slots i 6= kL, ∀k to transmit more data symbols while ensuring that the rate is higher

than that of the existing OFDM approach. If yes, how can the transmission strategy be designed
and to what extent can the rate be increased?
IV. P ROPOSED OFTN T RANSMISSION A PPROACH
This section describes the proposed OFTN transmission and data detection approach for SIMO
wireless systems to address the above objective. The extension of the proposed design for different
system setups will be discussed in the next section.
A. Transmission Scheme
One can observe from the discussions in the above section that the OFDM transmission scheme
transmits only one symbol in Ts seconds. However, as our objective is to transmit more than one
data symbol in Ts seconds, we may need to find a way to incorporate additional data symbols
per Ts . In this regard, we split the overall bandwidth B into L equivalent sub-bands where each
of them has bandwidth of BL ,

B
L

(i.e., traditional frequency division multiple access where the

symbol period of each of the sub-bands is the same as the delay spread of the channel)6 . With the
bandwidth BL , in fact, the existing OFDM transmission approach is able to transmit NL ,

Ns
L

data symbols in Tb seconds. Our objective will now be to transmit more than NL data symbols
in the same bandwidth which is detailed as follows.
For better exposition, we assume that the transmitted symbols, d0 , d1 , d2 , · · · are mapped into

period T̃s , and the transmitter and receiver filters employ bandwidth BL =

1
Td

and have a support

duration of Q̃ measured in terms of T̃s 7 . As in the above section, we start with the SISO system,
and the result is then extended to SIMO. With the above settings, we will have the following
6
7

This can be performed by employing an appropriate scheduling algorithm.
Note that this can be interpreted as oversampling of the transmitted signal P times where, unlike the traditional oversampling

with zero padding, the oversampled signals also contain independent information symbols which will be clear in the sequel.

transmitted data signals at time instant lT̃s
Q̃

xl , x(lT̃s ) =

2
X

k=− Q̃
2

g̃(lT̃s − k T̃s )dl+k .

(21)

We would like to emphasize here that the bandwidth of the transmitted signal x(lT̃s ) is determined
by the bandwidth of the filter g̃(t) which is still BL . For better understanding of the proposed
transmission, we assume that g̃(t) is a rectangular pulse shaping filter of width Td . With g̃(t) =
rect(t) (i.e., P = Q̃) and after appropriate time shifting, x(lT̃s ) can be re-expressed as
l
X

xl =

(22)

dk .

k=l−P +1

From this expression, one can understand that the only difference between the existing OFDM
transmission scheme and the proposed scheme is that the existing approach sets dk = 0, k 6= l

(i.e., zero padding approach) or dk = f (dl−P , dl , dl+P , · · · ), k 6= l (i.e., correlated dk in the
interpolation approach). In the following, we show that it is possible to introduce independent dk

and decode (estimate) them reliably, where the number of introduced independent symbols will
depend on N and the number of independent samples of h which will be clear in the sequel.
The received signal is given by the convolution of the transmitted signal and the channel which
can be expressed as
rl =

P
−1
X

(23)

hm xl−m + w̃l

m=0

where w̃l is the additive noise. Once the received signal is passed through the receiver matched
filter (i.e., a rect function) and after doing some mathematical manipulations, we will have
yl =

l
X

(rk + w̃k ) =

k=l−P +1

l
X

P
−1
X

k=l−P +1 m=0

=hT Adl + wl = h̄T RT Adl + wl
where wl =

Pl

k=l−Q̃+1

hm xk−m + wl
(24)

w̃k which is assumed to be a ZMCSCG random variable with variance

σ 2 , dl = [dl , dl−1 , · · · , dl−3P +4 , dl−3P +3 ]T , the last equality follows from (7) and A = A0 A1

with A0 ∈ <P ×2P −1 (A1 ∈ <2P −1×3P −2 ) is a convolution matrix obtained by utilizing a P row
vector of ones (i.e., A0 = convmtx(1TP , 2P − 1) and A1 = convmtx(1TP , 3P − 2)). As we can

see from (24), yl is a linear combination of 3P − 2 symbols (i.e., dl−3Q̃+3 to dl ) which has the
same dimension as the convolution of the transmit and receive filters, and the channel vector.

B. Number of Reliably Decoded Symbols
When the receiver has N antennas, one can stack yl of all antenna elements to get
¯ Ad + w
yl = H̄T R1/2 Adl + wl = Q1/2 H̄I R1/2 Adl + wl = H̄
l
l

(25)

where Q ∈ C N ×N is the spatial channel correlation matrix of all receive antennas, H̄ = [h̄1 , h̄2 , · · · , h̄N ]T =

(Q1/2 H̄I )T , each entry of H̄I is assumed to be an i.i.d ZMCSCG random variable with unit
¯ , Q1/2 H̄ R1/2 and w = [w , w , · · · , w ] with h̄ as h̄ of (7) corresponding to
variance, H̄
I
l
1l
2l
Nl
n

the nth antenna. One can model different propagation characteristics of a wireless environment

by setting different values of Q and R. For instance, a communication environment having few
scatterers can be represented by setting a very low rank R. Furthermore, when the antenna spacing
is not sufficiently large to experience independent channel coefficients, the spatial correlation
matrix Q will be different from a diagonal matrix. For these reasons, we believe that (25) can
be utilized for different wireless environments just by tuning Q and R.
Since (25) can be considered as a typical MIMO system, the decodability of dl with a linear
¯ A. And since,
decoder (estimator) and symbol-by-symbol basis entirely depends on the rank of H̄
¯ A) ≤ P , the number of independent symbols of d cannot be more than P . When the
rank(H̄
l

maximum number of independent transmitter symbols (S) is less than P , one can set the first S

coefficients of dl to be the desired symbols while leaving the others 0. By doing so, A becomes an
¯ A) = rank(H̄
¯ ) = min(rank(Q), rank(R)).
upper triangular matrix where in such a case rank(H̄
Therefore, we can express the maximum number of independent transmitted symbols as
S = min(rank(Q), rank(R)) ≤ P.

(26)

From this expression, we can understand that the maximum number of transmitted symbols
depends on the ranks of the channel covariance matrices Q and R. Furthermore, when these
two matrices have full ranks (i.e., favorable propagation conditions), we can reliably estimate P
independent data symbols of dl . The advantage of the proposed transmission approach compared
to that of the existing one (i.e., OFDM) is that it allows the transmission of S ≤ P symbols

in (3P − 2)T̃s interval whereas, the existing approach transmits only 3 symbols in the same

duration. Thus, in a favorable propagation environment, the proposed transmission scheme is
attractive especially when P and N are large.
To show the impact of N on the decodability of dli pictorially, we set Q = I, R = I and
considered binary phase shift keying (BPSK) transmitted symbols dli . Fig. 1 shows the eye

Eye Diagram (N=7)

Eye Diagram (N=8)

1

10 0

1

N=6
N=7
N=8

−0.5

BER

0

−1
−0.5

10 -1

0.5
Amplitude

Amplitude

0.5

0

−0.5

0
Time

0.5

−1
−0.5

10 -2

10 -3

0
Time

0.5

10 -4

0

5

10

15

20

25

30

35

SNR (dB)

Fig. 1: Eye diagram of proposed design Fig. 2: The uncoded BER of proposed dewhen P = 8 and γ → ∞.

sign for P = 8 and N = [5, 6, 7].

diagram of the estimate of dli obtained by applying ZF receiver for different values of N and
γ → ∞ (SNR) which is defined as γ =

E|dli |2
σ2

while setting P = 8. One can notice from this

figure that the eye diagram of the estimated data resembles that of the true one (i.e., with values
±1) when N = P . However, when N < P , the shape of the eye diagram is not improved even

if the SNR goes to infinity. This is due the fact that in such a case the received signal does not
have sufficient degrees of freedom (DoF) to reliably detect the transmitted symbols. This implies
that the maximum number of transmitted symbols in Td seconds will be min(N, P ). To verify
this, we have also plotted the uncoded bit error rate (BER) of the proposed transmission approach
for different values of N as shown in Fig. 2. From this figure, one can also observe that when
N < P , there is a BER floor and increasing SNR does not help avoid this undesirable effect8 .
However, the BER decreases consistently with SNR when N ≥ P .
V. DATA D ETECTION AND P ERFORMANCE A NALYSIS
As can be noticed from the discussions of the previous section, we are able to transmit S
symbols in (3P − 2)T̃s . This section discusses some of the possible data detection techniques to

recover dl from yl . In this regard, we discuss three possible approaches, namely MRC, ZF and
zero-forcing successive interference cancellation (ZF-SIC) as follows.
8

Note that the eye diagram of Fig. 1 (right side) is not exactly the same as that of the true one. This is because some of

the channel realizations of H̄ are not well conditioned which creates some non-zero residual interference on the estimate of dli .
However, such residual interference will not fundamentally affect the achievable BER (see Fig. 2).

From the existing OFDM transmission, one of the most widely used simple receiver designs
is the MRC. For the proposed OFTN transmission, one can design the MRC beamformer as
1/2 †
B = A−1 (R1/2 )† H̄H
).
I (Q

(27)

And for the ZF beamformer, we can design B as
B = (Q1/2 H̄I R1/2 A)† .

(28)

As can be seen from these two expressions, both the MRC and ZF beamformings have matrix
inversion operations. However, for the MRC approach A−1 , R† and Q† can be computed offline
as they are constant design coefficients. For this reason, the MRC beamforming does not need
to perform inverse operation when the channel is updated. However, the ZF beamformer indeed
requires such an operation whenever the channel is varied (i.e., once per channel coherence time).
For this reason, MRC beamforming is computationally less expensive than the ZF counterpart.
With these beamforming matrices, we can recover dl as
(29)

d̂l = Byl .

The main drawback of the MRC and ZF beamforming approaches is both of them require to
perform the matrix inverse operation on A which is an upper triangular matrix. For example,
when S = P = 4, we will get (30)

A(4×4)



1 2 3 4







 0 1 2 3 

.
=

 0 0 1 2 


0 0 0 1

(30)

Consequently, the recovered symbols d̂l experience quite different desired signal power which
consequently leads to unbalanced SINR and lower total achievable rate. To alleviate this drawback,
we propose a combined ZF and SIC (i.e., ZF-SIC) based data detection approach which employs
two steps: In the first step, we compute
¯ H H̄
¯ )−1 H̄
¯ H y = Ad + w
d̂l = (H̄
l
l
l

(31)

¯ H H̄
¯ )−1 H̄
¯ H w̄ . Now since the Sth element of d̂ is interference free, it can be
where w̄l = (H̄
l
l
decoded independently. In the second step, we subtract dlS from the (S − 1)th element of d̂l and

decode (i.e., dl(S−1) ). This is because the (S − 1)th element of d̂l contains only one interference

term which is known (i.e., dlS ). Finally, we repeat these two steps until all symbols are decoded.

The remaining question is to assess whether the proposed transmission strategy will always
achieve better performance than that of the OFDM or not? Addressing this question for a general
Q and R appears to be intractable. In the following, we analyze the average performances of the
proposed and OFDM transmission schemes for some practically relevant parameter settings.
Theorem 1: When Q = I and R = I, the proposed ZF-SIC based approach (assuming no error
propagation) and existing OFDM approach can have the following upper-bound average rates
REx = log2 (1 + γ Ex ), γ Ex =

LN
σ2

˜ log2 (1 + αγ P r ), γ P r =
RP r =Lα̃

N −L+1
σ2

˜ =
where RP r (REx ) is the rate achieved by the proposed (existing) approach, α̃
3L−2
.
L2

(32)
L
3L−2

and α =

˜ is introduced as the contribution of any L transmitted symbols in BL
The coefficient α̃

bandwidth will span up to (3L − 2)T̃s durations due to the channel delay, transmit filter and

receive filter, and L is introduced to ensure that the proposed OFTN approach will utilize a total
bandwidth B = LBL . And α is the power normalization factor which is introduced to maintain

that both the OFDM and OFTN approaches use the same average transmit power. Furthermore,
if

N
σ2

 1 and N  L (i.e., at high SNR and large N ), we will have RP r − REx ≥ 0 when
1

N ≥ max{L, 2 L2 −3L+2

[(3L−2) log2 (L)−L2 log2 (α)] 2

σ }.

(33)

Proof: See Appendix A.
The proposed design attempts to model the discrete time equivalent impulse response of the
channel using (7). In fact, the number of multipath channel coefficients definitely increases with
bandwidth. For instance, in the ultra wideband communication, one can have L in the order of
200 − 500, where the channel coefficients often tend to be sparse [17]. The considered model

captures the channel characteristics of such an environment by selecting Q and R appropriately.
This justifies the need to examine the performance of the proposed design for arbitrary correlation
matrices. This analysis is not trivial mainly because the SNR expression incorporates an inverse
matrix containing random variables which is difficult to handle. However, when N is large, by
utilizing the law of large numbers, most of the random components become deterministic which
enables us to study the system performance for large N in the following lemma.

Lemma 1: When N  L (i.e., for a large antenna array), the existing and proposed approaches

achieve the following upper-bound average rates
R

Ex

Pr
RZF

Pr
RZF
−SIC

Ns
1 X
N
=
log2 (1 + γiEx ), γiEx = 2 fiT RRT fi∗
Ns i=1
σ

S
N
11X
log2 (1 + α̃γiZF ), γiZF = 2
≈
α̃ S i=1
σ E(i,i)

S
11X
N
≈
log2 (1 + α̃γiZF −SIC ), γiZF −SIC = 2
α̃ S i=1
σ C(i,i)

(34)

Pr
Pr
where REx , RZF
and RZF
−SIC are the rates achieved by the existing, proposed ZF and proposed

ZF-SIC approaches, respectively, α̃ =

3L−2
,
LS

C = (RH R)−1 and E = (AH RH RA)−1 . As

expected, when Q = I, R = I and N  L, the rate expressions of (34) is almost the same as

that of (32) since N ≈ N − L.
Proof: See Appendix B.

From the results of Theorem 1 and Lemma 1, the following points can be highlighted:
1) In (32), as γ Ex ≈ Lγ P r for L ≥ 3 and medium to large N , one can interpret this equation

that the existing OFDM approach attempts to ensure that each transmitted symbol has a
maximum SNR (i.e., diversity gain achieving strategy). However, the proposed approach
increases the number of transmitted symbols while reducing the SNR achieved by each
symbol (i.e., multiplexing gain achieving strategy). Nevertheless, with the proposed approach,
it is always possible to decrease S while increasing the SNR, for instance by reducing S to

half, the SNR can be doubled. To show this fact, we have plotted the achievable γEqt = αγ P r
for different values of S and N in Fig. 3. As can be seen from this figure, the SNR of each
symbol can be increased just by decreasing S. This demonstrates that there exists a multipath
diversity-multiplexing tradeoff for SIMO systems. Having said this, however, exploiting this
tradeoff for general coding schemes (for example with Alamouti space time block coding
(STBC) [36]) is an interesting open problem which is left for future research.
We would like to emphasize here that the existence of a diversity-multiplexing tradeoff has
been known in the existing transmission scheme when both the transmitter and receiver are
equipped with multiple antennas (i.e., MIMO) [16]. However, to the best of our knowledge,
we are not aware of any existing work demonstrating the diversity-multiplexing tradeoff
which is exhibited by exploiting multipath components for a SIMO system.

2) From (33), one can notice that the proposed approach achieves better performance at large
N . Furthermore, for fixed L and σ 2 , increasing N increases the gap between the proposed
and existing designs. This shows that the proposed design is practically useful especially for
massive MIMO systems. Nevertheless, the proposed approach can still be customized for
small N by optimizing S. However, optimizing S for a given N (possible small) is still an
open research topic and is left for future research direction.
3) As the proposed algorithm does not require IFFT operation at the transmitter, it does not
suffer from high PAPR which merits the benefit of proposed design compared to OFDM.
Furthermore, the main computational load of the proposed algorithm is the matrix inversion
operation of (28) which is comparable with the OFDM transmission.
4) The noise samples obtained from (25) are i.i.d which consequently avoids the need of prewhitened matched filter which has been employed in the existing FTN transmission schemes
[4], [37]. The disadvantage of the pre-whitening operation is that it often amplifies adjacent
channel interference signals which consequently worsens the achievable BER of the FTN
transmission. Furthermore, the pre-whitening phase enables each data symbol to experience
ISI symbols which consequently requires complicated decoding operation especially when
the number of ISI symbols is large which is not desirable in practice [4].
5) One can also observe from (25) that the desired information is obtained every 3P − 2

sampling periods. This behavior will also allow to use low cost analog to digital converter
(ADC) which is one of the most expensive parts of a receiver unit especially in massive
MIMO regime [38]–[40].

Although the proposed ZF-SIC based data detection approach maintains the desired symbols of
dl to have balanced signal power, the latter approach may suffer from error propagation. In other
words, if dlS is decoded incorrectly, it will affect the SNR of the subsequently decoded symbols.
One way of addressing this issue is by enabling appropriate precoding operation at the transmitter.
In fact A of (31) is known a priori to the transmitter and it is an upper triangular matrix which is
invertible. Due to this reason, it is possible to employ ZF precoding (ZF-P) while ensuring Adl
becomes a diagonal matrix. Since such an approach ensures ISI free received signal, the receiver
can decode the transmitted symbols in symbol by symbol basis. From this discussion, one can
understand that the performance of the ZF-SIC (without error propagation) can be achieved by

30

7
ZF (N=16)
ZF-SIC (N=16)
ZF (N=32)
ZF-SIC (N=32)

6

25

Rate (bps/Hz)

5

log(γ Eqt)

Existing OFDM
Proposed (ZF)
Proposed (ZF-SIC)

4
3

20

15

10

2
5

1

0

0
1

2

3

4

5

6

7

3

8

4

5

6

7

8

Normalized number of antennas (N 0)

Number of multiplexed symbols (S)

Fig. 3: Logarithm of SNR log2 (γEqt ) when Fig. 4: Rate of proposed and existing deL = 8 and γEqt = αγ P r with σ 2 = 10 dB.

signs when L = 8 and γ = 5.3 dB.

employing a ZF-P operation at the transmitter which helps reduce the complexity of receiver.
Complexity: One can understand from the result of this section that the main complexity of the
OFDM approach arises due to the fast Fourier transform (FFT) operation which has complexity
Ns log(Ns ). On the other hand, for the proposed OFTN, the main complexity arises to compute
¯ H H̄
¯ )−1 in (29) of the revised manuscript which has a complexity O(L2.373 ) [41]. In practice
(H̄
since Ns is selected to be much greater than L, the complexity of OFDM and proposed OFTN
approaches are comparable.
Note that for a fixed bandwidth, one may think of increasing the number of sub-carriers by
decreasing the sub-carrier spacing ∆f . In an OFDM system, the symbols which are sent on
carrier frequencies with frequency spacing ∆f , the time-frequency product ∆t∆f may need to
be constant to achieve interference free sub-channels where ∆t is the OFDM duration9 . When
we employ the Nyquist transmission, the time-bandwidth product is roughly one [42]. Now if
we would like to decrease ∆f , the OFDM transmission duration ∆t may need to be increased
to maintain ∆t∆f ≈ 1. Thus, for the fixed OFDM duration ∆t, it is practically not possible to
decrease the sub-carrier spacing beyond some value.

VI. A DAPTIVE P OWER A LLOCATION WITH E XACT AVERAGE R ATES
Up to now we consider the upper-bound rate expressions obtained by existing OFDM and
proposed OFTN approaches by assuming equal power allocation policy. In fact, for the given
9

It is also possible to use the term time-bandwidth instead of time-frequency.

L, one can apply the law of large numbers to show that |hoi |2 and H̄H
I H̄I are almost constant

when N  L (i.e., massive MIMO setup) [34]. In such a case, equal power allocation yields

close to optimal rate (i.e., power adaptation results in negligible average rate gain). For a small to
medium N , however, appropriate power adaptation may help improve the average rates of both
the OFDM and OFTN approaches since the effect of fading is not negligible. Here we consider
the optimal power adaptation strategy which employs the exact average rates as discussed in
[43]. For a fading channel with Gaussian transmitted signal and SISO system, the optimal power
allocation that maximizes the average rate is formulated as


Z
Γ(ψ)
Rav = max
log2 1 +
ψ ρ(ψ)dγ
Pav
Z
s.t
Γ(ψ)ρ(ψ)dψ ≤ Pav

(35)

where Pav is the average power budget, ψ is the received signal SNR, Γ(ψ) is the allocated power,
and ρ(ψ) is the probability density function (PDF) of ψ. The optimal power allocation Γ(ψ) that
maximizes Rav satisfies (see (4)-(6) of [43])

1
1
Γ(ψ)  ψ0 − ψ ,
=
 0,
Pav

ψ ≥ ψ0

where ψ0 is the cut-off SNR which is selected to ensure [44]

Z ∞
1
1
−
ρ(ψ)dψ = 1.
ψ0 ψ
ψ0
Once ψ0 is computed, the optimal average rate is given as
 
Z ∞
ψ
Rav =
log2
ρ(ψ)dψ.
ψ0
ψ0

(36)

ψ < ψ0

(37)

(38)

One can observe from (38) that the optimal average rate depends on ρ(ψ)10 .
Ex
) and
Now if we adopt the steps (35) - (38) to compute the average rates of the OFDM (Rav

Pr
proposed OFTN (Rav
) transmissions, we will have
Ex
Rav
= Rav |ρ(ψ)=ρ(ψ)Ex
10

(39)

Note that to realize the power adaptation policy, the instantaneous received signal SNR ψ may need to be available at the

transmitter which can be obtained from the feedback SNR in frequency division duplex (FDD) system or from the estimated CSI
in time division duplex (TDD) system.

where ρ(ψ)Ex is a Gamma distribution with parameters k = N L and θ =

1
.
σ2

This is due to the

fact that 2σ 2 γ Ex in Theorem 1 is a Chi-square (χ2 ) distribution with 2N L degrees of freedom.
Furthermore, for the proposed OFTN approach with ZF-SIC (no error propagation), we will have
Pr
Rav
=

1
Rav |ρ(ψ)=ρ(ψ)P r
α

where ρ(ψ)P r is a Gamma distribution with parameters k =

(40)
N −L+1
2

and θ =

2α
σ2

when S = L.

VII. OFTN V S FTN, AND E FFECT OF P ULSE S HAPING F ILTER
A. OFTN Vs FTN for SISO Systems
The current paper considers a FTN transmission scheme for wideband channels having multipath components. This will raise a question what is the key difference between the proposed
and existing FTN transmission approaches just for SISO additive white Gaussian noise (AWGN)
channel which is the focus of many existing works [4]. In fact, at the transmitter side, one can still
use the methodology used in Section V to realize FTN transmission (for example with precoding).
Furthermore, by stacking the first Q̃ = L received samples taken every T̃s and after filtering and
doing some manipulations, we can express the filtered signal as
˜l
yl = dl + A0 w̃

(41)

˜ l noise vectors where each of its entry is i.i.d ZMCSCG random variable each with
where w̃
variance σ 2 and A0 is as defined in (24). As we can see (41) contains correlated noise vectors
which will share the same drawbacks as that of the existing FTN transmission schemes. This
shows that in the case of AWGN or flat fading channels, the proposed design does not have any
advantages compared to the existing FTN transmission. This also confirms that the gain achieved
by the proposed transmission does not come just because of rectangular pulse shaping filter11 .
We would like to emphasize here that the current paper focuses on a typical OFTN transmission
and the corresponding data estimation technique for a wideband SIMO system. The technique
used in this paper, however, may not be necessarily optimal. Nevertheless, as will be demonstrated
in the simulation section, the proposed technique achieves significantly better transmission rate
than that of the existing OFDM technique.
11

We would like to point out here that one can also enable the benefits of the other sampling time received signals of the SIMO

channel. However, such an approach will create correlated noise and increase the requirement of ADC speed which will introduce
additional complexity to the receiver. Nevertheless, this approach will definitely help to improve performance, and analyzing such
a case is beyond the scope of the paper is left for future.

B. Effect of Pulse Shaping Filter
As can be seen from the discussions in section IV, the proposed approach employs a rectangular
pulse shaping filter g(.). When g(t) is different from rect(x), (22) will be slightly modified which
consequently yields different yl . In such a case, however, it is difficult to express yl in closed
form. Nevertheless, as the channel coefficients h1 , h2 , · · · hP are not related to g(t), we believe that
employing different g(t) would change the gain achieved by the proposed transmission scheme

only slightly. Furthermore, one can still employ rect g(t) just by utilizing appropriate guard band
like in the OFDM technique.
VIII. S IMULATION R ESULTS
This section provides simulation results, we set B = 1.25 MHz, S = L = P = 8, Td =

L
B

and the carrier frequency is set to 2 GHz. For such settings, we compare the rates obtained
by the proposed and existing OFDM transmission approaches. For the OFDM approach, the
number of sub-carriers (i.e., FFT size) is set to 128 as in LTE system which is used in 1.4 MHz
channel bandwidth transmission [27], [28]. All results are obtained by averaging 10000 channel
realizations, the SNR is defined as γ =

E{|di |2 }
,
σ2

and N = 2N0 where N0 is an integer which is

chosen to attain the desired N .
A. Comparison of Existing OFDM and Proposed Approaches
This subsection compares the upper-bound rates achieved by the existing OFDM and proposed
approaches. Towards this end, we assume that Q = I and R = I. Fig. 4 shows the rates achieved
by the proposed and existing designs for different values of N when γ = 5.3 dB. As we can see
from this figure, the proposed design with ZF-SIC detection approach achieves better performance
than that of the existing OFDM approach approximately when N ≥ 12, and the rate gap between

the proposed and existing approach increases as N increases12 . Furthermore, the ZF-SIC data
detection approach yields better performance than that of the ZF data detection which is expected.
In the next simulation, we examine the effect of γ on the performances of the proposed and
existing OFDM transmissions as shown in Fig. 5. One can observe from this figure that the
proposed ZF-SIC approach achieves better performance than that of the existing approach when
N ≥ 8 at γ = 15.3 dB. Furthermore, the required N such that the proposed approaches (ZF
12

Recall that the performance of ZF-SIC is the same as that of ZF-P.

40

40
Existing OFDM
Proposed (ZF)
Proposed (ZF-SIC)

35

35
30

SNR=15.3dB
25

Rate(bps/Hz)

Rate (bps/Hz)

30

Theory (SNR=5.3dB)
Simulation (SNR=5.3dB)
Theory (SNR=15.3dB)
Simulation (SNR=15.3dB)

20
15

25
20
15

10
5

10

SNR=-4.6dB

0
3

4

5

6

7

Normalized number of antennas (N 0)

8

5
3

4
5
6
7
Normalized number of antennas (N )

8

0

Fig. 5: Effect of SNR on the performances Fig. 6: Comparison of Theoretical and Simof the Existing OFDM and proposed ap- ulation rates of the proposed design with
proaches when L = 8.

ZF-SIC detection approach.

and ZF-SIC) achieve better performance than that of the existing approach increases as the SNR
decreases. This demonstrates that the proposed approach can still be applied for existing LTE
systems where each BS is possibly equipped with few antennas in the high SNR regime.
B. Validation of Theoretical Results
This subsection validates the theoretical results derived in Section V. In this aspect, we demonstrate the theoretical upper-bound rate derived in Theorem 1 for different N0 while setting the
SNR as γ = {5.3, 15.3} dB in Fig. 6. As can be seen from this figure, the results given in

Theorem 1 fit well with the simulation results for all antenna array sizes. In Theorem 1, it is
stated that the proposed design will achieve better data rate than OFDM transmission when the
number of antennas is selected as in (33) which is demonstrated in Fig. 7. For the settings of
this figure, one can obtain N of (33) as N = [16, 8, 8] for γ = [−4.6, 5.3, 15.3] dB. As can be
seen from this figure, the required number of antennas to achieve better performance decreases
as SNR increases which demonstrates the benefit of the proposed design even for small antenna
array sizes (i.e., N = P = 8) at the high SNR regime. Furthermore, the approximated expression
of (33) can be regarded as a good estimation for practical system design.

35

40
Existing OFDM (SNR=-4.6dB)
Proposed ZF-SIC (SNR=-4.6dB)
Existing OFDM (SNR=5.3dB)
Proposed ZF-SIC (SNR=5.3dB)
Existing OFDM (SNR=15.3dB)
Proposed ZF-SIC (SNR=15.3dB)

Rate (bps/Hz)

30
25

Proposed MRC (SNR=5.3dB)
Proposed ZF (SNR=5.3dB)
Proposed MRC (SNR=15.3dB)
Proposed ZF (SNR=15.3dB)

30
25

Rate (bps/Hz)

35

20
15

20
15
10

10

5

5
0

0
3

4

5

6

7

8

Normalized number of antennas (N0)

3

4

5

6

7

8

Normalized number of antennas (N 0)

Fig. 7: Validation of the required number of Fig. 8: Comparison of MRC and ZF beamantennas expression of (33).

formings of the proposed design.

C. Effect of Beamforming Scheme
One can notice from the above results that the ZF-SIC beamformer achieves better performance
than that of ZF. However, when N and L are large, one may think of employing a simple MRC
beamforming approach to reduce complexity. Thus, it is desirable to compare the performance
of ZF and MRC beamforming approaches. In this regard, we have used R = I, Q = I and
γ = {5.3, 15.3} dB. One can observe from Fig. 8 that ZF beamformer achieves significantly better

performance than that of the MRC. Furthermore, increasing the SNR level of the transmission
does not help improve the upper-bound rate of MRC. The potential reason for this is that the MRC
beamformer does not yield a strictly diagonal H̄H
I H̄I which will ultimately create non-vanishing
ISI terms which are dominant at the very high SNR regime. Consequently, increasing the SNR
does not provide meaningful rate improvement.
D. Desired Rate Regions
One can observe from Fig. 4 - Fig. 7 that the proposed OFTN approach achieves better
performance in the region where the average rate is higher than a certain threshold. This raises
a concern whether the performance improvement is in a practically relevant rate region or not.
In general, it is not trivial to provide exact answer for this concern for a general setup. Here, we
try to provide some insight for the setup Q = I, R = I and L = P = Q̃ = S. For this setup,

one can compare the rates of OFDM and OFTN approaches for fixed N or SNR. In fact, the rate

gap between OFTN and OFDM approaches increases with N for fixed SNR. Furthermore, as
the proposed OFTN transmission packs L symbols in LT̃s seconds (i.e., one symbol period with
Nyquist transmission)13 , the actual constellation of these L symbols is dictated by the quantity
1 Pr
R .
L

These observations motivate us to compare

1 Pr
R
L

and

1 Ex
R
L

(termed as normalized rate)

for different SNR and L by setting the minimum N (i.e., N = L) as shown in Fig. 9. One can
observe from this figure that the proposed approach achieves better performance for the region
R ≥ Rth and most importantly Rth decreases with L. For instance, when L = 12, Rth ≈ 1 (i.e.,

BPSK transmission) and when L = 4, Rth ≈ 3 (i.e., 8 PSK transmission)14 . This phenomena
renders the proposed OFTN approach advantageous compared to OFDM as L is related to the

transmission bandwidth. And for fixed SNR and N , it is always possible to choose L optimally
while ensuring that the OFTN approach reaps the maximum benefit compared to OFDM. One
simple, albeit suboptimal, strategy would be to choose L as the maximum value that ensures
the OFTN approach to have a normalized rate at least 1 (i.e., BPSK constellation) which can be
computed by gradually decreasing L starting from L = N 15 .
From this discussion, one can understand that the proposed OFTN approach is flexible and can
be applied for high and low data rate applications.
E. Tightness of the Upper-bound Rates and Effect of Power Adaptation
In this simulation, we examine the tightness of the rate expressions REx and RP r , and the
effect of power adaptation on the performances of OFDM and proposed OFTN approaches. As
stated in Section VI, power adaptation may help improve the average rates of both the OFDM
and OFTN approaches for a small to medium N since the effect of fading is not negligible. For
this reason, we employ L = 8, R = I, Q = I, γ = 5.3 dB and small to medium N , and compare
Ex
Pr
{Rav
, Rav
} (i.e., optimal power adaptation), and {REx , RP r } (i.e., equal power allocation) as

shown in Fig. 10 where we compute the integral terms (37) and (38) numerically. One can
observe from this figure that the rate gap between the proposed approach and that of OFDM
increases as N increases for fixed L. Furthermore, this gap is almost similar with and without
power adaptation. In addition, the upper-bound rate expressions presented in Theorem 1 is tight
13

Note that the ratio

L
3L−2

is just due to the rate inefficiency because of padding 2L − 2 zero samples every 3L − 2 sampling

periods.
14
Note that the region with a normalized rate less than one is not desirable as it does not represent a realistic signal constellation.
15
Recall that a slightly better performance optimization can be achieved by further optimizing S for fixed L.

3.5

25
Existing OFDM (L=4)

Existing OFDM (L=8)

2.5

20

Proposed ZF-SIC (L=8)
Existing OFDM (L=12)

2

Exact rate with optimal power (Existing OFDM)
Upper-bound rate with equal power (Existing OFDM)
Exact rate with optimal power (Proposed ZF-SIC)
Upper-bound rate with equal power (Proposed ZF-SIC)

Proposed ZF-SIC (L=4)

Rate (bps/Hz)

Normalized rate (bps/Hz)

3

8 PSK

Proposed ZF-SIC (L=12)

1.5
1

15

10

0.5
BPSK
0
-24

5
-18

-12

-6

0

6

12

18

24

SNR (dB)

8

16

24

32

40

48

56

64

Number of antennas (N)

Fig. 9: Comparison of the normalized rates Fig. 10: Average rates achieved by OFDM
of OFTN and OFDM approaches for differ- and OFTN approaches with optimal and
ent SNR and L.

equal power allocations when γ = 5.3 dB.

and power adaptation provides almost no gain both for OFDM and proposed OFTN approaches
Ex
Pr
≈ RP r .
since Rav
≈ REx and Rav

F. Effect of Channel Correlations
This simulation examines the effect of temporal and spatial channel correlations on the performances of the proposed design. For the spatial channel correlation matrix, we employ the well
known exponential correlation model as Qi,j = ρ|i−j| , where 0 ≤ ρ < 1. The exponential model

is physically reasonable in a way that the correlation between two transmit antennas decreases
as the distance between them increases [45]. Furthermore, this model is a widely used antenna
correlation model for urban area communications where traffic is usually congested [46]. Fig. 11
shows the performances of the proposed designs for ρ = {0.1, 0.5, 0.95} with γ = 5.3 dB and

R = I. One can observe from this figure that the performance of the proposed design degrades
when ρ increases for all N . However, for fixed ρ, increasing N not only helps increase the upperbound rate but it also help decrease the rate gap between the correlated channel and the i.i.d (i.e.,
ρ = 0). The main reason for this observation is that for fixed ρ, increasing N will increase the
number of eigenvalues of Q. And for this simulation since S = P , the minimum normalized
power gain is the same as the P th maximum normalized eigenvalues of Q which increases with
N . To show this, we plot the normalized eigenvalues of Q for ρ = 0.5 and different values of N
in Fig. 12. As can be seen from this figure

λP
λ1

increases as N increases. This also validates the

1

30
Proposed ZF-SIC (ρ=0)
Proposed ZF-SIC (ρ=0.5)
Proposed ZF-SIC (ρ=0.95)

25

N=8
N=32
N=64

0.9

20

0.7
λk/λ1

Rate (bps/Hz)

0.8

15

10

0.6
0.5
0.4
0.3

5

0.2
0
3

4

5

6

7

0.1
1

8

Normalized number of antennas (N 0)

8

15

22
29
36
43
Eigenvalue index (k)

50

57

64

Fig. 11: Performance of the proposed ZF- Fig. 12: Normalized Eigenvalue distribution
SIC detection for different ρ.

λi
λ1

(decreasing order) when ρ = 0.5.

result of Lemma 1 which states that the performance of the proposed design does not depend on
Q when N  L and rank(Q) ≈ P .

A number of approaches have been proposed to model the temporal channel correlation matrix

R. One most commonly employed model is the COST 207 for urban environment where R is a
diagonal matrix with its elements decay exponentially [47], [48]. Under this model, Ri,i can be
chosen as (see (5), (7) and (19) of [47])
Ri,i = κ

Z

τi

e−Bx dx

(42)

τ(i−1)

where κ is a normalization factor which depends on the carrier frequency fc and is chosen to
ensure R0,0 = 1 and B is the transmission bandwidth. With this model, the diagonal elements of
R becomes [1, 0.3679, 0.1353, 0.0498, 0.0183, 0.0067, 0.0025, 0.0009]. As we can see, the diagonal
elements of R decrease with i. This will ultimately decrease the effective rank of R. For this
reason, one can expect reduced data rate when R is designed as in (42) compared to that of
R = I. The simulation results demonstrating this fact has been omitted for conciseness.
From the result of this section, we can understand that the performance of the proposed OFTN
transmission approach varies from one scenario to another. Furthermore, one can improve the
performance of the proposed design by selecting S adaptively based on the structures of R and
Q. Nevertheless, how to design optimal S for different scenario is left for future research

IX. O PEN R ESEARCH T OPICS
The paper focuses on the OFTN transmission for SIMO system by assuming typical multipath
channel models. In this section, we first discuss some of the open research topics for the SIMO
system. Then, we discuss how one can utilize the proposed OFTN transmission for other single
user and multiuser systems, and the associated research problems.
A. Single User SIMO Systems
In a typical OFDM system, the achievable rate can be computed by assuming equal power allocation, or by optimizing the powers allocated to all sub-carriers. In fact, applying a suitable power
allocation strategy will definitely help improve the achievable rate. However, this improvement
is achieved at the expense of additional computational cost. In the current paper, we compare the
performance of the proposed transmission with that of the OFDM with equal power allocation
just for simplicity. The comparison of the OFDM transmission and the proposed approach under
optimized power allocation is an interesting open research problem.
In general, the covariance matrices Q and R of (25) depend on different factors including the
deployment environment, carrier frequency and number of scatterers [49]. Furthermore, there are
different multipath channel covariance matrix models for practically relevant wireless environments like Extended Pedestrian A (EPA), Extended Vehicular A (EVA), or Extended Typical Urban
(ETU). In the current paper, however, we have provided numerical results for the most commonly
adopted channel covariance matrices in the literature such as R = I (i.e., multipath coefficients
have equal average power), to simplify the analysis and to get an insight about the performance of
the proposed approach. Furthermore, the current paper assumes that the channel state information
(CSI) is known to the receiver perfectly. In practice, however, the CSI may need to be estimated,
and often known imperfectly. Thus, the performance comparison of the proposed approach and
OFDM with perfect and imperfect channel model parameters of the above environments is still
an open problem.
The current paper assumes that the receive antennas are collocated (i.e., single cell setup).
However, it may be the case that the antennas may be placed at different locations like in the
distributed (network) MIMO system [50]. In general, the proposed approach can still be applied
for the latter system. And in such a system, the elements of R and Q will be different from that

of the collocated case. Nevertheless, comparing the proposed approach with that of OFDM for
network MIMO system case is an interesting questions which needs further research.
B. Single User MISO and MIMO Systems
In this subsection, we discuss the potential extensions and open research issues of the proposed
data transmission and detection approach for MISO and MIMO channels. In fact, if the transmitter
is equipped with multiple antennas, multiplexing gain can be achieved. In this regard, one may
think of two possibilities: the first possibility is to simply transmit independent information from
each antenna. The second possibility is to utilize appropriate precoding operation. The former
possibility can be realized like in the aforementioned sections which consequently help scale
the data rate with the number of transmitter antennas M in a favorable propagation condition16 .
However, the receiver requires at least N = M S antennas to reliably estimate the transmitted
signal with linear receivers. Although, the latter possibility is well understood from existing
communication system, the beamforming design with the proposed OFTN transmission is still
not understood and it is left for future work.
The remaining issue is whether the proposed approach can be extended for the MISO channel?
As stated previously, when the transmitter has multiple antennas, each antenna can transmit independent information. However, for each transmit antenna, S receive antennas are required. This
leads to a fundamental question on how to realize the proposed OFTN transmission and detection
(estimation) approach for MISO systems. From the existing OFDM transmission approach, it is
well known that the rates achieved by interchanging the roles of the transmitter and receiver are
the same (i.e., there is a so called uplink-downlink rate duality (see for example [51]–[53])).
However, showing the rate duality of MISO and SIMO channels under the proposed OFTN
transmission and data detection is still not understood well and it is left for future research.
C. Multiuser Systems
The proposed data transmission and detection approach is designed for a single user multiantenna system with multipath components. In fact, multipath components exist almost in all
wireless wideband channels which justifies that the proposed transmission approach can potentially
16

Here a favorable propagation condition means that the spatial channel covariance matrix between any transmitter and all

receiver antennas is close to I (or well conditioned).

be applied for different system setups. For example, it is applicable to a wideband multiuser
(massive) MIMO system where the BS is equipped with several antennas whereas, each UE is
equipped with single antennas. For such a system, each UE can apply the proposed transmission
approach. Consequently, the rate achieved by each UE increases with P . Nevertheless, examining
the performance of the proposed design for different system setups is one of the interesting
research directions. The performance can be studied by taking into account the effects of pulse
shaping filter g(t) and error propagation in ZF-SIC receivers.
X. C ONCLUSIONS
In this paper, we propose a novel and simple OFTN transmission and detection approach
for a wideband SIMO system with L multipath components. The proposed OFTN transmission
exploits the multiplexing gain obtained by the inherent characteristics of multipath components of
wideband channels. By doing so, the proposed design achieves significantly higher transmission
rate than the conventional Nyquist transmission approach (i.e., OFDM). The superiority of the
proposed approach has been shown theoretically and demonstrated through extensive computer
simulations. It is also shown that the capacity gap between the proposed approach and that of
the OFDM increases as the number of receiver antennas increases for fixed L which implies
that the new approach may have practical significance for large antenna array receiver systems.
The proposed method can be extended to (massive) MIMO systems. Having said this, however,
extending the proposed approach for a MISO system while achieving the same capacity as that
of the SIMO system is a non trivial fundamental challenge.
A PPENDIX A
P ROOF OF Theorem 1
For the existing OFDM transmission, by utilizing (7) and (11), the channel gain of the sth
sub-carrier of each receive antenna, gs of (16) can be expressed as
gs = fsH R1/2 h̄ = h̄T (R1/2 )T fs∗

(43)

where fs is an L sized truncated DFT matrix with overall size Ns [18], [29]. Using (25) and after
some manipulations, E{|g̃s |2 } becomes
E{|g̃s |2 } = E{|Q1/2 H̄I (R1/2 )T fs∗ |2 } = N fsT Rfs∗ ,

⇒ γEx =

E{|g̃s |2 }
N T ∗
=
f Rfs
σ2
σ2 s

(44)

where the second equality uses tr{Q} = N for any Q. When R = I, we will get γEx of (32).
For the proposed approach, when the ZF beamforming with SIC is used, we will have
1

γi =

−1
σ 2 E{([H̄H
I H̄I ] )i,i }

.

(45)

To compute the expected value of an inverse function, we consider the following lemma.
Lemma B.1: Let A = ZH Z be a non-singular Hermitian matrix of size K × K and B = A−1 ,

where Z ∈ C N ×K and K ≤ N . We partition A and B as




a11 aH
b
b
12
21
 , B =  11

A=
a21 A22
b21 B22

(46)

where a11 (b11 ) is a scalar value, and the remaining terms are appropriate dimension vectors or
matrices. If a11 6= 0 and A22 is non-singular, we can express b11 as
1
−1
= a11 − aH
21 A22 a21 .
b11

(47)

And, if each element of Z is taken from i.i.d ZMCSCG random variable with variance 1, then
1
∼ χ2N −K+1 .
b11

(48)

Proof: Equation (47) can be proved by applying the well known Schur Complement theorem.
The detailed derivation can also be found from Theorem A5.2 of [54].
−1
To prove (48) we note that both a11 and aH
21 A22 a21 are strictly non-negative real values. And
−1
when a11 − aH
21 A22 a21 is a non-negative real valued term, by applying Theorem 3.2.10 of [54],

the probability density function of

1
b11

can be expressed as W1 (N − K + 1, 1), where W(.) (., .)

denotes a real valued Wishart distribution. It follows

1
∼ W1 (N − K + 1, 1) ∼ χ2N −K+1
b11

(49)

where the second distribution is due to the fact that W1 (N − K + 1, 1) has the same distribution

as that of Chi-square (χ2 ) distribution with N − K + 1 degrees of freedom (see Corollary 3.2.2

of [54]). As expected when N = K,
Z = H̄I of (45), we will get γi =

1
b11

is a χ2 distribution with 1 degree of freedom. By setting

N −L+1
.
σ2

As discussed in Section IV, the proposed design transmits L symbols in 3L − 2 sampling

periods. On the other hand, the design splits the overall bandwidth B into L sub-bands. For these
)( L1 ) to ensure that the proposed design uses
reasons, we may need to incorporate α = ( 3L−2
L

the same transmit power as that of OFDM, while simultaneously transmitting

1
α

more symbols

compared to OFDM. Upon doing so, the proposed design can achieve RP r of (32).
Next we show (33) of Theorem 1. If γ0 =

N
σ2

 1 and N  L = 3, we will have

L2
(log2 (α) + log2 (γ0 ))
REx = log2 (γ0 ) + log2 (L), RP r =
3L − 2


L2
L2
RP r − REx ≥ 0 ⇒
− 1 log2 (γ0 ) ≥ log2 (L) −
log2 (α)
3L − 2
3L − 2


1
1
2
2
2 ( L2 −3L+2 ((3L−2) log2 (L)−L log2 (α)))
2 [ L2 −3L+2 ((3L−2) log2 (L)−L log2 (α))]
⇒N ≥σ e
= max L, σ e
(50)
where the last inequality follows from the fact that N ≥ L.
A PPENDIX B
P ROOF OF L EMMA 1
For the existing OFDM transmission, one can utilize the result of Appendix A. For the proposed
approach, when N → ∞, we obtain (Q1/2 H̄I )H (Q1/2 H̄I ) ≈ tr(Q)I = N I almost surely [19]. It

follows that

¯ A)H (H̄
¯ A))−1 ≈ N −1 (AH RA)−1 = N −1 E, (H̄
¯ H H̄
¯ )−1 ≈ N −1 C
((H̄

(51)

where C = R−1 and E = (AH RA)−1 which are deterministic matrices. Thus, for the ZF
detection, we will achieve
γi =

1
¯ A)H (H̄
¯ A)]−1 ) }
σ 2 E{([(H̄
i,i

≈

N
.
σ 2 E(i,i)

(52)

And for the ZF-SIC detection approach, we will have
γi =

1
¯H

¯ ]−1 ) }
σ 2 E{([H̄ H̄
S,S

≈

N
σ 2 C(i,i)

.

(53)

R EFERENCES
[1] E. Hossain, M. Rasti, H. Tabassum, and A. Abdelnasser, “Evolution towards 5G multi-tier cellular wireless networks: An
interference management perspective,” IEEE Wireless Commun. Mag., pp. 118 – 127, Jun. 2014.
[2] A. Osseiran, F. Boccardi, V. Braun, K. Kusume, P. Marsch, M. Maternia, O. Queseth, M. Schellmann, H. Schotten, H. Taoka,
H. Tullberg, M. A. Uusitalo, B. Timus, and M. Fallgren, “Scenarios for 5G mobile and wireless communications: The vision
of the METIS project,” IEEE Commun. Mag., pp. 26 – 35, May 2014.
[3] S. Sugiura, “Frequency-domain equalization of faster-than-Nyquist signaling,” IEEE Wireless Commun. Letters, vol. 2, no.
5, pp. 555 – 558, Oct. 2013.

[4] J. B. Anderson, F. Rusek, and V. Owall, “Faster-than-Nyquist signaling,” IEEE Proceedings., pp. 1817 – 1830, Aug. 2013.
[5] T. Bogale and L. Le, “Massive MIMO and mmWave for 5G wireless HetNet: Potential benefits and challenges,” IEEE
Vehic. Techno. Mag., vol. 11, no. 1, pp. 64 – 75, Mar. 2016.
[6] 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.
[7] J. E. Mazo., “Faster-than-Nyquist signaling,” Bell Syst. Technical J., vol. 54, no. 8, pp. 1451 – 1462, Oct. 1975.
[8] A. D. Liveris and C. N. Georghiades, “Exploiting faster-than-Nyquist signaling,” IEEE Trans. Commun., vol. 51, no. 9, pp.
1502 – 1511, Sep. 2003.
[9] M. E. Hefnawy and H. Taoka, “Overview of faster-than-Nyquist for future mobile communication systems,” in Proc. IEEE
Veh. Technol. Conf (VTC-Spring)., 2013, pp. 1 – 5.
[10] M. E. Hefnawy and G. Kramer, “Impact of spectrum sharing on the efficiency of faster-than-Nyquist signaling,” in Proc.
IEEE Wireless Commun. and Networking. Conf., 2014, pp. 648 – 653.
[11] M-S. Baek, N-H. Hur, and H. Lim, “Novel interference cancellation technique based on matrix computation for FTN
communication system,” in Proc. IEEE Military Commun. Conf., 2014, pp. 830 – 834.
[12] F. Schaich and T. Wild, “A reduced complexity receiver for multi-carrier faster-than-Nyquist signaling,” in Proc. IEEE
Global Telecommun. Workshop, 2013, pp. 235 – 240.
[13] A. Prlja and J. B. Anderson, “Reduced-complexity receivers for strongly narrowband intersymbol interference introduced
by faster-than-Nyquist signaling,” IEEE Trans. Commun., vol. 60, no. 9, pp. 2591 – 2601, Sep. 2012.
[14] F. Rusek, “On the existence of the Mazo-limit on MIMO channels,” IEEE Trans. Wireless Commun., vol. 8, no. 3, pp.
1118 – 1121, Mar. 2009.
[15] F. Rusek, “A first encounter with faster-than-Nyquist signaling on the MIMO channel,” in Proc. IEEE Wireless Commun.
and Networking. Conf., 2007, pp. 1094 – 1098.
[16] L. Zheng and D. N. C. Tse, “Diversity and multiplexing: A fundamental tradeoff in multiple-antenna channels,” Trans.
Commun., vol. 49, no. 5, pp. 1073 – 1096, 2003.
[17] D. Tse and P. Viswanath, Fundamentals of Wireless Communication, Cambridge University Press, 2005.
[18] D. P. Palomar, A unified framework for communications through MIMO channels, Ph.D. thesis, Technical University of
Catalonia (UPC), Barcelona, Spain, 2003.
[19] H. Q. Ngo, E. G. Larsson, and T. L. Marzetta, “Energy and spectral efficiency of very large multiuser MIMO systems,”
IEEE Trans. Commun., vol. 61, no. 4, pp. 1436 – 1449, Apr. 2013.
[20] J. Hoydis, C. Hoek, T. Wild, and S. T. Brink, “Channel measurements for large antenna arrays,” in Proc. IEEE International
Symposium Wireless Communication Systems (ISWCS), 2012, pp. 811 – 815.
[21] X. Gao, F. Tufvesson, O. Edfors, and F. Rusek, “Measured propagation characteristics for very-large MIMO at 2.6 GHz,”
in Proc. Asilomar Conf. Sign. Syst. and Comp., 2012, pp. 295 – 299.
[22] X. Gao, O. Edfors, F. Rusek, and F. Tufvesson, “Massive MIMO performance evaluation based on measured propagation
data,” IEEE Trans. Wireless Commun., vol. 14, no. 7, pp. 3899 – 3911, Jul. 2015.
[23] X. Gao, O. Edfors, F. Rusek, and F. Tufvesson, “Linear pre-coding performance in measured very-large MIMO channels,”
in Proc. IEEE Veh. Technol. Conf (VTC-Fall)., 2011, pp. 1 – 5.
[24] A. Manolakos, M. Chowdhury, and A. Goldsmith, “Energy-based modulation for noncoherent massive SIMO systems,” Jul.
2015, http://arxiv.org/abs/1507.04978.
[25] “Jensen’s inequality,” www.econ.hku.hk/˜wsuen/teaching/micro/jensen.pdf.

[26] T. E. Bogale, Transceiver design for single-cell and multi-cell downlink multiuser MIMO systems, Ph.D. thesis, University
Catholique de Louvain (UCL), Louvain la neuve, Belgium, 2013.
[27] H. Zarrinkoub, Understanding LTE with MATLAB: From Mathematical Modeling to Simulation and Prototyping, John
W&S, Ltd, 2014.
[28] J. Zyren, “Overview of the 3GPP long term evolution physical layer,” in Freescale semiconductor, 2007, pp. 1 – 24.
[29] 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.
[30] A. Goldsmith, Wireless Communications, Cambridge University Press, New York, NY, USA, 2005.
[31] Y. Liu, Z. Tan, H. Hu, L. J. Cimini, and G. Y. Li, “Channel estimation for OFDM,” IEEE Commun. Surv. Tutor., vol. 16,
no. 4, pp. 1891 – 1908, 2014.
[32] 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 (E-UTRA).
[33] ETSI TR 125 943, “Universal mobile telecommunications system (UMTS); deployment aspects,” 2010, 3GPP TR 25.943
version 9.0.0 Release 9.
[34] 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.
[35] J. G. Andrews, A. Ghosh, and R. Muhamed, Fundamentals of WiMAX: Understanding Broadband Wireless Networking,
Pearson Edu., 2007.
[36] S.M. Alamouti, “A simple transmit diversity technique for wireless communications,” IEEE J. Select. Areas in Commun.,
vol. 16, no. 8, pp. 1451 – 1458, 1998.
[37] F. Rusek and J. B. Anderson, “On information rates for faster than Nyquist signaling,” in Proc. IEEE Global Telecommun.
Conf., 2006, pp. 1 – 5.
[38] V. Venkateswaran and A-J. V. Veen, “Analog beamforming in MIMO communications with phase shift networks and online
channel estimation,” IEEE Trans. Signal Process., vol. 58, no. 8, pp. 4131 – 4143, Aug. 2010.
[39] X. Zhang, A. F. Molisch, and S-Y. Kung, “Variable-phase-shift-based RF-Baseband codesign for MIMO antenna selection,”
IEEE Trans. Signal Process., vol. 53, no. 11, pp. 4091 – 4103, Nov. 2005.
[40] T. E. Bogale, L. B. Le, A. Haghighat, and L. Vandendorpe, “On the number of RF chains and phase shifters, and scheduling
design with hybrid analog-digital beamforming,” IEEE Trans. Wireless Commun., vol. 15, no. 5, pp. 3311–3326, May 2016.
[41] Francois Le Gall, “Powers of tensors and fast matrix multiplication,” in Proceedings of the 39th International Symposium
on Symbolic and Algebraic Computation, New York, NY, USA, 2014, ISSAC, pp. 296–303, ACM.
[42] T. Moazzeni, “On the compactness of ofdm and hermite signals,” IEEE Commun. Letters, vol. 20, no. 7, pp. 1313–1316,
July 2016.
[43] A. J. Goldsmith and S-G. Chua, “Variable-rate variable-power MQAM for fading channels,” IEEE Trans. Commun., vol.
45, no. 10, pp. 1218 – 1230, Oct. 1997.
[44] A. J. Goldsmith and P. P. Varaiya, “Capacity of fading channels with channel side information,” IEEE Trans. Info. Theory,
vol. 43, no. 6, pp. 1986 – 1992, Nov. 1997.
[45] S. L. Loyka, “Channel capacity of MIMO architecture using the exponential correlation matrix,” IEEE Commun. Letters,
vol. 5, no. 9, pp. 369 – 371, 2001.
[46] M. Ding, Multiple-input multiple-output wireless system designs with imperfect channel knowledge, Ph.D. thesis, Queens
University Kingston, Ontario, Canada, 2008.

[47] M. Patzold, A. Szczepanski, and N. Youssef, “Methods for modeling of specified and measured multipath power-delay
profiles,” IEEE Trans. Veh. Technol., vol. 51, no. 5, pp. 978 – 988, Sep. 2002.
[48] D. R. Morgan, “Analysis and realization of an exponentially-decaying impulse response model for frequency-selective fading
channels,” IEEE Signal Process. Letters, , no. 15, pp. 441 – 444, Sep. 2008.
[49] T. E. Bogale, L. B. Le, and X. Wang, “Hybrid analog-digital channel estimation and beamforming: Training-throughput
tradeoff,” IEEE Trans. Commun., Dec. 2015.
[50] T. E. Bogale and L. Vandendorpe, “Weighted sum rate optimization for downlink multiuser MIMO coordinated base station
systems: Centralized and distributed algorithms,” IEEE Trans. Signal Process., vol. 60, no. 4, pp. 1876 – 1889, Dec. 2011.
[51] S. Vishwanath, N. Jindal, and A. Goldsmith, “Duality, achievable rates, and sum-rate capacity of Gaussian MIMO broadcast
channels,” IEEE Trans. Info. Theory, vol. 49, no. 10, pp. 2658 – 2668, Oct. 2003.
[52] P. Viswanath and D. Tse, “Sum capacity of the vector Gaussian broadcast channel and uplink-downlink duality,” IEEE
Trans. Info. Theory, vol. 49, no. 8, pp. 1912 – 1921, Aug. 2003.
[53] T. E. Bogale and L. Vandendorpe, “Linear transceiver design for downlink multiuser MIMO systems: Downlink-interference
duality approach,” IEEE Trans. Signal Process., vol. 61, no. 19, pp. 4686 – 4700, Oct. 2013.
[54] R. J. Muirhead, Aspects of Multivariate Statistical Theory, John Wiley and Sons Inc., New Jersey, 1982.

