Abstract

single-carrier (SC) modulation with frequency-domain equalization (FDE) combined with iterative (turbo) FDE schemes has been shown to be suitable for cyclic-prefix (CP)-assisted block transmission within broadband wireless systems. Its excellent performance in severely time-dispersive channels, makes it a promising candidate for future broadband wireless systems. In this paper we investigated the impact of the correlation factor estimation in the performance of iterative block decision feedback equalization (IB-DFE) receivers. Since the correlation factor represents a key parameter to ensure the good performance of these receivers, reliable estimates of this factor are needed in the feedback loop. We present several methods to estimate the correlation coefficient. We also propose a technique to compensate the inaccuracy of the estimation.

1. Introduction

Future wireless systems are required to support high quality of service at high data rates, for which we can have severely time-dispersion effects associated to the multipath propagation. It is known that block transmission techniques, with appropriate cyclic extensions and employing FDE techniques, are excellent candidates for severely time-dispersive channels [1, 2]. The most popular techniques based on this concept are orthogonal frequency division multiplexing (OFDM) and single-carrier with frequency-domain equalization (SC-FDE). However, the OFDM signals have high envelope fluctuations and a high peak-to-mean envelope power ratio (PMEPR) leading to amplification difficulties. An alternative to OFDM modulations are SC modulations using FDE and CP-assisted block transmission techniques, especially for the uplink of broadband wireless systems [1, 2]. The performances and overall implementation complexities are similar to SC-FDE and OFDM [1, 2].

Typically the receiver for SC-FDE schemes is a linear FDE. However, it is known that nonlinear equalizers outperform linear equalizers [3]. iterative block decision feedback equalizer (IB-DFE) [4] is a promising iterative FDE technique for SC-FDE that was first proposed in [5] and extended to diversity scenarios [6] and layered space-time schemes [7]. IB-DFE receivers can be regarded as iterative DFE receivers with the feedforward and the feedback operations implemented in the frequency domain. Since the feedback loop takes into account not just the hard decisions for each block, but also the overall block reliability we have small error propagation. Consequently, the IB-DFE techniques offer much better performances than the noniterative methods, with performances that can be close to the matched filter bound (MFB) [8] and can be regarded as low complexity turbo equalization schemes [9, 10], since the feedback loop uses the equalizer outputs instead of the channel decoder outputs. True turbo FDE schemes can also be designed based on the IB-DFE concept [11, 12].

In this paper, we study the impact of the correlation coefficient estimation on the asymptotic performance of IB-DFE schemes. Several methods to estimate the correlation coefficient are derived as well as a proposed technique to compensate the inaccuracy of the estimations.

This paper is organized as follows: conventional and turbo IB-DFE receivers are described in Section 2. Various methods to estimate the correlation coefficient are presented in Section 3. The compensation technique is also included. A set of performance results is presented in Section 4, and Section 5 is concerned with the conclusion of this paper.

2. System Description

2.1. Transmitted and Received Signals

We consider an SC-FDE modulation where the data is transmitted in blocks of ??-useful modulation symbols {????;??=0,,??-1}, resulting from a direct mapping of the original data bits into a selected signal constellation, for example, QPSK. A cyclic prefix with length longer than the channel impulse response is posteriorly appended, resulting the transmitted signal {????;??=-????,,??-1}, where ???? denotes the number of guard samples. This signal is transmitted over a time-dispersive channel. The received signal is sampled at the receiver where the CP samples are removed, leading in the time domain the samples {????;??=0,,??-1}. The corresponding frequency-domain block, obtained after an appropriate size-?? discrete fourier transform (DFT) operation, is {????;??=0,1,,??-1}, where ???? can be written as ????=????????+????,(1) where ???? denotes the overall channel frequency response for the ??th frequency of the ??th time block, and ???? represents channel noise term in frequency-domain.

2.2. IB-DFE Basic Receiver Structure

For the sake of simplicity, we will assume in this section perfect synchronization conditions.

