Recently, a closed-form approximated expression was derived by the same author for the achievable residual intersymbol interference (ISI) case that depends on the step-size parameter, equalizer's tap length, input signal statistics, signal to noise ratio (SNR), and channel power and is valid for fractional Gaussian noise (fGn) input where the Hurst exponent is in the region of 0.5≤H<1. But this expression was obtained for the blind adaptive case and cannot be applied to the nonblind adaptive version. Up to now, the achievable residual ISI for the non-blind adaptive case could be obtained only via simulation. In this paper, we derive a closed-form approximated expression (or an upper limit) for the residual ISI obtained by non-blind adaptive equalizers valid for fractional Gaussian noise (fGn) input where the Hurst exponent is in the region of 0.5≤H<1. This new obtained expression depends on the step-size parameter, equalizer's tap length, input signal statistics, SNR, channel power, and the Hurst exponent parameter. Simulation results indicate that there is a high correlation between the calculated results (obtained from the new obtained expression for the residual ISI) and those obtained from simulating the system.
1. Introduction
We consider a nonblind deconvolution problem in which we observe the output of an unknown, possibly nonminimum phase, linear system (single-input-single-output (SISO) FIR system) from which we want to recover its input (source) using an adjustable linear filter (equalizer) and training symbols. During transmission, a source signal undergoes a convolutive distortion between its symbols and the channel impulse response. This distortion is referred to as ISI [1, 2]. It is well known that ISI is a limiting factor in many communication environments where it causes an irreducible degradation of the bit error rate thus imposing an upper limit on the data symbol rate. In order to overcome the ISI problem, an equalizer is implemented in those systems [1–12].
In this paper, we consider the nonblind adaptive equalizer where training sequences are needed to generate the error that is fed into the adaptive mechanism which updates the equalizer’s taps [9–12]. The nonblind adaptive approach yields in most cases a better equalization performance considering convergence speed and equalization quality compared with the blind adaptive version [6]. In addition, the blind adaptive version has a higher computational cost compared with its nonblind approach [6].
The equalization performance from the residual ISI point of view depends on the channel characteristics, on the added noise, on the step-size parameter used in the adaptation process, on the equalizer’s tap length and on the input signal statistics [13, 14]. Fast convergence speed and reaching a residual ISI where the eye diagram is considered to be open (for the communication case) are the main requirements from a blind or nonblind equalizer. Fast convergence speed may be obtained by increasing the step-size parameter. But, increasing the step-size parameter may lead to a higher residual ISI which might not meet the system’s requirements any more. Recently [2], a closed-form approximated expression was derived for the achievable residual ISI case that depends on the step-size parameter, equalizer’s tap length, input signal statistics, SNR, Hurst exponent, and channel power. But this expression is valid only for the blind adaptive case and cannot be used for the nonblind version.
Up to now, the achievable residual ISI for the nonblind adaptive case (for the noisy or noiseless case) could be obtained only via simulation. Thus, the system designer had to spend a lot of time in simulating the whole system in order to find the best values for the step-size parameter and equalizer’s tap length that meet the system’s requirements from the residual ISI point of view. In this paper, we derive a closed-form approximated expression (or an upper limit) for the residual ISI obtained by nonblind adaptive equalizers that depends on the step-size parameter, equalizer’s tap length, input signal statistics, SNR, channel power, and Hurst exponent parameter. This expression is valid for fGn input where the Hurst exponent is in the region of 0.5≤H<1. Please note, H=1 is the limit case, which does not have much practical sense [15–17]. It should be pointed out that a white Gaussian process is a special case (H=0.5) of the fractional Gaussian noise (fGn) model [18]. FGn with H∈(0.5,1) corresponds to the case of long-range dependency (LRD) [18]. Thus, the new obtained expression for the achievable residual ISI is not only valid for the special case of white Gaussian process but also covers those cases that correspond to the case of LRD.
The paper is organized as follows. After having described the system under consideration in Section 2, the closed-form approximated expression (or upper limit) for the residual ISI is introduced in Section 3. In Section 4, simulation results are presented, and the conclusion is given in Section 5.
2. System Description
The system under consideration is illustrated in Figure 1, where we make the following assumptions.
The input sequence x[n] belongs to a two independent quadrature carriers case constellation input with variance σx2, where xr[n] and xi[n] are the real and imaginary parts of x[n], respectively, and σxr2 is the variance of xr[n].
The unknown channel h[n] is a possibly nonminimum phase linear time-invariant filter in which the transfer function has no “deep zeros"; namely, the zeros lie sufficiently far from the unit circle.
The equalizer c[n] is a tap-delay line.
The noise w[n] consists of w[n]=wr[n]+jwi[n], where wr[n] and wi[n] are the real and imaginary parts of w[n], respectively, and wr[n] and wi[n] are independent. Both wr[n] and wi[n] are fractional Gaussian noises (fGn) with zero mean. Note that σwr2=E[wr2[n]], σwi2=E[wi2[n]], for m≠k:E[wr[n-k]wr[n-m]] = (σwr2/2)[(|m-k|-1)2H-2(|m-k|)2H+(|m-k|+1)2H] and E[wi[n-k]wi[n-m]] = (σwi2/2)[(|m-k|-1)2H-2(|m-k|)2H+(|m-k|+1)2H], where E[·] denotes the expectation operator on (·) and H is the Hurst exponent.
The variance of w[n] is denoted as E[w[n]w*[n]]=σw2, where σw2=2σwi2=2σwr2 and (·)* is the conjugate operation on (·).
Block diagram of a baseband communication system.
The transmitted sequence x[n] is sent through the channel h[n] and is corrupted with noise w[n]. Therefore, the equalizer’s input sequence y[n] may be written as
(1)y[n]=x[n]*h[n]+w[n],
where “*” denotes the convolution operation. The equalized output signal can be written as
(2)z[n]=x[n]+p[n]+w~[n],
where p[n] is the convolutional noise, namely, the residual intersymbol interference (ISI) arising from the difference between the ideal equalizer’s coefficients and those chosen in the system and w~[n]=w[n]*c[n]. The ISI is often used as a measure of performance in equalizers’ applications, defined by
(3)ISI=∑m~|s~m~|2-|s~|max2|s~|max2,
where |s~|max is the component of s~, given in (4), having the maximal absolute value. Consider that
(4)s~[n]=c[n]*h[n]=δ[n]+ζ[n],
where δ is the Kronecker delta function and ζ[n] stands for the difference (error) between the ideal and the actual value used for c[n].
Next, we turn to the adaptation mechanism of the equalizer which is based on training symbols [9–12, 19]:
(5)c_eq[n+1]=c_eq[n]-μ(z[n]-x[n])y_*[n],
where μ is the step-size parameter, c_eq[n] is the equalizer vector where the input vector is y_[n]=[y[n]⋯y[n-N+1]]T, and N is the equalizer’s tap length. The operator ()T denotes for transpose of the function (). Please note that for the nonblind adaptive case, during the training period, a known data sequence is transmitted. A replica of this sequence is made available at the receiver in proper synchronism with the transmitter, thereby making it possible for adjustments to be made to the equalizer coefficients in accordance with the adaptive filtering algorithm employed in the equalizer design [19].
3. Residual ISI for Fractional Gaussian Noise Input
In this section, a closed-form approximated expression (or an upper limit) is derived for the residual ISI valid for the fGn input case.
Theorem 1.
Noted the following assumptions.
The convolutional noise p[n] is a zero mean, white Gaussian process with variance σp2=E[p[n]p*[n]]. The real part of p[n] is denoted as pr[n] and E[pr2[n]]=mp.
The source signal x[n] is a rectangular Quadrature Amplitude Modulation (QAM) signal (where the real part of x[n] is independent of the imaginary part of x[n]) with known variance and higher moments.
The convolutional noise p[n] and the source signal are independent.
The gain between the source and equalized output signal is equal to one.
The convolutional noise p[n] is independent of w~[n].
The added noise is fGn with zero mean.
The channel h[n] has real coefficients.
The Hurst exponent is in the range of 0.5≤H<1.
The residual ISI expressed in dB units may be defined as
(6)ISI=10log10(mp)-10log10(σxr2),
where mp is defined by
(7)mp=σw~r2(-μ(Nσx2∑k=0k=R-1hk2[n]+Nσx2SNR)2)-1μ(Nσx2∑k=0k=R-1hk2[n]+Nσx2SNR)2×((Nσx2∑k=0k=R-1hk2[n]+Nσx2SNR)22(Nσx2∑k=0k=R-1hk2[n]+Nσx2SNR)-μ(Nσx2∑k=0k=R-1hk2[n]+Nσx2SNR)2)-1)=σw~r2μ(Nσx2∑k=0k=R-1hk2[n]+(Nσx2/SNR))2-μ(Nσx2∑k=0k=R-1hk2[n]+(Nσx2/SNR)),σw~r2≅σxr2SNR∑k=0k=R-1hk2[n][1+(N-1)H(2H-1)],
and R is the channel length, σw~r2 is the variance of w~r[n] (w~r[n] is the real part of w~[n]), and SNR is given by SNR=σxr2/σwr2=σx2/σw2.
Comments. It should be pointed out that assumptions (1)–(5) from above are precisely the same assumptions made in [2, 14].
Proof.
Let us first recall the expression for the adaptation mechanism of the equalizer given in (5). Then, we substitute (2) into (5) and obtain
(8)c_eq[n+1]=c_eq[n]-μ(p[n]+w~[n])y_*[n].
Next, we recall from [14] the expression for E[Δ(pr2)], where pr is the real part of p[n] and Δ(pr2)=pr2[n+1]-pr2[n]:
(9)E[Δ(pr2)]≅≅-2E[pr(μPr(z)∑m=0m=N-1y[n-m]y*[n-m])]+++E[(-μPr(z)∑m=0m=N-1y[n-m]y*[n-m])2],
where Pr(z) is the real part of P(z) and is given in our case as
(10)P(z)=z[n]-x[n]=p[n]+w~[n]⟹Pr(z)=pr+w~r[n].
According to [13, 14], when the equalizer has converged, we may assume that E[Δ(pr2)]≅0. Therefore, by setting E[Δ(pr2)]=0 into (9), we obtain
(11)-2μmpE[∑m=0m=N-1y[n-m]y*[n-m]]+μ2(mp+σw~r2)×E[(∑m=0m=N-1y[n-m]y*[n-m])2]=0⇓mp=σw~r2(-μE[(∑m=0m=N-1y[n-m]y*[n-m])2])-1μE[(∑m=0m=N-1y[n-m]y*[n-m])2]×××××××(2E[∑m=0m=N-1y[n-m]y*[n-m]]-μE----------×[(∑m=0m=N-1y[n-m]y*[n-m])2])-1).
In [13], the expression E[(∑m=0m=N-1y[n-m]y*[n-m])2] was approximated as (E[∑m=0m=N-1y[n-m]y*[n-m]])2. It should be pointed out that this approximation fits the MPSK case where QPSK is a special case of it. However, satisfying results were obtained in [13] for the 16QAM and 64QAM cases in spite of the fact that the above mentioned approximation was applied. Thus, it makes sense to use the same approximation also here for our case. The expression E[∑m=0m=N-1y[n-m]y*[n-m]] may be written as
(12)E[∑m=0m=N-1y[n-m]y*[n-m]]=Nσx2∑k=0k=R-1hk2[n]+Nσx2SNR.
Next, we turn to find a closed-form approximated expression for σw~r2. The real part of w~[n], namely, w~r[n], may be expressed as
(13)w~r[n]=∑k=0k=N-1ck[n]wr[n-k].
Thus, the variance of w~r[n] is given by:
(14)σw~r2=E[∑k=0k=N-1ck[n]wr[n-k]∑m=0m=N-1cm[n]wr[n-m]]=∑k=0k=N-1∑m=0m=N-1ck[n]cm[n]E[wr[n-k]wr[n-m]],
which can be also written as
(15)σw~r2=σwr2∑k=0k=N-1ck2[n]+∑k=0,k≠mk=N-1∑m=0,k≠mm=N-1ck[n]cm[n]E[wr[n-k]wr[n-m]].
According to [2], expression (15) can be approximately written as
(16)σw~r2≅σxr2SNR∑k=0k=R-1hk2[n][1+N-1H(2H-1)],
by using assumption (4) from the system description section, assumptions (4) and (6)–(8) from this section, the Holder inequality [20], and the following approximation [21]:
(17)0.5[(|m-k|-1)2H-2(|m-k|)2H+(|m-k|+1)2H]≃≃H(2H-1)|m-k|2H-2.
Now, by substituting (12) and (16) into (11) we obtain (7). This completes our proof.
4. Simulation
In this section, we test our new proposed expression for the residual ISI for the 16QAM case (a modulation using ±{1,3} levels for in-phase and quadrature components) with the algorithm described in (5) for different SNR, step-size, and equalizer’s tap length values and for two different channel types. The following two channels were considered. Channel1 (initial ISI=0.44): the channel parameters were determined according to [22]
(18)hn=(0.4n-10forn<0;-0.4forn=0××0.84·0.4n-1forn>0).
Channel2 (initial ISI=0.88) the channel parameters were determined according to
(19)hn=(0.4851,-0.72765,-0.4851).
The equalizer was initialized by setting the center tap equal to one and all others to zero.
In the following, we denote the residual ISI performance according to (6) with (7) as “Calculated ISI.” Figure 2 to Figure 8 show the ISI performance as a function of the iteration number of our proposed expression (6) with (7) for the achievable residual ISI compared with the simulated results for two different channels and equalizer’s tap length and various values for H, SNR and step-size parameter. According to Figures 2, 3, 5, 6, 7, and 8, a high correlation is observed between the simulated and calculated results even for H=0.9. According to Figure 4, the calculated ISI may be considered as an upper limit for the simulated results for H>0.5.
A comparison between the simulated and calculated residual ISI for the 16QAM source input going through channel1 for SNR=10 [dB]. The averaged results were obtained in 100 Monte Carlo trials. The equalizer’s tap length and step-size parameter were set to 13 and 0.0006 respectively.
A comparison between the simulated and calculated residual ISI for the 16QAM source input going through channel1 for SNR=8 [dB]. The averaged results were obtained in 100 Monte Carlo trials. The equalizer’s tap length and step-size parameter were set to 13 and 0.0006 respectively.
A comparison between the simulated and calculated residual ISI for the 16QAM source input going through channel1 for SNR=10 [dB]. The averaged results were obtained in 100 Monte Carlo trials. The equalizer’s tap length and step-size parameter were set to 27 and 0.0006 respectively.
A comparison between the simulated and calculated residual ISI for the 16QAM source input going through channel2 for SNR=12 [dB]. The averaged results were obtained in 100 Monte Carlo trials. The equalizer’s tap length and step-size parameter were set to 13 and 0.0006 respectively.
A comparison between the simulated and calculated residual ISI for the 16QAM source input going through channel2 for SNR=12 [dB]. The averaged results were obtained in 100 Monte Carlo trials. The equalizer’s tap length and step-size parameter were set to 27 and 0.0006 respectively.
A comparison between the simulated and calculated residual ISI for the 16QAM source input going through channel2 for SNR=12 [dB]. The averaged results were obtained in 100 Monte Carlo trials. The equalizer’s tap length and step-size parameter were set to 27 and 0.0004 respectively.
A comparison between the simulated and calculated residual ISI for the 16QAM source input going through channel2 for SNR=12 [dB]. The averaged results were obtained in 100 Monte Carlo trials. The equalizer’s tap length and step-size parameter were set to 27 and 0.0002 respectively.
5. Conclusion
In this paper, we proposed a closed-form approximated expression (or an upper limit) for the residual ISI obtained by nonblind adaptive equalizers valid for the fGn input case where the Hurst exponent is in the region of 0.5≤H<1. This new obtained expression depends on the step-size parameter, equalizer’s tap length, input signal statistics, SNR, channel power, and the Hurst exponent parameter. According to simulation results, a high correlation is obtained between the calculated and simulated results for the residual ISI for some cases, while for others the new obtained expression is a relative tight upper limit for the averaged residual ISI results.
PinchasM.Two blind adaptive equalizers connected in series for equalization performance improvement201341647110.4236/jsip.2013.41008PinchasM.Residual ISI obtained by blind adaptive equalizers and fractional noise2013201311972174MR306270110.1155/2013/972174PinchasM.BobrovskyB. Z.A maximum entropy approach for blind deconvolution20068610291329312-s2.0-3374571209910.1016/j.sigpro.2005.12.009PinchasM.BobrovskyB. Z.A novel HOS approach for blind channel equalization2007638758862-s2.0-3394760615910.1109/TWC.2007.04404PinchasM.2009VDM Verlagsservice Gesellschaft mbHPinchasM.2012Bentham Science Publishers10.2174/97816080535201120101ImG. H.ParkC. J.WonH. C.A blind equalization with the sign algorithm for broadband access20015270722-s2.0-003524815310.1109/4234.905939GodardD. N.Self recovering equalization and carrier tracking in two-dimenional data communication system19802811186718752-s2.0-0019079180ReuterM.ZeidlerJ. R.Nonlinear effects in LMS adaptive equalizers1999476157015792-s2.0-003265473210.1109/78.765126MakkiA. H. I.DeyA. K.KhanM. A.Comparative study on LMS and CMA channel equalizationProceedings of the International Conference on Information Society (i-Society '10)June 20104874892-s2.0-80053365617TucuE.AkirF.ZenA.New step size control technique for blind and non-blind equalization algorithms201322144http://www.academypublisher.com/proc/iscsct10/papers/iscsct10p256.pdfPinchasM.A closed approximated formed expression for the achievable residual intersymbol interference obtained by blind equalizers2010906194019622-s2.0-7674913943810.1016/j.sigpro.2009.12.014PinchasM.A new closed approximated formed expression for the achievable residual ISI obtained by adaptive blind equalizers for the noisy caseProceedings of the IEEE International Conference on Wireless Communications, Networking and Information Security (WCNIS '10)June 2010Beijing, China26302-s2.0-7795767900010.1109/WCINS.2010.5541879BeranJ.199461New York, NY, USAChapman and HallMonographs on Statistics and Applied ProbabilityMR1304490LiM.ZhaoW.On 1/f noise2012201223673648MR3007794ZBL1264.94060LiM.Fractal time series—a tutorial review20102010262-s2.0-7795148927610.1155/2010/157264157264MR2570932ZBL1191.37002LiM.ZhaoW.Quantitatively investigating locally weak stationarity of modified multifractional Gaussian noise2012391246268627810.1016/j.physa.2012.07.043MalikG.SappalA. S.Adaptive equalization algorithms: an overview201123SpiegelM. R.1968New York, NY, USAMcGraw-HillSchaum's Outline SeriesLiM.ZhaoW.On bandlimitedness and lag-limitedness of fractional Gaussian noise201339291955196110.1016/j.physa.2012.12.035MR3038268ShalviO.WeinsteinE.New criteria for blind deconvolution of nonminimum phase systems (channels)19903623123212-s2.0-002540498510.1109/18.52478