Open Access Research

Factor graph based detection approach for high-mobility OFDM systems with large FFT modes

Pello Ochandiano1, Henk Wymeersch2, Mikel Mendicute1*, Lorena Martínez1 and Iker Sobrón1

Author Affiliations

1 Electronics and Computing Department, University of Mondragon, Loramendi, 4, 20500 Mondragon, Spain

2 Chalmers University of Technology, SE-412 96 Gothenburg, Sweden

For all author emails, please log on.

EURASIP Journal on Wireless Communications and Networking 2012, 2012:331  doi:10.1186/1687-1499-2012-331


The electronic version of this article is the complete one and can be found online at: http://jwcn.eurasipjournals.com/content/2012/1/331


Received:19 September 2011
Accepted:23 September 2012
Published:31 October 2012

© 2012 Ochandiano et al; licensee Springer.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

In this article, a novel detector design is proposed for orthogonal frequency division multiplexing (OFDM) systems over frequency selective and time varying channels. Namely, we focus on systems with large OFDM symbol lengths where design and complexity constraints have to be taken into account and many of the existing ICI reduction techniques can not be applied. We propose a factor graph (FG) based approach for maximum a posteriori (MAP) symbol detection which exploits the frequency diversity introduced by the ICI in the OFDM symbol. The proposed algorithm provides high diversity orders allowing to outperform the free-ICI performance in high-mobility scenarios with an inherent parallel structure suitable for large OFDM block sizes. The performance of the mentioned near-optimal detection strategy is analyzed over a general bit-interleaved coded modulation (BICM) system applying low-density parity-check (LDPC) codes. The inclusion of pilot symbols is also considered in order to analyze how they assist the detection process.

Keywords:
Factor graph (FG); Sum-product algorithm (SPA); OFDM Doppler; Intercarrier interference (ICI)

Introduction

Orthogonal frequency-division multiplexing (OFDM) is one of the most popular modulation schemes suitable for high-data rate wireless communication systems. The fact that the OFDM symbol period is much longer than a data symbol in a single-carrier system, makes it robust against intersymbol interference (ISI) caused by frequency-selective channels. For this reason it has been adopted by most of the last-generation wireless communication systems such as, IEEE’s 802.16 framework, the long-term evolution (LTE) project, or the second generation digital video broadcasting (DVB) specifications. Mobility support is one of the key features of these communication systems, which try to deal with the challenge of enabling mobile broadband services at high vehicular speed. Either the second generation DVB standards, with both terrestrial (DVB-T2) and handheld (DVB-NGH) versions, or the IEEE 802.16m are good examples of the mobility requirement of new wireless communication standards.

For conventional OFDM receivers, it is assumed that the channel remains static in an OFDM symbol period. In this case, equalization can be drastically simplified by turning the frequency selective channel into several parallel flat-fading channels. However, when the channel varies within an OFDM block, the subcarriers are not longer orthogonal and the system performance drops down severely. In fact, it is well known that one of the main drawbacks of OFDM is its susceptibility to the loss of orthogonality amongst subcarriers due to Doppler frequency shifts (mobile reception) or oscillator offsets, which lead to the so-called intercarrier interference (ICI). However, since it is shown that ICI can be modeled as a near-Gaussian random process that leads to an error floor [1-3], the use of advanced forward error correction (FEC) techniques such as low-density parity-check codes (LDPC) or turbo codes allows error-free reception at high-mobility scenarios without any particular reception strategy for ICI compensation.

With the aim of increasing spectral efficiency or enhancing the robustness in single-frequency networks (SFN), there is a growing trend toward using large OFDM symbol lengths (which means densely spaced subcarriers) in most of the new wireless communication systems (e.g., 32 K subcarriers in DVB-T2). As the delay spread increases, symbol duration must also be increased in order to maintain a near-flat channel in every subcarrier. Nevertheless, since the ICI power depends on the subcarrier spacing, the effect of time-varying channels becomes critical when large OFDM symbols are used and it is necessary to develop appropriate signal processing techniques to combat the mobility-induced ICI problem. Hence, from an ICI cancelation perspective, the most challenging scenario is a high-mobility scenario employing a large FFT mode.

This article focuses on systems with very large OFDM symbol lengths, and proposes a novel FG approach for ICI compensation based on the Forney observation model [4]. The proposed detector approaches the free-ICI case in high mobility scenarios and presents an intrinsic parallel structure suitable for high speed detection. We also incorporate the turbo principle applied to numerous fields [5], by which means the extrinsic information generated by the decoder is used as a priori information in the detection process. In fact, the receiver becomes a triple iterative scheme where iterations are performed in the detector, in the decoder and between both of them. Moreover, when LDPC codes are employed, a combined FG can be drawn which accomplishes both the detection and decoding process in a fully parallel manner [6], leading to potential hardware or DSP code saving.

