O EURASIP Journal on

Advances in Signal Processing

a SpringerOpen Journal

RESEARCH Open Access

FBMC receiver for multi-user asynchronous transmission on fragmented spectrum

Jean-Baptiste Doré*, Vincent Berg, Nicolas Cassiau and Dmitri Kténas

Abstract

Relaxed synchronization and access to fragmented spectrum are considered for future generations of wireless networks. Frequency division multiple access for filter bank multicarrier (FBMC) modulation provides promising performance without strict synchronization requirements contrary to conventional orthogonal frequency division multiplexing (OFDM). The architecture of a FBMC receiver suitable for this scenario is considered. Carrier frequency offset (CFO) compensation is combined with intercarrier interference (ICI) cancellation and performs well under very large frequency offsets. Channel estimation and interpolation had to be adapted and proved effective even for heavily fragmented spectrum usage. Channel equalization can sustain large delay spread. Because all the receiver baseband signal processing functionalities are proposed in the frequency domain, the overall architecture is suitable for multiuser asynchronous transmission on fragmented spectrum.

Keywords: FBMC; FDM; FDMA; Equalization; Synchronization; Channel estimation; Asynchronous processing; Multiuser; Fragmented spectrum; Carrier frequency offset

1 Introduction

The current generation of cellular physical layers such as used by Long Term Evolution (LTE) and LTE Advanced have been optimized to deliver high-bandwidth pipes to wireless users but require strict synchronism and orthogonality between users within a single cell. The advent of the smart phone and the expected explosion of MachineType Communication (MTC) are posing new and unexpected challenges. Fast dormancy necessary for these types of equipment to save battery power has resulted in significant control signalling growth. Furthermore, as the availability of large amounts of contiguous spectrum is getting more and more difficult to guarantee, the aggregation of non-contiguous frequency bands has to be considered. Therefore relaxed synchronization and access to fragmented spectrum are key parameters for future generations of wireless networks [1]. This requirement of spectrum agility has encouraged the study of alternative multicarrier waveforms such as filter bank multicarrier (FBMC) to provide better adjacent channel leakage performance without compromising spectral efficiency [2]. frequency division multiple access (FDMA)

Correspondence: jean- baptiste.dore@cea.fr CEA Leti, 17 rue des Martyrs, Grenoble 38054, France

has already been considered for FBMC and provides promising performance without strict synchronization [3] between users. However, efficient frequency domain FBMC receiver algorithms have to be considered to fully benefit from the scheme. A first attempt at a practical implementation in the frequency domain for FBMC has been proposed in [4]. Synchronization and channel estimation are also based on the use of a training sequence but performance is limited by intercarrier interference (ICI) in the presence of large carrier frequency offset (CFO). This paper presents an architecture of a FBMC receiver suitable for asynchronous multiuser FDMA with very relaxed synchronization constraints: timing synchronization, CFO compensation and channel equalization are addressed.

Correction of CFO has been largely documented for orthogonal frequency division multiplexing (OFDM) [5] and FBMC [6]. A coarse CFO compensation in the time domain is necessary prior to a phase tracking in the frequency domain in order to limit ICI. As a different level of frequency offset is associated with each user, these solutions require as many receivers as users and are not suitable for practical implementation. Alternative signal processing methods for CFO compensation in the frequency domain should be developed. Iterative ICI cancellation has been proposed through iterative inversion

Springer

© 2014 Doré 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 work is properly credited.

techniques [7] or turbo processing [8] but is extremely complex to implement. The authors of [9] proposed an interesting method in the context of an orthogonal frequency division multiple access (OFDMA) uplink network. This method compensates CFO after fast Fourier transform (FFT) using a circular convolution. A good trade-off between performance and complexity is demonstrated. A low complexity CFO compensation method for FBMC - FDMA based on this method has been investigated in this paper.

Channel estimation is a necessary process before demodulation as the radio channel is frequency selective. For multicarrier modulation systems, the estimation is usually performed by sending a training data sequence, called pilot tones, on a set of carriers known to the receiver. The channel is then estimated at the pilot frequencies using the classical least squares (LS) or minimum mean-square error (MMSE) estimators [10]. When pilot tones are not available on every carrier, a process of interpolation is required to recover the complete channel response [11]. The interpolation may be done in the time domain [12]. A time-domain channel impulse response is obtained using an inverse Fourier transform of the channel estimated at the pilot frequencies. A filter may then be applied to reduce noise effects and border effects. An interpolated frequency domain channel is then processed by Fourier transform. Another solution may be to interpolate the channel in the frequency domain. Many algorithms, such as linear interpolation, low-pass filtering and spline cubic interpolation have been proposed [13,14]. For most interpolation schemes, the channel is poorly interpolated on the carrier located on the edges of the frequency band. This effect may be neglected when the number of carriers per contiguous frequency band is large but may lead to significant performance degradation of the overall system when the multicarrier modulation is applied to a fragmented spectrum. A robust interpolation scheme for the complete spectrum, including the edges is therefore critical. A solution has already been proposed in [12] using interpolation in the time domain. Performance is increased at the price of higher complexity. Furthermore, instability issues are not addressed and may lead to noise enhancement. A stable and robust scheme based on interpolation filters in the frequency domain is proposed in this paper.

Conventional equalization techniques for FBMC are suitable for multiuser reception. Good spectral containment introduced by the matched filter of the receiver helps to avoid distortion from non-synchronous users. Channel equalization may therefore be independently processed per user. However, if the conventional polyphase network FBMC implementation is employed, then equalization has to be carried out in the time domain. To cope with this issue, frequency-spreading FBMC is

considered [15]: at the cost of a significantly larger FFT (the size of the FFT is multiplied by the overlapping ratio), equalization may be efficiently done using a one-tap complex coefficient per subcarrier. A high-performance equalization scheme is described in this paper. Performance of the receiver is evaluated and discussed.

The paper is organized as follows: In Section 2, the overall context for a multiuser asynchronous environment is presented and the benefits FBMC waveforms have in this context are developped. Then in Section 3, the complete architecture of a FBMC receiver adapted to these scenarios is described and performance results are given. In Section 4, the main features and results are summarized and some perspectives are provided.

2 Asynchronous multiuser context

2.1 FBMC

