Get 20M+ Full-Text Papers For Less Than $1.50/day. Start a 14-Day Trial for You or Your Team.

Learn More →

Decoding Schemes for FBMC with Single-Delay STTC

Decoding Schemes for FBMC with Single-Delay STTC Decoding Schemes for FBMC with Single-Delay STTC //// Hindawi Publishing Corporation Home Journals About Us About this Journal Submit a Manuscript Table of Contents Journal Menu Abstracting and Indexing Aims and Scope Article Processing Charges Articles in Press Author Guidelines Best Paper Awards Bibliographic Information Contact Information Editorial Board Editorial Workflow Free eTOC Alerts Reviewers Acknowledgment Society Affiliation Subscription Information Open Special Issues Published Special Issues Special Issue Guidelines Abstract Full-Text PDF Full-Text HTML Linked References How to Cite this Article Complete Special Issue EURASIP Journal on Advances in Signal Processing Volume 2010 (2010), Article ID 689824, 11 pages doi:10.1155/2010/689824 Research Article <h2>Decoding Schemes for FBMC with Single-Delay STTC</h2> Chrislin Lélé and Didier Le Ruyet (EURASIP Member) Electronics and Communications Laboratory, Conservatoire National Des Arts Et Métiers (CNAM), 75141 Paris, France Received 5 June 2009; Accepted 28 December 2009 Academic Editor: Markku Renfors Copyright © 2010 Chrislin Lélé and Didier Le Ruyet. This is an open access article distributed under the Creative Commons Attribution License , which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. Abstract Orthogonally multiplexed Quadrature Amplitude Modulation (OQAM) with Filter-Bank-based MultiCarrier modulation (FBMC) is a multicarrier modulation scheme that can be considered an alternative to the conventional orthogonal frequency division multiplexing (OFDM) with cyclic prefix (CP) for transmission over multipath fading channels. However, as OQAM-based FBMC is based on real orthogonality, transmission over a complex-valued channel makes the decoding process more challenging compared to CP-OFDM case. Moreover, if we apply Multiple Input Multiple Output (MIMO) techniques to OQAM-based FBMC, the decoding schemes are different from the ones used in CP-OFDM. In this paper, we consider the combination of OQAM-based FBMC with single-delay Space-Time Trellis Coding (STTC). We extend the decoding process presented earlier in the case of 𝑁 𝑡 = 2 transmit antennas to greater values of 𝑁 𝑡 . Then, for 𝑁 𝑡 ≥ 2 , we make an analysis of the theoretical and simulation performance of ML and Viterbi decoding. Finally, to improve the performance of this method, we suggest an iterative decoding method. We show that the OQAM-based FBMC iterative decoding scheme can slightly outperform CP-OFDM. 1. Introduction Orthogonal Frequency Division Multiplexing (OFDM) is an efficient Multicarrier Modulation (MCM) capable of fighting against multipath fading channels. Its robustness to multi-path propagation effects comes from the insertion of a CP and is, therefore, obtained at the price of a reduced spectral efficiency. Furthermore, the rectangular shape of OFDM symbols leads to a s i n ( 𝑥 ) / 𝑥 frequency spectrum. Studies have been conducted in order to find better MCM schemes with respect to the frequency and/or time-frequency localization criteria. As suggested in [ 1 – 3 ], OFDM/OQAM also called as OQAM-based Filter Bank Multicarrier (FBMC) is an MCM scheme which may be the appropriate alternative. In OFDM/OQAM each subcarrier is modulated with Offset Quadrature Amplitude Modulation (OQAM). This principle has been introduced in [ 4 , 5 ], but it is only recently [ 1 ] that FBMC has been presented as a viable alternative to OFDM. Compared to OFDM that transmits complex-valued symbols at a given symbol rate, OQAM-based FBMC transmits real-valued symbols at twice this symbol rate. Therefore, a similar spectral efficiency is achieved by both systems. In practice, OQAM-based FBMC may provide a higher useful bit rate since it operates without the addition of a CP. Furthermore, with a pulse shaping that can be optimized according to given channel characteristics, its performance can be improved. However, all the interesting features of OQAM-based FBMC come at the price of a relaxation of the orthogonality conditions that only hold in the real field. At the receive side the data is carried only by the real component of the signal (assuming a 0 or 𝜋 / 2 phase modulation term). Thus, the imaginary part appears as an interference term. This interference term is a source of problem in the presence of the complex-valued channel as it destroys the real orthogonality. Therefore, when combining OQAM-based FBMC with MIMO technique such as Space-Time Block Codes (STBC) or Space-Time Trellis Coding (STTC) [ 6 , 7 ], the decoding process cannot be done in the same way as with CP-OFDM modulation. In the case of a single delay STTC chain with 2 transmit and 1 receive antennas, refrence [ 8 ] proposed a simple preprocessing to cancel this imaginary interference component. In this paper, we extend the proposed method in [ 8 ] to 𝑁 𝑡 transmit antennas and introduce an iterative decoding method. In Section 2 , we give a short description of the discrete-time OQAM modulation. Then, in Section 3 , we provide an overview of the STTC single delay detection. In Section 4.1 , we provide a theoretical performance analysis of ML and Viterbi decoding. Section 5 is devoted to the iterative decoding method in order to improve the performance of the previous decoding method. Simulation results are presented in Section 6 . Conclusions and perspectives are given in Section 7 . In the rest of the paper, FBMC will be used to denote OQAM-based FBMC. 2. The FBMC Modulation Using the baseband discrete-time model, we can write at the transmit antenna 𝑖 the OQAM-based FBMC signal as follows [ 1 ]: 𝑠 𝑖 [ 𝑚 ] = 𝑀 − 1  𝑘 = 0  𝑛 ∈ 𝐙 𝑑 𝑘 , 𝑛 , 𝑖 𝑔 [ 𝑚 − 𝑛 𝑁 ] 𝑒 𝑗 ( 2 𝜋 / 𝑀 ) 𝑘 ( 𝑚 − 𝐷 / 2 ) 𝑒 𝑗 𝜙 𝑘 , 𝑛                                                          𝑔 𝑘 , 𝑛 [ 𝑚 ] , ( 1 ) where 𝑀 = 2 𝑁 is the even number of subcarriers, 𝐹 0 = 1 / 𝑇 0 = 1 / 2 𝜏 0 is the subcarrier spacing, 𝜙 𝑘 , 𝑛 is an additional phase term, 𝑔 is the pulse shape, and 𝐷 is the delay parameter associated to the length of the pulse shape. The transmitted symbols 𝑑 𝑘 , 𝑛 , 𝑖 are real-valued data transmitted by antenna 𝑖 . They are obtained from a 2 2 𝐾 -QAM constellation, taking the real and imaginary parts of these complex-valued symbols of duration 𝑇 0 = 2 𝜏 0 , where 𝜏 0 denotes the time offset between the two parts [ 1 – 3 , 9 ]. For a given subcarrier 𝑘 and symbol time index 𝑛 , the real and imaginary parts are driven by the phase term 𝜙 𝑘 , 𝑛 given by 𝜙 𝑘 , 𝑛 = 𝜙 0 + 𝜋 2 ( 𝑛 + 𝑘 ) ( m o d 𝜋 ) , ( 2 ) where 𝜙 0 can be arbitrarily chosen. Here, we set 𝜙 0 = 0 and 𝑔 is assumed to be real valued. Assuming a distortion-free channel, a perfect reconstruction of real symbols is obtained owing to the following real orthogonality condition: ℜ   𝑔 𝑘 , 𝑛 ∣ 𝑔 𝑝 , 𝑞   = ℜ  ∞  𝑚 = − ∞ 𝑔 𝑘 , 𝑛 [ 𝑚 ] 𝑔 ∗ 𝑝 , 𝑞 [ 𝑚 ]  = 𝛿 𝑘 , 𝑝 𝛿 𝑛 , 𝑞 , ( 3 ) where 𝛿 𝑘 , 𝑝 = 1 if 𝑘 = 𝑝 and 𝛿 𝑛 , 𝑞 = 0 if 𝑛 ≠ 𝑞 . However, in practice for transmission over a realistic channel, the orthogonality property is lost, leading to intersymbol and intercarrier interferences. It has been shown in previous studies [ 8 ] that, when combining FBMC with single delay STTC in presence of 2 transmit and one receive antennas, specific processing should be done in order to remove the interference terms. In this paper, we will extend this method for 𝑁 𝑡 ≥ 2 antennas. 3. Single-Delay STTC in FBMC with 𝑁 𝑡 Transmit Antennas 3.1. Transmission Model Let us first assume that only the 𝑖 th antenna is transmitting. At the receiver side, the demodulated signal 𝑦 𝑘 , 𝑛 at the frequency 𝑘 and time instant 𝑛 ( 𝑛 𝜏 0 ) can be written as 𝑦 𝑘 , 𝑛 = 𝐻 𝑘 , 𝑛 , 𝑖 𝑑 𝑘 , 𝑛 , 𝑖 + 𝑗 𝐼 𝑘 , 𝑛 , 𝑖 + 𝜐 𝑘 , 𝑛 , ( 4 ) where (i) 𝐻 𝑘 , 𝑛 , 𝑖 is the channel coefficient between transmit antenna 𝑖 and the receiver, at subcarrier 𝑘 and time instant 𝑛 , (ii) 𝜐 𝑘 , 𝑛 is the noise component at subcarrier 𝑘 and time instant 𝑛 , (iii) 𝐼 𝑘 , 𝑛 , 𝑖 = ( − 𝑗 )  ( 𝑘  , 𝑛  ) ≠ ( 𝑘 , 𝑛 ) 𝐻 𝑘  , 𝑛  , 𝑖 𝑑 𝑘  , 𝑛  , 𝑖 ∞  𝑚 = − ∞ 𝑔 𝑘 , 𝑛 [ 𝑚 ] 𝑔 ∗ 𝑘  , 𝑛  [ 𝑚 ] . ( 5 ) We assume that we have a prototype filter well localized in time and frequency. This implies that in the previous equation the main contribution comes from the closest neighborhood, that is, 𝑔 𝑘 , 𝑛 [ 𝑚 ] 𝑔 ∗ 𝑘  , 𝑛  [ 𝑚 ] takes a significant value only for | 𝑘 − 𝑘 ′ | ≤ 1 and | 𝑛 − 𝑛 ′ | ≤ 1 . Moreover, if we assume that the channel is constant over a set of at least three consecutive subcarriers and a set of at least three consecutive time indexes, then we can rewrite the previous expression as in [ 10 ]: 𝐼 𝑘 , 𝑛 , 𝑖 ≈ 𝐻 𝑘 , 𝑛 , 𝑖 ( − 𝑗 )  ( 𝑘  , 𝑛  ) ≠ ( 𝑘 , 𝑛 ) 𝑑 𝑘  , 𝑛  , 𝑖 ∞  𝑚 = − ∞ 𝑔 𝑘 , 𝑛 [ 𝑚 ] 𝑔 ∗ 𝑘  , 𝑛  [ 𝑚 ]                                                                                  𝑢 𝑘 , 𝑛 , 𝑖 . ( 6 ) Thus, the demodulated signal can be approximated by 𝑦 𝑘 , 𝑛 ≈ 𝐻 𝑘 , 𝑛 , 𝑖  𝑑 𝑘 , 𝑛 , 𝑖 + 𝑗 𝑢 𝑘 , 𝑛 , 𝑖  + 𝜐 𝑘 , 𝑛 . ( 7 ) Throughout the remainder of the paper, we will consider ( 7 ) as the expression of the signal at the output of the demodulator. 3.2. Problem Statement Let us consider the single delay STTC scheme with 𝑁 𝑡 antennas as shown in Figure 1 . The real data to be transmitted is modulated by an FBMC modulator and transmitted by the first antenna. The same stream of data is delayed by 2 𝑛 𝑖 real data before being modulated by FBMC modulator and transmitted by the 𝑛 𝑖 th antenna. The delay 2 𝑛 𝑖 is chosen to have the same delay as with a CP-OFDM system although a delay of 𝑛 𝑖 could also be chosen. We denote by 𝑎 𝑘 , 𝑛 the real data from the main stream of data at frequency 𝑘 and time index 𝑛 . Thus, at a given subcarrier 𝑘 the transmission is given at antenna 𝑖 by 𝑑 𝑘 , 𝑛 , 𝑖 = 𝑎 𝑘 , 𝑛 − 2 𝑖 . At the receiver side, the demodulated signal can be written as Figure 1: FBMC Single-delay STTC transmitter. 𝑦 𝑘 , 𝑛 = 𝑁 𝑡 − 1  𝑖 = 0 𝐻 𝑘 , 𝑛 , 𝑖  𝑑 𝑘 , 𝑛 , 𝑖 + 𝑗 𝑢 𝑘 , 𝑛 , 𝑖  + 𝜐 𝑘 , 𝑛 , ( 8 ) where 𝜐 𝑘 , 𝑛 is the noise component at the subcarrier 𝑘 and time instant 𝑛 . As the same stream of data is transmitted over the 𝑁 𝑡 antennas, we have 𝑢 𝑘 , 𝑛 , 𝑖 = 𝑢 𝑘 , 𝑛 − 2 𝑖 , 0 = 𝑏 𝑘 , 𝑛 − 2 𝑖 . In the remainder of the paper, we will assume a channel constant over time, that is, ( 𝐻 𝑘 , 𝑛 , 𝑖 = 𝐻 𝑘 , 𝑖 ); we get 𝑦 𝑘 , 𝑛 = 𝑁 𝑡 − 1  𝑖 = 0 𝐻 𝑘 , 𝑖 ( 𝑎 𝑘 , 𝑛 − 2 𝑖 + 𝑗 𝑏 𝑘 , 𝑛 − 2 𝑖 )                              𝑥 𝑘 , 𝑛 − 2 𝑖 + 𝜐 𝑘 , 𝑛 . ( 9 ) The problem is to recover from 𝑦 𝑘 , 𝑛 the data 𝑎 𝑘 , 𝑛 . The presence of the term 𝑏 𝑘 , 𝑛 − 2 𝑖 makes the decoding process from 𝑦 𝑘 , 𝑛 difficult. Some processing should be carried out in order to recover the real data. 4. Interference Cancelation Method 4.1. Cancelation Procedure For the case 𝑁 𝑡 = 2 , it has been shown in [ 8 ] that if we define 𝑧 𝑘 , 𝑛 + 2 as 𝑧 𝑘 , 𝑛 + 2 = 𝐻 ∗ 𝑘 , 1 𝑦 𝑘 , 𝑛 + 𝐻 ∗ 𝑘 , 0 𝑦 𝑘 , 𝑛 + 2 , ( 1 0 ) then we have ℜ  𝑧 𝑘 , 𝑛 + 2  = ℜ  𝐻 ∗ 𝑘 , 1 𝑦 𝑘 , 𝑛 + 𝐻 ∗ 𝑘 , 0 𝑦 𝑘 , 𝑛 + 2  = | | 𝐻 𝑘 , 1 | | 2 𝑎 𝑘 , 𝑛 − 2 + 2 ℜ  𝐻 ∗ 𝑘 , 1 𝐻 𝑘 , 0  𝑎 𝑘 , 𝑛 + | | 𝐻 𝑘 , 0 | | 2 𝑎 𝑘 , 𝑛 + 2 + 𝑤 𝑘 , 𝑛 + 2 , ( 1 1 ) with 𝑤 𝑘 , 𝑛 + 2 = ℜ { 𝐻 ∗ 𝑘 , 1 𝜐 𝑘 , 𝑛 + 𝐻 ∗ 𝑘 , 0 𝜐 𝑘 , 𝑛 + 2 } . Let 2 𝐿 𝑓 denotes the frame length, for 𝑒 ∈ { 0 , 1 } . If we denote by 𝑡 𝑒 =  ℜ { 𝑧 𝑘 , 𝑒 } ℜ { 𝑧 𝑘 , 𝑒 + 2 } ⋯ ℜ { 𝑧 𝑘 , 𝑒 + 2 ( 𝐿 𝑓 − 1 ) }  𝑇 , 𝑎 𝑒 =  𝑎 𝑘 , 𝑒 𝑎 𝑘 , 𝑒 + 2 ⋯ 𝑎 𝑘 , 𝑒 + 2 ( 𝐿 𝑓 − 1 )  𝑇 , 𝑤 𝑒 =  𝑤 𝑘 , 𝑒 𝑤 𝑘 , 𝑒 + 2 ⋯ 𝑤 𝑘 , 𝑒 + 2 ( 𝐿 𝑓 − 1 )  𝑇 ⎡ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎣ | | 𝐻 𝑘 , 0 | | 2 0 ⋯ ⋯ 0 2 ℜ { 𝐻 𝑘 , 0 𝐻 ∗ 𝑘 , 1 } | | 𝐻 𝑘 , 0 | | 2 ⋱ ⋯ ⋮ | | 𝐻 𝑘 , 1 | | 2 2 ℜ  𝐻 𝑘 , 0 𝐻 ∗ 𝑘 , 1  | | 𝐻 𝑘 , 0 | | 2 0 ⋮ 0 ⋱ ⋱ ⋱ ⋮ ⋮ ⋱ ⋱ ⋱ 0 0 ⋯ | | 𝐻 𝑘 , 1 | | 2 2 ℜ  𝐻 𝑘 , 0 𝐻 ∗ 𝑘 , 1  | | 𝐻 𝑘 , 0 | | 2 ⎤ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎦                                                                                                                                                𝐺 2 , ( 1 2 ) ( ( ⋅ ) 𝑇 denotes the transpose operation and ( ⋅ ) 𝐻 the transpose conjugate one) then we have 𝑡 𝑒 = 𝐺 2 𝑎 𝑒 + 𝑤 𝑒 . ( 1 3 ) In this last equation, the imaginary interference term is canceled. Thus the decoding process can be easily carried out by using either Maximum Likelihood (ML) decoding, Viterbi decoding, or linear equalization such as Zero Forcing (ZF) or Minimum Mean Square Error (MMSE) decoding. More generally with 𝑁 𝑡 ≥ 2 , let us note and compute 𝑧 𝑘 , 𝑛 + 2 𝑁 𝑡 − 2 = 𝑁 𝑡 − 1  𝑝 = 0 𝐻 ∗ 𝑘 , 𝑁 𝑡 − 1 − 𝑝 𝑦 𝑘 , 𝑛 + 2 𝑝 = 𝑁 𝑡 − 1  𝑝 = 0 𝑁 𝑡 − 1  𝑖 = 0 𝐻 ∗ 𝑘 , 𝑁 𝑡 − 1 − 𝑝 𝐻 𝑘 , 𝑖 𝑥 𝑘 , 𝑛 + 2 𝑝 − 2 𝑖 + 𝑁 𝑡 − 1  𝑝 = 0 𝐻 ∗ 𝑘 , 𝑁 𝑡 − 1 − 𝑝 𝜐 𝑘 , 𝑛 + 2 𝑝                                    𝑛 𝑘 , 𝑛 + 2 𝑁 𝑡 − 2 = 𝑁 𝑡 − 1  𝑖 = 1 𝑖 − 1  𝑝 = 0 𝐻 ∗ 𝑘 , 𝑁 𝑡 − 1 − 𝑝 𝐻 𝑘 , 𝑖 𝑥 𝑘 , 𝑛 + 2 𝑝 − 2 𝑖                                                            𝐵 𝑘 , 𝑛 + 𝑁 𝑡 − 1  𝑖 = 0  𝑝 = 𝑖 𝐻 ∗ 𝑘 , 𝑁 𝑡 − 1 − 𝑝 𝐻 𝑘 , 𝑖 𝑥 𝑘 , 𝑛 + 2 𝑝 − 2 𝑖                                                            𝐴 𝑘 , 𝑛 + 𝑁 𝑡 − 1  𝑖 = 0 𝑁 𝑡 − 1  𝑝 = 𝑖 + 1 𝐻 ∗ 𝑘 , 𝑁 𝑡 − 1 − 𝑝 𝐻 𝑘 , 𝑖 𝑥 𝑘 , 𝑛 + 2 𝑝 − 2 𝑖                                                                𝐶 𝑘 , 𝑛 + 𝑁 𝑡 − 1  𝑝 = 0 𝐻 ∗ 𝑘 , 𝑁 𝑡 − 1 − 𝑝 𝜐 𝑘 , 𝑛 + 2 𝑝 . ( 1 4 ) Moreover 𝐴 𝑘 , 𝑛 is given by 𝐴 𝑘 , 𝑛 = 𝑥 𝑘 , 𝑛 𝜇 𝑘 , ( 1 5 ) and details for this equation are given in Appendix A.1 . The expression of 𝐵 𝑘 , 𝑛 is given by 𝐵 𝑘 , 𝑛 = 𝑁 𝑡 − 1  𝑞 = 1 𝑥 𝑘 , 𝑛 − 2 𝑞 𝛾 𝑞 , ( 1 6 ) where 𝛾 𝑞 are real-valued quantities which depend only on the channel coefficients as shown in Appendix A.2 . The expression of 𝐶 𝑘 , 𝑛 is given by 𝐶 𝑘 , 𝑛 = 𝑁 𝑡 − 1  𝑞 = 1 𝑥 𝑘 , 𝑛 + 2 𝑞 𝛽 𝑞 , ( 1 7 ) where 𝛽 𝑞 are real-valued quantities which depend only on the channel coefficients as shown in Appendix A.3 . Therefore, 𝑧 𝑘 , 𝑛 + 2 𝑁 𝑡 − 2 = 𝑁 𝑡 − 1  𝑞 = 1 𝛾 𝑞 𝑥 𝑘 , 𝑛 − 2 𝑞 + 𝜇 𝑘 𝑥 𝑘 , 𝑛 + 𝑁 𝑡 − 1  𝑞 = 1 𝛽 𝑞 𝑥 𝑘 , 𝑛 + 2 𝑞 + 𝑁 𝑡 − 1  𝑝 = 0 𝐻 ∗ 𝑘 , 𝑁 𝑡 − 1 − 𝑝 𝜐 𝑘 , 𝑛 + 2 𝑝 . ( 1 8 ) Thus, by noting that 𝑡 ( 1 ) 𝑘 , 𝑛 + 2 𝑁 𝑡 − 2 = ℜ { 𝑧 𝑘 , 𝑛 + 2 𝑁 𝑡 − 2 } , we have 𝑡 ( 1 ) 𝑘 , 𝑛 + 2 𝑁 𝑡 − 2 = 𝑁 𝑡 − 1  𝑞 = 1 𝛾 𝑞 𝑎 𝑘 , 𝑛 − 2 𝑞 + 𝜇 𝑘 𝑎 𝑘 , 𝑛 + 𝑁 𝑡 − 1  𝑞 = 1 𝛽 𝑞 𝑎 𝑘 , 𝑛 + 2 𝑞 + ℜ ⎧ ⎨ ⎩ 𝑁 𝑡 − 1  𝑝 = 0 𝐻 ∗ 𝑘 , 𝑁 𝑡 − 1 − 𝑝 𝜐 𝑘 , 𝑛 + 2 𝑝 ⎫ ⎬ ⎭                                                  𝑤 𝑘 , 𝑛 + 2 𝑁 𝑡 − 2 . ( 1 9 ) For 𝑒 ∈ { 0 , 1 } , we note 𝑡 𝑒 = [ 𝑡 𝑘 , 𝑒 𝑡 𝑘 , 𝑒 + 2 ⋯ 𝑡 𝑘 , 𝑒 + 2 ( 𝐿 𝑓 − 1 ) ] 𝑇 , 𝑤 𝑒 = [ ℜ { 𝑤 𝑘 , 𝑒 } ℜ { 𝑤 𝑘 , 𝑒 + 2 } ⋯ ℜ { 𝑤 𝑘 , 𝑒 + 2 ( 𝐿 𝑓 − 1 ) ] 𝑇 , and 𝐺 𝑁 𝑡 = ⎛ ⎜ ⎜ ⎜ ⎜ ⎜ ⎜ ⎜ ⎜ ⎜ ⎜ ⎜ ⎜ ⎜ ⎜ ⎝ 𝛽 𝑁 𝑡 − 1 0 ⋯ ⋯ ⋯ ⋯ ⋯ ⋯ ⋯ ⋯ 0 𝛽 𝑁 𝑡 − 2 𝛽 𝑁 𝑡 − 1 0 ⋯ ⋯ ⋯ ⋯ ⋯ ⋯ ⋯ ⋮ ⋮ ⋱ ⋱ ⋱ ⋯ ⋯ ⋯ ⋯ ⋯ ⋯ ⋮ 𝛽 1 ⋱ ⋱ ⋱ ⋱ ⋯ ⋯ ⋯ ⋯ ⋯ ⋮ 𝜇 𝑘 ⋱ ⋱ ⋱ ⋱ ⋱ ⋯ ⋯ ⋯ ⋯ ⋮ 𝛾 𝑁 𝑡 − 1 ⋱ ⋱ ⋱ ⋱ ⋱ ⋱ ⋯ ⋯ ⋯ ⋮ ⋮ ⋱ ⋱ ⋱ ⋱ ⋱ ⋱ ⋱ ⋯ ⋯ ⋮ 𝛾 1 ⋱ ⋱ ⋱ ⋱ ⋱ ⋱ ⋱ ⋱ ⋯ ⋮ 0 ⋱ ⋱ ⋱ ⋱ ⋱ ⋱ ⋱ ⋱ ⋱ ⋮ ⋮ ⋱ ⋱ ⋱ ⋱ ⋱ ⋱ ⋱ ⋱ ⋱ 0 0 ⋯ 0 𝛾 1 ⋯ 𝛾 𝑁 𝑡 − 1 𝜇 𝑘 𝛽 1 ⋯ 𝛽 𝑁 𝑡 − 2 𝛽 𝑁 𝑡 − 1 ⎞ ⎟ ⎟ ⎟ ⎟ ⎟ ⎟ ⎟ ⎟ ⎟ ⎟ ⎟ ⎟ ⎟ ⎟ ⎠ . ( 2 0 ) We have: 𝑡 𝑒 = 𝐺 𝑁 𝑡 𝑎 𝑒 + 𝑤 𝑒 . ( 2 1 ) There is no imaginary interference in ( 21 ) and consequently Maximum Likelihood (ML) [ 11 ] or linear equalizers can be used to estimate 𝑎 𝑘 , 𝑛 . The computation of 𝑧 𝑘 , 𝑛 from 𝑦 𝑘 , 𝑛 according to ( 14 ) is referred to as Preprocessing1 as shown in Figure 2 . We will now provide a theoretical performance analysis of this scheme. Figure 2: FBMC Single-delay STTC receiver. 4.2. A Theoretical Performance Analysis Let us consider that the noise 𝜐 𝑘 , 𝑛 is an AWGN noise with 𝐸 { | 𝜐 𝑘 , 𝑛 | 2 } = 𝑁 0 . It is worth noticing that ℜ { 𝑤 𝑘 , 𝑛 } is Gaussian noise as it is the result of the real part of a linear transformation of Gaussian noise. However this noise is colored. For example, when 𝑁 𝑇 = 2 , we have (i) 𝐸 { 𝑤 𝑘 , 𝑛 𝑤 ∗ 𝑘 , 𝑛 + 2 } = 𝐸 { 𝑤 𝑘 , 𝑛 + 2 𝑤 ∗ 𝑘 , 𝑛 } = 𝑁 0 ( | 𝐻 𝑘 , 0 | 2 + | 𝐻 𝑘 , 1 | 2 ) ℜ { ( 𝐻 𝑘 , 0 ) ∗ 𝐻 𝑘 , 1 } / 2 , (ii) 𝐸 { 𝑤 𝑘 , 𝑛 𝑤 ∗ 𝑘 , 𝑛 } = 𝑁 0 ( | 𝐻 𝑘 , 0 | 2 + | 𝐻 𝑘 , 1 | 2 ) / 2 = 𝑈 0 / 2 , (iii) for 𝑞 ≠ { 0 , 1 } , 𝐸 { 𝑤 𝑘 , 𝑛 𝑤 ∗ 𝑘 , 𝑛 + 2 𝑞 } = 0 . Let us recall that if the noise was white the ML performance would have been obtained by the Viterbi decoder. Therefore, the performance of Viterbi decoding in this present case is suboptimal. In [ 12 ] the authors evaluate the loss of performance of Viterbi decoding in presence of correlated noise. The optimal performance using an ML decoding is very complex to implement since it requires an exhaustive search over all the possible transmitted sequences. Another alternative could be to perform a whitening followed by a Viterbi decoding. However, such Viterbi decoding will be more complex since the whitening will increase the number of states. Indeed, the noise 𝑤 𝑒 is colored with a correlation matrix 𝑅 . Since 𝑅 is a positive Hermitian matrix, its eigenvalues are real and positive. We have 𝑅 = 𝑄 ⎡ ⎢ ⎢ ⎢ ⎢ ⎣ 𝜆 0 0 ⋯ 0 0 ⋱ ⋱ ⋮ ⋮ ⋱ ⋱ 0 0 ⋯ 0 𝜆 𝐿 𝑓 − 1 ⎤ ⎥ ⎥ ⎥ ⎥ ⎦                                    Λ 𝑄 𝐻 , ( 2 2 ) with 𝑄 being a unitary matrix, that is, 𝑄 𝑄 𝐻 = 𝐼 𝐿 𝑓 . We denote Λ 1 / 2 = ⎡ ⎢ ⎢ ⎢ ⎢ ⎣ 𝜆 1 / 2 0 0 ⋯ 0 0 ⋱ ⋱ ⋮ ⋮ ⋱ ⋱ 0 0 ⋯ 0 𝜆 1 / 2 𝐿 𝑓 − 1 ⎤ ⎥ ⎥ ⎥ ⎥ ⎦ . ( 2 3 ) Therefore, the whitening process can be done by computing 𝑦 𝑒 = Λ − 1 / 2 𝑄 𝐻 𝑧 𝑒 = Λ − 1 / 2 𝑄 𝐻 𝐺 2                  𝐻 𝑎 𝑒 + Λ − 1 / 2 𝑄 𝐻 𝑤 𝑒                  𝜇 𝑒 = 𝐻 𝑎 𝑒 + 𝜇 𝑒 . ( 2 4 ) It can easily be proved that 𝜇 𝑒 is AWGN. As we will see in the simulation results section, the presence of the colored noise will lead to a degradation of performance. Let us now present an iterative decoding approach which should improve the performance compared to that of the previous decoding strategy. 5. Iterative Method 5.1. Iterative Procedure In this section we propose an iterative decoding procedure for FBMC single-delay STTC decoding. At the output of the Preprocessing1 block (see Figure 3 ), we can perform a decoding procedure (ML, Viterbi, or linear decoding) to derive an estimate value ̂ 𝑎 ( 1 ) 𝑘 , 𝑛 of 𝑎 𝑘 , 𝑛 . From ( 6 ) and using this estimate ̂ 𝑎 ( 1 ) 𝑘 , 𝑛 , we can compute an estimate ̂ 𝑢 ( 1 ) 𝑘 , 𝑛 of 𝑢 𝑘 , 𝑛 by Figure 3: Receiver decoding processing for FBMC modulation in the case of single delay STTC transmission. ̂ 𝑢 ( 1 ) 𝑘 , 𝑛 =  ( 𝑝 , 𝑞 ) ≠ ( 0 , 0 ) ̂ 𝑎 ( 1 ) 𝑘 + 𝑝 , 𝑛 + 𝑞 ∞  𝑚 = − ∞ 𝑔 𝑘 , 𝑛 [ 𝑚 ] 𝑔 𝑘 + 𝑝 , 𝑛 + 𝑞 [ 𝑚 ]                                            𝛾 𝑝 , 𝑞 . ( 2 5 ) It is worth noticing that for a well-localized prototype filter in time and frequency domain it is enough to consider the previous sum only for 𝑝 , 𝑞 ∈ { 1 , − 1 } , that is, ̂ 𝑢 ( 1 ) 𝑘 , 𝑛 ≈  | | 𝑝 | | = 1 , | | 𝑞 | | = 1 ̂ 𝑎 ( 1 ) 𝑘 + 𝑝 , 𝑛 + 𝑞 𝛾 𝑝 , 𝑞 . ( 2 6 ) This approximation is justified in [ 10 ]. 𝛾 𝑝 , 𝑞 can be computed off-line since the prototype filter response is known. Then in ( 9 ) we can remove the contribution of the 𝑢 𝑘 , 𝑛 components by computing 𝑦 ( 2 ) 𝑘 , 𝑛 = 𝑦 𝑘 , 𝑛 − 𝑁 𝑡 − 1  𝑖 = 0 𝐻 𝑘 , 𝑖 ̂ 𝑢 ( 1 ) 𝑘 , 𝑛 − 2 𝑖 = 𝑁 𝑡 − 1  𝑖 = 0 𝐻 𝑘 , 𝑖 𝑎 𝑘 , 𝑛 − 2 𝑖 + 𝑁 𝑡 − 1  𝑖 = 0 𝑗 𝐻 𝑘 , 𝑖  𝑢 𝑘 , 𝑛 − 2 𝑖 − ̂ 𝑢 ( 1 ) 𝑘 , 𝑛 − 2 𝑖  + 𝜐 𝑘 , 𝑛 . ( 2 7 ) If we assume a perfect cancelation of the 𝑢 𝑘 , 𝑛 terms, that is, 𝑢 𝑘 , 𝑛 = ̂ 𝑢 ( 1 ) 𝑘 , 𝑛 , then we have 𝑦 ( 2 ) 𝑘 , 𝑛 = 𝑁 𝑡 − 1  𝑖 = 0 𝐻 𝑘 , 𝑖 𝑎 𝑘 , 𝑛 − 2 𝑖 + 𝜐 𝑘 , 𝑛 . ( 2 8 ) The operation of estimating 𝑢 𝑘 , 𝑛 and canceling its contribution to the signal 𝑦 𝑘 , 𝑛 is referred to as “Interference estimation + Interference cancelation” as depicted in Figure 3 . Thus, we can perform from 𝑦 ( 2 ) 𝑘 , 𝑛 a new decoding (Decoder 2 block) to obtain a new estimate ̂ 𝑎 ( 2 ) 𝑘 , 𝑛 of 𝑎 𝑘 , 𝑛 . In the same manner, we can use either a Viterbi/ML decoding or a linear decoder. From ̂ 𝑎 ( 2 ) 𝑘 , 𝑛 and ( 19 ) we can also compute 𝑡 ( 2 ) 𝑘 , 𝑛 + 2 by 𝑡 ( 2 ) 𝑘 , 𝑛 + 2 𝑁 𝑡 − 2 = 𝑁 𝑡 − 1  𝑞 = 1 𝛾 𝑞 ̂ 𝑎 ( 2 ) 𝑘 , 𝑛 − 2 𝑞 + 𝜇 𝑘 ̂ 𝑎 ( 2 ) 𝑘 , 𝑛 + 𝑁 𝑡 − 1  𝑞 = 1 𝛽 𝑞 ̂ 𝑎 ( 2 ) 𝑘 , 𝑛 + 2 𝑞 . ( 2 9 ) 𝑡 ( 2 ) 𝑘 , 𝑛 + 2 can also be rewritten as 𝑡 ( 2 ) 𝑘 , 𝑛 + 2 𝑁 𝑡 − 2 = 𝑁 𝑡 − 1  𝑞 = 1 𝛾 𝑞 𝑎 𝑘 , 𝑛 − 2 𝑞 + 𝜇 𝑘 𝑎 𝑘 , 𝑛 + 𝑁 𝑡 − 1  𝑞 = 1 𝛽 𝑞 𝑎 𝑘 , 𝑛 + 2 𝑞 + 𝑁 𝑡 − 1  𝑞 = 1 𝛾 𝑞  ̂ 𝑎 ( 2 ) 𝑘 , 𝑛 − 2 𝑞 − 𝑎 𝑘 , 𝑛 − 2 𝑞  + 𝜇 𝑘  ̂ 𝑎 ( 2 ) 𝑘 , 𝑛 − 𝑎 𝑘 , 𝑛  + 𝑁 𝑡 − 1  𝑞 = 1 𝛽 𝑞  ̂ 𝑎 ( 2 ) 𝑘 , 𝑛 + 2 𝑞 − 𝑎 𝑘 , 𝑛 + 2 𝑞  = 𝑁 𝑡 − 1  𝑞 = 1 𝛾 𝑞 𝑎 𝑘 , 𝑛 − 2 𝑞 + 𝜇 𝑘 𝑎 𝑘 , 𝑛 + 𝑁 𝑡 − 1  𝑞 = 1 𝛽 𝑞 𝑎 𝑘 , 𝑛 + 2 𝑞 + n o i s e c o m p o n e n t . ( 3 0 ) 𝑡 ( 2 ) 𝑘 , 𝑛 + 2 𝑁 𝑡 − 2 is a new version of the 𝑡 ( 1 ) 𝑘 , 𝑛 + 2 𝑁 𝑡 − 2 signal which is obtained from the estimates of the Decoder 2 block output. Thus, this last equation can be used to perform another estimation ̂ 𝑎 ( 3 ) 𝑘 , 𝑛 of 𝑎 𝑘 , 𝑛 in the same manner as we compute ̂ 𝑎 ( 1 ) 𝑘 , 𝑛 . We expect to improve the estimation of 𝑎 𝑘 , 𝑛 since the noise component in ( 30 ) should be less correlated than the one in ( 19 ). Again from ̂ 𝑎 ( 3 ) 𝑘 , 𝑛 we can derive an estimate ̂ 𝑢 ( 2 ) 𝑘 , 𝑛 of 𝑢 𝑘 , 𝑛 as in ( 25 ). Therefore, we can repeat another decoding process as already presented. We can run this decoding process as many times as necessary. The process of computing 𝑡 ( 2 ) 𝑘 , 𝑛 + 2 𝑁 𝑡 − 2 from the ̂ 𝑎 ( 2 ) 𝑘 , 𝑛 is referred to as Preprocessing2; see Figure 3 . Let us have a look at the convergence of this iterative method. 5.2. A Convergence Analysis of the Iterative Procedure Let us consider the function 𝑃 𝑒 = 𝐶 1 ( S N R ) that we obtain when considering the perfect cancelation of the interference term by using ( 28 ) and the function 𝑃 𝑒 = 𝐶 2 ( S N R ) obtained using ( 19 ). 𝑃 𝑒 is the real symbol error probability and S N R = 2 𝜎 2 𝑎 / 𝑁 0 = 1 / 𝑁 0 assuming that the real symbol power 𝜎 2 𝑎 is fixed at 1/2. These functions are illustrated in Figure 4 for a given channel realization. Let us note that 𝐶 1 is Δ dB better than 𝐶 2 , that is, Figure 4: Convergence illustration. 𝐶 1 ⎛ ⎜ ⎝ 1  1 + 𝛼 Δ  𝑁 0 ⎞ ⎟ ⎠ = 𝐶 2  1 𝑁 0  , ( 3 1 ) with Δ = 1 0 l o g 1 0 ( 1 + 𝛼 Δ ) . At the first iteration, when using ( 19 ) for decoding, we obtain at S N R = 1 / 𝑁 0 a symbol probability of error 𝑃 𝑒 1 = 𝐶 2 ( 1 / 𝑁 0 ) . This first iteration is summarized by the point 𝐴 1 ( 1 / 𝑁 0 , 𝑃 𝑒 1 ) in Figure 4 . Now, from this probability of error we can derive the degradation that we obtain when applying interference cancelation. Indeed, the cancelation of the interference will add some noise to the current noise component. This additional noise component is given by the cancelation error 𝑛 + = 𝑁 𝑡 − 1  𝑖 = 0 𝑗 𝐻 𝑘 , 𝑖  𝑢 𝑘 , 𝑛 − 2 𝑖 − ̂ 𝑢 ( 1 ) 𝑘 , 𝑛 − 2 𝑖  = 𝑁 𝑡 − 1  𝑖 = 0 𝑗 𝐻 𝑘 , 𝑖  ( 𝑝 , 𝑞 ) ≠ ( 0 , 0 )  𝑎 𝑘 + 𝑝 , 𝑛 − 2 𝑖 + 𝑞 − ̂ 𝑎 ( 1 ) 𝑘 + 𝑝 , 𝑛 − 2 𝑖 + 𝑞  × ∞  𝑚 = − ∞ 𝑔 𝑘 , 𝑛 − 2 𝑖 [ 𝑚 ] 𝑔 𝑘 + 𝑝 , 𝑛 − 2 𝑖 + 𝑞 [ 𝑚 ] . ( 3 2 ) Using the current observation ̂ 𝑎 𝑘 , 𝑛 = 𝑎 𝑘 , 𝑛 w i t h p r o b a b i l i t y 1 − 𝑃 𝑒 1 , ̂ 𝑎 𝑘 , 𝑛 ≠ 𝑎 𝑘 , 𝑛 w i t h p r o b a b i l i t y 𝑃 𝑒 1 ( 3 3 ) and considering that [ 10 ]  ( 𝑝 , 𝑞 ) ≠ ( 0 , 0 ) | | | | ∞  𝑚 = − ∞ 𝑔 𝑘 , 𝑛 − 2 𝑖 [ 𝑚 ] 𝑔 𝑘 + 𝑝 , 𝑛 − 2 𝑖 + 𝑞 [ 𝑚 ] | | | | 2 = 1 , ( 3 4 ) we have 𝐸  | | 𝑛 + | | 2  = 𝑃 𝑒 1 𝑁 𝑡 − 1  𝑖 = 0 | | 𝐻 𝑘 , 𝑖 | | 2              𝛼 ℎ . ( 3 5 ) Therefore, the symbol probability of error is given at second iteration by 𝑃 𝑒 2 = 𝐶 1 ⎛ ⎜ ⎝ 1 𝑁 1 0 ⎞ ⎟ ⎠ = 𝐶 1  1 𝑁 0 + 𝑃 𝑒 1 𝛼 ℎ  = 𝐶 1 ⎛ ⎜ ⎝ 1 𝑁 0  1 + 𝛼 ℎ 𝐶 2  1 / 𝑁 0  / 𝑁 0  ⎞ ⎟ ⎠ , ( 3 6 ) where 1 / 𝑁 1 0 is the SNR at the input of Decoder 2. 𝐶 2 ( 1 / 𝑁 0 ) is a 𝑄 -function that is exponentially decreasing as SNR increases; thus, 𝛼 ℎ 𝐶 2 ( 1 / 𝑁 0 ) / 𝑁 0 decreases as SNR increases since the exponential function overwhelms the polynomial function. Then, there is a noise power 𝑁 𝑎 0 such that, for 𝑁 0 < 𝑁 𝑎 0 𝛼 ℎ 𝐶 2  1 / 𝑁 0  𝑁 0 < 𝛼 Δ , ( 3 7 ) and thus, 1 𝑁 0  1 + 𝛼 ℎ 𝐶 2  1 / 𝑁 0  / 𝑁 0  > 1  1 + 𝛼 Δ  𝑁 0 . ( 3 8 ) Therefore for 𝑁 0 < 𝑁 𝑎 0 , 𝐶 1 ⎛ ⎜ ⎝ 1 𝑁 0  1 + 𝛼 ℎ 𝐶 2  1 / 𝑁 0  / 𝑁 0  ⎞ ⎟ ⎠ < 𝐶 1 ⎛ ⎜ ⎝ 1  1 + 𝛼 Δ  𝑁 0 ⎞ ⎟ ⎠ = 𝐶 2  1 𝑁 0  , ( 3 9 ) that is, 𝑃 𝑒 2 < 𝑃 𝑒 1 . ( 4 0 ) For 𝑁 0 < 𝑁 𝑎 0 the output of the second iteration will give better performance than that of the first iteration. This second iteration is summarized by the point 𝐴 2 ( 1 / 𝑁 1 0 , 𝑃 𝑒 2 ) in Figure 4 . When recombining the signal at the input of Decoder 1 for the third iteration using ( 29 ), the noise component is now smaller than that in the previous case since 𝑃 𝑒 2 < 𝑃 𝑒 1 . Consequently, the third iteration performance is given by 𝐶 2 at S N R = 1 / 𝑁 2 0 with 𝑁 2 0 < 𝑁 1 0 . Thus, 𝐶 2 ( 1 / 𝑁 2 0 ) < 𝐶 2 ( 1 / 𝑁 0 ) , that is, the probability of error at the output of Decoder 1 for the third iteration 𝑃 𝑒 3 is less than that for 𝑃 𝑒 1 . This third iteration is summarized by the point 𝐴 3 ( 1 / 𝑁 2 0 , 𝑃 𝑒 3 ) in Figure 4 . Let us notice that 𝑃 𝑒 3 could be greater than 𝑃 𝑒 2 . The next iteration performance can be derived in the same manner since we just have to replace 𝑁 0 by 𝑁 2 0 . Thus, the probability of error at the output of a given decoder (Decoder 1 or Decoder 2) will always decrease or reach a fixed point. 6. Simulation Results In this section, we will evaluate the performance of the two decoding methods that we have presented. We consider a transmission scheme with two and three transmit antennas. For 𝑁 𝑡 = 2 , we have 𝑡 ( 1 ) 𝑘 , 𝑛 + 2 𝑁 𝑡 − 2 = | | 𝐻 𝑘 , 1 | | 2 𝑎 𝑘 , 𝑛 − 2 + 2 ℜ  𝐻 ∗ 𝑘 , 1 𝐻 𝑘 , 0  𝑎 𝑘 , 𝑛 + | | 𝐻 𝑘 , 0 | | 2 𝑎 𝑘 , 𝑛 + 2 + 𝑤 𝑘 , 𝑛 + 2 𝑁 𝑡 − 2 , ( 4 1 ) and for 𝑁 𝑡 = 3 , we get 𝑡 ( 1 ) 𝑘 , 𝑛 + 2 𝑁 𝑡 − 2 = | | 𝐻 𝑘 , 2 | | 2 𝑎 𝑘 , 𝑛 − 4 + 2 ℜ  𝐻 ∗ 𝑘 , 2 𝐻 𝑘 , 1  𝑎 𝑘 , 𝑛 − 2 +  2 ℜ  𝐻 ∗ 𝑘 , 2 𝐻 𝑘 , 0  + | | 𝐻 𝑘 , 1 | | 2  𝑎 𝑘 , 𝑛 + 2 ℜ  𝐻 ∗ 𝑘 , 1 𝐻 𝑘 , 0  𝑎 𝑘 , 𝑛 + 2 + | | 𝐻 𝑘 , 0 | | 2 𝑎 𝑘 , 𝑛 + 4 . ( 4 2 ) The simulation parameters we consider are given as follows: (i) no channel coding, (ii) QPSK modulation, (iii) Rayleigh channel per antenna, that is, flat over all the subcarriers. We assume that the channel coefficients are perfectly known by the receiver, (iv) number of subcarrier 𝑀 = 3 2 , (v) we used a truncation of the IOTA (Isotropic Orthogonal Transform Algorithm) prototype function [ 1 ]. Its duration is limited to 4 𝑇 0 , which leads to a nearly orthogonal prototype filter containing 𝐿 = 4 𝑀 = 1 2 8 taps. In this section, we give BER (Bit Error Rate) versus SNR simulation results, and consequently, we do not take into consideration the loss of efficiency due to the cyclic prefix in CP-OFDM modulation. In Figure 5 we show the performance of the FBMC decoding structure introduced in Figure 2 . For FBMC, we consider both ML and Viterbi decoding. ML decoding using an exhaustive search among all possible transmitted sequences of data outperforms Viterbi decoding by 1 dB. This is due to the fact that the noise is colored; thus, Viterbi decoding is suboptimal. We also give the CP-OFDM performance using a Viterbi decoding. We can see that CP-OFDM outperforms ML/FBMC by about 1 dB. Figure 5: Performance of single delay STTC with 2 transmit antennas and one receive antenna (FBMC and CP-OFDM modulation). In the rest of this section, we will focus on the iterative decoding performance. The simulation results are obtained using Viterbi decoding blocks implemented inside Decoder 1 and Decoder 2 blocks in Figure 3 . The Viterbi algorithm implemented in Decoder 1 is related to ( 41 ). For QPSK modulation, the Trellis is a 4 𝑁 𝑇 − 1 state Trellis with only two possible transitions per state since the detection is performed on real data. Whereas the Viterbi algorithm implemented in Decoder 2 is related to ( 28 ) and is a 2 𝑁 𝑇 − 1 state Trellis with two transitions per state, again detection is performed on real data. We also consider hard estimation of the data at the output of a given Viterbi decoder. For the CP-OFDM case with QPSK modulation, we have a 4 𝑁 𝑇 − 1 state Trellis with 4 transitions per state as the detection is performed on complex data. Therefore, this Viterbi algorithm is more complex compared with one of the two Viterbi algorithms used in the case of FBMC modulation. The two Viterbi algorithms used in FBMC taken together have a complexity comparable to the one used with CP-OFDM. However, the two Viterbi algorithms used in FBMC operate on a frame sequence which is two times longer than the one for CP-OFDM modulation. Then, in terms of complexity the proposed FBMC structure has a significantly higher complexity than that of CP-OFDM mainly due to the “Interference estimation + Interference cancelation” block. For uncorrelated Rayleigh channels, we plot the performance of this FBMC receiver structure for different iteration stages as well as the performance of CP-OFDM with ML decoding as a matter of comparison. Figures 6 and 7 provide the simulation results for 𝑁 𝑡 = 2 and 𝑁 𝑡 = 3 , respectively. For 𝑛 = 1 , we have a 2 dB degradation compared to CP-OFDM. For 𝑛 ≥ 2 (more than two-Viterbi decoding), we get closer to CP-OFDM. For 𝑛 = 5 or 6 , we almost reach the same performance as that of CP-OFDM. In Figure 6 we also plot the curve obtained when we assume perfect interference cancelation in the second iteration as mentioned in ( 28 ). In that case, there is a possible gain of 0 . 8 dB since the Viterbi structure with 2 states and two transitions per state (Decoder 2) provides better performance than the 4-state Viterbi decoder with 4 transitions per state implemented for CP-OFDM. Indeed, it is possible to show that the structures of the code related to these two Trellises have the same minimum distance. However, the performance gain is due to the distance distribution associated to the two Trellises. Figure 6: Performance of single delay STTC (iterative decoding) with 2 transmit antennas and one receive antenna (FBMC and CP-OFDM modulation). Figure 7: Performance of single delay STTC (iterative decoding) with 3 transmit antennas and one receive antenna (FBMC and CP-OFDM modulation). Moreover, let us evaluate this scheme in presence of a frequency selective channel. We consider the following channel parameters: (i) uncoded QPSK modulation, (ii) 𝑀 = 6 4 subcarriers, (iii) static channels (no Doppler), IOTA prototype, (iv) 3-tap channels between the transmit antennas and the receive: power profile: 0, −4, −10 (dB) Delay: 0, 1, 2 (number of samples), (v) OFDM Cyclic Prefix length: 4 samples, (vi) perfect channel estimation. As shown in Figure 8 , after one iteration ( 𝑛 = 1 ), we have about 2 dB degradation compared to CP-OFDM. For 𝑛 = 2 and 𝑛 = 3 , the loss is reduced to 0 . 7 dB, and for 𝑛 = 6 the degradation is about 0 . 3 dB compared to CP-OFDM. However, the iterative method has an inherent gain as FBMC does not use a CP contrary to CP-OFDM. Figure 8: Performance of single delay STTC (iterative decoding) with 2 transmit antennas and one receive antenna over frequency selective channels. 7. Conclusion In this paper, we have presented two general methods for data detection when combining FBMC and single delay STTC as well as the interference cancelation and the iterative methods. The interference cancelation method despite its simplicity has poorer performance compared to that of CP-OFDM. Thus, we have proposed an iterative decoding based on interference estimation and cancelation which does not require any channel coding or decoding block. We have shown that in the case of QPSK modulation and Rayleigh or frequency selective channels it is possible with this decoding method to perform as better as OFDM-STTC. Moreover if the iterative cancelation process is improved, then a potential gain can be achieved. This is obtained with a relatively higher complexity. In future work, we will look at FBMC with other STTC schemes and evaluate their performance under nonlocally flat channels. Appendix A. General Expression of 𝐴 𝑘 , 𝑛 , 𝐵 𝑘 , 𝑛 , and 𝐶 𝑘 , 𝑛 A.1. 𝐴 𝑘 , 𝑛 Value Let us compute 𝐴 𝑘 , 𝑛 . ( i ) Case 𝑁 𝑡 is even, that is, 𝑁 𝑡 = 2 𝑈 𝑡 , 𝐴 𝑘 , 𝑛 = 𝑥 𝑘 , 𝑛 ⎛ ⎜ ⎝ 𝑈 𝑡 − 1  𝑖 = 0 𝐻 ∗ 𝑘 , 2 𝑈 𝑡 − 1 − 𝑖 𝐻 𝑘 , 𝑖 + 2 𝑈 𝑡 − 1  𝑖 = 𝑈 𝑡 𝐻 ∗ 𝑘 , 2 𝑈 𝑡 − 1 − 𝑖 𝐻 𝑘 , 𝑖 ⎞ ⎟ ⎠ . ( A . 1 ) Using the relation 𝑞 = 2 𝑈 𝑡 − 1 − 𝑖 , we have 𝐴 𝑘 , 𝑛 = 𝑥 𝑘 , 𝑛 ⎛ ⎜ ⎝ 𝑈 𝑡 − 1  𝑖 = 0 𝐻 ∗ 𝑘 , 2 𝑈 𝑡 − 1 − 𝑖 𝐻 𝑘 , 𝑖 + 𝑈 𝑡 − 1  𝑞 = 0 𝐻 ∗ 𝑘 , 𝑞 𝐻 𝑘 , 2 𝑈 𝑡 − 1 − 𝑞 ⎞ ⎟ ⎠ = 𝑥 𝑘 , 𝑛 ⎛ ⎜ ⎝ 𝑈 𝑡 − 1  𝑖 = 0  𝐻 ∗ 𝑘 , 2 𝑈 𝑡 − 1 − 𝑖 𝐻 𝑘 , 𝑖 + 𝐻 ∗ 𝑘 , 𝑖 𝐻 𝑘 , 2 𝑈 𝑡 − 1 − 𝑖  ⎞ ⎟ ⎠ = 𝑥 𝑘 , 𝑛 ⎛ ⎜ ⎝ 2 𝑈 𝑡 − 1  𝑖 = 0 ℜ   𝐻 ∗ 𝑘 , 2 𝑈 𝑡 − 1 − 𝑖 𝐻 𝑘 , 𝑖   ⎞ ⎟ ⎠                                                            𝜇 𝑘 . ( A . 2 ) ( i i ) Case 𝑁 𝑡 is odd, that is, 𝑁 𝑡 = 2 𝑈 𝑡 + 1 , 𝜇 𝑘 = ⎛ ⎜ ⎝ 𝑈 𝑡 − 1  𝑖 = 0 𝐻 ∗ 𝑘 , 2 𝑈 𝑡 − 𝑖 𝐻 𝑘 , 𝑖 + 𝐻 ∗ 𝑘 , 𝑈 𝑡 𝐻 𝑘 , 𝑈 𝑡 + 2 𝑈 𝑡  𝑖 = 𝑈 𝑡 + 1 𝐻 ∗ 𝑘 , 2 𝑈 𝑡 − 𝑖 𝐻 𝑘 , 𝑖 ⎞ ⎟ ⎠ . ( A . 3 ) Again using 𝑞 = 2 𝑈 𝑡 − 𝑖 , we have 𝜇 𝑘 = ⎛ ⎜ ⎝ 2 𝑈 𝑡 − 1  𝑖 = 0 ℜ   𝐻 ∗ 𝑘 , 2 𝑈 𝑡 − 𝑖 𝐻 𝑘 , 𝑖   + | | 𝐻 𝑘 , 𝑈 𝑡 | | 2 ⎞ ⎟ ⎠ , ( A . 4 ) and we get 𝐴 𝑘 , 𝑛 = 𝑥 𝑘 , 𝑛 𝜇 𝑘 . ( A . 5 ) A.2. 𝐵 𝑘 , 𝑛 Value Let us now compute 𝐵 𝑘 , 𝑛 ; setting 𝑞 = 𝑝 − 𝑖 , we get 𝐵 𝑘 , 𝑛 = 𝑁 𝑡 − 1  𝑖 = 1 𝑖 − 1  𝑝 = 0 𝐻 ∗ 𝑘 , 𝑁 𝑡 − 1 − 𝑝 𝐻 𝑘 , 𝑖 𝑥 𝑘 , 𝑛 + 2 𝑝 − 2 𝑖 = 𝑁 𝑡 − 1  𝑖 = 1 𝑖  𝑞 = 1 𝑥 𝑘 , 𝑛 − 2 𝑞 𝐻 ∗ 𝑘 , 𝑁 𝑡 − 1 + 𝑞 − 𝑖 𝐻 𝑘 , 𝑖 . ( A . 6 ) This last equation is the sum over a triangular set of index; therefore, the sum can be taken either from lines or from columns where the total is the same. Therefore, 𝐵 𝑘 , 𝑛 = 𝑁 𝑡 − 1  𝑞 = 1 𝑥 𝑘 , 𝑛 − 2 𝑞 𝑁 𝑡 − 1  𝑖 = 𝑞 𝐻 ∗ 𝑘 , 𝑁 𝑡 − 1 + 𝑞 − 𝑖 𝐻 𝑘 , 𝑖 . ( A . 7 ) Taking 𝑚 = 𝑖 − 𝑞 , we get 𝐵 𝑘 , 𝑛 = 𝑁 𝑡 − 1  𝑞 = 1 𝑥 𝑘 , 𝑛 − 2 𝑞 𝑁 𝑡 − 1 − 𝑞  𝑚 = 0 𝐻 ∗ 𝑘 , 𝑁 𝑡 − 1 − 𝑚 𝐻 𝑘 , 𝑚 + 𝑞                                          𝛾 𝑞 . ( A . 8 ) ( i ) Case 𝑁 𝑡 − 𝑞 is even, that is, 𝑁 𝑡 − 𝑞 = 2 𝑈 𝑞 ; then, 𝛾 𝑞 = 𝑁 𝑡 − 1 − 𝑞  𝑚 = 0 𝐻 ∗ 𝑘 , 𝑁 𝑡 − 1 − 𝑚 𝐻 𝑘 , 𝑚 + 𝑞 = 𝑈 𝑞 − 1  𝑚 = 0 𝐻 ∗ 𝑘 , 2 𝑈 𝑞 + 𝑞 − 1 − 𝑚 𝐻 𝑘 , 𝑚 + 𝑞 + 2 𝑈 𝑞 − 1  𝑚 = 𝑈 𝑞 𝐻 ∗ 𝑘 , 2 𝑈 𝑞 + 𝑞 − 1 − 𝑚 𝐻 𝑘 , 𝑚 + 𝑞 = 2 𝑈 𝑞 − 1  𝑚 = 0 ℜ  𝐻 ∗ 𝑘 , 2 𝑈 𝑞 + 𝑞 − 1 − 𝑚 𝐻 𝑘 , 𝑚 + 𝑞  . ( A . 9 ) ( i i ) Case 𝑁 𝑡 − 𝑞 is odd, that is, 𝑁 𝑡 − 𝑞 = 2 𝑈 𝑞 + 1 ; then, 𝛾 𝑞 = 𝑁 𝑡 − 1 − 𝑞  𝑚 = 0 𝐻 ∗ 𝑘 , 𝑁 𝑡 − 1 − 𝑚 𝐻 𝑘 , 𝑚 + 𝑞 = 𝑈 𝑞 − 1  𝑚 = 0 𝐻 ∗ 𝑘 , 2 𝑈 𝑞 + 𝑞 − 𝑚 𝐻 𝑘 , 𝑚 + 𝑞 + 𝐻 ∗ 𝑘 , 𝑈 𝑞 + 𝑞 𝐻 𝑘 , 𝑈 𝑞 + 𝑞 + 2 𝑈 𝑞  𝑚 = 𝑈 𝑞 + 1 𝐻 ∗ 𝑘 , 2 𝑈 𝑞 + 𝑞 − 𝑚 𝐻 𝑘 , 𝑚 + 𝑞 = 2 𝑈 𝑞 − 1  𝑚 = 0 ℜ  𝐻 ∗ 𝑘 , 2 𝑈 𝑞 + 𝑞 − 𝑚 𝐻 𝑘 , 𝑚 + 𝑞  + | | | 𝐻 ∗ 𝑘 , 𝑈 𝑞 + 𝑞 | | | 2 . ( A . 1 0 ) A.3. 𝐶 𝑘 , 𝑛 Value Let us now compute 𝐶 𝑘 , 𝑛 ; setting 𝑞 = 𝑝 − 𝑖 , we get 𝐶 𝑘 , 𝑛 = 𝑁 𝑡 − 2  𝑖 = 1 𝑁 𝑡 − 1 − 𝑖  𝑞 = 1 𝐻 ∗ 𝑘 , 𝑁 𝑡 − 1 − 𝑞 − 𝑖 𝐻 𝑘 , 𝑖 𝑥 𝑘 , 𝑛 + 2 𝑞 . ( A . 1 1 ) This last equation is the sum over a triangular set of index; therefore, the sum can be taken either from lines or from columns where the total is the same. Therefore, 𝐶 𝑘 , 𝑛 = 𝑁 𝑡 − 1  𝑞 = 1 𝑥 𝑘 , 𝑛 + 2 𝑞 𝑁 𝑡 − 1 − 𝑞  𝑖 = 0 𝐻 ∗ 𝑘 , 𝑁 𝑡 − 1 − 𝑞 − 𝑖 𝐻 𝑘 , 𝑖                                        𝛽 𝑞 . ( A . 1 2 ) ( i ) Case 𝑁 𝑡 − 𝑞 is even, that is, 𝑁 𝑡 − 𝑞 = 2 𝑈 𝑞 ; then, 𝛽 𝑞 = 𝑁 𝑡 − 1 − 𝑞  𝑖 = 0 𝐻 ∗ 𝑘 , 𝑁 𝑡 − 1 − 𝑞 − 𝑖 𝐻 𝑘 , 𝑖 = 𝑈 𝑞 − 1  𝑖 = 0 𝐻 ∗ 𝑘 , 2 𝑈 𝑞 − 1 − 𝑖 𝐻 𝑘 , 𝑖 + 2 𝑈 𝑞 − 1  𝑚 = 𝑈 𝑞 𝐻 ∗ 𝑘 , 2 𝑈 𝑞 − 1 − 𝑖 𝐻 𝑘 , 𝑖 = 2 𝑈 𝑞 − 1  𝑖 = 0 ℜ  𝐻 ∗ 𝑘 , 2 𝑈 𝑞 − 1 − 𝑖 𝐻 𝑘 , 𝑖  . ( A . 1 3 ) ( i i ) Case 𝑁 𝑡 − 𝑞 is odd, that is, 𝑁 𝑡 − 𝑞 = 2 𝑈 𝑞 + 1 ; then, 𝛽 𝑞 = 𝑁 𝑡 − 1 − 𝑞  𝑖 = 0 𝐻 ∗ 𝑘 , 𝑁 𝑡 − 1 − 𝑞 − 𝑖 𝐻 𝑘 , 𝑖 = 𝑈 𝑞 − 1  𝑖 = 0 𝐻 ∗ 𝑘 , 2 𝑈 𝑞 − 𝑖 𝐻 𝑘 , 𝑖 + 𝐻 ∗ 𝑘 , 𝑈 𝑞 𝐻 𝑘 , 𝑈 𝑞 + 2 𝑈 𝑞  𝑚 = 𝑈 𝑞 + 1 𝐻 ∗ 𝑘 , 2 𝑈 𝑞 − 𝑖 𝐻 𝑘 , 𝑖 = 2 𝑈 𝑞 − 1  𝑚 = 0 ℜ  𝐻 ∗ 𝑘 , 2 𝑈 𝑞 − 𝑖 𝐻 𝑘 , 𝑖  + | | | 𝐻 ∗ 𝑘 , 𝑈 𝑞 | | | 2 . ( A . 1 4 ) Acknowledgments The authors would like to thank Pr. M. Bellanger for helpful discussions. This work was supported in part by the European Commission under Project PHYDYAS (FP7-ICT-2007-1-211887). <h4>References</h4> B. Le Floch, M. Alard, and C. Berrou, “Coded orthogonal frequency division multiplex,” Proceedings of the IEEE , vol. 83, pp. 982–996, 1995. H. Boelcskei, “Orthogonal frequency division multiplexing based on offset QAM,” in Advances in Gabor Analysis , Birkhäuser, Boston, Mass, USA, 2003. P. Siohan, C. Siclet, and N. Lacaille, “ Analysis and design of OFDM/OQAM systems based on filterbank theory ,” IEEE Transactions on Signal Processing , vol. 50, no. 5, pp. 1170–1183, 2002. R. W. Chang, “Synthesis of band-limited orthogonal signals for multi-channel data transmission,” Bell Labs Technical Journal , vol. 45, pp. 1775–1796, 1966. B. R. Saltzberg, “Performance of an efficient parallel data transmission system,” IEEE Transactions on Communication Technology , vol. 15, no. 6, pp. 805–811, 1967. V. Tarokh, N. Seshadri, and A. R. Calderbank, “Space-time codes for high data rate wireless communication: performance criterion and code construction,” IEEE Transactions on Information Theory , vol. 44, no. 2, pp. 744–765, 1998. V. Tarokh, H. Jafarkhani, and A. R. Calderbank, “ Space-time block coding for wireless communications: performance results ,” IEEE Journal on Selected Areas in Communications , vol. 17, no. 3, pp. 451–460, 1999. M. Bellanger, “ Transmit diversity in multicarrier transmission using OQAM modulation ,” in Proceedings of the 3rd International Symposium on Wireless Pervasive Computing (ISWPC '08) , pp. 727–730, Santorini, Greece, May 2008. B. Hirosaki, “Orthogonally multiplexed QAM system using the discrete Fourier transform,” IEEE Transactions on Communications Systems , vol. 29, no. 7, pp. 982–989, 1981. C. Lélé, P. Siohan, R. Legouable, and J.-P. Javaudin, “ Preamble-based channel estimation techniques for OFDM/OQAM over the powerline ,” in Proceedings of IEEE International Symposium on Power Line Communications and Its Applications (ISPLC '07) , pp. 59–64, Pisa, Italy, March 2007. G. D. Forney Jr., “Maximum-likelihood sequence estimation of digital sequences in the presence of intersymbol interference,” IEEE Transactions on Information Theory , vol. 18, no. 3, pp. 363–378, 1972. J.-D. Wang and H. Y. Chung, “ Trellis coded communication systems in the presence of colored noise: performance analysis, simulation, and the swapping technique ,” in Proceeingds of IEEE Conference Record on Global Telecommunications Conference, and Exhibition. ‘Communications for the Information Age’ (GLOBECOM '88) , vol. 2, pp. 1160–1165, Hollywood, Fla, USA, November-December 1988. // http://www.deepdyve.com/assets/images/DeepDyve-Logo-lg.png EURASIP Journal on Advances in Signal Processing Hindawi Publishing Corporation