In SC-FDE to deal with the channel effects, we can employ a linear FDE. However, the performance can be improved if the linear FDE is replaced with an IB-DFE [5], whose basic structure is depicted in Figure 1. In this case, for the ??th iteration, the frequency-domain block at the output of the equalizer is {?????(??);??=0,1,,??-1}, with ?????(??)=????(??)????-????(??)?????(??-1),(2) where {????(??);??=0,1,,??-1} are the feedforward coefficients, and {????(??);??=0,1,,??-1} are the feedback coefficients. {?????(??-1);??=0,1,,??-1} is the DFT of the block {????(??-1);??=0,1,,??-1}, with ???? denoting the “hard decision” of ???? from the previous FDE iteration. It can be shown that the optimum coefficients ???? and ???? that maximize the overall SNR (instantaneous signal-to-noise ratio) associated to the samples ????? are [6]????(??)???=????(??)?????-1,(3)????(??)=????*????????+1-??(??-1)?2?||????||2;(4) respectively, where ??=??[|????|2]/??[|????|2] (common to all data blocks), and ?? selected to guarantee that ?1/????-1??=0????(??)????=1.

As we can see from (3) and (4), the correlation factor ??(??-1) is a key parameter for the good performance of IB-DFE receivers, since it supplies a blockwise reliability measure of the estimates, associated to the previous iteration, that are employed in the feedback loop. This is done in the feedback loop by taking into account the “hard-decisions” for each block plus the overall block reliability, which reduces error propagation problems. The correlation factor ??(??-1) is defined as ????(??-1)=???????(??-1)??*??????||????||2?=????????(??-1)??*??????||????||2?,(5) where the block {????(??-1);??=0,1,,??-1} denotes the data estimates associated to the previous iteration, that is, the hard decisions associated to the time-domain block at the output of the FDE, {????(??);??=0,1,,??-1} = IDFT {?????(??);??=0,1,,??-1}.

It should be noted that (2) can be written as ?????(??)=????(??)????-?????(??)??(??-1)??,Block,(6) where ?????(??)=????(??)/????(??-1) and ??(??-1)??,Block=????(??-1)?????(??-1) (as stated before, ??(??-1) can be considered as the blockwise reliability of the estimates {?????(??-1);??=0,1,,??-1}).

To improve the IB-DFE performance, it is possible to use “soft decisions,” ????(??), instead of “hard decisions,” ????(??). Under these conditions, the “blockwise average” is substituted by “symbol averages” [11]. This can be done by using {??(??-1)??,Symbol;??=0,1,,??-1} instead of {??(??-1)??,Block;??=0,1,,??-1} = DFT {??(??-1)??,Symbol;??=0,1,,??-1}, where ??(??-1)??,Symbol denotes the average symbol values conditioned to the FDE output from previous iteration, ????(??-1). To simplify the notation, we will replace ??(??-1)??,Symbol with ????(??-1) in the following equations.

For QPSK constellations, the conditional expectations associated to the data symbols for the ??th iteration are given by????(??)???=tanh????(??)2????+??tanh????(??)2?=????????????+??????????????,(7) with the loglikelihood ratios (LLRs) of the “in-phase bit” and the “quadrature bit,” associated to ?????? and ??????, given by??????(??)=2??2????????(??),??(8)????(??)=2??2????????(??),(9) respectively, with??2??=12???||????-????(??)||2?˜12????-1???=0||????(??)-????(??)||2(10) with the signs of ?????? and ?????? defining the hard decisions ??????=±1 and ??????=±1, respectively. In (7), ?????? and ?????? denote the reliabilities related to the “in-phase bit" and the “quadrature bit” of the ??th symbol and are given by??????(??)=???????????????????||??????||2?=|||||???tanh????(??)2?|||||,??????(??)=???????????????????||??????||2?=|||||???tanh????(??)2?|||||.(11) Therefore, the correlation coefficient employed in the feedforward coefficients will be given by ??(??)=12????-1???=0???????(??)+??????(??)?.(12) Obviously, for the first iteration ??????(0)=??????(0)=0, and consequently ????=0.

