0 EURASIP Journal on

Advances in Signal Processing

a SpringerOpen Journal

RESEARCH

Open Access

Preamble-based frequency-domain joint CFO and STO estimation for OQAM-based filter bank multiearrier

Stijn Van Caekenberghe1, André Bourdoux2, Liesbet Van der Perre2 and Jérôme Louveaux3*

Abstract

Filter bank multicarrier systems, similarly to orthogonal frequency division multiplexing (OFDM), are very sensitive to carrier frequency offset (CFO) and symbol timing offset (STO). In this paper, a low-complexity preamble-based joint CFO and STO technique is presented. It is based on a relatively long preamble in order to improve the CFO estimation performance as well as avoid interference coming from the data following this preamble. After CFO and STO correction, the preamble can be reused to estimate the channel. Unlike most current techniques, the CFO and STO estimation occurs in the frequency domain. This allows for a low-complexity estimation with respect to time-domain techniques and, as will be shown by simulations, provide even better performance in a reasonable range. The drawback however is that the estimation range is shorter. Specifically, for large STOs (and to a smaller extent large CFOs), the performance decreases below time-domain estimations. Two versions of the STO estimation technique will be presented, the second one being an approximation of the first one, making it less complex yet also less precise. The performance is assessed by means of computer simulations, testing for both large and small STOs, and compared with existing techniques.

Keywords: Preamble-based estimation; Carrier frequency offset; Symbol timing offset; Filter bank multicarrier

1 Introduction

Filter bank multicarrier (FBMC) is a family of multicarrier modulation techniques that use discrete Fourier transform (DFT)-modulated filter bank in order to obtain a better spectral containment than the traditional orthogonal frequency division multiplexing (OFDM). There exists different versions of FBMC such as filtered multitone (FMT) and offset QAM (FBMC/OQAM). We focus on the latter in this paper. The FBMC/OQAM offers, at the expense of an increased complexity, several advantages over OFDM. The first one is the gain in spectral efficiency related to the removal of the cyclic extension. But the major advantage is the possibility to have several coexisting systems with very little guard bands, which is a very desirable property in wireless communications where spectrum is expensive and should be used as efficiently as possible. For this reason, FBMC has been strongly

•Correspondence: jerome.louveaux@udouvain.be

3Université catholique de Louvain, Place du Levant, 2, Louvain-la-Neuve

B-1348, Belgium

Full list of author information is available at the end of the article

considered for cognitive applications recently [1] as well as several other applications such as professional mobile radio (PMR) or 5G mobile networks. Just like OFDM, FBMC is highly sensitive to carrier frequency offset (CFO) and symbol timing offset (STO). A good estimation and correction technique is therefore essential.

There has been a lot of literature on CFO and STO estimation for OFDM, but most of these techniques cannot be directly applied to FBMC/OQAM due to the removal of the cyclic extension and due to the particular structure of the OQAM. Hence, a good amount of research has been devoted recently to specific techniques for the synchronization in FBMC/OQAM. The literature focused on blind estimation methods initially. In [2], a blind joint CFO and STO estimation has been presented based on the cyclostationarity of the FBMC/OQAM signal. In [3], the CFO estimation is further improved by using the conjugate second-order cyclostationarity statistics. Then, a frequency-domain implementation is proposed in [4]. In [5], a blind CFO estimator is obtained based on the

Springer

© 2014 Van Caekenberghe et al.; licensee Springer. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original workis properly credited.

maximum likelihood (ML) principle for low signal-to-noise ratio (SNR). In [6], a blind closed-loop method is proposed for the tracking of the STO based on the ML estimation, with several approximation to obtain a computationally efficient algorithm. More recently, pilot-based (or preamble-based) synchronization has received more attention. In [7], a periodic preamble is considered, and both STO and CFO estimators are designed based on a least-square (LS) approach. This is a timedomain approach and exhibits a stable performance independently of the actual STO, making it very well suited as a coarse alignment algorithm. It also provides good robustness against multipath channels but has rather high complexity. In [8], the same authors develop a joint STO and CFO estimator for short nonperiodic preamble based on the ML principle. A closed-form approximate expression of the CFO estimation is presented that provides accurate performance for moderate values. Another timedomain technique that should be mentioned is presented in [9] for FMT. And in [10], LS-based CFO and STO estimation is investigated for a short preamble, designed specifically for low latency and simplified channel estimation. In [11], a CFO estimation is derived for scattered pilots based on the ML principle and taking into account mobility as well as channel dispersion. In [12] and [13], a frequency-domain approach is considered for various pilot schemes inspired from WiMAX but using the auxiliary pilots [14] or pair of pilots [15] (POP) principle. Due to this frequency-domain approach, it leads to lower complexity algorithms but the performance suffers for large values of the CFO and/or STO and it is more suited to a tracking scenario or for refining the estimation.