The rest of this article is organized as follows. Section 2 gives a brief overview of the state of the art. Section 3 describes the system model in frequency selective channels affected by ICI. Section 4 reviews the principles behind the BP detection both for Forney and Ungerboeck approaches and proposes how to deal with pilot tones in the FG. Section 5 faces up to the complexity analysis. Section 6 analyzes the turbo reception scheme. Section 7 evaluates the proposed detector performance for several mobile scenarios, and finally, Section 8 summarizes the main conclusions of the article.

Notation: We denote the kth element of s by sk, and sk,l the element in the kth row and lth column of the matrix S. The conjugate transpose is denoted as (·)Hand ℜ{·} refers to the real part of a complex number. IN is the identity matrix and ∝ operator means that two terms are equivalent up to irrelevant additive and multiplicative values.

Background on ICI suppressing

In recent years, a wide range of OFDM receiver designs have been proposed in the literature for time-varying multipath channels. We can make a difference between linear equalizers based on zero-forcing (ZF) [7] or minimum mean-square error criterion [8-10], and non-linear equalizers based on ICI cancelation or decision-feedback [11-16]. Linear equalization requires the inversion of the frequency channel matrix, which is prohibitively complex for large OFDM symbols. Although many approaches have been proposed for solving this problem [7,8], in general terms, they show poor performance at high-mobility scenarios.

Non-linear equalizers outperform linear approaches at the cost of higher complexity. Since there are similarities between the ICI and the multiuser interference (MUI), well-known interference cancelation techniques such as, successive interference cancelation (SIC) and parallel interference cancelation (PIC) can be directly applied to ICI suppressing. While the former present higher computational complexity and processing delays, the latter suffer from performance loss at high Doppler frequencies. In [17], a combination of the MMSE method and PIC is presented, which outperforms the conventional PIC schemes. Recently, space alternating generalized expectation maximization (SAGE) has been proposed for ICI cancelation [18,19].

Other research studies explore optimal or near-optimal solutions, e.g., near-maximum likelihood (ML) approaches [20,21], sphere decoding (SD) [22], or maximum a posteriori (MAP) equalization employing Bahl Cocke Jelinek Raviv (BCJR) algorithm [23]. Apart from the frequency domain, time domain solutions have also been proposed [24-26].

Most of the techniques mentioned above, either do not show good performance in high mobility scenarios or they are not suited for large FFT modes in terms of complexity and design constraints. Peng and Ryan suggest a parallel implementation of a two-stage equalizer [27] suitable for large FFT sizes, which has low complexity comparing with other proposals when the OFDM symbol length is very large. Recently, a time domain per sub-block equalizer has been presented [28] to tackle the same problem formulated in this article. One of the drawbacks of this proposal is that, pilot tones, used in most of wireless communication systems can not be used for channel estimation.

Factor graphs (FG) have received the attention of many researchers and have been applied to solve many technical problems in communications [29,30]. Essentially, FGs represent graphically the factorization of a function, and have been shown to be a good alternative to solve complex inference problems. Many recent articles address the use of FG-based algorithms (also known as belief propagation (BP) algorithms) for detection or equalization in channels with memory [31-33]. Recently, this framework has been extended to ICI channels: a set of detection algorithms have been presented based on the Ungerboeck approach for MAP symbol detection [34], while in [35], progressive parallel ICI cancelation has been proposed based on FGs.

System model

We consider a single user bit interleaved coded modulation (BICM)-OFDM system over a frequency selective time varying channel. OFDM modulation divides a data stream with sample period T into N substreams, which are modulated on different subcarriers, equally spaced at a distance of Δf = 1/(TN). N coded symbols s = [s1,…,sN]T, equally likely and chosen from a complex constellation χare transformed into the time domain applying the inverse discrete Fourier transform (IDFT) rule. A cyclic prefix (CP) of length G is added at the output of the IDFT block. No interblock interference (IBI) is considered since it is assumed that the CP is longer than the channel delay spread.

According to the COST 207 model, the channel taps are considered independent random processes with Rayleigh statistics, and the Doppler frequency Fd is modeled with the widely used Jakes’ Doppler spectrum. The mobile scenarios considered in this article are characterized by the normalized Doppler frequency fd = FdΔf. The received signal after CP removal can be expressed in matrix form as [8]

r = H ¯ F H s + z , (1)

where H ¯ is a N×N time-domain channel matrix, Fstands for the standard N-point DFT matrix, and vectors r and z are the received signal and the additive white Gaussian noise (AWGN) discrete samples, respectively. We assume that z is a complex Gaussian noise vector with zero mean and covariance E zz H = σ n 2 I N . Applying the discrete Fourier transform (DFT) leads to the received signal in frequency domain