We may note that the receiver that employs “blockwise reliabilities” is referred as IB-DFE with “hard decisions,” while the receiver that employs “symbol reliabilities” is referred as IB-DFE with “soft decisions.” The feedforward coefficients used in both types of IB-DFE receivers are still given by (4), but the feedback loop of the IB-DFE with “hard decisions” uses the estimated data block, weighted by a reliability coefficient common to the entire block, while for IB-DFE with “soft decisions,” the feedback loop uses a different reliability coefficient for each symbol (changes from bit to bit).

With a conventional IB-DFE receiver, the log-likelihood values are computed on a symbol-by-symbol basis (i.e., we do not need to perform the channel decoding in the feedback loop). As an alternative, we can define a turbo IB-DFE that employs the channel decoder outputs instead of the uncoded “soft decisions” in the feedback loop. The main difference between conventional IB-DFE and turbo IB-DFE is in the decision device; in the first case, the decision device is a symbol-by-symbol soft decision (for QPSK constellation this corresponds to the hyperbolic tangent, as in (7)); for the turbo IB-DFE, a soft-in, soft-out (SISO) channel decoder is employed in the feedback loop. The SISO block can be implemented as defined in [13] and provides the LLRs of both the “information bits” and the “coded bits.” The input of the SISO block is LLRs of the “coded bits” at the FDE output, given by ??????(??) and ??????(??).

3. Correlation Coefficient Estimation

As shown in the previous section, the correlation factor is a key parameter for the good performance of IB-DFE receivers, since it supplies a blockwise reliability measure of the estimates employed in the feedback loop. However, contrarily to the channel frequency response, it changes from block to block and iteration to iteration. Therefore, it cannot be computed using reference blocks and needs to be obtained from the equalized output.

In this section, we present several methods to estimate the correlation coefficient. The correlation coefficient ??(??) was defined in (5), and for the sake of simplicity, the iteration number ?? will be ignored in the following equations. Therefore, (5) can be rewritten as ?????=??????*??????||????||2?=??????????*??????||????||2?.(13) Let us consider the transmitted symbols {????;??=0,,??-1} corresponding to a QPSK constellation, under a Gray mapping rule. Consequently, the ???? may be written as ????=??????+??????=±??±????,(14) in which ??????=Re{????} is the “in-phase” component of ????, and ??????=Im{????} is the “quadrature” component of ????. Also, ?? is given by ??=??/2, where ?? denotes the minimum Euclidean distance between two constellation symbols. Under these conditions, ??[|????|2]=??2/4. The time-domain samples at the FDE's output are given by ????=??????+??????=????+????,(15) where ???? denotes the global error term, which is Gaussian-distributed, with zeromean. The symbol estimates are then given by ????=????+??????+????????,(16) where ?????? denotes the error coefficient in ??????, and ?????? denotes the error coefficient in ??????. Clearly, if ?????? and ?????? have no errors, then ?????? and ?????? are null. On the other hand, if there are errors, then the value of ?????? and/or ?????? will be ±??. Consequently, ?????? and ?????? are random variables, with values 0, +?? and -?? with probabilities 1-2????, ???? and ???? respectively, where ???? denotes the bit error rate. Therefore, it can be shown that ??=1-2????.(17) Naturally, in practice, we do not know the transmitted symbols {????;??=0,,??-1}.

3.1. Method I: Estimation Based on the BER Estimate

The total variance of the overall noise plus residual ISI, ??2????, is given by ??2????=12???||????-????||2?.(18) The exact value of ??2???? can not be used in practice, because we do not know the transmitted symbols {????;??=0,,??-1}. Therefore, we may use an approximated value of ??2????, given by ???2????=12????-1???=0||????-????||2.(19) Provided that the numbers of decision errors in ???? are small and ?? is high, we have ??2????˜???2????.(20) For a QPSK constellation, the estimated BER, denoted by ?????, can be approximated by ??????1˜??????????,(21) and, from (17), the estimated value of ?? will be?????=1-2??.(22)