In this paper, we are interested in low-complexity synchronization methods using closed-form expressions while still providing accurate estimations. Because it is easier to implement in many system architectures, we focus on a frequency-domain implementation, i.e., working with the demodulated symbols after the receiver's analysis filter bank. As opposed to the literature described above, we do not focus on a particularly short preamble [8,10] or scattered pilots [11,13], but we instead consider a specific preamble designed to alleviate the interference structure of the OQAM modulation without requiring the use of auxiliary pilots or POP. This preamble is relatively long and might not be appropriate for low-latency applications but is able to provide efficient synchronization. In particular, the length helps improving the CFO accuracy. For this preamble, we design a specific STO estimation and also show that accurate CFO can be obtained with a simple adaptation of a known technique. As with other frequency-domain methods, the best performance is obtained for offsets (both CFO and STO) which are not too large. So, it might be necessary to perform a very low complexity coarse estimation before

applying the filter bank, in order to ensure that the STO is within reasonable range. Especially large STOs degrade the performance of this estimation technique. This will be illustrated in the simulation results.

In this paper, we will focus on the OQAM flavor of FBMC. The OQAM modulation sends symbols on the real and imaginary part alternatively with T/2 spacing. Because of this structure, it is frequent to perform fractionally spaced equalization at the receiver, using T/2 spacing at the output of the analysis filter bank [16]. The STO estimation method proposed here will be using this double sampling rate at the receiver. Other flavors of FBMC, such as FMT, do not necessarily have this double sampling rate. The method can be generalized to those cases as well, but it requires that the double sampling rate be introduced at the receiver, at least for the duration of the preamble.

The rest of the paper is organized as follows. In Section 2, the FBMC/OQAM system is described. In Section 3, the preamble is introduced, and we explain how the STO and CFO can be estimated using this received preamble. The simulation results of the CFO and STO estimation will be presented in Section 4, and the performance of the proposed method is compared with the LS approach of [7].

2 FBMC/OQAM system model

Consider an FBMC/OQAM system with M subcarriers, as shown in Figure 1. At the input of the transmitter, QAM symbols are converted to OQAM, which is represented by the C2R block on the figure. The QAM symbols have a duration T, with 1/T being the subcarrier spacing. The sampling rate is M/T at the output of the transmitter. For the description of the FBMC/OQAM, we use a formalism based on real symbols similar to the one used in [1]. The purely real OQAM symbol for subcarrier k at sampling instant nM/2 will be denoted by dR[nM/2]. The alternatively real and imaginary symbols to be transmitted are denoted as (see Figure 1)

M = dR M M

n— n — ok n—

2 k 2 2

M n — 2

k+n mod 2

The prototype filter is denoted by a[m]. The output of the transmitter s[m] can be written as

M-1 to s[m] = Y^ Y dk

k=0 n= to

n — a m — n —

¿m km (3)

The prototype filter used in this paper is a root raised cosine filter as defined in [17]. In the z-domain, this filter will be called A(z), with polyphase filters Ak(z) (as shown

Transmitter

Figure 1 The system model.

Receiver

in Figure 1). This filter is used in a DFT-modulated filter bank. The length of the prototype filter is KM, with K being the overlapping factor and M the number of subcarriers. In the frequency domain, neighboring subchannels will overlap, and this causes interference. When we send an impulse on one subcarrier, contributions will be received on the neighboring subcarriers. Other subcarriers only have negligible contributions. By using OQAM, this interference can easily be removed.

The received symbols after the analysis receive filter bank are denoted by xk [nM/2] for subcarrier k at sampling instant nM/2. After multiplication with 6^[«M/2] (the complex conjugate of 6k[nM/2]), the real part is taken to recover the estimation of the initial real symbol. The obtained value is denoted by xR[nM/2]. Now, for timing estimation purposes, it is worthwhile to look at xk[m] in between the symbol instants or, in other words, neglecting the downsampling that occurs in the receiver's synthesis filter bank. For an ideal channel, it can be written as

xk [m] = Y^ dk

' M' n — 2

m — n — 2

œ + Y dk+1 M ' n — 2 a k+i,k m M ' — n — 2

n=—œ

œ + dk—1 " M ' n — 2 a k—1,k m M ' — n — 2

ak,k' [m] = (a[m] eM km^ * ^a[m] e

'— k'r

" M km

is the convolution of the prototype filters on subcarriers k and k! (we use * to denote a convolution). The second and third terms of (4) are the interference terms from neighboring subcarriers. Non-neighboring subcarriers have

negligible interference, thanks to the spectral containment of the prototype, i.e., ak,k+w[m] ^ ak+w,k[m] ^ 0 for integer w > 1.

When taking into account the influence of CFO 0, STO S, and channel impulse response c[m], the received signal s[m] can be written as

S[m] = (s[m + S] * c[m])ej2n0(m+S)/M + n[m], (6)

where n[m] is the additive noise.

3 Joint CFO and STO estimation

The preamble suggested in this paper has a duration of four multicarrier symbols, i.e., 4T. The nth preamble symbol on the kth subcarrier will be denoted by pk [nM/2] in the transmitter and the corresponding received samples by yk [m] in the receiver (similarly to dk [m] and xk [m] for data symbols). The preamble can now be defined as

if n e {0,4} and k is even 0 otherwise

