Poveda Poveda Hector 2011
Poveda Poveda Hector 2011
THÈSE
présentée
pour obtenir le titre de
Notations xv
Introduction 1
i
CONTENTS
ii
CONTENTS
Appendices
iii
CONTENTS
References 161
iv
Acronyms and Abbreviations
v
Acronyms and Abbreviations
EM Expectation-Maximization
ENC low-data rate ENCoder
EQU channel EQUalizer
ESE Elementary Signal Estimator
FCC Federal Communications Commission
FDMA Frequency Division Multiple Access
FFT Fast Fourier Transform
GPRS General Packet Radio Service
GRV Gaussian Random Variable
GSM Global System for Mobile Communications
HSDPA High-Speed Downlink Packet Access
HSPA+ High-Speed Packet Access Plus
HSUPA High-Speed Uplink Packet Access
IBI Inter Block Interference
IC Interference Cancellation
ICI Inter Carrier Interference
IDFT Inverse Discrete Fourier Transform
IDMA Interleaved Division Multiple Access
IEEE Institute of Electrical and Electronics Engineers
IEKF Iterative Extended Kalman Filter
IFFT Inverse Fast Fourier Transform
ISI Inter Symbol Interference
ITU International Telecommunication Union
KF Kalman Filter
LLR Logarithm Likelihood Ratio
LMS Least Mean Squares
LO Local Oscillator
LS Least-Squares
LTE Long Term Evolution
LTE-A LTE Advanced
M-PSK M-ary Phase-Shift Keying
MAI Multiple Access Interference
MAP Maximum A Posteriori
MC-CDMA Multi Carrier CDMA
MC-DS-CDMA Multi Carrier Direct Sequence CDMA
MIMO Multiple-Input Multiple-Output
ML Maximum Likelihood
MMSE Minimum Mean Square Error
MMSE-SD Minimum Mean Square Error Successive Detection
MRC Maximum Ratio Combining
NAMTS Nippon Automatic Mobile Telephone System
vi
Acronyms and Abbreviations
vii
List of Figures
ix
LIST OF FIGURES
x
LIST OF FIGURES
xi
List of Tables
xiii
Notations
xv
Notations
xvi
Introduction
For more than 25 years, a great deal of interest has been paid to mobile communi-
cation systems and the design of new schemes to transmit and receive information.
Engineers and researchers have worked much and exchanged on that topic.
Several PhD in the signal group at the UMR CNRS 5218 IMS have dealt with
channel estimation and symbol detection in multicarrier systems involving sev-
eral users. In [Jamo 07b], the channel was assumed to be an autoregressive (AR)
process. Jamoos suggested estimating both the channel and its AR parameters
by using training sequences and optimal filtering. To avoid an approach dedi-
cated to the non-linear estimation issue, mutually interactive optimal filters were
used. One filter aimed at estimating the channel, whereas the other updated the
estimation of the AR parameters. H∞ and Kalman filtering were tested and the
resulting approaches were studied in a multicarrier direct-sequence code division
multiple access system (MC-DS-CDMA). In [Grol 07], the channel was modeled
by a sinusoidal stochastic process or a low-pass filter version of an AR process.
Then, both the symbols and the Rayleigh fading channel were jointly estimated in
MC-DS-CDMA systems by using a Rao-Blackwellized particle filter cross-coupled
with a Kalman filter.
Despite the great success of the CDMA as a multiple access technique, the cur-
rent trend is to use the orthogonal frequency division multiple access (OFDMA)
and the orthogonal frequency division multiplexing-interleaved division multi-
ple access (OFDM-IDMA). In that case, the input data stream is split into a
number of streams that are transmitted in parallel over a large number of orthog-
onal subcarriers. The frequency-selective fading over the entire bandwidth of the
transmitted signal has hence the advantage of being converted in frequency-flat
fading over each subcarrier. It should be noted that these schemes are partic-
ularly well adapted to mobile wireless communication to provide high-data rate
services, e.g. [3GPP 09], [IEEE 06].
However, there are two unwanted phenomena:
1/ the channel remains unknown and hence has to be estimated.
1
Introduction
2
Introduction
• the minimum mean square error successive detector (MMSE-SD) [Hou 08]
to increase the data rate for an OFDMA system [Pove 09a], [Pove 09b],
[Pove 11c].
• a statistical test, e.g. the binary hypothesis test (BHT) or the cumulative
sum (CUSUM) test, to detect the beginning and the end of the disturbances
induced by another system [Pove 11b], [Pove 11d]. This happens with cog-
nitive radio (CR) systems where the spectrum that is not used by licensed
(primary) users must be detected in order to be shared with unlicensed
(secondary) users1 .
• a modified version of the iterative IDMA receiver [Ping 02a] to design a new
receiver scheme for an OFDM-IDMA system [Pove 12].
• a space time block code (STBC) and a new variant of the iterative IDMA
receiver [Pove 11a].
3
Introduction
In chapter 2, we propose and compare several CFO and channel estimation meth-
ods based on optimal filtering techniques [Pove 10]. As the relation between the
CFO and the received signal is non-linear, a linearized version of the Least-
Squares (LS) approach can be first used to estimate the CFO. More generally, we
suggest using an extended Kalman filter (EKF). The EKF consists in analytically
propagating the estimation through the system dynamics, by means of a first-
order Taylor expansion of the functions defining the state-space representation
of the system. However, as the approximation may not be sufficient to describe
the non-linearity, the EKF may sometimes diverge. To solve this problem, a
second-order linearization can be considered and leads to the second-order EKF
(SOEKF) [Bar 01]. Another solution is to use the iterative extended Kalman
filter (IEKF). In the IEKF, the measurement model is linearized around the up-
dated state vector, instead of the predicted state vector. Then, the process is
iterated until the state vector estimate does not change much.
As the above approaches require the computation of the Jacobians and the Hes-
sians matrices for the first-order and the second-order linearizations respectively,
we also look at the sigma-point Kalman filter (SPKF) [VdMe 04a], namely the
unscented Kalman Filter (UKF) and the central difference Kalman filter (CDKF).
In that case, the state distribution is approximated by a Gaussian distribution,
which is characterized by the so-called sigma points. Then, the sigma points are
propagated through the non-linear system. A weighted combination of the re-
sulting values makes it possible to estimate the mean and the covariance matrix
of the transformed random variable. The UKF is based on the unscented trans-
formation whereas the CDKF is based on the second-order Sterling polynomial
interpolation formula.
Nevertheless, Kalman algorithms require an accurate system model as well as
perfect knowledge of the noise statistics. Therefore, we also propose to evaluate
4
Introduction
the performance of the H∞ filtering, which has been very popular in the field of
control and is more and more used in signal processing (see for instance [Shen 99],
[Laba 05] and [Gire 09]). This approach is designed to be robust against uncer-
tainties. In addition, no Gaussian assumption on the additive noise and the
model noise in the state-space representation of the system is required. Here, we
analyze the relevance of the “extended H∞ filter”, based on Taylor expansion like
the EKF and the “unscented H∞ filter” in which an unscented transformation is
embedded into the “extended H∞ filter” [Li 10].
Then, the above techniques are combined with different processing blocks to de-
sign new receiver schemes:
Firstly, the optimal filtering techniques are combined with the so-called
MMSE-SD [Hou 08] in an iterative way to obtain a CFO estimator that does
not require pilot subcarriers [Pove 09a], [Pove 09b], [Pove 11c].
Secondly, we use the SPKF in a dynamic bandwidth allocation context. The
SPKF is combined with a statistical test, either the BHT or the CUSUM test.
This approach allows the CFOs and the channels to be jointly estimated and the
beginning and the end of a narrow-band interference (NBI) produced by a CR
system to be detected [Pove 11b], [Pove 11d].
5
Introduction
In each case, simulation results confirm the efficiency of the proposed architec-
tures in terms of CFO estimation, channel estimation, computational complexity
and BER performance.
• In the first one, we recall the main results about the Kalman filter in the
linear case.
• In the fourth appendix, we introduce the EKF, the SOEKF and the IEKF.
• Finally, the last appendix deals with the H∞ filtering for non-linear cases.
The “extended H∞ filter” and the “unscented H∞ filter” are presented.
6
Chapter 1
Multicarrier High-Data Rate Mobile
Wireless Systems
Contents
1.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . 8
1.2 Wireless Communication Systems . . . . . . . . . . . 8
1.2.1 Conventional Systems . . . . . . . . . . . . . . . . . . 9
1.2.2 Cognitive Radio Systems . . . . . . . . . . . . . . . . 9
1.3 Propagation Channel Model . . . . . . . . . . . . . . . 11
1.3.1 Rician and Rayleigh Fading Channels . . . . . . . . . 13
1.3.2 Frequency-Selective and Time-Varying Fading Channels 15
1.4 Multiple Access Techniques for Mobile Systems . . . 19
1.4.1 About Multicarrier Multiple Access Techniques . . . . 22
1.4.2 OFDMA System . . . . . . . . . . . . . . . . . . . . . 27
1.4.3 OFDM-IDMA System . . . . . . . . . . . . . . . . . . 36
1.5 Conclusions . . . . . . . . . . . . . . . . . . . . . . . . . 43
7
1. Multicarrier High-Data Rate Mobile Wireless Systems
1.1 Introduction
The purpose of this chapter is first to give general information about wireless
communication systems. Then, we present the unwanted phenomena that may
appear when transmitting wireless information. More particularly, the channel
plays a key role. Its statistical properties may also vary and depend on the exis-
tence of a line-of-sight (LoS) between the transmitter and the receiver. Notions of
frequency-selective fading or flat-fading channels are recalled. In the second part
of this chapter, we introduce the various access schemes. After recalling what
FDMA, time division multiple access (TDMA) and code division multiple access
(CDMA) are all about, OFDMA and OFDM-IDMA are presented.
8
1.2 Wireless Communication Systems
9
1. Multicarrier High-Data Rate Mobile Wireless Systems
the same frequency spectrum. However, secondary users should transmit the sig-
nal with a low power that does not exceed an interference threshold [Pali 10]. On
the other hand, in overlay systems the secondary users have to find out the bands
that are not used by the primary users (PUs). Then, secondary users can utilize
these unused portions of the spectrum. For this purpose, the secondary users
need information about the spectrum allocation of the PUs by regularly perform-
ing spectrum sensing techniques. Various spectrum sensing techniques have been
proposed [Arsl 07]. Some of them aim at identifying the characteristics of the
transmission whereas others distinguish the signal type. The most common are
the matched filter, the cyclostationary feature detection and the energy detection
[Arsl 07]. For more details about CR systems, the reader is referred to [Pali 10].
Understanding that CR systems have the potential to exploit the underutilized
conventional system frequency bands by spectrum sharing, the Federal Commu-
nications Commission (FCC) released a second report in which unlicensed devices
are allowed to operate in the unused portions of the TV-spectrum in 2008.
CR systems will be considered in this PhD dissertation, in chapter 2, section 2.5.
In the following section, we present one of the more challenging tasks when de-
signing a wireless system, taking into account the influence of the propagation
channel.
10
1.3 Propagation Channel Model
• reflection: it arises when the plane waves are incident upon a surface with
dimensions that are very large compared to the wavelength.
• scattering: it happens when the plane waves are incident upon an object,
the dimensions of which are of the order of a wavelength or less, and causes
the energy to be redirected in many directions.
11
1. Multicarrier High-Data Rate Mobile Wireless Systems
the other two phenomena are statistically described. In the following, let us focus
our attention on the multipath effect.
The received signal is a superposition of several delayed and attenuated copies
of the transmitted signal. By considering a stationary propagation, the channel
propagation model h(τ ) at time τ can be mathematically expressed as:
L−1
X
h(τ ) = A(l)δ(τ − τl ) (1.1)
l=0
where L is defined as the number of paths, A(l) and τl = lTsamp are the amplitude
and the time delay associated to the lth path, respectively. In addition, Tsamp is
the sampling time.
Thus, the coefficients of the discrete-time channel impulse response (CIR) can be
stored in the following vector:
Figure 1.3 shows an example of a CIR, where τmax represents the maximum chan-
nel delay spread in seconds.
h(τ )
0 1 2 3 4 5 6 7 8
Tsamp τ
τmax
12
1.3 Propagation Channel Model
The fading that modifies the transmitted signal depends on the characteristics
of the channel and the nature of the transmitted signal. Moreover, the statis-
tical model of the channel is slightly different if the transmitter is in LoS with
the receiver. The multipath fading channel can be hence classified by looking at
its probability density distribution and its frequency response in the transmitted
signal bandwidth.
These classifications of fading are investigated in the next subsections.
13
1. Multicarrier High-Data Rate Mobile Wireless Systems
The received signal distribution can be rewritten by using the Rice factor defined
by:
Krice = µ2r /2σr2 (1.4)
Krice = 0
Krice = 3
1.8 Krice = 7
Krice = 10
1.6
1.4
1.2
Pdf (x)
0.8
0.6
0.4
0.2
0
0 0.5 1 1.5 2 2.5 3 3.5 4
x
14
1.3 Propagation Channel Model
plane waves arrive from many directions with equal probability without a direct
LoS component. Then, each coefficient of the CIR follows a Rayleigh distribution
[Jake 74]. Therefore, given (1.3) and under these assumptions, the magnitude of
the received signal has the following Rayleigh distribution:
2x − Ωx2
Pdf (x) = e r x≥0 (1.7)
Ωr
where Ωr = E{|αr (t)|2 } = 2σr2 is the average envelope power [Stub 02].
Figure 1.5 shows examples of Rayleigh distribution for different values of Ωr .
In the following subsection, we present the frequency-selective and the time-
varying channels.
Ωr = 0.5
Ωr =1
0.8 Ωr =2
Ωr =4
0.7
0.6
0.5
Pdf (x)
0.4
0.3
0.2
0.1
0
0 1 2 3 4 5 6 7
x
15
1. Multicarrier High-Data Rate Mobile Wireless Systems
i.e.
τmax << Ts (1.10)
At the receiver, without the noise contribution, the signal is an attenuated copy
of the transmitted signal:
r(n) = hs(n) (1.11)
where s(n) is the received symbol, h is the CIR coefficient and r(n) is the received
symbol.
On the other hand, if the multipaths are spread outside the symbol duration, i.e.
the maximum channel delay spread is higher than the symbol time, the fading
channel is considered as a frequency-selective fading channel. In this case, the
transmitted signal bandwidth is higher than the channel coherence bandwidth:
W >> Wc (1.12)
i.e.
τmax >> Ts (1.13)
16
1.3 Propagation Channel Model
expressed as the convolution between the transmitted signal and the CIR:
L
X
r(n) = h(l)s(n − l) (1.14)
l=0
The received signal is spread in time and this leads to the so-called inter sym-
bol interference (ISI). Figure (1.6) shows the frequency response of a frequency-
selective fading channel.
In both cases, estimating and correcting the CIR, also called channel equaliza-
tion, is usually performed using a known signal at the receiver, called pilot signal.
Given (1.11) and (1.14), we can expect a higher complexity of the equalization
when dealing with frequency-selective channel.
−10
Power (dB)
frequency response of
the transmitted signal
−20
−30
−0.5 −0.25 0 0.25 0.5
Normalized frequency
17
1. Multicarrier High-Data Rate Mobile Wireless Systems
where the Doppler shift f DOP = v fcc is the maximum measure in Hertz of a relative
frequency shift between the transmitted and the received signal, v is the relative
speed between the transmitter and the receiver, fc is the frequency carrier and c
is the speed of the light.
The Doppler shift is caused not only by the transmitter-receiver relative motion,
but also by the movement of surrounding objects. It leads to frequency offsets at
the receiver.
If the channel coherence time is higher than the transmitted symbol time:
Ts < Tc (1.16)
Ts > Tc (1.17)
In both cases, at the receiver, the frequency offset caused by the Doppler shift
has to be estimated in order to recover the original transmitted signal.
The design of a mobile wireless system implies the above technical challenges.
In the following, we present some schemes that have been proposed during the
last years to take into account the impact of the fading effects.
18
1.4 Multiple Access Techniques for Mobile Systems
Power
user 1 user 2
Frequency
FD
M
A
Fr
ma
e
du
ra
tio
n
User 1
Time bandwidth
Available bandwidth
19
1. Multicarrier High-Data Rate Mobile Wireless Systems
At the beginning of the ’90s, the 2nd generation mobile wireless systems (2G)
were developed, such as the digital AMPS in the United States, the global sys-
tem for mobile communications (GSM) in Europe and the personal digital cellular
(PDC) system in Japan [Hanz 98]. These systems employed TDMA as their mul-
tiple access scheme. When using TDMA, the whole bandwidth is assigned and
the time-domain transmission frame is divided into time slots, each assigned to
one user to transmit the data information [Stub 02]. See figure 1.8. TDMA is
used in the evolution of the 2G GSM standard, namely the general packet radio
service (GPRS) and the enhanced data rates for GSM Evolution (EDGE) sys-
tems.
Power
user 1 user 2
user U
Frequency
U
TD
se
M
r 3
A
tim
Fr
am
e
sl
e
ot
du
ra
tio
n
Available
Time bandwidth
Nevertheless in FDMA and TDMA, the number of frequency bands or time slots
is fixed for a given system, and one frequency band/time slot is assigned to one
user during the whole period of communications. This guarantees the service
quality for real-time and constant-bit-rate voice telephony. However, as the num-
ber of services is increasing from simple voice to multimedia traffic with different
requirements, fixed frequency band or time slot assignments have shown their
limitations, especially with the increasing number of users in the system. For
that reason, CDMA scheme, based on spread spectrum technology, has emerged.
In CDMA systems, the relatively narrow-band users information is spread into a
much wider spectrum using a high chip rate spreading code. When using different
20
1.4 Multiple Access Techniques for Mobile Systems
codes, multiple user information can be transmitted on the same frequency band
at the same time. See figure 1.9. The spreading code of each user is orthogonal1
to the codes of all other users to minimize the multiple access interference (MAI)
produced by other users [Schu 05]. The 2G American system IS-95 was the first
mobile cellular communication system to use CDMA technology followed by the
CDMA 2000 technology [Stub 02].
Frequency
C
D dur
M at
A io
Fr n
am
e
Available bandwidth
Time
The digital 2G has shown higher transmission capacity and better voice quality
than the analog 1G. Like the 1G, 2G was primarily designed to support voice
communication. In the last releases of these standards, capabilities were intro-
duced to support data transmission. However, the data rates were generally lower
than those supported by the existent bandwidth. Thus, an initiative of the inter-
national telecommunication union (ITU) laid the way for evolution to 3rd genera-
tion mobile wireless systems (3G). Requirements such as a specific high-data rate
and support for vehicular mobility were established. Both the GSM and CDMA
camps formed their own separate 3G partnership projects, the 3GPP and 3GPP2,
respectively. Within the 3GPP evolution track, the 2G GSM/GPRS/EDGE fam-
ily is based on TDMA and FDMA, whereas the 2G IS-95/CDMA 2000 is based
on TDMA and CDMA in the 3GPP2. See table 1.2.
1
P∞
Two vectors x(n) and y(n) are orthogonal if their dot product hx(n), y(n)i= −∞ x(n)y(n)
is equal to zero.
21
1. Multicarrier High-Data Rate Mobile Wireless Systems
The 3G in the 3GPP was first referred to the universal mobile telecommunica-
tion system wideband CDMA (UMTS WCDMA). Then, it evolved in the high-
speed downlink and uplink packet access (HSDPA and HSUPA) enhancements
and in the high-speed packet access plus (HSPA+) enhancement. Meanwhile in
the 3GPP2, the 3G was known as the CDMA evolution-data optimized (CDMA
EVDO) [Sesi 09]. All these systems use CDMA as their multiple access tech-
nique. The number of orthogonal codes used for uplink transmission1 from one
user is limited to a few codes or complex techniques are used to limit the uplink
signal peakiness and to improve the low noise amplifier efficiency2 . In addition,
CDMA usually employs a rake receiver [Proa 95] technique to suppress multipath
effects and the related cost increases with the number of paths. Therefore, the
complexity of the receiver can be very high for a high-data rate mobile wireless
systems.
To avoid the ISI, a solution could be to divide the entire bandwidth W in several
subbands of size ∆f << Wc and to send a large number of narrow-band signals
over several parallel subcarriers in the frequency band assigned to the transmis-
sion. This is the concept of a multicarrier modulation.
The preferred case of multicarrier modulation is the one that uses an orthogonal
1
In the uplink transmission, the transmitter is the user and the receiver is the base station.
2
For example, when using CDMA for an uplink transmission, a complex successive inter-
ference cancellation method is required to improve the performance of the receiver.
22
1.4 Multiple Access Techniques for Mobile Systems
basis, namely the orthogonal frequency division multiplexing (OFDM)1 [Bing 90].
In this scheme, the input symbols are transmitted at the same time over orthogo-
nal subcarriers [Wein 71]. See figure 1.10. The main idea of OFDM is to convert a
frequency-selective channel in the time domain into a collection of frequency-flat
channels in the frequency domain.
OFDM increases robustness against multipath distortions, making the system ro-
bust against ISI. In addition, OFDM systems use a cyclic prefix (CP) to combat
interblock interference (IBI). The CP consists in prefixing the OFDM symbol
with the end of it2 . The above advantages allow the channel equalization to be
easily performed in the frequency and time domains through a bank of one-tap
multipliers [Cimi 85]. Furthermore, OFDM exploits the spectral diversity and
allows an independent selection over each subcarrier of resources, such as power,
constellation size and necessary bandwidth [Kell 00], in order to maximize the
1
Orthogonality between two frequencies fk , fk′ is defined as ej2πfk t , ej2πfk′ t = δk,k′ , where
δk,k′ = 0 if k = k ′ , and δk,k′ = 1 if k 6= k.
2
The CP is detailed in subsection 1.4.2.1.
23
1. Multicarrier High-Data Rate Mobile Wireless Systems
Power
pilots
user
Frequency
O
FD
M
O
FD
sy
m
M
bo
fra
ld
m
ur
e
at
du
io
ra
n
tio
n
Subcarrier spacing
Time
link efficiency.
OFDM has been adopted in several communication standards such as digital
audio broadcasting (DAB) [ETSI 95], DVB-T [ETSI 97], and the WLAN IEEE
802.11a [IEEE 99].
The OFDM concept has been extended to multiuser communication scenarios. In
the following subsections three different multiple access schemes based on OFDM
are presented:
• OFDMA,
• and OFDM-IDMA.
This scheme was originally suggested for cable TV (CATV) networks [Sari 98]
and in the uplink communication of the Interaction Channel for Digital Terres-
trial Television (DVB-RCT) [ETSI 01]. The institute of electrical and electronics
engineers (IEEE) 3G project, called the mobile worldwide interoperability of mi-
crowave access (WiMAX) [IEEE 06] uses ODFMA technology. The last releases
of the 3GPP and 3GPP2, known as the long term evolution (LTE) and ultra
24
1.4 Multiple Access Techniques for Mobile Systems
mobile broadband (UMB) respectively, are based on OFDMA1 . Note that LTE
uses OFDMA in the downlink communication2 [3GPP 09].
In this scheme, unlike the OFDM case where all subcarriers are assigned to a
single user, subcarriers are divided in several mutually exclusive subchannels and
they are exclusively assigned to a particular user in an OFDMA network. See
figure 1.11.
OFDMA inherits from OFDM the flexibility for simultaneous transmissions and
frequency allocation algorithms aim at exploiting the spectral diversity to allo-
cate the communication resources to the different users. In addition, OFDMA
inherits the ability to compensate channel distortions in the frequency domain
without computationally demanding time domain equalizers.
In subsection 1.4.2, more details about the OFDMA system are presented.
pilots
Power user 1
user 2
user 3
user 4
Frequency
O
FD
M
O
FD
A
sy
M
m
A
bo
fra
ld
m
ur
e
at
du
io
ra
n
tio
n
Subcarrier spacing
Time
25
1. Multicarrier High-Data Rate Mobile Wireless Systems
munication [3GPP 09]. In this multiple access scheme, the symbols pass through
a discrete Fourier transform before going through the standard OFDMA mod-
ulation. This is often viewed as a DFT-coded OFDM. Thus, SC-FDMA in-
herits all the advantages of OFDMA over other well-known techniques such as
TDMA and CDMA. SC-FDMA brings additional benefit of low peak-to-average
power ratio (PAPR) compared to OFDMA, making it suitable for uplink trans-
missions [Sesi 09]. The SC-FDMA transceiver has similar structure as a typical
OFDMA system except the addition of a new DFT block before subcarrier map-
ping. Hence, SC-FDMA can be considered as an OFDMA system with a DFT
mapper.
When using OFDMA and SC-FDMA, an MAI free transmission can be achieved
by allocating different subcarriers to different users. Thus, if there are more than
one user in the system, the entire bandwidth has to be shared by all the users.
This may limit the data rate.
In 2002, Ping et al. [Ping 02a] proposed the interleave-Division Multiple Access
(IDMA) for asynchronous1 spread-spectrum mobile systems. Whereas transmit-
ters can be defined by using different orthogonal codes in CDMA systems, they
are distinguished by a different chip-level interleaver in IDMA systems. While
CDMA allows the MAI to be suppressed using the different codes, IDMA requires
an iterative process. Like CDMA, the entire bandwidth can be allocated to a sin-
gle user, achieving a very high single-user capacity when using IDMA.
IDMA systems require a method to suppress multipath effects and avoid the ISI.
One solution may be to use an ISI cancellation method as in [Ping 06]. However,
the corresponding computational cost of this method increases linearly with the
number of paths and may be unsuitable for high-data rate systems.
To solve this constraint, a scheme that combines OFDM and IDMA has been
proposed by Mahafeno et al. [Maha 06]. This architecture combines most of
the advantages of the OFDM and the IDMA and avoids their individual disad-
vantages. When an OFDM-IDMA is considered, ISI is resolved by an OFDM
1
The users of an IDMA system do not have to be time-synchronized one to another.
26
1.4 Multiple Access Techniques for Mobile Systems
This subsection describes the OFDMA transmitter model. After channel coding
and modulation, the symbols are grouped into blocks of length Ku < K, where
PU
u=1 Ku = K.
27
1. Multicarrier High-Data Rate Mobile Wireless Systems
Power user
user 2
user 3
user 4
Frequency
O
FD
M
O
FD
A
sy
M
m
A
bo
fra
ld
m
ur
e
at
du
io
ra
n
t
io
n
Subcarrier spacing
Time
• Interleaved CAS: the subcarriers of each user are uniformly spaced over the
signal bandwidth at a distance K from each other. This method can exploit
the channel frequency diversity. See figure 1.13.
Power user
user 2
user 3
user 4
Frequency
O
FD
M
O
FD
A
sy
M
m
A
bo
fra
l
m
du
e
ra
du
tio
ra
n
tio
n
Subcarrier spacing
Time
28
1.4 Multiple Access Techniques for Mobile Systems
• Generalized CAS: each user can select the best subcarriers that are currently
available, e.g. those with the highest signal-to-noise ratios (SNRs). In this
allocation strategy, there is no rigid association between subcarriers and
users; the generalized CAS allows dynamic resource allocation and provides
more flexibility than the other CAS. When using generalized CAS, a priori
knowledge of the propagation channel is necesssary, i.e. the transmitter
needs a channel feedback. See figure 1.14.
Power user
user 2
user 3
user 4
Frequency
O
FD
M
O
FD
A
sy
M
m
A
bo
fra
ld
m
ur
e
at
du
io
ra
n
tio
n
Subcarrier spacing
Time
Then, the CAS maps the Ku symbols (M-PSK symbols) of each block to the
subcarriers assigned to the uth user. This operation is easily performed by ex-
tending the block with the insertion of K −Ku zeros. The resulting column vector
of dimension K is defined as follows:
∆
Su = [Su (0), Su (1), . . . , Su (k), . . . Su (K − 1)]T (1.19)
according to the frequency allocation of each user, Su (k) is non-zero if the kth
subcarrier is allocated to the uth user. Su (k) is the symbol associated to the kth
subcarrier. The resulting block has a time duration T = KTs . In the following
let T be the OFDMA symbol time duration.
29
1. Multicarrier High-Data Rate Mobile Wireless Systems
Then, the IFFT output samples are gathered into a vector, resulting in the
OFDMA symbol:
∆
Xu = [Xu (0), Xu (1), . . . , Xu (n), . . . , Xu (K − 1)]T
(1.21)
= FH Qu Su
In addition, Qu is the CAS diagonal matrix where the kth coefficient of the main
diagonal is given by: (
1 if Su (k) 6= 0
Qu (k, k) = (1.23)
0 elsewhere
Adding a Cyclic Prefix (CP)
At the receiver due to the multipaths, the lth copy of a determined OFDMA
symbol2 arrives with a time delay τl , where l ∈ {0, . . . Lu − 1} and Lu is number
of paths of the channel. See section 1.3. Then, some samples of the OFDMA
symbol interfere with the adjacent OFDMA symbol causing IBI.
The next step of the OFDMA transmitter consists in introducing the cyclic prefix
(CP), that is a repetition of the Ng last samples of the vector Xu at the beginning
of it. The CP makes the OFDMA symbol robust to IBI. By inserting a CP of
1
It should be noted that FFH = IK .
2
The users are assumed to be time-synchronized one another.
30
1.4 Multiple Access Techniques for Mobile Systems
Figure 1.15: Cyclic prefix, where m denotes the OFDMA symbol index
size Ng ≥ max
u
(Lu ), the interferences between adjacent OFDMA symbols are
avoided. To maintain acceptable data throughput, the CP length must be chosen
just greater than the CIR length. See figure 1.15 which shows the advantage of
adding the CP.
Thus, a CP of size Ng is inserted to Xu such that:
31
1. Multicarrier High-Data Rate Mobile Wireless Systems
CP-OFDMA modulation
symbols to be
transmitted Su Xu XCP
u
ADD P/S
CAS IFFT
CP
Figure 1.16: Simplified OFDMA transmitter for the uth user, where Xu is the
OFDMA symbol
1
∆f = T
Power
User 1
User 2
User 3
32
1.4 Multiple Access Techniques for Mobile Systems
LO
Tolerance= fufc (ppm) fuLO (kHz) fuDOP + fuLO (kHz) ǫu
4 10.40 11.12 0.9966
3 7.80 8.52 0.7636
2 5.20 5.92 0.5306
1 2.60 3.32 0.2977
0.5 1.30 2.02 0.1812
Table 1.4: CFO values for different values of LO tolerance, using the IEEE 802.16
Wireless MAN system parameters, vu = 300 km/h and fc =2.6 GHz
After down-conversion and low-pass filtering, the received signal passes into an
1
analog-to-digital (A/D) converter, where it is sampled at the frequency Fs = Ts
.
The signal received by the BS is a superposition of the contributions from the
U active users. After the time synchronization, the received OFDMA symbol of
size KT can be expressed as:
U
X
RCP = RuCP + BCP (1.26)
u=1
∆
RuCP = [Ru (K − Ng ), . . . , Ru (K − 1), Ru (0), . . . , Ru (K − 1)]T (1.27)
At the BS, the signal transmitted by each user has been affected each propagation
channel, thus:
LX
u −1
ǫu n
RuCP (n) = ej2π K hu (l)Xu (n − l) − Ng ≤ n ≤ K − 1 (1.28)
l=0
where hu (l) is the lth coefficient of the CIR, ǫu = K(fuDOP + fuOSC )Ts is the carrier
frequency offset (CFO) normalized to the subcarrier spacing, f DOP is the Doppler
shift and fuOSC is the frequency difference in Hertz between the carrier frequencies
of the LOs . Tables 1.4 and 1.5 show examples of the CFO magnitudes.
Then, at the OFDMA uplink receiver a CFO estimation/correction is necessary
(See chapter 2). Let us consider that this step is performed perfectly. Then,
33
1. Multicarrier High-Data Rate Mobile Wireless Systems
Table 1.5: CFO values for different relative speeds, using the IEEE 802.16 Wire-
less MAN system parameters and fuLO = 0 Hz
removing the CP, the contribution from each user to the nth sample of the received
OFDMA symbol can be expressed as:
LX
u −1
Then, the contribution from each user to the received OFDMA symbol at the
BS satisfies:
∆
Ru = [Ru (0), Ru (1), . . . , Ru (n), . . . , Ru (K − 1)]T
(1.30)
= hu Xu
1
Due to the CP, at the reception, the convolution between the transmitted signal and the
channel can be expressed a circular convolution.
34
1.4 Multiple Access Techniques for Mobile Systems
Finally, given (1.30) the received OFDMA symbol can be expressed as:
∆
R = [R(0), R(1), . . . , R(n) . . . , R(K − 1)]T
U
X
= Ru + B
u=1
(1.32)
XU
= hu Xu + B
u=1
where b = FB.
Then, given (1.21), (1.22) and (1.23), (1.33) can be rewritten as:
U
X
r= Fhu F Qu Su + b
| {z }
u=1
Hu
(1.34)
U
X
= Hu Qu Su + b
u=1
where Hu = diag {[Hu (0), hu (1), . . . , Hu (K − 1)]} is the channel frequency re-
sponse matrix, the entries of which are the channel frequency responses associated
to the kth subcarrier:
L−1
X kl
Hu (k) = hu (l)e−j2π K (1.35)
l=0
where Su (k) 6= 0 if the kth subcarrier is allocated to the uth user and
P −j2π nk
b(k) = K−1
n=0 B(n)e
K .
35
1. Multicarrier High-Data Rate Mobile Wireless Systems
After subcarrier selection, the symbols r(k) pass through channel equalizers
(EQU) to obtain the transmitted symbols. It should be noted that due to the
diagonal scheme of the matrix Hu the channel equalization can “easily” be per-
formed. The OFDMA uplink receiver block scheme is shown in figure 1.18.
CFO R r subcarrier
Removing selection
estimation/ FFT
CP Qu
correction
RCP CP-OFDMA demodulation
Time
synchronization EQUs
36
1.4 Multiple Access Techniques for Mobile Systems
the subscript u denotes the information associated to the uth user. In addition,
the overall subcarriers are numbered from 0 to K − 1, i.e k ∈ {0, . . . K − 1}.
IDMA transmitter
The bits of the uth user to be transmitted are encoded by a low-rate encoder
(ENC), which consists of an association of channel coding1 and spreading code of
Q
size S. The encoded bits pass through the interleaver u and then are modulated.
The symbols are divided into blocks of length K; the resulting IDMA modulated
block is:
∆
Su = [Su (0), Su (1), . . . , Su (k), . . . Su (K − 1)]T (1.37)
CP-OFDM Modulation
An IFFT is performed over the IDMA modulated symbol as in (1.20) resulting
in the OFDM-IDMA symbol:
∆
Xu = [Xu (0), Xu (1), . . . , Xu (k), . . . Xu (K − 1)]T (1.38)
37
1. Multicarrier High-Data Rate Mobile Wireless Systems
IDMA transmitter
bits to be
transmitted Q Su CP-OFDM
Xu
ENC u modulation modulation
Q
Figure 1.19: OFDM-IDMA transmitter for the uth user, where u represents
the interleaver
where B = [B(0), B(1), . . . , B(n) . . . , B(K −1)]T is zero-mean AWGN vector with
variance σB2 IK and hu is the channel circulant matrix as is expressed in (1.31).
Thus, after the FFT step the received symbol can be expressed as:
∆
r = [r(0), r(1), . . . , r(k), . . . , r(K − 1)]T
= FR
(1.40)
U
X
= Fhu Xu + b
u=1
38
1.4 Multiple Access Techniques for Mobile Systems
39
1. Multicarrier High-Data Rate Mobile Wireless Systems
Q
where ⇒u represents the interleaving process. As there is any a priori information
at the first iteration, the initial values of all guESE (k) are set to zeros.
On the one hand, the a posteriori probability DEC is a standard function that
depends on the ENC, such as convolutional or turbo codes [Ping 01], [Berr 93].
The reader may refer to [Ping 01] and [Berr 93] for more details about the DECs.
On the other hand, let us present details about the ESE. The ESE generates coarse
estimates of the symbol Su (k). The symbol is treated as a random variable and
the a priori LLR guESE (k) about Su (k) is obtained by updating gDEC
u (k) as shown
in (1.45). Given (1.41) we have [Maha 06]:
ESE !
egu (k) − 1
∆ guESE (k)
µu (k) = E{Su (k)} = gESE (k) = tanh ∀u, k (1.46)
eu +1 2
∆
υu (k) = Var{Su (k)} = 1 − (µu (k))2 ∀u, k (1.47)
where E{.} and Var{.} represents the mean and the variance of {.}, respectively.
For the initialization process, guESE (k) = 0 for i = 1, where i denotes the iteration
number with i ∈ {1, . . . , Iidma }, where Iidma is the maximum number of iterations.
Let us recall the expression of the received symbol r(k) done for the OFDMA
uplink receiver in (1.36):
U
X
r(k) = Hu (k)Su (k) + b(k) (1.48)
u=1
40
1.4 Multiple Access Techniques for Mobile Systems
and r(k) can be considered. Thus, ςu (k) can be completely characterized by its
mean and variance:
where
U
X
E{r(k)} = Hu (k)µu (k) ∀u, k (1.52)
u=1
U
X
Var{r(k)} = σB2 + (Hu (k))2 υu (k) ∀u, k (1.53)
u=1
Algorithm 1 ESE
Require: r(k), σB2 , Hu (k), Iidma
Initialization step:
guESE (k) = 0 ∀u, k
Iterative process
for i = 1 to Iidma do ESE
µu (k) = E{Su (k)} = tanh gu 2 (k) ∀u, k
2
υu (k) = V ar{Su (k)} = 1 − (µu (k)) ∀u, k
P
E{r(k)} = Uu=1 Hu (k)µu (k) ∀u, k
P
V ar{r(k)} = σB2 + Uu=1 (Hu (k))2 υu (k) ∀u, k
E{ςu (k)} = E{r(k)} − Hu (k)µu (k) ∀u, k
V ar{ςu (k)} = V ar{r(k)} − (Hu (k))2 υu (k) ∀u, k
r(k)−E{ςu (k)}
u (k) = 2Hu (k) × V ar{ςu (k)}
gESE ∀u, k
end for
41
1. Multicarrier High-Data Rate Mobile Wireless Systems
2
σB Hu
R
Figure 1.20: OFDM-IDMA uplink receiver [Maha 06], where ESE is the elemen-
tary signal estimator and DEC is the a posteriori probability decoder
1
In this PhD dissertation we consider a parallel processing of the users in the IDMA receiver.
42
1.5 Conclusions
1.5 Conclusions
In this chapter, our purpose was to present the multicarrier high-data rate mo-
bile wireless systems. We recall various properties of the system such as: 1/ the
influence of the propagation channel leading to various kinds of channels (slow
or fast fading channel, flat-fading or frequency-selective channel), 2/ the various
multiple access techniques that have been derived: FDMA, TDMA, CDMA, SC-
FDMA, OFDMA and OFDM-IDMA.
Then, we focus our attention on two systems. We first present in detail the
OFDMA uplink transmitter and receiver, by giving the expression of the input
and the output of each block in the processing chain. We also look at the OFDM-
IDMA system and present the uplink receiver proposed in [Maha 06]
In both cases, we point out the importance of the estimation of the CFO and the
channel.
In chapter 2, we are going to propose several methods based on Kalman Filtering
and H∞ filtering to carry out the frequency synchronization and channel equal-
ization.
In chapter 3, we will propose some modifications of the conventional OFDM-
IDMA uplink receiver to counteract the effect of the multiple CFOs and the
channel in the received OFDM-IDMA signal.
43
Chapter 2
Frequency Synchronization and Channel
Equalization of OFDMA Uplink Systems
Contents
2.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . 47
2.2 State of the Art on the OFDMA Uplink Systems -
Carrier Frequency Offset (CFO) and Channel Esti-
mation . . . . . . . . . . . . . . . . . . . . . . . . . . . . 47
2.3 A Joint CFO/Channel Estimator for OFDMA Systems 53
2.3.1 State-Space Representation of the System . . . . . . . 54
2.3.2 Kalman Filtering in the Non-Linear Case . . . . . . . 56
2.3.3 H∞ Filtering in the Non-Linear Case . . . . . . . . . . 57
2.3.4 Simulation Results . . . . . . . . . . . . . . . . . . . . 58
2.3.5 Conclusions . . . . . . . . . . . . . . . . . . . . . . . . 65
2.4 An OFDMA Non-Pilot Aided Iterative CFO Estimator 66
2.4.1 Problem Formulation . . . . . . . . . . . . . . . . . . 67
2.4.2 User Detection - The Minimum Mean Square Error
Successive Detector (MMSE-SD) . . . . . . . . . . . . 68
2.4.3 CFO Estimation with the MMSE-SD Preambles . . . 71
2.4.4 Simulation Results . . . . . . . . . . . . . . . . . . . . 72
2.4.5 Conclusions . . . . . . . . . . . . . . . . . . . . . . . . 78
2.5 An OFDMA Robust CFO Estimator . . . . . . . . . . 79
2.5.1 System Description . . . . . . . . . . . . . . . . . . . . 80
2.5.2 Joint Disturbance Detection and CFO Estimation . . 84
45
2. Frequency Synchronization and Channel Equalization of OFDMA
Uplink Systems
46
2.1 Introduction
2.1 Introduction
As mentioned in chapter 1, in OFDMA uplink systems, estimating CFO and
channel plays a key role. This chapter deals with the frequency synchronization
and the channel equalization, and presents our contributions in that field.
The chapter is organized as follows: a state of the art dedicated to the CFO
and channel estimation in OFDMA uplink systems is given. We first focus our
attention on the estimation of the CFO for systems with interleaved and gener-
alized CAS [More 04], [Zhao 06], [Jing 08]. The estimation can be done in the
frequency domain or in the time domain. Then, we present the existing methods
dealing with the joint estimation of the CFO and the channel [Pun 04a], [Pun 06],
[Sezg 08]. At that stage, we suggest a recursive estimation of the CFO and the
channel by using optimal filtering, i.e. Kalman filtering and H∞ filtering (and
their variants). These approaches are derived:
• in combination with a statistical test, i.e. the BHT and the CUSUM test,
to design a robust CFO estimator that detects the beginning and the end
of the disturbances induced by another system.
47
2. Frequency Synchronization and Channel Equalization of OFDMA
Uplink Systems
0
10
User 1
User 2
−1
10
guard interval
Power (dB)
−2
10
−3
10
−4
10
−10 −8 −6 −4 −2 0 2 4 6 8 10
Frequency (MHz)
If the CFOs are smaller than the guard intervals, the received signals can be sepa-
rated at the BS by passing through a bank of digital band-pass filters (BPFs), each
one selecting a subband. Therefore, CFO-estimation methods used in OFDM or
OFDMA downlink transmission can be used. For example, the method proposed
in [Moos 94] to estimate the CFO in an OFDM transmission can be performed.
1
Indeed, in an OFDMA downlink communication, the receiver is the user. Thus, one user
has to estimate only a single CFO and a single channel. Then, the CFO estimation can be
performed by using the methods proposed for OFDM single-user communications [Moos 94],
[VdBe 99].
2
The control data are composed by known training sequences called “preambles” and/or
pilot subcarriers.
48
2.2 State of the Art on the OFDMA Uplink Systems - Carrier
Frequency Offset (CFO) and Channel Estimation
The CFO is estimated by calculating a phase shift in the frequency domain be-
tween two successive identical OFDM symbols. The repetitive OFDM symbol
remains identical after passing through the channel except for a phase shift pro-
duced by the CFO. In [Huan 10a], the authors use control data having a tile
structure in the frequency domain. Then, a high-resolution subspace method al-
lows the CFO to be estimated. In [VdBe 99], Van de Beek et al. take advantage
of the CP and estimate the CFOs by calculating the phase shifts in the time
domain between the Ng last samples of the received OFDMA symbols and the
CP of each user.
In OFDMA uplink systems with interleaved or generalized CAS, it is impossible
to take advantage of the frequency guard intervals. In those systems the CFO
estimation is a challenging task. When using control data, several methods have
been proposed. In [More 04], the CFO is estimated in the ML sense. However,
it requires an exhaustive grid search over all the possible values of the CFO,
leading to a high computational cost. For that reason, Morelli also proposes a
reduced-complexity frequency estimator, that compares the phase shifts of sev-
eral identical received OFDMA symbols in the frequency domain. Nevertheless,
the author assumes that all the active users, except the first one, are already
synchronized in frequency. In [Zhao 06], Zhao et al. use an EKF to estimate the
CFO. In addition, Jing et al. propose an extended H∞ filter to combat intercar-
rier interference for an OFDM system [Jing 08]. However, in both approaches,
the channels are preliminary estimated and this assumption cannot be necessar-
ily satisfied in real cases. Indeed, during the uplink CFO-synchronization stage,
channel state information1 (CSI) is not available, and hence has to be estimated
either jointly with the CFO or after CFO compensation.
In order to correct the MAI, the CFO has to be compensated. In the single-user
detector method, shown in figure 2.2, the CFOs are compensated by multiplying
the complex envelope of the signal before the CP-OFDMA demodulation at the
receiver. However, the number of FFTs increases linearly with the number of
users as well as the computational complexity of the receiver. As an alternative,
new correction methods have been proposed. More particularly in [Dai 07], the
1
The CSI refers to known channel properties of a communication link. It is usually associ-
ated to the CIR.
49
2. Frequency Synchronization and Channel Equalization of OFDMA
Uplink Systems
CFO is compensated by mitigating the effects of the major lobes and the side
lobes of the frequency response of the signals received from each user. The CFO
estimation can be obtained by calculating the phase shifts of two identical re-
ceived OFDMA symbols in the time domain. In [Choi 00], CFO compensation is
based on a circular convolution after the CP-OFDMA demodulation1 , as shown
in figure 2.3. The authors suggest using the method presented in [Moos 94] to
estimate the CFOs. However, both methods only work for subband CAS. Con-
sequently, in [Huan 05], the authors extend the concept proposed in [Choi 00]
and use it in a system with interleaved or generalized CAS. In that case, the
ǫ̂1
to the 1st EQU
BPF r1 BPF Ẑ1
circular
Q1 convolution Q1
R r ǫ̂U
FFT
to the U th EQU
rU ẐU
BPF circular BPF
QU convolution QU
Figure 2.3: CFO correction by means of a circular convolution after the FFT
1
Due to the CP, the received signal without CFO correction after the CP-OFDMA demod-
PU
ulation can be expressed as: r = u=1 FEu ⊗ Zu + b, where Zu = Fhu Xu , Eu is the CFO
matrix and ⊗ represents a circular convolution.
50
2.2 State of the Art on the OFDMA Uplink Systems - Carrier
Frequency Offset (CFO) and Channel Estimation
CFOs are corrected via a circular convolution after the CP-OFDMA demodula-
tion, but an iterative interference-cancellation (IC) step is required, as shown in
algorithm 2. The authors suggest again using the method presented in [Moos 94].
In [Cao 04a], the CFO is corrected by means of a linear detection, whereas the
CFO is estimated using a high-resolution subspace method, i.e. the multiple sig-
nal classification (MUSIC) method [Ther 92]. More recently, Hou et al. have
proposed a minimum mean square error successive detector (MMSE-SD) to sup-
press the MAI, but the CFO is assumed to be known [Hou 08].
In the above approaches, channel estimation is not included. The joint esti-
mation of the CFO and the channel has been also investigated. In [Pun 04a] and
[Pun 06], the authors study how to obtain the joint ML estimations of the chan-
nels and the CFOs of multiple users. Thus in [Pun 04a], a conventional EM is
first proposed: during the E-step, the received signals transmitted by each user,
namely the “complete data”, are estimated. During the M-step, all the CFOs and
the channels are jointly estimated by using these complete data. To simplify the
optimization issue, the value of the channel is replaced by its expression depend-
ing on the CFO in the criterion to be minimized. Therefore, only the CFO of
each user has to be estimated. Even if the criterion is explicitly given, the authors
do not mention the estimation method they used. For instance, exhaustive grid
search could be considered as suggested by the same authors in [Pun 04b]. To
reduce the computational cost, the authors in [Pun 04a] use the SAGE algorithm.
51
2. Frequency Synchronization and Channel Equalization of OFDMA
Uplink Systems
52
2.3 A Joint CFO/Channel Estimator for OFDMA Systems
[Mova 08], a blind method to estimate CFOs based on a linear precoder is pro-
posed. Using two OFDMA symbols, the idea is to find a time correlation using
a precoder which gives a second-order moment based CFO estimation for each
user.
Several approaches using ML estimation of the CFO and channel estimation in
multicarrier systems have been proposed. However, only few approaches are
based on optimal recursive filtering to perform the CFO and channel estimation
in OFDMA uplink systems. One feature is that this multiparameter estimation
solution can be implemented recursively. In particular, each updated estimate
is computed from the previous estimate and the current observation. Therefore,
only the previous estimation has to be stored. The optimal filtering is hence
computationally more efficient than algorithms that compute the estimation us-
ing the entire past observed data at each step of the estimation process.
In the following, different optimal filtering techniques aiming at estimating the
CFO and the channel are proposed for conventional and CR environments.
Due to the propagation channels and the radio frequency (RF) effects, the uplink
signals at the BS are plagued by synchronization and estimation errors. There-
fore at the BS, synchronization and equalization process has to be done. These
processes usually require control data. Commonly, the synchronization process is
performed before the CIR estimation. However, if a time-synchronized OFDMA
uplink system is considered, the CFO and the channel can be jointly estimated.
Then, for a time-synchronized system the nth sample of the received OFDMA
symbol at the BS after CP removing is:
U LX
X u −1
ǫu n
R(n) = ej2π K hu (l)Xu (n − l) + B(n) 0 ≤ n ≤ K − 1 (2.1)
u=1 l=0
53
2. Frequency Synchronization and Channel Equalization of OFDMA
Uplink Systems
nh io
where Eu = diag 1, ej2πǫu /K , . . . , ej2π(K−1)ǫu /K (See chapter 1 in section 1.4.2.2).
In order to restore orthogonality among each user subcarrier, both the CFO-
synchronization error vector ǫ and the CSI vector h have to be estimated at the
BS.
PU
Let us define the vector ǫ of size U and the vector h of size L = u=1 Lu by
storing the CFOs and the CSIs of the U users in the system as follows:
ǫ = [ǫ1 , ǫ2 , . . . , ǫu , . . . , ǫU ] (2.3)
h = [h1 , h2 , . . . , hu , . . . , hU ] (2.4)
54
2.3 A Joint CFO/Channel Estimator for OFDMA Systems
Let us define the vector x(n) of size U = U + 2L, which is the state vector of the
system (2.1)-(2.2):
h iT
x(n) = ǫ(n) Re {h(n)} Im {h(n)} (2.5)
Y(ǫ, h) = C(ǫ, h) + B
(2.6)
= [Y(0, ǫ, h), Y(1, ǫ, h), . . . , Y(n, ǫ, h), . . . , Y(K − 1, ǫ, h)]
where C(ǫ, h) is a 2 × K matrix that stores the real and the imaginary parts of
the contributions of alls users to the received OFDMA symbol:
nP oT
U
Re u=1 Ru
C(ǫ, h) =
nP
oT
Im U
u=1 Ru (2.7)
= [C(0, ǫ, h), C(1, ǫ, h), . . . , C(n, ǫ, h), . . . , C(K − 1, ǫ, h)]
and " #
Re{B(n)}T
B= = [B(0), B(1), . . . , B(K − 1)] (2.8)
Im{B(n)}T
B(n) is an AWGN vector with zero-mean and covariance matrix σ2B I2 .
Now, let us introduce the state-space representation of the system (2.1)-(2.2) to
estimate the CFO and the channel:
State equation:
Measurement equation:
55
2. Frequency Synchronization and Channel Equalization of OFDMA
Uplink Systems
where w(n) is an AWGN matrix with zero-mean and covariance matrix σw2 IU . It
should be noted that σw2 is very low and can even be equal to 0. In this latter
case, the CFO and the channel are assumed to be constant during one OFDMA
symbol.
56
2.3 A Joint CFO/Channel Estimator for OFDMA Systems
Non-linear Estimation
No model noise Model noise
Linearized RLS Kalman filter
Linearization Gaussian approximation
EKF UKF
SPKF
SOEKF CDKF
IEKF
• the unscented H∞ filter [Li 10], which is based on the unscented transfor-
mation embedded into the “extended H∞ filter” architecture.
57
2. Frequency Synchronization and Channel Equalization of OFDMA
Uplink Systems
The reader is referred to table 2.2. In addition, appendix C deals with a compar-
ison between the equations that define Kalman filtering and H∞ filtering.
No statistical
No linearization
assumptions
EKF x x
SPKF X x
Extended
x X
H∞ filter
Unscented
X X
H∞ filter
In the following subsection, simulation results show the efficiency of the proposed
optimal filtering algorithms.
• and the methods presented in [Pun 04a] and [Pun 06] where a grid search
approach is used to update the CFO estimation.
Simulation protocol:
An OFDMA IEEE 802.16 WirelessMANTM uplink system composed of U = 4
users sharing K = 128 subcarriers and with a cyclic prefix Ng = K/8 ≥ max
u
(Lu )
is considered. A transmission over a Rayleigh slow-fading frequency-selective
58
2.3 A Joint CFO/Channel Estimator for OFDMA Systems
2
σB
When using H∞ filtering: Ξ = 102 , V = I
2 2
and W = IU σw2 . One assumes
that there is a state noise w(n) with very small variance, e.g. σw2 = 10−3 . For
the CFO, the initialization parameters of the algorithm is ǫ̂u (0) = 0 ∀u.
The grid search algorithms proposed in [Pun 04a] and [Pun 06] are based on
Iem = 20 iterations and a grid search precision equal to 10−4 ; this means that
2
b = 10−4
+ 1 = 20001 values of CFO are studied in the grid search algorithm.
The estimations are obtained by using control data of one OFDMA symbol.
We focus our attention on the first user in the system. First of all, our approaches
provide similar results when UKF and CDKF are used. Therefore, in the follow-
ing, we make reference to SPKF filtering performance.
Figure 2.4 and 2.5 show the results in terms of CFO and channel MMSE consid-
ering a perfect knowledge of the noise statistics. When the noise characteristics
are available, the EKF provides quite similar results in comparison with the “ex-
tended H∞ filter” and there is no real difference between the SPKF and the
“unscented H∞ filter”. For that reason, in figure 2.4 and 2.5, we only show the
results of the Kalman approaches. In addition, the algorithms in [Pun 04a] and
[Pun 06] provide similar results after a given number of iterations. Thus, figure
2.4 and 2.5 show only the performance of the EM-based algorithm presented in
[Pun 04a].
The IEKF gives the best estimation performance. However, due to the calcula-
tion of the Jacobian matrix at each iteration, the computational complexity is
the highest among the proposed approaches. The SPKF algorithms give better
CFO and channel estimations than the SOEKF and the EKF.
59
2. Frequency Synchronization and Channel Equalization of OFDMA
Uplink Systems
−2
10
EKF
SOEKF
SPKF
IEKF
EM
−3
10
MMSE
−4
10
0 2 4 6 8 10 12 14 16
SNR (dB)
0.11
EKF
SOEKF
0.1 SPKF
IEKF
EM
0.09
0.08
0.07
0.06
MMSE
0.05
0.04
0.03
0.02
0.01
0 2 4 6 8 10 12 14 16
SNR (dB)
Figure 2.6 shows that the IEKF converges faster than the other Kalman ap-
proaches. In addition, the SPKF seem to require less observations than the EKF
60
2.3 A Joint CFO/Channel Estimator for OFDMA Systems
−0.1
CFO value
SOEKF
EKF
−0.15
−0.2
−0.25
20 40 60 80 100 120
Sample Index
Figure 2.6: Recursive CFO estimation for a joint CFO/channel estimation using
optimal filtering, when the noise statistics are available
61
2. Frequency Synchronization and Channel Equalization of OFDMA
Uplink Systems
0.1
SNR = 10 dB
0.08
5 dB noise variance error
EKF
SPKF
0.06
extended H∞ filter
unscented H∞ filter
CFO value
0.04
0.02
0
CFO real value
−0.02
20 40 60 80 100 120
Sample Index
Figure 2.7: Comparison between the Kalman filtering approaches and the H∞
filtering approaches in terms of convergence speed, when the noise statistics are
not available
One of the goals in the design of wireless communication systems is to reduce the
energy consumption of the system. In current wireless communication systems,
only a few Go/s are dedicated to the channel estimation and the synchroniza-
tion. So, in the following let us have a look at the computational cost of the EM
algorithm and the SPKF algorithm. These approaches provide similar results in
terms of BER, but the EM-based algorithm have a higher computational com-
plexity due to the exhaustive grid search over the possible range of CFOs and the
iterative estimation.
Table 2.4 shows the number of arithmetic operations performed by the EM; N ob
62
2.3 A Joint CFO/Channel Estimator for OFDMA Systems
EM algorithm
E-step
additions and
Iem Nob U + Iem U (2Lu + 6) 1.1200 × 104
subtractions
multiplications
Iem Nob U + Iem U (U + 2Lu ) 1.1040 × 104
and divisions
M-Step
add./sub. for
Iem bN ob U Lu + Iem bU Lu (1 + Lu ) − Iem bU 6.2271 × 107
the grid search
mult./div. for Iem bU Nob (3U Lu + Lu )
8.0040 × 107
the grid search + Iem bU (L2u + 15U + 2Lu + 1)
other
Iem N ob U Lu + Iem U L2u 1.1366 × 106
add./sub.
other
Iem N ob U (3U Lu + Lu ) + Iem U (L2u + 15U + Lu ) 3.9329 × 106
mult./div.
Total arithmetic operations 1.4740 × 108
of arithmetic operations performed by the SPKF. The first step of the algorithm,
corresponding to the selection of the sigma points requires the Cholesky decom-
position1 of a matrix of size U × U. As the observation vector length is equal to
2, the measurement update also requires the inversion of a 2 × 2 matrix.
Table 2.6 shows the number of Giga-operations per second (Go/s) performed by
the EM for different grid search precisions.
Table 2.7 shows the number of Go/s for different numbers of users in the system
with a channel composed of 3 multi-paths. In addition, table 2.7 shows the num-
ber of Go/s for different number of channel multipaths with 4 users in the system.
1
When using the SPKF and by considering U the length of the state vector, a U×U Cholesky
decomposition is necessary to obtain the “square root” of the error covariance matrix.
63
2. Frequency Synchronization and Channel Equalization of OFDMA
Uplink Systems
SPKF algorithm
step 1: calculation of the sigma points
additions and
2N ob U 2 (4L2u + 4Lu + 1) 2.0070 × 105
subtractions
multiplications
2N ob U 2 (4L2u + 4Lu + 1) 2.0070 × 105
and divisions
step 2: estimation update
additions and 2N ob U (8U 2 L3u + 12U 2 L2u + 6U 2 Lu
6.2290 × 106
subtractions 12U L2u + U 2 + 12U Lu + 3U + 2Lu +1)
multiplications N ob U (16U 2 L3u + 24U 2 L2u + 12U 2 Lu
6.1286 × 106
and divisions 20U L2u + 2U 2 + 20U Lu + 5U + 4Lu + 2)
step 3: measurement update
additions and N ob U (44U Lu + 19U + 13)
1.8348 × 105
subtractions +N ob U (28U L2u + 96Lu )
multiplications 2N ob U (12U Lu + 7U +27) + 2N ob
2.9651 × 105
and divisions +2N ob U (4U L2u + 50Lu )
Total arithmetic operations 1.3239 × 107
Table 2.5: Number of arithmetic operations performed by the SPKF joint chan-
nel/CFO estimator
64
2.3 A Joint CFO/Channel Estimator for OFDMA Systems
Table 2.7: Number of Go/s performed by the EM and the SPKF joint chan-
nel/CFO estimators
2.3.5 Conclusions
The architectures proposed in this section require control data to perform the
estimation, but they do not need an initialization step. The robustness to un-
certainties of the H∞ filtering has been confirmed, but the selection of the noise
attenuation level remains quite difficult. When the noise characteristics are avail-
able, the best estimation performance have been obtained by the IEKF among
the proposed estimators. However, the computational complexity is the high-
est among the proposed algorithms. SPKF gives the best compromise between
computational complexity and estimation performance. They have lower compu-
tational cost than the grid search methods. This advantage is crucial, because
one of the goals in the design of wireless systems is to reduce the energy con-
65
2. Frequency Synchronization and Channel Equalization of OFDMA
Uplink Systems
: data sub-carriers
: pilot sub-carriers
Tc : coherence time of the channel
1 st OFDMA Symbol 2 nd OFDMA Symbol 3 rd OFDMA Symbol (M − 1)th OFDMA Symbol M th OFDMA Symbol
time (s)
Tofdma << Tc
66
2.4 An OFDMA Non-Pilot Aided Iterative CFO Estimator
CSI for the whole frame and to provide an initial estimation of the CFO. See
figure 2.9.
: data sub-carriers
time (s)
Tofdma << Tc
The idea of the proposed receiver is to correct the MAI and to estimate the
signals sent from all users in the system by using an MMSE-SD. These estimated
signals called the “MMSE-SD preambles” play the role of the preambles and
are used to estimate the CFOs by means of the optimal filter. The iterative
architecture is shown in figure 2.10.
The frequency synchronization in an uplink communication can be seen as a
multiparameter estimation problem. A solution is to distinguish one user from
the others before starting the synchronization process. After the estimation of the
U “MMSE-SD preambles”, each one is inserted in U optimal filtering processes,
i.e. U Kalman filters or U H∞ filters. See figure 2.11. Finally, the CFOs between
the BS and each user are estimated independently, thanks to the optimal filter
and an additional MAI cancellation process.
67
2. Frequency Synchronization and Channel Equalization of OFDMA
Uplink Systems
^
A ^
1,1 1,1
R Optimal filter
MMSE-SD ^ ^
iteration 1 A U,1 based
U,1
CFO estimation
^ ^
A 1,2 1,2
MMSE-SD Optimal filter
^
A based ^
iteration 2 U,2
U,2
CFO estimation
^ ^
A 1,i 1,i
MMSE-SD Optimal filter
^
A based ^
iteration i U,i
U,i
CFO estimation
Figure 2.10: Proposed OFDMA non-pilot aided iterative receiver using optimal
filtering
68
2.4 An OFDMA Non-Pilot Aided Iterative CFO Estimator
MMSE-SD iteration i
S^ 1,i ^ ^
R A 1,i 1,i
W1,i Q1 sign G1,i
+ -
^
Y2,i S^ 2,i ^
A 2,i Kalman 2,i
W2,i Q2 sign G2,i based
+ -
CFO
estimation
+ -
^
S ^
A
^
YU,i U,i U,i U,i
WU,i QU sign GU,i
Figure 2.11: Iterative MMSE-SD using a Kalman filtering based estimator and
a MAI suppression for BPSK modulation
iteration of the MMSE-SD, where 1 ≤ i ≤ Isd and Isd denotes the max-
imum iteration number. The user’s decoding order is denoted as ũ where
ũ ∈ {1, . . . , U }; when ũ = 1, it is associated with the maximum signal to
interference-plus-noise ratio1 (SINR) among the users whereas ũ = U rep-
resents the user with the lowest SINR. In addition the so-called (ũ + 1)th
order MMSE-SD residual at the ith iteration Yũ+1,i ∀ũ 6= 1 is the differ-
ence between the received signal R and the components transmitted by the
detected users2 (namely those corresponding to the ũ highest SINRs) and
Y1,i = R ∀i. See figure 2.11.
The IC steps can be summarized as follows:
( P
U
ũ=1 Êũ,i Hũ Qũ if ũ = 1
with Ĝũ,i = (2.14)
Ĝũ−1,i (IK − Qũ−1 ) if 2 ≤ ũ ≤ U
1
As the BS calculates the different SINRs, the signal can be decoded in an ordered way.
2
Given (1.21), (1.22) and (1.23), the received OFDMA symbol can be expressed as
PU
R = u=1 Gu Su + B, where Gu = Eu FH Hu Qu , See chapter 1 in section 1.4.2.2.
69
2. Frequency Synchronization and Channel Equalization of OFDMA
Uplink Systems
where Hũ = FH Hũ , Hũ is the channel frequency response matrix associated
to the ũth user1 , Qũ is the CAS matrix associated to the ũth user defined in
(1.23), and Ŝũ,i is the vector at the ith iteration that contains the estimated
symbols sent by the ũth user:
h i
Ŝũ,i = Ŝũ,i (0), Ŝũ,i (1), . . . , Ŝũ,i (k), . . . , Ŝũ,i (K − 1) (2.15)
nh io
Êũ,i = diag 1, ej2πǫ̂ũ,i−1 /K , . . . , ej2π(K−1)ǫ̂ũ,i−1 /K and ǫ̂ũ,i is the estimation
of the CFO associated with the ũth user at the ith iteration. A new iteration
begins when all the users have been processed and when the estimation ǫ̂ũ,i ,
using the recursive estimator approach proposed in subsection 2.4.3, has
been done.
Then, (2.16) is used to decode the selected user and to obtain the estimated
signal sent by uth user4 .
70
2.4 An OFDMA Non-Pilot Aided Iterative CFO Estimator
In the next subsection, those estimated signals the “MMSE-SD preambles” are
used to estimate the CFO of each user.
Measurement equation:
ǫũ,i n
R̂ũ,i (n) = Âũ,i (n)ej2π K + Bũ (n)
(2.21)
= f (ǫũ,i (n)) + Bũ (n)
where R̂ũ,i (n) is the estimated received signal from each user (see subsection
2.4.3.1). In addition, w and Bũ are white zero-mean gaussian noises with variances
assumed to be σw2 and σũ2 respectively. It should be noted that σw2 is very low and
can even be equal to 0. In this latter case, the CFO is assumed to be constant
during one OFDMA symbol.
Let us introduce a third state-space equation that is used in the general case
for H∞ filtering to focus on a linear combination of the state-vector components.
Here, as the dimension of the state vector is equal to 1, one has:
where L is equal to 1.
As the state-space model is non-linear due to (2.21), the study of several kinds
of algorithms is proposed:
71
2. Frequency Synchronization and Channel Equalization of OFDMA
Uplink Systems
• the EKF approaches including the EKF, the SOEKF and the IEKF, pre-
sented in appendix D,
• the SPKF, including the UKF and the CDKF, presented in appendix E,
(est)
MAI estimation: R̂ũ,i (n) = Âũ,i (n)ej2πnǫ̂ũ,i (n−1)/K (2.23)
U
X (est)
MAI correction: R̂ũ,i (n) = R(n) − R̂ũ,i (n) (2.24)
j=1,j6=ũ
After some recursions, the algorithm makes its possible to estimate the value of
the uth user CFO.
72
2.4 An OFDMA Non-Pilot Aided Iterative CFO Estimator
When using H∞ filtering: Ξ = 102 , V = σũ2 and W = σw2 . One assumes that
there is a state noise w(n) with very small variance, e.g. and σw2 = 10−3 . It is
assumed that the channel estimation and a CFO pre-estimation have been per-
formed using a preamble. The algorithm estimates the “MMSE-SD preambles”
of the 4 users, and then estimates the CFOs for each OFDMA symbol so that the
coherent detection can be robust against variations over all the OFDMA frame.
1.2
Eb/N0 = 10 dB Doppler User 4
0.8
Eb/N0 = 0 dB
Estimated doppler
0.6
Doppler User 1
0.4
Doppler User 3
0.2
0 Doppler User 2
−0.2
0 50 100 150 200 250 300 350 400 450 500
Sample index
Figure 2.12: CDKF based approach for CFO estimation with a known preamble
In this first case, it is assumed that the MMSE-SD preambles have been perfectly
estimated at the receiver. The CFO estimation algorithms for users with different
CFO are tested. The users’ CFO estimation errors are considered fixed during
one OFDMA symbol. 500 montecarlo runs are performed.
Based on preliminary tests, it was noticed that the Kalman filtering approaches,
“the extended H∞ filter” and the “unscented H∞ filter” provide very similar
results. In the following, we choose the SPKF and the unscented H∞ filter to
perform the CFO estimation, in order to avoid the linearization step. Therefore,
in figure 2.12, we only show the results obtained when using a CDKF. According
73
2. Frequency Synchronization and Channel Equalization of OFDMA
Uplink Systems
0.8
0.6
Iteration 5
Estimated doppler
CFO
0.4 Iteration 4
0.2
0
Iteration 3
−0.2 Iteration 2
Iteration 1
−0.4
−0.6
Figure 2.13: MMSE-SD combined with an UKF based approach for CFO esti-
mation
to figure 2.12, our approach makes it possible to accurately estimate the CFO
recursively. It should be noted that this is the recursive CFO estimation over
only one OFDMA symbol.
In this second case, the receiver estimates the MMSE-SD preambles. The simu-
lation results deal with the first user of the system. The users’ CFO errors vary
between the different OFDMA symbols (but they are considered fixed during one
OFDMA symbol). They are modeled as independent zero-mean Gaussian ran-
2
dom variables with a variance of σcf o . Figure 2.13 shows the results in terms of
CFO estimation over one OFDMA symbol. The proposed algorithm provides an
estimation of the CFO in an iterative way. As expected, the first iteration leads
to poor performance, but iterating our approach especially up to 5 times leads to
good performance.
In figure 2.14, the performance of our algorithm in terms of BER are shown; the
estimation is done by using either the Kalman filtering with an error of 3 dB over
the variance or the H∞ filtering. In the following, we choose “the unscented H∞
74
2.4 An OFDMA Non-Pilot Aided Iterative CFO Estimator
−1
10
BER
H∞ filter - iteration 1
H∞ filter - iteration 2
H∞ filter - iteration 3
H∞ filter - iteration 4
H∞ filter - iteration 5
−2
10 No CFO
KF - iteration 5
CFO perfectly estimated
4 6 8 10 12 14 16
Eb/N0 (dB)
filter” to perform the CFO estimation, in order to avoid the linearization step. In
addition, we show the results when a single-user detector is used by considering
a perfectly estimated CFOs. According to figure 2.14, the results of the proposed
method are better after performing the fourth iteration. A gain of around 1 dB for
a BER of 10−2 is obtained when using the proposed method with the “unscented
H∞ filter” at the fifth iteration in comparison with the single-user detector. In
[Huan 05] and [Cao 04a], the algorithms tend to the theoretical value. However,
in [Huan 05] more subchannels than users are required. In addition, in [Cao 04a]
a training sequence is necessary.
In figure 2.15, the results in terms of MMSE are shown. Based on previous tests,
it was seen that the results in [Mova 08] are better than the results in [Cao 04b]
in terms of CFO estimation performance. So, we propose a comparative study
with the multiuser interference resilient (MUI) approach [Cao 04b]. Note that
the algorithm presented in [Cao 04b] needs two OFDMA symbols for the estima-
tion. A gain of 4 dB for an MMSE of 10−4 is obtained between the MUI method
and the proposed method at the third iteration when using a SPKF. In addition,
figure 2.16 compares the SPKF and the “unscented H∞ filter” in terms of con-
vergence speed, for a Eb /No = 10 dB, for different errors over the noise variance
75
2. Frequency Synchronization and Channel Equalization of OFDMA
Uplink Systems
−3
10
KF iteration 3
KF iteration 4
KF iteration 5
MUI
−4
10
MMSE
−5
10
−6
10
4 6 8 10 12 14 16
Eb/N0 (dB)
Figure 2.15: CFO estimation performance for a non-pilot aided CFO estimation
when using the SPKF and for different values of Ξ when using the H∞ filtering.
Firstly, one considers an error over the variance of 3dB and 5dB; secondly, no
error over the variance is assumed. It can be seen that the convergence speed of
the SPKF is affected by the variance noise error. For the lowest values of Ξ, the
convergence may be slightly faster. If the value of Ξ increases, the “unscented
H∞ ” filter performance tend to be the same as the SPKF ones.
In this third case, we propose to analyze the influence of the CSIs variations over
the receiver performance. The CFO estimation is performed by a SPKF, provided
that the noise characteristics are known. The simulation results are focused on
the first user of the system. It is assumed that the CSIs have small variations in
time over the OFDMA frame. ĥu denotes the estimated channel impulse response
for each user using the preamble and hũ is the value of the channel response for
each user. In a first test, the following normalized random error over each channel
is considered |ĥu − hũ | < ILu ×1 10−2 , where ILu ×1 is a Lu × 1 vector composed of
ones. In the second test |ĥu − hũ | < ILu ×1 10−1 . The results obtained at the fifth
iteration for the different BER with a different error over the channels are shown
76
2.4 An OFDMA Non-Pilot Aided Iterative CFO Estimator
in figure 2.17. One can notice that the performance of the architecture are not
0.4
0.35
0.3
0.15
0.1
0.05
0
0 50 100 150 200 250 300 350 400 450 500
Sample Index
(a)
0.4
0.35
CFO H∞ filter Ξ = 20
H∞ filter Ξ = 102
Estimated doppler
0.3
SPKF with a
an estimation error over the
SPKF with a noise variance of 5 dB
an estimation error over the
noise variance of 3 dB
0.25
SPKF with a
an estimation error over the
0.2 noise variance of 0 dB
(b)
Figure 2.16: Comparison between the SPKF and the “unscented H∞ filter” in
terms of convergence speed
77
2. Frequency Synchronization and Channel Equalization of OFDMA
Uplink Systems
−1
10
Theoretical value
Channel perfectly estimated
Channel estimation error < 10−2
Channel estimation error < 10−1
BER
−2
10
4 6 8 10 12 14 16
Eb /No (dB)
Figure 2.17: BER performance when the channel is estimated and a non-pilot
aided CFO estimation is considered
Now let us look at the complexity of the frequency synchronization scheme. Con-
cerning with the MMSE-SD step, at each iteration and for each user, the MMSE-
SD is based on the inversion of a matrix of size K × K in order to calculate the
suppression weight matrix Wũ,i . This matrix inversion requires O(K 3 ) opera-
tions. Thus, the overall complexity of the MMSE-SD step is O(Isd U K 3 ).
Then, concerning the estimation step, as the dimension of the state vector is equal
to 1, the optimal filtering step requires the computation of O(1) operations for
each user and for each iteration. Thus, the overall complexity of the estimation
step is O(Isd U ).
Therefore, we can conclude that the overall complexity of the proposed iterative
architecture is O(Isd U K 3 ).
2.4.5 Conclusions
The non-pilot aided MAI suppression scheme for an OFDMA uplink transmission
followed by the CFO estimation algorithm based on optimal filtering can estimate
78
2.5 An OFDMA Robust CFO Estimator
and detect each user’s CFO and frame respectively, with no need of pilot sub-
carriers. This leads to a maximum transmission data rate.
In addition, simulation results confirm that the proposed scheme can effectively
suppress the MAI caused by a relatively large CFO, sufficiently robust to CFO
variations. The decoding of OFDMA signal is an ordered serial processing that
combines interference suppression and interference cancellation techniques. The
iterative decoding applies simple hard interference cancellation techniques, result-
ing in moderate complexity. The CFO estimation performance of the different
approaches give quite similar results. On the one hand, one advantage of the
SPKF is that it does not need the calculations of Jacobians or Hessians. In those
cases, the computational complexity is lower. On the other hand, when using the
H∞ filter, no information about the noise statistics is required, but the selection
of the noise attenuation level may be a problem.
79
2. Frequency Synchronization and Channel Equalization of OFDMA
Uplink Systems
80
2.5 An OFDMA Robust CFO Estimator
Figure 2.18: PU and CR system spectrum, where PUA and PUB denote the
frequency bands used by other PUs
hupu=[hpu pu pu pu pu pu
u (0), hu (1), . . . , hu (l ), . . . , hu (Lu − 1)] (2.25)
where Lpu
u is the length of the maximum channel delay spread.
To help the reader, the PU system parameters are summarized in table 2.8.
At the receiver, due to the propagation conditions, time offset and CFO are
induced in the signal. The nth sample of the mth PU-OFDMA received symbol
after time synchronization and cyclic prefix removal can be written as follows
(See chapter 1 in section 1.4.2.2):
pu −1
KX
ǫu n kpu n
pu j2π K
Ru (m , n) = e pu
Supu (mpu , k pu )Hupu (k pu )ej2π K pu (2.26)
kpu =0
1 pu
It should be noted that Tof dma is the time duration of one PU-OFDMA symbol and is not
the time duration of one PU-OFDMA frame.
81
2. Frequency Synchronization and Channel Equalization of OFDMA
Uplink Systems
Description Parameter
Bandwidth W pu
Symbol time T = W1pu
pu
Number of subcarriers K pu
Number of emitted OFDMA symbols M pu
Subcarrier index k ∈ {0, . . . K pu − 1}
pu
Concerning the CR-signal, the analog baseband signal emitted I(t) by the
CR-antenna, without cyclic prefix, at time t is given by:
cr −1
KX
cr (t−mcr T cr
I(t) = S cr (mcr , k cr )ej2πfk of dm ) (2.28)
kcr =0
where Skcrcr (mcr ) is the emitted CR-symbol with a time duration T cr at the k cr th
subcarrier associated to the mcr th CR-OFDM symbol, where M cr is the number
of emitted CR-OFDM symbols. Moreover, Skcrcr (mcr ) = 0 if the spectrum sensing
decides that the k cr th subcarrier is busy and mcr Tof
cr cr cr
dm ≤ t < (m + 1)Tof dm .
cr cr
Tof dm is the CR-OFDM symbol time and fk is the baseband frequency of the
k cr th subcarrier.
At the PU-BS antenna, the CR-received signal after the convolution with the
82
2.5 An OFDMA Robust CFO Estimator
Description Parameter
Bandwidth W cr
Symbol time T = W1cr
cr
Number of subcarriers K cr
Number of emitted OFDM symbols M cr
Subcarrier index k cr ∈ {0, . . . K cr − 1}
OFDM emitted symbol index mcr ∈ {0, . . . M cr − 1}
cr
Emitted symbol at the k th subcarrier
Skcrcr (mcr )
associated to the mcr th CR-OFDM symbol
cr cr cr
OFDM symbol time Tof dm = K T
Baseband frequency of the subcarrier cr k
fk = T cr = K crkT cr
of dm
Description Parameter
Bandwidth W cr = W pu /α
Symbol time T cr = αT pu
pu pu
OFDM symbol time Tof dm = γTof dma = α(K cr /K pu )Tof
cr
dma
where hcr
lcr are the coefficients of the CR-channel impulse response,
lcr ∈ {0, 1, . . . , Lcr − 1} and Lcr is the length of the maximum channel delay
spread.
Once again, to help the reader, the CR system parameters are summarized in
table 2.9. In addition, table 2.10 shows the relations between the PU and the
CR-system parameters.
The received signal is sampled at the sampling frequency:
1
fs = (2.30)
T pu
83
2. Frequency Synchronization and Channel Equalization of OFDMA
Uplink Systems
Then, using (2.29), (2.30) and tables 2.8-2.10, the unknown downsampled
CR-NBI1 produced by a fault of spectrum sensing detection can be expressed
as follows:
LX −1 cr KX −1 cr
n cr n cr pu
V (mcr , n) = Z(mcr Tof
cr
dm + )= hcr (lcr ) S cr (mcr , k cr )ej2πfk ( fs −l T )
fs lcr =0 kcr =0
cr −1
KX cr −1
LX
kcr lcr T pu kcr nT pu
cr cr cr
= S (m , k ) hcr (lcr )e−j2π K cr T cr ej2π K cr T cr
kcr =0 lcr =0
cr −1
KX
kcr n
= S cr (mcr , k cr )Hkcrcr ej2π αK cr
kcr =0
(2.31)
PLcr −1 cr cr
cr cr cr cr −j2π kαKlcr
where H (k ) = lcr =0 h (l )e represents the unknown CR-channel
frequency response associated to the k cr th subcarrier.
Thus, the PU-OFDMA received signal at the BS satisfies:
(
f (mpu , n, ǫ, hpu ) + V (mcr , n) + B(mpu , n) if n1 ≤ n ≤ n2
R(mpu , n) =
f (mpu , n, ǫ, hpu ) + B(mpu , n) otherwise
(2.32)
PU
where f (mpu , n, ǫu , hupu ) = u=1 Ru (mpu , n) and B(mpu , n) is a zero-mean com-
plex AWGN with variance σB2 . In addition, n1 ∈ {0, 1, . . . , K pu (1 − γ)} is the
first sample of the CR-NBI and n2 = n1 + γK pu − 1 is the last sample of the
CR-NBI, with 0 ≤ γ ≤ 1. The time-domain representation of the signals is shown
in figure 2.19. To have a clearer idea, a time-frequency domain representation of
the received signal at the BS is shown in figure 2.20.
1
CRs may likely use wide bands. However, due to the fault of spectrum sensing, a CR-NBI
is produced over some subcarriers allocated to the PU.
84
2.5 An OFDMA Robust CFO Estimator
85
2. Frequency Synchronization and Channel Equalization of OFDMA
Uplink Systems
Secondly, let Y(n) be the observation column vector that stores the real and the
imaginary parts of the PU-OFDMA received signal R(mpu , n):
(
C(n, x(n)) + V(n) + B(n) if n1 ≤ n ≤ n2
Y(n) = (2.35)
C(n, x(n)) + B(n) otherwise
where
h
C(n, x(n)) = Re {f (n, ǫ, hpu )} Im {f (n, ǫ, hpu )} ]T ,
h
B(n) = Re {B(n)} Im {B(n)} ]T is a 2 × 1 white Gaussian noise vector with co-
h
variance matrix (σB2 /2)I2 and V(n) = Re {V (n)} Im {V (n)} ]T is the
2 × 1 zero-mean Gaussian CR-NBI vector with covariance matrix (σV2 /2)I2 . The
CR-NBI signal can be considered as a zero-mean Gaussian vector due to the large
number of subcarriers K cr [Huan 10b].
The SPKF provide “accurate” estimations as long as the received signal is only
disturbed by the AWGN [Pove 10], whereas significant degradations are expected
in the presence of the CR-NBI. To maintain the estimation performance, it is
proposed to jointly detect the CR-NBI and estimate the PU-CFOs.
The approach operates for the nth sample in three steps:
1. estimating the state vector x̂(n|n) by means of SPKF from the observations
{Y(0), Y(1), . . . , Y(K pu − 1)}
86
2.5 An OFDMA Robust CFO Estimator
2
In the following, let us define the SNR expressed in dB as SNR= 10 log( σσ2u ),
B
pu 2
where E{|R(m , n)| } = σu2 is the mean power of the received signal from
the uth user. In addition, we define the signal-to-interference ratio (SIR)
2
expressed in dB as SIR = 10 log( σσ2u ).
V
Remark: If the CR-NBI is assumed to appear at the nth sample, the esti-
mated state vector at time n − 1 and its covariance matrix are kept stored:
x̂nu = x̂(n − 1|n − 1) and Pnu = P(n − 1|n − 1). They will be used as a new
starting point for the CFO and channel estimation once the CR-NBI is considered
to have disappeared. During the CR-NBI perturbation, SPKF provides an esti-
mation of the state vector that is not reliable, but useful to detect the presence
of the CR-NBI.
In the next subsections, two tests to detect the CR-NBI are presented.
40
35
30
25
n1 : beginning n2 : end
||Ỹ(n)||2
15
0
50 100 150 200 250 300 350 400 450 500
Sample Index
87
2. Frequency Synchronization and Channel Equalization of OFDMA
Uplink Systems
To decide whether there is CR-NBI in the OFDMA symbol or not, the following
binary hypothesis is tested:
(
H0 : C(n, x(n)) + V(n) + B(n)
(2.36)
H1 : C(n, x(n)) + B(n)
where λ(n) is a threshold at the nth sample set by the practitioner and cdf (.) is
the cumulative density function (CDF) of (.).
Using ||Ỹ(n)|| the properties and the CDF of the Rayleigh distribution1 , one has:
n o 2
λ(n) = tr PYY (n) ln( ) (2.38)
Pf a
88
2.5 An OFDMA Robust CFO Estimator
BHT inputs:
2
SPKF x̂(n) Measurement Ŷ(n) 2
||Ỹ(n)|| OK OK
Y(n) estimation
||.|| < λ(n) T update
+ -
NOT OK NOT OK
2
T : ||Ỹ(n − 1)|| < λ(n − 1)
update
Figure 2.22: SPKF-BHT algorithm
the mean of the sequence. When using the CUSUM, if the mean increases much,
the cumulative sum also increases. It has been used in various applications such
as information sciences [Huan 11] and biomedicine [Yang 10].
To decide wheter the CR-NBI is present or not, we now suggest using the upper
control limit of the CUSUM test defined as follows:
C + (n) = max(0, C + (n − 1)) + ||Ỹ(n)||2 − tr {PỸỸ (n)}
(2.39)
= max(0, C + (n − 1)) + ỸH (n)Ỹ(n) − E{ỸH (n)Ỹ(n)}
Given (2.39) and fixing a threshold, the CUSUM test makes it possible to detect
when the CR-NBI begins, but does not give information when it ends. See figure
2.23.
Thus, it is proposed to look at the slope of the curve C + (n):
δ(C + (n))
= C + (n) − C + (n − 1) (2.40)
δn
89
2. Frequency Synchronization and Channel Equalization of OFDMA
Uplink Systems
δ(C + (n))
= ||Ỹ(n)||2 − C + (n − 1) > ||Ỹ(n)||2 (2.45)
δn
Otherwise if there is no CR-NBI distortion, given (2.39), (2.40) and (2.43) one
has:
δ(C + (n)) n o
= ||Ỹ(n)||2 − tr PYY (n) + C + (n − 1)
δn (2.46)
≈ C + (n − 1)
δ(C + (n))
Therefore, one has to compare δn
to a threshold in order to detect the end of
the CR-NBI in the PU-OFDMA symbol. The CR-NBI is assumed to be present
δ(C + (n))
provided that δn
is higher than the threshold.
Concerning the choice of the threshold, we suggest the one defined in (2.38).
2000
1500
n2 : end
n1 : beginning of the CR-NBI
C + (n)
of the CR-NBI
1000
500
0
50 100 150 200 250 300 350 400 450 500
Sample Index
90
2.5 An OFDMA Robust CFO Estimator
In this subsection, the goal is to reduce the gap between the CR-NBI probability
of detection Pcr 2
d = P (||Ỹ(n)|| > λ(n)|H0 ) and the PU probability of detection
Ppu 2
d = P (||Ỹ(n)|| < λ(n)|H1 ), in order to obtain better estimation performance
when using the SPKF-BHT and the SPKF-CT. If the CR-NBI probability of
pu
detection Pcr
d increases, the PU probability of detection Pd decreases. Therefore,
there are less samples available to perform the CFO estimation. This hence
may lead to poor estimation of the CFO. The PU probability of detection Pdpu
determines the values that are not disturbed by the CR-NBI and that are going
to be used in the estimation task.
When looking at figure 2.21, the values of ||Ỹ(n)||2 are not always higher than the
threshold between n1 and n2 sample indexes. In addition, the values of ||Ỹ(n)||2
may sometimes exceed the threshold outside the interval [n1 , n2 ]. To improve the
CR-NBI and the PU detections when using the BHT, we propose to take into
account the evolution of ||Ỹ(n)||2 over several samples. Therefore, we suggest
defining a new residual value ||Ỹmean (p)||2 as follows:
1 β(p+1)−1
X
||Ỹmean (p)||2 =
||Ỹ(q)||2 (2.47)
j pu k β q=βp
where p ∈ {0, . . . P − 1}, P = Kβ and ⌊.⌋ is the integer part of {.}.
j k
K pu K pu K pu
If β
6= β
, i.e. β
∈
/ Z, the last sample of the new residual value is:
KX −1 pu
2 1
||Ỹmean (P)|| = pu ||Ỹ(q)||2 (2.48)
K − βP q=βP
If the value of ||Ỹmean (p)||2 exceeds the value of the threshold, the algorithm
decides that the CR-NBI is present during the corresponding samples. See algo-
rithm 3.
Based on the same idea, to improve the detection performance when using the
j pu
k
CUSUM test, we propose to keep only P − P − Kβ (i.e P or P + 1) samples
as follows:
+
Cmean (p) = C + (βp − 1) (2.49)
+ K pu
with p ∈ {1, . . . P − 1}, Cmean (0) = C + (0) and if β
∈
/ Z then
+
Cmean (P) = C + (K pu − 1). See figure 2.24.
91
2. Frequency Synchronization and Channel Equalization of OFDMA
Uplink Systems
+ K pu
Figure 2.24: Cmean (p + 1) and C + (n) with β ∈Z
+
δ(Cmean (p)) C + (p)−C + (p−1)
If the value of δp
= β
exceeds the value of the threshold, the
algorithm decides that the CR-NBI is present in the corresponding samples.
In the next subsection, simulation results show the relevance of our algorithms
and their improvements to jointly detect the CR-NBI and estimate the PU-CFOs.
92
2.5 An OFDMA Robust CFO Estimator
mation bits. The carrier frequency is at fc = 2.6GHz and the bandwidth is set
to W pu = 5MHz. A CR-NBI affects K nbi = 40 subcarriers of the PU-OFDMA
system bandwidth. It is supposed a transmission over a normalized Rayleigh
slow-fading frequency-selective channel composed of Lpu
u = 7 ∀u multipaths.
For very low values of SIR, as ||Ỹ(n)||2 >> 0 during the CR-NBI presence, its
detection can be easy. In addition for high values of SIR, i.e SIR> 0dB the dif-
ference between σB2 and (σB2 + σV2 ) may be negligible and would not affect the
Kalman filter estimation performance. Therefore the case when SIR= −10 dB is
studied. The CR-NBI can begin anywhere in the PU-OFDMA received symbol.
For comparison reasons with [More 08], the CR-NBI is supposed to affect only
the subcarriers assigned to one PU.
SPKF-BHT
SPKF-CT
0.8
CR-NBI Probability of detection
0.75
0.7
0.65
1 2 3 4 5 6 7 8
β
Figure 2.25: CR-NBI detection, SNR= 0 dB, SIR= −10 dB, Pf a = 0.05 and
γ = 0.5
93
2. Frequency Synchronization and Channel Equalization of OFDMA
Uplink Systems
This subsection presents how to select a “well suited” value of β, which is the
parameter presented in subsection 2.5.2.3, to reduce the gap between the CR-NBI
pu
probability of detection Pcr
d and the PU probability of detection Pd .
Firstly the impact of β in the CFO estimation performance is shown. The SNR
is set to 0 dB and the Pf a is set to 0.05.
Let us consider that the CR-NBI disturbs half of the PU-OFDMA symbol
(γ = 0.5). Figure 2.25 shows the CR-NBI probability of detection Pcr
d for dif-
ferent values of β. We clearly see that if β increases, the CR-NBI probability of
detection Pcr
d increases and the CFO estimation should normally be improved.
However, when looking at figure 2.26 we can see that the PU probability of de-
tection Ppu
d does not increase with the value of β as the CR-NBI probability of
detection Pcr
d . For both proposed algorithms, there are “well suited” values of β
in terms of PU probability of detection. For the SPKF-BHT, the “well suited”
value of β is between 4 and 6 whereas the “well suited” value of β is between 2
and 3 when using the SPKF-CT. The CR-NBI probability of detection Pcr
d and
the PU probability of detection Ppu
d has been increased, then better estimation
performance should be obtained by using this optimal β parameter.
1
SPKF-BHT
SPKF-CT
0.995
PU Probability of detection
0.99
0.985
0.98
1 2 3 4 5 6 7 8
β
Figure 2.26: PU detection, SNR= 0 dB, SIR= −10 dB, Pf a = 0.05 and γ = 0.5
94
2.5 An OFDMA Robust CFO Estimator
Tables 2.11 and 2.12 show the CFO estimation performance in terms of MMSE
for different values of β and duration of the CR-NBI over the symbol. It should
be noted that the duration of the CR-NBI over the PU-OFDMA symbol does not
change the influence of β in the algorithm. For different durations of the CR-NBI
over the symbol, we can see that better performance in terms of CFO estimation
are obtained for values β between 4 and 6 for the SPKF-BHT. In addition, it can
be seen that the “well suited” value of β in terms of probability of detection is
between 2 and 3 for the SPKF-CT.
Table 2.12: SPKF-CT CFO estimation performance for different values of γ and
β, Pf a = 0.05
95
2. Frequency Synchronization and Channel Equalization of OFDMA
Uplink Systems
0.85
0.8
SIR = −10 dB
SNR = 0 dB
0.75
CR-NBI probability of detection
0.7
0.65
0.6
0.55
Improved SPKF-CT β = 1
Improved SPKF-BHT β = 1
0.5
Improved SPKF-CT β = 2
Improved SPKF-BHT β = 6
0.45
Figure 2.27: CR-NBI detection, SNR= 0 dB, SIR= −10 dB and γ = 0.5
Figure 2.29 shows the estimation performance in terms of MMSE for different
durations of the CR-NBI over the PU-OFDMA symbol. The SNR is again set to
96
2.5 An OFDMA Robust CFO Estimator
0.995
0.99
PU probability of detection
0.985
0.98
SIR = −10 dB
SNR = 0 dB
0.975
0.97
Improved SPKF-CT β = 1
Improved SPKF-BHT β = 1
0.965
Improved SPKF-CT β = 2
Improved SPKF-BHT β = 6
0.96
0.955
0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09
Probability of false alarm
MMLE[11]
−1
10 SPKF-BHT β = 1
SPKF-BHT β = 6
MMSE
−2
10
−3
10
−4
10
0.3 0.4 0.5 0.6 0.7 0.8
Pf a=0.05 γ
SIR= −10dB
SNR= 0 dB
MMLE[11]
−1 SPKF-CT β = 1
10
SPKF-CT β = 2
MMSE
−2
10
−3
10
−4
10
0.3 0.4 0.5 0.6 0.7 0.8
γ
Figure 2.29: Robust CFO estimation performance, SIR= −10 dB, SNR= 0 dB
and Pf a = 0.05
0 dB. We clearly see that the proposed algorithms provide better CFO estimation
performance since γ < 0.6 in comparison with [More 08].
Figure 2.30 shows the CFO estimation performance in terms of MMSE. The CR-
NBI disturbs half of the PU-OFDMA symbol (γ = 0.5).
97
2. Frequency Synchronization and Channel Equalization of OFDMA
Uplink Systems
Normal KF
SIR = −10 dB
Pf a = 0.05
MMLE[12]
−3
10
MMSE
−4
10
SPKF-CT β = 1
SPKF-BHT β = 1
SPKF-CT β = 2
−5 SPKF-BHT β = 6
10
0 2 4 6 8 10 12 14 16 18
SNR
Figure 2.30: Robust CFO estimation performance, SIR= −10 dB, Pf a = 0.05
and γ = 0.5
Table 2.13 shows the estimation performance in terms of MMSE for different
values of SIR. It should be noted that the proposed algorithms are not disturbed
by the power of the CR-NBI.
MMSE
SNR= 0 dB
SIR −30 dB −20 dB −10 dB 0 dB
SPKF 0.43520 0.08150 0.003557 0.00172
SPKF-CT β = 1 0.0028500
SPKF-BT β = 1 0.0012670
SPKF-CT β = 2 0.0008217
SPKF-BT β = 6 0.0003958
Since ĥpu is determined by the values ǫ̂, one can conclude that the proposed
estimators achieve accurate channel estimation. Figure 2.31 confirms the channel
98
2.6 Conclusions
−3
x 10
6.5
SPKF-CT β = 1
6 SPKF-BHT β = 1
SPKF-CT β = 2
5.5 SPKF-BHT β = 6
4.5
MMSE
3.5
2.5
SIR = −10 dB
Pf a = 0.05
2
1.5
0 2 4 6 8 10 12 14 16 18
SNR
Figure 2.31: Robust channel estimation performance, SIR= −10 dB, Pf a = 0.05
and γ = 0.5
2.6 Conclusions
In this chapter we proposed to jointly estimate the CFOs and the channels for
an OFDMA uplink system by using a recursive approach. On the one hand, the
EKF, the SOEKF, the IEKF and the SPKF including the UKF and the CDKF
were considered. On the other hand, to relax the assumptions on the model and
measurement noise, the “extended H∞ filter” and the “unscented H∞ filter” filter
were studied.
Then, the above techniques have been combined with an MMSE-SD to obtain a
non-pilot aided CFO estimator for an OFDMA uplink. In addition, as the SPKF
give the best compromise between computational complexity and estimation per-
formance, it has been chosen to be combined with two statistical tests: the BHT
and the CUSUM test, to design two OFDMA robust CFO estimators.
Both SPKF-based algorithms are used to estimate the CFOs and the channels of
different users in a PU-OFDMA system disturbed by a localized in time CR-NBI.
99
2. Frequency Synchronization and Channel Equalization of OFDMA
Uplink Systems
The proposed estimators operate in the time domain unlike other approaches that
estimate the CFO in presence of an NBI in the frequency domain [More 08].
By using a SPKF and without increasing the computational complexity, we take
advantage of the innovation to detect a sharp variation of the measurement noise
covariance matrix when there is CR-NBI. The SPKF enable the state vector to
be estimated and the innovation energy to be calculated. Then, a test detects
the variation of the measurement-noise covariance matrix by comparing the in-
novation energy with a threshold set by the practitioner.
When using the BHT, better performance are obtained in terms of the CR-NBI
probability of detection. However, when using the CUSUM, better performance
in terms of PU probability of detection are obtained. To reduce the gap between
the CR-NBI probability of detection and the PU probability of detection, a tech-
nique is then proposed. It consists in improving both probabilities of detection
by watching the evolution of the innovation energy.
Our method outperforms other existing methods [More 08]. However, the pro-
posed estimators are sensitive to the duration in time of the CR-NBI.
The next chapter still deals with the constraints addressed in this chapter, i.e.
frequency synchronization and channel equalization, but in the context of the
uplink OFDM-IDMA systems.
100
Chapter 3
Frequency Synchronization and Channel
Equalization of OFDM-IDMA Uplink
Systems
Contents
101
3. Frequency Synchronization and Channel Equalization of
OFDM-IDMA Uplink Systems
3.1 Introduction
The OFDM-IDMA technique initially proposed by Mahafeno et al. [Maha 06]
has been introduced in chapter 1.
Firstly, OFDM-IDMA can be seen as a special case of multicarrier CDMA
(MC-CDMA). IDMA inherits many advantages from CDMA, in particular, its
robustness against different kinds of interferences. In addition, the autocorrela-
tion properties of the IDMA transmitted signal make it possible to take advantage
of the multipath fading channel diversity1 . The key principle of IDMA is that
the spreading code is the same for all users, unlike CDMA2 . The user distinction
is done by using different interleavers. They are generated independently and
randomly, allowing BER to be near the theoretical values of single-user systems.
Secondly, OFDM-IDMA can be seen as a special case of OFDMA. OFDM-IDMA
inherits the robustness to ISI from OFDMA. In addition, the multicarrier scheme
allows the frequency allocation algorithms to improve the spectral efficiency by
exploiting the spectral diversity. The main difference between these techniques
is that the entire bandwidth can be allocated to a single user for all the OFDM-
IDMA frame duration. This may lead to a very high single-user capacity. Indeed,
if we assume an OFDMA system with U users and K subcarriers all used, each
K
user can transmit U
symbols at the same time. Now if one also consider an
OFDM-IDMA system with a spreading factor3 S = log2 (U ), each user can trans-
K K
mit S
= log2 (U )
symbols at the same time. We clearly see that log2 (U ) < U ; thus
the users in the OFDM-IDMA system can transmit more symbols at the same
time than in an OFDMA system.
Furthermore, an OFDM-IDMA scheme may improve the duty cycle of a single-
user transmission4 in comparison with an OFDMA system. Indeed, in an OFDMA
1
In chapter 1 in section 1.3, we show that the received signal is a superposition of several
delayed and attenuated copies of the transmitted signal, due to the multipath fading channel.
This diversity may be exploited by a combination of the different copies using a rake receiver
[Proa 95].
2
In chapter 1 in section 1.4.3.1, we present that the ENC of the OFDM-IDMA transmitter
incorporates a spreading code.
3
S = log2 (U ) is the lowest spreading factor to allow U users transmit at the same time.
4
The duty cycle of a single-user transmission is defined as the ratio between the duration
of one user frame and the time between two successive frames.
102
3.1 Introduction
system, the maximum number of users allowed by the system to transmit at the
same time is U = K. In that case, each user transmits one symbol at the same
time. In an OFDM-IDMA system, the maximum number of users allowed by
the system to transmit at the same time is U = 2S . In addition, to transmit
one symbol at the same time, it is necessary that K = S. Then, the maximum
number of users allowed by an OFDM-IDMA system to transmit one symbol at
the same time is U = 2K . We clearly see that 2K > K; the OFDM-IDMA enables
more users to perform a transmission at the same time.
Today, the OFDM-IDMA scheme is not used as the multiple access technique by
mobile wireless standards, but for the above reasons, it can be seen as a promising
technique for the future mobile wireless standards.
To become one of the multiple access schemes of future mobile standards, several
constraints have to be addressed, such as the channel equalization. Nevertheless,
few papers deal with this issue in IDMA systems. Thus, Zhou et al. develop an
iterative channel estimation process that uses pilot subcarriers [Zhou 07]. The
authors propose two different approaches: the first one is based on an unbiased LS
estimator whereas the second is based on the maximum ratio combining (MRC)
method. In [Song 09], the authors propose an EM estimation method to obtain
the CSIs. More specifically, in [Suya 08], Suyama et al. insert a least mean
square (LMS) estimator in the iterative process of an OFDM-IDMA receiver. In
addition, the authors in [Rehm 08] present a low-complexity uplink channel es-
timation. By taking advantage of the channel coding, the transmitter creates a
training sequence in order to estimate the channel through a linear detection at
the receiver.
However, in the above approaches, a perfectly frequency synchronized system
is considered. In wireless communication systems with multicarrier modulation,
this assumption is usually not true. Therefore, a CFO estimation/correction has
to be considered at the receiver.
Several approaches of the frequency synchronization in multicarrier systems have
been proposed in the last years. More particularly, in MC-CDMA systems, Chien
et al. present a blind recursive estimation of the CFOs by using an EM algo-
rithm [Chie 07]. In [Thia 07], the authors take advantage of the knowledge of the
103
3. Frequency Synchronization and Channel Equalization of
OFDM-IDMA Uplink Systems
second-order statistics of the received signal. They suggest a two-stage CFO esti-
mation. The first one involves a grid search to obtain the initialization of the CFO
estimations, whereas the second consists of a CFO estimation based on the LMS
algorithm. In [Chen 10], the CFOs are estimated by searching the peaks of the
correlation function of the received signal. In addition, the reader can look at the
frequency synchronization techniques for OFDMA systems, [Pun 04a], [Pun 06]
and [Sezg 08] presented in chapter 2.
In the above methods, the CFO correction can be performed before or after the
CP-OFDM demodulation. The CFO correction methods can take advantage of
the orthogonal codes in MC-CDMA systems, whereas the CFO correction meth-
ods usually utilize the CAS1 in OFDMA systems. The CFO estimation/correction
issue in the particular OFDM-IDMA scheme becomes a challenging task2 . In
[Cao 04a], Cao et al. propose to correct the CFO in an OFDMA uplink system
after the CP-OFDMA demodulation through a linear detection3 , without taking
advantage of the CAS. This CFO correction method can be completely applied
to OFDM-IDMA systems. However, CFOs and CSIs are assumed to be perfectly
known. Although frequency synchronization has been addressed in MC-CDMA
and OFDMA, it has never been studied in an OFDM-IDMA system to our knowl-
edge.
Given the above considerations, our contribution in this chapter is twofold:
104
3.2 OFDM-IDMA Modified Receiver
where we recall that hu (l) is the lth coefficient of the CIR, ǫu is the normalized
CFO, Xu (n) is the nth sample of the transmitted OFDM-IDMA symbol and B(n)
is a zero-mean AWGN with variance σB2 .
Then, the received OFDM-IDMA symbol can be expressed as follows:
∆
R = [R(0), R(1), . . . , R(n), . . . , R(K − 1)]T
U
X (3.2)
= Eu hu Xu +B
u=1
| {z }
Ru
nh io
where Eu = diag 1, ej2πǫu /K , . . . , ej2π(K−1)ǫu /K is the CFO matrix, hu is the
channel circulant matrix defined in (1.31), Xu is the transmitted OFDM-IDMA
symbol defined in (1.21) and B is a zero-mean AWGN vector with covariance
matrix σB2 IK (See chapter 1 in section 1.4.3.2).
Given (1.20), it should be noted that (3.1) can be rewritten as follows:
U K−1
X ǫu n X k′ n
R(n) = ej2π K Hu (k ′ )Su (k ′ )ej2π K + B(n) 0 ≤ n ≤ K − 1 (3.3)
u=1 k′ =0
105
3. Frequency Synchronization and Channel Equalization of
OFDM-IDMA Uplink Systems
PL−1 k′ l
where Hu (k ′ ) = l=0 hu (l)e−j2π K is the channel frequency response associated
to the k ′ th subcarrier and Su (k ′ ) is the symbol transmitted over the k ′ th subcar-
rier.
Given (3.3), although the effect of the CFO of one user or another can be elim-
inated, they cannot be suppressed all at the same time before the CP-OFDM
demodulation. Whatever the strategy chosen by the receiver to compensate the
CFO, there are still some interferences. Therefore, we suggest introducing the
CFO correction after the CP-OFDM demodulation.
After the CP-OFDM demodulation, the received symbol r(k) can be expressed
as follows:
K−1
X kn
r(k) = R(n)e−j2π K + b(k)
n=0
K−1
X X U K−1
X (k′ −k)n
j2π ǫK
un
= e Hu (k ′ )Su (k ′ )ej2π K + b(k)
n=0 u=1 k′ =0
U K−1
(3.4)
X X K−1
X (k′ −k+ǫu )n
j2π
= e K Hu (k ′ )Su (k ′ ) + b(k)
u=1 k′ =0
n=0
XU K−1
X
= F(k ′ − k + ǫu )Hu (k ′ )Su (k ′ ) + b(k)
u=1 k′ =0
where1
sin(xπ) jπ x(K−1)
F(x) = e K (3.5)
sin( xπ
K
)
and
K−1
X kn
b(k) = B(n)e−j2π K (3.6)
n=0
U
X U K−1
X X
r(k) = F(ǫu )Hu (k)Su (k) + F(ǫu + k ′ − k)Hu (k ′ )Su (k ′ ) +b(k)
u=1 u=1 k′ =0 (3.7)
k′ 6=k
| {z }
ICI(k)
1
PK−1 nx 1−ej2πx sin(xπ) jπ x(K−1)
Let us recall that n=0 ej2π K = j2π x = sin( xπ e K .
1−e K K )
106
3.2 OFDM-IDMA Modified Receiver
If we look at the received signal from the uth user point of view, (3.7) can be
expressed as a summation of a useful term and an interference term as follows:
U
X
r(k) = F(ǫu )Hu (k)Su (k) + F(ǫu′ )Hu′ (k)Su′ (k) + ICI(k) + b(k)
u′ =1
u′ 6=u
U
X (eq) (3.8)
= Hu(eq) (k)Su (k) + Hu′ (k)Su′ (k) + ICI(k) + b(k)
u′ =1
u′ 6=u
where
Hu(eq) (k) = F(ǫu )Hu (k) ∀u ∈ {1, . . . , U } (3.9)
U
X (eq)
ςu (k) = Hu′ (k)Su′ (k) + ICI(k) + b(k) (3.10)
u′ =1
u′ 6=u
2
and ICI(k) is a zero-mean ICI with variance σICI = E{|ICI(ǫ, k)|2 }, caused by
the CFOs [Moos 94]. This interference can be considered as a Gaussian vector due
to the large number of subcarriers [Huan 10b]. It should be noted that the ICI
is included in ςu (k), the characteristics of which are used by the IDMA receiver
to estimate the transmitted bits (See chapter 1 in section 1.4.3.2).
In the following, we will show that the transmitted bits can be recovered by taking
into account the ICI characteristics in the IDMA receiver. In addition, as r(k)
(eq)
in (3.8) is a function of ǫu′ and Hu′ (k), it is necessary to know the CFOs and
the CSIs values at the IDMA receiver. In the next subsection, the estimations of
these vectors are presented.
ǫ = [ǫ1 , ǫ2 , . . . , ǫu , . . . , ǫU ] (3.11)
h = [h1 , h2 , . . . , hu , . . . , hU ] (3.12)
107
3. Frequency Synchronization and Channel Equalization of
OFDM-IDMA Uplink Systems
Let Y be the 2 × K observation matrix that stores the real and the imaginary
parts of the received OFDM-IDMA symbol R:
Y(ǫ, h) = C(ǫ, h) + B
(3.14)
= [Y(0, ǫ, h), Y(1, ǫ, h), . . . , Y(n, ǫ, h), . . . , Y(K − 1, ǫ, h)]
where C(ǫ, h) is a 2 × K matrix that stores the real an the imaginary parts of
the contribution of all users to the received OFDM-IDMA symbol:
nP oT
U
Re u=1 Ru
C(ǫ, h) =
nP
oT
Im U
u=1 Ru (3.15)
= [C(0, ǫ, h), C(1, ǫ, h), . . . , C(n, ǫ, h), . . . , C(K − 1, ǫ, h)]
and " #
Re{B(n)}T
B=
Im{B(n)}T (3.16)
= [B(0), B(1), . . . , B(K − 1)]
2
σB
B(n) is an AWGN vector with zero-mean and covariance matrix I.
2 2
Now let us introduce the state-space representation of the system to estimate the
CFO and the channel:
State equation:
108
3.2 OFDM-IDMA Modified Receiver
Measurement equation:
where w(n) is an AWGN matrix with zero-mean and covariance matrix σw2 IU . It
should be noted that σw2 is very low and can even be equal to 0.
The SPKF detailed in appendix E is used to estimate the CFOs and the channels,
after the time synchronization and CP removal.Then, the OFDM demodulation
is performed without any CFO correction and the resulting samples are used in
the modified IDMA receiver presented in the next subsection.
109
3. Frequency Synchronization and Channel Equalization of
OFDM-IDMA Uplink Systems
Then, given (3.20) and by replacing Hu(eq) (k) by its estimate, one has:
U
X
E(r(k)) = Ĥu(eq) (k)µu (k) ∀u, k
u=1
U
(3.21)
X
V ar(r(k)) = σZ2 + (Ĥu(eq) (k))2 υu (k) ∀u, k
u=1
with σZ2 = σICI
2
+ σB2 .
The outputs of the modified ESE are inserted in the user decoding process (See
1st user
gESE
1 (k) estimated
bits
g1ESE(k)
Ĥu
(eq) 2 + σ2
σz2 = σB ICI
ĥ, ǫ̂
R SPKF initialization step
Figure 3.1: Proposed OFDM-IDMA uplink receiver, where the user decoding is
the same as in figure 1.20
110
3.2 OFDM-IDMA Modified Receiver
chapter 1 in section 1.4.3.2). Then, at the ith iteration for i < Iidma , the outputs
of the decoding process are the values of guESE (k), that are reinserted in the ESE.
The ESE process is again performed. Finally, the decoding process produces
hard decisions to obtain the estimations of the transmitted bits during the final
iteration Iidma . See algorithm 4, where the modified ESE is summarized. Figure
3.1 shows the modified OFDM-IDMA uplink receiver.
111
3. Frequency Synchronization and Channel Equalization of
OFDM-IDMA Uplink Systems
random variables with a normal distribution with zero-mean and variance 0.1.
The estimations are obtained by using control data of one OFDM-IDMA symbol.
Figure 3.2 shows the convergence speed of the SPKF for a different number of
users in the system. We can see that the number of users in the system does
not really affect the estimation performance of the SPKF. This is a key point
in OFDM-IDMA systems that allow a large number of users to transmit at the
same time. However, the computational complexity of the SPKF increases with
the number of users in the system. Then, we fix the number of user to U = 4.
Figure 3.3 shows the CFO estimation performance in terms of MMSE of one user
in the system. The SPKF allows the CFO and the channel of the users to be
estimated. The MMSE for an AWGN channel and for channels with different
length are presented. The channel estimation for different channel lengths give
the same estimation performance, around 10−2 for a SNR= 0dB. Previous studies
in [Ping 02a] and [Maha 06] show the performance of the IDMA receiver for dif-
ferent number of multipaths, for different number of users and different number
of iterations. It has been already proved that the performance of the IDMA re-
ceiver in terms of BER are better when the number of multipaths and iterations
increases and also when the number of users decreases [Ping 02a], [Maha 06]. For
0.2
U=4
0.15 U=8
U=16
CFO value
0.1
0.05
−0.05
0 100 200 300 400 500 600 700 800 900 1000
Sample Index
112
3.2 OFDM-IDMA Modified Receiver
SPKF Lu = 3
10
−2 SPKF Lu = 1
SPKF AWGN
MMSE
−3
10
−4
10
0 2 4 6 8 10 12 14 16 18
SNR (dB)
this reason, we keep the number of users equal to U = 4, and the number of iter-
ations to 3. The following results are the average of 5000 OFDM-IDMA symbols.
Figure 3.4 shows the BER for a transmission over a normalized Rayleigh slow-
fading frequency selective channel channel composed of ∀u Lu = 1 and Lu = 3
multipaths. The efficiency of the algorithm is confirmed. Indeed, we can see
that for Lu = 1 and a BER = 10−1 the difference between the results of the
proposed scheme and the perfect case is lower than 1dB. In addition, we can
notice that for Lu = 3 and a BER= 10−2 the difference between the results of the
proposed scheme and the perfect case is less than 1 dB. In both cases (Lu = 1 and
Lu = 3), that the proposed scheme gives a gain of around 3dB for BER= 10−1
in comparison with the linear detection method presented in [Cao 04a]. Figure
3.5 shows the BER for a transmission over a normalized Rayleigh slow-fading
frequency selective channel channel composed of ∀u Lu = 5 multipaths. We can
see that the difference between the results of the proposed scheme and the perfect
case is lower than 1 dB. In addition, the proposed scheme gives a gain of around
2dB for BER= 10−2 in comparison with the linear detection method presented
in [Cao 04a]. It should be noticed that the performance for Lu = 5 are better
than for Lu = 3 and Lu = 1. The IDMA receiver takes advantage of the diversity
introduced by the propagation channel.
113
3. Frequency Synchronization and Channel Equalization of
OFDM-IDMA Uplink Systems
−1 Lu = 1
10
Lu = 3
−2
10
BER
−3
10
No CFO - Channel Perfectly estimated
SPKF-CFO estimated - Channel perfectly estimated
SPKF-CFO/channel estimated
Linear detection-CFO and Channel perfectly estimated
−4
10
0 2 4 6 8 10 12 14
SNR (dB)
0
10
−1
10
BER
−2
10
−3
10
No CFO - Channel Perfectly estimated
SPKF-CFO estimated - Channel perfectly estimated
SPKF-CFO/channel estimated
Linear detection-CFO and Channel perfectly estimated
−4
10
0 2 4 6 8 10 12 14
SNR (dB)
114
3.3 Space-Time Block Code OFDM-IDMA Receiver
3.2.4 Conclusions
In this section, a modified OFDM-IDMA receiver is proposed. The receiver is
robust against ISI due to the OFDM modulation and it is robust against the
MAI due to the IDMA scheme. Unlike common approaches, the IDMA scheme
allows the CFO correction step to be skipped. After the OFDM demodulation,
the CFOs and the channels are jointly estimated by using a SPKF. Then, the
OFDM demodulation is performed without any CFO correction. Finally, the
modified IDMA receiver accurately estimates the transmitted bits of the uth
user, by taking into account the estimations of the SPKF and the interference
produced by the CFO. The proposed IDMA receiver is able to cancel the ICI
produced by the CFOs , as well as the MAI produced by the other users in the
system. Simulation results clearly show the efficiency of the proposed OFDM-
IDMA receiver. In addition, we show that the proposed receiver outperforms
existent CFO correction methods such as [Cao 04a].
In the next section, we show that the proposed scheme can be extended to the
multiple-input multiple-output (MIMO) case.
115
3. Frequency Synchronization and Channel Equalization of
OFDM-IDMA Uplink Systems
Su,1 Xu,1
CP-OFDM
bits to be modulation
transmitted Su
IDMA Alamouti
transmitter encoder
Xu,2
CP-OFDM
Su,2 modulation
Figure 3.6: OFDM-IDMA-STBC transmitter for the uth user, the IDMA trans-
mitter is the same as in figure 1.19, Su,p are the IDMA symbols and Xu,p are the
CP-OFDM-IDMA symbols
∆
Su = [Su (0), Su (1), . . . , Su (K − 1)]T (3.23)
116
3.3 Space-Time Block Code OFDM-IDMA Receiver
∆
Xu,1 = FH Su,1 = FH [Su (0), −Su∗ (1), . . . Su (K − 2), −Su∗ (K − 1)]T
∆
(3.24)
Xu,2 = FH Su,2 = FH [Su (1), Su∗ (0), . . . Su (K − 1), Su∗ (K − 2)]T
hu,p = [hu,p (0), hu,p (1), . . . , hu,p (l), . . . , hu,p (Lu,p − 1)] (3.25)
ǫ = [ǫ1 , ǫ2 , . . . , ǫu , . . . , ǫU ]
h = [h1 , h2 , . . . , hu , . . . , hS ] (3.26)
∆
Ru = [Ru (0), . . . , Ru (n) . . . , Ru (K − 1)]
= Eu,1 hu,1 Xu,1 + Eu,2 hXu,2 (3.27)
h i
where Eu,p = diag 1, ej2πǫu,p /K , . . . , ej2π(K−1)ǫu,p /K , ǫu,p is the normalized CFO
and hu,p is the channel circulant matrix defined as follows (See chapter 1 in
117
3. Frequency Synchronization and Channel Equalization of
OFDM-IDMA Uplink Systems
section 1.4.2.2):
hu,p (0) 0
. . . hu,p (Lu,p − 1) . . . hu,p (1)
... ...
hu,p (1) hu (0) 0 hu,p (2)
.. ... ... ... ... ..
. .
hu,p = hu,p (Lu − 1) hu,p (Lu,p − 2) . . .
0
...
0
(3.28)
... ... ..
0 hu,p (Lu − 1) hu,p (0) .
.
.. . .. . .. ... ... ..
.
0 0 . . . hu,p (Lu,p − 2) . . . hu,p (0)
where B(n) is a zero-mean AWGN noise vector with covariance matrix IK σB2 . As
in section 3.2, the CP-OFDM is performed without any CFO correction. Then,
after some manipulations (See section 3.2), the received samples r(k) and r(k +1)
can be expressed can be expressed as:
U
X (eq) (eq)
r(k) = {Hu,1 (k)Su (k) + Hu,2 (k)Su (k + 1)} + ICI(k) + b(k)
u=1
XU
(eq) (eq)
r(k + 1) = {−Hu,1 (k + 1)Su∗ (k + 1) + Hu,2 (k + 1)Su∗ (k)}
u=1
+ ICI(k + 1) + b(k + 1)
(eq)
where Hu,p (k) = F(ǫu,p )Hu,p (k) for p ∈ {1, 2} and ICI(k) is a zero-mean Gaus-
2
sian ICI with variance σICI = E{|ICI(ǫ, k)|2 } caused by all the CFOs. In addi-
tion,
PK−1 j2πnk
b(k) = n=0 B(n)e− K , F(ǫu,p ) is defined in (3.5) and
Ls,p −1
X j2πlk
Hu,p (k) = hu,p (l)e− K forp ∈ {1, 2} (3.30)
l=0
In the following, we show that the transmitted bits can be recovered by taking
into account the ICI characteristics in the IDMA receiver. As it is necessary to
118
3.3 Space-Time Block Code OFDM-IDMA Receiver
know the values of CFOs and the CSIs at the IDMA receiver, the estimation of
these vectors is presented in the next subsection.
Now let us introduce the state-space representation1 of the system estimate the
CFO and the channel:
State equation:
Measurement equation:
where w(n) is an AWGN matrix with zero-mean and covariance matrix σw2 IU ,
Y(n) is the observation 2 × 1 vector that stores the real and the imaginary parts
of the STBC-OFDM-IDMA symbol R(n), B(n) is an AWGN matrix with zero-
mean and covariance matrix σB2 I2 .
The SPKF detailed in appendix E is used to estimate the CFOs and the chan-
nels. After the time synchronization and CP removing, the CFOs and channels
are estimated. Then, the OFDM demodulation is performed without any CFO
correction and the obtained samples are inserted in the MRC receiver presented
in the next subsection.
119
3. Frequency Synchronization and Channel Equalization of
OFDM-IDMA Uplink Systems
and Ŝu (k + 1) are the decisions obtained by the MRC of two successive received
symbols r(k) and r(k + 1), and the channel frequency responses Ĥu(eq) [Alam 98]:
U
X (eq)∗ (eq) (eq) (eq)∗
Ŝu (k) = {Ĥu,1 (k)Ĥu′ ,1 (k) + Ĥu,2 (k + 1)Ĥu′ ,2 (k + 1)}Su (k)
u′ =1
XU
(eq)∗ (eq) (eq) (eq)∗
+ {Ĥu,1 (k)Ĥu′ ,2 (k) − Ĥu,2 (k + 1)Ĥu′ ,1 (k + 1)}Su (k + 1)
u′ =1
(eq)∗ (eq)
+ Ĥu,1 (k){ICI(ǫ, k) + b(k)} + Ĥu,2 (k + 1){ICI ∗ (ǫ, k + 1) + b∗ (k + 1)}
(3.34)
U
X (eq) (eq)∗ (eq)∗ (eq)
Ŝu (k + 1) = {Ĥu,1 (k + 1)Ĥu′ ,1 (k + 1) + Ĥu,2 (k)Ĥu′ ,2 (k)}Su (k + 1)
u′ =1
XU
(eq)∗ (eq) (eq) (eq)∗
+ {Ĥu,2 (k)Ĥu′ ,1 (k) − Ĥu,1 (k + 1)Ĥu′ ,2 (k + 1)}Su (k)
u′ =1
(eq)∗
+ Ĥu,2 (k){ICI(ǫ, k) + b(k)}
(eq)
+ Ĥu,1 (k + 1){ICI ∗ (ǫ, k + 1) + b∗ (k + 1)}∗
(3.35)
Then, the values Ŝu (k) and Ŝu (k + 1) are inserted in the ESE presented in the
next subsection.
Then, given the outputs of the MRC Ŝ(k) and Ŝ(k + 1), the values of Ĥu(eq) (k),
the additive-noise variance σB2 , the ICI variance σICI
2
and guESE (k), the IDMA
reception is performed.
It should be noted that (3.34) and (3.35) can be rewritten as:
(eq) (eq)
Ŝu (k) = {||Ĥu,1 (k)||2 + ||Ĥu,2 (k + 1)||2 }Su (k) + ςu (k)
(eq) (eq)
Ŝu (k + 1) = {||Ĥu,1 (k + 1)||2 + ||Ĥu,2 (k)||2 }Su (k + 1)
+ςu (k + 1) (3.36)
120
3.3 Space-Time Block Code OFDM-IDMA Receiver
and the distortions ςu (k) and ςu (k +1) contained in the output values of the MRC
can be expressed as:
U
X (eq)∗ (eq) (eq) (eq)∗
ςu (k) = {Ĥu,1 (k)Ĥu′ ,1 (k) + Ĥu,2 (k + 1)Hu′ ,2 (k + 1)}Su (k)
u′ =1
u′ 6=u
U
X (eq)∗ (eq) (eq) (eq)∗
+ {Ĥu,1 (k)Hu′ ,2 (k) − Ĥu,2 (k + 1)Hu′ ,1 (k + 1)}Su (k + 1) (3.37)
u′ =1
u′ 6=u
(eq)∗
+ Ĥu,1 (k){ICI(ǫ, k) + b(k)}
(eq)
+ Ĥu,2 (k + 1){ICI ∗ (ǫ, k + 1) + b∗ (k + 1)}
U
X (eq) (eq)∗ (eq)∗ (eq)
ςu (k + 1) = {Ĥu,1 (k + 1)Hu′ ,1 (k + 1) + Ĥu,2 (k)Ĥu′ ,2 (k)}Su (k + 1)
u′ =1
u′ 6=u
U
X (eq)∗ (eq) (eq) (eq)∗
+ {Ĥu,2 (k)Hu′ ,1 (k) − Ĥu,1 (k + 1)Hu′ ,2 (k + 1)}Su (k)
u′ =1
u′ 6=u
(eq)∗
+ Ĥu,2 (k){ICI(ǫ, k) + b(k)}
(eq)
+ Ĥu,1 (k + 1){ICI ∗ (ǫ, k + 1) + b∗ (k + 1)}
(3.38)
Then, assuming that Su are independent and identically distributed, and applying
the central limit theorem, a Gaussian approximation for ςu (k) and Ŝu (k) can
be considered. Thus, they can be completely characterized by their means and
variances. The modified ESE process for the uth user at the ith iteration generates
the value of the a posteriori LLR as follows:
For the initialization process guESE (k) = 0 for i = 1, where i denotes the iteration
number with i ∈ {1, . . . , Iidma }, with Iidma the maximum number of iterations.
The outputs of the modified ESE are inserted in the user decoding process pre-
sented in chapter 1 in section 1.4.3.2. Then, at the ith iteration for i < Iidma ,
the outputs of the decoding process are the values of guESE (k), that are reinserted
in the ESE. The ESE process is again performed. Finally, the decoding process
produces hard decisions to obtain the estimations of the transmitted bits during
121
3. Frequency Synchronization and Channel Equalization of
OFDM-IDMA Uplink Systems
the final iteration Iidma . Figure 3.7 shows the proposed STBC-OFDM-IDMA
scheme. 1st user
gESE
1 (k) estimated
bits
g1ESE(k)
Ŝ1 (k)
r(k) user U th user
modified
MRC estimated
ESE gESE
U (k) decoding
ŜU(k) bits
ESE(k)
OFDM gU
demodulation
(eq) σz2 = σB
2 + σ2
ICI
Ĥu
ĥ, ǫ
R SPKF initialization step
Figure 3.7: Proposed STBC-OFDM-IDMA uplink receiver, where the user de-
coding is the same as in figure 1.20
122
3.3 Space-Time Block Code OFDM-IDMA Receiver
estimated. The MMSE for an AWGN channel and for channels with different
length are presented. The channel estimation for different channel lengths give
the same estimation performance, around 10−2 for a SNR= 0dB. The following
−2
10
SPKF Lu = 3
SPKF Lu = 1
SPKF AWGN
MMSE
−3
10
−4
10
0 2 4 6 8 10 12 14 16
SNR (dB)
results are the average of 5000 OFDM-IDMA symbols. Figure 3.9 and 3.10 show
the BER for a transmission over a normalized Rayleigh slow-fading frequency
selective channel with Lu,p = 1 and Lu,p = 3 multipaths ∀u, p respectively. The
efficiency of the OFDM-IDMA-STBC proposed scheme is confirmed, we can see
that for a BER of 10−2 there are less than 1 dB between the curve without CFO
and the curve of the proposed estimator. In addition, we can see a gain of around
3dB of the proposed scheme in comparison with the linear detection method
presented in [Cao 04a].
123
3. Frequency Synchronization and Channel Equalization of
OFDM-IDMA Uplink Systems
−1
10
BER
−2
10
No CFO - Channel perfectly estimated
SPKF-CFO estimated - Channel perfectly estimated
SPKF-CFO/channel estimated
Linear detection - CFO and Channel perfectly estimated
0 2 4 6 8 10 12 14
SNR (dB)
−1
10
−2
10
BER
−3
10
3.4 Conclusions
In this section, two OFDM-IDMA receivers are proposed. A joint CFO/channel
estimation is performed using a SPKF. Then, the estimated parameters are in-
124
3.4 Conclusions
125
Conclusions and Perspectives
Estimating the CFOs and channels plays a key role in the design of mobile wire-
less systems. In this PhD, our purpose was to create new receiver schemes that
take into account both phenomena. We have proposed to jointly estimate the
CFOs and the channels for an OFDMA uplink system by using a recursive ap-
proach [Pove 10]. On the one hand, the EKF, the SOEKF, the IEKF and the
SPKF including the UKF and the CDKF were considered. On the other hand,
to relax the Gaussian assumptions on the model noise and measurement noise,
the “extended H∞ filter” and the “unscented H∞ filter” were studied.
From the estimation point of view, we first compared the different Kalman tech-
niques. We can draw the following conclusions:
1/ The SOEKF gives slightly “better” estimates of the CFO and the channels
than the EKF in terms of estimation accuracy. Nevertheless, the computational
cost of the SOEKF is higher than the EKF one, due to the calculation of the
Hessian matrix for the second-order expansion.
2/ Concerning the SPKF, the UKF and the CDKF give the same results in terms
of estimation accuracy. In addition, their computational costs are quite similar.
They converge faster than the EKF and the SOEKF. They have the advantage
of avoiding the calculation of Jacobians and Hessians matrices, but the Cholesky
factor of the autocorrelation matrices must be computed.
3/ The IEKF gives the “best” estimation performance in terms of estimation
accuracy and convergence speed. It needs less samples to converge to the true
value than the other Kalman approaches. However, due to the calculation of the
Jacobian matrix at each iteration, the computational complexity is the highest
among the proposed approaches.
Let us now look at the relevance of the H∞ approaches.
1/ the “unscented H∞ filter” gives better estimation performance than the “ex-
tended H∞ filter”. In addition, its convergence is faster than the one of the
“extended H∞ filter”.
2/ when the noise characteristics are available, the “extended H∞ filter” and the
127
Conclusions and Perspectives
EKF give similar results. In addition, there is no real difference between the “un-
scented H∞ filter” and the SPKF. If the noise variances are a priori unknown,
the H∞ approaches converge faster than the Kalman algorithms. Nevertheless,
the noise attenuation level and the weighting matrices have to be chosen by the
practitioner. If the noise attenuation level is set to a value that is too small, a
solution to the H∞ issue may not exist. If the noise attenuation level is set to a
high value to guarantee the existence of a solution, the resulting H∞ filter “looks
like” the Kalman filter.
3/ The computational complexity of the “extended H∞ filter” is slightly higher
than the EKF while the computational cost of the “unscented H∞ filter” is slightly
higher than the SPKF.
Then, the above techniques have been compared with the EM-based algorithms
proposed in [Pun 04a] and in [Pun 06] in terms of estimation accuracy and com-
putational complexity. Even if the EM-based estimator gives better estimation
performance, its computational complexity is much higher than the costs of the
SPKF and the “unscented H∞ filter”. Our methods are therefore very useful for
practical cases.
Based on the above considerations, we think that the SPKF gives the best com-
promise in terms of estimation performance, computational complexity and a
priori information to be tuned.
From the digital communication point of view, the above optimal filtering tech-
niques have been combined with an MMSE-SD to obtain a non-pilot aided CFO
estimator for an OFDMA uplink [Pove 09a], [Pove 09b], [Pove 11c]. The MAI
suppression scheme followed by the CFO estimation algorithm based on optimal
filtering can jointly estimate and detect each user’s CFO and symbols respectively.
According to the simulations we carried out, the proposed scheme can effectively
suppress the MAI caused by a relatively large CFO, and it is sufficiently robust
to CFO variations. The iterative decoding applies simple hard interference can-
cellation techniques, resulting in a moderate complexity.
Then, the SPKF has been chosen to be combined with two statistical tests: the
128
Conclusions and Perspectives
BHT and the CUSUM test [Pove 11b], [Pove 11d] to detect the beginning and
the end of a CR-NBI disturbing the estimation of the CFOs and the channels of
different users in a PU-OFDMA system. The proposed estimator operates in the
time domain unlike other approaches such as [More 08] that estimate the CFO
in presence of an NBI in the frequency domain. By using a SPKF and without
increasing the computational complexity, we take advantage of the innovation in
the Kalman filter to detect the CR-NBI. When using the BHT, better perfor-
mance are obtained in terms of the CR-NBI probability of detection. However,
when using the CUSUM, better performance in terms of PU probability of detec-
tion are obtained. The results show that the SPKF without the statistical tests is
not able to estimate the CFO of each user. The proposed estimator outperforms
the existing CFO estimation methods such as [More 08] in presence of NBI and
it is robust to the NBI power. Nevertheless, the architecture is sensitive to the
duration in time of the CR-NBI.
We have also proposed to combine the SPKF with a modified version of an itera-
tive IDMA receiver to design a new OFDM-IDMA receiver [Pove 12]. Simulation
results show that the SPKF provides “accurate” estimations of the CFO and the
channel of each user. In addition, the modified OFDM-IDMA receiver accurately
estimates the transmitted bits of each user, by taking into account the estima-
tions of the SPKF and the interference produced by the CFO. According to the
simulation tests we carried out, our receiver outperforms CFO correction meth-
ods proposed for other multicarrier systems such as [Cao 04a].
Our last contribution consists in combining the SPKF with a STBC and a new
variant of an iterative IDMA receiver to take advantage of the space diversity
[Pove 11a]. The simulation results confirm that the SPKF provides accurate es-
timations of the CFO and the channel. In addition, the STBC-OFDM-IDMA
receiver can accurately estimate the transmitted bits of each user. We show
that the proposed receiver outperforms existing CFO correction methods such as
[Cao 04a].
129
Conclusions and Perspectives
1
Due to the large number of independent modulated subcarriers that are added up coher-
ently, the transmitted signal exhibits a high PAPR.
130
Appendix A
Kalman Filter - Linear Case
Firstly presented in [Kalm 60] by Rudolf E. Kalman, the Kalman filter (KF) is
based on a state-space representation of the system. It has been used in a wide
range of applications, from speech enhancement to time-varying autoregressive
parameters tracking, from biomedical applications to mobile communications.
It uses two equations to describe the system. When the state-space equations
are linear and the noises are additive white zero-mean Gaussian processes, they
satisfy:
State equation:
x(n) = Φ(n − 1)x(n − 1) + Γw(n) (A.1)
Measurement equation:
y(n) = Ψ(n)x(n) + b(n) (A.2)
where x(n) is the state vector1 of size U at time n and y(n) is the measurement
vector of size K at time n. The model noise w and the observation noise b are
uncorrelated white zero-mean Gaussian vectors with covariance matrices assumed
to be Q = σw2 IU and R = σb2 IK respectively. In addition, Φ(n−1) is the transition
matrix of size U × U from time n − 1 to n, Γ is the input gain matrix of size
U × U and Ψ(n) is the measurement matrix K × U at time n.
The KF operates in two steps: the prediction step and the filtering step.
In the prediction step, the KF uses the estimated state at the previous instant to
estimate the current state, without taking into account the current observation.
This is the so-called a priori estimation of the state vector defined as follows:
131
A. Kalman Filter - Linear Case Appendix A
Given (A.1) and (A.3) and as w is a zero-mean AWGN, the a priori estimation
of the state vector can be expressed as follows:
where x̂(n − 1|n − 1) is the estimation of the state vector at time n − 1 given the
set of observations [y(0), y(1), . . . , y(n − 1)].
Then, let us introduce the a priori estimation error:
Given (A.1) and (A.4), the a priori estimation error of the state vector can be
rewritten as:
x̃(n|n − 1) = Φ(n − 1)x(n − 1) + Γw(n) − Φ(n − 1)x̂(n − 1|n − 1)
(A.7)
= Φ(n − 1)x̃(n − 1|n − 1) + Γw(n)
Then, the covariance matrix of x̃(n|n − 1) satisfies:
h i
P(n|n − 1) = Φ(n − 1)E x̃(n − 1|n − 1)x̃H (n − 1|n − 1) ΦH (n − 1)
h i
+ ΓE w̃(n)w̃H (n) ΓH (A.8)
= Φ(n − 1)P(n − 1|n − 1)ΦH (n − 1) + ΓQΓH
132
A. Kalman Filter - Linear Case
At that stage, the filtering step uses the current observation to correct the a
priori estimated state vector, in order to obtain the a posteriori estimation. To
estimate the current state and by following the kind of update equation obtained
in the RLS algorithm, one has:
where K(n) is the Kalman gain and ỹ(n) is the innovation defined as:
and it satisfies:
K(n) = {Pxy (n)}{Pyy (n)}−1 (A.12)
and Pxy (n) is the cross-covariance matrix between y(n) and the state prediction
x̂(n|n − 1) defined as:
Finally, by combining (A.9) and (A.12) it can be easily shown that the error
estimation covariance matrix is updated as follows:
133
A. Kalman Filter - Linear Case Appendix A
where P{x(n)|Y(n − 1)} = P{x(n)|y(0), y(1), . . . , y(n − 1)} is the a priori prob-
ability density of the state vector. As it is Gaussian, it can be completely char-
acterized by its mean and its covariance matrix. Then, by using the Bayes rule
the a posteriori probability density can be obtained as follows:
1
The Riccati equation is defined as: P = A{I − PBH (BPBH + R)−1 B}PAH + Q.
134
Appendix B
H∞ Filter - Linear Case
The H∞ approach was introduced in the field of control in 1981 [Zame 81]. A
solution of the H∞ estimation problem can be based on polynomial decompo-
sition techniques [Grim 90]. Nevertheless, they lead to equations with high-
computational cost that cannot be used in practical situations. As an alternative
solution, state-space approaches have emerged. In [Gero 99], the author presents
a solution of the H∞ estimation constraint, based on the resolution of a convex
optimization problem under linear matrix inequality constraints approaches. In
[Shak 92], Shaked et al. present a state-space approach solution based on the
resolution of a quadratic Riccati-type equation. This algorithm is easy to imple-
ment and has a lower computational cost than the above approaches.
The H∞ filter has been widely used in the field of control. For the last 12 years,
several studies have been conducted by the signal processing community.
In [Shen 99], instead of using a Kalman filtering, Shen et al. suggest using an
H∞ filter to enhance a speech signal disturbed by an additive noise and recorded
from a single microphone. For this purpose, the signal is assumed to be modeled
by an AR process. However, the AR parameters are unknown and hence need
to be estimated. Shen et al. propose to estimate them directly from the noisy
observations by using a second H∞ filter. Therefore, the resulting AR parameter
estimates are biased, as pointed out by Labarre et al. in [Laba 05], who inves-
tigated the relevance of the H∞ filtering for speech enhancement. To avoid the
problem of biased AR parameter estimates that can be obtained in [Shen 99],
Labarre et al. suggest estimating both the AR model and its parameters. Al-
though this leads to a non-linear estimation issue, they have developed a structure
based on two mutually interactive standard H∞ filters [Laba 07]. The first one
aims at estimating the AR model, while the second one updates the estimation
of the AR parameters. In addition, in [Jamo 07a], the authors take advantage
135
B. H∞ Filter - Linear Case Appendix B
of the two mutually interactive H∞ filter based approach to jointly estimate the
fading channel and its AR parameters. However, the authors do not obtain better
performance than a KF based method.
Let us now look more carefully at the H∞ filtering approach. Given the state-
space model in (A.1) and (A.2), let us introduce a third state-space equation to
focus on a linear combination of the state-vector components:
where N ob denotes the number of available observations, e(n) = z(n) − ẑ(n), and
V and W are positive weighting matrices tuned by the practitioner to achieve
performance requirements.
1
V− 2 (n)b(n)
1
e(n) = z(n) − ẑ(n)
W− 2 (n)w(n) Transfer operator
P(0)e(0)
Figure B.1: Transfer operator from disturbances to estimation error for H∞ -norm
based estimation
1
Namely w(n), b(n) and the initial conditions on the state vector.
136
B. H∞ Filter - Linear Case
J ∞ < Ξ2 (B.3)
where
" # " #
V 0 Ψ(n) h i
M= + P∞ (n|n − 1) ΨH (n) LH
0 −Ξ2 IU L
" # (B.5)
Ψ(n)P∞ (n|n − 1)ΨH (n) + V Ψ(n)P∞ (n|n − 1)LH
=
LP∞ (n|n − 1)ΨH (n) LP∞ (n|n − 1)LH − Ξ2 IU
When using the H∞ filter, the state vector can be estimated recursively as follows:
137
B. H∞ Filter - Linear Case Appendix B
It should be noted that the purpose of the H∞ filter is to minimize the peak error
power in the frequency domain whereas the Kalman filter minimizes the average
error power [Grim 90].
138
Appendix C
Kalman Filter vs H∞ Filter
In this appendix, our purpose is to compare the H∞ filter and the KF in the linear
case when estimating model parameters, by comparing the Ricatti equations of
both algorithms. To our knowledge, there is no work dealing with this kind of
comparison.
Let us first assume that:
and
L = IU (C.2)
and
Pyy∞ (n) = Pyy (n) (C.4)
Then, let us compare P∞ (n + 1|n) and P(n + 1|n). For this purpose, the matrix
M defined in (B.6) must be inverted in (B.4). Among the approaches that could
be considered such as the matrix inversion lemma, we suggest using the one based
on the Schur complement [Brez 88] for the sake of simplicity1 . Indeed by defining
" # " #" #" #
−1 −1
1 A B I K −A B A 0 I K 0
M= , then M−1 =
C D 0 IU 0 S−1 −CA−1 IU
where A is matrix of size K × K, B is matrix of size K × U , C is matrix of size U × K, D is
matrix of size U × U , S = D − CA−1 B is the Schur complement of A in M and 0 is a zero
matrix.
139
C. Kalman Filter vs H∞ Filter Appendix C
Υ(n) = P(n|n − 1)
(C.1),(C.4)
= P(n|n) − Ξ2 IU (C.6)
(A.15)
Then given the previous footnote, M−1 can be expressed as the product of three
matrices, the coefficients of which are defined from the coefficients of M and the
Schur complement Υ(n) as follows:
" #" #" #
−1 IK −KH (n) {Pyy (n)}−1 0 IK 0
M =
0 IU 0 Υ−1 −K(n) IU
= W(n)X(n)Y(n) (C.7)
At that stage, given (C.7), we can rewrite the Riccati recursion (B.4) as:
140
C. Kalman Filter vs H∞ Filter
Hence, the solution of the Riccati equation when using the H∞ filter satisfies:
where
141
C. Kalman Filter vs H∞ Filter Appendix C
2 Ξ
where {− λuλ−Ξ
u
2 }u∈{1,2,...,U } are the eigenvalues of Q (n).
It should be noted that when Ξ tends to +∞, QΞ (n) tends to be a zero matrix and
P∞ (n + 1|n) tends to P(n + 1|n).
If Ξ2 > λu , for u ∈ {1, 2, . . . , U }, QΞ (n) is a positive-definite matrix. Therefore,
the solution to the Ricatti equation for the H∞ filter can be seen as an upper
bound of the Kalman a priori error covariance matrix.
142
Appendix D
Extended Kalman Filter Approaches
Measurement equation:
where ∇Φn |x̂(n−1|n−1) is the Jacobian matrix of Φn evaluated for x̂(n − 1|n − 1).
Given (D.3), (D.1) can be rewritten as follows:
143
D. Extended Kalman Filter Approaches Appendix D
Given (A.3) and (D.4), the a priori estimation of the state vector is defined as
follows:
x̂(n|n − 1) = Φn (x̂(n − 1|n − 1))
+ ∇Φn |x̂(n−1|n−1) E [x(n − 1)|y(0), y(1), . . . , y(n − 1)]
(D.5)
− ∇Φn |x̂(n−1|n−1) x̂(n − 1|n − 1)}
+ ΓE [w(n)|y(0), y(1), . . . , y(n − 1)]
Thus, given (D.4) and (D.7), the a priori error estimation can be defined as:
x̃(n|n − 1) = x(n) − x̂(n|n − 1)
= Φn (x̂(n − 1|n − 1)) + ∇Φ|x̂(n−1|n−1) x(n − 1)
− ∇Φn |x̂(n−1|n−1) x̂(n − 1|n − 1) + Γw(n) − Φn (x̂(n − 1|n − 1))
= ∇Φn |x̂(n−1|n−1) x̃(n − 1|n − 1) + Γw(n)
(D.8)
and its corresponding covariance matrix P(n|n − 1) is defined by:
It should be noted the a priori estimation of the state vector (D.7) does not
depend on the Jacobian matrix ∇Φn |x̂(n−1|n−1) . However, ∇Φn |x̂(n−1|n−1) is re-
quired to calculate the a priori error estimation covariance matrix P(n|n − 1).
Now, let us focus our attention on the a posteriori estimation of the state vector.
For this purpose, it is necessary to calculate the first-order Taylor expansion of
(D.2) around x̂(n|n − 1) as follows:
144
D. Extended Kalman Filter Approaches
where ∇Ψn |x̂(n|n−1) is the Jacobian matrix of Ψn evaluated for x̂(n|n − 1).
Then, given (D.2) and (D.10), the innovation can be expressed as:
When looking at the linear case and more particularly (A.10), the innovation is
defined as follows:
A similarity can be noticed between (D.11) and (D.12). The EKF can be hence
easily derived. The way the state-vector estimation can be updated and the
definition of the Kalman gain can be obtained similarly as in appendix A, by
replacing Ψ(n) by ∇Ψn |x̂(n|n−1) .
Thus, to obtain the a posteriori estimation of the state vector, one has:
y(n) K(n)
145
D. Extended Kalman Filter Approaches Appendix D
and Pxy (n) the cross-covariance matrix between y(n) and the state prediction
x̂(n|n − 1) defined as:
The EKF may sometimes diverge if the first-order approximation is not sufficient
to describe the non-linearity. Therefore, a second-order approximation can be
used and leads to the SOEKF [Bar 01]. It is still based on a Taylor expansion,
but a second-order expansion is considered.
Given the state-space equations in (D.1) and (D.2), let us look at the second-order
Taylor expansion of (D.1) around x̂(n − 1|n − 1):
146
D. Extended Kalman Filter Approaches
∂ 2 Φn,u
Φn,u = tr{ |x̂(n−1|n−1)
∂x2 (n)
× {x(n − 1) − x̂(n − 1|n − 1)}{x(n − 1) − x̂(n − 1|n − 1)}H }
(D.20)
= tr{{x(n − 1) − x̂(n − 1|n − 1)}{x(n − 1) − x̂(n − 1|n − 1)}H
∂ 2 Φn,u
× 2 |x̂(n−1|n−1) }
∂x (n)
Remark 2: one can approximate the expectation of Φn,u given the set of obser-
vations {y(0), . . . , y(n − 1)} by:
mean
Φn,u = E{Φn,u |y(0), . . . , y(n − 1)}
∂ 2 Φn,u (D.21)
≈ tr{P(n − 1|n − 1) |x̂(n−1|n−1) }
∂x2 (n)
mean
The above approximation of Φn,u has the advantage of being “easily” computed
∂ 2 Φn,u
by using P(n − 1|n − 1) and |
∂x2 (n) x̂(n−1|n−1)
.
Given (D.19) and (D.21), one has:
mean
Φn = E{Φn |y(0), . . . , y(n − 1)}
U
X (D.22)
mean
= ΘΦ
u Φn,u
u=1
By following the same development as in (D.5) and (D.6), the a priori estimation
of the state vector can be obtained as follows:
1 mean
x̂(n|n − 1) = Φn (x̂(n − 1|n − 1)) + Φn (D.23)
2
mean
By assuming that Φn ≈ Φn and by combining (D.18) and (D.23), the a priori
estimation error x̃(n|n − 1) can be expressed as follows:
147
D. Extended Kalman Filter Approaches Appendix D
∂ 2 Ψn,k
Ψn,k = tr{ |x̂(n|n−1)
∂x2 (n)
× {x(n) − x̂(n|n − 1)}{x(n) − x̂(n|n − 1)}H }
(D.28)
= tr{{x(n) − x̂(n|n − 1)}{x(n) − x̂(n|n − 1)}H
∂ 2 Ψn,k
|x̂(n|n−1) }
∂x2 (n)
Then, given (D.2) and (D.26) the innovation can be expressed as follows:
1
y(n) − Ψn (x̂(n|n − 1)) ≈ ∇Ψn |x̂(n|n−1) {x(n) − x̂(n|n − 1)} + Ψn + b(n)
2
(D.29)
Unlike the Kalman filter in the linear case or the EKF, the innovation depends
on Ψn . Therefore, to obtain the a posteriori estimation of the state vector, the
following relation can be considered:
1
x̂(n|n) = x̂(n|n − 1) + K(n){y(n) − Ψn (x̂(n|n − 1)) − Ψn } (D.30)
2
where K(n) is the Kalman gain defined in (D.14) and the error estimation co-
variance matrix is updated as in (D.17).
148
D. Extended Kalman Filter Approaches
However, in practical case x(n) − x̂(n|n − 1) is not available and hence Ψn cannot
be computed. For this reason, instead of Ψn one suggests using:
K
X mean
Ψn = ΘΨ
k Ψn,k (D.31)
k=1
mean ∂2Ψ
n,k
where Ψn,k = tr{P(n|n − 1) ∂x2 (n) |x̂(n−1|n−1) }.
As consequence, the state vector is updated as follows:
1
x̂(n|n) = x̂(n|n − 1) + K(n){y(n) − Ψn (x̂(n|n − 1)) − Ψn } (D.32)
2
The IEKF is another alternative that aims at improving the EKF [Gelb 74]. The
purpose is to linearize the measurement model around the updated state x̂(n|n),
instead of the predicted state x̂(n|n − 1). Indeed, the Bayesian approach makes
it possible to separate the Kalman filter in two steps: the prediction step and the
filtering step. The nearer the linearization point is to the real value, the lower
the linearization error should be. In addition, after the filtering step, x̂(n|n) is
expected to have a lower error variance than x̂(n|n − 1). For the above reasons,
using the IEKF can be of interest.
Given the state-space equations in (D.1) and (D.2), the a priori estimation of
the state vector and its covariance matrix can be obtained as in (D.7) and (D.9)
respectively.
To obtain the a posteriori estimation of the state vector, an iterative process
is introduced. At the i + 1th iteration, with i ∈ {0, 1, 2, . . . , Ikf }, where Ikf is
T : |x̂i (n|n) − x̂i+1 (n|n)| < αcvg Iu×1
x̂(n|n − 1) x̂(n|n)
x̂(n − 1|n − 1) = x̂0 (n|n) x̂i+1 (n|n) = x̂i+1 (n|n)
(D.7) (D.33) i ≥1 T
OK OK
y(n)
Not OK Not OK
Ki (n) i = i+1 i = i+1
149
D. Extended Kalman Filter Approaches Appendix D
the number of iterations, the (i + 1)th a posteriori estimate of the state vector
satisfies:
where
Ki (n) = {Pxy yy
i (n)}{Pi (n)}
−1
(D.34)
Pyy H
i (n) = ∇Ψn |x̂i (n|n) P(n|n − 1)∇Ψn |x̂i (n|n) + R (D.35)
Pxy H
i (n) = P(n|n − 1)∇Ψn |x̂i (n|n) (D.36)
and x̂0 (n|n) = x̂(n|n − 1). Figure D.2 shows the block scheme of the IEKF.
The estimation-error covariance matrix has to be updated as follows:
Then, having a new value of the estimate, the procedure is repeated until the
difference between two consecutive estimates is lower than a predefined value
αcvg :
|x̂i (n|n) − x̂i+1 (n|n)| < αcvg IU ×1 αcvg > 0 (D.38)
150
Appendix E
Sigma-Point Kalman Filter
When dealing with the SPKF, the state distribution is still approximated by a
Gaussian distribution, which is now characterized by a set of points lying along
the main eigenaxes of the random variable covariance matrix, as shown in figure
E.1.
−3
x 10
5
4
Probability density
0
10
5 10
0 5
0
−5
−5
x2 −10 −10 x1
These so-called sigma-points are propagated through the non-linear system (D.2),
as shown in figure E.2. A weighted combination of the resulting values makes
it possible to estimate the mean and the covariance matrix of the transformed
random vector, i.e. the Gaussian random variable that undergoes the non-linear
transformation.
On the one hand, the UKF is based on the unscented transformation [Wan 01].
151
E. Sigma-Point Kalman Filter Appendix E
sigma-points
Non-linear
transformation
y(n) = Ψn (x(n))
When the density is odd, the weights are chosen to provide the exact second-order
Taylor expansion around the mean of the random variable. The UKF was pro-
posed by their authors as an alternative to the EKF to avoid the linearization step.
Nevertheless, Lefebvre et al. [Lefe 02] showed that the sigma-point approach cor-
responds to a weighted statistical linear regression (WSLR). In [VdMe 01], a
square-root UKF is presented. It has better numerical properties and guarantees
the positive semi-definiteness of the underlying state covariance.
On the other hand, the CDKF is based on the second-order Sterling polynomial
interpolation1 formula [Ito 00]. The linearization is based on Sterling polynomial
interpolation around the last available estimate of the state vector. The CDKF
can be seen as a Bayesian approach of a second-order divided difference filter
(DDF). According to Nörgaard et al. [Norg 00], if the last available estimate is
“far” from the real value, the Sterling polynomial interpolation approximation of
the non-linear function is better than a Taylor expansion approximation.
Let us consider the state-space representation done in (D.1) and (D.2). Thus,
as the size of the state vector is U , there are 2U + 1 sigma-points stored in the
following matrix:
h i
X(n − 1|n − 1) = X0 (n − 1|n − 1) Xa (n − 1|n − 1) Xb (n − 1|n − 1) (E.1)
1
Sterling polynomial interpolation [Boas 64] is the approximation in an infinite sum of terms
of a function yi = f (xi )|i∈Z = f (x0 + hu) around a value x0 , where u = xi −xh
0
and h is the
distance between two points chosen arbitrarily. Thus, the second-order Sterling formula of the
2
function f (xi ) is yi = f (xi ) ≈ y0 + u2 (y1 − y−1 ) + u2 (y1 + y−1 − 2y0 ).
152
E. Sigma-Point Kalman Filter
where
q
⋆
Xb (n − 1|n − 1) = x̂(n − 1|n − 1) − Λ P(n − 1|n − 1) m ∈ {U, . . . , 2U }
m
(E.2)
q q
( P(n − 1|n − 1))m is the mth column of the matrix square-root, (.) is the so-
⋆ ⋆
called matrix square-root that can be obtained by using the lower triangular
Cholesky decomposition and Λ is a scalar scaling factor that determines the
spread of the sigma-points. These points are deterministically chosen so that
they completely capture the true mean and the covariance matrix of the Gaussian
random variable.
The difference between CDKF and UKF stands in the way the mean and the
covariance matrix of the transformed Gaussian random variable are calculated.
The CDKF uses only a single scalar scaling parameter Λ, as opposed to three
√
required by the UKF. When using the CDKF, Λ = 3 [VdMe 04b]. If the UKF
is used, one has:
• Λ = ι2 (U + κ) + U , with 10−4 ≤ ι ≤ 1,
Let us define the weight scalar parameters for the UKF [Wan 01] as follows:
η
w0(c0) =
η+U
η
w0(c1) = + (1 − ι2 + η) (E.3)
η+U
1
wm(c0) = wm(c1) =
2(U + Λ)
153
E. Sigma-Point Kalman Filter Appendix E
while the a priori estimation error covariance matrix when using the UKF satis-
fies:
2U
X h i⋆
P(n|n − 1) = wm(c1) Xm (n|n − 1) − x̂− (n|n − 1) +Q (E.7)
m=0
where:
X0 (n|n − 1) = x̂(n|n − 1)
q
⋆
Xa (n|n − 1) = x̂(n|n − 1) + Λ P(n|n − 1) m ∈ {1, . . . , U } (E.10)
q m
⋆
Xb (n − 1|n − 1) = x̂(n|n − 1) − Λ P(n|n − 1) m ∈ {U, . . . , 2U }
m
154
E. Sigma-Point Kalman Filter
Once the a priori sigma points are calculated, they are propagated through the
non-linear function:
Y(n|n − 1) = Ψn (X(n|n − 1)) (E.11)
The mean and the covariance of y(n) are approximated by using the posterior
sigma points as follows:
2U
X
ŷ− (n) = wmc0 Ym (n|n − 1) (E.12)
m=0
On the one hand, when using the UKF, the covariance matrix Pyy (n) of the
innovation ỹ(n) = y(n) − ŷ− (n) and the cross-covariance matrix Pxy (n) between
y(n) and the state prediction error satisfy:
2U
X h i⋆
yy
P (n) = wmc1 Ym (n|n − 1) − ŷ− (n) +R (E.15)
m=0
2U
X h i h i
xy
P (n) = wmc1 Xm (n|n − 1) − x̂− (n) × Ym (n|n − 1) − ŷ− (n) (E.16)
m=0
On the other hand, when using the CDKF, the covariance matrix of the innova-
tion and the cross-covariance matrix between y(n) and the state prediction error
satisfy:
U
X
yy
P (n) = {wmc1 [Ym (n|n − 1) − Ym+U (n|n − 1)]⋆
m=0 (E.17)
+ wmc2 [Ym (n|n − 1) + Ym+U (n|n − 1) − 2Y0 (n|n − 1)] } + R ⋆
q
Pxy (n) = ⋆ wmc1 P(n|n − 1) [Y1:U (n|n − 1) − YU +1:2U (n|n − 1)]H (E.18)
Finally, the estimation error covariance matrix is updated as follows:
155
Appendix F
Extended and Unscented H∞ Filter
P∞ (n + 1|n) = ∇ΦH ∞ H
n+1 |x̂(n|n) P (n|n)∇Φn+1 |x̂(n|n) + ΓWΓ
H
(F.1)
h i
H
= ∇Φn+1 |x̂(n|n) P∞ (n|n − 1){IU − ∇ΨH
n |x̂(n|n−1) L
" #
−1 ∇Ψn |x̂(n|n−1)
×M P∞ (n|n − 1)}∇ΦH
n+1 |x̂(n|n)
L
+ΓWΓH (F.2)
where ∇ΦH
n+1 |x̂(n|n) is the Jacobian matrix of Φn+1 evaluated for x̂(n|n), ∇Ψn |x̂(n|n−1)
is the Jacobian matrix of Ψn evaluated for x̂(n|n − 1) and
" # " #
V 0 ∇Ψn |x̂(n|n−1) h
H
i
M= + P∞ (n|n − 1) ∇ΨH
n |x̂(n|n−1) L
0 −Ξ2 IU L
" H
#
∇Ψn |x̂(n|n−1) P∞ (n|n − 1)∇ΨH ∞
n |x̂(n|n−1) + V ∇Ψn |x̂(n|n−1) P (n|n − 1)L
=
LP∞ (n|n − 1)∇ΨH n |x̂(n|n−1) LP∞ (n|n − 1)LH − Ξ2 IU
(F.3)
157
F. Extended and Unscented H∞ Filter Appendix F
Remark 1: when Ξ2 tends to +∞, the “extended H∞ filter” reduces to the EKF
[Hass 99].
Remark 2: the computational cost of the “extended H∞ filter” is slightly higher
than the EKF.
158
F. Extended and Unscented H∞ Filter
where ŷ− (n) is obtained as in (E.12) by performing the same process as in the
UKF, but by replacing the matrix Q by W. The gain K∞ (n) is obtained as in
(F.5), but Pxy∞ (n) and Pyy∞ (n) are now approximated as follows:
2U
X h i⋆
Pyy∞ (n) ≈ wmc1 Ym (n|n − 1) − ŷ− (n) +V (F.9)
m=0
and
2U
X h i h i
Pxy∞ (n) ≈ wmc1 Xm (n|n − 1) − x̂− (n) × Ym (n|n − 1) − ŷ− (n) (F.10)
m=0
where the wmc1 , Xm (n|n − 1), x̂− (n) and Ym (n|n − 1) are obtained as in (E.3),
(E.5), (E.6) and (E.11) respectively, by the using the UKF process.
Considering L = IU , given (F.6), (F.7), (F.9) and (F.10) the matrix M in (F.3)
can be rewritten as:
" #
Pyy∞ (n) {Pxy∞ (n)}H
M= (F.11)
Pxy∞ (n) P∞ (n|n − 1) − Ξ2 IU
159
References
[3GPP 09] 3GPP Standard. Technical Specification Group Radio Access Net-
work; Evolved Universal Terrestrial Radio Access (E-UTRA); Physi-
cal layer procedures (Release 8). 3GPP TS 36.213 V8.8.0, September
2009. 1, 25, 26
[Arsl 07] H. Arslan. Cognitive Radio, Software Defined Radio, and Adaptive
Wireless Systems. Springer, 2007. 10
[Bian 05] H. Bian, Z. Jin, and W. Tian. Study on GPS Attitude Determination
System Aided INS Using Adaptive Kalman Filter. Measurement
Science and Technology, Vol. 6, No. 10, pp. 2072–2079, September
2005. 80
161
References References
[Chen 10] C. J. Chen, T. T. Lin, and T. C. Chen. Blind Joint Frequency Off-
set and Channel Estimation for MC-CDMA Systems over Multipath
Fading Channels. Proceedings of the Mechanical and Electronics En-
gineering (ICMEE ’10), pp. 137–141, August 2010. 104
162
References References
[Cibl 04] P. Ciblat and E. Serpedin. A Fine Blind Frequency Offset Estimator
for OFDM/OQAM Systems. IEEE Transactions on Signal Process-
ing, Vol. 52, No. 1, pp. 291–296, January 2004. 52
[Dini 10] D. H. Dini and D. P. Mandic. Analysis of the Widely Linear Complex
Kalman Filter. Proceedings of the UDRC Sensor Signal Processing
for Defense (SSPD ’10), pp. 1–4, September 2010. 55
[ETSI 01] ETSI Standard. Interaction Channel for Digital Terrestrial Televi-
sion (RCT) Incorporating Multiple Access OFDM. ETSI DVB RCT,
March 2001. 24
163
References References
[ETSI 95] ETSI Standard. Radio Broadcasting Systems: Digital Audio Broad-
casting to Mobile, Portable and Fixed Receivers. ETS 300 401,
European Telecommunications Standard, 1995. 24
[ETSI 97] ETSI Standard. Digital Video Broadcasting (DVB-T); Frame Struc-
ture, Channel Coding, Modulation for Digital Terrestrial Television.
ETS 300 744, European Telecommunications Standard, 1997. 24
[Fu 06] X. Y. Fu, H. Minn, and C. Cantrell. Two Novel Iterative Joint
Frequency-Offset and Channel Estimation Methods for OFDMA Up-
link. Proceedings of the IEEE Global Communications Conference
(GLOBECOM ’06), pp. 1–6, November 2006. 2, 52
[Gelb 74] A. Gelb. Applied Optimal Estimation. MIT Press, 1974. 149
164
References References
[Hayk 96] S. Haykin. Adaptative Filter Theory. Prentice Hall, 1996. 56, 143
165
References References
[Huan 11] Y. Huang, H. Li, K. A. Campbell, and Z. Han. Defending False Data
Injection Attack on Smart Grid Network Using Adaptive CUSUM
Test. Proceedings of the IEEE Conference on Information Sciences
and Systems (CISS ’11), pp. 1–6, March 2011. 89
[IEEE 06] IEEE Standard. Part 16: Air Interface for Fixed and Mobile Broad-
band Wireless Access System; Amendment 2: Physical and Medium
Access Control Layers for Combined Fixed and Mobile Operation in
Licensed Bands. IEEE Std. 802.16e-2005, February 2006. 1, 3, 24
[IEEE 99] IEEE Standard. Part 11: Wireless LAN Medium Access Con-
trol (MAC) and Physical Layer (PHY) Specifications; Higher-Speed
Physical Layer Extension in the 5 GHz Band. IEEE Std. 802.11a,
December 1999. 3, 24
[Ito 00] K. Ito and K. Xiong. Gaussian Filters for Nonlinear Filtering Prob-
lems. IEEE Transactions on Automatic Control, Vol. 5, No. 5,
pp. 910–927, May 2000. 56, 130, 152
[ITU 08] ITU Report. ITU Corporate Annual Report 2008. pp. 38–40, June
2008. 9
166
References References
[Kang 07] T. H. Kang and R. A. Iltis. Iterative Decoding, Offset and Chan-
nel Estimation for OFDM using the Unscented Kalman Filter.
Proceedings of the IEEE Asilomar Conference on Signals, Systems
and Computers (ACSSC ’07), pp. 1828–1732, November 2007. 54,
55
167
References References
168
References References
Source and Channel Coding, April 2006. 4, 26, 40, 41, 42, 43, 102,
104, 112
[Mito 00] J. Mitola III. An Integrated Agent Architecture for Software Defined
Radio. PhD thesis, under the supervision of G. Q Maguire, Royal
Institute of Technology (KTH), May 2000. 3, 9
[Mito 99] J. Mitola III and G. Q. Maguire Jr. Cognitive Radio: Making Soft-
ware Radios More Personal. IEEE Personal Communications, Vol. 6,
No. 4, pp. 13–18, August 1999. 8, 9
[More 04] M. Morelli. Timing and Frequency Synchronization for the Uplink
of an OFDMA System. IEEE Transactions on Communications,
Vol. 52, No. 2, pp. 296–306, February 2004. 2, 47, 49
169
References References
170
References References
[Pove 09a] H. Poveda, G. Ferré, and E. Grivel. Estimation Aveugle Itérative par
Filtrage H∞ des Décalages Doppler d’un Système OFDMA Entrelacé
en Liaison Montante. Actes de la Conférence Groupe de Recherche et
d’Etudes du Traitement du Signal (GRETSI ’09), September 2009.
In French. 3, 5, 128
[Pove 10] H. Poveda, G. Ferré, and E. Grivel. Joint Channel and Frequency
Offset Estimation Using Sigma-Point Kalman Filter for an OFDMA
Uplink System. Proceedings of the EURASIP European Signal Pro-
cessing Conference (EUSIPCO ’10), pp. 855–859, August 2010. 4,
54, 86, 127
[Pove 11b] H. Poveda, G. Ferré, and E. Grivel. Frequency Offset and Channel
Estimations for an OFDMA Uplink Primary User System Despite
a Cognitive Radio Interference. To be submitted IEEE Transactions
on Wireless Communications, 2011. 3, 5, 80, 129
171
References References
172
References References
[Sari 98] H. Sari and G. Karam. Orthogonal Frequency Division Multiple Ac-
cess and its Application to CATV Networks. European Transactions
on Telecommunications, Vol. 45, pp. 507–516, November 1998. 24
[Schu 05] H. Schulze and C. Lüders. Theory and Applications of OFDM and
CDMA. Wiley, 2005. 21
[Sesi 09] S. Sesia, I. Toufik, and M. Baker. LTE, The UMTS Long Term
Evolution: From Theory to Practice. Wiley, 2009. 22, 26
173
References References
[Shen 99] X. Shen and L. Deng. A Dynamic System Approach to Speech En-
hancement Using the H∞ Filtering Algorithm. IEEE Transactions
on Speech and Audio Processing, Vol. 7, pp. 391–399, July 1999. 5,
135
174
References References
[Ther 92] C. W. Therrien. Discrete Random Signals and Statistical Signal Pro-
cessing. Prentice-Hall, 1992. 51
[VdMe 01] R. van der Merwe and E. A. Wan. The Square-Root Unscented
Kalman Filter for State and Parameter-Estimation. Proceedings of
the IEEE International Conference on Acoustics, Speech, and Signal
Processing (ICASSP ’01), Vol. 6, pp. 3461–3464, May 2001. 152
[VdMe 04a] R. van der Merwe. Sigma-Point Kalman Filters for Probabilistic
Inference in Dynamic State-space Models. PhD thesis, under the
supervision of E. A. Wan, OGI School of Science and Engineering at
Oregon Health and Science University, April 2004. 4, 56, 134
[VdMe 04b] R. van der Merwe, E. A. Wan, and S. J. Julier. Sigma-Point Kalman
Filters for Nonlinear Estimation and Sensor-Fusion: Applications to
Integrated Navigation. Proceedings of the AIAA Guidance, Navi-
gation and Control Conference (GNC ’04), pp. 1–30, August 2004.
153, 154
[Wan 01] E. A. Wan and R. van der Merwe. Chapter 7, The Unscented
Kalman Filter. Adaptive and Learning Systems for Signal Process-
ing, Communications and Control, Wiley, 2001. 56, 151, 153
175
References References
176
References References
[Zhao 06] P. K. Zhao, L. L. Kuang, and J. H. Lu. Carrier Frequency Offset Es-
timation Using Extended Kalman Filter in Uplink OFDMA Systems.
Proceedings of the IEEE International Conference on Communica-
tions (ICC ’06), Vol. 6, pp. 2870–2874, June 2006. 2, 47, 49, 72
177