Mattera etal. EURASIP Journal on Advances in Signal Processing (2015) 2015:23 DOI 10.1186/s13634-015-0191-4

0 EURASIP Journal on

Advances in Signal Processing

a SpringerOpen Journal

RESEARCH Open Access

Analysis of an FBMC/OQAM scheme for asynchronous access in wireless communications

Davide Mattera1, Mario Tanda1* and Maurice Bellanger2

Abstract

The OFDM/OQAM transceiver belongs to the filter-bank-based multicarrier (FBMC) family and, unlike OFDM schemes, it is particularly able to meet the requirements of the physical layer of cognitive radio networks such as high level of adjacent channel leakage ratio and asynchronous communications. The paper proposes and analyzes a new implementation structure, named frequency spreading, for the OFDM/OQAM transceiver. On flat channels, it is equivalent to the standard one in terms of input-output relations, though more complex. On multipath channels, it offers a crucial advantage in terms of equalization, which is performed in the frequency domain, leading to high performance and no additional delay. With its flexibility and level of performance, the analyzed scheme has the potential to outperform OFDM in the asynchronous access context and in cognitive radio networks.

Keywords: OFDM/OQAM; Multicarrier systems; Prototype filter; FBMC; Asynchronous access; Cognitive radio

1 Introduction

The cognitive radio transmission context exhibits a number of specific features which make it significantly different from the conventional transmission environment [1,2]. First, the available bandwidth is likely to be fragmented, i.e., it is made of non-adjacent spectrum chunks that have to be exploited jointly for high speed data communications. Then, the sections of the spectrum that are not available might be occupied by a primary user and a high level of protection must be provided. Specifically, the transmission system must guarantee a high level of adjacent channel leakage ratio (ACLR)a. Next, the transmission band is likely to be changing on short notice or even without notice. On the exploitation side, the total bandwidth available might be dedicated to a single user requiring high bit-rates or it can be dynamically shared by several users in proportion to their instantaneous capacity needs. If opportunistic operation is contemplated, these users have the freedom to show up and disappear as they wish. In such conditions, a rigid communication procedure, where each user must be aligned before the transmission can start, is inadequate. In fact, asynchronous operation is necessary to reach an acceptable

Correspondence: tanda@unina.it

1 Dipartimento di Ingegneria Elettrica e delleTecnologie dell'Informazione, Universita degli Studi di Napoli Federico II, via Claudio 21, 80125 Napoli, Italy Full list of author information is available at the end of the article

level of spectral efficiency. Clearly, to cope with such a context, an appropriate physical layer is required.

The spectrum granularity offered by multicarrier transmission techniques has proven its efficiency for spectrum exploitation, and the most popular technique, orthogonal frequency division multiplexing (OFDM), has been widely used in communications for more than a decade now. However, for the cognitive radio context as described above, it lacks flexibility and it is likely to lead to poor spectral efficiency, even with the introduction of additional processing [3]. Therefore, an enhanced multicarrier technique is needed, as pointed out in [4], where it is shown that a filter-bank-based multicarrier (FBMC) physical layer can meet the ACLR requirements [2,4]. In particular, FBMC/OQAM can overcome the limits of OFDM provided that we impose a constraint on the cut-off frequency of the prototype filter, which cannot exceed the sub-channel spacing as pointed out in [5]. High stop-band attenuation filters have been proposed that do not satisfy this constraint, such as isotropic orthogonal transform algorithm (IOTA) [6] and Hermite filters. These filters are associated with single-tap equalizers as mentioned in [7]: such a scenario does not allow the exploitation of the potential advantages ofFBMC/OQAM systems for cognitive and, therefore, they cannot compete with CP-OFDM. On the contrary, the high performance equalization objective can be met if the prototype filter (employed in the

Springer

© 2015 Mattera 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/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited.

FBMC/OQAM scheme) is designed using the frequency sampling technique introduced in [8] and developed in [5]. With this simple approach, the coefficients are derived from a few samples of the filter frequency response, which makes the implementation of the filter bank in the frequency domain practical.

When comparison between such FBMC/OQAM system and OFDM system is considered, it appears that a key advantage of OFDM with cyclic prefix (CP) is the capability to achieve perfect channel equalization, as long as the channel impulse response remains shorter than the guard time provided by the CP. Thus, in order to be accepted, the FBMC/OQAM approach must have a high performance equalization capability, particularly in the asynchronous context, characterized by the fact that the system must compensate simultaneously the timing offset, the frequency offset, and the channel distortion.

In the absence of CP, the equalization capability of the FBMC system rests on the sub-channel equalizers, which cannot be single-tap as in OFDM, but must be multitap to reach similar performance. However, the use of multitap equalizers implies an increase of the receiver latency; this motivates the search for an equalizer structure that does not introduce such a disadvantage. The main contribution of the paper lies in the proposal of a new transceiver structure for FBMC/OQAM systems that is able to provide satisfactory performance without increasing the transceiver delay and accepting the presence of significant timing and frequency offsets among the users that are performing the multiple access, as it is common in a cognitive radio scenario.

The concept has been presented in [9,10], along with preliminary performance results, under the name frequency spreading (FS)-FBMC, but a rigorous analysis of the corresponding scheme is still missing. An objective of the present paper is to provide such an analysis and prove the equivalence of FS-FBMC with the conventional polyphase network (PPN)-FFT scheme in both transmitter and receiver. This equivalence is important because it opens the way to mixed implementations. For example, in uplink transmission, the distant user can be equipped with the conventional IFFT-PPN transmitter, while the high performance but more complex FS-FBMC receiver is implemented at the base station.

Many advances in the applications of FBMC to various scenarios will be able to take advantages from the proposed transceiver structure. In particular, the capability to use multiple antennas at the transmitter and/or at the receiver, which significantly increases bandwidth efficiency, can be easily carried out along the lines introduced in [11-13], which however do not take into account the frequency-despreading structure at each receiver; further works is needed to define the details of the MIMO

extension of the proposed structure. Alternative structures are also under consideration for achieving the same goal of operating on multipath channel with a minimum implementation complexity. For example, the fast-convolution structure [14], which is currently under study for its extension on the multipath channel, is superior to the proposed structure in terms of computational complexity. It is equivalent in terms of flexibility (e.g., it shows a similar capability to easily compensate a time-offset in the frequency-domain as suggested in [15]), while the fast-convolution transceiver latency is larger than that achieved by the frequency-despreading system [14].

The organization of the paper is as follows. In Section 2, the FS-FBMC scheme for the transmitter is described and the proof of the equivalence with the standard FBMC, namely the IFFT-PPN cascade, is provided. Section 3 is dedicated to the receiver structures and, again, the proof of the equivalence between FS-FBMC and standard FBMC, namely the cascade of PPN and FFT, is provided; moreover, in Section 3, it is shown that the FS-FBMC structure is computationally more complex while in Section 4, it is shown that on a multipath channel, it offers a crucial advantage in terms of equalization, which is performed in the frequency domain, just like OFDM, leading to high performance and no additional delay. In Section 5, the performance of the proposed scheme is illustrated and contrasted with the results obtained for OFDM and standard FBMC when the sub-channel equalizer has a single-tap. In Section 6, the main aspects of FS-FBMC are summarized and the potential impact is discussed.

Notation: j = sf—1, superscripts (•)*, (-)T, and (-)H denote the complex conjugation, the transpose, and the conjugate transpose, respectively, is the real part, log is the base-2 logarithm, <g> is the linear convolution, S [//] is the Kronecker delta, ceil[x] is the smallest integer larger than or equal to x, and (•) denotes the time average, i.e.,

<x[n]) = 2N+T EN=-Nx[«] and modM№ =

1 - qM with q such that modM(W) e {0,1,...,M - 1}. Moreover, we denote with DFT[x] the vector x whose //th component can be written as Tx/ = N ^¿io1 xie-j~Nkl and with IDFT[x] the vector x whose //th component can be written as x/ = "£.4=0 xiejNkl where xi is the lth component of the N x 1 input vector x. Finally, lowercase boldface letters denote column vectors, x the componentwise product between two vectors and, finally, 0 denotes the null vector.

2 The transmitter with standard and frequency-spreading structures

Let us consider an FBMC system using offset QAM modulation, often designated by OFDM-OQAM [16]. We

assume that the number M of subcarriers be even; the low-pass transmitted signal s(t) can be written as

s(t) = sR(t) + js1 (t - T/2)

Nh+Ns-1

sR(t) = j Tt+ f)g(t - nT) (2)

n=0 keA

Nb+Ns-1

s (t) = £ Tt+ % (t - nT)

n=0 keA

s[i] = sR (iTs) + js1 ((i - M/2)Ts).

g[i] = g(iTs).