y = F H ¯ F H s + Fz = Hs + z f , (2)

where H = F H ¯ F H is the frequency-domain channel matrix, and zf = Fzis the frequency domain noise vector which has the same statistics as time domain noise z. Note that y = Fr.

In time invariant channels (fd = 0), H ¯ is circulant and therefore Hbecomes diagonal, leading to the low complexity per subcarrier equalization that characterizes OFDM systems. Nevertheless, when the channel is time variant, the off-diagonal elements of H are not longer zero, and the ICI comes up making the implementation of conventional equalizers prohibitively complex. The ICI power PICI is expressed by [2]

P ICI = N 1 N 2 N 2 k = 1 N 1 J 0 2 Π f d N N k , (3)

where J0(·) is the zero-order Bessel function of the first kind. From (3) it is clear that the ICI power is a function of the Doppler frequency, the subcarrier spacing and the OFDM block length, and does not depend on the signal constellation.

Since most of the energy in His concentrated around the main diagonal, the matrix can be approximated by its banded version, where q defines how many off-diagonals above and below the main diagonal are taken into consideration (Figure 1). For example, more than 64 % of the energy clusters in the neighboring three subcarriers (q = 1) for fd = 0.13. As a result of increasing the vehicular speed, or the OFDM block size, more and more energy leaks into the off-diagonal elements of H, and it is necessary to define a proper value of q to get a trade-off between performance and computational complexity. For many applications, it is shown that q = 1 is enough to achieve good BER performance in relatively high mobility scenarios. In this article, q = 1 is further assumed.

thumbnailFigure 1. Frequency-domain system input-output relation after removing the CP.

The proposed BP detector’s performance is evaluated considering the 64800 bits-long LDPC code adopted by DVB [36] as channel code. Ideal channel state information (CSI) and perfect synchronization are assumed at the receiver. Nevertheless, we consider the presence of pilot carriers in order to describe how they can assist the detection process. The PP1 DVB-T2 pilot pattern [36] has been adopted for this purpose. The structure of the receiver is shown in Figure 2.

thumbnailFigure 2. Block diagram of the BICM iterative receiver chain including the proposed BP detector.

Factor graph based approach for ICI mitigation

Factor graph & sum-product algorithm

A FG is a bipartite graph G = {V,E} which represents the factorization of a complex function into the product of local functions, where V represents the set of vertices of the graph and E is the set of edges. Let s = {s1,…,sD} represent a set of variables and f (s) a multivariate function. We consider the factorization of f (s) into K factors:

f s = k = 1 K f k ψ k , (4)

where ψk⊆ {s1,…,sD} is the kth variable subset, and fk(·) is a real-valued function. Variable nodes sk are associated with vertex called variable nodes, factors fkare represented with vertex called function nodes and there is an edge e = (sn,fk) ∈ E for a factor fk and every variable snϑ (fk), being ϑ(v) the set of neighbors of a given node v.

Assuming that fk(ψk) is a probability mass function (pmf), and the graph has no cycles, the application of the sum-product algorithm (SPA) to the FG leads to the exact computation of the marginal pmfs, avoiding the cumbersome evaluation of the joint pmf. Denoting by μ s n f k s n the message sent from a variable node snto a factor node fk(ψk), and by μ f k s n s n the message in the opposite direction (log domain is assumed), the message computation in variable and factor nodes, respectively, becomes

μ s n f k s n = h ϑ s n f k μ h s n s n , (5)

μ f k s n s n = ln s n exp ln f k ( ψ ) + z ϑ f k s n μ z f k ( z ) , (6)

where s n indicates the sum over all variables excluding sn(summary operation). After a fix number of iterations between variable and factor nodes, marginal a posteriori pmfs P (sn|y) are computed, necessary to evaluate the MAP symbol detection rule given by

ŝ n = arg max s n P ( s n | y ) , (7)

as follows

P s n | y = exp h ϑ s n μ h s n s n exp s n χ h ϑ s n μ h s n s n . (8)

The order in which the messages are updated in the SPA is defined by the message-passing schedule. Different schedules can be considered, and the performance of the algorithm varies depending on it. In this article, in order to parallelize the algorithm as much as possible, the so-called flooding schedule is adopted: all variable and function nodes pass new messages to their neighbors in each iteration.

If the FG is not cycle-free, the algorithm does not converge to the exact marginal pmfs and it is hard in general to state whether it converges to reasonable distributions or not. Consequently, approximate marginal a posteriori distributions are computed. However, for many applications, it has been shown that loopy FGs provide excellent empirical results and turn out to be a feasible alternative for designing near optimal algorithms [37].

MAP symbol detection based on Forney approach