The power of one nonzero symbol is G. The sign of a nonzero symbol can be chosen arbitrarily to improve the peak-to-average power ratio (PAPR) of the preamble but should be the same for symbols on the same subcarrier. On odd subcarriers, the preamble only has zeros. This is to avoid the interference on even subcarriers which can not easily be mitigated before estimation of the channel, the STO and the CFO. On even subcarriers, the preamble has exactly two nonzero symbols spaced 2T from each other. This relatively large spacing, while still reasonable, allows high-precision CFO estimation and also alleviates the OQAM interference issues, making it possible to estimate the STO via the early-late tracking technique presented below. The tail of the preamble only consists of zeros to avoid the interference coming from subsequent data symbols.

n=—oo

n=—oo

The received preamble is processed right after the IFFT in the receiver, i.e., the subchannel processing blocks in Figure 1. According to (4), and assuming that the channel is approximately flat inside each subcarrier, the received symbols on subcarrier k (denoted by yk [m] instead of xk [m] when they correspond to the preamble) can be written as

yk [m] = Yh pk

Ck a k,k

+ vk [m]

= ± vG Ck (ak,k [m] +akM [m — 2M]) + vk [m],

where vk [m] is the additive noise sample, and where Ck is the channel coefficient on subcarrier k. The channel is assumed to be constant on the duration of the four preamble symbols. We assume additive white Gaussian noise (AWGN) with variance an2. In case of CFO $ and STO 5, this becomes

yk [m] = 4GCk (ak,k,$,s [m] +ak,k,$,5 [m — 2M]) + vk [m]

ak,k',$,5[m] = (a[m + 5] e"M-(k+$)(m+5)^ * (a[m] e'

'— k>m M km

3.1 STO estimation

The STO estimator is based on the observation of the amplitude of the received preamble symbols |yk [nM/2]| on all subcarriers k for the first part of the preamble n = 0,1,... ,4 (the second part n = 5,6,7 is potentially corrupted by intersymbol interference from the data symbols that follow). Note that even though the preamble is nonzero only for n = 0 and n = 4, all samples contain some information for the purpose of timing estimation, and we can thus take advantage of the structure of OQAM working at T/2 to utilize the overall information here.

In order to understand the derivation of the STO estimator below, it is interesting to investigate the amplitude of the received preamble |yk[m]| on the different subcarriers k for all sample instants m. As an example, the amplitude |yo[m]| for subcarrier k = 0 is illustrated in Figure 2 for an ideal channel in the absence of noise. Note the raised cosine filter shape caused by the root raised cosine prototype filter in the filter bank. The STO can be estimated by looking at the difference in amplitude between the received preamble symbols |yk[M/2]| and |yk[3M/2]|, similarly to the way it is done for early-late tracking, and as it is illustrated in Figure 2. For instance, when the STO increases, the amplitude of yk [M/2] will decrease while the amplitude of yk [3M/2] will increase. To cope with frequency selective channels and to increase

the precision, |yk [M/2]| and |yk [3M/2]| are combined for all even subcarriers k.

The estimation method proposed here is using four amplitude samples per subcarrier: |yk [0]|, |yk [M/2]|, |yk[3M/2]|, and |yk[2M]|. It is based on the early-late principle [18] and can be derived by using a few approximations and assumptions:

• The four amplitude samples are modeled as linearly dependent on the STO, using a first-order approximation around 5 = 0. In particular, the samples |yk [0] | and |yk [2M] |, which have a zero slope around 5 = 0 (see Figure 2), are assumed to be roughly independent of the STO. This approximation is obviously valid only for small STO and makes the method less accurate at high STO. This effect can be partly compensated by using the overall reference functionas defined and explained below, which provides a reasonable range to the method.

• The noise variance is assumed to be constant on all subcarriers (before applying any equalization coefficient). This is usually a valid assumption.

• The combination across all subcarriers is performed using maximum ratio combining (MRC), which requires knowledge of the channel coefficients amplitudes. To this end and based on the approximation described above, the samples |yk [0] | and |yk [2M] | are used as estimations of the channel amplitudes.

• The channel coefficients are assumed to be constant on the duration of the preamble, which is the case for most applications.

The expression of the estimator is derived below. Based on the linear approximation described above, the amplitude sample |yk [0] | can be written as

|yk [0]| = VG|Ck ||a k,k,$,5 [0] +aw,5 [-2M]|

+ nk,0 (11)

^ Vgc| + nk,0 (12)

since ak,k,0,0[—2M] = 0 and ak,k,0,0[0] = 1 due to the normalization of the prototype, and where nk,i = nk [iM/2] denotes the contribution of additive noise on the amplitude samples of interest®. Similarly,

|yk rn^^VGck |+nk,4.

For the middle points, performing a linear approximation around 5 = 0, we get

|yk[m/2]| ^ VG|Ck| (|ak,k,$,5=0[M/2]

+ a k,k,$,5=0[—3M/2]|— Sk,$5) + nk,1,

Figure 2 Amplitude of a received preamble for k = 0 with G = 1. The time axis has a sampling rate of M/T (with M = 512 in this case), while the processed preamble symbolsy0[nM/2] (crosses in the figure) have a rate of 2/T. Ifan STO 5 is applied, the amplitude ofy0[nM/2] will be different (circles in the figure).

where is the slope of the amplitude with respect to the STO

