Scholarly article on topic 'Partial interference cancellation with maximum likelihood sequence detection in FBMC spatial multiplexing system'

Partial interference cancellation with maximum likelihood sequence detection in FBMC spatial multiplexing system Academic research paper on "Electrical engineering, electronic engineering, information engineering"

Share paper

Academic research paper on topic "Partial interference cancellation with maximum likelihood sequence detection in FBMC spatial multiplexing system"

Zakaria and Le Ruyet EURASIP Journal on Advances in Signal Processing (2016) 2016:47 DOI 10.1186/s13634-016-0339-x

EURASIP Journal on Advances in Signal Processing

RESEARCH Open Access

Partial interference cancellation with maximum likelihood sequence detection in FBMC spatial multiplexing system

Rostom Zakaria* and Didier Le Ruyet


Intrinsic interference in filter bank multicarrier (FBMC) modulation prevents the maximum likelihood (ML) detection in spatial multiplexing (SM) system. This intrinsic interference is caused by the transmultiplexer impulse response in time-frequency domain. Solutions based on interference cancellation are not always effective because they may introduce error propagation. In this paper, we propose to study some receivers based on reducing the interference concerned by the cancellation. These solutions can be seen as a trade-off between the whole Viterbi (ML) detection and the whole interference cancellation. The principle of the proposed receivers is to partially estimate and cancel the interference, and the rest of the interference is processed by a low complexity Viterbi detector. We show in this work that the receiver performance depends on the transmultiplexer impulse response and on the choice of the partial Viterbi detector.

1 Introduction

Nowadays, the most well-known multicarrier modulation is the orthogonal frequency division multiplex (OFDM) modulation. Thanks to the cyclic prefix (CP) insertion, OFDM technique provides a high robustness against the multipath fading channels. However, since OFDM uses CP and also rectangular pulse shape filtering that causes some disadvantages concerning the spectrum control and efficiency, to tackle this issue, filter bank multicarrier systems with offset quadrature amplitude modulation (FBMC/OQAM) was proposed as an alternative to CP-OFDM. Many works on the comparison between FBMC/OQAM and OFDM can be found in the literature such as [1]. One of the features of FBMC/OQAM is the possibility of shaping subcarrier signals with well-localized prototype filter in time and frequency axis [1]. That is, sidelobes steeply decay and allow a better spectrum control and offer an important resistance against time and frequency misalignment. Thus, FBMC/OQAM enables asynchronous multiple access and reduces intercarrier interference since each subcarrier signal is spectrally confined in a band and has negligible


LEATITIA/CEDRIC Laboratory, CNAM, 292 Rue Saint Martin, 75141 Paris, France

interference to other bands. In cognitive radios, the filtering capability of FBMC systems makes them the perfect choice for filling in the spectrum holes [1]. In the other hand, thanks to its spectral confinements, FBMC-based systems can reduce their guard band at the frequency boundaries, thereby increasing the spectral efficiency. Moreover, the absence of CP and any guard interval in time domain also contributes to increase the spectral efficiency. Many research works have been carried out in the issues related to doubly spread channels, and it has been shown that FBMC modulations are far better choices when compared to OFDM [2-4]. For instance, the author in [1] clearly shows the outperformance of the FBMC over OFDM in terms of the resilience to the doubly dispersive channels under the same assumptions of channel mobility and symbol density.

In FBMC/OQAM, each subcarrier is modulated with an OQAM transmitting pulse-amplitude modulation (PAM) symbols at each half a period (T/2). The orthogonality condition is considered only in the real field [5]. Indeed, the data at the receiver side is carried only by the real components of the signal, and the imaginary parts appear as interference terms. This interference is caused by the data symbols transmitted in the neighborhood area in the time-frequency domain. Since the

Springer Open

© 2016 Zakaria and Le Ruyet. Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a linkto the Creative Commons license, and indicate if changes were made.

interference is orthogonal to the useful data symbols, data detection is easily performed when the channel is flat fading or lowly selective [6]. However, in some situations like when the channel is highly selective, the intrinsic interference may cause performance degradation. Several works based on sub-channel equalization were carried out to tackle to this issue as in [7, 8] where multi-tap equalization is used. In [9], the authors have performed a minimum mean square error (MMSE) per subcarrier equalization. Whereas in [10, 11] the authors propose to implement a per-subcarrier decision-feedback equalizer (DFE). The extension of the latter to the multiple-input multiple output (MIMO) spatial multiplexing (SM) system was addressed in [12]. The per-subchannel equalization in FBMC with MIMO systems was also studied in [13, 14]. Some other works [15-20] address the design of MIMO precoding and decoding techniques for the FBMC modulation. For low MIMO selective channels, linear equalizations such as zero forcing (ZF) or MMSE can be applied to FBMC as shown in [21].

The presence of the inherent interference causes problems when combining FBMC/OQAM with some MIMO techniques such as Alamouti space-time block coding (STBC) scheme and SM systems with maximum likelihood (ML) detection. Regarding the Alamouti scheme, its application in a straightforward manner to the FBMC makes an inherent interference appear that cannot be easily removed [22]. Many works have been carried out on this topic such as [22] where the authors show that Alamouti coding can be performed but only when it is combined with code division multiple access (CDMA). A pseudo-Alamouti scheme was introduced in [23] where it is combined with single-carrier FBMC using the cyclic prefix (CP). Another solution was proposed by Renfors et al. in [24] where the Alamouti coding is performed in a block-wise manner. We have proposed in [25, 26] an iterative Alamouti scheme for FBMC based on intrinsic interference cancellation. On the other side, ML detection in SM system, which is supposed to offer a diversity order equal to the number of the receive antennas [27], cannot be applied straightforwardly with FBMC due to the presence of the interference.

Full ML receivers, in principle, offer the best possible performance but require an impractically high complexity when the impulse response is long. Moreover, the intrinsic interference in FBMC is two-dimensional which further complicates the detection task in maximum likelihood sense. Interference cancellation approaches generally offer the possibility of removing interference with low complexity increase and without enhancing the level of noise already present in the received signal [28]. Interference cancellation schemes are essentially based on using preliminary decisions to estimate and cancel the interference. In [29], we have addressed the possibility to cancel the