The generation of the sequence sR [i] is equivalent to the generation of the sequence of M x 1 vectors d<nR) whose kth component d^l is equal to sR[nM + k] for k e {0,1,..., M — 1}. In the following, we consider two implementation structures and their implementation complexities: though the standard implementation structure based on an IFFT over M points exhibits a reduced computational complexity, the frequency-spreading structure based on an IFFT over a larger number of points provides useful insights into the structure of the transmitted signal.

2.1 Standard transmitter structure

The standard structure, often named the polyphase network, for the implementation of the OFDM/OQAM transceiver has been first proposed in [17,18]. To make clear its comparison with the proposed alternative, we briefly

recall its derivation here. The kth component d«'k =

sR[nM + k] of d"R) can be written as

Nb+Ns—1

dnkk = E bn-kg[l + « — "'M k e {0,1,..., M — 1}

where T is the multicarrier symbol interval, A С {0,1,..., M — 1} is the set of active subcarriers whose size is Mu, the sequences aRk and aInk indicate the real and imaginary parts of the complex data symbols transmitted on the kth subcarrier during the nth QAM symbol, Nb is the number of training symbols, Ns is the number of payload symbols, while g(t) is the prototype filter. It is assumed that the data symbols aRk and aInk are statistically independent with zero-mean and variance a The discrete-time low-pass version s[j] = s{t)\t=iTs of

the transmitted signal (Ts = T/M is the sampling interval) can be written as

bR = n,k

jk<«) jMk

k e {0,1,..., M - 1}

which is the IDFT of the sequence )kttRk with respect to

the index k. If we define the vector as the M x 1 vector whose kth component (for k e {0,1,..., M — 1}) is b^l in Equation 8, we can compactly write

bnR) = IDFT |w x aiR1] (9)

where IDFT[-] denotes the IDFT operator on the input vector and, for k e A, the kth component wk of the M x 1-vector w is

In the next subsections, we report the derivation of an efficient generation procedure for the signal sR(t). An analogous derivation can be straightforwardly obtained for the signal s1 (t). Since the continuous time signal is generated by D/A conversion, we consider the generation of its discrete-time samples

Nb+Ns—1

sR[j] = sR(iTs) = J2 J2alkeik{ M'+ % [i — nM] n=0 keA

where we have used Equation 2 and the following definition

Wk = j

and the kth component of the vector a(nR) is the symbol aRk in Equation 2 while, for k e A, wk = 0 and the components of aiR) are irrelevant. Note that Equation 9 is only defined for n e {0,1,..., Nb + Ns — 1}, but we can straightforwardly extend it to any n provided that we assume that

anR) = bnR) = 0

n e {0,1,...,Nb + Ns - 1}.

We can compactly write Equation 7 as

Nb+Ns-1

dnR) = £ b

X gn — n

where the vector gn is defined so that its kth component

gn,k is

gn,k = g[k + nM]

ke {0,1,..., M - 1}.

Therefore, Equation 13 defines the polyphase components of g[•]. The prototype filter g(t) satisfies the following property

g (t) = 0 t e [0, KT)

where K is the overlap parameter. The property in Equation 14 implies that the vector gn is nonnull only for n e {0,1,..., K — 1}. Consequently, Equation 12 can be rewritten as

More specifically, let us first define the vectors h® such that

= [hm,iM hm,iM+1 hm,iM+2 ■ ■ ■ hm,iM+M-1] ^

i e {0,1, ■■■,K - 1}.

= go X bf-' + gi X b^-i + ■■■ + gK-1 X h{n-_{K-1). From Equations 22 and 23, it follows that

Analogously, the generation of the sequence s(I)[i] in Equation 4 is equivalent to the generation of the sequence of vectors d^"1 defined as the output of the PPN:

dtf = go X bni) + g1 x +...+ gK-1 x b^^ (16)

h(0)T h(1)T h(K-i)T ^'it) ■■■ 11™

b« = IDFT

w x a,

Using Equations 23 and 24, Equation 19 can be re-written as

d(nK) = hf + h- + hi-2 + --- + hi/-^K)-1). (25)

We have now to study the structure of the vector hn in Equation 24 for a general time-step n in order to simplify its generation. From Equation 20, it follows that its kth component hn/ can be written as