Decoding Schemes for FBMC with Single-Delay STTC

Loading next page...
 
/lp/hindawi-publishing-corporation/decoding-schemes-for-fbmc-with-single-delay-sttc-axyKC6H7l1

References

References for this paper are not available at this time. We will be adding them shortly, thank you for your patience.

Publisher
Hindawi Publishing Corporation
Copyright
Copyright © 2010 Chrislin Lélé and Didier Le Ruyet.
ISSN
1687-6172
eISSN
1687-6180
Publisher site
See Article on Publisher Site

Abstract

Decoding Schemes for FBMC with Single-Delay STTC //// Hindawi Publishing Corporation Home Journals About Us About this Journal Submit a Manuscript Table of Contents Journal Menu Abstracting and Indexing Aims and Scope Article Processing Charges Articles in Press Author Guidelines Best Paper Awards Bibliographic Information Contact Information Editorial Board Editorial Workflow Free eTOC Alerts Reviewers Acknowledgment Society Affiliation Subscription Information Open Special Issues Published Special Issues Special Issue Guidelines Abstract Full-Text PDF Full-Text HTML Linked References How to Cite this Article Complete Special Issue EURASIP Journal on Advances in Signal Processing Volume 2010 (2010), Article ID 689824, 11 pages doi:10.1155/2010/689824 Research Article <h2>Decoding Schemes for FBMC with Single-Delay STTC</h2> Chrislin Lélé and Didier Le Ruyet (EURASIP Member) Electronics and Communications Laboratory, Conservatoire National Des Arts Et Métiers (CNAM), 75141 Paris, France Received 5 June 2009; Accepted 28 December 2009 Academic Editor: Markku Renfors Copyright © 2010 Chrislin Lélé and Didier Le Ruyet. This is an open access article distributed under the Creative Commons Attribution License , which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. Abstract Orthogonally multiplexed Quadrature Amplitude Modulation (OQAM) with Filter-Bank-based MultiCarrier modulation (FBMC) is a multicarrier modulation scheme that can be considered an alternative to the conventional orthogonal frequency division multiplexing (OFDM) with cyclic prefix (CP) for transmission over multipath fading channels. However, as OQAM-based FBMC is based on real orthogonality, transmission over a complex-valued channel makes the decoding process more challenging compared to CP-OFDM case. Moreover, if we apply Multiple Input Multiple Output (MIMO) techniques to OQAM-based FBMC, the decoding schemes are different from the ones used in CP-OFDM. In this paper, we consider the combination of OQAM-based FBMC with single-delay Space-Time Trellis Coding (STTC). We extend the decoding process presented earlier in the case of 𝑁 𝑡 = 2 transmit antennas to greater values of 𝑁 𝑡 . Then, for 𝑁 𝑡 ≥ 2 , we make an analysis of the theoretical and simulation performance of ML and Viterbi decoding. Finally, to improve the performance of this method, we suggest an iterative decoding method. We show that the OQAM-based FBMC iterative decoding scheme can slightly outperform CP-OFDM. 1. Introduction Orthogonal Frequency Division Multiplexing (OFDM) is an efficient Multicarrier Modulation (MCM) capable of fighting against multipath fading channels. Its robustness to multi-path propagation effects comes from the insertion of a CP and is, therefore, obtained at the price of a reduced spectral efficiency. Furthermore, the rectangular shape of OFDM symbols leads to a s i n ( 𝑥 ) / 𝑥 frequency spectrum. Studies have been conducted in order to find better MCM schemes with respect to the frequency and/or time-frequency localization criteria. As suggested in [ 1 – 3 ], OFDM/OQAM also called as OQAM-based Filter Bank Multicarrier (FBMC) is an MCM scheme which may be the appropriate alternative. In OFDM/OQAM each subcarrier is modulated with Offset Quadrature Amplitude Modulation (OQAM). This principle has been introduced in [ 4 , 5 ], but it is only recently [ 1 ] that FBMC has been presented as a viable alternative to OFDM. Compared to OFDM that transmits complex-valued symbols at a given symbol rate, OQAM-based FBMC transmits real-valued symbols at twice this symbol rate. Therefore, a similar spectral efficiency is achieved by both systems. In practice, OQAM-based FBMC may provide a higher useful bit rate since it operates without the addition of a CP. Furthermore, with a pulse shaping that can be optimized according to given channel characteristics, its performance can be improved. However, all the interesting features of OQAM-based FBMC come at the price of a relaxation of the orthogonality conditions that only hold in the real field. At the receive side the data is carried only by the real component of the signal (assuming a 0 or 𝜋 / 2 phase modulation term). Thus, the imaginary part appears as an interference term. This interference term is a source of problem in the presence of the complex-valued channel as it destroys the real orthogonality. Therefore, when combining OQAM-based FBMC with MIMO technique such as Space-Time Block Codes (STBC) or Space-Time Trellis Coding (STTC) [ 6 , 7 ], the decoding process cannot be done in the same way as with CP-OFDM modulation. In the case of a single delay STTC chain with 2 transmit and 1 receive antennas, refrence [ 8 ] proposed a simple preprocessing to cancel this imaginary interference component. In this paper, we extend the proposed method in [ 8 ] to 𝑁 𝑡 transmit antennas and introduce an iterative decoding method. In Section 2 , we give a short description of the discrete-time OQAM modulation. Then, in Section 3 , we provide an overview of the STTC single delay detection. In Section 4.1 , we provide a theoretical performance analysis of ML and Viterbi decoding. Section 5 is devoted to the iterative decoding method in order to improve the performance of the previous decoding method. Simulation results are presented in Section 6 . Conclusions and perspectives are given in Section 7 . In the rest of the paper, FBMC will be used to denote OQAM-based FBMC. 2. The FBMC Modulation Using the baseband discrete-time model, we can write at the transmit antenna 𝑖 the OQAM-based FBMC signal as follows [ 1 ]: 𝑠 𝑖 [ 𝑚 ] = 𝑀 − 1  𝑘 = 0  𝑛 ∈ 𝐙 𝑑 𝑘 , 𝑛 , 𝑖 𝑔 [ 𝑚 − 𝑛 𝑁 ] 𝑒 𝑗 ( 2 𝜋 / 𝑀 ) 𝑘 ( 𝑚 − 𝐷 / 2 ) 𝑒 𝑗 𝜙 𝑘 , 𝑛                                                          𝑔 𝑘 , 𝑛 [ 𝑚 ] , ( 1 ) where 𝑀 = 2 𝑁 is the even number of subcarriers, 𝐹 0 = 1 / 𝑇 0 = 1 / 2 𝜏 0 is the subcarrier spacing, 𝜙 𝑘 , 𝑛 is an additional phase term, 𝑔 is the pulse shape, and 𝐷 is the delay parameter associated to the length of the pulse shape. The transmitted symbols 𝑑 𝑘 , 𝑛 , 𝑖 are real-valued data transmitted by antenna 𝑖 . They are obtained from a 2 2 𝐾 -QAM constellation, taking the real and imaginary parts of these complex-valued symbols of duration 𝑇 0 = 2 𝜏 0 , where 𝜏 0 denotes the time offset between the two parts [ 1 – 3 , 9 ]. For a given subcarrier 𝑘 and symbol time index 𝑛 , the real and imaginary parts are driven by the phase term 𝜙 𝑘 , 𝑛 given by 𝜙 𝑘 , 𝑛 = 𝜙 0 + 𝜋 2 ( 𝑛 + 𝑘 ) ( m o d 𝜋 ) , ( 2 ) where 𝜙 0 can be arbitrarily chosen. Here, we set 𝜙 0 = 0 and 𝑔 is assumed to be real valued. Assuming a distortion-free channel, a perfect reconstruction of real symbols is obtained owing to the following real orthogonality condition: ℜ   𝑔 𝑘 , 𝑛 ∣ 𝑔 𝑝 , 𝑞   = ℜ  ∞  𝑚 = − ∞ 𝑔 𝑘 , 𝑛 [ 𝑚 ] 𝑔 ∗ 𝑝 , 𝑞 [ 𝑚 ]  = 𝛿 𝑘 , 𝑝 𝛿 𝑛 , 𝑞 , ( 3 ) where 𝛿 𝑘 , 𝑝 = 1 if 𝑘 = 𝑝 and 𝛿 𝑛 , 𝑞 = 0 if 𝑛 ≠ 𝑞 . However, in practice for transmission over a realistic channel, the orthogonality property is lost, leading to intersymbol and intercarrier interferences. It has been shown in previous studies [ 8 ] that, when combining FBMC with single delay STTC in presence of 2 transmit and one receive antennas, specific processing should be done in order to remove the interference terms. In this paper, we will extend this method for 𝑁 𝑡 ≥ 2 antennas. 3. Single-Delay STTC in FBMC with 𝑁 𝑡 Transmit Antennas 3.1. Transmission Model Let us first assume that only the 𝑖 th antenna is transmitting. At the receiver side, the demodulated signal 𝑦 𝑘 , 𝑛 at the frequency 𝑘 and time instant 𝑛 ( 𝑛 𝜏 0 ) can be written as 𝑦 𝑘 , 𝑛 = 𝐻 𝑘 , 𝑛 , 𝑖 𝑑 𝑘 , 𝑛 , 𝑖 + 𝑗 𝐼 𝑘 , 𝑛 , 𝑖 + 𝜐 𝑘 , 𝑛 , ( 4 ) where (i) 𝐻 𝑘 , 𝑛 , 𝑖 is the channel coefficient between transmit antenna 𝑖 and the receiver, at subcarrier 𝑘 and time instant 𝑛 , (ii) 𝜐 𝑘 , 𝑛 is the noise component at subcarrier 𝑘 and time instant 𝑛 , (iii) 𝐼 𝑘 , 𝑛 , 𝑖 = ( − 𝑗 )  ( 𝑘  , 𝑛  ) ≠ ( 𝑘 , 𝑛 ) 𝐻 𝑘  , 𝑛  , 𝑖 𝑑 𝑘  , 𝑛  , 𝑖 ∞  𝑚 = − ∞ 𝑔 𝑘 , 𝑛 [ 𝑚 ] 𝑔 ∗ 𝑘  , 𝑛  [ 𝑚 ] . ( 5 ) We assume that we have a prototype filter well localized in time and frequency. This implies that in the previous equation the main contribution comes from the closest neighborhood, that is, 𝑔 𝑘 , 𝑛 [ 𝑚 ] 𝑔 ∗ 𝑘  , 𝑛  [ 𝑚 ] takes a significant value only for | 𝑘 − 𝑘 ′ | ≤ 1 and | 𝑛 − 𝑛 ′ | ≤ 1 . Moreover, if we assume that the channel is constant over a set of at least three consecutive subcarriers and a set of at least three consecutive time indexes, then we can rewrite the previous expression as in [ 10 ]: 𝐼 𝑘 , 𝑛 , 𝑖 ≈ 𝐻 𝑘 , 𝑛 , 𝑖 ( − 𝑗 )  ( 𝑘  , 𝑛  ) ≠ ( 𝑘 , 𝑛 ) 𝑑 𝑘  , 𝑛  , 𝑖 ∞  𝑚 = − ∞ 𝑔 𝑘 , 𝑛 [ 𝑚 ] 𝑔 ∗ 𝑘  , 𝑛  [ 𝑚 ]                                                                                  𝑢 𝑘 , 𝑛 , 𝑖 . ( 6 ) Thus, the demodulated signal can be approximated by 𝑦 𝑘 , 𝑛 ≈ 𝐻 𝑘 , 𝑛 , 𝑖  𝑑 𝑘 , 𝑛 , 𝑖 + 𝑗 𝑢 𝑘 , 𝑛 , 𝑖  + 𝜐 𝑘 , 𝑛 . ( 7 ) Throughout the remainder of the paper, we will consider ( 7 ) as the expression of the signal at the output of the demodulator. 3.2. Problem Statement Let us consider the single delay STTC scheme with 𝑁 𝑡 antennas as shown in Figure 1 . The real data to be transmitted is modulated by an FBMC modulator and transmitted by the first antenna. The same stream of data is delayed by 2 𝑛 𝑖 real data before being modulated by FBMC modulator and transmitted by the 𝑛 𝑖 th antenna. The delay 2 𝑛 𝑖 is chosen to have the same delay as with a CP-OFDM system although a delay of 𝑛 𝑖 could also be chosen. We denote by 𝑎 𝑘 , 𝑛 the real data from the main stream of data at frequency 𝑘 and time index 𝑛 . Thus, at a given subcarrier 𝑘 the transmission is given at antenna 𝑖 by 𝑑 𝑘 , 𝑛 , 𝑖 = 𝑎 𝑘 , 𝑛 − 2 𝑖 . At the receiver side, the demodulated signal can be written as Figure 1: FBMC Single-delay STTC transmitter. 𝑦 𝑘 , 𝑛 = 𝑁 𝑡 − 1  𝑖 = 0 𝐻 𝑘 , 𝑛 , 𝑖  𝑑 𝑘 , 𝑛 , 𝑖 + 𝑗 𝑢 𝑘 , 𝑛 , 𝑖  + 𝜐 𝑘 , 𝑛 , ( 8 ) where 𝜐 𝑘 , 𝑛 is the noise component at the subcarrier 𝑘 and time instant 𝑛 . As the same stream of data is transmitted over the 𝑁 𝑡 antennas, we have 𝑢 𝑘 , 𝑛 , 𝑖 = 𝑢 𝑘 , 𝑛 − 2 𝑖 , 0 = 𝑏 𝑘 , 𝑛 − 2 𝑖 . In the remainder of the paper, we will assume a channel constant over time, that is, ( 𝐻 𝑘 , 𝑛 , 𝑖 = 𝐻 𝑘 , 𝑖 ); we get 𝑦 𝑘 , 𝑛 = 𝑁 𝑡 − 1  𝑖 = 0 𝐻 𝑘 , 𝑖 ( 𝑎 𝑘 , 𝑛 − 2 𝑖 + 𝑗 𝑏 𝑘 , 𝑛 − 2 𝑖 )                              𝑥 𝑘 , 𝑛 − 2 𝑖 + 𝜐 𝑘 , 𝑛 . ( 9 ) The problem is to recover from 𝑦 𝑘 , 𝑛 the data 𝑎 𝑘 , 𝑛 . The presence of the term 𝑏 𝑘 , 𝑛 − 2 𝑖 makes the decoding process from 𝑦 𝑘 , 𝑛 difficult. Some processing should be carried out in order to recover the real data. 4. Interference Cancelation Method 4.1. Cancelation Procedure For the case 𝑁 𝑡 = 2 , it has been shown in [ 8 ] that if we define 𝑧 𝑘 , 𝑛 + 2 as 𝑧 𝑘 , 𝑛 + 2 = 𝐻 ∗ 𝑘 , 1 𝑦 𝑘 , 𝑛 + 𝐻 ∗ 𝑘 , 0 𝑦 𝑘 , 𝑛 + 2 , ( 1 0 ) then we have ℜ  𝑧 𝑘 , 𝑛 + 2  = ℜ  𝐻 ∗ 𝑘 , 1 𝑦 𝑘 , 𝑛 + 𝐻 ∗ 𝑘 , 0 𝑦 𝑘 , 𝑛 + 2  = | | 𝐻 𝑘 , 1 | | 2 𝑎 𝑘 , 𝑛 − 2 + 2 ℜ  𝐻 ∗ 𝑘 , 1 𝐻 𝑘 , 0  𝑎 𝑘 , 𝑛 + | | 𝐻 𝑘 , 0 | | 2 𝑎 𝑘 , 𝑛 + 2 + 𝑤 𝑘 , 𝑛 + 2 , ( 1 1 ) with 𝑤 𝑘 , 𝑛 + 2 = ℜ { 𝐻 ∗ 𝑘 , 1 𝜐 𝑘 , 𝑛 + 𝐻 ∗ 𝑘 , 0 𝜐 𝑘 , 𝑛 + 2 } . Let 2 𝐿 𝑓 denotes the frame length, for 𝑒 ∈ { 0 , 1 } . If we denote by 𝑡 𝑒 =  ℜ { 𝑧 𝑘 , 𝑒 } ℜ { 𝑧 𝑘 , 𝑒 + 2 } ⋯ ℜ { 𝑧 𝑘 , 𝑒 + 2 ( 𝐿 𝑓 − 1 ) }  𝑇 , 𝑎 𝑒 =  𝑎 𝑘 , 𝑒 𝑎 𝑘 , 𝑒 + 2 ⋯ 𝑎 𝑘 , 𝑒 + 2 ( 𝐿 𝑓 − 1 )  𝑇 , 𝑤 𝑒 =  𝑤 𝑘 , 𝑒 𝑤 𝑘 , 𝑒 + 2 ⋯ 𝑤 𝑘 , 𝑒 + 2 ( 𝐿 𝑓 − 1 )  𝑇 ⎡ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎣ | | 𝐻 𝑘 , 0 | | 2 0 ⋯ ⋯ 0 2 ℜ { 𝐻 𝑘 , 0 𝐻 ∗ 𝑘 , 1 } | | 𝐻 𝑘 , 0 | | 2 ⋱ ⋯ ⋮ | | 𝐻 𝑘 , 1 | | 2 2 ℜ  𝐻 𝑘 , 0 𝐻 ∗ 𝑘 , 1  | | 𝐻 𝑘 , 0 | | 2 0 ⋮ 0 ⋱ ⋱ ⋱ ⋮ ⋮ ⋱ ⋱ ⋱ 0 0 ⋯ | | 𝐻 𝑘 , 1 | | 2 2 ℜ  𝐻 𝑘 , 0 𝐻 ∗ 𝑘 , 1  | | 𝐻 𝑘 , 0 | | 2 ⎤ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎦                                                                                                                                                𝐺 2 , ( 1 2 ) ( ( ⋅ ) 𝑇 denotes the transpose operation and ( ⋅ ) 𝐻 the transpose conjugate one) then we have 𝑡 𝑒 = 𝐺 2 𝑎 𝑒 + 𝑤 𝑒 . ( 1 3 ) In this last equation, the imaginary interference term is canceled. Thus the decoding process can be easily carried out by using either Maximum Likelihood (ML) decoding, Viterbi decoding, or linear equalization such as Zero Forcing (ZF) or Minimum Mean Square Error (MMSE) decoding. More generally with 𝑁 𝑡 ≥ 2 , let us note and compute 𝑧 𝑘 , 𝑛 + 2 𝑁 𝑡 − 2 = 𝑁 𝑡 − 1  𝑝 = 0 𝐻 ∗ 𝑘 , 𝑁 𝑡 − 1 − 𝑝 𝑦 𝑘 , 𝑛 + 2 𝑝 = 𝑁 𝑡 − 1  𝑝 = 0 𝑁 𝑡 − 1  𝑖 = 0 𝐻 ∗ 𝑘 , 𝑁 𝑡 − 1 − 𝑝 𝐻 𝑘 , 𝑖 𝑥 𝑘 , 𝑛 + 2 𝑝 − 2 𝑖 + 𝑁 𝑡 − 1  𝑝 = 0 𝐻 ∗ 𝑘 , 𝑁 𝑡 − 1 − 𝑝 𝜐 𝑘 , 𝑛 + 2 𝑝                                    𝑛 𝑘 , 𝑛 + 2 𝑁 𝑡 − 2 = 𝑁 𝑡 − 1  𝑖 = 1 𝑖 − 1  𝑝 = 0 𝐻 ∗ 𝑘 , 𝑁 𝑡 − 1 − 𝑝 𝐻 𝑘 , 𝑖 𝑥 𝑘 , 𝑛 + 2 𝑝 − 2 𝑖                                                            𝐵 𝑘 , 𝑛 + 𝑁 𝑡 − 1  𝑖 = 0  𝑝 = 𝑖 𝐻 ∗ 𝑘 , 𝑁 𝑡 − 1 − 𝑝 𝐻 𝑘 , 𝑖 𝑥 𝑘 , 𝑛 + 2 𝑝 − 2 𝑖                                                            𝐴 𝑘 , 𝑛 + 𝑁 𝑡 − 1  𝑖 = 0 𝑁 𝑡 − 1  𝑝 = 𝑖 + 1 𝐻 ∗ 𝑘 , 𝑁 𝑡 − 1 − 𝑝 𝐻 𝑘 , 𝑖 𝑥 𝑘 , 𝑛 + 2 𝑝 − 2 𝑖                                                                𝐶 𝑘 , 𝑛 + 𝑁 𝑡 − 1  𝑝 = 0 𝐻 ∗ 𝑘 , 𝑁 𝑡 − 1 − 𝑝 𝜐 𝑘 , 𝑛 + 2 𝑝 . ( 1 4 ) Moreover 𝐴 𝑘 , 𝑛 is given by 𝐴 𝑘 , 𝑛 = 𝑥 𝑘 , 𝑛 𝜇 𝑘 , ( 1 5 ) and details for this equation are given in Appendix A.1 . The expression of 𝐵 𝑘 , 𝑛 is given by 𝐵 𝑘 , 𝑛 = 𝑁 𝑡 − 1  𝑞 = 1 𝑥 𝑘 , 𝑛 − 2 𝑞 𝛾 𝑞 , ( 1 6 ) where 𝛾 𝑞 are real-valued quantities which depend only on the channel coefficients as shown in Appendix A.2 . The expression of 𝐶 𝑘 , 𝑛 is given by 𝐶 𝑘 , 𝑛 = 𝑁 𝑡 − 1  𝑞 = 1 𝑥 𝑘 , 𝑛 + 2 𝑞 𝛽 𝑞 , ( 1 7 ) where 𝛽 𝑞 are real-valued quantities which depend only on the channel coefficients as shown in Appendix A.3 . Therefore, 𝑧 𝑘 , 𝑛 + 2 𝑁 𝑡 − 2 = 𝑁 𝑡 − 1  𝑞 = 1 𝛾 𝑞 𝑥 𝑘 , 𝑛 − 2 𝑞 + 𝜇 𝑘 𝑥 𝑘 , 𝑛 + 𝑁 𝑡 − 1  𝑞 = 1 𝛽 𝑞 𝑥 𝑘 , 𝑛 + 2 𝑞 + 𝑁 𝑡 − 1  𝑝 = 0 𝐻 ∗ 𝑘 , 𝑁 𝑡 − 1 − 𝑝 𝜐 𝑘 , 𝑛 + 2 𝑝 . ( 1 8 ) Thus, by noting that 𝑡 ( 1 ) 𝑘 , 𝑛 + 2 𝑁 𝑡 − 2 = ℜ { 𝑧 𝑘 , 𝑛 + 2 𝑁 𝑡 − 2 } , we have 𝑡 ( 1 ) 𝑘 , 𝑛 + 2 𝑁 𝑡 − 2 = 𝑁 𝑡 − 1  𝑞 = 1 𝛾 𝑞 𝑎 𝑘 , 𝑛 − 2 𝑞 + 𝜇 𝑘 𝑎 𝑘 , 𝑛 + 𝑁 𝑡 − 1  𝑞 = 1 𝛽 𝑞 𝑎 𝑘 , 𝑛 + 2 𝑞 + ℜ ⎧ ⎨ ⎩ 𝑁 𝑡 − 1  𝑝 = 0 𝐻 ∗ 𝑘 , 𝑁 𝑡 − 1 − 𝑝 𝜐 𝑘 , 𝑛 + 2 𝑝 ⎫ ⎬ ⎭                                                  𝑤 𝑘 , 𝑛 + 2 𝑁 𝑡 − 2 . ( 1 9 ) For 𝑒 ∈ { 0 , 1 } , we note 𝑡 𝑒 = [ 𝑡 𝑘 , 𝑒 𝑡 𝑘 , 𝑒 + 2 ⋯ 𝑡 𝑘 , 𝑒 + 2 ( 𝐿 𝑓 − 1 ) ] 𝑇 , 𝑤 𝑒 = [ ℜ { 𝑤 𝑘 , 𝑒 } ℜ { 𝑤 𝑘 , 𝑒 + 2 } ⋯ ℜ { 𝑤 𝑘 , 𝑒 + 2 ( 𝐿 𝑓 − 1 ) ] 𝑇 , and 𝐺 𝑁 𝑡 = ⎛ ⎜ ⎜ ⎜ ⎜ ⎜ ⎜ ⎜ ⎜ ⎜ ⎜ ⎜ ⎜ ⎜ ⎜ ⎝ 𝛽 𝑁 𝑡 − 1 0 ⋯ ⋯ ⋯ ⋯ ⋯ ⋯ ⋯ ⋯ 0 𝛽 𝑁 𝑡 − 2 𝛽 𝑁 𝑡 − 1 0 ⋯ ⋯ ⋯ ⋯ ⋯ ⋯ ⋯ ⋮ ⋮ ⋱ ⋱ ⋱ ⋯ ⋯ ⋯ ⋯ ⋯ ⋯ ⋮ 𝛽 1 ⋱ ⋱ ⋱ ⋱ ⋯ ⋯ ⋯ ⋯ ⋯ ⋮ 𝜇 𝑘 ⋱ ⋱ ⋱ ⋱ ⋱ ⋯ ⋯ ⋯ ⋯ ⋮ 𝛾 𝑁 𝑡 − 1 ⋱ ⋱ ⋱ ⋱ ⋱ ⋱ ⋯ ⋯ ⋯ ⋮ ⋮ ⋱ ⋱ ⋱ ⋱ ⋱ ⋱ ⋱ ⋯ ⋯ ⋮ 𝛾 1 ⋱ ⋱ ⋱ ⋱ ⋱ ⋱ ⋱ ⋱ ⋯ ⋮ 0 ⋱ ⋱ ⋱ ⋱ ⋱ ⋱ ⋱ ⋱ ⋱ ⋮ ⋮ ⋱ ⋱ ⋱ ⋱ ⋱ ⋱ ⋱ ⋱ ⋱ 0 0 ⋯ 0 𝛾 1 ⋯ 𝛾 𝑁 𝑡 − 1 𝜇 𝑘 𝛽 1 ⋯ 𝛽 𝑁 𝑡 − 2 𝛽 𝑁 𝑡 − 1 ⎞ ⎟ ⎟ ⎟ ⎟ ⎟ ⎟ ⎟ ⎟ ⎟ ⎟ ⎟ ⎟ ⎟ ⎟ ⎠ . ( 2 0 ) We have: 𝑡 𝑒 = 𝐺 𝑁 𝑡 𝑎 𝑒 + 𝑤 𝑒 . ( 2 1 ) There is no imaginary interference in ( 21 ) and consequently Maximum Likelihood (ML) [ 11 ] or linear equalizers can be used to estimate 𝑎 𝑘 , 𝑛 . The computation of 𝑧 𝑘 , 𝑛 from 𝑦 𝑘 , 𝑛 according to ( 14 ) is referred to as Preprocessing1 as shown in Figure 2 . We will now provide a theoretical performance analysis of this scheme. Figure 2: FBMC Single-delay STTC receiver. 4.2. A Theoretical Performance Analysis Let us consider that the noise 𝜐 𝑘 , 𝑛 is an AWGN noise with 𝐸 { | 𝜐 𝑘 , 𝑛 | 2 } = 𝑁 0 . It is worth noticing that ℜ { 𝑤 𝑘 , 𝑛 } is Gaussian noise as it is the result of the real part of a linear transformation of Gaussian noise. However this noise is colored. For example, when 𝑁 𝑇 = 2 , we have (i) 𝐸 { 𝑤 𝑘 , 𝑛 𝑤 ∗ 𝑘 , 𝑛 + 2 } = 𝐸 { 𝑤 𝑘 , 𝑛 + 2 𝑤 ∗ 𝑘 , 𝑛 } = 𝑁 0 ( | 𝐻 𝑘 , 0 | 2 + | 𝐻 𝑘 , 1 | 2 ) ℜ { ( 𝐻 𝑘 , 0 ) ∗ 𝐻 𝑘 , 1 } / 2 , (ii) 𝐸 { 𝑤 𝑘 , 𝑛 𝑤 ∗ 𝑘 , 𝑛 } = 𝑁 0 ( | 𝐻 𝑘 , 0 | 2 + | 𝐻 𝑘 , 1 | 2 ) / 2 = 𝑈 0 / 2 , (iii) for 𝑞 ≠ { 0 , 1 } , 𝐸 { 𝑤 𝑘 , 𝑛 𝑤 ∗ 𝑘 , 𝑛 + 2 𝑞 } = 0 . Let us recall that if the noise was white the ML performance would have been obtained by the Viterbi decoder. Therefore, the performance of Viterbi decoding in this present case is suboptimal. In [ 12 ] the authors evaluate the loss of performance of Viterbi decoding in presence of correlated noise. The optimal performance using an ML decoding is very complex to implement since it requires an exhaustive search over all the possible transmitted sequences. Another alternative could be to perform a whitening followed by a Viterbi decoding. However, such Viterbi decoding will be more complex since the whitening will increase the number of states. Indeed, the noise 𝑤 𝑒 is colored with a correlation matrix 𝑅 . Since 𝑅 is a positive Hermitian matrix, its eigenvalues are real and positive. We have 𝑅 = 𝑄 ⎡ ⎢ ⎢ ⎢ ⎢ ⎣ 𝜆 0 0 ⋯ 0 0 ⋱ ⋱ ⋮ ⋮ ⋱ ⋱ 0 0 ⋯ 0 𝜆 𝐿 𝑓 − 1 ⎤ ⎥ ⎥ ⎥ ⎥ ⎦                                    Λ 𝑄 𝐻 , ( 2 2 ) with 𝑄 being a unitary matrix, that is, 𝑄 𝑄 𝐻 = 𝐼 𝐿 𝑓 . We denote Λ 1 / 2 = ⎡ ⎢ ⎢ ⎢ ⎢ ⎣ 𝜆 1 / 2 0 0 ⋯ 0 0 ⋱ ⋱ ⋮ ⋮ ⋱ ⋱ 0 0 ⋯ 0 𝜆 1 / 2 𝐿 𝑓 − 1 ⎤ ⎥ ⎥ ⎥ ⎥ ⎦ . ( 2 3 ) Therefore, the whitening process can be done by computing 𝑦 𝑒 = Λ − 1 / 2 𝑄 𝐻 𝑧 𝑒 = Λ − 1 / 2 𝑄 𝐻 𝐺 2                  𝐻 𝑎 𝑒 + Λ − 1 / 2 𝑄 𝐻 𝑤 𝑒                  𝜇 𝑒 = 𝐻 𝑎 𝑒 + 𝜇 𝑒 . ( 2 4 ) It can easily be proved that 𝜇 𝑒 is AWGN. As we will see in the simulation results section, the presence of the colored noise will lead to a degradation of performance. Let us now present an iterative decoding approach which should improve the performance compared to that of the previous decoding strategy. 5. Iterative Method 5.1. Iterative Procedure In this section we propose an iterative decoding procedure for FBMC single-delay STTC decoding. At the output of the Preprocessing1 block (see Figure 3 ), we can perform a decoding procedure (ML, Viterbi, or linear decoding) to derive an estimate value ̂ 𝑎 ( 1 ) 𝑘 , 𝑛 of 𝑎 𝑘 , 𝑛 . From ( 6 ) and using this estimate ̂ 𝑎 ( 1 ) 𝑘 , 𝑛 , we can compute an estimate ̂ 𝑢 ( 1 ) 𝑘 , 𝑛 of 𝑢 𝑘 , 𝑛 by Figure 3: Receiver decoding processing for FBMC modulation in the case of single delay STTC transmission. ̂ 𝑢 ( 1 ) 𝑘 , 𝑛 =  ( 𝑝 , 𝑞 ) ≠ ( 0 , 0 ) ̂ 𝑎 ( 1 ) 𝑘 + 𝑝 , 𝑛 + 𝑞 ∞  𝑚 = − ∞ 𝑔 𝑘 , 𝑛 [ 𝑚 ] 𝑔 𝑘 + 𝑝 , 𝑛 + 𝑞 [ 𝑚 ]                                            𝛾 𝑝 , 𝑞 . ( 2 5 ) It is worth noticing that for a well-localized prototype filter in time and frequency domain it is enough to consider the previous sum only for 𝑝 , 𝑞 ∈ { 1 , − 1 } , that is, ̂ 𝑢 ( 1 ) 𝑘 , 𝑛 ≈  | | 𝑝 | | = 1 , | | 𝑞 | | = 1 ̂ 𝑎 ( 1 ) 𝑘 + 𝑝 , 𝑛 + 𝑞 𝛾 𝑝 , 𝑞 . ( 2 6 ) This approximation is justified in [ 10 ]. 𝛾 𝑝 , 𝑞 can be computed off-line since the prototype filter response is known. Then in ( 9 ) we can remove the contribution of the 𝑢 𝑘 , 𝑛 components by computing 𝑦 ( 2 ) 𝑘 , 𝑛 = 𝑦 𝑘 , 𝑛 − 𝑁 𝑡 − 1  𝑖 = 0 𝐻 𝑘 , 𝑖 ̂ 𝑢 ( 1 ) 𝑘 , 𝑛 − 2 𝑖 = 𝑁 𝑡 − 1  𝑖 = 0 𝐻 𝑘 , 𝑖 𝑎 𝑘 , 𝑛 − 2 𝑖 + 𝑁 𝑡 − 1  𝑖 = 0 𝑗 𝐻 𝑘 , 𝑖  𝑢 𝑘 , 𝑛 − 2 𝑖 − ̂ 𝑢 ( 1 ) 𝑘 , 𝑛 − 2 𝑖  + 𝜐 𝑘 , 𝑛 . ( 2 7 ) If we assume a perfect cancelation of the 𝑢 𝑘 , 𝑛 terms, that is, 𝑢 𝑘 , 𝑛 = ̂ 𝑢 ( 1 ) 𝑘 , 𝑛 , then we have 𝑦 ( 2 ) 𝑘 , 𝑛 = 𝑁 𝑡 − 1  𝑖 = 0 𝐻 𝑘 , 𝑖 𝑎 𝑘 , 𝑛 − 2 𝑖 + 𝜐 𝑘 , 𝑛 . ( 2 8 ) The operation of estimating 𝑢 𝑘 , 𝑛 and canceling its contribution to the signal 𝑦 𝑘 , 𝑛 is referred to as “Interference estimation + Interference cancelation” as depicted in Figure 3 . Thus, we can perform from 𝑦 ( 2 ) 𝑘 , 𝑛 a new decoding (Decoder 2 block) to obtain a new estimate ̂ 𝑎 ( 2 ) 𝑘 , 𝑛 of 𝑎 𝑘 , 𝑛 . In the same manner, we can use either a Viterbi/ML decoding or a linear decoder. From ̂ 𝑎 ( 2 ) 𝑘 , 𝑛 and ( 19 ) we can also compute 𝑡 ( 2 ) 𝑘 , 𝑛 + 2 by 𝑡 ( 2 ) 𝑘 , 𝑛 + 2 𝑁 𝑡 − 2 = 𝑁 𝑡 − 1  𝑞 = 1 𝛾 𝑞 ̂ 𝑎 ( 2 ) 𝑘 , 𝑛 − 2 𝑞 + 𝜇 𝑘 ̂ 𝑎 ( 2 ) 𝑘 , 𝑛 + 𝑁 𝑡 − 1  𝑞 = 1 𝛽 𝑞 ̂ 𝑎 ( 2 ) 𝑘 , 𝑛 + 2 𝑞 . ( 2 9 ) 𝑡 ( 2 ) 𝑘 , 𝑛 + 2 can also be rewritten as 𝑡 ( 2 ) 𝑘 , 𝑛 + 2 𝑁 𝑡 − 2 = 𝑁 𝑡 − 1  𝑞 = 1 𝛾 𝑞 𝑎 𝑘 , 𝑛 − 2 𝑞 + 𝜇 𝑘 𝑎 𝑘 , 𝑛 + 𝑁 𝑡 − 1  𝑞 = 1 𝛽 𝑞 𝑎 𝑘 , 𝑛 + 2 𝑞 + 𝑁 𝑡 − 1  𝑞 = 1 𝛾 𝑞  ̂ 𝑎 ( 2 ) 𝑘 , 𝑛 − 2 𝑞 − 𝑎 𝑘 , 𝑛 − 2 𝑞  + 𝜇 𝑘  ̂ 𝑎 ( 2 ) 𝑘 , 𝑛 − 𝑎 𝑘 , 𝑛  + 𝑁 𝑡 − 1  𝑞 = 1 𝛽 𝑞  ̂ 𝑎 ( 2 ) 𝑘 , 𝑛 + 2 𝑞 − 𝑎 𝑘 , 𝑛 + 2 𝑞  = 𝑁 𝑡 − 1  𝑞 = 1 𝛾 𝑞 𝑎 𝑘 , 𝑛 − 2 𝑞 + 𝜇 𝑘 𝑎 𝑘 , 𝑛 + 𝑁 𝑡 − 1  𝑞 = 1 𝛽 𝑞 𝑎 𝑘 , 𝑛 + 2 𝑞 + n o i s e c o m p o n e n t . ( 3 0 ) 𝑡 ( 2 ) 𝑘 , 𝑛 + 2 𝑁 𝑡 − 2 is a new version of the 𝑡 ( 1 ) 𝑘 , 𝑛 + 2 𝑁 𝑡 − 2 signal which is obtained from the estimates of the Decoder 2 block output. Thus, this last equation can be used to perform another estimation ̂ 𝑎 ( 3 ) 𝑘 , 𝑛 of 𝑎 𝑘 , 𝑛 in the same manner as we compute ̂ 𝑎 ( 1 ) 𝑘 , 𝑛 . We expect to improve the estimation of 𝑎 𝑘 , 𝑛 since the noise component in ( 30 ) should be less correlated than the one in ( 19 ). Again from ̂ 𝑎 ( 3 ) 𝑘 , 𝑛 we can derive an estimate ̂ 𝑢 ( 2 ) 𝑘 , 𝑛 of 𝑢 𝑘 , 𝑛 as in ( 25 ). Therefore, we can repeat another decoding process as already presented. We can run this decoding process as many times as necessary. The process of computing 𝑡 ( 2 ) 𝑘 , 𝑛 + 2 𝑁 𝑡 − 2 from the ̂ 𝑎 ( 2 ) 𝑘 , 𝑛 is referred to as Preprocessing2; see Figure 3 . Let us have a look at the convergence of this iterative method. 5.2. A Convergence Analysis of the Iterative Procedure Let us consider the function 𝑃 𝑒 = 𝐶 1 ( S N R ) that we obtain when considering the perfect cancelation of the interference term by using ( 28 ) and the function 𝑃 𝑒 = 𝐶 2 ( S N R ) obtained using ( 19 ). 𝑃 𝑒 is the real symbol error probability and S N R = 2 𝜎 2 𝑎 / 𝑁 0 = 1 / 𝑁 0 assuming that the real symbol power 𝜎 2 𝑎 is fixed at 1/2. These functions are illustrated in Figure 4 for a given channel realization. Let us note that 𝐶 1 is Δ dB better than 𝐶 2 , that is, Figure 4: Convergence illustration. 𝐶 1 ⎛ ⎜ ⎝ 1  1 + 𝛼 Δ  𝑁 0 ⎞ ⎟ ⎠ = 𝐶 2  1 𝑁 0  , ( 3 1 ) with Δ = 1 0 l o g 1 0 ( 1 + 𝛼 Δ ) . At the first iteration, when using ( 19 ) for decoding, we obtain at S N R = 1 / 𝑁 0 a symbol probability of error 𝑃 𝑒 1 = 𝐶 2 ( 1 / 𝑁 0 ) . This first iteration is summarized by the point 𝐴 1 ( 1 / 𝑁 0 , 𝑃 𝑒 1 ) in Figure 4 . Now, from this probability of error we can derive the degradation that we obtain when applying interference cancelation. Indeed, the cancelation of the interference will add some noise to the current noise component. This additional noise component is given by the cancelation error 𝑛 + = 𝑁 𝑡 − 1  𝑖 = 0 𝑗 𝐻 𝑘 , 𝑖  𝑢 𝑘 , 𝑛 − 2 𝑖 − ̂ 𝑢 ( 1 ) 𝑘 , 𝑛 − 2 𝑖  = 𝑁 𝑡 − 1  𝑖 = 0 𝑗 𝐻 𝑘 , 𝑖  ( 𝑝 , 𝑞 ) ≠ ( 0 , 0 )  𝑎 𝑘 + 𝑝 , 𝑛 − 2 𝑖 + 𝑞 − ̂ 𝑎 ( 1 ) 𝑘 + 𝑝 , 𝑛 − 2 𝑖 + 𝑞  × ∞  𝑚 = − ∞ 𝑔 𝑘 , 𝑛 − 2 𝑖 [ 𝑚 ] 𝑔 𝑘 + 𝑝 , 𝑛 − 2 𝑖 + 𝑞 [ 𝑚 ] . ( 3 2 ) Using the current observation ̂ 𝑎 𝑘 , 𝑛 = 𝑎 𝑘 , 𝑛 w i t h p r o b a b i l i t y 1 − 𝑃 𝑒 1 , ̂ 𝑎 𝑘 , 𝑛 ≠ 𝑎 𝑘 , 𝑛 w i t h p r o b a b i l i t y 𝑃 𝑒 1 ( 3 3 ) and considering that [ 10 ]  ( 𝑝 , 𝑞 ) ≠ ( 0 , 0 ) | | | | ∞  𝑚 = − ∞ 𝑔 𝑘 , 𝑛 − 2 𝑖 [ 𝑚 ] 𝑔 𝑘 + 𝑝 , 𝑛 − 2 𝑖 + 𝑞 [ 𝑚 ] | | | | 2 = 1 , ( 3 4 ) we have 𝐸  | | 𝑛 + | | 2  = 𝑃 𝑒 1 𝑁 𝑡 − 1  𝑖 = 0 | | 𝐻 𝑘 , 𝑖 | | 2              𝛼 ℎ . ( 3 5 ) Therefore, the symbol probability of error is given at second iteration by 𝑃 𝑒 2 = 𝐶 1 ⎛ ⎜ ⎝ 1 𝑁 1 0 ⎞ ⎟ ⎠ = 𝐶 1  1 𝑁 0 + 𝑃 𝑒 1 𝛼 ℎ  = 𝐶 1 ⎛ ⎜ ⎝ 1 𝑁 0  1 + 𝛼 ℎ 𝐶 2  1 / 𝑁 0  / 𝑁 0  ⎞ ⎟ ⎠ , ( 3 6 ) where 1 / 𝑁 1 0 is the SNR at the input of Decoder 2. 𝐶 2 ( 1 / 𝑁 0 ) is a 𝑄 -function that is exponentially decreasing as SNR increases; thus, 𝛼 ℎ 𝐶 2 ( 1 / 𝑁 0 ) / 𝑁 0 decreases as SNR increases since the exponential function overwhelms the polynomial function. Then, there is a noise power 𝑁 𝑎 0 such that, for 𝑁 0 < 𝑁 𝑎 0 𝛼 ℎ 𝐶 2  1 / 𝑁 0  𝑁 0 < 𝛼 Δ , ( 3 7 ) and thus, 1 𝑁 0  1 + 𝛼 ℎ 𝐶 2  1 / 𝑁 0  / 𝑁 0  > 1  1 + 𝛼 Δ  𝑁 0 . ( 3 8 ) Therefore for 𝑁 0 < 𝑁 𝑎 0 , 𝐶 1 ⎛ ⎜ ⎝ 1 𝑁 0  1 + 𝛼 ℎ 𝐶 2  1 / 𝑁 0  / 𝑁 0  ⎞ ⎟ ⎠ < 𝐶 1 ⎛ ⎜ ⎝ 1  1 + 𝛼 Δ  𝑁 0 ⎞ ⎟ ⎠ = 𝐶 2  1 𝑁 0  , ( 3 9 ) that is, 𝑃 𝑒 2 < 𝑃 𝑒 1 . ( 4 0 ) For 𝑁 0 < 𝑁 𝑎 0 the output of the second iteration will give better performance than that of the first iteration. This second iteration is summarized by the point 𝐴 2 ( 1 / 𝑁 1 0 , 𝑃 𝑒 2 ) in Figure 4 . When recombining the signal at the input of Decoder 1 for the third iteration using ( 29 ), the noise component is now smaller than that in the previous case since 𝑃 𝑒 2 < 𝑃 𝑒 1 . Consequently, the third iteration performance is given by 𝐶 2 at S N R = 1 / 𝑁 2 0 with 𝑁 2 0 < 𝑁 1 0 . Thus, 𝐶 2 ( 1 / 𝑁 2 0 ) < 𝐶 2 ( 1 / 𝑁 0 ) , that is, the probability of error at the output of Decoder 1 for the third iteration 𝑃 𝑒 3 is less than that for 𝑃 𝑒 1 . This third iteration is summarized by the point 𝐴 3 ( 1 / 𝑁 2 0 , 𝑃 𝑒 3 ) in Figure 4 . Let us notice that 𝑃 𝑒 3 could be greater than 𝑃 𝑒 2 . The next iteration performance can be derived in the same manner since we just have to replace 𝑁 0 by 𝑁 2 0 . Thus, the probability of error at the output of a given decoder (Decoder 1 or Decoder 2) will always decrease or reach a fixed point. 6. Simulation Results In this section, we will evaluate the performance of the two decoding methods that we have presented. We consider a transmission scheme with two and three transmit antennas. For 𝑁 𝑡 = 2 , we have 𝑡 ( 1 ) 𝑘 , 𝑛 + 2 𝑁 𝑡 − 2 = | | 𝐻 𝑘 , 1 | | 2 𝑎 𝑘 , 𝑛 − 2 + 2 ℜ  𝐻 ∗ 𝑘 , 1 𝐻 𝑘 , 0  𝑎 𝑘 , 𝑛 + | | 𝐻 𝑘 , 0 | | 2 𝑎 𝑘 , 𝑛 + 2 + 𝑤 𝑘 , 𝑛 + 2 𝑁 𝑡 − 2 , ( 4 1 ) and for 𝑁 𝑡 = 3 , we get 𝑡 ( 1 ) 𝑘 , 𝑛 + 2 𝑁 𝑡 − 2 = | | 𝐻 𝑘 , 2 | | 2 𝑎 𝑘 , 𝑛 − 4 + 2 ℜ  𝐻 ∗ 𝑘 , 2 𝐻 𝑘 , 1  𝑎 𝑘 , 𝑛 − 2 +  2 ℜ  𝐻 ∗ 𝑘 , 2 𝐻 𝑘 , 0  + | | 𝐻 𝑘 , 1 | | 2  𝑎 𝑘 , 𝑛 + 2 ℜ  𝐻 ∗ 𝑘 , 1 𝐻 𝑘 , 0  𝑎 𝑘 , 𝑛 + 2 + | | 𝐻 𝑘 , 0 | | 2 𝑎 𝑘 , 𝑛 + 4 . ( 4 2 ) The simulation parameters we consider are given as follows: (i) no channel coding, (ii) QPSK modulation, (iii) Rayleigh channel per antenna, that is, flat over all the subcarriers. We assume that the channel coefficients are perfectly known by the receiver, (iv) number of subcarrier 𝑀 = 3 2 , (v) we used a truncation of the IOTA (Isotropic Orthogonal Transform Algorithm) prototype function [ 1 ]. Its duration is limited to 4 𝑇 0 , which leads to a nearly orthogonal prototype filter containing 𝐿 = 4 𝑀 = 1 2 8 taps. In this section, we give BER (Bit Error Rate) versus SNR simulation results, and consequently, we do not take into consideration the loss of efficiency due to the cyclic prefix in CP-OFDM modulation. In Figure 5 we show the performance of the FBMC decoding structure introduced in Figure 2 . For FBMC, we consider both ML and Viterbi decoding. ML decoding using an exhaustive search among all possible transmitted sequences of data outperforms Viterbi decoding by 1 dB. This is due to the fact that the noise is colored; thus, Viterbi decoding is suboptimal. We also give the CP-OFDM performance using a Viterbi decoding. We can see that CP-OFDM outperforms ML/FBMC by about 1 dB. Figure 5: Performance of single delay STTC with 2 transmit antennas and one receive antenna (FBMC and CP-OFDM modulation). In the rest of this section, we will focus on the iterative decoding performance. The simulation results are obtained using Viterbi decoding blocks implemented inside Decoder 1 and Decoder 2 blocks in Figure 3 . The Viterbi algorithm implemented in Decoder 1 is related to ( 41 ). For QPSK modulation, the Trellis is a 4 𝑁 𝑇 − 1 state Trellis with only two possible transitions per state since the detection is performed on real data. Whereas the Viterbi algorithm implemented in Decoder 2 is related to ( 28 ) and is a 2 𝑁 𝑇 − 1 state Trellis with two transitions per state, again detection is performed on real data. We also consider hard estimation of the data at the output of a given Viterbi decoder. For the CP-OFDM case with QPSK modulation, we have a 4 𝑁 𝑇 − 1 state Trellis with 4 transitions per state as the detection is performed on complex data. Therefore, this Viterbi algorithm is more complex compared with one of the two Viterbi algorithms used in the case of FBMC modulation. The two Viterbi algorithms used in FBMC taken together have a complexity comparable to the one used with CP-OFDM. However, the two Viterbi algorithms used in FBMC operate on a frame sequence which is two times longer than the one for CP-OFDM modulation. Then, in terms of complexity the proposed FBMC structure has a significantly higher complexity than that of CP-OFDM mainly due to the “Interference estimation + Interference cancelation” block. For uncorrelated Rayleigh channels, we plot the performance of this FBMC receiver structure for different iteration stages as well as the performance of CP-OFDM with ML decoding as a matter of comparison. Figures 6 and 7 provide the simulation results for 𝑁 𝑡 = 2 and 𝑁 𝑡 = 3 , respectively. For 𝑛 = 1 , we have a 2 dB degradation compared to CP-OFDM. For 𝑛 ≥ 2 (more than two-Viterbi decoding), we get closer to CP-OFDM. For 𝑛 = 5 or 6 , we almost reach the same performance as that of CP-OFDM. In Figure 6 we also plot the curve obtained when we assume perfect interference cancelation in the second iteration as mentioned in ( 28 ). In that case, there is a possible gain of 0 . 8 dB since the Viterbi structure with 2 states and two transitions per state (Decoder 2) provides better performance than the 4-state Viterbi decoder with 4 transitions per state implemented for CP-OFDM. Indeed, it is possible to show that the structures of the code related to these two Trellises have the same minimum distance. However, the performance gain is due to the distance distribution associated to the two Trellises. Figure 6: Performance of single delay STTC (iterative decoding) with 2 transmit antennas and one receive antenna (FBMC and CP-OFDM modulation). Figure 7: Performance of single delay STTC (iterative decoding) with 3 transmit antennas and one receive antenna (FBMC and CP-OFDM modulation). Moreover, let us evaluate this scheme in presence of a frequency selective channel. We consider the following channel parameters: (i) uncoded QPSK modulation, (ii) 𝑀 = 6 4 subcarriers, (iii) static channels (no Doppler), IOTA prototype, (iv) 3-tap channels between the transmit antennas and the receive: power profile: 0, −4, −10 (dB) Delay: 0, 1, 2 (number of samples), (v) OFDM Cyclic Prefix length: 4 samples, (vi) perfect channel estimation. As shown in Figure 8 , after one iteration ( 𝑛 = 1 ), we have about 2 dB degradation compared to CP-OFDM. For 𝑛 = 2 and 𝑛 = 3 , the loss is reduced to 0 . 7 dB, and for 𝑛 = 6 the degradation is about 0 . 3 dB compared to CP-OFDM. However, the iterative method has an inherent gain as FBMC does not use a CP contrary to CP-OFDM. Figure 8: Performance of single delay STTC (iterative decoding) with 2 transmit antennas and one receive antenna over frequency selective channels. 7. Conclusion In this paper, we have presented two general methods for data detection when combining FBMC and single delay STTC as well as the interference cancelation and the iterative methods. The interference cancelation method despite its simplicity has poorer performance compared to that of CP-OFDM. Thus, we have proposed an iterative decoding based on interference estimation and cancelation which does not require any channel coding or decoding block. We have shown that in the case of QPSK modulation and Rayleigh or frequency selective channels it is possible with this decoding method to perform as better as OFDM-STTC. Moreover if the iterative cancelation process is improved, then a potential gain can be achieved. This is obtained with a relatively higher complexity. In future work, we will look at FBMC with other STTC schemes and evaluate their performance under nonlocally flat channels. Appendix A. General Expression of 𝐴 𝑘 , 𝑛 , 𝐵 𝑘 , 𝑛 , and 𝐶 𝑘 , 𝑛 A.1. 𝐴 𝑘 , 𝑛 Value Let us compute 𝐴 𝑘 , 𝑛 . ( i ) Case 𝑁 𝑡 is even, that is, 𝑁 𝑡 = 2 𝑈 𝑡 , 𝐴 𝑘 , 𝑛 = 𝑥 𝑘 , 𝑛 ⎛ ⎜ ⎝ 𝑈 𝑡 − 1  𝑖 = 0 𝐻 ∗ 𝑘 , 2 𝑈 𝑡 − 1 − 𝑖 𝐻 𝑘 , 𝑖 + 2 𝑈 𝑡 − 1  𝑖 = 𝑈 𝑡 𝐻 ∗ 𝑘 , 2 𝑈 𝑡 − 1 − 𝑖 𝐻 𝑘 , 𝑖 ⎞ ⎟ ⎠ . ( A . 1 ) Using the relation 𝑞 = 2 𝑈 𝑡 − 1 − 𝑖 , we have 𝐴 𝑘 , 𝑛 = 𝑥 𝑘 , 𝑛 ⎛ ⎜ ⎝ 𝑈 𝑡 − 1  𝑖 = 0 𝐻 ∗ 𝑘 , 2 𝑈 𝑡 − 1 − 𝑖 𝐻 𝑘 , 𝑖 + 𝑈 𝑡 − 1  𝑞 = 0 𝐻 ∗ 𝑘 , 𝑞 𝐻 𝑘 , 2 𝑈 𝑡 − 1 − 𝑞 ⎞ ⎟ ⎠ = 𝑥 𝑘 , 𝑛 ⎛ ⎜ ⎝ 𝑈 𝑡 − 1  𝑖 = 0  𝐻 ∗ 𝑘 , 2 𝑈 𝑡 − 1 − 𝑖 𝐻 𝑘 , 𝑖 + 𝐻 ∗ 𝑘 , 𝑖 𝐻 𝑘 , 2 𝑈 𝑡 − 1 − 𝑖  ⎞ ⎟ ⎠ = 𝑥 𝑘 , 𝑛 ⎛ ⎜ ⎝ 2 𝑈 𝑡 − 1  𝑖 = 0 ℜ   𝐻 ∗ 𝑘 , 2 𝑈 𝑡 − 1 − 𝑖 𝐻 𝑘 , 𝑖   ⎞ ⎟ ⎠                                                            𝜇 𝑘 . ( A . 2 ) ( i i ) Case 𝑁 𝑡 is odd, that is, 𝑁 𝑡 = 2 𝑈 𝑡 + 1 , 𝜇 𝑘 = ⎛ ⎜ ⎝ 𝑈 𝑡 − 1  𝑖 = 0 𝐻 ∗ 𝑘 , 2 𝑈 𝑡 − 𝑖 𝐻 𝑘 , 𝑖 + 𝐻 ∗ 𝑘 , 𝑈 𝑡 𝐻 𝑘 , 𝑈 𝑡 + 2 𝑈 𝑡  𝑖 = 𝑈 𝑡 + 1 𝐻 ∗ 𝑘 , 2 𝑈 𝑡 − 𝑖 𝐻 𝑘 , 𝑖 ⎞ ⎟ ⎠ . ( A . 3 ) Again using 𝑞 = 2 𝑈 𝑡 − 𝑖 , we have 𝜇 𝑘 = ⎛ ⎜ ⎝ 2 𝑈 𝑡 − 1  𝑖 = 0 ℜ   𝐻 ∗ 𝑘 , 2 𝑈 𝑡 − 𝑖 𝐻 𝑘 , 𝑖   + | | 𝐻 𝑘 , 𝑈 𝑡 | | 2 ⎞ ⎟ ⎠ , ( A . 4 ) and we get 𝐴 𝑘 , 𝑛 = 𝑥 𝑘 , 𝑛 𝜇 𝑘 . ( A . 5 ) A.2. 𝐵 𝑘 , 𝑛 Value Let us now compute 𝐵 𝑘 , 𝑛 ; setting 𝑞 = 𝑝 − 𝑖 , we get 𝐵 𝑘 , 𝑛 = 𝑁 𝑡 − 1  𝑖 = 1 𝑖 − 1  𝑝 = 0 𝐻 ∗ 𝑘 , 𝑁 𝑡 − 1 − 𝑝 𝐻 𝑘 , 𝑖 𝑥 𝑘 , 𝑛 + 2 𝑝 − 2 𝑖 = 𝑁 𝑡 − 1  𝑖 = 1 𝑖  𝑞 = 1 𝑥 𝑘 , 𝑛 − 2 𝑞 𝐻 ∗ 𝑘 , 𝑁 𝑡 − 1 + 𝑞 − 𝑖 𝐻 𝑘 , 𝑖 . ( A . 6 ) This last equation is the sum over a triangular set of index; therefore, the sum can be taken either from lines or from columns where the total is the same. Therefore, 𝐵 𝑘 , 𝑛 = 𝑁 𝑡 − 1  𝑞 = 1 𝑥 𝑘 , 𝑛 − 2 𝑞 𝑁 𝑡 − 1  𝑖 = 𝑞 𝐻 ∗ 𝑘 , 𝑁 𝑡 − 1 + 𝑞 − 𝑖 𝐻 𝑘 , 𝑖 . ( A . 7 ) Taking 𝑚 = 𝑖 − 𝑞 , we get 𝐵 𝑘 , 𝑛 = 𝑁 𝑡 − 1  𝑞 = 1 𝑥 𝑘 , 𝑛 − 2 𝑞 𝑁 𝑡 − 1 − 𝑞  𝑚 = 0 𝐻 ∗ 𝑘 , 𝑁 𝑡 − 1 − 𝑚 𝐻 𝑘 , 𝑚 + 𝑞                                          𝛾 𝑞 . ( A . 8 ) ( i ) Case 𝑁 𝑡 − 𝑞 is even, that is, 𝑁 𝑡 − 𝑞 = 2 𝑈 𝑞 ; then, 𝛾 𝑞 = 𝑁 𝑡 − 1 − 𝑞  𝑚 = 0 𝐻 ∗ 𝑘 , 𝑁 𝑡 − 1 − 𝑚 𝐻 𝑘 , 𝑚 + 𝑞 = 𝑈 𝑞 − 1  𝑚 = 0 𝐻 ∗ 𝑘 , 2 𝑈 𝑞 + 𝑞 − 1 − 𝑚 𝐻 𝑘 , 𝑚 + 𝑞 + 2 𝑈 𝑞 − 1  𝑚 = 𝑈 𝑞 𝐻 ∗ 𝑘 , 2 𝑈 𝑞 + 𝑞 − 1 − 𝑚 𝐻 𝑘 , 𝑚 + 𝑞 = 2 𝑈 𝑞 − 1  𝑚 = 0 ℜ  𝐻 ∗ 𝑘 , 2 𝑈 𝑞 + 𝑞 − 1 − 𝑚 𝐻 𝑘 , 𝑚 + 𝑞  . ( A . 9 ) ( i i ) Case 𝑁 𝑡 − 𝑞 is odd, that is, 𝑁 𝑡 − 𝑞 = 2 𝑈 𝑞 + 1 ; then, 𝛾 𝑞 = 𝑁 𝑡 − 1 − 𝑞  𝑚 = 0 𝐻 ∗ 𝑘 , 𝑁 𝑡 − 1 − 𝑚 𝐻 𝑘 , 𝑚 + 𝑞 = 𝑈 𝑞 − 1  𝑚 = 0 𝐻 ∗ 𝑘 , 2 𝑈 𝑞 + 𝑞 − 𝑚 𝐻 𝑘 , 𝑚 + 𝑞 + 𝐻 ∗ 𝑘 , 𝑈 𝑞 + 𝑞 𝐻 𝑘 , 𝑈 𝑞 + 𝑞 + 2 𝑈 𝑞  𝑚 = 𝑈 𝑞 + 1 𝐻 ∗ 𝑘 , 2 𝑈 𝑞 + 𝑞 − 𝑚 𝐻 𝑘 , 𝑚 + 𝑞 = 2 𝑈 𝑞 − 1  𝑚 = 0 ℜ  𝐻 ∗ 𝑘 , 2 𝑈 𝑞 + 𝑞 − 𝑚 𝐻 𝑘 , 𝑚 + 𝑞  + | | | 𝐻 ∗ 𝑘 , 𝑈 𝑞 + 𝑞 | | | 2 . ( A . 1 0 ) A.3. 𝐶 𝑘 , 𝑛 Value Let us now compute 𝐶 𝑘 , 𝑛 ; setting 𝑞 = 𝑝 − 𝑖 , we get 𝐶 𝑘 , 𝑛 = 𝑁 𝑡 − 2  𝑖 = 1 𝑁 𝑡 − 1 − 𝑖  𝑞 = 1 𝐻 ∗ 𝑘 , 𝑁 𝑡 − 1 − 𝑞 − 𝑖 𝐻 𝑘 , 𝑖 𝑥 𝑘 , 𝑛 + 2 𝑞 . ( A . 1 1 ) This last equation is the sum over a triangular set of index; therefore, the sum can be taken either from lines or from columns where the total is the same. Therefore, 𝐶 𝑘 , 𝑛 = 𝑁 𝑡 − 1  𝑞 = 1 𝑥 𝑘 , 𝑛 + 2 𝑞 𝑁 𝑡 − 1 − 𝑞  𝑖 = 0 𝐻 ∗ 𝑘 , 𝑁 𝑡 − 1 − 𝑞 − 𝑖 𝐻 𝑘 , 𝑖                                        𝛽 𝑞 . ( A . 1 2 ) ( i ) Case 𝑁 𝑡 − 𝑞 is even, that is, 𝑁 𝑡 − 𝑞 = 2 𝑈 𝑞 ; then, 𝛽 𝑞 = 𝑁 𝑡 − 1 − 𝑞  𝑖 = 0 𝐻 ∗ 𝑘 , 𝑁 𝑡 − 1 − 𝑞 − 𝑖 𝐻 𝑘 , 𝑖 = 𝑈 𝑞 − 1  𝑖 = 0 𝐻 ∗ 𝑘 , 2 𝑈 𝑞 − 1 − 𝑖 𝐻 𝑘 , 𝑖 + 2 𝑈 𝑞 − 1  𝑚 = 𝑈 𝑞 𝐻 ∗ 𝑘 , 2 𝑈 𝑞 − 1 − 𝑖 𝐻 𝑘 , 𝑖 = 2 𝑈 𝑞 − 1  𝑖 = 0 ℜ  𝐻 ∗ 𝑘 , 2 𝑈 𝑞 − 1 − 𝑖 𝐻 𝑘 , 𝑖  . ( A . 1 3 ) ( i i ) Case 𝑁 𝑡 − 𝑞 is odd, that is, 𝑁 𝑡 − 𝑞 = 2 𝑈 𝑞 + 1 ; then, 𝛽 𝑞 = 𝑁 𝑡 − 1 − 𝑞  𝑖 = 0 𝐻 ∗ 𝑘 , 𝑁 𝑡 − 1 − 𝑞 − 𝑖 𝐻 𝑘 , 𝑖 = 𝑈 𝑞 − 1  𝑖 = 0 𝐻 ∗ 𝑘 , 2 𝑈 𝑞 − 𝑖 𝐻 𝑘 , 𝑖 + 𝐻 ∗ 𝑘 , 𝑈 𝑞 𝐻 𝑘 , 𝑈 𝑞 + 2 𝑈 𝑞  𝑚 = 𝑈 𝑞 + 1 𝐻 ∗ 𝑘 , 2 𝑈 𝑞 − 𝑖 𝐻 𝑘 , 𝑖 = 2 𝑈 𝑞 − 1  𝑚 = 0 ℜ  𝐻 ∗ 𝑘 , 2 𝑈 𝑞 − 𝑖 𝐻 𝑘 , 𝑖  + | | | 𝐻 ∗ 𝑘 , 𝑈 𝑞 | | | 2 . ( A . 1 4 ) Acknowledgments The authors would like to thank Pr. M. Bellanger for helpful discussions. This work was supported in part by the European Commission under Project PHYDYAS (FP7-ICT-2007-1-211887). <h4>References</h4> B. Le Floch, M. Alard, and C. Berrou, “Coded orthogonal frequency division multiplex,” Proceedings of the IEEE , vol. 83, pp. 982–996, 1995. H. Boelcskei, “Orthogonal frequency division multiplexing based on offset QAM,” in Advances in Gabor Analysis , Birkhäuser, Boston, Mass, USA, 2003. P. Siohan, C. Siclet, and N. Lacaille, “ Analysis and design of OFDM/OQAM systems based on filterbank theory ,” IEEE Transactions on Signal Processing , vol. 50, no. 5, pp. 1170–1183, 2002. R. W. Chang, “Synthesis of band-limited orthogonal signals for multi-channel data transmission,” Bell Labs Technical Journal , vol. 45, pp. 1775–1796, 1966. B. R. Saltzberg, “Performance of an efficient parallel data transmission system,” IEEE Transactions on Communication Technology , vol. 15, no. 6, pp. 805–811, 1967. V. Tarokh, N. Seshadri, and A. R. Calderbank, “Space-time codes for high data rate wireless communication: performance criterion and code construction,” IEEE Transactions on Information Theory , vol. 44, no. 2, pp. 744–765, 1998. V. Tarokh, H. Jafarkhani, and A. R. Calderbank, “ Space-time block coding for wireless communications: performance results ,” IEEE Journal on Selected Areas in Communications , vol. 17, no. 3, pp. 451–460, 1999. M. Bellanger, “ Transmit diversity in multicarrier transmission using OQAM modulation ,” in Proceedings of the 3rd International Symposium on Wireless Pervasive Computing (ISWPC '08) , pp. 727–730, Santorini, Greece, May 2008. B. Hirosaki, “Orthogonally multiplexed QAM system using the discrete Fourier transform,” IEEE Transactions on Communications Systems , vol. 29, no. 7, pp. 982–989, 1981. C. Lélé, P. Siohan, R. Legouable, and J.-P. Javaudin, “ Preamble-based channel estimation techniques for OFDM/OQAM over the powerline ,” in Proceedings of IEEE International Symposium on Power Line Communications and Its Applications (ISPLC '07) , pp. 59–64, Pisa, Italy, March 2007. G. D. Forney Jr., “Maximum-likelihood sequence estimation of digital sequences in the presence of intersymbol interference,” IEEE Transactions on Information Theory , vol. 18, no. 3, pp. 363–378, 1972. J.-D. Wang and H. Y. Chung, “ Trellis coded communication systems in the presence of colored noise: performance analysis, simulation, and the swapping technique ,” in Proceeingds of IEEE Conference Record on Global Telecommunications Conference, and Exhibition. ‘Communications for the Information Age’ (GLOBECOM '88) , vol. 2, pp. 1160–1165, Hollywood, Fla, USA, November-December 1988. //

Journal

EURASIP Journal on Advances in Signal ProcessingHindawi Publishing Corporation

Published: Feb 18, 2010

There are no references for this article.