interference before applying a local ML detection and proposed some receiver schemes. However, the obtained performance was limited and far from the optimum due to the error propagation. We have also proposed in [26, 30] to modify the FBMC/OQAM modulation by transmitting QAM symbols at each period of T instead of transmitting PAM symbols at each T/2. We referred to this modulation as FBMC/QAM. This modulation allows to reduce the intrinsic interference power but at the expense of the orthogonality [26]. Error propagation problem in iterative inter-symbol interference (ISI) cancellation has been addressed in many research fields. In the general case, the authors in [31] have established conditions under which the cancellation scheme is effective. The idea is to consider the interference as the sum of two terms. The first one is cancelled by using tentative decisions, and the second uncancelled one is considered by a maximum likelihood sequence equalizer (MLSE).

In this work, in order to counteract the error propagation in SM-FBMC and make the cancellation scheme effective, we were inspired by [31]. We apply the theory developed in [31] and propose some receiver structures that satisfy the conditions of the interference cancellation effectiveness. The receivers are based on partial interference cancellation followed by a Viterbi detector instead of an ML detector. The tentative detector is first used to only partially cancel the intrinsic interference. We have presented a part of this work in [32] where we have only focused on FBMC/OQAM with one prototype filter. In this paper, both FBMC/OQAM and FBMC/QAM modulations are addressed, and their performances are compared. Furthermore, each of the both modulations are analyzed with two different prototype filters.

The paper is organized as follows. We start in Section 2 by giving an overview on FBMC/OQAM modulation highlighting the issue of FBMC when it is combined with SM-ML detection. Then, in Section 3, we give a background on the use of tentative decisions to cancel the interference, and we present the principle of partial interference cancellation. Sections 4 and 5 are devoted to the analysis and adaptation of the partial interference cancellation to the FBMC/OQAM and FBMC/QAM, respectively. Simulation results of the different proposed receivers are presented in Section 6. Then, we finish by a conclusion in Section 7.

2 The FBMC/OQAM modulation

2.1 System model

At the transmitter side, the discrete time of the equivalent baseband FBMC signal is written as follows [5]:

s[m] = J^ J2ak'"g[m - nM/2] e!Mk(m-2), (1)

k=0 kgZ

where M is an even number of subcarriers, g [m] is the prototype filter taking values in real field, and j is the delay term which depends on the length (Lg) ofg[m]. For the sake of simplification, we can set, as in [5], D = Lg — 1 and Lg = KM, where K is the overlapping factor. The transmitted symbols are real-valued symbols which are the real or the imaginary parts of QAM symbols. The additional phase term is given by $k,n = 2 (n + k) — nnk. We can rewrite (1) in a simple manner M—1

s[m] [m], (2)

k=0 neZ

where gk,n[m] are the shifted versions of g[m] in time and frequency. In the case of a perfect channel without noise, the demodulated symbol over the k'th subcarrier and the n'th instant is determined using the inner product of s[m] and gk',n [m]

rk',n' = (s,gk',n') = Y^ s[m]gkn[m]


= J2 J2 J2ak,ngk,n[m] gkn[m].

m=—<x k=0 neZ

The transmultiplexer impulse response can be derived assuming null data except at one time-frequency position (k0, n0) where a unit impulse is applied. Then, Eq. (3) becomes

rk',n'= gk0,n0[m] g*k',n'[m]


= J2 g[m] g [m — AnM/2] e^^ —m)


Jn(Ak+k0)An-j n (Ak+An) _ r(k0) (4)

e e = 1 Ak,An (4)

where An = n' — n0, Ak = k' — k0 and r(£fAn is then the transmultiplexer impulse response coefficients. It is worth noticing that the impulse response of the transmultiplexer depends on k0. Indeed, because the multiplicative factor ejirkoAn = (—1)koAn, the sign of some impulse response coefficients (with An odd) depends on the parity of k0.

The prototype filters are designed in a way that they provide a well-localized spectrum and spread over only a few adjacent subcarriers. Several pulse shaping prototype filters g[m] can be used according to their properties. A survey on prototype filters can be found in [33]. In this paper, we consider two different pulse shapes with overlapping factor set to K = 4. The first one is that which was adopted in the European Physical Layer for Dynamic Access (PHYDYAS) project on FBMC [34] and is referred to as the PHYDYAS prototype filter. The second is the one based on the so-called isotropic orthogonal

transform algorithm (IOTA) function [3]. IOTA filter is designed according to the energy concentration criterion, it is obtained by orthogonalizing the Gaussian function to prevent interference to neighboring points in the lattice [2]. Thus, it targets to preserve the optimum concentration property of Gaussian filter and fulfills the Nyquist criterion and can yield isotropic response in time and frequency [33]. As for PHYDYAS filter, it is designed according to the criterion of rapid decaying of the side-lobes. It is shown in [35] that, for an overlapping factor of K = 4, the minimum stop-band attenuation is 39.86 dB, and the sidelobes fall at \a)\—-3 where w is the angular frequency. Although the decaying rate and the stop-band attenuation increase with increasing the overlapping factor K, the choice of K = 4 is considered as a good trade-off between the stop-band attenuation and the system delay caused by the time spread of the filter [36]. Tables 1 and 2 depict the transmultiplexer impulse response coefficients (TAkAn) for even subcarriers k of the PHYDYAS and IOTA filter, respectively. We should note that the trans-multiplexer impulse response is independent from the subcarrier number (M) as long as M is sufficiently high, M > 8. All the prototype filters g[m] are designed to satisfy the real orthogonality condition given by [5]

Re £ [m] gk,n[m] = Sk,k &n,n'. (5)