3.2. Method II: Estimation Based on the LLR

This technique consists in the estimation of the correlation coefficient with resort to the log-likelihood ratio (LLRs), as derived in Section 2.2. However, as in method I, an approximated value for the variance of channel and interference noise is considered, instead of the optimum value ???2????=12????-1???=0||????-????||2.(23) The calculation of the LLRs of the “in-phase bit” and the “quadrature bit”, associated to ?????? and ??????, are based in (8) and (9), respectively, but with ???2???? in place of ??2????.

Consequently, the reliabilities related to the “in-phase bit” and the “quadrature bit” of the ??th symbol, are based on (11), respectively, but with ??????? and ??????? in place of ?????? and ??????.

Therefore, the estimated value of the correlation coefficient will be given by1??=2????-1???=0???????+???????.(24)

3.3. Method III: Estimation Based on the MSE

This estimation method is based on the same technique employed in method I, given by (22). However, unlike method I, the BER calculation employs the mean-squared error (MSE) of the equalized samples, instead of an approximated value for the variance of the channel and interference noise. Therefore, in method III, the estimated bit error rate denoted by ????? is given by ?????????????=??||??||2||??·????||2??2MSE????.(25) Once again, assuming a QPSK constellation with a gray mapping and ????=±1±??, then (25) can be written as ??????????||??||=?????2MSE?????,(26) where ??2MSE corresponds to the MSE and is defined as??2MSE=1??2??-1???=0T??,(27) withT???||???=????-????||2?.(28) By combining (1) and (6), and assuming ??[|????·?*??|]=0, ??[|????·??*??|]=0, and ??[|????·?*??|]=0, we can rewrite (28) asT???||???=????-????||2??||??=????????????+????????-???????????-????||2??||??=????????-??????2||-12?·2·??2???||??+????||2?·2·??2???||??+????||2???2·?1-??2?·2·??2??,(29) where the variance of the transmitted frequency-domain data symbols is given by??2??=???||????||2?2,(30) the variance of the channel noise is obtained by??2??=???||????||2?2,(31) the expected value of the error term for the ??th frequency-domain “hard decision” estimate is defined as???||???||2?=??2???1-??2?,(32) and the estimated value of ?? can be expressed by?????=1-2??,(33) with ????? given by (26).

4. Performance Results

In the following, we present a set of results obtained during the research of the impact of the correlation coefficient estimation on the performance of IB-DFE receivers, for 10 iterations, for a given channel realization. We consider uncoded transmissions of FFT blocks with ??=1024 data symbols, selected from a QPSK constellation under a gray mapping rule. Similar results were observed for other values of ??, provided that ??»1.

4.1. Results without Correlation Coefficient Compensation

Figure 2 presents the evolution of the correlation coefficient for values of ????/??0. The results show the differences between the ideal value of the correlation coefficient (obtained with (13)), and the estimated value for methods I, II, and III (obtained by (22), (24) and (33), resp.). As we can see, for methods I and II, the ?? curve suffered a deviation from the curve corresponding to the evolution of the optimum (true) correlation coefficient. This is caused by the use of the approximated value ??????? in the correlation coefficient calculation of both methods, which leads to an optimistic ???? derivation for method I, while as for method II results in optimistic versions of the LLRs. Consequently, these effects will originate optimistic estimates of the transmitted symbols in the feedback loop. In the case of method III, the ?? curve is very close to the optimum (true) ?? curve. This means that the reliability measure estimation, with resort to the MSE, used in the feedback loop is very accurate. It is important to note that the curve related with method III is presented for comparison purposes only, since the variance of the channel and interference noise estimation, ???????, is not considered in the estimation process.