In this article, we are interested in evaluating the joint a posteriori pmf of the transmitted information sequence S conditioned on an observation sequence y(Figure 1). Based on the Bayesian estimation theory, it is well known that the transmission of a random signal over an AWGN-affected random channel can be modeled as

P ( s | y ) P ( y | s ) P ( s ) . (9)

For the system model assumed in this article (2), and assuming q = 1, the ψ variable subset described in (4) is defined as ψk = {sk−1sksk + 1}. Based on the Forney observation model [4], the likelihood function P(y|s) is factorized as follows

P y | s k = 1 N f k s k 1 , s k , s k + 1 , (10)

where implicitly sk = 0 for k ≤ 0 and k > N and

f k s k 1 , s k , s k + 1 exp y k n = k 1 k + 1 H kn s n 2 σ n 2 . (11)

Part of the corresponding FG is depicted in Figure 3, where it is clear that the shortest cycle in the graph is length 4 (the graph has girth 4). Although it has been shown that length 4 cycles degrade severely the performance of the algorithm, it is known that in particular cases they do not pose significant problems in terms of BER. Throughout this article, the performance of the proposed BP algorithm is compared with the optimal forward backward (FB) algorithm with the aim of determining the performance loss due to the presence of cycles in the graph.

thumbnailFigure 3. Factor Graph for the Forney approach-based MAP symbol detection.

MAP symbol detection based on Ungerboeck approach

Two approaches have been proposed for MAP sequence detection until now: Forney approach [4] and the Ungerboeck approach [38]. The former uses the output of a whitened match filter for computing the branch metrics of the Viterbi algorithm (VA) leading to a Markovian channel model, while the latter can work directly with the output of a filter matched to the received pulse. Both strategies are equivalent and only differ in the expression of the VA branch metrics. Nevertheless, when the MAP symbol detection strategy is adopted, while the Forney observation model can be directly applied, the extension of the Ungerboeck observation model is not trivial. In [39], this problem is solved and it is proven that both models are also equivalent for MAP symbol detection. Based on the keen observation made in [39], a BP detection algorithm based on the Ungerboeck approach [40] was proposed for dispersive channels providing impressive complexity reduction comparing with [33]. Similar low-complexity FGs have been derived for channels with memory in [32]. Recently, a novel BP detector employing the Ungerboeck observation model has been presented [34] for ICI channels which shows a good performance-complexity trade-off at low Doppler frequencies.

Based on the Ungerboeck observation model the likelihood function P(y|s) is factorized as follows

P ( y | s ) exp y H s 2 σ n 2 exp 2 s H x s H G s σ n 2 , (12)

where x = HHy and G is the Hermitian matrix defined as G = HHH. The same factorization can be expressed in scalar form as

P ( y | s ) k = 1 N exp 1 σ n 2 2 s k x k 1 2 g k , k s k 2 i = 1 2 q g k , k i s k s k i . (13)

By defining the next functions

f k s k = exp 1 σ n 2 2 s k x k 1 2 g k , k s k 2 , (14)

u k , i s k , s k i = exp 1 σ n 2 2 g k , k i s k s k i , (15)

the likelihood function can be expressed in a compact form as

P y | s k = 1 N f k s k i = 1 2 q u k , i s k , s k i , (16)

where implicitly sk = 0 for k ≤ 0 and k > N. Part of the FG derived from (16) is depicted in Figure 4, where it can be appreciated that has girth 6. This fact seems indicate that the algorithm based on the Ungerboeck approach may perform better than the algorithm based on the Forney observation model (Figure 3). However, numerical results presented in Section 5 show that for high Doppler frequencies, the performance of the Ungerboeck approach drops down severely. Note from (16) that there is still a fundamental complexity difference with respect to the factorization in (11). The complexity of both approaches will be further analyzed in Section 5.

thumbnailFigure 4. Factor Graph for the Ungerboeck approach-based MAP symbol detection.

Pilot-assisted graph-based detection

Many wireless communication systems use training data sequences for channel estimation. For instance, DVB standards consider different pilot patterns in order to adjust the transmission requirements to different channel conditions. Due to the frequency correlation introduced by the ICI, BP detection has to be performed before pilot removal. This fact makes the FG longer since variable and function nodes corresponding to pilot carriers have to be taken into account, but, at the same time, there is known information that we can use in order to assist the detection process.

Figure 5 shows how to deal with pilot carriers in the message passing process, where sn and fk are the variable and function nodes corresponding to a given pilot symbol, respectively. On one hand, outgoing messages from fkare more accurate representations of marginal distributions since known information corresponding to snhas been used in (11). On the other hand, messages transmitted by variable node sn are uniform distributions of the corresponding pilot symbol, which are useless for data symbol detection, thus, we can consider variable node sn as an idle node and disconnect from their neighbors. As we can appreciate in Figure 5, as a consequence of disconnecting variable nodes corresponding to pilot symbols, some cycles are straightforwardly removed from the FG. As an example, considering PP1 pilot pattern of DVB-T2 where pilot carriers represent the 8% of the total carriers, 8% of length 4 cycles are directly removed, which yields an additional performance improvement.