Sk$ = -

Similarly,

\Ук [3M/2] |

3 \акш [M/2] +ак,к,ф,& [-3M/2]

VG|Ck\ (|йк,к,ф,8=о[3М/2]

+ аш,8=о[-М/2]\+ S^8)

+ Пк 33.

and optimizing the weights to minimize the estimation variance under the constraint of an unbiased estimator, it can be shown that the MRC weights corresponding to the different subcarriers must be proportional to |Ck |. Hence, the overall MRC estimate can be written as

1 M—1

5 = -- J2 C|(|yk[3M/2]| — |yk[M/2]|) (19)

Due to the symmetry of the prototype, it is easy to show that the slopes at M/2 and 3M/2 are exactly opposite to each other and that the linearization points at M/2 and 3M/2 have the same amplitude:

akM,5 [m/2] + ak,k,$,5 [—3M/2]| = |ak,k,$,5=0[3M/2] +

ak,k,$,5=0 [—M/2] |. Hence, based on the linearization and on the early-late principle, a first quantity proportional to the STO can easily be obtained from the samples at subcarrier k:

5 k = |yk [3M/2]| — |yk [M/2] | (17)

= 25Sk$ Ck|VG + (nk,3 — nk,1). (18)

Now, one such quantity can be obtained for each subcarrier k. All theses quantities can then be combined using MRC to form an estimate of the STO. It can be shown that for an ideal channel and for the prototype filter used here, the slopes Sk$ are identical for all subcarriers k. Based on this, assuming identical noise variances on all subcarriers

with some normalization coefficient Anorm. In practice, the channel amplitudes are not yet available, so the values \ук [0]\ and \ук [2M]\ are used as estimates of the channel amplitude inside each subcarrier. The estimation is then normalized in order to be independent of the channel coefficient. Finally, only even subcarriers are taken into account as no symbols are sent on odd subcarriers in the chosen preamble. In the end, the estimation is based on the following quantity:

z(8, ф) =yt -y;

J2 \У2к' [M/2] \\у2к'

к'=0_

J2 \У2к' [о]\2 к/=0

J2 \y2k'[3M/2]||y2k'[2M]| n = -. (22)

/T M/2-1 V '

J2 \y2k' [2M]\2 k=0

Note that this quantity is a function of both the STO S and the CFO $ as emphasized in the notation. It is represented in Figure 3 as a function of the STO when there is no CFO ($ = 0), for a prototype filter with overlapping factor K = 4 and for an ideal channel in the absence of noise. It appears clearly that it is approximately linear on a significant range of STO values and can therefore be used efficiently to perform the STO estimation. In theory, the function can even be used if it is not linear, as long as it is a known one-to-one relationship with the true STO. In this paper, both methods are considered. We start with the more general one, assuming a known one-to-one relationship between the STO and the value of the quantity (20). In order to analyze this relationship, we define the so-called reference function. This reference function will be denoted by z(S, $) and is defined as the value of z(S, $) for an ideal channel and in the absence of noise (the effect of noise will be investigated in more detail in Section 3.1.3). In other words, z(S, $) represents the actual measured value computed with (20) to (22), while z(S, $) represents the theoretical value that would be obtained on an ideal channel and in the absence of noise. If a reasonable estimate $ of the CFO has been obtained (for instance

using the technique explained in the next subsection), the STO can be estimated as

S = argmin \z(A, $) - z(S, $))\ (23)

In the second part, we consider a linear approximation of the reference function which provides a simpler but less precise estimation.

3.1.1 General version

Let us first analyze the reference function z(S, $). As previously stated, this reference function is defined as z(S, $) on an ideal channel and in the absence of noise. Figure 3 illustrates this reference function z(S, $) for three values of the CFO. It is unbiased and exhibits a very good linearity except for large STO (close to ±M/2). The slope of the curve however depends on the CFO. This is further illustrated in Figure 4 which represents z(1, $) as a function of the CFO $. A larger slope is of course preferable as it makes the estimate less sensitive to additional noise. So, the estimation method performs better when the CFO is small although the difference is not very large, as can be seen in Figure 3.

The principle of the estimation, as described in (23) is to compute a reference function in advance and identify which value of the STO corresponds to the observed value of the quantity (20). Note that z(S, $) does not have to be recalculated for each estimation. It can be precal-culated and stored in memory. Therefore, in a practical implementation, the minimization of (23) does not require a long search over a large set of values; it simply corresponds to a look-up table. The estimation method is thus

x 10-3

Figure 4 The reference function z(1, 0) in function of the CFO 0 in the case of an STO & = 1. The larger the CFO in absolute value, the smaller z(1,0).

of low complexity; it amounts to the computation of one closed-form expression (20) followed by a look-up table. Regarding the memory needed, the STO is discrete, but the CFO is not. The reference function should be precal-culated for a number of CFOs and interpolated for the others. The larger that number, the more precise the reference function (and hence the STO estimation) will be and the larger the memory usage as well. As can be seen in Figure 4, z(8, 0) = z(S, -0) which can help reduce the memory usage.

3.1.2 Linear approximation

In order to reduce the memory usage even more, the reference function can be approximated linearly:

~z(S, 0) = z(0,0) + z(1,0)8 (24)

It is clear from Figure 3 that this approximation is quite accurate for moderate values of the STO. For large STOs, the approximation error becomes more significant however. Using this approximation, the complexity of the STO estimation reduces even further:

. y\ - n - z (0,0)

8 =-, (25)

z (1,0)

3.1.3 Effect of the noise

When AWGN is added to the channel, all the amplitude samples \yk[iM/2]\ are corrupted by noise. Now, since the noise on the initial yk[iM/2] samples is Gaussian,

the probability density function of the amplitude samples \yk [iM/2] \ is a Rice distribution. In particular, it also means that the average effect of the noise is not zero. On average, the respective contributions of the noise on y^ and y\ do not cancel each other, and the estimate z(8,0) deviates from the reference function z(8,0). The overall effect is illustrated in Figure 5 which represents the average value of the estimate z(8,0) in the presence of noise as a function of the STO 8 and when the CFO 0 = 0. Two SNR cases (15 and 25 dB, respectively) are presented, and the result is compared to the reference function z(8,0) in the absence of noise. Once again, the effect is negligible for small STOs and more significant at high STOs. This generates an estimation error that gets larger for higher STOs. However, it is interesting to observe that the average effect of the noise at high (positive) STO is to decrease the estimate z(8, 0), which is the opposite of the nonlinear behavior of the reference function z(8, 0) that tends to deviate above the linear slope. The overall result is that the average estimate z(8, 0) exhibits an even better linear behavior than the reference function z(8,0) as can be seen on Figure 5. In order to explain this, a complete analytical derivation of the noise distribution for z(8,0) would be long and tedious, so we restrict ourselves to a qualitative justification which is provided in the Appendix.

3.1.4 Effect of the multipath channel

The frequency selectivity of the channel also has an influence on z(8, 0), not only for large STO but for the entire range. For instance, the bias z(0,0) might not be zero

Figure 5 Effect of the noise on the reference function. The mean of z(S,0) on a 25- and 15-dB AWGN channel, as well as the reference function z(S,0).

anymore depending on the channel impulse response taps. The longer the channel, the larger the divergence with the reference function z(S, <) can be.

To improve the estimation, it is possible to use some basic information about the channel. The idea is to assume some statistical channel model and try to take its effect into account in the reference function. A new reference function zmult(5, <) is used in that case, that is simply replacing z(S, <) which was calculated for an ideal channel. This new reference function zmult(<5, <) is defined as the expectation of z(S, <) in (20) in the absence of noise and averaged over the possible realizations of the channels, according to the chosen model. In practice, it is difficult to obtain the true expectation; so, the practical computation of zmult(<5, <) comes down to computing it for a certain number of realizations and compute the average.

Just as previously, this new reference function is computed in advance without the knowledge of the true channel realization, but some channel model needs to be available. Obviously, the accuracy of the model has a direct impact on the performance of this method. Several results are presented below in the simulation section.

3.1.5 Complexity

Even though a detailed complexity analysis would depend on the chosen implementation, and hence is outside the scope of this paper, a few comments can be made on the issue of complexity. As mentioned above, the proposed method relies on a closed-form expression, and does not require a min or max search over a potentially

large number of candidates, which helps reduce the complexity significantly. The method also assumes that the frequency-domain samples of the preamble are available, so the method is for instance very well suited to an architecture where the analysis filter bank is implemented separately and applies to all received symbols, including the preamble.

3.2 CFO estimation

The CFO estimation used here is a direct application of the one presented in [19] for OFDM. Similar CFO estimation methods have also been used for FBMC/OQAM systems in [12,13,20] although for different preamble schemes. The CFO < is estimated by looking at the phase difference between the received preamble symbols yk[0] and yk [2M] on each even subcarrier k. The estimated CFO will be denoted by <f>:

1 /M/2-1 \

4nz I T> y*2k' [0] y2k' [2M] j (26)

With the preamble considered in this paper, the distance of 2T between yk [0] and yk [2M] is quite large. This improves the precision of the estimation but also limits the range of CFOs that can be estimated correctly. More precisely, this only allows correct estimation of CFOs in the range of < e [-0.25,0.25]. A CFO of < = 0.30 would be estimated as <p = -0.20. Because of the noise, the practical range of this estimator is of course much

smaller than [-0.25,0.25] and depends on the SNR of the channel.

To cope with this problem, a heuristic adjustment has been used. It is taking into account the sign of the phase difference betweenyk[M/2] andyk[3M/2]. This phase difference will be denoted by 0s. When the CFO is large and there is a risk of ambiguity, 0s is taken into account. When the CFO is small on the other hand, 0s is neglected since it is more susceptible to noise in this case than 0. Hence, the estimated CFO 0 is

0.5 + $ if 101 > 0.15 and sign(0) = -1 and sign (<f>s) = 1 - 0.5 if \0p\ > 0.15 and sign(0) = 1 and sign (fa) = -1 0 otherwise

/M/2-1

& = Z J2 y*2k> [M/2] y2k' [3M/2]

0s = ^ I

\ k'=0

