Hindawi Publishing Corporation EURASIP Journal on Advances in Signal Processing Volume 2010, Article ID 741429, 18 pages doi:10.1155/2010/741429

Research Article

Pilot-Based Synchronization and Equalization in Filter Bank Multicarrier Communications

Tobias Hidalgo Stitz, Tero Ihalainen, Ari Viholainen, and Markku Renfors (EURASIP Member)

Department of Communications Engineering, Tampere University of Technology, P.O. Box 553, 33101 Tampere, Finland Correspondence should be addressed to Tobias Hidalgo Stitz, tobias.hidalgo@tut.fi Received 16 June 2009; Revised 21 October 2009; Accepted 30 December 2009 Academic Editor: Pierre Siohan

Copyright © 2010 Tobias Hidalgo Stitz et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This paper presents a detailed analysis of synchronization methods based on scattered pilots for filter bank based multicarrier (FBMC) communications, taking into account the interplay of the synchronization, channel estimation, and equalization methods. We show that by applying pilots designed specifically for filter banks, the carrier frequency offset (CFO), fractional time delay (FTD), and channel response can be accurately estimated. Further, a novel joint FTD and channel estimation scheme, based on iterative interference cancelation, permits extending the FTD estimation range well beyond the limit imposed by the pilot separation. The channel parameter estimation and compensation are successfully performed totally in the frequency domain, in a subchannel-wise fashion, which is appealing in spectrally agile and cognitive radio scenarios. The performance evaluation is done in a hypothetical WiMAX scenario in which an FBMC system would substitute OFDM maintaining as much physical layer compatibility as possible.

1. Introduction

Research on increasing data transmission rates and the development of applications exploiting these improved rates, or requiring even more bandwidth, have been fueling each other for decades already, and it seems that this trend will continue. The challenge is therefore to boost high data rates in wireless communications.

Multicarrier (MC) techniques have proven to be a means to adequately overcome many challenges of wide bandwidth transmission while providing also high spectral efficiency. Using MC communications, a frequency selective channel can be divided into several parallel subchannels with flat or mildly selective fading, facilitating channel estimation and equalization (the terms "subchannel", "subband" and "subcarrier" are used interchangeably throughout the whole text, referring to the frequency subband centered at the subcarrier frequency). Timing synchronization is also easier and limited narrowband interference can be easily mitigated.

In fact, the division of the whole bandwidth into many subchannels provides scalability and flexibility when configuring the communication link [1].

The flagship of MC techniques is orthogonal frequency division multiplexing (OFDM) with cyclic prefix (CP) [1]. A correctly chosen CP elegantly turns severely frequency selective channels into flat fading ones at subcarrier level, enabling very simple subcarrier-wise equalization (one complex coefficient per subcarrier). OFDM has become the MC technique of choice in a number of communication systems, including emerging cellular standards such as WiMAX [2] (based on the IEEE 802.16e [3] standard) and E-UTRA LTE (Evolved UMTS Terrestrial Radio Access Long Term Evolution) [4]. These two exploit the flexibility of OFDM to also provide multiple access in time and frequency through orthogonal frequency division multiple access(OFDMA).

Filter bank based multicarrier (FBMC) modulation can be considered an evolved OFDM due to the advantages it presents and also because it can easily be implemented

building on the core blocks of OFDM, that is, the IFFT/FFT pair [5]. The filter banks (FBs) display the following two highly desirable properties. Firstly, their subchannels can be optimally designed in the frequency domain to have desired spectral containment [6]. Secondly, FBMC systems do not make use of the CP, which is pure redundancy, enabling a more efficient use of channel resources. Basically, the subchannel filters are designed with the Nyquist pulse shaping principle, which means that the consecutive symbol waveforms are overlapping in time. Compared to OFDM, where the adjacent subbands are only attenuated 13 dB, the excellent spectral containment allowed by the FBs is crucial for avoiding distortion from asynchronous signals in adjacent bands. Further savings in spectral resources appear at the edges of the transmission band, where ideally only one subchannel can be used as guard band to the next transmission band. An OFDM system with slowly decaying sidelobes would need more subbands or very sharp additional filtering to provide similar attenuation outside the transmission band. In the context of cognitive radio [7, 8], where a secondary user scans the spectrum assigned to primary users for a transmission opportunity (spectral sensing), the high spectral definition of FBMC is especially valuable. The performance of Haykin's method of choice for radio scene analysis, the multitaper method [9], can be well approximated using FBs with greatly reduced complexity [10].

In order to synchronize and equalize the channel, there are methods to perform timing, frequency offset and channel estimation. In the MC context, timing estimation involves locating the start of a transmission burst, finding the first multicarrier symbol, and estimating the fractional time delay (FTD), which is a fraction of the multicarrier symbol period. Since OFDM is the most widespread MC technique, many methods for synchronization and channel estimation can be found for it. They are based on scattered pilots and training sequences [11, 12], on exploiting the redundancy of the CP [13] or even blind methods that do not make use of overheads in the signal [12].

Due to the different nature of the FB waveforms, some of the OFDM methods can be applied to FBMC and others, such as the CP-based ones, cannot. Although there exist preamble based [14], training sequence based [15] and blind [16] approaches, we focus our interest here towards scattered-pilot-based approaches. However, when using efficient filter bank implementations, as for example the FBMC/OQAM (offset quadrature amplitude modulation) [5], the application of scattered pilots is not as straightforward as in OFDM due to the complementary interference that a low-rate symbol suffers from adjacent symbols in time and frequency. This has motivated different proposals for generating pilots [17, 18] compatible with the filter bank class of interest. In this paper, we build our synchronization and channel estimation subsystem based on pilots and auxiliary pilots [19] similarly to [18] because of good estimation performance.

Here it is interesting to mention that the authors of [20] also present an approach to perform full synchronization in both FBMC/OQAM and OFDM/QAM systems based on scattered pilots. However, in that publication the time

domain properties of the regularly scattered pilots are exploited and synchronization takes place in the time domain instead of subcarrier-wise, as is our goal in this paper. Both approaches complement each other and could be used together, for example, for performing coarse acquisition and fine tuning or to improve their respective performances.

In order to obtain results of practical relevance, we develop a testbed system that aims to maintain a certain resemblance and even compatibility up to a degree with the WiMAX profile. Also the channel model and distortion parameters are chosen to be typical of WiMAX communications. We use rather basic synchronization parameter estimation, and compensation methods, many of which are known from the literature in the OFDM context, and adapt them to the FBMC/OQAM system model. Our focus is on understanding the interplay of different synchronization and channel equalization functions. We also explore the possibilities to enhance the overall performance by iterating the estimation and compensation tasks. It is important to underline, that all the necessary synchronization, channel estimation and equalization operations take place after the analysis bank at the receiver, at the low sampling rate. This makes it possible to utilize the filter bank itself for efficient implementation of the needed frequency channelization selectivity for all signal processing functions. Further, the channelization can be done in a dynamically adjustable manner, efficiently suppressing immediately adjacent and even (narrow) in-band interference components.

The paper is structured as follows: first, we review the concept of efficient filter bank based multicarrier communications. In Section 2, we describe how to implement the filter banks and perform efficient subchannel-wise equalization. In Section 3, we first discuss the synchronization requirements in spectrally agile radios. Then we present the method for obtaining pilots for FBMC/OQAM and study the effect of the channel distortions on these pilots. In the end of this section, we present solutions for correcting the distortions, including a novel, iterative interference canceling, joint FTD, and channel estimating algorithm. In Section 4, we introduce the design of the WiMAX-like simulation testbed based on FBMC. Section 5 presents the estimation and correction performance of the studied methods. Finally, Section 6 draws the conclusions from this research.

2. FBMC and Efficient Subchannel-Wise Equalization

2.1. Filter Bank Based Multicarrier Communications. Exponentially modulated filter banks (EMFB) [21], modified DFT (MDFT) filter banks [22], and OFDM with offset QAM (OFDM/OQAM or FBMC/OQAM) [5], among others, are complex filter bank structures that can produce complex I/Q baseband signals for transmission, making them suitable for FBMC systems in spectrally efficient radio communications. In FBMC communications, the filter banks are used in the transmultiplexer (TMUX) configuration [23], with the synthesis filter bank (SFB) in the transmitter and the analysis filter bank (AFB) in the receiver. Figure 1 shows the filter banks in this configuration as fundamental part

s[m] »

lM/2 > Bo(z2) 4-

Bi(z2) >

■>Lm/2> bm-i(z2)4

—>§H>

IM-1,n —X8H>

Subch. proc

o,n . ^

Subch. «ôiv proc.

Subch. proc

Re d0,n R2C0

°M-1,,

dM-1,n »

R2CM-1 -t>

OQAM Transform Polyphase p/s

modulation block filtering conversion

S/P polyphase Transform

conversion filtering block

OQAM demodulation (a) (b)

Figure 1: (a) Synthesis and (b) analysis filter banks for complex FBMC transmultiplexer (TMUX) with per-subchannel processing.

of a complete FBMC/OQAM transmission/reception system. This FB technique builds on uniform modulated filter banks [23], in which a prototype filter p[m] of length Lp is shifted in frequency to generate subbands which cover the whole system bandwidth. The output of such a synthesis filter bank can be expressed by

= X X dk,n6knp

k=0 n=-x

,j(2n/M)k(m-n(M/2)-((Lp-1)/2))

0knn = ej(n/2)(k+n) = j