where the /th component (// e X) of the vector a^ is the K,k = g /] E cRejM ^ k e{0,1,..., KM - 1} symbol aIn k in Equation 3. Again, we have assumed

o — bm — 0

PI — — v

*n — "n

n e {0,1, ■■■, Nb + Ns - 1}.

2.2 Frequency-spreading structure

In the present subsection, we derive an alternative structure for implementing the OFDM/OQAM transmitter, named the frequency-spreading structure. By using Equations 7 and 8, in Appendix A it is shown that the //th component d^/ of d(n" can be written as:

d^l = hnM + hn-H+M + hn—2,/+2M + hn—3,/+3M

k'eA K-1

k'=-(K-1) K-1

k'=-(K-1)

i 2n. kk e> km kk

i 2*. kk1 e> km kk

lk'e.4

(R). e'KM Kkk

IKM-1 \

km (26) m=0 J

where we have assumed that the prototype filter has been designed according to a frequency-sampling approach so that its Fourier transform

G(F) = g[i] e-i2nFi

+ ■■■ + hn-(K-1),k+(K-1)M

satisfies the property

where we have defined

hm,p = g [p] E C ep e{0,1,..., KM -1}.

c(R) a J fafyp n' e{0,1,..., Nb + Ns - 1} and / e X Cn',//' [ 0, n' e {0,1,...,Nb + Ns - 1} and V/

Equation 19, which has to be applied for / e{0,1,..., M—1}, is similar to that present in the classic fast-convolution procedure, usually named overlap-and-add structure; here, it operates on the sequence of KM x 1 vectors hn, defined as follows:

G(F)| k = 0 k e {K,K + 1, ■■■,KM - K}.

The Fourier coefficients in Equation 26 are given by

Vk. (29)

Gk =-G(F)i k

k KM ( If= KM

Note that the 2K — 1 nonnull values {G/}K=-(/K—1) are the free parameters of the prototype filter when it is designed according to the frequency-sampling procedure used in [5].

Furthermore, in Equation 26, we have also introduced the following definition of the K-times upsampled version of c(R in Equation 21:

hm = [hm,0 hm,1 ■ ■ ■ hmKM-1]T

c(KR) A C m = kK nm \ 0 otherwise .

From Equation 26, it follows that

hn = IDFT[G] xIDFT[cn] = IDFT[zn]

where the vectors G and cn in Equation 31 are defined as follows:

G = [Go Gi G2 ... Gkm-\]t

= [G0 G1 . . . GK-1 01 x [KM — 2K +1] GKM-(K-i) . . . GKM—1] T

cn —

C^O °1 x(K—1) cR O1 x(K—1) C(R O1 x(K—1)

c(nM—1 01 x (K—1)

and the vector zn in Equation 31 for the well-known property of the DFT and IDFT operators is the circular convolution of the two KM x 1 vectors G and cn. We have used the property in Equation 28 to write Equation 33. In Appendix B, it is shown that the components of the vector zn in Equation 31 can be written as

zn,pK+k = j

aR,pGk + /'aR,modM(p+1)Gk—K (1 — ¿[k] )

wherep e {0,1,..., M} and k e {0,1,..., K — 1}.

Therefore, the frequency-spreading structure requires:

1. to use the input symbols aRk to calculate (for

t e {0,1,..., KM — 1}) the components zn,t of zn according to Equation 35. Note that the symbol aR,p is spread over 2K — 1 components of the vector zn and for this reason, the structure is named FS-FBMC; in fact, each component of zn is dependent on two adjacent symbols and each symbol anp, according to Equation 35, not only determines the component pK of the frequency-domain vector zn but also spreads its effect, weighted by the frequency response of the prototype filter, on the different components of the same vector ranging from pK — (K — 1) up to pK + K — 1;

2. to determine hn starting from zn by performing the IDFT over KMpoints in the right-hand side of Equation 31;

3. to evaluate d^ by the overlap-and-add processing defined in Equations 24 and 25.

2.3 Complexity comparison of the two structures

The standard transmitter structure requires to calculate (a) the IFFT over M samples according to the definition in Equation 9 then to calculate the vector d.^ according to the PPN (Equation 15). The frequency-spreading structure requirements have been just summarized.

In a structure with a single processor, the complexity comparison is equivalent to the count of the number of flops required by the two structures. The number of complex flops for calculating the IFFT over M samples can be

written as 1.5M log(M) while the number of real flops can be writtenb as 4M log(M) — 6M + 8 as in the split-radix.

The number of complex multiplications for calculating IFFT over M samples can be written as 0.5M x(log M — 1) while the number of real multiplications can be written as Mlog(M) — 3M + 4 [19] by removing most of the trivial operations and using three real multiplications per complex multiplication.

To calculate the vector dn ) according to the PPN (Equation 15), the following number of real-valued flops are necessary: 2M real-valued multiplications for each of the K terms and M complex-valued additions for each of the K — 1 couples of vectors to be summed.

With the frequency-spreading structure, we need 2KM — M real-valued multiplications for calculating zn according to Equation 35 while the IFFT for calculating hn by using Equation 31 requires 4KMlog (KM — 6KM + 8 real-valued flops. Finally, (K — 1)M complex-valued additions are needed by the overlap-and-sum structure in Equation 25.

Therefore, for generating each vector dn}, with the standard structure, we need C(ST) real-valued flops or C(mT) real-valued multiplications where

C(ST) — 4M log(M) — 6M + 8 + 2MK + 2M (K — 1)

(9) (15)

— 8 + M [4 log(M) + 4K — 8] C((ST) — M log(M) — 3M + 4 + 2MK — M [log(M) + 2K — 3] + 4

while with the frequency-spreading structure, we need c(FT) real-valued flops or c(T> real-valued multiplications where

C(FT) — 2KM — M + 4KM log(KM) — 6KM + 8 + 2(K — 1)M

(35) (31) (25)

— 8 + M

4K [log(KM) — 0.5] — 3

C(T) — 2KM — M + KM log(KM) — 3KM + 4

K [log(K) + log(M) — 1] — 1

Note that C(FT> ~ MK

4log(M + 4log (K ) — 2

4log(M + 4K — 8

— KC

The first approximation is obvious while the second approximation holds provided that M is sufficiently large and K sufficiently small (e.g., for K = 4 and M = 1024, we

have a normalized approximation error of 6.1%). Moreover, Cf)/C(FT) is around 58% for K = 2, 35% for K = 4, and 24% for K = 8, mainly independently of M e {512,1024,2048,4096}. Therefore, the frequency-spreading structure is about K times more complex when implemented using a structure with a single processor.

3 The receiver for standard and frequency-spreading structures

According to Equations 4 and 5, assuming perfect synchronization on the flat channel (i.e., the channel with unit response in the frequency domain) and neglecting the presence of the noise, the received signal r[i] can be written as follows

Nb+Ns-1

r[i] = a ^RjM'+f)g [i - nM]

n=0 /eA

+ jaIn,k(-1)//ei//(M'+ % [i - nM - M/2] J = aaRnkAMi+ n)g[i - nM] +o(n) [i] (39)

jaan,/

i M i+ n

Mi+ i)g [i - nM - M/2] + o(I) [i]

where 0 < a < 1 denotes the channel gain and o(n) [i] denotes terms present in the transmitted signal that do not depend on ann,/ while o(I)[i] denotes terms present in the transmitted signal that do not depend on aInk, with n e {0,1,..., Nb + Ns} and / e X. We consider at the receiver the decision variables:

n,/ = m DR

an,/( = m ) Dn,/

D(R = 1 J2 r[i] e-/§)g[i - nM]

D(Ii = —j

J2 r[i] —fe^iMi+ f)

x g[i - nM - M/2]

where A is a proper constant amplitude, defined so that R = aR, when o(n) [i] = 0 and a^ = aIn/ when o(I) [i] =

M- / - M- /

n,kn,k

0. Consequently,

A = a£.

with Eg = Y^,K==0—1 g2[i]. Equation 41 is motivated by a proper design of the prototype filter g[-] that guarantees a negligible projection of the interference terms o(n) [i] on the matched filter, as usual in FBMC receiver [16,20,21]:

In fact, according to Equation 39, o(n)[i] denotes the additive signals present in the received signal r[i] that do not depend on the useful symbol ann,k; the condition in Equation 44 therefore implies that such additive signals do not interfere with the useful signal when the matched-filter projection (designed according to the useful term) is performed (i.e., the result of the matched-filter projection is independent of the interfering signals); on the other hand, the matched-filter projection is optimum (in the maximum-likelihood sense) when the interference signals are not present and only the noisy version of the useful term is taken into account. Consequently, the decision variable in Equation 41 operating on the flat channel without synchronization error is optimum (in the maximum-likelihood sense) for estimating statistically independent information symbols; analogously, the same optimality holds for the decision variable in Equation 42. Obviously, the receiver implemented according to Equation 41 has to be modified in order to operate on a multipath channel. However, before discussing such modifications, we first need to describe the two structures implementing Equation 41. Consequently, we recall the standard structure for implementing Equation 41 in Subsection 3.1 and we introduce an alternative structure in Subsection 3.2; moreover, we compare their complexities in Subsection 3.3.

3.1 The standard receiver structure

In the present subsection, we briefly recall the derivation of the standard receiver structure. From Equations 41 and 14, it follows that

dR = A j\ E r[nM + i] e-j/ M ig[i]

_ i=0 M-1

+ r [(n + 1)M + i] e-jkMig[M + i] i=0 M-1

+ ^ r [(n + 2)M + i] e-jkMig [2M + i] + ... i=0

M-1 n 1

+ ^ r [(n + K - 1)M + i] e-jkMig [(K - 1)M + i] 1. i=0 j

Let us introduce the vector rn whose ith component rn,i (i e {0,1,..., M - 1}) is defined as follows

rn,i = r[nM + i]. (46)

By using also Equation 13, let us consider the vector

o(n)[i] ej M i+ f)g [i - nM] - 0.

Jg,n) a

= rn+W x ,

e {0,1,...,K - 1}

and let us denote with r+i its ith component for i e {0,1,..., M— 1}. Using such definitions, from Equation 45, it follows that

AjkD(R = £ r^fe-** q + £ r^e'*q

q=0 M-1

, V r(gn) e-jkM q+ + V rgn) e-jkM q

+ / , rn+2,qe M + ••• + Z^ rn+(K-1),qe

n(K) n(K) n(K) Dn,0 Dn,1 • • • Dn*-1

A DR X w = M DFT

+ M DFT

r(g,n) rn+1

+ M DFT

(g,n) n+(K-1)

= MJ2 DFT [tn+^ x gi]

or, equivalently,

DnR) = M— x DFT nA

Analogously, in Appendix C, it is shown that

tj\ w*

D(I = M— x DFT n jA

r,„ ) o X ;

D(I) d(D

Dn,0 Dn,1

D(R = J E r[nM + ¿] e-jkM' j £ Gk'eJ

r I i=0

k'=-(K-1)

E G*' E r[nM + i] e 'km

kmm i(kK+k!)

k'=-(K-1)

where we have taken into account the properties of the prototype filter (see Equations 14, 28, and 29) and the fact that it is real and, consequently, G—k = G*k.

Let us introduce the KM x 1 vector r(nK) whose kth component r^k is defined as

rK = r[nM + k]

k e {0,1,..., KM - 1}

Defining the vector D^ =

taining the decision variables on the vector a^K it follows from Equation 48 that

and the vector Rn, whose kth component is denoted with Rn,k, defined as the DFT over KM samples of the vector r(K).

Rn = DFT

Using such definitions, Equation 54 can be re-written as

DR = Aj-kKM J2 G*' Rn,,

modKM(kK+k0

k e {0,1,M - 1}.

k'=-(K -1)

Let us define the (2K - 1) x 1 vector R^f (k e{0,1,...,

M -1}) ,

a(FS) A r

Rn,k = LRn,modKM(kK-(K-1)) . . . Rn,modKM(K ■ ■ ■

n,modKM (kK + (K

-1))]1

and the vector

where the vector Dn) contains the decision variables on the vector a?1:

G(K) = [G— (k—1)... G—1 G0 G1... GK—1]r. (59)

Taking into account Equations 58 and 59, Equation 57 becomes

G(K )H r(FS) G Rn,k

while the ith component r^i (i e {0,1,..., M — 1}) of the vector r^ is defined as follows

n — r [nM + M/2 + i]. (53)

3.2 The frequency-despreading receiver structure

In the present subsection, we derive the receiver counterpart of the frequency-spreading transmitter described in Subsection 2.2. From Equation 41, it follows that for k e {0,1,...,M — 1}

Analogously, in Appendix D, it is shown that the decision variable on a(p, can be written as

D(I) _ KM -k-1

g(k r r:

where we define the (2K — 1) x 1 vector R^k^ (k e {0,1,..., M — 1}) as follows ,

R(I,FS) A

Rn,k =

n,modKM (kK-(K-1)) . . . n,modKM(kK)

n,modKM(kK+(K-1))

The vector R

R(I) = DFT r(IK

whose /th component is denoted with n( / and used in Equation 62 is defined as the DFT over KM samples of the vector tn^ whose /th component r^/P is defined as

r^f) = r [nM + M/2 + /] / e{0,1,..., KM - 1}.

n,/ (64)

Therefore, Equations 60 and 61 define the following receiver structure, illustrated in Figure 1:

1. Collect the samples of the received signal r [•] to form the KM x 1 vector rin Equation 55 (and half

period later the vector r^/^ in Equation 64);

2. Calculate the vector Rn in Equation 56 by FFT over KM points (and half period later the vector R^*1 in Equation 63);

3. Perform M different projections of the vector Rn according to Equations 60 and 58 in order to obtain each decision variable DR for / e {0,1,..., M — 1} (and half period later according to Equations 61 and 62 to obtain Dn/). Thus, to obtain the datum a^/ that in the transmitter has been spread over 2K — 1 components of the vector zn (see the first point of the sentence after Equation 35), the same components of the vector Rn are exploited by using as weights the (conjugate) Fourier coefficients with / e {-(K - 1),..., K - 1}. For this reason,

such a structure is called frequency-despreading receiver: in fact, it collects all the components of the frequency-domain vector Rn dependent on the useful symbol a/k, due to the spreading performed at the transmitter, and weights them according to the frequency response of the prototype filter, achieving the despreading of the useful symbol.

The importance of the proposed structure is not limited by the assumption in Equation 28 because, when it is necessary to introduce a possible mismatch (i.e., to use at the transmitter a prototype filter that does not satisfy Equation 28), it can be managed with very marginal performance loss (i.e., the frequency despreading receiver can be still employed at the receiver, with its advantages considered in the paper and without appreciable disadvantages due to the presence of a mismatch).

3.3 Complexity comparison of the two structures

The standard receiver structure (see Equation 50) requires to calculate (a) MK multiplications, (b) K — 1 additions of M x 1 vectors, and (c) an FFT over M samples, which is the same requirement of the standard transmitter structure. Therefore, the overall number of real-valued flops Cjn and the overall number of real-valued multiplications

C® are equal to their transmitter counterpart CjT and

C® in Equation 36. The frequency-despreading structure

Figure 1 Frequency despreading. The frequency despreading structure at the receiver side. On a flat channel Hm = 1 for any m while on a multipath channel, the values of Hm depends on the channel frequency response and are used to equalize the channel distortions.

is described in Equations 60, 58, and 56. With reference to such a structure, the overall number of real-valued flops CjR and the overall number of real-valued multiplications Cm can be written as

CjK = 4KM log (KM) - 6KM + 8

FFT in Equation 56

2(2K - 1)

+ 2(2K - 2)

multiplications in Equation 60 additions in Equation 60

ke{0,1....M-1} in Equation 60

= 8 + M [4K log (KM) + 2K - 6]

CFm = KMlog(KM) - 3KM + 4 + 2(2K-1) •

FFT in Equation 56 multiplications in Equation 60

ke(0,1....,M—1) in Equation 60

= 4 + M [K log(KM) + K - 2]. Note that

CfR = MK [4 log(M) + 4 log(K) + 2] + 8 - 6M

~ MK [4 log(M) + 4K - 8] ~ KCf (66)

where the first and the second approximations are obvious while the third one holds provided that M is sufficiently large and K is sufficiently small (e.g, for M = 1024 and K = 4, 4 log(M) + 4log(K) + 2 = 50 while 4log(M) + 4K - 8 = 48). Therefore, under the same assumptions used at the transmitter side, we can obtain the following approximation: CjR = KCjR; moreover, C0/C0 is about 50% for K = 2, 30% for K = 4, and 21% for K = 8, independently of M e (512,1024,2048,4096). Thus, also for the receiver case, the frequency-spreading structure has a computational complexity about K times larger.

The complexities of the different considered structures are summarized in Table 1. Note that the analysis shows that the complexity increase is due to the fact that the frequency-spreading structure requires to evaluate the DFT over KM points instead of the DFT over M points required by the standard structure. Since many solutions

exist for implementing FFT processing, the result of the comparison depends on the particular solution for its implementation; however, a K-fold increase of the implementation costs can be considered an upper-bound; not always, however, the installed processing power can be optimized to the effective needs and, consequently, the increase of the actual costs (needed to perform the FFT over a larger number of points) can be much smaller depending on the effective implementation details.

4 Adapting the frequency-despreading structure to the multipath channel

In the present section, we first define the adaptation of the frequency-despreading structure to the multipath channel, then we recall a standard approach to adapt the standard structure to the multipath channel and we finally compare their performance.

4.1 Frequency-despreading structure operating on multipath channel

When the multipath nature of the channel is taken into account, the received signal r[i] can be written as

r[i] = J2 h[t] s[i - I] +n[i]

where the complex-valued sequence h[i] oflength Lh + 1 models the multipath channel.

In Equation 67, the noise term n[i] is assumed to be the ith sample of the low-pass equivalent of a white Gaussian random process with null average and power spectral density equal to N0/2 in the signal bandwidth; the analog low-pass equivalent has null average and power spectral density equal to 2N0. The low-pass equivalent is sampled with sampling period Ts; assuming an ideal antialiasing filter with bandwidth equal to 1/(2Ts), the random variable n[i] ideally sampled at the filter output has null average and variance VAR[n[i] ] = 2N0/Ts. Moreover, the sequence n [ i] is zero-mean, white, Gaussian and circularly symmetric; thus,

n[n] n*[n - m]

n[n] n[n - m]

Table 1 Structure complexity

Structure Number of flops Number of multiplications

Standard transmitter 8 + M [4 log (M) + 4K - 8] M [log(M) + 2K - 3] + 4

Standard receiver 8 + M [4 log (M) + 4K - 8] M [log(M) + 2K - 3] + 4

FS transmitter 8 + M 4K [log(KM) - 0.5] - 3 M K [log(K) + log (M) - 1] - 1 + 4

FS receiver 8 + M [4K log (KM) + 2K - 6] 4 + M [K log(KM) + K - 2]

The number of flops and the number of multiplications of the two transceiver structures for each half multicarrier symbol period.

On a multipath channel, an equalization stage has to be included in the structures already described in Section 3. In particular, when the frequency-despreading receiver (see Equation 57) is considered, the following modified structure is considered

— 1rkKMY, F(k) Rn,modKM(kK+k0 k e{0,1,..., M—1}

k'——(K —1)

where the complex-valued coefficients F^ in Equation 69 replace the coefficients Gy in Equation 54; therefore, the coefficients F^ can be set to Gk> obtaining the structure for a flat channel; on a multipath channel, we can set them in order to equalize the channel improving the receiver performance. The value of a to be used in Equation 43

is given by a = Pr/Ps where Pr = (e \h[i] ®s[i] represents the average power of the useful component of

the received signal r[i] and Ps = [E [|s[ i] |2] represents the average power of the transmitted signal s[i].

The decision variable D(Rk in Equation 69 can be written as follows:

d(R) _ j(R,k)an j .

Dn,k — j0,0 an,k + jn,k + Vn,k

m?—1)m « 4)*}

— ^{(—1)( (v2(—1,qfk

fk —

'(k) (k) (k) Ak) Ak)

F—(K—1) . . . F—1 F0 F1 ... FK—1

influence the useful coefficient and the interference power; we can, therefore, set the vector fk in order to equalize the effects of the multipath channel in Equation 67.

The noise term nn,k in Equation 70 is a zero-mean complex-valued Gaussian random variable with variance

2||G(K°|a2 where y is defined as the signal-to-noise ratio per

subcarrier, i.e., y — N where Es is defined as the energy of the useful term of the received signal in a multicarrier symbol period that is dedicated to each active subcarrier.

On a flat channel (h[n] — S[n]), we denote with v^^'^ the vector v^' in Equation 71; its k'th component v^qk^'^ (k e {(—K — 1),..., —1,0,1,..., K — 1}) canbe written as

v(FLAT,k) n,q,k'

KM+nM/2—1

l|G(K )|

11 r—q(—1)k«j Kn ^ g[i] e-

i—nM/2

Then, the properties of the prototype filter and the choice fk — G(K) for setting the receiver coefficients guarantee that the useful coefficient /0R,k'1 — 1 and the interference coefficients are practically null; in practice, the achieved signal-to-interference-ratio is around 65 dB, i.e.,

i,(FLAT,k)" (K )

G(K ^ 1

H v™" G(K^|« 1

V(m, q) — (0,0) (74)

where l00k) represents the coefficient of the useful term aRk, nn,k describes the effect of the background noise, and In,k describes the intersymbol and intercarrier interferences of the symbols aR—m k—q ((m, q) = (0,0)) and a1 , on the useful symbol aR,. Such interferences

n— m,k—q n,k

would be negligible on a flat channel but they become significant on the multipath channel (see Equation 67).

It is easy to prove that /n,k can be written as a linear combination of the interfering symbols; the coefficients of such a combination can be written as

The satisfaction of the condition in Equation 74 is equivalent to the condition in Equation 44 and it concerns the design of the prototype filter. Moreover, it guarantees the optimality (in the maximum-likelihood sense) of the receiver structure on the flat channel.

When all the components of the vector VmnLq'1^ are nonnull, the vector Vmna can be written as v(k) _ v(FLAT,k) w(k)

m,q = m,q x m,q

where the k'th component w()qk' of the vector w(;q is defined as

m,q,k'

KM+mM/2—l—\ qK+k' .

E g [i] e—' 'KM'1

i—mM/2—l

KM+mM/2—\ k+U .

E g[i] e—' 'km

i—mM/2

—j2n

kK+k M KM l

where iC'* and are the coefficients of the symbols of an and aj , respectively. This shows that the

n—m,k—q n—m,k—q L 1

coefficients F® of the receiver structure, collected in the (2K — 1) x 1 vector

When v^f'^ = 0, then the denominator of Equation 76 is null and Equation 75 is still valid provided that we replace with unit both the denominator of Equation 76 and the same quantity in Equation 73.

From Equations 71, 74, and 75, it follows that the choice

fk — G(KVw!

(k)* m,q

where we denote with / the component-wise division of the two vectors and (m, q) = (0,0), impliesc that

HI v(mq fk

However, a single vector fk has to be chosen and, consequently, Equation 78 can be satisfied for a single value of (m, q). A robust choice usually done for setting fk lies in

using Equation 77 with (m, q) = (0,0), which implies that the useful coefficient is practically set to unit:

f^ ^ 1

Such a derivation shows that also the optimum choice of fk may not able to achieve the optimum performance achieved on a flat channel, unless the vector wOq exhibits a weak dependence on (o, q) so that a single vector fk can approximately satisfy the condition in Equation 77 for any value of (m, q) or at least for the most important values of (o, q). Such a weak dependence is present when Lh ^ M; in such a case, the quantity in Equation 76 can be approximated as follows:

KM+mM/2-e-1 qK+k,

E g [i] ~KMT' i=mM/2-l

KM+mM/2-1 qK+k,.

E g[i] e-j

i=mM/2

Vl < Lh

and, consequently, the vector wOq can be approximated by the following vector h, which is independent of (o, q):

hk = [Hki<-(k-1) ... HkK-1 HkK HkK+1 • • • HkK+k-1]T.

where Ho denotes the multipath channel frequency response at the frequency F = o/(KM). Consequently, the choice in Equation 77 with the use of the expression hk for wOq, i.e.,

fk = G(K)/h* ^ F(kk) =

of the transceiver. It may appear that the FS equalizer be equivalent to a multitap sub-channel equalizer following the standard structure and therefore that the considered comparison be unfair. However, they are not equivalent for two reasons: (a) because the delay introduced by the two structures is different and obtaining the minimum delay is important in a transceiver, like the OFDM/OQAM one, with an already larger delay in comparison with the OFDM system; (b) the PPN-FFT scheme performs equalization after sampling rate reduction which introduces an interpolation operation. The distinction of the two structures in terms of sampling rate reduction lies in the fact that the FS structure performs equalization in its internal behavior and therefore before sampling rate reduction while the PPN structure performs equalization after sampling rate reduction and consequently needs to use the single-tap equalizer to not increase the transceiver delay.

4.3 Comparing the signal-to-interference-and-noise ratios of the two structures

The equivalence of the two structures described in Section 3 on a flat channel implies that the behavior of the single-tap equalizer can be described with the same relations introduced in Subsection 4.1 provided that the coefficients expression in Equation 82 is replaced by the following one

p(k) _ Fk, =

guarantees to the structure the same performance achieved on flat channel provided that Lh ^ M (see Figure 1 for a scheme of the FS structure in Equation 69 when Equation 82 is chosen).

4.2 Recalling the standard approach

The standard approach consists in including an equalizer stage in cascade with the structure described in Subsection 3.1. The effects of the multipath channel can be equalized by using a single-tap structure [7]. In this case, the kth entry of the DFT output (see Equations 50 or 51) is multiplied by 1/H(Fk) with Fk = M in the standard structure for k e (0,1,..., M - 1).

More sophisticated multitap structures could be used and have also been proposed with reference to the standard structure [22]. Since they operate in the frequency domain and at twice the multicarrier symbol rate, they introduce an additional delay proportional to the number of taps. We consider the single-tap equalizer in both structures since it maintains limited the overall latency

In other terms, the choice in Equation 83 makes the FS receiver equivalent to the PPN structure equipped with the single-tap equalizer of coefficient 1/Hk. Therefore, by comparing Equations 82 and 83, we can obviously note that the advantage of the frequency-despreading equalizer lies in its capability of using the coefficient Hk+k' instead of the constant term Hk. Since the FS structure in Equation 69 first extracts the DFT of the input vector r(nK) according to Equation 56 and subsequently uses the coefficient F^ to equalize the channel effect at frequency (kK + k')/(KM) = k/M + kk/(KM), the FS structure with the choice in Equation 82 uses the right coefficient (i.e., HkK+k ) to equalize the channel response at frequency (kK + k')/(KM) while the FS structure with the choice in Equation 83, which is equivalent to the PPN structure equipped with the single-tap equalizer, always uses the same coefficient HkK to equalize the channel responses at the different frequencies (kK + k')/(KM) for k' e {-(K - 1),..., -1,0,1,...,K - 1}. In other terms, differently from the PPN structure equipped with the single-tap equalizer, the FS structure is able to equalize with different coefficients the different parts of the subcarrier band. Since the effect of an offset nT in timing synchronization, perfectly compensated however in the frequency domain,

can be obtained by setting h[£] = S[£ — nT], our analysis shows that a performance improvement of the FS equalizer, which for nT = 0 (flat channel) is equivalent to the single-tap equalizer, appears when larger values of nT determine faster variations of the channel frequency response so that non-negligible variations appear within the subcarrier band; in such a case, the FS structure is able to use different coefficients to equalize each part of the subcarrier band and can therefore achieve improved performance in comparison with the PPN structure equipped with single-tap equalizer. Such a superior capability of the FS structure is irrelevant in the presence of a flat channel; therefore, the two structures show the same performance on the flat channel or on channels where the variations on the subcarrier band (of length 1/M) are marginal.

Furthermore, an unified expression of the signal-to-interference-and-noise ratio (SINR) on the //th subcarrier, denoted as SINR/, of the two structures on the //th subcarrier can be written as follows:

SINR/ =

114II2 Y II G*>||2

K=Kmin,K=0

E jv/f) + "E E1 ^ {<f

• "=" L ' "="min q=1 L '

where the SINR of the frequency-despreading receiver is determined by the use of the vector f> described in Equation 82 whereas the single-tap equalizer is determined by the use of the vector f/ in Equation 83. Note

that the signal-to-interference ratio SIR/ can be obtained by employing Equation 84 without the first term at the denominator.

In order to plot the resulting SINRs for the two structures, we have to set the prototype filter g[ •]: we set it to the filter designed in [5]. For the two structures, we report the SIR and the SINR for a choice of the subcarrier index k = 56, M = 1024, and by using 21 subcarriers and 21 time-interval values to approximate the sum in Equation 84. The channels employed to determine the receiver SINR are 1,000 random realizations of the ITU-R Vehicular B [23] and the SINR values, evaluated from Equation 84 and reported in Figures 2, 3, and 4, have been normalized to \HkKI2 in order to separate the effects due to the equalizer capability from that due to the fading of the kth subcarrier.

Figure 5 shows that the frequency-despreading receiver is able to get the same SIR performance that it achieves on the flat channel on a large fraction of the considered different channels, randomly selected and ordered for decreasing values of output SIR achieved by the frequency-despreading receiver. The multipath channels where it cannot achieve the SIR achieved on the flat channel are those with very fast variations of the channel frequency response on the considered subcarrier. The single-tap shows poorer performance mainly because, on many channels realizations, the channel frequency response is not flat in the subcarrier band of size 1/M. From the SIR evaluation, it follows that, when we set the signal-to-noise ratio y between the SIR achieved by the

---Frequency-spreading equalizer

-Single-tap equalizer

100 200 300 400 500 600 700 800 900 1000

Channel Number

Figure 2 Samples of output SINR when y is 30 dB. The normalized output SINRs for the two structures on 1,000 different channel realizations according to the ITU-R Vehicular B model when y is 30 dB.

P4 22 Z

- Frequency-spreading equalizer ■ Single-tap equalizer

100 200 300 400 500 600 700 800 900 1000

Channel number

Figure 3 Samples of output SINR when y is 20 dB. The normalized output SINRs for the two structures on 1,000 different channel realizations according to the ITU-R Vehicular B model when y is 20 dB.

single-tap equalizer and that achieved by the frequency-spreading one, the SIR limitation implies also a SINR limitation and, therefore, a performance advantage of the frequency-spreading structure is present. In fact, from Figures 2, 3, and 4, we can notice a significant advantage of the frequency-despreading equalizer when we set the value of y at 30 dB; such advantage is reduced but it is still significant at y equal to 20 dB while the noise

limitation becomes dominant for y equal to 10 dB and the frequency-despreading structure cannot take advantage from its superior interference-rejection capability.

5 Performance comparison of the two structures

In the present section, we assess via computer simulations the equalization performance achieved by using the frequency-spreading structure and compare it with

Channel Number

Figure 4 Samples of output SINR when y is 10 dB. The normalized output SINRs for the two structures on 1,000 different channel realizations according to the ITU-R Vehicular B model when y is 10 dB.

100 200 300 400 500 600 700 800 900 1000

Channel Number

Figure 5 Samples of output SIR. Samples of the output SIRs for the two structures on 1,000 different channel realizations according to the ITU-R Vehicular B model.

that achieved by the standard structure. We have also included in the performance comparison the classical OFDM system that is often considered for opportunistic transmissions because it is a classical scheme employing the multicarrier approach where many practical difficulties have already been resolved; this has a strong impact on the overall cost. However, for the cognitive radio context, it lacks flexibility and it is likely to lead to poor spectral efficiency. Since the latencyd of the FBMC receiver is K times larger than that of the OFDM receiver with the same number of subcarriers, we have set the number of subcarriers in the OFDM transceiver K times larger in order to compare two structures with the same latency. Moreover, with such a choice, the OFDM receiver and the FS-FBMC receiver perform the FFT procedure on the same size, though FS-FBMC has still to perform it to a rate 2K times larger. In particular, we have used 2,048 subcarriers for OFDM while we have used only 512 subcarriers for FBMC transceiver and we have used K e {2,3,4} in order to verify the effect of the overlap factor.

A number of 104 Monte Carlo trials has been performed under the following conditions:

1. The considered FBMC and OFDM systems have a bandwidth 1 = 11.2 MHz;

2. The transmitted symbols are the real and imaginary parts of 64-QAM symbols;

3. The considered multipath fading channel model is the ITU-R Vehicular B [23];

4. The used prototype filter is that proposed in [5]. Actually, any type of prototype filter can be implemented with an extended FFT, due to the equivalence between time and frequency domains. However, in order to be practical, the number of frequency domain filter coefficients must be the smallest possible, which is the case of the used filter;

5. The channel is fixed in each run but it is independent from one run to another;

6. The residual timing offset (RTO) and the normalized residual carrier frequency offset (RCFO) are controlled as simulation parameters;

7. Both systems exploit a one-tap subcarrier equalizer with perfect knowledge of the channel and of the residual timing error, i.e., when simulating the presence of the timing offset nT, we have used exp(-j2nnTk/M)/H(k/M) instead of 1/H(k/M) as coefficient of the single-tap equalizer in the standard structure and exp(-j2nnTk/KM)/H(k/KM) as coefficient of the single-tap equalizer for the frequency-despreading structure;

8. The effect of the RCFO on the phase of each decision variable in the frequency domain, which increases [20] linearly with time, is not compensated; therefore, the BER is dependent on the specific multicarrier symbol interval considered for equalization. In order to maintain sufficiently limited the effects of such nonideal receiver behavior, we evaluate the BER on the data transmitted in one of the first multicarrier symbol intervals, the eighth one;

9. In order to use the same bandwidth in both FBMC and OFDM, which exhibits a larger spectral leakage, we have set the percentage of active subcarriers in OFDM transceiver as 82% of the overall number of subcarriers while we have set to 89% this percentage in OFDM/OQAM transceiver;

10. The length of the cyclic prefix is 1/8 of the OFDM multicarrier symbol period (note that since in the FBMC system the cyclic prefix is not used, in the considered case an increase of the bit-rate nearly equal to 11.1% with respect to the OFDM system is obtained).

Note that, in consequence of the choices reported at the points 9 and 10, the data rate of the OFDM system is about 82% of the data rate of the FBMC system.

We first consider a simulation scenario where no timing or frequency offset is introduced; Figure 7 reports the results of the experiment. We can verify that the two FBMC structures, which are equivalent in the flat channel, perform differently on multipath channel because of the difference between the equalization mechanisms; more specifically, the frequency-spreading structure provides improved performance with respect to the standard structure and this improvement increases as the overlapping factor increases. Moreover, both structures exhibit a performance floor at large values of Eb/N0. However, while the performance achieved by the standard FBMC receiver is very poor, in the range 0 dB < Eb/N0 < 20 dB,

the performance of the FS-FBMC system is practically equivalent to that of the OFDM transceiver employing 2,048 subcarriers for a sufficiently large value of the overlap parameter K. Since the channel coherence time may impose a smaller number of subcarriers also to the OFDM transceiver, it is interesting to test the OFDM system with only 512 subcarriers; the results not reported for clarity in Figure 7 show a poor performance, practically equivalent to the standard structure, since it is unable to equalize the considered channel (the length of the cyclic prefix (512/8) Ts is smaller than the ITU-R Vehicular B channel length).

We have also performed other simulation experiments to verify the performance on the less hostile ITU-R Vehicular A channel: here, the condition Lh ^ M is better satisfied. In fact, the corresponding results, shown in Figure 6, report that the two structures and the OFDM system are practically equivalent on channel A. Only for larger values of Eb/N0 we can observe some difference; in particular, we note that the single-tap equalizer for K = 2 provides the worst performance; moreover, the three dashed-line curves, corresponding to the single-tap equalizer for K = 3 and K = 4 as well as the frequency spreading structure for K = 2. Only for K = 3 and K = 4 the frequency despreading structure behaves practically equivalent to the OFDM system. Therefore, we can conclude that similar performance is achieved by OFDM and FS-FBMC transceivers also if the latter uses only 512 subcarriers while the former uses 2,048 subcarriers.

—e— 1-tap K=4

— 1-tap K=3

—i— 1-tap K=2

—A— FS K=4

—B— FS K=3

—*— FS K=2

Figure 6 Comparison wrt Eb/N0 (ITU-R Vehicular B channel). The effect of the noise on the performance of the OFDM transceiver (with 2,048 subcarriers) and of the two structures for FBMC transceiver (with 512 subcarriers) on ITU-R Vehicular B channel. Three possible values of the overlap parameters are considered.

-e- l-tap K=4 : - V - l-tap K=3 t— l-tap K=2 ■A— FS K=4 e— FS K=3

Figure 7 Comparison wrt Eb/N0 (ITU-R Vehicular A channel). The effect of the noise on the performance of the OFDM transceiver (with 2,048 subcarriers) and of the two structures for FBMC transceiver (with 512 subcarriers) on ITU-R Vehicular A channel. Three possible values of the overlap parameters are considered.

In the second experiment, where we have imposed the value of the RTO, we have not introduced frequency offset and we have set Eb/N0 equal to 20 dB. The results are reported in Figure 8 and show the superior performance of the proposed FS-FBMC system. The presence of

the cyclic prefix explains the threshold effect on the performance of the OFDM receiver; apart from this effect, we can verify that the OFDM performance worsens and becomes poor for larger values of the magnitude of the timing offset. The performance of the standard FBMC

-200 -100 0 RTO

Figure 8 Comparison wrt RTO. The effect of a timing offset on the performance of the OFDM transceiver (with 2,048 subcarriers) and of the two considered structures for FBMC transceiver (with 512 subcarriers) on ITU-R Vehicular B channel. Three possible values of the overlap parameters are considered.

structure is very poor independently of the value of the overlap parameter K also in the absence of the timing offset and saturates for larger values of the magnitude of the RTO. Regarding the sensitivity of the performance to residual timing offset, it increases as the overlap parameter K decreases.

Furthermore, in the third experiment, we have imposed the value of the RCFO, we have not introduced timing offset, and we have set again Ey/N0 equal to 20 dB. The results reported in Figure 9 show that the performance degradation of the FS-FBMC structure due to the RCFO is more pronounced for higher values of the overlap parameter K; however, such a degradation is lower than that of OFDM system. Moreover, when the CFO is larger than 5/100, the performance of the OFDM and of the FS-FBMC structures is equivalent.

In order to better simulate the working conditions of a cognitive radio scenario, we next consider an uplink scenario where all active subcarriers are block-wise shared among four users whose delays are mutually independent and uniformly distributed within {—M/2, -M/2 + 1,..., M/2 — 1}. Only a single subcarrier is left as guard between adjacent sets of active subcarriers whereas a maximum value of the normalized CFO equal to 0.1 is admitted on each user. The value of K is set to 4 as, in the previous experiments, it resulted to be the best choice for K.

The considered uplink scenario is quite general and can be encountered in a number of cognitive network

architectures, in particular when the cognitive cellular networks are taken into account [24]. For example, consider the case where the four terminals are located in an area where a large number of cells have been deployed by different operators for the local coverage enhancement of a potentially crowded zone interested in opportunistic multimedia downloads at low cost. Each relay is provided of a proper backhaul connection and cooperates with the other relays (and with the terminals) by means of a control channel, designed in dedicated, common, or underlay fashion [25]. The cooperation is aimed at determining according to an optimization procedure (e.g., that proposed in [26]), the subcarriers available for the transmission of each active terminal, and its transmitted power (to be minimized). Such a cooperation obviously does not include a timing alignment procedure not only because it would make much more complex the access control but mainly because it is impossible, i.e., the distances among relays imply that the transmission of each terminal is received by different relays. Therefore, in absence of a timing alignment procedure, the delay and the CFO of the user of interest are perfectly compensated and, then, all the other users remain asynchronous; moreover, we first consider the case where each signal, coming from the other users, at the receiver arrives with the same power of the useful signal.

The results, reported in Figure 10, show the BER versus Eb/N0 for the user of interest. It is important to note that the OFDM system performance is equivalent on the

W 10 PQ

-0.25 -0.2 -0.15 -0.1 -0.05

0.05 0.1 0.15 0.2

Figure 9 Comparison wrt CFO. The effect of a carrier offset on the performance of the OFDM transceiver (with 2,048 subcarriers) and of the two considered structures for FBMC transceiver (with 512 subcarriers) on ITU-R Vehicular B channel. Three possible values of the overlap parameters are considered.

□ Single-tap > FS 0 OFDM

Figure 10 Comparison wrt Eb/N0 (powers equal to that of the user of interest). The effect of the noise on the performance of the OFDM transceiver (with 2,048 subcarriers) and of the two structures for FBMC transceiver (with 512 subcarriers) on both the considered channels for K = 4 in an uplink scenario with the powers of the interfering users equal to that of the user of interest. The continuous line refers to channel ITU-R Vehicular A while the dashed line refers to the ITU-R Vehicular B channel.

two channels and practically equivalent to the result of the FS system on channel ITU-R Vehicular B. The FS receiver performance is similar to that of the single-tap receiver on channel ITU-R Vehicular A and they achieve practically the same results achieved on the ideal flat channel. The single-tap receiver achieves very poor performance on channel ITU-R Vehicular B in comparison with the proposed frequency-despreading structure. Therefore, it is interesting to note that the FBMC system can provide comparable performance with the OFDM system only thanks to the receiver structure here proposed. We have also repeated the experiment when all the signals transmitted by interfering users arrive at the receiver with a power ten times larger than that of the useful signal; we have reported in Figure 11 the corresponding results. We note that the performance of FS structure is practically not affected by the larger powers of the interfering users whereas the OFDM system is severely affected and it exhibits a performance gap with respect to the FBMC system equipped with the frequency despreading receiver structure.

It is interesting to note that, differently from what happens on channel ITU-R Vehicular A, on channel ITU-R Vehicular B the performance of the frequency despread-ing receiver is not equivalent to its performance on the flat channel. This is due to the fact that the powers of the interference terms are increased since the conditions that guarantee the optimality of the frequency despreading

structure are not satisfied. Therefore, the interferences limit the performance for larger values of Eb/N0. When the performance achieved by the frequency despreading structure is not satisfactory, transceiver performance can be improved by using another procedure for the design the coefficients F^P , more sophisticated than that in Equation 82. An alternative approach would require to introduce a time-domain filtering of each component of the DFT output, which would however also increase the receiver latency (however, the introduced latency would not be worse than that required by a lighter time-domain filtering approach following the PPN structure).

6 Conclusions

We have considered the FS-FBMC multicarrier scheme to meet requirements of cognitive radio such as high level of adjacent channel leakage and asynchronous communications. This structure is based on an FFT whose size is the length of the prototype filter and, then, the data samples are spread over several carriers. It has been shown that the FS-FBMC scheme is equivalent to the standard FBMC scheme in terms of input-output relation when operating on the flat channel and computationally more complex. However, it offers a crucial advantage in multi-path channel where the FS structure can perform equalization in the frequency domain. The benefit of improved equalization with no additional delay for the user might well overweight the increase in complexity. In fact, the

□ Single-tap : > FS 0 OFDM

Figure 11 Comparison wrt Eb/N0 (powers ten times larger than that of the user of interest). The effect of the noise on the performance of the OFDM transceiver (with 2,048 subcarriers) and of the two structures for FBMC transceiver (with 512 subcarriers) on both the considered channels for K = 4 in an uplink scenario with the powers of the interfering users ten times larger than that of the user of interest. The continuous line refers to channel ITU-R Vehicular A while the dashed line refers to the ITU-R Vehicular B channel.

simulation results clearly show that the standard structure is not able to achieve satisfactory performance without using a multitap equalizer and, therefore, is not able, in multipath channel, to achieve the minimum transceiver latency achieved by the FS structure. Moreover, the equivalence proven in the paper shows that the increased complexity may be needed in the receiver only. Furthermore, it has been shown that the FS-FBMC transceiver can achieve, in a wide range of values of Eb/N0, a performance equivalent to that of an OFDM system with the same FFT size and similar latency. Therefore, in a cognitive radio transmission context, the considered FS-FBMC system can assure the following advantages in comparison with an OFDM transceiver:

• smaller guard bands in the frequency domain, which means improved performance of the multiplexing scheme. Such a performance advantage can be quantified in terms of the overall bit-rate resulting from the use of OFDM/OQAM or equivalently can be quantified by comparing the complexity of the two transceivers since OFDM may need a heavy digital filtering of the signal to be transmitted in order to maintain similar overall bit-rate, as pointed out in [2];

• robustness to residual timing and carrier frequency offsets;

• capability to fully exploit fragmented spectrum;

• more flexibility in multiuser exploitation since asynchronous users can be accommodated.

The main disadvantage of the FS-FBMC receiver is related to the higher rate at which the FFT has to be performed when M/K subcarriers are employed in the frequency-spreading structure in comparison with the classical OFDM transceiver employing M carriers. Here, however, a pipelining approach to the FFT implementation may significantly reduce the importance of such a disadvantage.

Endnotes

aThe adjacent channel leakage ratio (ACLR) is the ratio of the filtered mean power centered on the assigned channel frequency to the filtered mean power on the adjacent channel frequency.

bWe use the approximate expression 4M log(M)-6M + 8 as a first approximation and to better appreciate the behavior of the complexity; however, the exact calculation of the number of flops for IFFT is given [27] by 394Mlog(M) - 124M - 2log(M) + 2f (-l)log(M)(l -3/8log(M)) + 8.

c Also if a component of the vector w^ is null the condition existing on a flat channel can be still restored provided that we accept that the definition (77), which does not apply to such a case, is obviously meant in the sense that the corresponding component of the vector f> is null.

dThe delay of each symbol transmitted with OFDM/OQAM system is equal to the length of the

prototype filter while the delay of each symbol transmitted with OFDM system is equal to the length of the multicarrier symbol period. Thus, the FBMC and the OFDM transceivers have the same delay if the OFDM one has K times more sub-carriers than OFDM/OQAM.

Appendix A

By using Equations 7, 8, and 21, the kth component df of dLR) can be written as:

Moreover, zn,i in Equation B.1 for I — pK + k with p e {0,1,..., M - 1} and k e {1,2,..., K - 1}

— Y^ AR) r

zn,pK+k — cn,mGmodkm(pK+k-mK)

— E c<nRRm [GkS [p - m] + Gk-KS [modM(p + 1 - m)]]

(R)^ (R) cn,pGk + Cn,modM (p+1)Gk-K.

Nb+Ns-1

dR = T, M kk' g [k + (n - n')M]

n'=0 k'eA

k €{0,1,...,M - 1}

= E E c% eiM kk g [k + (n - n')M]

n'=n-(K-1) k'eA

= E cei2Mkk' g[k] + E cRX,eiM(k+M)k' g[k + M] k'eA k'eA

+ E cn\k-eig[k + 2M]

+ E cn)u>ei2l(k+3M)k'g[k + 3M] +...

+ E cf-K-1),k'ei2l(k+(K-1)M)k'g[k + (K - 1)M]

Using Equation 20 in A.1, Equation 19 directly follows. AppendixB

The Ith component zn,i of the vector zn can be calculated by circular convolution of the two vectors G and cn

ZnX = E c(KR GmodKM(i-m) i e {0,1,..., KM - 1}

m=0 M-1

= El cnRiGmodKM (i-qK) = El ^nqGmodKMiA-qK). q=0 qeA

Using Equation 33, zn,i in Equation B.1 for i = pK and p e{0,1,..., M - 1} can be written as

zn,pK ^ ' CnRm GmodKM (pK-mK) m—0 M-1

— Go E c(RRs [p - m] — cnpGo. (B.2)

Consequently, taking into account Equations B.2, B.3, and 21, for p e {0,1,..., M - 1} and k e{0,1,..., K - 1}, we can compactly write

Zn,pK+k = cVpGk + c® odM(p+1)Gk-K(1 - S[k] )

= falpGk + jmodM(p+1)4modM(p+1)Gk-K(1 - 5[k] ) = faRn,pGk + f+KmodM(p+1)Gk-K(1 - Sk] )

From Equation B.4, Equation 35 directly follows.

In the derivation, we have used the property that jmodM (p+1) = jp+1 that holds for p e {0,1,..., M - 1} under the condition, always satisfied, that M is multiple of 4. In fact, the property is trivial for p e {0,1,..., M - 2} since modM(p + 1) = p + 1 while, for p = M - 1, it is equivalent to the property jM = 1, which holds only for M integer multiple of 4.

Appendix C

From Equations 42 and 14, it follows that

J2r[nM + M/2 + i] e-jkM'g[i] i—o

+ r[(n + 1)M + M/2 + i] e-jk M ig[M + i] + ...

+ r [(n +K-1) M+M/2 + i] e-jkMig[ (K-1)M+i] 1.

By using Equations 13 and 53, let us introduce the vector

(I,g,n) A (I) Tn+i — Tn+i x ■

e {0,1,...,K - 1} (C.2)

and let us denote with r^ff! its ith component for i e {0, l,..., M - l}. Using such definitions, from Equation C.l, it follows that

M-l M-l

V+ <k = E ^ e-ik M q + E e-jkM q

q=0 q=0

+ E e-ikMq + •••

q=0 M-l

, V r(I

+ rn+o

n+(K-1),q e

From Equations 52 and C.3, it follows that

jA D« x w = M DFT [r(lg'n)] + M DFT ^fl] + ... + M DFT [rl+K-D ]

= M£ DFT [r<+ x ge]

The last equality can be equivalently re-written as in Equation 5l.

Appendix D

From Equations 42 and l4, it follows that for k e {0, l,..., M - l}

D% = Ajj-k ^ r[nM + M/2 + i] e

i=0 K -1

E Gk'e

k '=-(IK-1)

KAA ik

J2 G* J2 r[nM + M/2 + i]

e-jKM m+k ' )

k'=-(K -1)

where we have taken into account that the prototype filter is real (G-k = Gk). Using Equations 63 and 64, Equation D.l can be re-written as

_ • k-1

->(I) _ j j-k^M X^ C* J?(I)

n,k = Aj ¿-^ Gk Rn,modKM (kK+k* )

k'=-(K -1)

, ) k g {0,1,M - 1}.

Taking into account Equations 62 and 59, Equation D.2 can be written as in Equation 61.

Competing interests

The authors declare that they have no competing interests. Author details

1 Dipartimento di Ingegneria Elettrica e delleTecnologie dell'Informazione, Université degli Studi di Napoli Federico II, via Claudio 21,80125 Napoli, Italy. 2CNAM-Electronique, 292 rue Saint-Martin, 75141, Paris cedex 03, France.

Received: 31 July 2014 Accepted: 31 December 2014 Published online: 10 March 2015

References

1. M Fitch, M Nekovee, S Kawade, K Briggs, R MacKenzie, Wireless service provision in TV white space with cognitive radio technology: a telecom operator's perspective and experience. IEEE Commun. Mag. 49(3), 64-73 (2011)

2. D Noguet, M Gautier, V Berg, Advances in opportunistic radio technologies for TVWS. EURASIP J. Wireless Commun. Netw. 2011(1), 170 (2011)

3. W Jiang, M Schellmann, in Proc. of IEEE Global Communications Conference (GLOBECOM2012). Suppressing the out-of-band power radiation in multicarrier systems: a comparative study (Anahem CA, USA, 2012)

4. H Bogucka, A. M Wyglinski, S Pagadarai, A Kliks, Spectrally agile multicarrier waveforms for opportunistic wireless access. IEEE Commun. Mag. 49(6), 108-115 (2011)

5. M Bellanger, in Proc. of2001 IEEE International Conference on Acoustics, Speech, and Signal Processing (ICASSP'01). Specification and design of a prototype filter for filter bank based multicarrier transmissions (Salt Lake City, UT, USA, 2001)

6. B Le Floch, M Alard, C Berrou, Coded orthogonal frequency division multiplex. Proc. IEEE. 83(6), 982-996 (1995)

7. B Farhang-Boroujeny, OFDM versus filter bank multicarrier. IEEE Signal Process. Mag. 28(3), 92-112 (2011)

8. KW Martin, Small sidelobe filter design for multitone data-communication applications. IEEE Trans. Circuits Syst.CAS-II. 45(8), 1155-1161 (1988)

9. M Bellanger, in Proc. of the 5th Symposium on Communications, Control and Signal Processing (ISCCSP 2012). FS-FBMC: an alternative scheme for filter bank based multicarrier transmission (Rome, Italy, 2012)

10. D Mattera, MTanda, M Bellanger, in 2012 International Symposium on Wireless Communication Systems. Frequency-spreading implementation of OFDM/OQAM systems (Paris, France, 2012)

11. E Kofidis, A. A Rontogiannis, in 2010 IEEE21st International Symposium on Personal Indoor and Mobile Radio Communications. Adaptive BLAST decision-feedback equalizer for MIMO-FBMC/OQAM systems (Istanbul, Turkey, 2010)

12. M Newinger, LG Baltar, AL Swindlehurst, JA Nossek, in 2010IEEE21st International Symposium on Personal Indoor and Mobile Radio Communications. MISO broadcasting FBMC system for highly frequency selective channels (Istanbul, Turkey, 2010)

13. T Ihalainen, A Ikhlef, J Louveaux, M Renfors, Channel equalization for multi-antenna FBMC/OQAM receivers. IEEE Trans. Vehicular Technol. 60(3), 2070-2085 (2011)

14. M Renfors, J Yli-Kaakinen, FJ Harris, Analysis and design of efficient and flexible fast-convolution based multirate filter banks. IEEE Trans. Signal Process. 62(15), 3768-3783 (2014)

15. M Renfors, J Yli-Kaakinen, in The Tenth International Symposium on Wireless Communication Systems 2013. Timing offset compensation in fast-convolution filter bank based waveform processing, (2013), pp. 1-4

16. P Siohan, C Siclet, N Lacaille, Analysis and design of OFDM/OQAM systems based on filterbank theory. IEEE Trans. Signal Process. 50(5), 1170-1183 (2002)

17. M Bellanger, J Daguet, TDM-FDM transmultiplexer: Digital polyphase and FFT. IEEE Trans. Commun. 22(9), 1199-1205 (1974)

18. B Hirosaki, An orthogonally multiplexed QAM system using the discrete Fourier transform. IEEE Trans. Commun. 29(7), 982-989 (1981)

19. H Sorensen, M Heideman, C Burrus, On computing the split-radix FFT. IEEE Trans. Acoust. Speech Signal Process. 34(1), 152-156 (1986)

20. H Saeedi-Sourck, Y Wu, JWM Bergmans, S Sadri, B Farhang-Boroujeny, Sensitivity analysis of offset QAM multicarrier systems to residual carrier frequency and timing offsets. Signal Process. 91(7), 1604-1612 (2011)

21. H Lin, M Gharba, P Siohan, Impact of time and carrier frequency offsets on the FBMC/OQAM modulation scheme. Signal Process. 102(9), 151-162 (2014)

22. T Ihalainen, TH Stitz, M Rinne, M Renfors, Channel equalization in filter bank based multicarrier modulation for wireless communications. EURASIP J. Adv. Signal Process. 2007(Article ID 49389), 18(2007)

23. International Telecommunication Union, in Recommendation ITU-R M. 1225: Guidelines for evaluation of radio transmission technologies for IMT-2000, (ITU, 1997)

24. X Hong, J Wang, C-X Wang, J Shi, Cognitive radio in 5G: a perspective on energy-spectral efficiency trade-off. IEEE Commun. Mag. 52(7), 46-53 (2014)

KM — 1

25. DL Wasden, H Moradi, B Farhang-Boroujeny, Design and implementation of an underlay control channel for cognitive radios. IEEE J. Select. Areas Commun. 30(10), 1875-1889 (2012)

26. M Shaat, F Bader, in Proceedings of the Fifth International Conference on Cognitive Radio Oriented Wireless Networks Communications (CROWNCOM). An uplink resource allocation algorithm for OFDM and FBMC based cognitive radio systems (Cannes, France, 2010)

27. FG Johnson, M Frigo, A modified split-radix FFT with fewer arithmetic operations. IEEE Trans. Signal Process. 55(1), 111-119 (2007)

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