[m=—^ I

Let us consider the SISO FBMC transmission. When passing through the radio channel and adding the zero-mean and circular white Gaussian noise contribution b[m], Eq. (3) becomes [37]

+ J2 hk ,n ak,n gk,n [m] gkn [m]

(k,n)=(k',n') m=—<xi

+ Y^ b[m] gkn [m],


where h^n is the channel coefficient at subcarrier k and time index n' and the term Ik'n is defined as an intrinsic interference. Since the prototype filters have a finite

Table 1 FBMC/OQAM transmultiplexer impulse response (main part) using PHYDYAS filter

n — 3 n — 2 n — 1 n n + 1 n + 2 n + 3

k — 1 0.043/ 0.125/ 0.206/ 0.239/ 0.206/ 0.125/ 0.043/

k —0067j 0 —0.564/ 1 0.564/ 0 0.067/

k + 1 0.043/ —0.125/ 0.206/ —0.239/ —0.206/ —0.125/ 0.043/

Table 2 FBMC/OQAM transmultiplexer impulse response using IOTA filter

Rq - 3 no - 2 n0 - 1 n0 n0 + 1 n0 + 2 n0 + 3

k0 - 2 -0.0016/ 0 -0.0381/ 0 0.0381 / 0 0.0016/

k0 - 1 0.0103/ 0.0381/ 0.228/ 0.4411 / 0.228/ 0.0381 / 0.0103/

k0 -0.0182/ 0 -0.4411/ 1 0.4411 / 0 0.0182/

k0 + 1 +0.0103/ -0.0381/ 0.228/ -0.4411/ 0.228/ -0.0381/ 0.0103/

k0 + 2 -0.0016/ 0 -0.0381/ 0 0.0381 / 0 0.0016/

support in time and we assume that they are spectrally well-localized, then the energy of the impulse response is confined in a finite set around the considered symbol [6]. Consequently, we can consider that the intrinsic interference term depends only on this restricted set (denoted by ,n'). Moreover, for the sake of derivation simplicity, we consider that the channel is sufficiently low-selective such that we can assume that the channel coefficient hk,n is approximately constant at least over the summation set &k,n. Hence, we can write as in [37]

ak,vi + ^

ak,n X gk,n[m] gkn [m]

+ Yk,,

According to (5) and because ak,n is real-valued, the intrinsic interference I^n = ju^n is pure imaginary. Thus, we can write the demodulated signal as

rk',n' ^ hk',n' iak',n' + juk',n') + Yk',n'•

It is worth noting that the noise term yk'n defined in Eq. (6) is colored, and its autocorrelation is given by

E {Yk',n' YkJ = a2 J2 gk,n [m] gin [m] m=-<x

- a2r(k)

= a 1 k>-k,n'-n-

The last equality is obtained according to Eq. (4).

2.2 Problem statement

In MIMO spatial multiplexing context, when Nt antennas are used to transmit Nt real data symbols and Nr antennas are used to collect the transmitted signals, the FBMC demodulated signal at the tth receive antenna and at a given time-frequency position (k, n) can be expressed by

= L C (al + j»l) + YÍ i=l

We should recall that the above expression is only valid under the assumption of the low-selectivity of the channels, as assumed for (7), between the transmit and receive

antennas. Finally, the matrix formulation of the system can be expressed as follows:

rk,n = Hk,n (ak ,n + juk,n ) + Y k ,n

where Hk,n is an (Nr x Nt) channel matrix with the element hk) at the tth row and ith column, rk,n =

,(1) J2)

r ' rv ' r

k,n k,n k,n

>) a(2) a(N)

ak,n akn ■■■, ak,n

' (1) (2)

is the received vector, ak,n =

is the transmitted symbol vector,

is the vector of the inter-

(1) (2) Yk,n, Yk,n ,

■ Yk,n

is the

ference values, and yk,n = noise vector.

Linear equalization as ZF and MMSE in FBMC/OQAM context can be performed as described in [21] where a virtually transmitted vector ck,n = ak,n + ;uk,n is considered instead of the effective one. The vector rk,n represents the input of the equalizer having as output the equalized virtually transmitted vector Ck,n

ck,n = G¡k,nrk„

where Gk,n is the equalization matrix based on the ZF or MMSE criterion. Then, the real part retrieval of ck,n yields the real equalized data vector ak,n [21].

As for maximum likelihood (ML) detection, the presence of the interference vector term uk,n in (11) prevents the application of ML separately at each time-frequency position (k, n). This is because the interference terms take their values in a large set and depend on the transmitted data symbols in the neighborhood around the considered position (k, n). The global ML detector is the one that considers all the transmitted data symbols within a frame. Obviously, such a receiver implementation is by far impractical due to its huge complexity.

3 Overview on partial interference cancellation

In [31], the authors have considered the channel model depicted in Fig. 1, where f0 (ak,ak-i, •••,ak^+i) is a function of S data symbols and represents the target response expected by the receiver. f1 (ak+y, ak+Y-1, •••, ak-x+1) is a function of y + X data symbols and represents a small channel perturbation.


It should be noted that, in general, both/0 and/i may be nonlinear functions. The samples of the signal at the input of the receiver are

rk = /0 {ak, ak—1..., a.k—h+1)

+/1 {ak+Y, ak+Y —1, •••, ak—x+1) + wk, (13)

where wk is the noise contribution. The receiver is composed by a tentative detector producing tentative decisions and a main Viterbi detector which assumes that the channel is described anly by /0. Before performing the main Viterbi detector, the tentative decisions are used only to cancel the remaining ISI (RISI) represented by/1. The receiver scheme is depicted in Fig. 2.

Given the correct data sequence {ak}, let us define a wrong sequence jake) J where the wrong symbols aare determined by a given error event e with a= ak + ek. That is, for a value of k, if ek = 0, then the symbol a(e in

the sequence ja^"1 J is correct. Let us define also /0 a

A(e) =

a^ ,a^ 1, •••,a


(e) k—«+1

—f0 (ak, ak—1, •••, ak—&+1) >

*(e) = A0e), A1e), •••, AK— 1

where K is assumed to be the total number of transmitted symbols. The authors in [31] classified the error events in terms of their distance d0(e) in the absence of RISI (f1= 0), which is given, in the presence of correlated noise, by


d0(e) =


Fig. 2 Receiver scheme with ISI cancellation using tentative decisions

where R = E{ww^}/oW is the normalized noise auto-covariance matrix, with oW is the noise variance and w = [w0,w1,•••,wK—1]T. The events whose distance d0(e) is the minimum are called "first-order" error events. Similarly, events whose distance is the second smallest are called "second-order" error events and so on [31].

The conditions for which RISI cancellation is satisfying are summarized as follows [31]:

1. Errors affecting the main (Viterbi detector) and the tentative detector must be statically independent.

2. The RISI (described by/1) must be small enough to guarantee that the main Viterbi detector can make relatively reliable decisions even when the tentative detector makes a decision error and such that the tentative detector also makes relatively reliable decisions in spite of the ISI.