thumbnailFigure 5. Pilot processing in the factor graph.

Complexity analysis

As it has been mentioned before, the implementation of the SPA is accomplished in log domain. This can be done using the Jacobian logarithm which only requires additions and the evaluation of a non-linear function.

The complexity of the proposed detection algorithm is mainly a function of the following parameters: constellation size M, number of total subcarriers N and the bandwidth parameter q. Table 1 describes the complexity of different algorithms considered in this article. As it can be appreciated, it is obvious the weakness of the proposed algorithm comparing with the Ungerboeck approach based algorithm presented in [34] in terms of complexity, since the Forney approach leads to exponential dependence on the bandwidth parameter q.

Table 1. Complexity analysis

We can adopt different strategies to reduce the complexity of the algorithm proposed in this article. Compared with FB algorithms, the complexity reduction is much easier in BP algorithms due to there are no constraints imposed by the Trellis structure. Reduced state sequence detection (RSSD) [41,42] can be an effective method for reducing the complexity with a minor performance loss.

Turbo reception

The turbo principle has been applied to numerous fields beyond the practical decoding of powerful channel codes such as turbo codes and LDPC codes. In general terms, it consists of iterative exchanging of soft information between different blocks in a communications receiver with the aim of improving the system performance. While the original turbo receiver is conceived as two independent processors exchanging extrinsic soft information [5] (classical approach), there is a generalized view of the turbo approach derived from [6] which is described as passing soft decision messages on cyclic graphical models (graphical approach).

Considering that LDPC codes are decoded using BP algorithms, the joint data detection and decoding process can be depicted by a higher-order FG following the recent trend derived from the graphical approach. A triple iterative scheme is proposed where messages are exchanged among internal LDPC decoder nodes (LDPC iterations), internal detector nodes (BP iterations), and between the detector and the decoder (turbo iterations) (Figure 6).

thumbnailFigure 6. High order factor graph performing joint data detection and decoding following the graphical approach for turbo reception.

In order to reach a good balance between the latency and system performance, we are strongly interested in analyzing the convergence behavior of the proposed receiver and finding out the best turbo schedule, i.e., how many BP and LDPC iterations perform in each turbo iteration. In this article, this task is faced up by extensive simulation results.

Numerical results

In this section, numerical results show the performance of the proposed BP detector for several mobile scenarios. The next table summarizes the main simulation parameters.

The mobile scenarios considered in this article are characterized by the normalized Doppler frequency. For the simulation parameters described in Table 2, fd = 0.13 corresponds to about 200 km/h of vehicular speed and fd = 0.5 represents the same vehicular speed for the 32 K OFDM mode. 6-taps typical urban (TU6) channel and 6-taps rural area (RA6) channel are used for modeling urban and rural environments, respectively. The channels of interest in this article correspond to Doppler frequencies beyond fd = 0.3.

Table 2. Simulation parameters

Forney approach versus Ungerboeck approach

The performance of the proposed detector is compared with the Ungerboeck approach-based BP detector presented in [34]. TU6 channel model has been simulated for QPSK constellation. Five iterations have been performed in the detector in both cases, and no turbo reception (no turbo iterations) is considered for simplicity. Min-sum (MS) algorithm is adopted for LDPC decoding with code rate (CR) 2/3. Both BP algorithms are compared with the optimal FB algorithm in order to quantize the performance loss of both BP algorithms with respect to the optimal MAP symbol detection, although it is well known that FB algorithms can not be considered for large FFT sizes due to its intrinsic serial structure.

As it is shown in Figure 7, both approaches show almost the same performance for low Doppler frequencies and they are very closed to the optimal FB. Consequently, and taking into account that the Ungerboeck approach is less complex than the proposed Forney approach, it is clear that for low mobile scenarios, the best option is the Ungerboeck approach-based BP detector [34]. However, it turns out that the performance of the Ungerboeck approach-based BP detector drops down severely at high-mobility scenarios. In Figure 8, it is shown that the Ungerboeck approach tends to an error-floor whereas the Forney approach presents a performance loss of about 0.6 dB at BER = 10−4 with respect to the optimal FB. A detailed observation of the messages into the FG at asymptotic regime (high SNR and high Doppler frequency) indicates that the factorization based on the Ungerboeck observation model (13) collapses when in the presence of very reliable input information factor nodes output uniform distributions. Therefore, the Ungerboeck approach is not a good option for the channels of interest in this article.