The threshold for using <j>s is set on \0\ = 0.15. This value was chosen to assure correct CFO estimation for CFOs in the range of 0 6 [—0.25,0.25] even when the SNR is low. It is the result of a trade-off but does not come from any specific theoretical justification.

Note that in the method proposed here, the CFO is estimated before the STO. Hence, the CFO estimation is sensitive to the actual STO (as it could not be compensated yet). This is mainly due to the interference between the preamble symbols. For 8 = 0, there is no interference from one preamble symbol to the other on yk [0] and yk [2M]. However, when the STO increases, the interference increases which modifies the observed phases and degrades the CFO estimation.

4 Simulation results

To assess the performance of the CFO and STO estimation, the technique presented in this paper is compared with the LS technique of [7]. Note that this LS technique is a time-domain algorithm and hence corresponds to a different implementation architecture, with different constraints on the complexity. It is however also a preamble-based technique, and it is an appropriate benchmark to evaluate the performance of the proposed method. The preamble used in [7] has a duration of 3T, so we have added a zero guard symbol of length T to reach the same total length 4T as the preamble used in this paper, without

any impact on the estimation technique. Both preambles have been normalized for equal transmitted power. About 104 trials were performed under the following conditions:

• The number of subcarriers was M = 512. The overlapping factor of the prototype filter was K = 4 ([17]).

• The (normalized) CFO was uniformly distributed in 0 6 [-0.25,0.25]. Notice that this is the maximum range that can be estimated correctly. In order to have a good estimation for CFOs on the edges of this range, the threshold to use the adjustment was set on \0\ = 0.15, as in (27).

• The STO was simulated in two ranges:

- STO 8 6 [-M/2,M/2] = [-256,256]

- STO 8 6 [-M/16,M/16] = [-32,32]

• The multipath channel has been modeled to consist of 17 independent Rayleigh fading taps h(l) with an exponentially decaying power delay profile. Specifically, E[\h(l)\2] = Ce-l/4, where the constant C is chosen for total unit energy ^\=0 E[\h(l) \2] = 1. The channel was different in each trial.

• The STO was estimated (with a granularity of one sample) using reference functions z(8,0), zmult(8,0), and their respective linear approximations. All of these reference functions were sampled in the CFO domain with a step size of 0 = 0.01.

When zmult(8,0) is used, it is computed based on the 17-tap channel model detailed above. Hence, in this case, the model used for computing zmult(8,0) is the same as the model used to generate the channels, but the true channel realization is of course not known and may be different from the ones used in computing the function zmult(8,0). As will be shown below, it proves that some basic knowledge on the channel (delay spread and power delay profile) can already help improve the method.

The first set of simulations were done for an STO uniformly distributed in the range [-M/16, M/16] = [-32,32]. The results are shown in Figure 6. On an AWGN channel, the STO was estimated with the reference function z(8,0) as well as its linear approximation. The linear approximation performs equally well since the linearity is quite good in the range 8 6 [-32,32]. On a multipath channel, the STO was estimated with the reference functions zmult(8,0) and their respective linear approximations. Again, the linear approximations are doing equally well. Note that the root-mean-square error (RMSE) of the STO estimation using zmult (8,0) is a lot lower than the RMSE using z(8,0). All techniques exhibit an error floor at high SNR in the presence of multipath

-B- AWGN - z

AWGN - z

approx

—'— Multipath z

— * - Multipath - z approx

—&— Multipath - z

- > - Multipath - z mult,appro)

-0- AWGN - LS

■ —1— Multipath - LS

0 5 10 15 20 25 30

SNR [dB]

Figure 6 RMSE of the STO S, with S s [-32,32] on both an AWGN channel and a multipath channel. On the AWGN channel, the STO was