Single-carrier (SC) FDMA has been chosen for the LTE to provide radio resource access for mobile users to the base station. This uplink multiuser access technique provides flexible resource allocation and is spectrally efficient. However, frequency offset between users should be strictly contained and received data signal should be aligned. If these synchronization requirements are not fulfilled, the orthogonality condition is not respected: intercarrier and intersymbol interference dramatically deteriorate performance. In order to guarantee mobility, the base station constantly monitors the time of arrival of received transmission signals. Signaling is sent to the mobile user via the downlink channel in order to synchronize this time of arrival at the base station. This constant exchange of control information introduces an overhead on the network that could be significant for asynchronous data communication services such as web access or machine-to-machine communications.

A multicarrier system can be described by a synthesis-analysis filter bank, i.e. a trans-multiplexer structure. The synthesis filter bank is composed of all the parallel transmit filters and the analysis filter bank consists of all the matched receive filters. The most widely used multicar-rier technique is OFDM, based on the use of inverse and forward FFT for the analysis and the synthesis filter banks. The prototype filter of OFDM is a rectangular window whose size is equal to the duration of the FFT. At the receiver, perfect signal recovery is possible under ideal channel conditions thanks to the orthogonality of the prototype filters. Nevertheless, under more realistic multipath channels, a data rate loss is induced by the mandatory introduction of a cyclic prefix (CP), longer than the impulse response of the channel. FBMC waveforms utilize a more advanced prototype filter design to better localize the subcarriers. The prototype filter used in this paper is based on the frequency sampling technique. This technique gives the advantage of using a closed-form

representation that includes only a few adjustable design parameters. The most significant parameter is the duration of the impulse response of the prototype filter also called overlapping factor, K. The impulse response of the prototype filter is given by [16]

k-1 , 2 k v

h(t) = Gp(0) + 2J2 (-1)k Gp(k) cos i KN (t + 1) J k=1 ^ '

where GP(0..3) =

1, 0.97195983, 1 - Gp(1)2

an overlapping factor of K = 4 and N is the number of carriers. In the following sections, the term carrier will refer to one of these N carriers and the term subcarrier will refer to one of the KN frequency domain FFT outputs at the receiver (see Section 3.1). The larger the overlapping factor K, the more localized the signal will be in frequency. In filter bank-based systems, transmit pulses are localized in time and in frequency. The orthogonality between carriers is maintained by introducing half a symbol period delay between the in-phase and the quadrature components of every complex symbol. The well-adjusted frequency localization of the prototype filter guarantees that only adjacent carriers interfere with each other. This justifies the use of FBMC waveforms in a non-synchronous context and particularly for the fragmented scenario. Nevertheless, adjacent carriers significantly overlap with this kind of filtering. In order to keep adjacent carriers orthogonal, real and pure imaginary values alternate on successive carrier frequencies and on successive transmitted symbols for a given carrier at the transmitter side. In order to maximize spectral efficiency of the offset QAM (OQAM) modulation, the symbol period T is halved.

FDMA access schemes have already been considered for FBMC modulations [3], but access on fragmented spectrum has not been envisaged yet. Fragmented spectrum may be viewed as the consequence of relaxed extension of channel aggregation for mobile communications. A user may access more than one frequency band being contiguous or not in frequency at any given time. The relaxed

synchronization between multiple users, i.e., users are not strictly synchronized in time, could lead to the generation of a heavily fragmented spectrum both in time and frequency. An example of the resulting spectrum is depicted in Figure 1. This phenomenon should be further exacerbated as messages get shorter as is often the case in machine-to-machine communications.

FBMC-FDMA access schemes appear therefore to be a very promising flexible multicarrier waveform. Efficient implementations of the FBMC-FDMA receivers should thus be considered in such multiuser asynchronous environments. A solution for this context is investigated in this paper.

2.2 Proposed burst structure

In order to keep a flexible frequency and time block allocation, a preamble-based burst approach is considered. Synchronization and channel estimation is performed using the training sequence. Its structure has been defined and is illustrated in Figure 2. It is composed of a preamble of duration P-FBMC symbols (P is set to 4 in Figure 2). The preamble has been designed to accurately detect the start of the burst and gives an estimate of the channel frequency response while preserving the localization properties of the FBMC signal. It is mainly composed of pilot carriers spaced every D active carriers for the whole duration of the preamble (D is set to 4 in Figure 2). The pilot carriers are designed so that the signal transmitted on each pilot carrier is constant for the duration of the preamble. Synchronization carriers are added on the first multicarrier symbol but are more sparsely distributed than pilot carriers. These are designed to accurately estimate the start of burst. By implementing at the receiver, all the baseband signal processing functions in the frequency domain, the proposed scheme may be extended without loss of generalization to the aforementioned multiuser asynchronous environment.

2.3 Notations

Bold letters denote vectors and matrices. Upper-case and lower-case letters denote frequency domain and time

Carriers A

^ Pilots carriers Q Synchronisation carriers ^ Data carriers

Figure 2 Considered FBMC burst structure.

domain variables, respectively. The following notations are used:

• (. )t Transpose

• (. )H Hermitian transform

• £[.] Expectation operator

• tr[.] Trace operator

• ||. || /-norm

F stands for the Nx N-DFT (discrete Fourier transform) matrix defined as:

1 2(N-1)

WNN1-1

(N-1)(N-1) N

where wN = ej n. Matlab notation was used to index the matrix. Therefore, A = B(:, 1 : U) means that A is built with the first U columns and all the rows of B.

3 Receiver architecture

3.1 Overview

A flexible architecture for multiuser asynchronous reception on fragmented spectrum is able to exploit the advantages of FBMC if the signal is efficiently demodulated in the frequency domain without a priori knowledge of the FFT timing alignment (i.e. the location of the FFT block, this property is called asynchronous FFT). A receiver architecture based on this assumption is depicted in Figure 3. An asynchronous FFT of size KN is processed every block of N/2 samples generating KN points, i.e. if rm is the mth received vector, a KN-point FFT is computed

for samples k = (n + m x N/2) with n = 0,1,..., NK — 1. These successive KN points are stored in a memory unit.