3. The distance of second-order and higher-order error events that could cause error propagation must be significantly larger than that of first-order error events.

4 Application to MIMO-SM FBMC/OQAM

Now that the conditions for effective RISI cancellation are summarized, we will attempt to apply them to FBMC. Hence, the problem is, essentially, how to select the functions /0 and /1 such that the three conditions cited above are fulfilled. The intrinsic interference in FBMC is seen as a two-dimensional intersymbol interference (2D-ISI). An extension of the works of Agazzi and Seshardi [31] to 2D-ISI channels was treated in [38] assuming that the noise samples are uncorrelated (which is not the case in FBMC). Hence, in general, the target response /0 may also represent a 2D-ISI channel. Then, a 2D-Viterbi algorithm is required to match with/0. Designing a 2D-Viterbi is quite challenging. Therefore, for simplicity reasons, we opted to set the additional constraint that the target response /0 must be one-dimensional and that f1 covers the rest of 2D-ISI.

Obviously, the receiver complexity depends essentially on the complexity of the Viterbi detector. Therefore, we have to choose a configuration with the least complex Viterbi detector that meets the conditions for effective RISI cancellation. We will select three configurations with different sizes of the target response/0. According to the decond condition,/0 must contain the largest coefficients (rAk An) in each configuration. Hence, from Table 1, the selected target responses are [32]

/(1>(ak ,n) = ak,n,

(2) (k) /0 (ak,n, ak,n—1) = ak,n + 1 0,—1 x ak,n—1,

,(3), \ _ r(k)

Jo (ak,n+1, ak,n, ak,n-1) — r 0,1 ak,n+1 + ak,n


The first configuration f0(1) ) corresponds to the whole ISI cancellation which has been studied in [29, 39]. Since, in FBMC, we have |r|2 — 2 [6], it is easy to calculate the power of the RISI (represented byfl) for each configuration.

Regarding the first condition, it is easily satisfied when the tentative detector is different from the main one (Viterbi) [31]. We recall that we consider the case of spatial multiplexing system. Then, we chose the MIMO-MMSE equalizer as the tentative detector.

The third condition concerns the spectrum distances d0(e) defined by (16). Hence, for each configuration ,

i e {1,2,3}^, we compare the non-minimum distances to the minimum one. Then, according to (16), we have first to determine the matrix R. Since we consider that the target response are one-dimensional and Viterbi algorithm is performed on each subcarrier "k", the noise vector w is then w — [ wk,n, Wk,n+1,..., Wk,n+K-1]r. Hence, the matrix R — E{ww^}/oW is only composed by the coefficients r^q, q e Z and is given by

R — —

E{wk,nwk,n} E{wknw\n+1} E{wk,nwk,n+2} E{wk,n+1w'k,n} E{wk,n+1w'*n+1} E{wk,n+1w'k,n+2} E{wk,n+2wk,n} E{wk,n+2W^,n+1} E{wk,n+2wk,n+2}

r(k) 1 0,-1 r(k) 1 0,-2

r(k) r0,1

r(k) r 0,-1

r(k) r0,2 r(k) r0,1

where the last equality is obtained thanks to Eq. (9).

In Tables 3 and 4, we summarize the values of the first- and second-order distances obtained by using (16) and also the power of the RISI for the three considered configurations and both the PHYDYAS and IOTA filters, respectively.

Let us denote by e1 and e2 the first- and second-order error events, respectively. We also define w(e) as the number of error position in the error event e. The values of w(e) are also shown in the table. We remark that the difference between the second-order and the first-order distances is almost the same for all the configurations (0.8 ± 0.03), so we consider (as considered in [31]) that the higher-order distances are sufficiently larger than the minimum distance for each configuration. Hence, condition 3 is fulfilled for the three configurations.

Now, we have only to determine the configuration(s) for which the second condition is satisfied. Unfortunately, the determination of the RISI power for which the cancellation starts to be effective (or equivalently, error propagation ceases) is not trivial and depends also on the noise variance o2 [31]. We will show by simulations that only the third configuration f0(3)) allows to obtain effective RISI cancellation.

As for the receiver complexity, it strongly depends on that of the Viterbi detector. When we consider a spatial multiplexing system with Nt transmit antennas, the Viterbi detector has to compute qixNt branch metrics, where q is the number of all possible symbols akn (constellation size) and i e {1,2,3} is the number of the taps inJ0(i). In order to reduce the receiver complexity, we can replace the Viterbi detection algorithm by the M-Algorithm [40] which keeps only a fixed number (J) of inner states instead of all the inner states (q(,-1)xNt). Hence, the M-algorithm has to compute only the J x qNt branch metrics. In the rest of the paper, we call the proposed receivers "PaIC/Viterbi" (for Partial Interference Cancellation with Viterbi detector) and can be followed by an index to indicate the considered configuration.

5 Application of PaIC/Viterbi receivers in FBMC/QAM

Aiming to reduce the intrinsic interference power, we consider in this section the FBMC system transmitting complex QAM symbols instead of real ones (OQAM) [30]. That is, instead of transmitting real symbols akn at each half a period (T/2), complex symbols sk,n are transmitted at each period (T). Therefore, the equivalent baseband

Table 3 Spectrum distances and RISI power for FBMC/OQAM with the PHYDYAS filter