Next we present the BER performance results. Since that the results were similar for all methods, as an example, we present in Figure 3 the results related with method I. These results are expressed as function of ????/??0, where ??0 is the one-sided power spectral density of the noise, and ???? is the energy of the transmitted bits (i.e., the degradation due to the useless power spent on the cyclic prefix is not included). These results also refer to a scenario where the IB-DFE suffers from error propagation.

Figure 4 presents the required BER to achieve ????/??0=6 (dB) as a function of the iterations number. Clearly, there is a performance degradation after a few iterations. Although the approximation in (20) may seem reliable, the true value of ??2???? is in fact higher than the estimated one due to decision errors. Consequently, the estimated BER, ?????, is lower than the true bit error rate, ????, leading to ??=??, that is, we are assuming that the estimates used in the feedback loop are more reliable than what they are in fact. It should be pointed out that while underestimating ?? (i.e., using ??<??) leads to a slower convergence of the IB-DFE, overestimating ?? (i.e., using ??>??) leads to a fast convergence but worse BER values.

4.2. Results with Correlation Coefficient Compensation
4.2.1. Technique Description

To overcome the problem of using an optimistic version of transmitted symbols estimates in the feedback loop, we propose a technique to compensate the inaccuracy of the correlation coefficient estimation.

The compensation factor denoted, by ??(??), can be expressed as ??(??)=??/??, where ?? is the estimated correlation coefficient derived from a given estimation method. The curves in Figure 5, obtained by simulation, show the relation between the correlation coefficient estimation, ??, and the compensation factor, ??(??), for methods I, II, and III. Clearly for each method, we can determine the value of the compensation factor ??(??) for each correspondent value of ??, by a simple interpretation of the plot presented in Figure 5. Thus, by knowing ??(??), we can determine a very precise estimated version of the optimum correlation factor given by ??˜??Est+Comp=??,??(??)(34) where ??Est+Comp denotes the compensated correlation coefficient. Therefore, ??Est+Comp can then be used in the derivation of the feedback and feedforward coefficients parameters of the IB-DFE receiver.

4.2.2. Results

We now present the performance results obtained regarding the correlation coefficient compensation. Figure 6 presents the evolution of the correlation coefficient as function of ????/??0. As we can see, for methods I and II the curves are now very close to the curve corresponding to the optimum (true) correlation coefficient evolution, when compared with the results obtained without compensation given by Figure 2. This demonstrates that the correlation coefficient inaccuracy, due to the use of optimistic estimates of the transmitted symbols in the feedback loop, can be efficiently avoided with resort to the proposed compensation technique. As expected for method III, the impact of the correlation coefficient with compensation in the evolution of the correlation coefficient is not so significative. This can be explained by the fact that correlation coefficient ?? obtained under method III was already close to the curve corresponding to the optimum (true) correlation coefficient evolution.

In Figure 7, we present the BER performance results. Once again, since that the results were similar to all methods, as an example, we present in Figure 3 the results related with method I. The improvements obtained with compensation are very significative, since that the performance curves obtained with the compensated correlation coefficient (denoted as “??Est1+Comp”, “??Est2+Comp,” and “??Est3+Comp” in the figures) are now very close to those obtained with optimum (true) correlation coefficient (denoted as “??True” in the figures), as the error propagation problem was corrected.

Figure 8 presents the required BER to achieve ????/??0=6 (dB) as a function of the iterations number. Clearly, there is a performance enhancement when compared with the correlation coefficient estimates without compensation for all methods.

5. Conclusions

In this paper, we studied the impact of the correlation factor estimation in the performance of IB-DFE receivers. Several methods to estimate the correlation coefficient were presented as well as a technique to compensate the estimation errors. It was shown that the inaccuracy due to the use of optimistic estimates of the transmitted symbols in the feedback loop can be efficiently avoided with resort to the proposed compensation technique.

Acknowledgment

This work was partially supported by Fundação para a Ciência e Tecnologia (pluriannual funding and ADCOD projects).