The detection of a start of burst is then achieved on the frequency domain (i.e. at the output of the FFT) using a priori information from the preamble. CFO is first estimated using the pilot subcarrier information of the preamble by computing the phase of the product between two consecutive FBMC symbols at the location of the pilot subcarriers. The propagation channel is assumed static for the duration of the burst. As described in [17], when large CFO correction is required, a first step in the estimation process consists of scanning the subcarriers around the pilot subcarrier locations to determine the subcarrier with the highest energy. A tracking algorithm of the CFO may complete the synchronization process when the duration of the burst is large and the accuracy of the preamble-based detection algorithm does not meet the required level [18]. CFO compensation is then performed in the frequency domain using a feed-forward approach.

The channel coefficients may be estimated on the pilot subcarriers of the preamble. Authors of [6,19] have already considered a similar approach by introducing a phase term to correct the CFO. In Section 3.2, this technique is completed by an efficient algorithm that compensates inter-carrier interference. The channel is then estimated on the pilot subcarriers before being interpolated on every active subcarrier. The use of a KN-point FFT makes the interpolation particularly specific to this receiver and a description of the proposed algorithm is detailed in Section 3.3.

Once the channel is estimated on all the active subcarriers, a one-tap per subcarrier equalizer is applied

Digital Front End

Memory

Synchro time FD

CFO Estimation

LLR calculation

Demapping

Prototype Filtering

Inner Receiver

Prototype Filtering

Outer decoder

Figure 3 Block diagram of the proposed receiver.

Equalizer

CFO Correction

Pilot Extraction

Channel Estimation

Channel Estimation

before filtering by the FBMC prototype filter (Section 3.4). Demapping and log-likelihood ratio (LLR) computation complete the inner receiver architecture. A soft-input forward error correction (FEC) decoder recovers finally the original message.

The asynchronous frequency domain processing of the receiver combined with the high stop-band attenuation of the FBMC prototype filter provides a receiver architecture that allows for multiuser asynchronous reception. FFT and Memory Unit are common modules, while the remaining of the receiver should be duplicated as many times as the number of parallel asynchronous users the system may tolerate.

3.2 Carrier frequency offset compensation 3.2.1 Problem formulation

When the received symbol rmFO is affected by CFO, the signal rmFO can be written as

rCFO =