First configuration f f01 M Second configuration (fo(2M Third configuration (/q(3)

First-order distance 2 1.8857 1.9189

Second-order distance 2V2 2.4936 2.6170

w(ei) 1 1 1

w(e2) 2 2 2

Power of the RISI 1 0.6819 0.3638

Table 4 Spectrum distances and RISI power for FBMC/OQAM with the IOTA filter

First configuration ( f0

Second configuration ( f0

Third configuration lf(_

First-order distance 2 1.8985 1.8871

Second-order distance 2V2 2.5395 2.6594

w(Cï) 1 1 1

W(€2) 2 2 2

Power of the RISI 1 0.8054 0.6109

discrete time FBMC/QAM signal can be deduced from Eq. (1) and given by

s[m] — J2 J2Sk'"g[m - M e2Mk{m-t)e)9k-n, (21)

k=0 neZ

where dk,n — $k,2n — nn + 2k. Consequently, the trans-multiplexer impulse response coefficients tpA can be also derived from Eq. (4) and expressed as

select the following three functions /0®, i = 1,2,3 that determine the Viterbi target responses

f0 (sk,n) — sk,n,

f02) {sk,n, sk-1,n) = sk-1,n + f—1,0 X sk,n—1,

f AkAn = J2 g[m - AnM] e> -m)en W+2*«).

It is worthy to note that, in contrast to fAkAn> the coefficients f Ak,An do not depend on the subcarrier index k. The impulse response coefficients f Ak,An for FBMC/QAM using the PHYDYAS and IOTA filters are depicted in Tables 5 and 6, respectively.

The main principle of the PaIC/Viterbi receiver remains the same as for conventional FBMC presented in Section 4. However, the configuration of the PaIC/Viterbi receiver has to be adapted to the new FBMC/QAM context. Indeed, we have seen that the function f0 that matches the Viterbi detector must contain the largest coefficients f Ak,An. According to the transmultiplexer impulse responses depicted in Tables 5 and 6, the three largest coefficients ({1, ±0.239;} for PHYDYAS and {1, ±0.4411;} for IOTA) in each table are aligned with the direction of the frequency axis. Therefore, we can

/0(3) (^+1,^ skn sk—1,n) = I11,0sk+1,n + sk,n + 1 — 1,0sk—1,n,

where/0^') obviously corresponds to the IIC-ML receiver [26] where the whole interference cancellation is performed. Hence, the advantage for PaIC/Viterbi receivers in FBMC/QAM lies in the fact that the Viterbi algorithms are performed in the frequency axis direction. Whereas in FBMC/OQAM, they are performed in the time axis direction. Consequently, from an implementation point of view, the PaIC/Viterbi receivers are less complicated for implementation in FBMC/QAM than in the conventional FBMC, because only one Viterbi algorithm has to be performed for each one received multicarrier symbol, whereas in FBMC/OQAM, we need to perform a Viterbi detector simultaneously for each subcarrier.

We have seen in Section 4 that the performance of PaIC/Viterbi receivers depend on the power of the residual ISI which is concerned by the interference cancellation and also depend on the values of the spectrum distances d0(e) given by Eq. (16). To calculate the values of d0(e), we need first to determine the normalized autocovariance matrix of the noise R. Since the Viterbi algorithm is performed in the direction of the frequency axis, the noise vector w is then w = [w0,n, w1,n, •••, wM—1,n]T. Hence,

Table 5 FBMC/QAM transmultiplexer impulse response using the PHYDYAS filter

n — 1 n n + 1

k — 1 0.125/ 0.239/' 0.125/

k 0 10 k + 1 — 0.125/ — 0.239/ — 0.125/

Table 6 FBMC/QAM transmultiplexer impulse response using the IOTA filter

n0 — 1 n0 n0 + 1

k0 — 1 0.0381/ 0.4411/ 0.0381/

k0 0 1 0 k0 + 1 —0.0381/ —0.4411/ —0.0381/

the matrix R = Ejww^}/&% is only composed by the coefficients rp,0, p e Z and is given by

E{wo,«w0,«} E{wi,„w0,«}

E{wo,nW*u E{wi,nwfn

E{W0,nWM-1,nl E{wi,„wM_i,K}

E{wM-i,«w0,«} E{wM_i,«w1,K} ■■■ E{wM_i,«wM_1,«} _

T_i,o r_2,0

r 1,0 1

r 2,0 r 1,0 1

r M_1,0 r M_2,0 rM_3,0

- ri_M,0 r2_M,0 r3_M,0

(for both FBMC/OQAM and FBMC/QAM) will be affected. Indeed, channel estimation error will lead additional noise to the detected signals. However, according to the considered system model given by Eqs. (8) or (10), the channel estimation error is uncorrelated with the intrinsic interference. Hence, since the proposed configurations /q \ i e {1,2,3} are only related to the FBMC intrinsic interference, then they will not be changed in the context of imperfect channel estimation. Nevertheless, the metric used in the Viterbi algorithm might be changed to improve the detection performance under imperfect channel knowledge assumption. It is worthwhile to note that a generalized robust (against channel estimation errors) ML decoder that yields near optimal BER performance was proposed in [41].

where M is the number of subcarriers. Then, the values of d0(e) can be calculated for each single error event e and for each one of the chosen receiver configurations (according to f(), i = 1,2,3). We remind that the error events are classified according to their spectrum distance d0(e). The error events corresponding to the smallest d0(e) are called first-order error events e\, and those corresponding to the second smallest d0(e) are called second-order error events e2. Tables 7 and 8 give the different values of d0(e) of the first- and second-order error events for each receiver configuration, respectively, for the PHY-DYAS and IOTA filters. We also show in the tables the corresponding power of the residual interference (<uRISI) and the number of error positions (w(e)) within the first and second error events.

It is worth noticing that in the third configuration using the IOTA filter, the values of the two smallest distances are practically the same. Moreover, the first- and second-order error events contain, respectively, 3 and 4 error positions (w(ei) = 3 and w(e2) = 4 ), whereas in the other cases, they only contain 1 and 2 errors, respectively. Consequently, even if the residual interference is perfectly removed, the bit error rate (BER) performance would be poor because the most likely error events are those that contain 3 or 4 errors.

Finally, we note that if we assume imperfect channel estimation, the performance of the proposed detectors

6 Simulation results

In the following, we provide the simulation results concerning the three configurations of the PaIC/Viterbi receivers for both FBMC/OQAM and FBMC/QAM schemes with the PHYDYAS and IOTA filters. Since the motivation of this work is to address the problem of optimum detection in spatial multiplexing with FBMC, we have considered the simple 2 x 2 spatial multiplexing scheme. We assume perfect channel knowledge at the receiver side, and the four Rayleigh sub-channels are spatially non-correlated and follow the Veh-A channel model [42]. The parameters of the Veh-A channel are [42]: delays = [0 300 700 1100 1700 2500] ns, powers = [0 _ 1 _ 9 _ 10 _ 15 _ 20] dB, and the velocity is v = 60 km/h. For the simulation settings, we considered that the sampling period is Ts = 100 ns, the carrier frequency is fc = 1 GHz, and the number of subcarriers is M = 512. Therefore, the maximum Doppler spread is given byfD = f, where c = 3 x 108m/s is the light speed, thus fD ^ 55.56 Hz. We should notice that with the above parameters, the approximation considered in Eq. (10) is still valid. The data symbols are offset quadrature phase-shift keying (OQPSK) modulated for FBMC/OQAM, and quadrature phase-shift keying (QPSK) modulated for FBMC/QAM. That is, the number of bits per data symbol is q = 2 for both FBMC/OQAM and FBMC/QAM. The system performance is assessed in terms of BER as function of the signal-to-noise ratio (SNR) per bit Eb/N0 and is compared

Table 7 Spectrum distances and RISI power for FBMC/QAM with the PHYDYAS filter

First configuration f f01 M Second configuration (f^j Third configuration (/q(3)

First-order distance 2 1.8857 1.7185

Second-order distance 2.5407 2.6520 1.9189

w(ei) 1 1 2

w(e2) 2 2 1

Power of the RISI .1771 0.1198 0.0626

Table 8 Spectrum distances and RISI power for FBMC/QAM with the IOTA filter

First configuration lf0

Second configuration (f-

Third configuration f

First-order distance 2 1.8985 1.7833

Second-order distance 2.3561 2.5395 1.7892

W(€2) 2 2 4

Power of the RISI 0.3956 0.2010 0.0065

to that of the conventional CP-OFDM with ML detector. The CP duration for CP-OFDM is set to TCP = 32Ts = 3.2 ¡.is. The considered SNR per bit for FBMC is given by

Eb/No = Nt

Es/q a«/2

where Es = E{\s[m] \2} is the transmitted signal energy in each transmit antenna, Nt = 2 is the number of transmit antennas, and on2 is the noise power in each receive antenna. As for CP-OFDM, we have

Eb/No = Nt

2Es MTs + Tcp 17Es

8a« '

Therefore, the performance curves of exhibit an SNR loss of about 10log10 (17/16) -

CP-OFDM * 0.26 dB.


Figure 3 depicts the performance of the MMSE equalizer (which is our tentative detector) and of the PaIC/Viterbi for the three considered configurations in FBMC/OQAM using the PHYDYAS filter. We clearly notice that PaIC/Viterbi-3 exhibits almost the same performance as CP-OFDM/ML, and that the RISI cancellation is effective.

Hence, the value of the RISI power given in Table 3 for the third configuration is sufficiently small so that condition 2 is satisfied.

As for the first and second configurations, the performance degradation compared to CP-OFDM/ML begins from about 12 dB. This relatively high degradation is due to the high values of the corresponding RISI powers causing error propagation.

Now, we consider only the third configuration fo^) since the RISI cancellation is effective, and we assess the BER performance when using M-Algorithm instead of the Viterbi one in order to reduce the complexity. In Fig. 4, we show the obtained performance of the M-Algorithm with two different values of J (J = 2 and J = 4). We notice that with J = 2, we have an SNR loss of about 2.5 dB compared to PaIC/Viterbi-3. This SNR loss is due to the suboptimality of the M-Algorithm when J is small. Moreover, PaIC/M-Algo with J = 2 provides a performance worse than the one provided by the tentative detector (MMSE) as long as the SNR is less than 12 dB. However, with J = 4, we can observe that PaIC/M-Algo exhibits almost the same performance as PaIC/Viterbi-3 but with much lower algorithm complexity (4 inner states instead of 16).

Tentative detector (MMSE) —e— PaIC/Viterbi-1 - - PaIC/Viterbi-2

Eb/N0 (dB)

Fig.3 Performance of PaIC/Viterbi receivers for 2 x 2 spatial multiplexing using FBMC/OQAM modulation with the PHYDYAS filter

Eb/N0 (dB)

Fig. 4 Performance of PaIC/M-Algo receiver for 2 x 2 spatial multiplexing using FBMC/OQAM modulation with the PHYDYAS filter

In the case where the IOTA filter is used, the BER performance curves are plotted in Fig. 5. We remark that no configuration among the three ones can provide the same BER performance as CP-OFDM/ML. We can observe, for all configurations, an SNR loss with respect to CP-OFDM/ML performance in the range where SNR is greater than 15 dB. This bad performance results are explained by the high residual interference powers presented in Table 4 that prevent the convergence to the optimal performance.


Figure 6 depicts the BER performance of the proposed receivers in FBMC/QAM using the PHYDYAS filter. We notice that all the three PalC/Viterbi configurations exhibit performance that converge to the CP-OFDM/ML performance. This is because the conditions of effective interference cancellation are fulfilled, namely the RISI powers are sufficiently small, and the first-order error events are separated enough from the other error events by their spectrum distance.

In the case where the IOTA filter is used, we can observe in Fig. 7 that PaIC/Viterbi-3 receiver with the IOTA filter suffers from an SNR loss of about 1.2 dB with respect to CP-OFDM/ML. This SNR loss is due to the errors caused by the Viterbi algorithm itself. Indeed, we have seen in Table 8 that the most likely error events are those that contain 3 or 4 errors, and the RISI power is negligible. That is, the BER performance curve is shifted up. In the other hand, since the effect of erroneous tentative decisions are practically negligible because the variance of the residual interference is very small (oRISI ^ 0.007), the tentative detector (MMSE equalizer) is not useful and can be withdrawn.

As for PaIC/Viterbi-2 and PaIC/Viterbi-1 receivers, we can notice that with the IOTA filter the BER performance

Fig. 6 Performance of PalC/Viterbi receivers for 2 x 2 spatial multiplexing using FBMC/QAM modulation with the PHYDYAS filter

does not reach the CP-OFDM/ML one. We obtain an SNR loss at BER = 10-2 with respect to CP-OFDM/ML of about 2.25 dB for PaIC/Viterbi-2, and 6.25 dB for PaIC/Viterbi-1. Hence, both latter configurations suffer from error propagation and the interference cancellation is not effective. This performance limitation is explained by the fact that the residual interference variances are not sufficiently small; according to Table 8, we have aRisi ^ 0.2 and aRISI ^ 0.4 for PaIC/Viterbi-2 and PaIC/Viterbi-1, respectively.

Because the orthogonality condition is lost in FBMC/QAM, we can observe the effect of the interference by comparing the BER-performance of the MMSE equalizer in both cases where the PHYDYAS or IOTA filter is used. We notice in Fig. 6 that the MMSE equalizer reaches the BER value of 10-1 at SNR = 15.5 dB, whereas in Fig. 7, we have BER = 10-1 at SNR = 9.5 dB. This is

Eb/N0 (dB)

Fig. 5 Performance of PalC/Viterbi receiver for 2 x 2 spatial multiplexing using FBMC/OQAM modulation with the IOTA filter

Tentative detector (MMSE)

- PaIC/Viterbi-1

- PaIC/Viterbi-2

- PaIC/Viterbi-3 CP-OFDM/ML

Eb/N0 (dB)

Fig. 7 Performance of PalC/Viterbi receivers for 2 x 2 spatial multiplexing using FBMC/QAM modulation with the IOTA filter

due to the fact that the BER floor in the IOTA case is much higher than the BER floor in the case where the PHYDYAS filter is used.

7 Conclusions

The intrinsic interference in FBMC is a 2D-ISI in the time-frequency plan. In order to avoid a full 2D-Viterbi detector, we have proposed in this paper to study a tradeoff between a whole interference cancellation receiver and a full 2D-Viterbi detector. The proposed receiver is composed by a tentative detector giving decisions which serve to partially cancel the interference, followed by a Viterbi detector matching to the non-cancelled interference. Three configurations were treated. The first one is called PaIC/Viterbi-1 and correspond to the whole interference cancellation. The second one is PaIC/Viterbi-2, where the Viterbi detector matches with the two largest coefficients of the transmultiplexer impulse response. The third one is PaIC/Viterbi-3 and the Viterbi detector matches with the three largest coefficients. We have shown that the performance of the PaIC/Viterbi receivers essentially depends on the power of the residual interference (RISI) which is not concerned by the Viterbi detector. When the RISI power is not sufficiently low, the receivers suffer from error propagation.

We have studied the proposed receivers in both FBMC/OQAM and FBMC/QAM. This latter can offer the best performance in some situation since the global intrinsic interference is reduced. For FBMC/OQAM, we have shown by simulations that only PaIC/Viterbi-3 with the PHYDYAS filter can provide satisfactory performance compared to CP-OFDM/ML one. Whereas for FBMC/QAM, all the three PaIC/Viterbi configurations with the PHYDYAS filter exhibit the same performance as CP-OFDM/ML. However, the interference coefficients of the IOTA filter do not allow any PaIC/Viterbi configuration to reach the optimal BER performance. Indeed, in spite of the lowest RISI power for PaIC/Viterbi-3 with the IOTA filter, the BER performance is limited because of the fact that the most likely error events in the Viterbi detector contain more than 3 errors.

Finally, since the Viterbi detection increases the computational complexity, we can replace the Viterbi detector by another one based on M-Algorithm in order to reduce the receiver complexity without performance deterioration.

As a future work, one can consider extending the present work to the context of imperfect channel state information. As we have aforementioned, in order to minimize the BER performance degradation, we should adapt the Viterbi algorithm and investigate an efficient Viterbi metric capable of achieving near optimal performance.

Competing interests

The authors declare that they have no competing interests.

Received: 8 October 2015 Accepted: 22 March 2016 Published online: 18 April 2016


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

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

3. M Alard, Construction of a multicarrier signal. Patent WO 96/35278 (1996)

4. T Strohmer, S Beaver, Optimal OFDM design for time-frequency dispersive channels. IEEE Trans. Commun. 51 (7), 1111-1122 (2003)

5. 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)

6. C Lélé, P Siohan, R Legouable, J-P Javaudin, in Power Line Communications and Its Applications, 2007. ISPLC '07. IEEE International Symposium on. Preamble-based channel estimation techniques for OFDM/OQAM over the powerline (IEEE, Piscataway, NJ, 2007), pp. 59-64

7. T Ihalainen, TH Stitz, M Renfors, in Circuits and Systems, 2005. ISCAS 2005. IEEE International Symposium on. Efficient per-carrier channel equalizer for filter bank based multicarrier systems, vol. 4 (IEEE, Piscataway, NJ, 2005), pp. 3175-3178

8. T Ihalainen, T Hidalgo Stitz, M Rinne, M Renfors, Channel equalization in filter bank based multicarrier modulation for wireless communications. EURASIP J. Appl. Signal Process. 2007(1), 140-140 (2007)

9. DS Waldhauser, LG Baltar, J Nossek, in Signal Processing Advances in Wireless Communications, 2008. SPAWC2008. IEEE 9th Workshop on. MMSE subcarrier equalization for filter bank based multicarrier systems (IEEE, Piscataway, NJ, 2008), pp. 525-529

10. Zs Kollar, G Peceli, P Horvath, in Advances in Cognitive Radio, 2011. COCORA. IEEE First International Conference on. Iterative decision feedback equalization for FBMC systems, vol. 7, (2011), p. 60056

11. LG Baltar, DS Waldhauser, JA Nossek, in Circuits and Systems, 2009. ISCAS 2009. IEEE International Symposium on. MMSE subchannel decision feedback equalization for filter bank based multicarrier systems, (2009), pp. 2802-2805

12. E Kofidis, AA Rontogiannis, in Personal Indoor and Mobile Radio Communications (PIMRC) 2010. IEEE 21st International Symposium on. Adaptive blast decision-feedback equalizer for MIMO-FBMC/OQAM systems, (2010), pp. 841-846

13. A Ikhlef, J Louveaux, in Communications, Computers and Signal Processing,

2009. PacRim 2009. IEEE Pacific Rim Conference on. Per subchannel equalization for MIMO FBMC/OQAM systems (IEEE, Piscataway, NJ, 2009), pp. 559-564

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

15. M Caus, AI Perez-Neira, Transmitter-receiver designs for highly frequency selective channels in MIMO FBMC systems. IEEE Trans. Signal Process. 60(12), 6519-6532 (2012)

16. I Estella, A Pascual-Iserte, M Payaro, in Future Network and Mobile Summit,

2010. OFDM and FBMC performance comparison for multistream MIMO systems (IEEE, Piscataway, NJ, 2010), pp. 1-8

17. M Caus, AI Perez-Neira, in Acoustics, Speech and Signal Processing (ICASSP), 2013 IEEE International Conference on. Multi-stream transmission in MIMO-FBMC systems (IEEE, Piscataway, NJ, 2013), pp. 5041-5045

18. M Payaro, A Pascual-Iserte, M Najar, in Wireless Conference (EW), 2010 European. Performance comparison between FBMC and OFDM in MIMO systems under channel uncertainty (IEEE, Piscataway, NJ, 2010),

pp. 1023-1030

19. M Caus, AI Perez-Neira, Multi-stream transmission for highly frequency selective channels in MIMO-FBMC/OQAM Systems. IEEE Trans. Signal Process. 62(4), 786-796 (2014)

20. Y Cheng, P Li, M Haardt, in IEEE International Conference on Acoustic, Speech and Signal Processing (ICASSP). Coordinated Beamforming in MIMO FBMC/OQAM systems (IEEE, Piscataway, NJ, 2014), pp. 484-488

21. M El Tabach, J-P Javaudin, M Helard, in Communications, 2007ICC'07. IEEE International Conference on. Spatial data multiplexing over OFDM/OQAM modulations (IEEE, Piscataway, NJ, 2007), pp. 4201-4206

22. C Lélé, P Siohan, R Legouable, The Alamouti scheme with CDMA-OFDM/ OQAM. EURASIP J. Adv. Sig. Proc. 2010,1-13 (2010)

23. H Lin, C Lélé, P Siohan, in Signal Processing Advances in Wireless Communications, 2009 SPAWC 09. IEEE 10th Workshop on.Apseudo Alamouti transceiver design for OFDM/OQAM modulation with cyclic prefix (IEEE, Piscataway, NJ, 2009), pp. 300-304

24. M Renfors, T Ihalainen, TH Stitz, in Wireless Conference (EW), 2010 European. A block-Alamouti scheme for filter bank based multicarrier transmission (IEEE, Piscataway, NJ, 2010), pp. 1031-1037

25. R Zakaria, D Le Ruyet, in Wireless Communication Systems (ISWCS 2013), Proceedings of the Tenth International Symposium on. On interference cancellation in Alamouti coding scheme for filter bank based multicarrier systems (IEEE, Piscataway, NJ, 2013), pp. 1-5

26. R Zakaria, D Le Ruyet, Intrinsic interference reduction in a filter bank-based multicarrier using QAM modulation. Phys. Commun. (0), (2013)

27. Rvan Nee, A van Zelst, G Awater, in VehicularTechnology Conference Proceedings, 2000. VTC2000-Spring Tokyo 2000 IEEE51st. Maximum likelihood decoding in a space division multiplexing system, vol. 1 (IEEE, Piscataway, NJ, 2000), pp. 6-10

28. A Gersho, TL Lim, Adaptive cancellation of intersymbol interference for data transmission. Bell Syst Tech. J. 60,1997-2021 (1981)

29. R Zakaria, D Le Ruyet, M Bellanger, in Wireless Conference (EW), 2010 European. Maximum likelihood detection in spatial multiplexing with FBMC (IEEE, Piscataway, NJ, 2010), pp. 1038-1041

30. R Zakaria, D Le Ruyet, Y Medjahdi, in Wireless Communication Systems (ISWCS), 2012 International Symposium on. On ISI cancellation in MIMO-ML detection using FBMC/QAM modulation (IEEE, Piscataway, NJ, 2012),

pp. 949 -953

31. OE Agazzi, N Seshadri, On the use of tentative decisions to cancel intersymbol interference and nonlinear distortion (with application to magnetic recording channels). IEEE Trans. Inf. Theory. 43(2), 394-408 (1997)

32. R Zakaria, D Le Ruyet, in Wireless Communication Systems (ISWCS), 2011 8th International Symposium on. Partial ISI cancellation with Viterbi detection in MIMO filter-bank multicarrier modulation (IEEE, Piscataway, NJ, 2011), pp. 322 -326

33. A Sahin, I Guvenc, H Arslan, A survey on multicarrier communications: prototype filters, lattice structures, and implementation aspects. IEEE Commun. Surv.Tutor. 16(3), 1312-1338 (2014)

34. M Bellanger, PHYDYAS team, FBMC physical layer: a primer (2010). Website: Accessed 2010

35. S Mirabbasi, K Martin, in Circuits and Systems, 2002. ISCAS 2002 IEEE International Symposium on. Design of prototype filter for near-perfect-reconstruction overlapped complex-modulated transmultiplexers, vol. 1 (IEEE, Piscataway, NJ, 2002), pp. I-821-I-824

36. MG Bellanger, in Acoustics, Speech, and Signal Processing, 2001 Proceedings. (ICASSP'01),2001 IEEEInternationalConferenceon. Specification and design of a prototype filter for filter bank based multicarrier transmission, vol. 4 (IEEE, Piscataway, NJ, 2001), pp. 2417-2420

37. D Katselis, E Kofidis, A Rontogiannis, S Theodoridis, Preamble-based channel estimation for CP-OFDM and OFDM/OQAM systems: a comparative study. IEEE Trans. Signal Process. 58(5), 2911-2916 (2010)

38. S Van Beneden, J Riani, JWM Bergmans, in Communications, 2006 ICC '06. IEEE InternationalConferenceon. Cancellation of linear intersymbol interference for two-dimensional storage systems, vol. 7 (IEEE, Piscataway, NJ, 2006), pp. 3173-3178

39. RZakaria, D Le Ruyet, in Personal Indoor and Mobile Radio Communications (PIMRC), 2011 IEEE22nd International Symposium on. On spatial data multiplexing over coded filter-bank multicarrier with ML detection (IEEE, Piscataway, NJ, 2011), pp. 1391-1395

40. J Anderson, S Mohan, Sequential coding algorithms: a survey and cost analysis. IEEE Trans. Commun. 32(2), 169-176 (1984)

41. B SimThian, A Goldsmith, in Communication, Control, and Computing (Allerton), 201149th Annual Allerton Conference on. Decoding for MIMO systems with correlated channel estimation errors, (2011), pp. 524-530

42. ITU-R M.1225, Guidelines for evaluations of radio transmission technologies for IMT-2000 (1997)

43. D Tse, P Viswanath, Fundamentals of Wireless Communication. (Cambridge University Press, New York, 2005)

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