thumbnailFigure 7. Ungerboeck approach versus Forney approach for fd = 0.16 with QPSK.

thumbnailFigure 8. Ungerboeck approach versus Forney approach for fd  = 0.4. TU6 channel and QPSK constellation are considered.

On the other hand, we can notice that the performance loss in comparison with the optimal FB is higher for high Doppler frequencies, but even for the highest Doppler frequency considered in this section (fd = 0.4), our proposal is closed to the optimal one. Hence, it is shown that, in this case, short cycles do not prevent from an acceptable performance of the SPA and the BP algorithm provides good ICI compensation.

BP detection versus MAP detection with ICI cancelation

The proposed detector is also compared with the MAP ICI canceler proposed in [27]. This detector presents a fully parallel architecture which accommodates hardware implementation and its complexity is much less than other proposals’ in the literature when the OFDM symbol length is very large. As a result, it consists of a good candidate for ICI compensation when very large OFDM symbol lengths are used. In general terms, it is composed of two stages: the first one estimates the transmitted data by means of a Viterbi-like algorithm and the second one performs the MAP detection suppressing the ICI previously reconstructed with the output of stage one. TU6 channel model and QPSK constellation is adopted, 5 BP iterations have been performed in the detector, and again no turbo reception (no turbo iterations) is considered for simplicity.

As it is shown in Figure 9, both schemes are very closed to each other for fd = 0.13. Nevertheless, the estimates produced by the first stage of the MAP ICI canceler degrade severely at high Doppler frequencies, thus it suffers from a performance loss which leads to an error-floor. Hence, although the MAP ICI canceler proposed in [27] is well suited to large FFT sizes from an implementation point of view, it is clearly shown that the proposed BP detector outperforms the MAP ICI canceler in terms of BER performance for the channels of interest in this article.

thumbnailFigure 9. MAP ICI canceler [27] vs BP detector for fd  = 0.13 and 0.3. TU6 channel and QPSK constellation are considered. MS algorithm is adopted for LDPC decoding with CR= 2/3.

BER performance analysis for the graphical approach

The performance of the proposed BP detector is evaluated for QPSK and 16QAM constellations and over two channel models. Turbo reception based on the graphical approach presented in Section 6 is considered and 2 BP iterations and 20 LDPC iterations are performed in each of the 5 turbo iterations. MS algorithm and CR= 2/3 are considered for LDPC decoding in QPSK simulations and SPA and CR= 1/2 in 16QAM simulations.

Figure 10 shows the system performance for the TU6 channel model. Regarding the QPSK constellation, it can be appreciated that the BP detector is able to remove the error-floor caused by the ICI even for the highest Doppler frequency considered in this article (fd = 0.5). We can also appreciate that the BER curve shifts to the right side as the Doppler frequency increases. This is caused by the residual ICI due to not considering the bandwidth parameter q higher that 1, which would enhance the BER performance with the cost of much higher complexity. It is considered that q = 1 gives a good trade-off between performance/complexity, since the performance loss comparing with the free-ICI case is about 0.8 dB for fd = 0.5. In the case of 16QAM constellation, the mentioned residual ICI, which can be modeled as a near-Gaussian noise, affects much more to the signal preventing from achieving the free-ICI performance. Still, the BP detector is able to remove the error-floor caused by the Doppler frequency and approach the free-ICI case up to 1.2 dB.

thumbnailFigure 10. BER performance of the proposed BP detector over TU6 channel model. Turbo reception based on the graphical approach is adopted. QPSK and 16QAM constellations are considered.

Figure 11 shows the performance of the BP detector for the RA6 channel, where it is shown that the behavior of the receiver changes comparing to the previous case: the BER performance improves significantly as the Doppler frequency increases. The reason behind this behavior is that the maximum frequency diversity order that a BICM-OFDM system can achieve is defined by min(dfreeL) [43], where dfreerefers to the Hamming distance of the channel code and L stands for the number of channel taps. Since the Hamming distance of the LDPC code is very high, the diversity order of the BICM-OFDM system is determined by the total number of channel taps, which is very low for RA6 (6-taps) comparing with TU6 (47-taps).

thumbnailFigure 11. BER performance of the proposed BP detector over RA6 channel model. Turbo reception based on the graphical approach is adopted. QPSK and 16QAM constellations are considered.

In the case of RA6 channel, the exploitation of the frequency diversity introduced by the ICI in the BP detector increases the diversity order of the system becoming the BER curve steeper as the Doppler frequency increases. Taking into account that high mobility scenarios usually take place in rural environments, it is concluded that the combination of the BICM-OFDM scheme with the proposed BP detector consists of an excellent solution for high mobility scenarios. For instance, it is shown that, the BP detector achieves free-ICI performance for fd = 0.22 with 16QAM constellation.