J2 (dirine

j(2nS'm+&)

di [k, k] — e

- e>2n N

The considered scenario is an uplink asynchronous access with multiple different CFOs. Since CFOs associated with different users are different, correction must be performed separately for each user. The correction of the CFO for one given user cannot then be realized in the time domain with reasonable complexity as described by the authors of [9]. For this reason, we propose a frequency domain processing. Thanks to the very good frequency localization of the FBMC carriers, the asynchronous time reception of users does not cause any frequency interference between users (providing a frequency guard band of only one carrier). Thus, after FFT, the correction of the CFO in the frequency domain is described without loss of generality for one user. In the following, subscript index i is dropped. After KN fast Fourier transform, (3) becomes

R — CR ej(2nSm+®) _i_ 7

Rm — CRme + 7m

where Zm e CKNx1 is the additive white Gaussian noise vector and C e CKN xKN is a Toeplitz matrix defined by

C — FdF

where rmFO e <CKNx1 is the received time domain vector, is the time domain vector associated with user i, S' is the CFO for user i relative to the carrier spacing (assuming N carriers), is a random phase for user i and zm is the noise vector d1 e

ClŒxIŒ

a diagonal matrix

defined for user i by

The coefficients of C measure the level of intercarrier interference (ICI) and may be written as

C [', k] —

sinc(KS - (' - k)) n (N-i)iKs-(i-k))

(K s-(i-k)^

From Equation 7, it can be noted that when KS = q, q e Z, C is a subdiagonal matrix. q is denoted interger part of the CFO. The fractional part is labeled e. The parameter S

may then be decomposed into its integer and its fractional part:

8 = K+e

q e Z and e e R with e e [— 1/(2K); 1/(2K)]. Regardless of the integer part, q, ICI is only present when e = 0. Assuming q is known, the required range of detection for S may be very small. For instance, in the case of FBMC with a prototype filter of duration K = 4, the maximum level of ICI is introduced when e = 12.5%.

3.2.2 Proposed correction scheme

As described in the previous section, CFO may be decomposed into integer and fractional parts. The integer part is easily corrected by a shift of q subcarriers at the output of the KN-point FFT. The phase term should then be compensated by a phase correction factor. A simple and efficient way to reduce ICI introduced by the factional part of the CFO, e, may be achieved by complex filtering of the received sequence. In order to derive the complex coefficients of filter W that mitigates ICI, two criteria are considered:

• Zero-forcing criterion (): using Equation 5, and omitting the phase term, which has been corrected, we define W such as:

WC = I

By substituting C with the expression in (6), (9) may be rewritten as

WFdFH = I (10)

As FFH = FHF = I and ddH = I, the filter W may be derived by

W = FdH FH (11)

Minimum mean square error criterion (): the filter may also be optimized to minimize the mean square error (MSE), taking the noise level at the receiver into account:

West = argminw WRm - Rm

By taking the derivative of the expectation of the trace, the minimization problem becomes:

E [tr (WCRmRmHCHWH + WZmZmHWH -WCRmRm"-RmRmH CHWH)] = 0

In the presence of the additive white Gaussian noise (AWGN), E [ZmZmH] = all and if Hr is defined by

E [RmRmH] = ÖR, (13) becomes CÖrChWh + olwh - CÖr = 0

When = oRl, the solution is straightforward, and can be expressed by:

oR + o i

-FdH F

For both cases, the matrix W is Toeplitz and is therefore characterized by only KN complex coefficients. In general, is not a diagonal matrix but rather a band diagonal matrix and depends on the considered FBMC prototype filter. A closed-form expression of W may be obtained if and only if is invertible. However, W may be derived by a non-linear optimization process or by a using a pseudo-inverse based for instance on singular value decomposition.

The correction of the fractional part of the CFO with either or filters requires a KN complex-tap filter. Nevertheless the complexity may be significantly reduced as ICI introduced by the CFO rapidly decreases with the index of the interfering subcarriers (since |e| < 1/(2K)). In Figure 4, the power of the W coefficients for a 256-point FFT and a CFO of 10% is plotted. In that example, almost all the power is located around the diagonal. We proposed to approximate matrix W by a sparse band Toeplitz matrix with 2Q + 1 terms on each column and centred on the diagonal. This is done by forcing all the other coefficients to 0. A trade-off between complexity and accuracy of the frequency offset correction should be found for parameter Q. The 2Q + 1 coefficients may be extracted from W through equation (11) or (15).

3.2.3 Performance and practical implementation of the CFO correction algorithm

The performance of the proposed algorithm is analysed and a practical implementation is derived in this section using the following FBMC parameters: N = 1,024 and K = 4. Figure 5 gives the architecture of the proposed implementation. The first step, performed by the Shift module, consists in correcting the integer part q. The result is then filtered by the 2Q + 1-tap filter W. A phase correction performed on each FBMC symbol completes the CFO correction process.

The relative mean square error (RMSE) for S values from 1% to 12% and for different values of Q is given in Figure 6. RMSE has been defined by

RMSE = 10log10 ■

WRm - Rm

E [jjRmii2]

(dB) (16)

RMSE is a measure of the signal to interference level as seen on the constellation. In the example, a QPSK modulation has been considered and RMSE decreases as Q becomes larger.

When Q = 0, the algorithm only performs phase correction. The benefit of ICI mitigation is clearly demonstrated. Parameter Q could be chosen as a function of the system signal-to-noise ratio (SNR) so as to limit the system by thermal noise but not by interference.

In practice, the estimation of the CFO is never perfect. As illustrated in Figure 6, a residual CFO of 1% generated by the correction mismatch exhibits a RMSE power level below -30 dB and does not degrade performance significantly. This result is interesting to note for practical implementations. Indeed, the determination of matrix W requires complex computation (matrix inversion or non-linear optimization). Therefore, a possible practical implementation is to pre-compute a set of filters and use a combination of them. Figure 7 depicts the implemented filtering architecture based on a cascade of P = 4 pre-computed filters. Using the MUX command, it is possible to apply or bypass the filter correction. When a negative CFO correction should be applied, the INV command is activated. As filters for negative CFO are derived from filters for positive CFO through a simple permutation and complex conjugation of its coefficients, the operation is performed with little complexity.

For instance, Wi can be optimized for an incremental drift of 1% by

Si = 2i x 1%, i = 0 ...P - 1 (17)

In case of Q = 3 and P = 4, 28 complex coefficients must be stored. A correction range of [ -15%; 15%] is then possible. As an example, a CFO of 11.7% is corrected by

setting MUX0 = 1, MUX1 = 1, MUX2 = 0 and MUX3 = 1. The effective CFO corrected is 1% + 2% + 8% = 11%.

The RMSE performance of the proposed scheme is illustrated in Figure 8. We considered P = 4 filters with Q = 3 and using the optimization process of (17) and perfect phase estimation. A step shape of the curve is observed and comes from the finite resolution of the compensation (i.e. a step of S = 1%). When S > 12.5% a shift of p = 1 subcarrier is performed. These results demonstrate the trade-off between performance and complexity of the algorithm combined with the proposed cascade filters architecture.

A technique to mitigate ICI generated by CFO has been introduced in this section. All the operations have been carried out in the frequency domain and the scheme has given satisfactory performance results while being adapted for implementation.

3.3 Channel estimation and interpolation

Another important function that should be performed in the frequency domain is to recover the channel coefficients on each active subcarrier. This operation is achieved after the channel coefficients have been estimated on the pilot subcarrier location.

3.3.1 Problem formulation

One of the main advantages of multicarrier modulation techniques over single-carrier modulation is a greatly simplified equalization process. In the case of OFDM, as long as the duration of the channel impulse response is shorter than the guard interval and the channel is constant over

Phase correction

Filtering W

O ■i—i (O

Figure 5 Block diagram of the proposed CFO correction scheme.

the duration of the OFDM symbol, a frequency-selective wideband channel converts to a number of subcarrier channel with flat fading. For FBMC, Hirosaki [20] showed that this property may be preserved if the equalizer at each subcarrier channel is fractionally spaced.

Therefore, under these assumptions, the received vector R may be written after Fourier transform as

R = diag (H) X + Z

where diag (H) is a NH x NH diagonal matrix of the channel frequency response, X the NH x 1 vector of symbols to transmit and Z the NH x 1 vector of additive white Gaussian noise. NH is the number of active subcarriers.

We labelled UP, the indices of the pilots (located at the central location, as a KN-point FFT is considered at the front end of the receiver) and UH, the indices of the active subcarriers. Np is the number of pilot subcarriers. The LS

channel estimation of P, the pilot vector, at indices UP is given by [10]

R (Up [0]) R (Up [1]) R (Up [Np - 1]) '

P [0] P [1]

P N -1

Accurate estimation, H, of the channel coefficients H may be derived from the observation of P using an interpolation filter W of size NH x Np.

H = WP

To construct the filter W, the following minimization problem should be solved:

West = arg min E W

Let to be defined as to= E WP- H

WP - H

Then, since P is a least square estimate of P, the following equation may be written:

P = P + Zp (23)

where Zp is a Np x 1 noise vector. Equation 22 may then be rewritten as

to = E =E =E

||W(P + Zp) - H||2]

tr [(W (P + Zp) - H) (W (P + Zp) - H)"]] tr [WPP"W" - WPH" - HP"W"+ HH" + WZPZP"W" ]]

Then, if h is the N x 1 vector of the channel impulse response, H and P may be expressed as

P = F (Up, :)h = F^h (25)

H = F (Uh,:)H = Fuhh

Finally (22) may be rewritten as

to = tr [WFup $kFupHWH - WFup $kFuh

Fuh $hFUH

-Fuh $hFuph WH-

+o2 WW' Zp

where is the time domain channel autocorrelation matrix of size N x N and aZ2 the noise power. By taking the partial derivative of to with respect to W and making it equal to zero, (27) becomes

dto _ 0 d w = 0

(Fup ^hFup

H + io 2

WH - Fup ^hFuh

where I is the identity matrix of size Np x Np. In many cases, the matrix A is ill-conditioned and impossible to invert. However, a pseudo-inverse using a singular value decomposition (SVD) may be computed to derive W, but this method may lead to an unstable result. The power distribution of matrix W coefficients is given in Figure 9 for a channel distribution with a square delay profile in the time domain (NH = 256 and Np = 64).

In this example, the power distribution of matrix W coefficients is mainly located around the diagonal. It

should be noted that the coefficients located at the center of the matrix have very similar values. On the contrary, at the edges of the matrix, the coefficients are significantly different. This property has been exploited to construct a new optimization criterion. The criterion imposes complexity constraints on matrix W coefficients as a function of their location within the matrix.

From a practical point of view, implementing a complex filter with a large set of complex coefficients may

VO MUXO INV1 MU)

i__L 1 J

Filtering W Filtering W -

n n

Filtering W2 - Filtering W3 -

1 ^ T T ^ 1

INV2 MUX2 INV3 MUX3

Figure 7 Block diagram of the proposed filtering approach based on a cascade of pre-computed filters.

Figure 8 RMSE of the proposed filtering approach based on a cascade of P = 4 pre-computed filters with Q = 3.

be extremely costly. The overall complexity should often be kept under control in order to fit implementation area constraints. The following constraints have thus been added to the minimization problem:

• W is a matrix with real Q coefficients per row instead of complex Np coefficients per row

• The pilot carrier distribution follows a regular pattern, i.e. the sampling of the channel is uniform®

• At least Q pilots are active per set of contiguous subcarriers

The structure of matrix W is divided into three subfilter blocks, a left sub-filter block, a middle sub-filter

block and a right sub-filter block, so that matrix W may be rewritten as

r Wl 0

0 Wm 0

0 W m 0

0 Wr J

where Wi is the Ni x Q matrix representing the left subfilter block, Wm is the Nm x Q matrix representing the middle sub-filter block and Wr is the Nr x Q matrix representing the right sub-filter block.

The minimization problem of (27) may then be decomposed into three constraint minimization problems function of the three sub-filter blocks Wl, Wm and Wr. Further constraints on the coefficients may be added on the minimization algorithm to guarantee stability. This optimization process may be computationally demanding but allows for a control of the complexity of the implemented interpolation process.

Filters Wl, W m and Wr are obtained through the following minimization problem by using a priori knowledge of

A and o|:

minimize E

minimize E

minimize E

WiP (1 : Q)

H (1 : Ni)

WmP (Q + 1: Np - Q) -

H (Nl + 1 : Nh - Nr)

WrP (Np - Q + 1: Np) -

H (Nh - nr + 1: Nh )

By applying (27) into (30), the minimization may be realized considering an a priori knowledge of the autocorrelation matrix of the time domain channel impulse

response. By applying this constraint, the level of implementation complexity may be traded off against the target performance of the channel estimation module.

The optimization process is not necessarily implemented on the real time system. A stable and efficient interpolation scheme adapted to the channel conditions is provided while complexity is kept under control.

3.3.2 Application to the proposed preamble structure The estimation of the channel coefficients is performed before applying the FBMC prototype filtering. The LS estimates of the channel coefficients are computed by applying a CFO phase correction on the received signal according to the location of the pilot symbol on the prototype filter. Performance may be further improved by considering all the positions of the prototype filter and by applying a maximum ratio combining algorithm after interpolation.

Therefore, by considering LS estimates at the center of the FBMC prototype filter, the structure of the pilot pattern is depicted in Figure 10. A pilot subcarrier is located every KD active subcarriers. The first pilot is located on subcarrier K — 1. No pilot subcarrier is active on the right edge of the spectrum, therefore K — 1 + D channel coefficients should be derived. The number of active subcarriers is equal to Nh = 2(K — 1) + NpKD , where Np is the number of active pilots.

Under these assumptions, and considering Q is even, Nr, Nm and Ni may be expressed as follows:

Nl = k -1 + KDQ

Nm = KD

Mr = K - 1 + KDDQ - [K^]

The channel taps are assumed to be uncorrelated and as a consequence is a diagonal matrix. The most difficult task of the problem formulation consists in defining the a priori time domain channel autocorrelation delay profile. Different profiles may be considered and Figure 11 gives examples of possible distributions. For instance, an a priori autocorrelation $h profile adapted to single frequency network is characterized by a 0-dB echo channel

and is given in Figure 11d. The most conservative shape, the rectangular distribution (see Figure 11a), has been considered in the following section of the paper.

The optimization process is computationally intensive, making it difficult for a real time implementation. In order to maximize the interpolation performance while maintaining a reasonable level of complexity, the interpolation architecture depicted in Figure 12 has been considered. A

set of x filters is optimized according to a target SNR, time domain channel duration, time domain profile, etc. The choice of the filter to apply is controlled by an entity that estimates the received channel conditions and decides which filter is the most suitable to the channel conditions measured at the receiver. Estimated pilot subcarriers are fed through the three filter blocks and generate three interpolated channel estimates. The choice of the channel

Figure 12 Proposed architecture for adaptive frequency domain channel interpolation.

estimate is controlled by a multiplexer according to the subcarrier index.

3.3.3 Performance results

The following parameters have been taken for a numerical evaluation:

• K = 4

• D = 4

• N = 1,024

The RMSE of the channel frequency response versus the channel delay spread expressed in number of time domain samplesb for a SNR of 15 dB is given in Figure 13. The RMSE is defined by

RMSE = 10log10 ■

WP - H

E [||H||2]

Filters have been first optimized for Q = 10, a channel duration of 64 samples and a SNR of 15 dB. The performance is estimated using a Monte Carlo simulation based on the draw of 1,000 randomly generated channels with rectangular time profiles.

The RMSE is below the noise level for the three filter blocks. However, only a marginal gain is achieved for the right filter block. This is explained by the pattern considered in the pilot scheme; on the right edge, KD + K — 1 channel coefficients have to be interpolated while no pilot subcarriers are located at the edge.

From a practical point of view, it is interesting to determine the number of coefficients required depending on the channel delay spread. In Figure 14, the variation of RMSE versus channel delay spread for filter block designed for a channel delay spread length of 64 time domain samples, SNR = 15 dB and for various interpolation length Q is given. When filter interpolation length Q is increased, performance improves. Right filter blocks designed with Q = 4 and Q = 6 exhibit poor performance due to the channel delay spread. On the other hand, for left filter blocks, the variation on performance is much smaller and less than 1 dB between the case Q = 4 and the case Q = 10. Increasing interpolation length Q also increases the number of coefficients to store. This number of coefficient exhibits polynomial complexity in O(Q2) making large numbers of Q costly for implementation.

The influence of the SNR on the channel estimation is given in Figure 15. Three filter blocks have been designed for a channel delay spread length of 16 time domain samples and SNR values of 15, 25 and 60 dB. The distance between the thermal noise and RMSE is given as a function of the SNR. For low SNR of 0 dB, a gain of 5 dB is given by a filter block optimized for a SNR of 15 dB. As the noise level is increased, the filter blocks optimized for a SNR of 25 dB become more optimal up to around 36 dB, when the filter blocks optimized for 60 dB become more appropriate. These results showed the importance of using filter blocks with key characteristics designed around the actual working SNR in order to maximize channel estimation performance.

-2 -4 -6 -8

111 _12 <0

-16 -18 -20 -22

0 20 40 60 80 100 120 140 160 180 200 Channel length

Figure 13 RMSE versus channel delay spread for a SNR of 15 dB. The filter blocks were designed for a channel delay spread length of 64 time domain samples and a SNR of 15 dB. Channel length in the figure stands for the channel delay spread expressed in number of time domain samples.

Figure 14 RMSE versus channel delay spread for a SNR of 15 dB for various interpolation filter lengths Q from 4 to 10. Channel length in the figure stands for the channel delay spread expressed in number of time domain samples.

To conclude this section, a set of filter blocks for a channel delay spread lengths of 16, 32 and 64 time domain samples at a SNR = 15 dB and with Q = 8 has been optimized. The RMSE across the active subcarrier band (this includes the left and right edges) as a function of the channel delay spread for a fragment of 262 active subcarriers has been plotted in Figure 16. When the channel

delay spread is short, the use of an adequate filter block improves the RMSE of the interpolated channel. This figure illustrates the gain offered by the filtering approach if the channel delay spread is known and/or estimated. A comparison between an interpolation structure based on the state of the art and the proposed filter block structure adapted to the edges of the active spectrum is also given.

Channel length

Figure 16 RMSE as a function of the channel delay spread for a SNR of 15 dB and for different interpolation filter blocks. Channel length in the figure stands for the channel delay spread expressed in number of time domain samples.

The proposed method gives a performance improvement of up to 7 dB on the RMSE of the estimated channel. Interpolation errors at the edge of the fragment dominate the RMSE. For large frequency fragment the gain becomes negligible. Performance improvement is mainly due to the small number of pilot subcarriers located in the frequency fragment making the proposed approach particularly adapted to multiuser asynchronous fragmented spectrum usage.

3.4 Inner receiver architecture 3.4.1 Proposed architecture

The proposed equalizer is placed directly after the FFT and before the FBMC prototype filtering. This strategy is suited for asynchronous processing of multiusers in the frequency domain using a shared FFT processor. Assuming the channel delay spread, L, small compared to the KN-point FFT, the symbols at the output of the FFT may be written as follows:

R(k, p) = H (k)X (k, p) + Z(k, p) + Z,f (k, p) (33)

where k is the subcarrier index, p the index of the FBMC symbol, H(k) the complex channel coefficients for subcarrier k, Z(k,p) the AWGN sample and Zf (k,p) the inter-symbol interference (ISI) contribution. The channel here is again assumed static for the duration of the burst. When the matrix is diagonal, the ISI term can

be omitted. This assumption is a good approximation when the delay spread of the channel is small compared to the size of the FFT. To illustrate the validity of this assumption, we depicted in Figure 17 the power of the channel matrix coefficients in decibels for a channel delay spread size of KN/64 and KN/16 time domain samples (KN is set to 512). An average over 1,000 channel realizations was done, each tap following a Rayleigh process.

When the channel is short, the power is located on the diagonal term validating the proposed assumption. When the channel is longer, power spreads over the diagonal term of the matrix, and as a result, performance of the one-tap equalizer starts showing its limitations.

The goal of the equalizer is to recover X(k,p) from the observation R(k,p) and the estimated channel coefficient H(k). The ZF or MMSE criteria are classically used and may be expressed as follows:

H * (k)

X (k, p) = ^TT) R(k, p) (34)

where Q(k) is the optimized factor according to the ZF or MMSE criterion:

Qzf(T) = |H (k) |2 (35)

9 o 2(k)

Qmmse(T) = |H (k) |2 + (36)

oX (k)

L= KN/64

50 100 150 200 250 300 350 400 450 500 L= KN/16

50 100 150 200 250 300 350 400 450 500 Figure 17 Average power of the channel matrix coefficients for different channel delay spread.

where aX(k) is the expectation of the power of X on subcarrier k and o| is the expectation of the power of Z on subcarrier k. For the following, noise power is assumed to be constant over all the subcarriers.

Once the signal is equalized, the samples are filtered by the FBMC prototype filter before down-sampling by a factor K. After payload extraction, OQAM inverse transform, LLR values are obtained. The LLR is defined as

LLR(b(m, n)) = log-

P b(m) = 1|Y (n)

P b(m) = 0|Y (n)

= log-

Y.a£Q.0n P Yn = a| Y (n)

Yn = a| Y (n)

where ^ (resp. is the set of symbols comprising a bit 0 (resp. 1) for bit n, Y(n) is the nth real symbol after filtering and inverse OQAM transform and b(m, n) is the mth bit mapped to the nth symbol. Y(n) is written as

Y (n) = T(n)Y (n) + Zy (n) (38)

where T(n) is a factor depending on the criterion of the equalizer and ZY(n) is the noise associated with symbol n of variance a2 :

ZY (n)

Table 1 Simulated system parameters

ZY (n)

G(p)2H(Sn(p))2

Q (Sn(p))2

filter coefficient

r(n) = £

G(p)H (Sn(p)) Q (Sn(p))

We assumed that the conditional pdf of Y(n) is Gaussian and is expressed by

Y (n)Y (n) = a

\[2noZ

Y(n) - T(n)a

By applying Bayes rules, assuming equal distribution of symbols and using; the max log approximation (log£i fa, = max,- log fa), equation (37) may be rewritten as

LLR(b(m, n)) a

V2K-1 G(p)2H(Sn(p))2

^p=0 Q(Sn(p))2

max f-|Y(n) - r(n^

Ue^m v

— max — |Y (n) — F(n)a\2 aenm V > j

Further, complexity simplification of (42) may be considered with negligible performance loss such as the techniques described in [21].

3.4.2 Simulations results

The proposed equalization scheme has been evaluated by simulation using the parameters derived from the LTE mode 10-MHz parameters and can be found in Table 1.

First, the effects of timing offset (i.e. location of the FFT) on the proposed equalizer has been simulated. The MSE of a QPSK constellation with ZF equalization has been

Parameter

where G(p) is the pth prototype expressed in the frequency domain, the noise variance and Sn(p) is a set of subcarrier indices. We group into this set the indices of subcarriers used for the computation of symbol n. T{n) may be expressed by

Prototype filter K G

Modulation N Fe

Active carriers Type

Nfec R

Decoding algorithm Channel estimation Q

{0.235147, V2/2,0.97196,1 0.97196, V2/2,0.235147,0}

1, 024 15.36 MHz 601 (9 MHz) QPSK

CC 1/2, k = 7 1,944 bits 2/3 SOVA

estimated using the proposed equalization technique for various time offset of the KW-point FFT location. The following channel has thus been considered:

hk (t) = 5 t -

where parameter k is the timing offset introduced by the FFT processor. An interpolation filter optimized for a SNR of 15 dB and a delay spread of 32 x 1/Fe was considered for channel estimation. The simulation results are depicted in Figure 18.

The MSE on the constellation depends on the timing offset. The worst-case MSE is found when k is equal to N/4; in this case, the receiver is affected the most by the interference from the previous or next multicarrier symbols. The best case comes as expected when the FFT is perfectly aligned, i.e. timing offset is close top x N/2 Wp. This confirms that the proposed frequency domain equalizer scheme combined with channel estimation does not require FFT synchronization and is therefore adapted to asynchronous multiuser reception.

The combined performance of the proposed receiver algorithm is then evaluated. A working SNR of 5 dB has been considered and the bit error rate (BER) after Viterbi decoding for various channel delay spread length L of time domain samples has been measured. The channel impulse response has been defined as

hL(t) = 1 + J2atSi t -

Figure 18 Effects of the timing offset on the RMSE on a QPSK modulation. k is the timing offset.

where L is the number of taps in the channel impulse response and ai are complex coefficients following a Rayleigh distribution. With these assumptions, the delay spread of the channel is equal to L/Fe. The resulting BER at the output of the receiver has been evaluated and averaged for 10,000 channel realizations. Figure 19 gives the

simulated performance of the receiver for various channel interpolation filters as a function of the channel delay spread.

The performance obtained with the interpolation filter optimized for a channel delay spread length of N/(4Fe) time domain samples is of particular interest as

Channel delay spread ( us )

0 3.3333 6.6667 10 13.3333 16.6667

::::::::::::::::::

..................

r-............."

- J......J..... ..................

m.......I...... ..................

: : : \'¿: : : : : : : : : : :: ::::::::::::::::::

— i............. ..................

>■ v

y-i \ • Interp filters - 1/32 Interp filters - 9/128

Interp filters - 1/8 Interp filters - 1/4

0 0.05 0.1 0.15 0.2 0.25

Channel delay spread (normalized by N/Fe)

Figure 19 Performance of the proposed receiver for various interpolation filters as a function of channel delay spread. Interpolation filters

were optimized for a SNR of 15 dB and a channel delay spread of respectively N/ (32Fe), 9N/ (128Fe), N/ (8Fe) and N/ (4Fe). Equivalent CP durations (guard interval, GI) have been added for reference.

it demonstrated the limitations of the proposed scheme. When the delay spread of the channel is greater than 0.15 x N/Fe, performance is limited by the proposed equalization scheme rather than the channel estimation processing as the one-tap equalizer becomes inefficient. On the other hand, for channel delay spread below 0.15 x N/Fe, one-tap equalization is sufficient and channel estimation performs beyond requirement when an appropriate interpolation filter is considered. These results compare with the performance of an OFDM receiver with a guard interval length of at least 1/8 x N. In Figure 19, the guard interval length for an equivalent OFDM system with a FFT size of N is given. Performance of the receiver may be further improved by considering multi-tap equalization.

The choice of the interpolation filter also impacted on performance. When a filter optimized for a channel exhibiting a large delay spread is applied to a channel with a short delay spread, a significant amount of noise is not filtered. Performance is then better if the interpolation filter is adapted to the actual channel delay spread.

The proposed frequency domain receiver architecture seems very attractive for asynchronous multiuser processing. Fair performance has been obtained by the proposed equalizer scheme; the receiver is robust for large delay spread environment. As a comparison, the normal LTE guard interval is set to 4.69A |is or approximately 1/14 x N.

4 Conclusion

This paper presented a novel architecture and algorithms for FBMC reception. All the baseband signal processing functions are implemented in the frequency domain and no strict synchronization requirement on the FFT, the first element of the receiver, is required. This asynchronous frequency domain processing of the receiver combined with the high stop-band attenuation of the FBMC prototype filter provides a receiver architecture that allows for multiuser asynchronous reception.

Particular attention has been paid to CFO compensation in order to relax synchronism requirements beyond one carrier spacing. An algorithm to mitigate ICI has been proposed and simulated. The performance of a practical reduced complexity implementation is simulated. Depending on the receiver target SNR, complexity may be traded off to keep RMSE introduced by CFO below thermal noise. Channel interpolation has also been carefully considered. As multiuser asynchronous FDMA generates heavily fragmented spectrum blocks, channel estimation should be optimized for the edges of the receiver active carrier bands. A performance improvement of up to 7 dB on the RMSE of the estimated channel may result on some simulated scenarios. Finally, a new equalizer scheme has been thoroughly presented. Its complexity is contained

while good performance for channel exhibiting large delay spread is achieved. As a comparison, using the 10-MHz LTE parameters, the receiver performs well for channels with delay spread of up to 8.3 | s. This compares with standard LTE that has been designed for channels with delay spread of up to 4.7 | s. Moreover, we have demonstrated that the proposed equalizer does not require FFT synchronization and is therefore adapted to asynchronous multiuser reception.

Further work should consider implementation performance and complexity estimation with a comparison to similar flexible systems using traditional OFDM techniques. Eventual cost overhead of FBMC implementation and finite precision effects in this context would further complete the study. At the system level, radio resource simulations are necessary to understand the benefits of the multiuser asynchronous approach. Its potential gains in system capacity and energy consumption compared to already existing solutions could be evaluated.

Endnotes

aThe proposed method also applies to non-uniform sampling, but the problem is more complex to formulate.

bChannel delay spread is the time delay between the arrival of the first received signal component and the last received signal component associated with a single transmitted pulse.

Competing interests

The authors declare that they have no competing interests. Acknowledgements

The research leading to these results has received funding from the European Commission's seventh framework program FP7 ICTCall8 under grant agreement 318555 also referred to as 5GNOW.

Received: 30 November 2013 Accepted: 13 March 2014 Published: 28 March 2014

References

1. G Wunder, M Kasparick, S ten Brink, F Schaich,T Wild, I Gaspar, E Ohlmer, S Krone, N Michailow, A Navarro, G Fettweis, D Ktenas, V Berg, M Dryjanski, S Pietrzyk, B Eged, 5GNOW: Challenging the LTE design paradigms of orthogonality and synchronicity, in IEEE77th VehicularTechnology Conference (VTCSpring), 2013 (IEEE, Piscataway, NJ, 2013)

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

3. Y Medjahdi, MTerre, D Le Ruyet, D Roviras, A Dziri, Performance analysis in the downlink of asynchronous OFDM/FBMC based multi-cellular networks. IEEE Trans. Wireless Commun. 10(8), 2630-2639 (2011)

4. TH Stitz, T Ihalainen, M Renfors, Practical issues in frequency domain synchronization for filter bank based multicarrier transmission, in 3rd International Symposium on Communications, Control and Signal Processing, 2008. ISCCSP 2008 (IEEE Signal Processing Society, Piscataway, NJ, 2008), pp. 411-416

5. M Speth, SA Fechtel, G Fock, H Meyr, Optimum receiver design for wireless broad-band systems using OFDM.I. IEEE Trans. Comm. 47(11), 1668-1677(1999)

6. TH Stitz, A Viholainen, T Ihalainen, M Renfors, CFO estimation and correction in a WiMAX-like FBMC system, in IEEE 10th Workshop on Signal Processing Advances in Wireless Communications, 2009. SPAWC'09 (IEEE Signal Processing Society, Piscataway, NJ, 2009), pp. 633-637

7. AF Molisch, MToeltsch, S Vermani, Iterative methods for cancellation of Intercarrler interference in OFDM systems. IEEE Trans. Veh. Tech. 56(4),

2158-2167 (2007)

8. S Salari, M Ardebilipour, M Ahmadian, J-P Cances, V Meghdadi, Turbo receiver design with carrier-frequency offset estimation for LDPC-coded MIMO OFDM systems, in The 9th International Conference on Advanced Communication Technology, vol. 3 (IEEE, Piscataway, NJ, 2007),

pp. 1911-1915

9. J Choi, Y-H Lee, C Lee, H-W Jung, Carrier frequency offset compensation

for uplink of OFDM-FDMA systems, in IEEEInternationalConferenceon Communications, 2000. ICC2000, vol. 1 (IEEE, Piscataway, NJ, 2000), pp. 425-4291

10. J-J van de Beek, O Edfors, M Sandell, SK Wilson, P Ola Borjesson, On channel estimation in OFDM systems, in IEEE45th Vehicular Technology Conference, /995, vol. 2 (IEEE, Piscataway, NJ, 1995), pp. 815-8192

11. P-YTsai, T-D Chiueh, Frequency-domain interpolation-based channel estimation in pilot-aided OFDM systems, in IEEE59th VehicularTechnology Conference, 2004. VTC2004-Spring, vol. 1 (IEEE, Piscataway, 2004),

pp. 420-4241

XG Doukopoulos, R Legouable, Robust channel estimation via FFT interpolation for multicarrier systems, in IEEE 65th VehicularTechnology Conference, 2007. VTC2007-Spring (Piscataway, NJ, 2007), pp. 1861-1865 S Coleri, M Ergen, A Puri, A Bahai, Channel estimation techniques based on pilot arrangement in OFDM systems. IEEE Trans. Broadcast. 48(3), 223-229 (2002)

M-H Hsieh, C-H Wei, Channel estimation for OFDM systems based on comb-type pilot arrangement in frequency selective fading channels. IEEE Trans. Consum. Electron. 44(1), 217-225 (1998) M Bellanger, FS-FBMC: An alternative scheme for filter bank based multicarrier transmission, in 5th International Symposium on Communications Control and Signal Processing (ISCCSP), 20/2 (IEEE Signal Processing Society, Piscataway, NJ, 2012), pp. 1-4

16. M Bellanger, D LeRuyet, D Roviras, M Terré, J Nossek, L Baltar, Q Bai, D Waldhauser, M Renfors,T Ihalainen, A Viholainen,TH Stitz, J Louveaux, A Ikhlef, V Ringset, H Rustad, M Najar, C Bader, M Payaro, D Katselis, E Kofidis, L Merakos, A Merentitis, N Passas, A Rontogiannis, STheodoridis, D Triantafyllopoulou, DTsolkas, D Xenakis, M Tanda, T Fusco, M Huchard, J Vandermot, A Kuzminskiy, F Schaich, P Leclair, A Zhao, FBMC physical layer: a primer (2010). http://www.ict-phydyas.org, Accessed December 2013

17. N Cassiau, D Kténas, J-B Doré, Time and frequency synchronization for downlink CoMP with FBMC, in The Tenth International Symposium on Wireless Communication Systems 20/3. ISWCS20/3 (IEEE, Piscataway, NJ, 2013), pp. 46-50

18. P Amini, B Farhang-Boroujeny, Packet format design and decision directed tracking methods for filter bank multicarrier systems. EURASIP J. Adv. Signal Process. 2010, 7-1711 (2010)

19. T Fusco, A Petrella, M Tanda, Data-aided symbol timing and CFO synchronization for filter bank multicarrier systems. IEEE Trans. Wireless Comm. 8(5), 2705-2715 (2009)

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

21. FTosato, P Bisaglia, Simplified soft-output demapper for binary interleaved COFDM with application to HIPERLAN/2, in IEEE International Conference on Communications, 2002. ICC2002, vol. 2, (2002), pp. 664-6682

doi:10.1186/1687-6180-2014-41

Cite this article as: Doré etal.: FBMC receiver for multi-user asynchronous transmission on fragmented spectrum. EURASIP Journal on Advances in Signal Processing 2014 2014:41.

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