estimated with both the reference function z(S, and its linear approximation. On the multipath channel, the STO was estimated with the reference functions z(S, and zmujt(5, and their respective linear approximations.

due to the impact of the channel on the reference function which is not perfectly known. The LS technique from [7] provides a much higher RMSE than the early-late techniques in this STO range, except for one particular case: at high SNR, in the case of multipath channel and if the reference function is used without taking into account the multipath model.

The second set of simulations was performed for an STO uniformly distributed in the range [-M/2, M/2] = [-256,256]. The results are shown in Figure 7. The LS technique has the same performance as in the previous set of simulations. Since it estimates in the time domain, its performance is independent of the actual STO range. The performance of the early-late technique on the other hand has degraded heavily. There is a clear RMSE difference now between the early-late technique with full reference functions and the early-late technique with linear approximated reference functions. This is because the linear approximation is less accurate for large STOs. The results also show that for SNRs lower than about 10 dB, the RMSE is lower when using the linearly approximated reference functions. This is due to the effect of the noise, as explained in the Appendix, that tends to compensate for the nonlinearity of the reference function and provide an overall better linearity. It even appears that this linearized version itself performs better for more practical values of the SNR, with a minimum around SNR = 10 dB. Note that when the full reference functions are used

(not their approximations), the early-late technique still performs better than the LS technique.

In Figure 8, the results of the CFO estimation are shown. Again, the RMSE of the LS technique is independent of the actual STO range. The frequency-domain technique presented in this paper, on the other hand, is highly sensitive to the STO range (remember that the CFO estimation is performed before the STO estimation here). For an STO 5 e [-256,256], the RMSE is higher than the RMSE of the LS technique. It might therefore be useful to have at least a coarse estimate of the STO before performing this CFO estimation. It is possible for instance to reduce the uncertainty on the STO to roughly [-128,128] by comparing the amplitudes of yk[-M/2], yk[0], and yk[M/2] beforehand. Note also the high RMSE when the SNR = 0 dB. It is caused by CFOs at the edges of the range <p e [-0.25,0.25] being estimated as CFOs at the opposite edges, causing a very large estimation error. Although the use of <ps corrects some of these errors, it is obviously not perfect, especially at low SNR.

5 Conclusion

The simulations have illustrated that the presented CFO and STO estimation technique outperforms current timedomain estimation techniques for small STOs. The low complexity of the technique makes it even more attractive. However, since the estimation is done in the frequency domain, the estimation error will increase when the actual

AWGN - z

AWGN - z

Multipath - z

' - * - Multipath - z r approx

. —&— Multipath - z

- f> - Multipath - z „ mult,appro)

-0- AWGN - LS

—h- Multipath - LS

0 5 10 15 20 25 30

SNR [dB]

Figure 7 RMSE of the STO &, with & s [-256,256] on both an AWGN channel and a multipath channel. On the AWGN channel, the STO was

estimated with both the reference function z(8,0) and its linear approximation. On the multipath channel, the STO was estimated with the reference functionsz(8,0) and zmujt(8,0) and their respective linear approximations.

STO and CFO increase. Hence, it is advisable to have a prior coarse estimation. This is not the case for timedomain estimation techniques. Having the possibility to reuse the preamble for channel estimation purposes is another advantage.

The focus of this paper was on FBMC/OQAM, since the double sampling rate of the analysis filter bank was used to estimate the STO. FMT does not have this double sampling rate. As stated before, by introducing this to FMT (at least for the duration of the preamble),

Figure 8 Results of the CFO estimation. RMSE of the CFO 0 for 8 6 [-32,32] and 8 6[-256,256] on both an AWGN channel and a multipath channel.

everything presented in this paper can be applied on FMT as well.

Endnote

aThis is the contribution of the noise taking into account the norm operation. It is no longer a Gaussian noise. More detail is provided in Section 3.1.3.

Appendix

Effect of the noise on the STO estimation

This section is aimed at explaining more precisely why the effect of the noise on average at high STO is to decrease the estimate z(5, and can thus help improve the linearity of the reference function. We define \yk,d [m]\ as the ideal value of \yk [m] \ in the absence of noise. As mentioned above, in the presence of noise, \yk[m]\ is a Rice distribution and its expected value is always larger than without noise: E[\yk[m]\] > \yk,id[m]\. In addition, it is easy to show based on the properties of the Rice distribution that the difference E[\yk[m]\] - \yk,d[m]\ is larger when the ideal amplitude \yk,d [m] \ is small and vice versa. Now, we are interested in the expected values of (21) and (22) in the presence of noise. The numerator and denominator are not independent, but for this qualitative inspection, we approximate the expectation as the ratio of expectations. The effect of the noise on the denominators is a fixed value equal to the sum of the noise variances. Now, it is easily seen that for very small STO, \y[0]\ and \y[2M]\ have similar distributions, as well as \y[M/2]\ and \y[3M/2]\. Hence, the effect of the noise is equal on average for y^ and y\, and it has a zero mean on the estimate z(S, $). For high (positive STO), \yk,d[M/2]\ becomes very small and on the contrary \yk,id [3M/2]\ is larger. Based on the above comments about Rice distribution, the average of the noise will be larger on y^ than on y\ (\yk[0]\ and \yk[2M]\ still have similar distributions). Hence the noise will have the tendency to decrease the estimate z(S, $). This is exactly what we observed on Figure 5.

In conclusion, if the general version of the algorithm is used, the noise generates an estimation error that gets larger with higher STO. If the linear version is used, however, the noise can be useful to improve the linearity of the overall curve. This effect was confirmed in the simulation results. In both cases in addition, the denominators in \y2k [0] \ and \y2k' [2M] \ can get small when the STO is large. This means that the overall estimation gets more sensitive to the noise (it increases the noise variance). Clearly, the larger \5\, the larger the STO estimation error. This was also confirmed in the simulations of Section 4 for 5 e[-M/2, M/2].

Competing interests

The authors declare that they have no competing interests.

Acknowledgements

This work was partially supported by the European project EMPhAtiC

(ICT-318362).

Author details

1 Broadcom Corp., Battelsesteenweg 455B, Mechelen B-2800, Belgium. 2IMEC,

Kapeldreef 75, Leuven B-3001, Belgium. 3Universite catholique de Louvain,

Place du Levant, 2, Louvain-la-Neuve B-1348, Belgium.

Received: 1 December 2013 Accepted: 26 June 2014

Published: 25 July 2014

References

1. MTanda,T Fusco, M Renfors, J Louveaux, M Bellanger, Deliverable

2.1—data-aided synchronization and initialization (single antenna).Tech. Rep., FP7-ICT PHYDYAS - PHYsical layer for DYnamic AccesS and cognitive radio (2010)

2. H Bolcskei, Blind estimation of symbol timing and carrier frequency offset in pulse shaping OFDM systems. Proc IEEE Int. Conf. Acoust. Speech Signal Process. 5,2749-2752 (1999)

3. P Ciblat, E Serpedin, A fine blind frequency offset estimator for OFDM/OQAM systems. IEEE Trans. Signal Process. 52(1), 291-296 (2004)

4. G Lin, L Lundheim, N Holte, New methods for blind fine estimation of carrier frequency offset in OFDM/OQAM systems, in Proceedings of the 7th IEEE Workshop on Signal Processing Advances in Wireless Communications (SPAWC) (Cannes, France, 2-5 July 2006), pp. 1-5

5. T Fusco, A Petrella, MTanda, Non-data-aided carrier-frequency offset estimation for pulse-shaping OFDM/OQAM systems. Signal Process. 88(8), 1958-1970 (2008)

6. V Lottici, M Luise, C Saccomando, F Spalla, Non-data-aided timing recovery for filter-bank multicarrier wireless communications. IEEE Trans. Signal Process. 54(11), 4365-4375 (2006)

7. T Fusco, A Petrella, MTanda, Data-aided symbol timing and CFO synchronization for filter bank multicarrier systems. IEEE Trans. Wireless Commun. 8, 2705-2715 (2009)

8. T Fusco, A Petrella, MTanda, Joint symbol timing and CFO estimation for OFDM/OQAM systems in multipath channels. EURASIP J. Adv. Signal Process. 2010, 897607 (2010)

9. ATonello, F Rossi, Synchronization and channel estimation for filtered multitone modulation, in Proceedings of the International Symposium on Wireless Personal Multimedia Communications (Abano Terme Italy, 2004), pp. 590-594

10. D Mattera, M Tanda, Data-aided synchronization for OFDM/OQAM systems. Signal Process. 92(9), 2284-2292 (2012)

11. V Lottici, R Reggiannini, M Carta, Pilot-aided carrier frequency estimation for filter-bank multicarrier wireless communications on doubly-selective channels. IEEE Trans. Signal Process. 58(5), 2783-2794 (2010)

12. TH Stitz, A Viholainen, T Ihalainen, M Renfors, CFO estimation and correction in a WiMAX-like FBMC system, in Proceedings of the 10th IEEE Workshop on Signal Processing Advances in Wireless Communications (SPAWC) (Perugia, Italy, 21-24 June 2009), pp. 633-637

13. TH Stitz, T Ihalai, A Viholainen, M Renfors, Pilot-based synchronization and equalization in filter bank multicarrier communications. EURASIP J. Adv. Signal Process. 2010, 741429 (2010)

14. JP Javaudin, D Lacroix, A Rouxel, Pilot-aided channel estimation for OFDM/OQAM, in Proceedings of the 57th IEEE Semiannual Vehicular Technology Conference (VTC), vol. 3 (Jeju, South Korea, 22-25 Apr 2003), pp. 1581-1585

15. J-P C Lele, R Javaudin, A Legouable, P Skrzypczak, Siohan, Channel estimation methods for preamble-based OFDM/OQAM modulations. Eur. Trans. Telecomm. 19(7), 741-750 (2008)

16. D Waldhauser, L Baltar, J Nossek, MMSE subcarrier equalization for filter bank based multicarrier systems, in IEEE 9th Workshop on Signal Processing Advances in Wireless Communications (SPAWC) (Receife, Brazil, 6-9 July 2008), pp. 525-529

17. MG Bellanger, Specification and design of a prototype filter for filter bank based multicarrier transmission, in Proceedings of the IEEE International Conference on Acoustics, Speech, and Signal Processing (Salt Lake City, UT, USA, 7-11 May 2001), pp. 2417-2420

18. M H Meyr, S Moeneclaey, Fechtel, Digital Communication Receivers: Synchronization, Channel Estimation, and Signal Processing. (Wiley, New York, 1998)

19. PH Moose, A technique for orthogonal frequency division multiplexing frequency offset correction. IEEE Trans. Commun. 42,2908-2914 (1994)

20. H Saeedi-Sourck, Yan Wu, JWM Bergmans, S Sadri, B Farhang-Boroujeny, Low-complexity carrier frequency offset estimation for multiuser offset QAM filter bank multicarrier systems uplink, in Proceedings of the 75th IEEE VehicularTechnology Conference (VTC) (Yokohama, Japan, 6-9 May 2012). 5 pages

/ \ doi:10.1186/1687-6180-2014-118

Cite this article as: Van Caekenberghe etal.: Preamble-based frequency-domain joint CFO and STO estimation for OQAM-based filter bank multicarrier. EURASIP Journal on Advances in Signal Processing 2014 2014:118. --1_I_/

Submit your manuscript to a SpringerOpen journal and benefit from:

7 Convenient online submission 7 Rigorous peer review 7 Immediate publication on acceptance 7 Open access: articles freely available online 7 High visibility within the field 7 Retaining the copyright to your article

Submit your next manuscript at 7 springeropen.com