Conclusions

A novel FG-based detector has been proposed for wireless mobile OFDM systems which enables exploitation of the frequency diversity available in the received signal affected by ICI. The implementation of this novel scheme for signal detection allows high speed joint detection and decoding providing excellent performance at high mobility scenarios comparing with other ICI compensation schemes in the literature. The performance of the proposed detector has been analyzed over TU6 and RA6 channel models using LDPC codes. Furthermore, the inclusion of pilot symbols has also been considered in order to analyze how they can assist the detection process. Compared to the traditional FB receiver, the proposed BP detector is better suited to large FFT sizes.

Competing interests

The authors declare that they have no competing interests.

Acknowledgements

The authors want to thank Department of Industry and Innovation and Department of Education, Universities and Research of the Basque Government, for its support and funding through IKERTU programme.

References

  1. Y Li, LJ Cimini, Bounds on the interchannel interference of OFDM in time-varying impairments. IEEE Trans. Commun 49, 401–404 (2001). Publisher Full Text OpenURL

  2. M Russell, G Stuber, Interchannel interference analysis of OFDM in a mobile environment, IEEE 45th Vehicular Technology Conference,, 820–824 (1995)

  3. T Wang, J Proakis, E Masry, J Zeidler, Performance degradation of OFDM systems due to Doppler spreading. IEEE Trans. Wirel. Commun 5, 1422–1432 (2006)

  4. GD Forney, Lower bounds on error probability in the presence of large intersymbol interference. IEEE Trans. Commun 20, 76–77 (1972). Publisher Full Text OpenURL

  5. M Tuchler, R Koetter, A Singer, Turbo equalization: principles and new results. IEEE Trans. Commun 50, 754–767 (2002). Publisher Full Text OpenURL

  6. N Wiberg, Codes and decoding on general graphs (Ph, 1996), . D. thesis, Linkoping Univ., Sweden

  7. CY Hsu, WR Wu, Low-complexity ICI mitigation methods for high-mobility SISO/MIMO-OFDM systems. IEEE Trans. Veh. Technol 58, 2755–68 (2009)

  8. K Fang, L Rugini, G Leus, Low-complexity block turbo equalization for OFDM systems in time-varying channels. IEEE Trans. Signal Process 56, 5555–5566 (2008)

  9. X Huang, HC Wu, Robust and efficient intercarrier interference mitigation for OFDM systems in time-varying fading channels. IEEE Trans. Veh. Technol 56, 2517–2528 (2007)

  10. L Rugini, P Banelli, G Leus, Low-complexity banded equalizers for OFDM systems in Doppler spread channels. EURASIP. J. Appl. Signal Process 2006, 13 (2006)

  11. X Cai, G Giannakis, Bounding performance and suppressing intercarrier interference in wireless mobile OFDM. IEEE Trans. Commun 51, 2047–2056 (2003). Publisher Full Text OpenURL

  12. WS Hou, BS Chen, ICI cancellation for OFDM communication systems in time-varying multipath fading channels. IEEE Trans. Wirel. Commun 4, 2100–2110 (2005)

  13. X Li, R Zhou, V Chakravarthy, S Hong, Z Wu, Total intercarrier interference cancellation for OFDM mobile communication systems, IEEE Consumer Communications and Networking Conference, ((Piscataway, NJ, 2010), pp. 1–5

  14. G Taubock, M Hampejs, P Svac, G Matz, F Hlawatsch, K Grochenig, Low-complexity ICI/ISI equalization in doubly dispersive multicarrier systems using a decision-feedback LSQR algorithm. IEEE Trans. Signal Process 59, 2432–2436 (2011)

  15. S Tomasin, A Gorokhov, H Yang, JP Linnartz, Iterative interference cancellation and channel estimation for mobile OFDM. IEEE Trans. Wirel. Commun 4, 238–245 (2005)

  16. Y Zhao, SG Haggman, Intercarrier interference self-cancellation scheme for OFDM mobile communication systems. IEEE Trans. Commun 49, 1185–1191 (2001). Publisher Full Text OpenURL

  17. SU Hwang, JH Lee, J Seo, Low complexity iterative ICI cancellation and equalization for OFDM systems over doubly selective channels. IEEE Trans. Broadcast 55, 132–139 (2009)

  18. E Panayirci, H Senol, HV Poor, Joint channel estimation, equalization, and data detection for OFDM systems in the presence of very high mobility. IEEE Trans. Signal Process 58, 4225–4238 (2010)

  19. H Wu, S Yang, J Ou, L Yang, Improved ICI mitigation scheme over time-varying channels for high-mobility OFDM systems. J. Converg. Inf. Technol 6, 264–272 (2011)

  20. ML Ku, WC Chen, CC Huang, EM-based iterative receivers for OFDM and BICM/OFDM systems in doubly selective Channels. IEEE Trans. Wirel. Commun 10, 1405–1415 (2011)

  21. S Ohno, Maximum likelihood inter-carrier interference suppression for wireless OFDM with null subcarriers. IEEE International Conference on Acoustics, Speech, and Signal Processing,, 849–852 (2005)

  22. Y Kou, WS Lu, A Antoniou, Application of sphere decoding in intercarrier-interference reduction for OFDM systems. IEEE Pacific Rim Conference on Communications, Computers and Signal Processing, 360–363 (2005)

  23. D Liu, M Fitz, Iterative MAP equalization and decoding in wireless mobile coded OFDM. IEEE Trans. Commun 57, 2042–51 (2009)

  24. E Peiker, WG Teich, J Lindner, Windowing in the receiver for OFDM systems in high-mobility scenarios. International Conference on Multimedia Communications, Services and Security, 57–65 (2009)

  25. P Schniter, Low-complexity equalization of OFDM in doubly selective channels. IEEE Trans. Signal Process 52, 1002–11 (2004). Publisher Full Text OpenURL

  26. A Stamoulis, S Diggavi, N Al-Dhahir, Intercarrier interference in MIMO OFDM. IEEE Trans. Signal Process 50, 2451–64 (2002). Publisher Full Text OpenURL

  27. F Peng, W Ryan, A low-complexity soft demapper for, OFDM fading channels with ICI. IEEE Wireless Communications and Networking Conference,, 1549–1554 (2006)

  28. P Baracca, S Tomasin, L Vangelista, N Benvenuto, A Morello, Per sub-block equalization of very long ofdm blocks in mobile communications. IEEE Trans. Commun 59, 363–368 (2011)

  29. F Kschischang, B Frey, HA Loeliger, Factor graphs and the sum-product algorithm. IEEE Trans. Inf. Theory 47, 498–519 (2001). Publisher Full Text OpenURL

  30. H Wymeersch, Iterative Receiver Design (Cambridge University Press, Cambridge, 2007)

  31. A Anastasopoulos, KM Chugg, G Colavolpe, G Ferrari, R Raheli, Iterative detection for channels with memory. Proc. IEEE 95, 1272–1294 (2007)

  32. G Colavolpe, On LDPC codes over channels with memory. IEEE Trans. Wirel. Commun 5, 1757–66 (2006)

  33. G Colavolpe, G Germi, On the application of factor graphs and the sum-product algorithm to ISI channels. IEEE Trans. Commun 53, 818–825 (2005)

  34. W Haselmayr, B Etzlinger, A Springer, Factor-graph-based detection algorithms for coded OFDM over time-varying channels. JNCW 2011 - NEWCOM++/COST 2100 joint Workshop, (2011)

  35. CW Huang, PA Ting, CC Huang, A novel message passing based MIMO-OFDM data detector with a progressive parallel ICI canceller. IEEE Trans. Wirel. Commun 10, 1260–1268 (2011)

  36. ETSI: Digital Video Broadcasting (DVB); Frame structure channel coding and modulation for a second generation digital terrestrial television broadcasting system (DVB-T2) (ETSI EN 302 755 V1), . 1.1 (2009)

  37. M Kaynak, T Duman, E Kurtas, Belief propagation over MIMO frequency selective fading channels. Joint International Conference on Autonomic and Autonomous Systems and International Conference on Networking and Services (ICAS/ICNS), 45–45 (2005)

  38. G Ungerboeck, Adaptive maximum-likelihood receiver for carrier-modulated data-transmission systems. IEEE Trans. Commun 22, 624–36 (1974)

  39. G Colavolpe, A Barbieri, On MAP symbol detection for ISI channels using the Ungerboeck observation model. IEEE Commun. Lett 9, 720–722 (2005). Publisher Full Text OpenURL

  40. D Fertonani, A Barbieri, G Colavolpe, Novel graph-based algorithms for soft-output detection over dispersive channels. IEEE Global Telecommunications Conference,, 1–5 (2008)

  41. PR Chevillat, E Eleftheriou, Decoding of trellis-encoded signals in the presence of intersymbol interference and noise. IEEE Trans. Commun 37, 669–676 (1989). Publisher Full Text OpenURL

  42. M Eyuboglu, S Qureshi, Reduced-state sequence estimation with set partitioning and decision feedback. IEEE Trans. Commun 36, 13–20 (1988). Publisher Full Text OpenURL

  43. E May, E Ayanoglu, Full frequency diversity codes for single input single output systems. IEEE Vehicular Technology Conference,, 1870–1874 (2004)