^{1}

^{1}

^{1}

To solve the problem of large steady state residual error of momentum constant modulus algorithm (CMA) blind equalization, a momentum CMA blind equalization controlled by energy steady state was proposed. The energy of the equalizer weights is estimated during the updating process. According to the adaptive filtering theory, the energy of the equalizer weights reaches to the steady state after the algorithm is converged, and then the momentum can be set to 0 when the energy change rate is less than the threshold, which can avoid the additional gradient noise caused by momentum and further improve the convergence precision of the algorithm. The proposed algorithm takes advantage of momentum to quicken the convergence rate and to avoid the local minimum in the cost function to some extent; meanwhile, it has the same convergence precision with CMA. Computer simulation results show that, compared with CMA, momentum CMA (MCMA) and adaptive momentum CMA (AMCMA) blind equalization, the proposed algorithm has the fastest convergence rate and the same steady state residual error with CMA.

Blind equalization can compensate and track the communication channel characteristic without training sequence, which can save the communication bandwidth and improve the communication quality [

Among the adaptive filtering algorithms based on stochastic gradient descent, the momentum gradient algorithm has been proved to be an effective algorithm for accelerating the convergence rate. Meanwhile, under the condition of non-convex cost function, the gradient momentum algorithm can avoid the shallow local minimum of the cost function to a certain extent, which can improve the global convergence performance [

The essence of blind equalization is to recovery the send signal x ( n ) without the information of send signal and the channel h ( n ) . According to the

communication signal transmission principle, the received signal y ( n ) can be given by

y ( n ) = h ( n ) x ( n ) + n ( n ) (1)

And the output x ˜ ( n ) can be given by

x ˜ ( n ) = w H ( n ) y ( n ) (2)

CMA blind equalization utilize the higher order statistics of the transmission signal indirectly, which cost function can be given by [

J ( n ) = 1 2 [ R − | x ˜ ( n ) | 2 ] 2 (3)

where r is the constant modulus, which can be computed according to transmission signal as follow

R = E [ | x ˜ ( n ) | 4 ] E [ | x ˜ ( n ) | 2 ] (4)

According to stochastic gradient descent algorithm, the blind equalizer weights updating formula of CMA blind equalization can be given by [

w ( n + 1 ) = w ( n ) + μ e ( n ) y * ( n ) x ˜ ( n ) (5)

where μ is the study step size and e ( n ) is the transient gradient error which can be given by

e ( n ) = R − | x ˜ ( n ) | 2 (6)

For the convergence rate of CMA blind equalization is very slow, it needs a large number of information symbols to capture the channel characteristics, which results in high bit error rate (BER) in receiver signal. Especially in high speed communication system or in time varying channel conditions, the performance of CMA blind equalization shows very poor for tracking the channel characteristics. Momentum gradient algorithm can improve the convergence rate significantly, and it can avoid falling into the shallow local minimum of the non-convex cost function to improve the global convergence ability. The momentum CMA blind equalization updates the equalizer weights by Equation (7), which a momentum term adds to the updating formula of CMA.

w ( n + 1 ) = w ( n ) + μ e ( n ) y * ( n ) x ˜ ( n ) + α [ w ( n ) − w ( n − 1 ) ] (7)

where 0 < α < 1 is the momentum factor which controls the contribution of momentum term in the algorithm and α [ w ( n ) − w ( n − 1 ) ] is the momentum term. If let Δ w ( n ) = w ( n + 1 ) − w ( n ) , then

Δ w ( n ) = μ e ( n ) y * ( n ) x ˜ ( n ) + α Δ w ( n − 1 ) (8)

Furthermore, Δ w ( n ) can be rewritten as follow

Δ w ( n ) = ∑ i = 1 n − 1 α n − i ∇ J ( i ) + α n ( w ( 1 ) − w ( 0 ) ) (9)

where ∇ J ( n ) is the transient gradient that is defined according to CMA cost function. To ensure the algorithm has robust convergence performance, it requires | α | < 1 . If α < 0 , the momentum term has no effect on the accelerating convergence of the algorithm. Therefore, the range of momentum factor is 0 < α < 1 . When the ideal equalization condition is reached, it expects Δ w ( n ) = 0 . Meanwhile, Equation (7) shows that when ideal equalization is reached, the momentum term incurs additional gradient noise, which results in big steady state residual error after the algorithm convergence.

Adaptive momentum algorithms update the weights of equalizer by the bigger momentum factor in the initial stage to accelerate the convergence rate, and the momentum factor is gradually reduced to obtain a smaller steady-state residual error in the iterative process, which can make a tradeoff between convergence rate and convergence precision. Typical adaptive momentum algorithms have the following forms

α ( n ) = β e ( n ) / ∑ j = 1 n e ( j ) (10)

where β is ratio coefficient as a variable parameter.

α ( n ) = { 2 3 E ( n ) ≥ 2 E ( n ) / 3 E ( n ) < 2 (11)

where E ( n ) is mean square error which can be estimated by

E ( n + 1 ) = λ E ( n ) + ( 1 − λ ) | x ˜ ( n ) − x ^ ( n ) | 2 (12)

where λ is the forgetting factor and x ^ ( n ) is the decision output of x ˜ ( n ) .

α ( n ) = γ [ 1 − e p − E ( n ) ] (13)

α ( n ) = γ [ 1 − e p − | e ( n ) | 2 ] (14)

In Equation (13) and Equation (14), γ is an adjustable parameter and p is a parameter which control the positive or negative of α ( n ) . Among the above adaptive momentum gradient algorithms, Equation (13) and Equation (14) have little difference. In the convergence direction of the algorithm, both E ( n ) and e ( n ) gradually reduce according to the statistical significance and tend to 0. If the parameters are set appropriate, the momentum factor α ( n ) must have a decreasing trend and if the algorithm obtains convergence result, the following result can be obtained

lim n → ∞ α ( n ) = γ ( 1 − e p ) (15)

Therefore, the adaptive momentum gradient algorithms according to Equation (13) and Equation (14) have the similar form with Equation (10). Furthermore, they have experience parameter p needing to set and more parameters lead the algorithm to have not to be used universally. By contrast, the adaptive momentum gradient algorithm given by Equation (11) needs no parameter to be set, however, the momentum factor has no theory basis in the algorithm and the performance shows difference under different modulate signal or different SNR conditions.

Based on the analysis of momentum gradient algorithm, the momentum factor should has a bigger value to accelerate the convergence rate in the initial stage of the algorithm, when the ideal equalization is reached, the momentum factor should be 0 to avoid the additional gradient noise. A point of view was proposed in [

Hereby an adaptive momentum CMA blind equalization controlled by energy steady state was proposed. According to adaptive filter theory, the weights of equalizer reach to steady state when the algorithm achieves convergence. The energy of the weights of the equalizer is defined as follow

σ ( n ) = ∑ j = 1 L | w j ( n ) | 2 (16)

where L is the length of the blind equalizer. Then the energy change rate of the weights of the blind equalizer can be defined as follow

Δ σ ( n ) = σ ( n ) − σ ( n − 1 ) = ∑ j = 1 L ( | w j ( n ) | 2 − | w j ( n − 1 ) | 2 ) (17)

The energy of the weights of the blind equalizer will reach to steady state, and then the energy change rate tends to 0. However, in the actual communication system, the received signal is often interfered by the noise, which results in the energyof the blind equalizer shows jitter downward trend. Here Δ σ ( n ) is estimated by sliding window as follow

Δ σ ( n + 1 ) = λ Δ σ ( n ) + ( 1 − λ ) Δ σ ( n − 1 ) (18)

where λ is the forgetting factor. When Δ σ ( n ) < 10 − 3 , it can judge the algorithm has converge to the steady state. Therefore, the updating method of the momentum factor of the adaptive momentum CMA blind equalization controlled by the energy change ratio can be given by

α ( n ) = { α ( 0 ) 0 Δ σ ( n ) ≥ 10 − 3 Δ σ ( n ) < 10 − 3 (19)

To verify the performance of the adaptive momentum CMA blind equalization controlled by energy steady state, computer simulation has done to compare with the statistical momentum CMA (SMCMA) and mean square error controlled momentum CMA (MSEMCMA). Because the adaptive momentum algorithm proposed in [

H ( z ) = 1 + 0.4 z − 12 (20)

The transmission signal is equal probability binary sequence and modulate by QPSK. The noise in the channel is adding Gaussian white noise and S N R = 20 d B . The length of the blind equalizer is 48 with the centre tap initializes to 1 and the others taps initialize to 0. The study step size μ = 0.001 . In SMCMA, the parameter of the momentum control function β = 0.75 and the initial momentum factor α ( 0 ) = 0.75 . In the estimation of the mean square error E ( n ) and the energy change rate Δ σ ( n ) , the forgetting factor λ = 0.99 . The initial energy change rate value of the blind equalizer is set Δ σ ( 0 ) = 1 . The performance comparison conducted by the residual inter-symbol interference (ISI) [

I S I ( n ) = ∑ i | C i ( n ) | 2 − | C i ( n ) | max 2 | C i ( n ) | max 2 (21)

where C ( n ) is the combined impulse response of the channel and the blind equalizer. The residual inter-symbol interference convergence result in the simulation achieved across 500 times Monte-Carlo is shown in

^{−3}, which shows 10^{−3} can be taken as the criterion of the momentum factor adjustment.

The simulation results show that ECMCMA can take full advantage of the momentum term to accelerate the convergence rate. Meanwhile, it can avoid the additional gradient noise introduced by momentum term after the algorithm convergence. Also, ECMCMA can be applied with the other improved MCMA to further improve the convergence performance.

Momentum CMA blind equalization can improve the convergence rate and avoid the local minimum of the non-convex cost function during the iterative process. However, the additional gradient noise is introduced by momentum term after the algorithm convergence, which results in the steady state residual error increase. Adaptive momentum CMA blind equalization reduces the momentum factor during the iterative process to obtain a compromise result between convergence rate and convergence precision. This work proposed an adaptive momentum CMA blind equalization controlled by energy steady state, in which the energy change rate of the blind equalizer is defined. The threshold is set according to the energy change rate which controls the scope of action of the momentum term in the algorithm. ECMCMA takes full advantage of the momentum term to accelerate the convergence rate of the algorithm and avoid the additional gradient noise introduced by momentum term. The computer simulation results show the effectiveness of ECMCMA.

This work was partially supported by National Science Foundation of P.R. China (Grant: 61201418), Fundamental Research Funds for the Central Universities (Grant: DC201502060302) and Liaoning Province High School Talent Support Program (Grant: LJQ2013126).

Xiao, Y., Dong, Y.H. and Sun, J.Y. (2017) An Adaptive Momentum CMA Blind Equalization Based on Error Energy. Int. J. Communications, Network and System Sciences, 10, 333-340. https://doi.org/10.4236/ijcns.2017.105B033