m is the sample index at the output of the SFB (at high rate), M is the number of subchannels in the filter bank, and dk,n are the real-valued data symbols in subchannel k, transmitted at a rate 2/T. The signaling interval is defined as T = 1/A f, where Af is the subcarrier spacing. The pair of symbols dk,n and dk,n+1 can be interpreted as carrying the in-phase and quadrature information of a complex-valued symbol transmitted at rate 1/T. Therefore, the filter bank presented in Figure 1 is critically sampled. The "C2Rk"-blocks indicate the conversion into real-valued data from the real and imaginary parts of the complex-valued input symbols ck,n and can be considered as introducing upsampling by 2. "R2Ck" carries out the inverse operation after the AFB in the receiver, effectively downsampling the signal by 2. In FBMC/OQAM, ckn belongs to a QAM alphabet and the real and imaginary parts are interleaved with a relative time offset

of T/2 (hence offset QAM) and C2R performs the following mapping:

dk,2n =

[Re [ck,n [Im[ck>n], k odd,

{Im[Ck,n], k even Re [ck,n], k odd.

Note that the signs of the sequences in (2)-(4) could be chosen arbitrarily, but the pattern of real and imaginary symbols after multiplication by 9k,n has to follow the above definitions to maintain orthogonality [5].

This type of filter bank pairs can be efficiently implemented using FFT and IFFT of size M aided by polyphase filtering structures. The different parts of the polyphase SFB structure of Figure 1(a) can be better identified by noting that ak[m] = p[m + kM] and rewriting (1) as

s[m] = dknO^nPtcnp

k=0 n = -x

ej(2n/M)km, (5)

|k,n = (-1) ■ e-j(2nk/M)((LP-1)/2).

Here, the factor (-1)kn centers the low-rate output signal of each subchannel k of the analysis filter bank around DC.

The synthesis and analysis banks in Figure 1, between the OQAM modulation and the subchannel processing parts, respectively, can also be understood from the classical FB presentation: each subchannel k with its own up/downsampling and subchannel filters. The synthesis/analysis subchannel

M-1 oo

filters are, respectively, obtained from the prototype filter p[m] as

gk [m] = p[m]ej(2n/M)k(m-(LP-1)/2), fk [m] = Lp - 1 - m],

where m = 0,1,...,Lp - 1 and (■)* indicates complex conjugation.

Further, the length Lp of the prototype filter p[m] depends on the size of the filter bank and the integer overlapping factor K as Lp = KM, where the factor K indicates the number of FBMC/OQAM symbol waveforms that overlap in time. In [24], other lengths of the prototype filter, close to Lp = KM, are explored. High values for K allow more freedom in designing the prototype filter, for example to achieve very high stopband attenuation. On the other hand, it increases the time required for processing each symbol. The prototype filter p[m] can be designed in such a way that the filter bank pair yields perfect reconstruction (PR) of the transmitted data in case of an ideal channel, that is, the received data sequence dk,n equals the transmitted data dk,n (except for the FB processing delay), if there is no additional processing involved. Methods to design PR prototype filters can be found in [25]. However, in practical communication systems, the channel will always introduce some distortion to the signal. Therefore, the design constraints can be somewhat relaxed and the prototype can be optimized to achieve nearly PR (NPR). The trade-off, when comparing with PR designs, is that for prototype filters ofthe same length, NPR designs can achieve higher stopband attenuations, or with fixed stopband attenuations, the NPR prototype filter can be shorter. This happens at the cost of allowing some marginal intersymbol (ISI) and intercarrier (ICI) interference from the filter bank, well under the noise level of the communication channel.

The NPR prototype can be designed using, among others, window-based techniques [26] or the frequency sampling approach [27, 28]. In the studies described here, the prototype is obtained using the latter method, in which the filter impulse response is obtained by the inverse Fourier transform of the desired frequency response at certain frequency locations. The design of such an NPR prototype filter under different optimization criteria has been addressed in [24].

2.2. Subchannel-Wise Equalization in FBMC. In OFDM, as long as the channel delay spread and the possible synchronization errors remain within the cyclic prefix time, equalization can simply be done with a complex coefficient multiplication at subcarrier level. This approach is also applicable to FBMC, if the ratio of channel delay spread in samples and number of subchannels is sufficiently low, since the frequency variation within a subchannel is then small enough that it can be considered flat fading. But as this is not the general case, more effective channel equalization methods have been developed for FBMC. The single complex coefficient is usually considered when the FB waveforms are well localized in time and frequency

domains to limit the effect on consecutive symbols and neighboring subchannels [5, 29, 30]. Longer finite impulse response (FIR) filters as subcarrier equalizers with cross-connections between the adjacent subchannels to cancel the ICI are studied in [31, 32]. To avoid the cross-connections between subchannels, an oversampled receiver filter bank structure with per-subcarrier FIR equalizers can restore the orthogonality of the subcarrier waveforms. This approach is followed in [33-37] and more recently with MMSE equalizers explicitly for FBMC/OQAM in [38].

The authors of this paper have worked on a low-complexity, subcarrier-wise FBMC equalizer using oversampled subcarrier signals [39, 40]. The equalizers were considered earlier for the exponential modulation type of filter banks, which are basically a generalization of the sine and cosine modulated filter banks to complex-valued signals. Nevertheless, these equalizers can readily be applied to FBMC/OQAM receivers by noting that now the filter bank is even stacked, as compared to the odd stacking of the EMFB-based system. We base the equalization in the current study on this approach because it yields good equalization performance with practical channels and subcarrier spacings of up to 100 kHz, which is more than enough for the WiMAX-like system under consideration.

The working principle is based on frequency sampling: assuming a roll-off factor of the prototype filter a = 1 or smaller, each subchannel overlaps only with the immediately neighboring subchannels. At the oversampled part of the receiver bank, before taking the real part of the subcarrier signals, the equalizer can perform equalization at a number offrequency points according to its complexity. For example, if the equalizer is a 3-tap complex FIR filter, 3 frequency points within the subchannel can be completely equalized, according to the zero-forcing (ZF) or the mean squared error (MSE) criterion. With the filter bank structure in Figure 1, all the subchannels alias to frequencies centered around DC, and a straightforward choice is to equalize at DC and at ±n/2, that is, the center of the subchannel and the passband edge frequencies, respectively. Also other frequencies can be used, as well as longer filters, but the described solution is the computationally least demanding.

In the 3-tap complex FIR equalizer of subchannel k at time n with (noncausal) transfer function

Wkn(z) = w(k-n1]z + wk0) + wkin z-1,

„(1),

the filter coefficients w^with d = -1,0,1, can be tuned in a way that the filter achieves at the mentioned frequency points the following target values:

Here, i <E {0,1,2}, where i = 0 corresponds to the lower subband edge, i = 1 to the subband center, and i = 2 to the upper subband edge frequency. Hk(i,)n is the channel frequency response in subchannel k and time n at frequency position given by i, and y and £ are scaling factors. If the equalizer

applies the ZF criterion, £ = 0 and y = 1. In the MSE case, £ is the noise-to-signal power ratio and the choice of

X2=o( |Hk,n 1 /( |Hk,n

removes the bias of the MSE solution.

Taking into account these assumptions, the equalizer coefficients can readily be derived from the target values in (9) above as [39]

-Xk0n(1 - j) +2xk1n -Xk2n(1+ j)

(0) (2) (0) _ xk,n + xk,n k,n 2 ,

xk0)( 1 + j) +2xk1) -xk2)(1 - j)

Furthermore, it is worth mentioning here that this equalizer can perform also limited FTD and CFO correction in addition to channel equalization, as will be discussed in the next section.

3. Pilot-Based Synchronization in FBMC

3.1. Synchronization in Spectrally Agile Radios. Traditional wireless communication systems are characterized by dedicated frequency bands and well-defined frequency channels. After analog and digital receiver front-ends, the signal contains only the transmissions allocated to that channel, which have well-controlled dynamics under the radio resource management functionalities of the wireless network.

The scarcity of frequency spectrum which can be used in wireless communications is a significant factor that has given raise to the concepts of flexible dynamic spectrum use and cognitive radio. In this context, the used frequency band is not anymore dedicated to a specific service and specific waveforms. The band cannot be expected anymore to be free of other waveforms, the utilized frequency spectrum may be noncontiguous, and the dynamics of signal power levels are not well-controlled anymore. This calls for efficient means to dynamically separate the used portions of the frequency spectrum from other portions that are considered interference. Filter banks are ideal for this purpose. It should be noted that in OFDM, the plain FFT processing does not provide effective filtering to signal elements that are not synchronized to the CP structure, and the frequency channelization selectivity has to be implemented in the analog and digital front-end.

As for the synchronization functionalities in case of spectrally agile radios, it is clear that synchronization parameter estimation cannot be implemented in time domain before major part of the selectivity is implemented and strong interferers are suppressed. As the filter bank itself can be used effectively for implementing the selectivity, the feasibility of time domain synchronization becomes questionable. Therefore, in the context of cognitive radio,

there is a strong motivation to develop synchronization methods which are operating in frequency domain, utilizing the subchannel signals only. However, for compensation of coarse synchronization errors, time domain methods are still favorable. In the following discussion, it is assumed that coarse timing and frequency offsets have been compensated in time domain. The required accuracy of coarse CFO and FTD estimates is an outcome ofthis study.

The primary use ofthese methods is for channel tracking. In normal tracking mode, with continuous flow of data packets, only small CFO and FTD values are expected, and it is enough to use only the estimation algorithms in conjunction with basic time domain compensation methods and some filtering to reduce the random variations of block-wise estimates. However, in advanced packet-based radio interfaces, there can be long gaps in the packet flow, especially when the terminal is in idle mode. This may result in significant drift of the CFO and FTD values, and it is advantageous to be able to compensate significant synchronization errors right away for the first received data packet.

For initial synchronization, the developed scheme could be a part of a search procedure, where different coarse CFO and FTD values are tested until synchronization can be established.

3.2. Signal Models. In order to perform channel equalization and synchronization with the presented equalizer, it has to be fed with channel estimates that provide information about the channel state and the possible synchronization errors. In this paper, we study only the effect of a linear multipath channel with additive white Gaussian noise (AWGN), fractional time delay, and carrier frequency offset. In this scenario, the baseband signal model at the receiver input can be expressed as

r(t) = (s(t) * h(t, t) * S(t - tftdT))ej2n(e/T)t + n(t),

where * represents the convolution operation, s(t) is the continuous time version of (1) (or (5)), h(t, t) is the time-varying transmission channel, and S is the Dirac delta. Moreover, tftd is the fractional time delay as a fraction of the signalling interval T, e is the CFO as a fraction of the subcarrier spacing Af, and n(t) is complexvaluedAWGN. At the receiver, r(t) is sampled at Ts = T/M into r[m] and then passes the analysis bank. Before the subchannel processing stage (estimation and synchronization/equalization in our case), the kth subchannel sequence yk,n can be expressed as

yk,n = [r[m] * fk [m]] M/2 = X vi,n * qi,k,n + nk,n, (13)

qi,k,n

gi[m]*^[n, m]*8

m-T ftdT

,j2n(£/M)m

*fk [m

Table 1: Interference weights for data surrounding the symbol in the center of the table

n = -3

n = -1

k = -2 k = -1 k=0 k = 1 k=2

0.0054 0

0.0054 0

0.0006 j0.0429 -0.0668 -j0.0429 0.0006

-0.0001 -0.1250 0.0002 -0.1250 -0.0001

-j0.2058

0.2393

j0.2058

0.5644

0.5644

j0.2058 0

0.2393 0

-j0.2058 0

-0.0001 -0.1250 0.0002 -0.1250 -0.0001

0.0006 -j0.0429 -0.0668 j0.0429 0.0006

0.0054 0

0.0054 0

Above, qj,k,n is the subchannel-dependent 2-dimensional impulse response, including the channel effects, from subchannel i to subchannel k, and (see Figure 1)

vi,n 6i,ndi,n

The discrete-time time varying channel is h[n, m], the delay tftd( T/Ts) = tftdM is assumed in this paper to be an integer number of samples, for simplicity, and M/2 represents the downsampling by M/2 of the preceding expression. The sampled and filtered noise is nk,n.

In an FB system with a sufficiently frequency selective prototype filter and roll-off a < 1, only adjacent subchannels overlap and have an effect on subchannel k of interest. This implies that the sum in (13), which is basically a 2-dimensional convolution, can be limited to run from i = k-1 to i = k + 1:

yk,n = X Vi,n * qi,k,n + nk,n-i=k-1

This also assumes a reasonably small frequency offset e, in which case the overlaps with subchannels i = k ± 2 remain limited.

3.3. Pilots for FBMC. In a multicarrier system with sufficient number of subcarriers, it is just intuitive to obtain information about a doubly selective (in frequency and time) transmission channel by sampling it in frequency and time directions at certain intervals. The samples are obtained by known data symbols (pilots) that are transmitted at given time and frequency locations and from which the channel information at these locations is recovered. This information is extended to cover the whole signal domain by interpolating between the pilots. In OFDM the application and exploitation of the pilots is straightforward: the channel state in subchannel k at time instant n is just the received symbol divided by the transmitted symbol (this implies ignoring the additive noise). In efficient modulated filter banks it is not that simple, as a closer look into (14) reveals. Indeed, if we first remove all the channel effects in that equation, that is, rFTD = 0, e = 0, and h[n, m] = S[n], qj,k,n becomes the time- and frequency invariant 2-dimensional impulse response of the TMUX, relating Vj,n with yk,n (or equivalently dj,n with the received kth subchannel signal

before the operation of taking its real part, In this case, (16) can be expressed as

yk,n = X X vi,iqi,k,n-l + nk,n

i=k-1 l=-x

k+1 x (17)

= Vk,n + X X Vi,iqi,k,n-l + nk,n-

i=k-1 l=-x

(i,l) = (k,n)

Equivalently, we can write

yk,n = 0k,nyk,n

= dk,nvk,n + Qi nX X v.qkn-l + 01 nVk,n

i=k-1 l=-x

(i,l) = (k,n)

= dk,n + X X ii kdi,iqi,k,n-i + Qtnnkn

i=k-1 l=-x

(i,l) = (k,n)

= dk,n + (

CR + juk,n

k,n^k,n-

(18) (19)

The elements nNpR and uk,n within the parentheses of (19) are obtained by, respectively, taking the real and imaginary part of the summation term in (18). They can be considered interference on the desired symbol dk,n. Now, even in a PR TMUX, in which the real-valued term nN,PR = 0, the imaginary-valued interference juk,n does not add up to 0. The orthogonality is only obtained after taking the real part, that is, this interference summation has only imaginary values [5]. In an NPR design, nNPR is a small real-valued contribution, generally well below the level of nk,n, and will be ignored from now on. Another property of the imaginary-valued interference jukn is that it is time-varying, since it depends on the data in the adjacent channels and on the symbols preceding and following dk,n.

For simplicity of notation, we define the noncausal TMUX response at subchannel k and instant n as tk,n. It is assumed that it is normalized at k = 0 and n = 0 as t0,0 = 1. Table 1 presents the interference weights that multiply neighboring symbols in the case of an NPR prototype designed with the frequency sampling method in [28] and with an overlapping factor K = 4. Thus, we obtain at the location (k0, n0) of interest

= 6k0,n0 yk0,n0 = dk0,n0 + juh,n0 +

where the imaginary interference is

uk0,n0 X dk,ntk0-k,n0-n,

(k,n)enk0in0

tkn = Im[ eyKn]

and Qk0>n„ is the set of subcarrier and time indices that are considered to contribute to the interference on the symbol at (k0, n0). Without channel distortions, t and q are related as

qi,k,n tk-i,n.

Table 1 shows that the interference from subchannels not adjacent to the subchannel of interest have negligible effect and that in time direction the interference goes K symbols from the symbol of interest in both directions. With this design, the residual distortion on the real part is less than -65 dB compared to the actual data.

It is clear that a pilot located at subcarrier kp and time np cannot be immediately recovered even if the channel at that subcarrier is a simple complex coefficient Hkp,np (flat fading), because the pilot information and the complex-valued, time-varying, and data-dependent interference will be mixed by the channel. In the following paragraphs we describe the method we use to handle the described complex-valued interference in order to make pilot-based channel sampling feasible and briefly point out several alternative approaches. The interested reader is referred to the cited references for detailed description of the methods. We first consider a channel with neither FTD nor CFO.

3.3.1. Auxiliary Pilots. This approach is based on an elegant idea presented in [18] utilizing an auxiliary pilot located at ka, na adjacently to the pilot kp, np and which cancels the interference Ukp,np of (21). The advantage of eliminating the interference is that the pilots can be used at the receiver in a similar fashion as OFDM pilots are used; the estimate of the channel at the pilot location is

Hkp,np =

®kp,npykp,np

dkp ,n

assuming that the channel is constant over the whole subchannel bandwidth. The auxiliary pilot is indeed one element of the sum (21) and has to be calculated online with data transmission every time a pilot is inserted since the interference terms vary with the data. In [18], the auxiliary pilot is chosen to cancel the interference from the 8 surrounding symbols. This is a good approximation in well time-frequency localized prototype filters, such as the ones based on the isotropic orthogonal transform algorithm (IOTA) function [29]. Here we chose to include more terms in the calculation of the interference because also symbols located further away can add significant interference. The auxiliary pilot that cancels the imaginary interference can be calculated as

dka,na

tkp-ka,n

kp-ka,np-na (k,n)enkp,„p

(k,n) = (kp,np)

(k,n) = (ka,na)

In typical filter bank designs it is wise to locate the auxiliary pilot immediately preceding or following the pilot, that is, na = np - 1 or na = np + 1. This way, the absolute value of the denominator tkp-ka,np-na is maximized (see Table 1) and the magnitude of the auxiliary pilot is minimized on the average, wasting less transmission energy on the pilot/auxiliary pilot pair and preventing possible strong effects on the peak-to-average-power ratio (PAPR) due to excessively strong auxiliary pilots. As an example, with the prototype filter presented above, it is sensible to use the shaded area of Table 1 for the computation of the interference and the auxiliary pilot. This leaves the residual imaginary interference below -38 dB level with respect to the data. Further, if the auxiliary pilot is chosen to precede or succeed the pilot, then the auxiliary pilot power is on the average 3.3 dB stronger than the data surrounding the pilot/auxiliary pilot pair. Note that the use of two real-valued symbols as pilots does not mean a penalty in overhead with respect to OFDM, since there the pilot is complex-valued.

3.3.2. Alternative Pilot Techniques for Channel Estimation. The pair of real pilots (POP) [17] method also uses two consecutive OQAM subsymbols to send known pilots, in the simplest case, similar pilots. Ignoring the noise and assuming that the channel remains unchanged during both subsymbols, an equation system yields the equalizer coefficient (the same for both pilots) that restores the pilots to their original phase. Its inverse is the estimate for the channel at the positions of the POP. This method places the computational complexity on the receiver part and has the advantage that it is independent from the prototype filter design, since the interference term is not used explicitly in the equations. However, if the noise is not negligible, it will be enhanced in a random fashion, depending on the data surrounding the pilots, which makes the performance unpredictable and generally worse.

The authors of [17] present also the interference approximation method (IAM). The philosophy here is that most of the symbols surrounding the pilot in time and frequency are fixed and known at the receiver. In this case, the interference Ukp,np can be approximatively calculated with (21), which

leads to the estimate Hkp,np = d^pnpykp,np/(dkp,np + jukp,np).

The more symbols in the shaded area of Table 1 are fixed, the better is the approximation of Ukp,np. Unfortunately, this approach leads to an unacceptable overhead if a good approximation of the interference, necessary for accurate channel estimates, is desired. Therefore, it is practical only in situations in which pilots are packed closely together, for example, in preambles, where the fixed symbols can be simultaneously utilized by nearby located pilots.

More sophisticated methods for scattered-pilot-based channel estimation in FBMC/OQAM are presented in [41]. By means of an orthonormal transformation of the data around the pilot, the imaginary-valued term can be nulled while avoiding the need of an auxiliary pilot with increased power. Nevertheless, it is still necessary to fix one of the surrounding symbols, which cannot be used for data transmission. Another approach in [41, 42] consists of iteratively

approximating Ukpnp at the receiver with help of the demodulated data. At first, it is considered that Ukp,np = 0, and the channel is estimated based on this assumption. The equalization and detection then yield estimates of the data surrounding the pilot, which, in turn, permit calculating a better estimate of Ukp,np. This technique converges in 3-4 iterations.

The methods presented above provide a way of sampling the time-varying frequency response of the transmission channel. The estimates required at the remaining subchannels and time-instants can be obtained by interpolating between the estimates obtained at the pilot locations. Note that for the 3-tap subcarrier equalizer presented in Section 2.2, the interpolation has to return estimates not only at the center frequencies of the subcarriers, but also at the subband edges in the middle of those frequencies. Studying how improved interpolation techniques [12] affect the FBMC channel estimation performance is an interesting research subject but goes beyond the scope of this paper.

3.4. Timing Estimation and Correction. Assuming that the coarse location of a transmission frame is obtained by other means, the effect of not synchronizing exactly with the MC symbols, hence introducing a fractional time delay, can be understood with help of (14). In the frequency domain we can write

Q,,k[e}a

= [ ( (Gi(eja) Hn(eja) e-jWTFTDM) *8(u -2ne)) Fkj]

= [ (G; (ej(w-2ne)) Hn{ej{"-2ne)) e-j(^-2ne)TFTD^ Fk

= [G;(ej(w-2ne)) Hn(ej(w-2ne)) F^ej^e-j(w-2ne)TFTD^ ,

which shows that FTD introduces a phase term that linearly varies with the frequency. This frequency-dependent phase term destroys the orthogonality of the subchannels and causes the appearance ofICI.

For estimation purposes, let us first assume that there is neither channel distortion nor CFO present. In this case, (16) can be approximated as

ykn - e-j2nkTFTD X Vi,n * ti-k^e-i=k-1

ej2nTFTD((k-i)/2) + n

where the term e±j2nTFTD(1/2) at i = k - 1 and i = k + 1, respectively, causes the interference uk,n not to be purely imaginary anymore (after compensation of the common phase rotation e-j2nkTTTD and multiplication by n), hence causing ICI. However, from (26) it seems straightforward to estimate the FTD from a transmitted pilot at kp, np, allowing for some uncertainty caused by the ICI. In practice, the estimate is obtained from the phase difference between two pilots separated by Ak subcarriers, since the signal can have

a constant random phase rotation that is eliminated when calculating the phase difference. Thus,

ykp,np) - ykp+Ak,np)

ykp,np) - ykp+Ak,np) + Ak(n/2)

where Z(-) is the phase of (•) and Ak(n/2) comes from the 6k,n that relate ykn and yk,n. This result limits the carrier separation for the pilots, if unambiguous phase differences are to be calculated. The phase difference remains below n if TFTDAk < 1/2. For example, if pilot subcarriers are separated by Ak = 10, only a delay of ±[M/20\ samples ([■ j rounds to the closest smaller integer) can be estimated without phase ambiguity.

The estimate can be improved by averaging techniques such as least squares linear curve fitting over all pilot pairs that are sufficiently closely placed within a frame. In case of CFO, consecutive symbols within a subcarrier are rotated with respect to each other, as will soon be discussed, forcing the averaging to be independently done over the subchannels for every MC symbol time. The CFO additionally increases the additive distortion. Finally, in presence of the transmission channel, the (flat-fading) channel coefficient within the subchannel of interest also multiplies the transmitted signal, introducing additional rotation. Since in general this rotation is different at different subcarriers, it would be expected that this method cannot be used anymore with practical transmission channels. Nevertheless, if the effect of the FTD on the phase is predominant, the phase difference introduced by the channel coefficients will add up as noise at the end. Consequently, we can consider AHkp +Ak,np ) = ¿(Hkp,np) +2nAkTFTD + np, where is an additive phase term, which depends on how much the channel is correlated in frequency direction at the pilot distance. Again, averaging over a sufficient number of pilots reduces the harmful effect of this additional phase term. In general, the FTD estimation can be seen as a problem of phase slope estimation over the active subcarriers.

To compensate for the FTD, an equalizer can be designed to reverse the effect of the frequency-varying phase. The 3-tap complex FIR equalizer presented in Section 2.2 is well suited for performing this task simply by including the FTD effect on the phase of the channel estimates at the subband edges. When computing these estimates for all subcarriers and time instants, the correct interpolation of the channel from the estimates at the pilots is not a trivial problem in the case of FTD due to the modulo 2n phase ambiguity. For example, if amplitude and phase are interpolated separately, in order to include the FTD estimation for compensating the phase, 2n hops can result in a very unsmooth 2-D channel interpolation. In our studies using very basic linear interpolation we obtained the smoothest results by first eliminating the phase slope from the center frequencies of the pilots. Then, 2-D complex interpolation to obtain the channel estimates for the data symbols both at the center and edge frequencies is performed, and finally the phase slope is

restored to the whole estimation matrix for final calculation of the equalizer coefficients.

3.5. Joint FTD and Channel Estimation Based on Iterative Interference Cancelation. The analysis above shows a clear interplay between the quality of the estimates and the presence of FTD and CFO. Therefore, it is desirable to apply the information of the estimated synchronization parameters to the estimated pilots in order to improve the quality of all the pilot-based estimates. One way to achieve this is by jointly estimating the subcarrier-wise channel coefficients and the phase slope. Here, we present a new approach that utilizes the FTD-induced linear phase slope within each pilot subcarrier. We assume a simplified signal model, where the amplitude response is assumed to be constant within each subchannel and the phase is assumed to be a linear function of frequency. This is motivated by the observation that, in many cases, a 1-tap equalizer is able to equalize the channel quite well in the absence of timing offsets. The joint channel-FTD estimation includes two parts: estimation of the complex channel coefficient for each pilot, and estimation ofthe phase slope jointly for all pilots participating in the estimation window QjE.

Simplifying the notation for the 1-tap equalizer coefficient at k, n in (11) as wk,n = wfl, the FTD-compensating target response of the 3-tap subcarrier equalizer presented in Section 2.2 can be written as

WKn(e-j(n/2)) = wk,nejV,

Wk,n(ej0) = wkn, Wk,n(ej(n/2)) = wk,ne-j f,

where the equalizer phase difference correction between the subchannel center frequency and the edge is

f = -ntFTD-

The subcarrier equalizer in the frequency sampling design can now be written as

Wk,n(z) = [0.5(1 - cos( f) - sin( f))z + cos( f) +0.5(1 - cos(f) + sin(f ))z-1]wk,n,

resulting in the following signal model for the subcarrier equalizer output:

yk,n(f) = cos( f)wk,nyk,n

+ 0.5(1 - cos(f) + sin(f))wk,nyk,n-i (31) + 0.5(1 - cos(f) - sin(f))Wk,nyk,n+1.

If small f is assumed we can approximate (31) as

yk,n(f) - Wk,nyk,n + 2Wk,n(yk,n-1 -yk,n+0. (32)

This can be rewritten as

dk,n + jUk,n = 0k,nyk,n(f)

wk,n6£nyk,n + f Wk,n0ln(yk,n-1 - yk,n+1)

= Wk,nyk,n + fWk,nAyk,n,

where we define

Ok,n(yk,n-1 - yk,n+1)

Using the auxiliary pilot scheme, for the pilot symbols the ideal output is dk,n (the subindices ( ■ )p are ignored here for readability). Now, the channel equalization problem can be formulated as

{f, Wk,n]

= arg min j X I dk,n

f,wk,n [k,nenjE

Wk,nyk,n + fWk,nA ykn

where QjE is the set of subcarrier symbols used in the joint estimation. The idea is to adjust f and Wk,n in such a way that the difference between the equalizer output and the known pilot dk,n is minimized in the least-squares sense. This is a nonlinear optimization problem, which can be solved for example by iterating the following two steps.

(1) Assuming that the phase slope is known from the previous iteration (0 in the beginning), [w'kn] are solved from pilots yk,n + f'-lAykn. Note that here i is the iteration index.

(2) Assuming that [w'kn] are known, the observation is a linear function of f and the optimum f can easily be calculated with the derivative of the expression in brackets in (35) with respect to f and setting the result to 0, yielding

Sk,nenjE dk,nRe[wkAyKn] - \wk\2 R^yk^Ayl^

Sk,nenjE I wk Ayk,n I

which, from (29), yields tFTD = -f/n. Note that Wk is here the average of Wk,n within subchannel k. This process converges typically in 5-10 iterations, depending on the FTD and QjE. Moreover, the estimates ft are quite accurate even for higher values of f, as will be seen in Section 5.

Since this method relies on iteratively canceling the ICI, it will hereafter be referred to as iterative interference cancelation (IIC) approach.

3.6. CFO Estimation and Correction. The effect of carrier frequency offset can be analyzed with help of (12), (14), and (25). We observe that before the AFB, the complex exponential that represents the frequency shift induces a linearly time-varying phase rotation of 2n(e/M) on consecutive high-rate samples. Recall that e is the normalized CFO with respect to the subcarrier separation Af. After the receiver bank and downsampling by M/2, the phase rotation between two consecutive symbols at rate 2/T is ne radians. This leads to ISI from symbols in the same subchannel and ICI from symbols in neighboring subchannels at earlier and later time-instants (the neighboring symbols located at the same timing instant have the same phase rotation as the symbol under study, so they do not contribute to ICI, if only CFO is present). The linearly time-varying phase rotation is the property exploited in the forthcoming pilot-based CFO estimation and permits also to implement the basic compensation scheme.

A further effect of the CFO that can be deduced from (25) is that the subchannel filter at the AFB is not frequency-aligned with the corresponding subchannel filter at the SFB. This distorts the effective impulse response of the synthesis-analysis cascade by a factor

fCFo( eit°) = p(ej2u/M) :

where w is the normalized angular frequency at subchannel sample rate and P(ej2w/M) is the prototype filter frequency response in the filter bank design. In filter banks with linearphase channel filters in the analysis (and synthesis) filter banks and in the zero-phase subchannel processing model, only the magnitude of the distorting frequency response is significant. The described distortion on the frequency response for the prototype filter introduced above is shown in Figure 2 for some selected CFO values.

Finally, there are some effects like aliasing effects close to the subband edges and the appearance of distortion from the subchannel located 2 subcarrier spacings away.

3.6.1. CFO Estimation. Based on the CFO effects described above, pilot-based frequency offset estimation is straightforward. If a pilot at kp, np is followed by another one An samples later at kp, np + An, the phase rotation between them will cover Ap = neAn radians. This yields an estimate

¿(ykp,np) - ¿{ykt,np +An)

¿{ykp,np) - ¿(ykp,np+An) + An(n/2)

where An(n/2) comes from the 6k,n that relate yk,n and yk,n.

The performance of the estimation can be improved by averaging over many pilots to reduce the effect of the distortion. For example, if pilots are all separated by the same

2 0 e d

ga - 2

/ / / ' ' /

/ / ' y

-0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8

Normalized frequency

- e = 0.2

--- e = 0.1

e = 0.05 e = 0.01

Figure 2: Subchannel frequency response distortion due to CFO for normalized frequency offsets e = 0.01, 0.05, 0.1, and 0.2.

number of samples, the following expression [43] can be used for estimating the CFO, while simultaneously weighting the subcarrier pilots according to their powers:

¿(Z(k,n)enCFO j Anyk,n+Anyk,t

This way, stronger subcarriers are favored over weaker ones. Again the subindices (•)p indicating pilot locations are not included in the equation for readability. Here, QCFo is the set of pilot locations used for CFO estimation. Note that the same problem of phase-ambiguity appears as in FTD estimation when eAn > 1.

3.6.2. CFO Compensation. If the CFO is moderate, for example, in CFO tracking scenarios, or the subchannel signal modulation is robust, then the offset compensation can be performed subchannel-wise, as sketched in Figure 3(a) with the multiplication by the complex exponential sequence that undoes the frequency shift: e-jnne. The performance of this CFO correction on a WiMAX-like FBMC system is presented in [19], although there no other synchronization tasks are considered. The frequency shifting correction will obviously not recover the information leaked to the neighboring channel nor get rid of the distortion coming from 2 subchannels away, since the frequency shift is circular when performed at subchannel level. Further, the amplitude distortion due to the subchannel filter misalignments, depicted in Figure 2, is also ignored when shifting the signal back to its original location subchannel-wise.

It is possible, however, to compensate for this amplitude distortion in a simple way with the help of the 3-tap equalizer presented above. Similarly as with the phase slope in the presence of FTD, the low complexity equalizer can approximatively correct the amplitude distortion given by

fk [m] -> |M/2 —xg>-> Wk [n]

,j2nme/M

ejnn(e-e)

Figure 3: Subchannel receiver signal model in the presence of CFO, with CFO compensation and subchannel equalization. (a) Basic model. (b) Equivalent form.

(37). This is achieved by setting the desired target responses to

FcFo(ejw)'

for a = -n/2,0, n/2.

Figure 3 presents the subchannel model with the effect of the first two mentioned distortion types, which are predominating. It includes also the basic subchannel-wise CFO correction that undoes the CFO phase rotation according to the estimated frequency shift e. In the equivalent model (b), the effect of inaccurate CFO estimation is represented by the residual frequency modulation ejnn(e-^.

When the frequency shift is too large, the degradation of the signal advocates for time domain solutions for correcting and possibly also for estimating the CFO. The compensation performance degrades faster than the estimation performance: At sufficiently high CFO, the subcarrier-wise compensation cannot recover enough of the signal of interest that has leaked to the neighboring subchannel and eliminate the interference. This logically affects higher order modulations most. However, since several pilots are used for CFO estimation, it is possible that the CFO estimate is still acceptable at the given CFO level. In this case, CFO estimation may still be done at subcarrier level, but compensation has to take place before the analysis filter bank.

It is also worth underlining that when multipath channel, FTD and CFO concur in a communication link, the low-complexity subband equalizer can compensate for the three distortions at the desired frequency locations by combining (multiplying) the target frequency responses needed for correcting each of the estimated distortions.

3.7. Synchronization Scheme with Iteration. In order to estimate and compensate the effect of CFO and FTD and the interplay between these effects shown in (14) and (25) we propose a receiver to jointly estimate and correct the channel impairments subchannel-wise. First, the pilots are recovered from the received signal and from them, the CFO is estimated and the frequency shift is reverted. Then, the synchronization block estimates the FTD and uses

the obtained information when interpolating the channel between the pilots. Finally, at the equalization stage, the channel estimates at the equalization frequency points are weighted by the frequency-dependent amplitude distortion caused by the CFO (Figure 2) and the 3-tap equalizer then equalizes the signal. It is also possible to iterate this loop, since the 3-tap equalizer is able to remove a lot of distortion from the pilots, allowing for better estimates after the iteration.

4. FBMC for WiMAX

Next we describe the parameters and adjustments required for an FBMC communication system that aims to maintain a certain degree of compatibility with WiMAX specifications. We have aimed to design a system taking as much parameters from the mobile WiMAX specifications [2] and the underlying 802.16e standard [3] as possible. For our setup, we have selected the time division duplexing (TDD) specification for the 10 MHz bandwidth at sampling rate fs = 1/Ts = 11.2 MHz. With a transform size M = 1024, the subcarrier spacing is Af = 10.94 kHz. The frame duration of 5 ms allows for transmission of 47 OFDMA symbols with the cyclic prefixes permitted in the specifications. If the whole frame was to be used for a downlink (DL) transmission, FBMC could fit 53 FBMC symbols in the same time because of the absence of the CP. In a more realistic scenario including the uplink (UL) subframe, one or two FBMC symbols would have to be sacrificed for guard times between the forward and reverse links. WiMAX provides for different data configuration modes. Here, we observe two of the possible configurations: the downlink partial usage of subcarriers (DL-PUSC) and the adaptive modulation and coding (AMC23). PUSC and AMC23 use 840 and 864 active consecutive subcarriers, respectively, and have a null subcarrier at the center ofthe transmission band. Compared to AMC23 in OFDMA without additional filtering, if an attenuation of 40 dB outside the 10 MHz transmission band is desired, the FBMC design with the prototype filter presented above could transmit on around 50 additional subcarriers thanks to its good spectral containment. The configurations also define the size of the transmission slots and how the pilots are located within the frame. Figure 4 shows these pilot configurations for the mentioned cases. Here, each OFDMA symbol corresponds to two consecutive subcarrier samples in the FBMC model.

The auxiliary pilot scheme introduced in Section 3.3.1 can directly be applied to the AMC23 pilot configuration. However, if the auxiliary pilot is to be placed in the advantageous position preceding or succeeding the actual pilot, the proximity of the pilots in DL-PUSC poses a computational problem. As shown in Figure 5, directly allocating the pilot and auxiliary pilot fixes a previous auxiliary pilot inside the range used for calculating the complex interference needed for obtaining the auxiliary pilot of interest. If the contribution of this previous auxiliary pilot is not to be ignored, the calculation requires optimization of the auxiliary pilots over the whole burst duration, taking into account their interdependencies (basically, the first

: <- Frequency ->

..... oooo#ooo»ooooooooo«ooo«o

tooooooooooototooooooooo TimeO o oo»ooo»oooooooooiooo«o * tooooooooooototooooooooo

..... otooooooootoooooooot

ooootooooooootoooooo Time oooooootooooooootooo * otooooooootoooooooot

Veh-A, 64-QAM, 3000 frames

t Pilot O Data

Figure 4: Pilot distribution in WiMAX (a) DL-PUSC and (b) AMC23 configurations.

QAM-symbol index

1 2 3 4 5

P P A P A

' Inconvenient pilot placing

OQAM-symbol subindex 12 12 12 12 12

A A

P A P P

Convenient pilot placing

Figure 5: Alternative pilot and data allocations for FBMC with DL-PUSC-like pilot pattern. P: Pilot. A: Auxiliary Pilot.

auxiliary pilot would need to be calculated, then the next one and so on successively). It is more practical to be able to calculate the interference-nulling auxiliary pilots only from the surrounding data. This can be achieved by switching the positions of the pilot/auxiliary pilot pair every other pilot. The lower part of Figure 5 shows this convenient allocation, where the interference window extends only over data and a known pilot.

5. Simulation Setup and Results

We have tested the filter bank based multicarrier system with WiMAX-like parameters and using the synchronization and channel equalization methods that have been described above. The prototype filter is the NPR optimized design with K = 4 described in Section 3.3. The FBMC signal is sent through a quasistatic channel, if not indicated otherwise, modeled according to the International Telecommunication Union ITU-R Vehicular-A Channel guidelines [44]. This means, that for the duration of a transmission burst, the channel remains constant in time unless a certain mobility is indicated. In this case, the fading varies with time

OFDM, DL-PUSC pilots -A- FBMC, DL-PUSC pilots, 3-tap equalizer -0- OFDM, PCI -e- FBMC, PCI, 3-tap equalizer

Figure 6: Comparison between OFDM-based WiMAX and FBMC-based WiMAX with 3-tap equalizers with respect to Eb/N0. 64-QAM transmission with perfect channel information (PCI) and pilot-based estimation in a quasistatic channel.

according to the mobile velocity. For each simulation, 3000 independent burst transmissions with independent channel realizations are performed. The estimation of the synchronization parameters and the channel state relies on scattered pilots obtained with the auxiliary pilot technique. The combined pilot/auxiliary pilot symbol is boosted on average by 4.5 dB with respect to the data. The overall channel response is obtained from the scattered channel estimates by triangulation-based linear interpolation between the pilots. It is important to note that all the processing is done subcarrier-wise, after the AFB.

Extensive results presented in [45] justify favoring the 3-tap equalizer over the single tap solution. With perfect synchronization in the Veh-A channel model, the differences are hardly visible, but in channel models with longer delay spreads and especially when correcting synchronization errors, the performance advantage of the 3-tap equalizer is overwhelming. With respect to the OFDM-based WiMAX, Figure 6 shows the BER performance in the synchronized case of FBMC and OFDM in a system with around 1000 subcarriers that transmit during a whole frame of 5 ms. The OFDM pilots are boosted by 4.5 dB with respect to the data, that is, by the same boost applied to the FBMC pilot/auxiliary pilot pair. The obtained BER curves are quite similar. The better performance of FBMC with perfect knowledge of the channel (PCI) is due to the power that goes into the OFDM CP of length 1/8 of the OFDM symbol [2]. This advantage is almost completely used up when channel estimation takes place because in OFDM all the pilot power is used but in FBMC part of it is dedicated to the auxiliary pilot, which does not contribute to the actual estimation.

CFO ( 1/subcarrier spacing) (a) BER versus CFO with fixed FTD tftd = 0 and 0.1(T/2)

FTD (1/OQAM subsymbol (T/2))

No marker: basic correction x CFO amplitude correction 0 Additional synchronization iteration

(b) BER versus FTD with fixed CFO £ = 0 and 0.04Af

Figure 7: BER performance for different synchronization options and equalization of a DL-PUSC zone consisting of 4 MC symbols within Vehicular-A channel at Eb/N0 = 20 dB, 64-QAM modulation.

In an actual communication system there are synchronization mismatches, and a receiver has to be able to cope with the joint effect of CFO and FTD. The approach detailed above in Section 4, including the iteration of the synchronization/equalization stage, is applied. In our simulations we have concluded that the described synchronization and equalization chain plus one iteration delivers nearly all the improvement, since the differences in offset estimates with additional iterations are marginal.

The first results in Figure 7 show the BER-performance of the mentioned system with 64-QAM at Eb/N0 = 20 dB in a quasistatic Vehicular-A channel. All the subcarriers are used,

Veh-A, Eb/N0 = 20 dB, 64-QAM, UL-AMC23 slots

+ 4 X 2 slots

^4 x 2 slots, CFO amplitude compensation, 1 iteration * 16 X 8 slots

O 16 X 8 slots, CFO amplitude compensation, 1 iteration

Figure 8: BER performance versus CFO and FTD for uplink transmission of 16 X 8 and 4 X 2 AMC23 slots in Vehicular-A channel at Eb/N0 = 20 dB using 64-QAM. Different synchronization options.

and the pilots are in the DL-PUSC configuration, simulating a downlink scenario. The transmitted burst consists of 4 FBMC symbols. Until otherwise mentioned, the FTD is calculated through averaging the phase differences between pilots. Figure 7(a) presents the BER results with changing CFO for two fixed FTD values, while Figure 7(b) presents the BER as a function of the FTD for two fixed CFO values. These figures show a clear improvement in performance with respect to CFO due to the CFO amplitude compensation. Further, the iteration of the synchronization and equalization stages also improves the BER, especially with respect to rFTD. Because improvements in 4-QAM are not that visible (since the modulation is very robust by itself), these results are not shown here. Nevertheless, a slight improvement can still be achieved.

Next, a situation similar to an uplink is tested. The pilots are now placed according to the AMC23 configuration, as in Figure 4(b), and the user sends its data in AMC23 slots, each slot consisting of 18 consecutive subcarriers and 3 MC symbols. Within each slot there are 6 pilot/auxiliary pilot pairs, as can be seen from the figure. We have studied two sizes for the transmitted bursts: 4 X 2 and 16 X 8 slots, where the first number indicates the number of slots in frequency direction and the second number in time direction. The smaller burst includes 48 pilot/auxiliary pilot pairs, the larger one has 768. We assume that other users are separated by the required guard band and that they do not produce multiuser interference.

Veh-A, Eb/N0 = 14 dB, 4-QAM, UL-AMC23 slots

---16 X 8 slots OFDM

A 4 X 2 slots, CFO amplitude compensation, 1 iteration O 16 X 8 slots, CFO amplitude compensation, 1 iteration

Figure 9: BER performance versus CFO and FTD for uplink transmission of 16 X 8 and 4 X 2 AMC23 slots in Vehicular-A channel at Eb/N0 = 14 dB using 4-QAM. OFDM- and FBMC-based WiMAX. FBMC Synchronization performed with CFO amplitude compensation in the 3-tap equalizer and one synchronization iteration.

Figures 8 and 9 present the BER performances for 64-QAM and 4-QAM, respectively. The degradation that occurs when using a smaller burst size can be observed. This is due to the lower number of pilots available to perform synchronization and estimation. We further see that the range in which the synchronization is successful is much smaller than in the downlink case, but this is mainly due to the larger pilot separation. The effect of this is that the phase ambiguity when estimating FTD or CFO appears for Tftd between 0.05 to 0.1(1/2) and £ between 0.15 to 0.2A/ (at around £ = 0.17, because only every 6th symbol is a pilot in time direction). Figure 9 also compares the BER performance of an OFDM-based WiMAX system using the same synchronization and estimation methods as its FBMC counterpart. It occupies 16 x 8 AMC23 slots and the BER results are similar to the FBMC case. The main difference is that in the CFO axis the performance degradation commences earlier. Due to the CP extension of the OFDM symbols, effectively slightly "separating" the pilots in time, the phase ambiguity appears approximately when £ = 0.14. Further, Figure 10 shows the effects on the FBMC system of user mobility and low SNR on the smaller burst with robust subchannel modulation. The performance is barely affected at this speed and only when the CFO approaches the nonambiguity limit the difference becomes visible. On the other hand, the stronger noise is evidently worsening the BER, but still allowing for a certain synchronization range. The following simulation results prove this by taking a closer

Veh-A, Eb/N0 = 5 and 14dB, 4-QAM, UL-AMC23 4 X 2 slots, CFO amplitude compensation + 1iteration

...... Eb/N0 = 5 dB, 0km/h

-A- Eb/No = 14 dB, 60 km/h O Eb/N0 = 14 dB, 0 km/h

Figure 10: BER performance versus CFO and FTD for uplink transmission of 4 X 2 AMC23 slots in Vehicular-A channel at Eb/N0 = 5 dB and Eb/N0 = 14 dB using 4-QAM. Static user and user with 60 km/h mobility. Synchronization with CFO amplitude distortion correction and one iteration.

look at the estimation performance of the synchronization parameters.

The CFO estimator performance is evaluated with help of the root mean square (RMS) error with respect to the actual CFO. The RMS error of the previously presented simulation scenarios is shown for a fixed tftd = 0.1(772) in Figure 11(a). The CFO estimation performance proves to be quite independent from the actual FTD, having only a minor degradation with increasing tftd, as the interested reader can verify in [45]. In the low SNR scenario, even with the small burst and with FTD present, the estimate for the frequency offset is acceptable if the CFO stays below 7% of the subcarrier spacing. Also the user with vehicular mobility can be synchronized well by the base station. The other considered cases have a low variation of the performance with increasing CFO, until the ambiguity limit is reached. Although not presented here, the RMS CFO error of e increases slightly if the synchronization procedure is iterated. Nevertheless, the estimates are still good enough for accurate synchronization within the WiMAX requirements for maximum CFO offset, specifically, ±2% of the subcarrier spacing [3].

Similar results, now for the FTD RMS estimation error at the fixed CFO £ = 0.1 A/, are summarized in Figure 11(b). Here the estimation error with iteration of the synchronization part is almost overlaying with the presented results. Also the case with the AMC23 mobile at 60 km/h has been omitted for readability, since the estimation RMS error is similar to the other AMC23 results. As in the CFO case, also here the FTD error is very much independent

RMS CFO error, Veh-A, fixed TFTD = 0.1(772)

Veh-A, Eb/No = 14 dB, 4-QAM, UL-AMC23 slots

CFO (1/subcarrier spacing (Д/))

■ ■ о ■ ■ 4 X 2 UL-AMC23 slots, Eb/N0 = 5dB, 0 km/h -a- 4 X 2 UL-AMC23 slots, Eb/N0 = 14dB, 60km/h -е- 4 X 2 UL-AMC23 slots, Eb/N0 = 14 dB, 0km/h

- 4 DL-PUSC MC symbols, Eb/N0 = 14 dB, 0 km/h

-*- 16 X 8 UL-AMC23 slots, Eb/N0 = 14 dB, 0 km/h

(a) RMS CFO error versus CFO at fixed FTD rFTD = 0.1(1/2)

RMS FTD error, Veh-A, fixed £ = 0.1Д/

FTD (1/OQAM subsymbol (T/2))

-в- 4 X 2 UL-AMC23 slots, Eb/N0 = 14 dB, 0km/h -*- 16 X 8 UL-AMC23 slots, Eb/N0 = 14 dB, 0 km/h

■ ■ о ■ - 4 X 2 UL-AMC23 slots, Eb/N0 = 5dB, 0 km/h -4 DL-PUSC MC symbols, Eb/N0 = 14 dB, 0km/h

■ - OFDM, 16 X 8 UL-AMC23 slots, Eb/N0 = 14 dB, 0 km/h

(b) RMS FTD error versus FTD at fixed CFO £ = 0.1Д/

Figure 11: RMS errors in CFO and FTD estimation in AMC23 and DL-PUSC transmission with different burst sizes. Eb/N0 values are referred to 4-QAM modulation.

of the CFO (see again [45]) and is almost constant within the nonambiguous estimation range, between 5 and 10 high rate samples. This error floor for the different scenarios is caused by the inherent delay the channel introduces due to the delay spread distribution. Its effect is a mild slope in the phase with respect to the subcarrier index, even in the

W 10"1 я 10

10"2 0.5

9-%;2 0.1 0.1 0

Conventional FTD estimation, 4 X 2 UL-AMC23 slots ■■*■■ Conventional FTD estimation, 16 X 8 UL-AMC23 slots о IIC, 4 X 2 UL-AMC23 slots * IIC, 16 X 8 UL-AMC23 slots

Figure 12: BER performance versus CFO and FTD of the IIC joint FTD and channel estimation method compared to classical FTD estimation. 4-QAM transmission of 4 X 2 and 16 X 8 AMC23 slots in Vehicular-A channel at Eb/N0 = 14 dB. Static user. No additional synchronization iteration.

absence of proper FTD. This mild slope is also detected by the estimation methods as part of the FTD, inducing the error floor when comparing it to the actual FTD.

Figures 11(a) and 11(b) include the RMS estimation error performance of the OFDM-based system. It can be seen that due to the higher available power boost for the OFDM pilots, the CFO estimates are better than in the FBMC system in the nonambiguous region. The FTD estimation error almost matches the error in the FBMC system because the mild slope of the inherent channel delay mentioned in the previous paragraph is detected here, too. The better CFO estimation performance does not translate into a better BER performance, as can be verified above in Figure 9.

The FTD estimation range can be extended further than the limit imposed by the frequency separation ofthe pilots by applying the interference minimization approach described in Section 3.5, yielding (36). The BER performance of the IIC technique for the AMC23 transmission configuration is shown in Figure 12, where the novel method is iterated 10 times. Figure 13 presents the corresponding RMS FTD estimation error. The extension into longer delay operating ranges is evident, and the BER performance for this Eb/N0 achieves the same level as the conventional technique when the larger bursts are compared. However, when a mobile user transmits using small bursts, and consequently few pilots, the BER and estimation performances of the IIC approach suffer greatly. If the interference cancelation method does not perform well enough in estimating the FTD because it has not enough pilots, it can, nevertheless, be applied for a coarse estimation within a wide FTD range to pin the delay

Veh-A, Eb/N0 = 14 dB, 4-QAM, UL-AMC23 slots

■* ■ Conventional FTD estimation, 4 X 2 and 16 X 8 UL-AMC23 slots O IIC, 4 X 2 UL-AMC23 slots * IIC, 16 X 8 UL-AMC23 slots

Figure 13: RMS error in FTD estimation versus CFO and FTD of the IIC FTD estimation method compared to classical FTD estimation. 4-QAM transmission of 4 X 2 and 16 X 8 AMC23 slots in Vehicular-A channel at Eb/N0 = 14 dB. Static user. No additional synchronization iteration.

to a value from where the conventional FTD estimation can take over the task.

Figure 14 presents the performance of the IIC method for AMC23 transmission of 4 x 2 slots under different channel conditions. We see that this method is more sensitive to the mobility conditions than the conventional FTD and channel estimation method (see Figure 10). The BER performance with low SNR is also worse in this small burst case. If only 5 iterations in the IIC are performed, then the BER worsens for longer fractional time delays. Although not presented here, 5 iterations suffice if the larger burst of 16 x 8 slots is sent, since the performance is practically the same as with 10 iterations.

6. Conclusions

In this contribution we proposed an integrated synchronization subsystem for filter bank based multicarrier communications. If the dimensioning of the filter bank leads to approximately flat fading subchannels, scattered pilots can be used to perform the synchronization and channel estimation tasks. Spectrally efficient filter banks with efficient implementations, such as the FBMC/OQAM transmultiplexer, can provide well-contained communication channels with high isolation to adjacent, nonsynchronous signals. However, the orthogonality conditions between subchannels for this kind of systems require special solutions in pilot design. The auxiliary pilot approach permits creating pilots that can take advantage of the basic synchronization procedures for OFDM to estimate the fractional time delay, carrier offsets, and communication channel. More important, it enables the design of FBMC-specific synchronization and channel estimation methods taking into account the interplay of the

Veh-A, Eb/No = 5 and 14 dB, 4-QAM, UL-AMC23 4 X 2 slots, CFO amplitude compensation

......Eb/N0 = 5 dB, 0 km/h, 10 IIC iterations

x Eb/N0 = 14 dB, 0 km/h, 5 IIC iterations -A- Eb/N0 = 14 dB, 60 km/h, 10 IIC iterations -Eb/N0 = 14 dB, 0km/h, 10 IIC iterations

Figure 14: BER performance versus CFO and FTD of the IIC joint FTD and channel estimation method in different channel scenarios. 4-QAM transmission of 4 X 2 AMC23 slots in Vehicular-A channel at Eb/N0 = 5dB and Eb/N0 = 14 dB and at different velocities. 5 and 10 iterations in the IIC technique. No additional synchronization iteration.

different stages involved in estimating and compensating for the distortions experienced by the signal. The utilized low-complexity subchannel-wise equalizer for filter banks can use these estimates for equalizing the channel and also to correct distortions caused by incorrect timing or frequency tuning in an efficient way. We have also presented a novel iterative interference cancelation (IIC) approach that estimates the fractional time delay and can return channel estimates by iterating an optimization routine specific to the low-complexity equalizer.

The mentioned synchronization apparatus was put to test in a 10 MHz WiMAX-like communications scenario, where the discussed tools succeed in estimating and compensating the channel distortions. Special emphasis is put on maintaining as many compatible physical layer elements with WiMAX as possible. As a matter of fact, the FBMC system can be updated to send over more subchannels and more MC symbols within the given time frame and system bandwidth, compared to the OFDM-based WiMAX, while achieving similar BER performance. For example, if a whole frame is used for downlink transmission and a 40 dB attenuation is assumed at the transmission band limits, the additional MC symbols and subcarriers in the FBMC-based system permit transmitting 16% more QAM symbols.

Clever reuse of the subchannel equalizer leads to improved BER performance, compared to basic techniques, especially for higher order modulations, for example, 64-QAM. In this case, performance degradation with significant

synchronization offsets is due to insufficient compensation by the synchronization subsystem, whereas with more robust modulations, such as 4-QAM, the bottleneck is in the estimation performance, that is, the performance after compensation is satisfactory as long as estimation is acceptable, which mainly depends on the pilot distribution. In this sense, the introduced joint FTD and channel estimation approach can outperform the conventional method under favorable channel conditions (good SNR and low mobility) allowing acceptable reception under severe FTD. In adverse conditions, it provides sufficiently good estimates for a wide range of FTD, making it suitable for coarse estimation and synchronization of robust data signals.

This can also be exploited when multiple-input multiple-output (MIMO) techniques are used to improve the throughput. The pilots of one transmit antenna have to coincide with "silence" on the other antennas to enable MIMO channel estimation similar as in the single stream case. If the pilot (including the "silent pilots") overhead is to remain the same, as for example in the optional 2 antenna WiMAX AMC transmission [3], the density of active pilots that can be used for actual estimation of the MIMO channel diminishes. In this scenario of greater separation of active pilots, the new FTD estimation with its wider estimation range can prove itself very practical.

Some of the simulation results of Section 5 were obtained for small- and medium-sized transmission bursts. Since the FTD and CFO values are relatively slowly changing in most wireless communication scenarios, it is always possible to refine the estimates by filtering them over multiple transmission bursts to reduce the random variations of block-wise estimates.

An attractive aspect of the presented method is that it operates completely in the frequency domain, after the analysis filter bank. This is especially useful in cognitive radio scenarios, in which high spectral containment and resolution permit secondary users to exploit appearing and vanishing spectral time-frequency holes in the communications of the primary system. Processing after the AFB means that the signal has been well delimited in frequency domain.

In an uplink multiuser scenario, a base station using the presented synchronization, estimation, and equalization methods for FBMC can independently synchronize different users in different bands transmitting in an unsynchronized way, that is, with different FTDs, as long as they are separated by a guard band of at least one subchannel. The bands assigned to each user are processed subcarrier-wise and independently. It is not necessary for the users to transmit synchronously because the guard subchannel assures that the signals remain orthogonal to each other, independently of their relative timings.

A possible direction of future work is to examine the behavior of smaller filter banks in the same bandwidth, with higher subcarrier spacings. This situation poses challenges to the estimation part, since the use ofscattered pilots is not that well justified when the subchannel channel is not flat fading anymore.

All in all, we conclude that FBMC can perform sufficiently well to offer an alternative to OFDM as modulation

and multiple access technique, offering several advantages especially as a cognitive radio physical layer.

Acknowledgments

This research was supported in part by the European Commission under Project PHYDYAS (FP7-ICT-2007-1-211887). The authors would like to acknowledge the valuable contributions of their colleagues in the PHYDYAS project. They also wish to thank the reviewers and the editor for their comments that helped to improve the manuscript.

References

[1] R. van Nee and R. Prasad, OFDM for Wireless Multimedia Communications, Artech House, Boston, Mass, USA, 2000.

[2] WiMAX Forum, "Mobile system profile," Release 1.0 Approved Specification, Revision 1.5.0, November 2007.

[3] "IEEE standard for local and metropolitan area networks. Part 16: air interface for fixed and mobile broadband wireless access systems. Amendment 2: physical and medium access control layers for combined fixed and mobile operation in licensed bands and corrigendum 1," IEEE Std 802.16e-2005 and IEEE Std 802.16-2004/Cor 1-2005 (Amendment and Corrigendum to IEEE Std 802.16-2004), pp. 1-822, 2006.

[4] TSG-RAN 3GPP TS 36.211, "Physical channels and modulation," Release 8, February 2009.

[5] P. Siohan, C. Siclet, and N. Lacaille, "Analysis and design of OFDM/OQAM systems based on filterbank theory," IEEE Transactions on Signal Processing, vol. 50, no. 5, pp. 1170-1183, 2002.

[6] H. S. Malvar, Signal Processing with Lapped Transforms, Artech House, Boston, Mass, USA, 1992.

[7] J. Mitola III and G. Q. Maguire Jr., "Cognitive radio: making software radios more personal," IEEE Personal Communications, vol. 6, no. 4, pp. 13-18, 1999.

[8] S. Haykin, "Cognitive radio: brain-empowered wireless communications," IEEE Journal on Selected Areas in Communications, vol. 23, no. 2, pp. 201-220, 2005.

[9] D. Thomson, "Spectrum estimation and harmonic analysis," Proceedings of the IEEE, vol. 70, no. 9, pp. 1055-1096, 1982.

[10] B. Farhang-Boroujeny and R. Kempter, "Multicarrier communication techniques for spectrum sensing and communication in cognitive radios," IEEE Communications Magazine, vol. 46, no. 4, pp. 80-85, 2008.

[11] J. H. Manton, "Optimal training sequences and pilot tones for OFDM systems," IEEE Communications Letters, vol. 5, no. 4, pp. 151-153, 2001.

[12] M. Pun, M. Morelli, and C. Jay Kuo, Multi-Carrier Techniques for Broadband Wireless Communications. A Signal Processing Perspective, Imperial College Press, London, UK, 2007.

[13] J.-J. van de Beek, M. Sandell, and P. O. Borjesson, "ML estimation of time and frequency offset in OFDM systems," IEEE Transactions on Signal Processing, vol. 45, no. 7, pp. 18001805, 1997.

[14] T. Fusco, A. Petrella, and M. Tanda, "Data-aided symbol timing and CFO synchronization for filter bank multicarrier systems," IEEE Transactions on Wireless Communications, vol. 8, no. 5, pp. 2705-2715, 2009.

[15] T. Hidalgo Stitz, T. Ihalainen, and M. Renfors, "Practical issues in frequency domain synchronization for filter bank based multicarrier transmission," in Proceedings of the 3rd International Symposium on Communications, Control, and

Signal Processing (ISCCSP '08), pp. 411-416, St. Julians, Malta, March 2008.

[16] H. Bolcskei, "Blind estimation of symbol timing and carrier frequency offset in wireless OFDM systems," IEEE Transactions on Communications, vol. 49, no. 6, pp. 988-999, 2001.

[17] C. Lélé, J.-P. Javaudin, R. Legouable, A. Skrzypczak, and P. Siohan, "Channel estimation methods for preamble-based OFDM/OQAM modulations," in Proceedings of the European Wireless Conference, Paris, France, April 2007.

[18] J.-P. Javaudin, D. Lacroix, and A. Rouxel, "Pilot-aided channel estimation for OFDM/OQAM," in Proceedings of the 57th IEEE Semiannual Vehicular Technology Conference (VTC '03), vol. 3, pp. 1581-1585, April 2003.

[19] T. Hidalgo Stitz, A. Viholainen, T. Ihalainen, and M. Renfors, "CFO estimation and correction in a WiMAX-like FBMC system," in Proceedings of IEEE Workshop on Signal Processing Advances in Wireless Communications (SPAWC '09), pp. 633637, Perugia, Italy, June 2009.

[20] B. Jahan, M. Lanoiselee, G. Degoulet, and R. Rabineau, "Full synchronization method for OFDM/OQAM and OFDM/QAM modulations," in Proceedings of IEEE International Symposium on Spread Spectrum Techniques and Applications (ISSSTA), pp. 344-348, Perugia, Italy, June 2008.

[21] J. Alhava and M. Renfors, "Exponentially-modulated filter bank-based transmultiplexer," in Proceedings of IEEE International Symposium on Circuits and Systems (ISCAS), vol. 4, pp. 233-236, Bangkok, Thailand, May 2003.

[22] S. Mirabbasi and K. Martin, "Overlapped complex-modulated transmultiplexer filters with simplified design and superior stopbands," IEEE Transactions on Circuits and Systems II, vol. 50, no. 8, pp. 456-469, 2003.

[23] P. P. Vaidyanathan, Multirate Systems and Filter Banks, Prentice-Hall, Englewood Cliffs, NJ, USA, 1993.

[24] A. Viholainen, T. Ihalainen, T. Hidalgo Stitz, M. Renfors, and M. Bellanger, "Prototype filter design for filter bank based multicarrier transmission," in Proceedings of European Signal Processing Conference (EUSIPCO), Glasgow, Scotland, August 2009.

[25] R. Bregovic, Optimal design of perfect-reconstruction and nearly perfect-reconstruction multirate filter banks, Ph.D. dissertation, Tampere University of Technology, TUT Publications 428, Tampere, Finland, 2003.

[26] P. Martin-Martin, R. Bregovic, A. Martin-Marcos, F. Cruz-Roldan, and T. Saramaki, "A generalized window approach for designing transmultiplexers," IEEE Transactions on Circuits and Systems I, vol. 55, no. 9, pp. 2696-2706, 2008.

[27] K. W. Martin, "Small side-lobe filter design for multitone data-communication applications," IEEE Transactions on Circuits and Systems II, vol. 45, no. 8, pp. 1155-1161, 1998.

[28] M. G. Bellanger, "Specification and design of a prototype filter for filter bank based multicarrier transmission," in Proceedings of IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP '01), vol. 4, pp. 2417-2420, May 2001.

[29] B. Le Floch, M. Alard, and C. Berrou, "Coded orthogonal frequency division multiplex," Proceedings of the IEEE, vol. 83, no. 6, pp. 829-996, 1995.

[30] A. Vahlin and N. Holte, "Optimal finite duration pulses for OFDM," IEEE Transactions on Communications, vol. 44, no. 1, pp. 10-14, 1996.

[31] S. D. Sandberg and M. A. Tzannes, "Overlapped discrete multitone modulation for high speed copper wire communications," IEEE Journal on Selected Areas in Communications, vol. 13, no. 9, pp. 1571-1585, 1995.

[32] L. Vandendorpe, L. Cuvelier, F. Deryck, J. Louveaux, and O. Van de Wiel, "Fractionally spaced linear and decision-feedback detectors for transmultiplexers," IEEE Transactions on Signal Processing, vol. 46, no. 4, pp. 996-1011, 1998.

[33] B. Hirosaki, "An analysis of automatic equalizers for orthogonally multiplexed QAM systems," IEEE Transactions on Communications Systems, vol. 28, no. 1, pp. 73-83, 1980.

[34] T. Wiegand and N. J. Fliege, "Equalizers for transmultiplexers in orthogonal multiple carrier data transmission," in Proceedings of European Signal Processing Conference (EUSIPCO), vol. 2, pp. 1211-1214, Trieste, Italy, September 1996.

[35] S. Nedic, "An unified approach to equalization and echo cancellation in OQAM-based multi-carrier data transmission," in Proceedings of the IEEE Global Telecommunications Conference, vol. 3, pp. 1519-1523, Phoenix, Ariz, USA, November 1997.

[36] K. van Acker, G. Leus, M. Moonen, O. van de Wiel, and T. Pollet, "Per tone equalization for DMT-based systems," IEEE Transactions on Communications, vol. 49, no. 1, pp. 109-119, 2001.

[37] A. M. Wyglinski, P. Kabal, and F. Labeau, "Adaptive filterbank multicarrier wireless systems for indoor environments," in Proceddings of IEEE Vehicular Technology Conference, pp. 336340, Vancouver, Canada, September 2002.

[38] D. S. Waldhauser, L. G. Baltar, and J. A. Nossek, "MMSE subcarrier equalization for filter bank based multicarrier systems," in Proceedings of IEEE Workshop on Signal Processing Advances in Wireless Communications (SPAWC '08), pp. 525529, July 2008.

[39] T. Ihalainen, T. Hidalgo Stitz, M. Rinne, and M. Renfors, "Channel equalization in filter bank based multicarrier modulation for wireless communications," EURASIP Journal on Advances in Signal Processing, vol. 2007, Article ID 49389, 18 pages, 2007.

[40] M. Renfors, T. Ihalainen, and T. Hidalgo Stitz, "Channel equalization," Patent EP1 716 681, November 2008.

[41] C. Lele, OFDM/OQAM modulation: channel estimation methods, and applications to multicarrier CDMA and multi-antenna transmission, Ph.D. dissertation, Conservatoire National des Arts et Meitiers, November 2008.

[42] C. Lele, R. Legouable, and P. Siohan, "Iterative scattered pilot channel estimation in OFDM/OQAM," in Proceedings of IEEE Workshop on Signal Processing Advances in Wireless Communications (SPAWC '09), pp. 176-180, June 2009.

[43] S. Kay, "Fast and accurate single frequency estimator," IEEE Transactions on Acoustics, Speech, and Signal Processing, vol. 37, no. 12, pp. 1987-1990, 1989.

[44] ITU-R, "Guidelines for evaluation of radio transmission technologies for IMT-2000," Recommendation M.1225, 1997.

[45] INFSO-ICT-211887 Project PHYDYAS, "Deliverable 2.2: synchronization and initialization with single antenna. Blind techniques," January 2009, http://www.ict-phydyas.org/ userfiles/file/Phydyas-D2-2.pdf.