A NEW ROBUST FREQUENCY DOMAIN ECHO CANCELLER WITH CLOSED-LOOP
LEARNING RATE ADAPTATION
Jean-Marc Valin and Iain B. Collings

arXiv:1602.08609v1 [cs.SD] 27 Feb 2016

CSIRO ICT Centre, Sydney, Australia (jmvalin@jmvalin.ca)
ABSTRACT
One of the main difficulties in echo cancellation is the fact that
the learning rate needs to vary according to conditions such as
double-talk and echo path change. Several methods have been
proposed to vary the learning. In this paper we propose a new
closed-loop method where the learning rate is proportional to
a misalignment parameter, which is in turn estimated based
on a gradient adaptive approach. The method is presented in
the context of a multidelay block frequency domain (MDF)
echo canceller. We demonstrate that the proposed algorithm
outperforms current popular double-talk detection techniques
by up to 6 dB.
Fig. 1. Block diagram of echo cancellation system.
1. INTRODUCTION
In any echo cancellation system, the presence of near end
speech (double-talk) tends to make the adaptive filter diverge.
To counter the effect, robust echo cancellers require adjustment of the learning rate to account for the presence of doubletalk in the signal.
Most echo cancellation algorithms attempt to explicitly
detect double-talk [1] conditions and then react by freezing
the coefficients of the adaptive filter (setting the learning rate
to zero). Reliable double-talk detection is a difficult problem
and sometimes it is not clear what should be considered as
double-talk, especially in an acoustic echo cancellation context with stationary background noise.
In previous work [2], we proposed a new approach to
make echo cancellation more robust to double-talk. Instead of
attempting to explicitly detect double-talk conditions, a continuous learning rate was used. The learning rate depends on
a misalignment estimate, which is obtained through a linear
regression. While the technique gives good results, the estimation of the misalignment remains a difficult problem.
In this paper, we propose a new approach where the misalignment is estimated in closed-loop based on a gradient adaptive approach. This closed-loop technique is applied to the
block frequency domain (MDF) adaptive filter [3] and shows
a significant improvement over previous approaches.
c 2007 IEEE. Personal use of this material is permitted. Permission
from IEEE must be obtained for all other uses, in any current or future media,
including reprinting/republishing this material for advertising or promotional
purposes, creating new collective works, for resale or redistribution to servers
or lists, or reuse of any copyrighted component of this work in other works.

In Section 2, we discuss the optimal learning rate in presence of noise. Section 3 introduces the multidelay frequency
domain (MDF) adaptive filter. In Section 4, we propose a
gradient adaptive technique for adjusting the learning rate of
the MDF algorithm. Experimental results and a discussion
are presented in Section 5 and Section 6 concludes this paper.
2. OPTIMAL LEARNING RATE IN THE PRESENCE
OF DOUBLE-TALK
In the acoustic echo cancellation context illustrated in Fig. 1,
the speech signal x (n) received from the far end is played
into a loudspeaker. The microphone signal d (n) captures
both the convoluted far end speech y (n) and the near end
noisy speech v (n). The adaptive filter attempts to estimate
the impulse response ĥ (n) to be as close as possible to the
real impulse response h (n). The estimated echo ŷ (n) is subtracted from d (n), so the output signal e (n) contains both
double-talk v (n) and residual echo r (n) = y (n) − ŷ (n).
The conventional approach to double-talk robustness consists of setting the learning rate to zero when double-talk is
detected. Double-talk detectors [1] are thus an important aspect of the approach. Unfortunately, they are sometimes unreliable, especially in acoustic echo cancellation context when
background noise is present. In this paper, we investigate continuous learning rates that do not depend on a binary doubletalk decision.
Whenever an adaptive filter is not perfectly adjusted, its

residual signal r (n) can be used to gain better information
about the exact (time-varying) filter weights h (n). However,
the amount of information about h (n) present in e (n) decreases with the amount of noise and near end speech v(n). In
the case of the normalised least mean square (NLMS) filter, it
means that the stochastic gradient becomes less reliable when
the noise increases or when the filter misalignment decreases
(as the filter converges). The theoretical optimal learning rate
is approximately equal to the residual-to-error ratio [2]:

E r2 (n)
µopt (n) ≈
(1)
E {e2 (n)}
where r (n) = y (n) − ŷ (n) is the (unknown) residual echo
and e (n) is the error signal.
One possible method to vary the learning rate µ (n) would
be to use the generalized normalized gradient descent (GNGD)
algorithm [4], which includes the NLMS learning rate:

µ0 E x2 (n)
(2)
µ (n) =
E {x2 (n)} + ǫ (n)
where ǫ (n) is adapted based on the NLMS stochastic gradient
behaviour. To examine ǫ (n) more closely, it is reasonable to
surmise that (2) eventually converges to the optimal learning
rate defined by (1). Assuming steady state behaviour (ǫ (n) is
stable) and µ0 = 1, we find that:


E r2 (n)
E r2 (n)
≈
(3)
E {r2 (n)} + γ (n) ǫ (n)
E {e2 (n)}
 2
 2
where
 2 γ = E r (n)2 /E x (n)
 . Knowing that we have
E e (n) = E r (n) + E v 2 (n) , we find the relation ǫ (n) ≈ E v 2 (n) /γ (n). In other words, the gradientadaptive parameter ǫ (n) is approximately proportional to the
variance of the near-field signal and independent of the farfield signal. Because ǫ (n) can only be adapted slowly
 over
time, it is clear that (2) implicitly assumes that E v 2 (n)
also varies slowly. While this is a reasonable assumption in
applications where the background noise is constant or slowly
varying, the assumption does not hold for acoustic echo cancellation, where double-talk can start or stop at any time.
In previous work [2], we proposed
to use (1) directly to

adapt the learning rate. While E e2 (n) can
 easily be estimated, the estimation of the residual echo E r2 (n) is difficult because one does not have access to the real filter coefficients. One reasonable assumption we make is that:



E r2 (n) = η (n) E y 2 (n) ≈ η (n) E ŷ 2 (n)
(4)
where η (n) is the normalised filter misalignment (or the inverse of the echo return loss enancement) and is easier to estimate because it is assumed to vary slowly as a function of
time. Although direct estimation of η (n) through linear regression can lead to good results, estimating η (n) remains a
difficult problem. In this paper we propose to apply a gradient
adaptive approach to the problem of estimating η (n).

3. THE MDF ALGORITHM
In this paper, we consider the special case of the multidelay
block frequency domain (MDF) adaptive filter [3]. The MDF
algorithm in matrix form is detailed here for the sake of clarity. Let N be the MDF block size, K be the number of blocks
and F denote the 2N × 2N Fourier transform matrix, we denote the frequency-domain signals for frame ℓ as:
T

(5)
e (ℓ) = F [01×N , e (ℓN ) , . . . , e (ℓN + N − 1)]
o
n
T
xk(ℓ) = diag F[x((ℓ−k−1)N ),. . . ,x((ℓ−k+1)N −1)]
(6)



X (ℓ) = x0 , x1 , . . . , xK−1



(7)

d (ℓ) = F [01×N , d (ℓN ) , . . . , d (ℓN + N − 1)]

T

(8)

The MDF algorithm is then expressed in matrix form as:
e (ℓ) = d (ℓ) − y (ℓ)

(9)

ŷ (ℓ) = G1 X (ℓ) ĥ (ℓ)

(10)

ĥ (ℓ + 1) = ĥ (ℓ) + G2 µ (ℓ) ∇ĥ (ℓ)

(11)

H
∇ĥ (ℓ) = Φ−1
xx (ℓ) X (ℓ) e (ℓ)

(12)

where Φxx (ℓ) is the diagonal normalisation matrix as computed in [3], G1 and G2 are the constraint matrices:


0N ×N 0N ×N
F−1
(13)
G1 = F
0N ×N IN ×N


IN ×N 0N ×N
F−1
(14)
G̃2 = F
0N ×N 0N ×N
o
n
(15)
G2 = diag G̃2 , G̃2 , . . . , G̃2
and µ (ℓ) is the 2KN × 2KN diagonal learning rate matrix:
n
T o
(16)
µ (ℓ) = diag µT (ℓ) , µT (ℓ) , . . . , µT (ℓ)
µ (ℓ) = [µ0 (ℓ) , µ1 (ℓ) , . . . , µ2N −1 (ℓ)]

T

(17)

If µk (ℓ) = µo , we have the standard MDF algorithm.
4. GRADIENT-ADAPTIVE LEARNING RATE
In [2], we proposed to use the frequency-dependent learning
rate:


2
(ℓ)
ŷ
k


(18)
µk (ℓ) = min η (ℓ)
2 , µ0 
|ek (ℓ)|
where η(ℓ) is a direct estimation of the normalised misalignment exploiting the non-stationarity of the signals and using
linear regression between the power spectra of the estimated
echo and the output signal. The motivation behind this formulation is that it factors the residual echo estimation into a

slowly-evolving (but unfortunately difficult to estimate) normalised misalignment η(ℓ) and a rapidly-evolving (but easy to
2

estimate) term far-end term yk (ℓ) . The learning rate can
thus react quickly to double-talk even if the estimation of the
residual echo (through the misalignment estimate) requires a
longer time period. A remaining problem with that approach
is that η(ℓ) is difficult to estimate and the algorithm does not
know whether its estimate of η(ℓ) is too low or too high. In
this sense, the η(ℓ) update in [2] is an open-loop estimate.
4.1. Adaptation algorithm
In this paper we bypass the difficulty of estimating η(ℓ) directly and instead propose a closed-loop gradient adaptive estimation of η(ℓ). The parameter η(ℓ) is no longer an estimate of the normalised misalignment, but is instead adapted
in closed-loop in such a way as to maximise convergence of
the adaptive filter. As with other gradient-adaptive methods
[5, 4] we compute the derivative of the mean square error
E (ℓ) = eH (ℓ) e (ℓ), this time with respect to η (ℓ − 1), using the chain derivation rule:
∂E (ℓ)
∂E (ℓ) ∂e (ℓ) ∂ŷ (ℓ) ∂ ĥ (ℓ) ∂µ (ℓ − 1)
=
∂η(ℓ−1)
∂e (ℓ) ∂ŷ (ℓ) ∂ ĥ (ℓ) ∂µ (ℓ − 1) ∂η (ℓ − 1)
(19)
= −e

≈−

H

H
(ℓ)G1 X(ℓ)G2 Φ−1
xx (ℓ)X (ℓ) e(ℓ)

ŷk (ℓ)

ŷk (ℓ)

2
2

|ek (ℓ)|
(20)

2

|ek (ℓ)|2

H

∇ĥ (ℓ) Φxx ∇ĥ (ℓ)

(21)

We propose to use a filtered version of the gradient with a
multiplicative update, similar to the general approch in [6],
and to drop the normalisation term Φxx . This results in the
new gradient-based learning rate adaptation rule:

n
o 
2
H
(ℓ)
ℜ
ψ
(ℓ)
∇
ĥ
(ℓ)
ŷ
k

n
o 
η (ℓ + 1) = η (ℓ) expρ

2P
∗
|ek (ℓ)|
k ℜ ψk (ℓ) ∇ĥk (ℓ)
(22)

ψ (ℓ + 1) = αψ (ℓ) + ∇ĥ (ℓ)

(23)

where ρ is the step size and α controls the gradient smoothing
(typically ρ = 1, α = 0.9).
An intuitive interpretation for the update equation in (22)
is that when the learning rate is too high (because η (ℓ)is too
high), the gradient oscillates, causing η (ℓ) to decrease. On
the other hand, when the learning rate is too low (perhaps because the echo path has changed) the gradient keeps pointing
in the same direction and η (ℓ) increases. Also, because of
|ŷ (ℓ)|2
the k 2 factor, we can apply (22) at each frame regardless
|ek (ℓ)|

of double-talk or even when no far-end speech is present (in
2

which case ŷk (ℓ) = 0 anyway).
The last aspect that needs to be addressed is the initial
condition. When the filter is initialised, all the weights are set
to zero, which means that ŷ (ℓ) = 0 and no adaptation can
take place in (18) and (22). In order to bootstrap the adaptation process, the learning rate µk (ℓ) is set to a fixed constant
(we use µk (ℓ) = 0.25) for a short time equal to twice the filter
length (only non-zero portions of signal x(n) are taken into
account). This procedure is only necessary when the filter is
initialised and is not required in case of echo path change.
4.2. Algorithm behaviour
Here we show that the adaptive learning rate described above
is able to deal with both double-talk and echo path change
without explicit modelling. From (18) we can see that when
2
double-talk occurs, the denominator |ek (ℓ)| rapidly increases,
causing an instantaneous decrease in the learning rate that
lasts only as long as the double-talk period lasts. In the case
of background noise, the learning rate depends on both the
presence of an echo signal as well as the misalignment estimate. As the filter misalignment becomes smaller, the learning
rate also becomes smaller.
One major difficulty involved in double-talk detection is
the need to distinguish between double-talk and echo path
change since both cause a sudden increase in the filter error
signal. This distinction is made implicitly in the gradientbased adaptation of η(ℓ). When the echo path changes, ∇ĥ (ℓ)
starts pointing steadily in the same direction, thus significantly increasing η(ℓ), which is a clear sign that the filter is no
longer properly adapted.
In gradient adaptive methods [5, 6], the implicit assumption is that both the near-end and the far-end signals are nearly
stationary. We have shown that the GNGD algorithm [4] only
requires the near-end signal to be nearly stationary. In our new
proposed method, both signals can be non-stationary, which
is a requirement for double-talk robustness.
5. RESULTS AND DISCUSSION
The proposed system is evaluated in an acoustic echo cancellation context with background noise, double-talk and a
change in the echo path (16 seconds into a 32-second recording). The two impulse responses used are 1024-sample long
and measured from real recordings in a small office room with
both the microphone and the loudspeaker resting on a desk.
The proposed algorithm is compared to our previous work
[2], to the normalised cross-correlation (NCC) method [1] and
to a baseline with no double-talk detection (no DTD). The
optimal threshold found for the NCC algorithm was 0.35 with
a learning rate µ = 0.25. It was found that choosing µmax =
0.75 as the upper bound on the learning rate gave good results
for the proposed algorithm. In practise, µmax has little impact

25

4.5

proposed
previous [2]
NCC [1]
no DTD

20

4

less double−talk

more double−talk

PESQ LQO-MOS

ERLE (dB)

3.5
15

10

3
2.5
2
proposed
previous [2]
NCC [1]
no DTD
no cancellation

5
1.5
0
-20

less double−talk

-15

-10

more double−talk

-5
0
5
10
Near-end-to-far-end ratio (dB)

15

20

1
-20

-15

-10

-5
0
5
10
Near-end-to-far-end ratio (dB)

15

20

Fig. 2. Steady-state ERLE (first two seconds of adaptation
not considered) as a function of the near-end to far-end ratio.

Fig. 3. PESQ objective listening quality measure (LQOMOS) as a function of the near-end to far-end ratio.

on the algorithm because the gradient-based adaptation would
compensate for a higher value of µmax by reducing η (ℓ). The
GNGD algorithm [4] is not included in the evaluation because
it is not intended for applications where the near-field signal is
highly non-stationary and it was not possible to obtain better
results than the baseline.
Fig. 2 shows the average steady-state (the first 2 seconds
of adaptation are not considered) ERLE for the test data with
different ratios of near-end signal and echo. Clearly, the proposed algorithm performs better than both our previous work
(2 dB average improvement) and the NCC algorithm (6 dB
average improvement). The perceptual quality of the output
speech signal is also evaluated by comparing it to the near
field signal v(n) using the Perceptual Evaluation of Speech
Quality (PESQ) ITU-T recommendations P.862 [7]. The perceptual quality of the speech shown in Fig. 3 is evaluated
based on the entire file, including the adaptation time. It is
again clear that the proposed algorithm performs better than
all other algorithms. It is worth noting that the reason why the
results in Fig. 2 improve with double-talk (unlike in Fig. 2) is
that the signal of interest is the double-talk v(n), so the higher
the double-talk the less (relative) echo in the input signal.

algorithm, we believe future work could apply it to other adaptive filtering algorithms, including the NLMS algorithm.
7. REFERENCES
[1] J. Benesty, D.R. Morgan, and J.H. Cho, “A new class of
doubletalk detectors based on cross-correlation,” IEEE
Trans. on Speech and Audio Processing, vol. 8, no. 2, pp.
168–172, 2000.
[2] J.-M. Valin, “On adjusting the learning rate in frequency
domain echo cancellation with double-talk,” IEEE Trans.
on Audio, Speech and Language Processing, 2007.
[3] J.-S. Soo and K.K. Pang, “Multidelay block frequency
domain adaptive filter,” IEEE Trans. ASSP, vol. 38, no.
2, pp. 373–376, 1990.
[4] D.P. Mandic, “A generalized normalized gradient descent
algorithm,” IEEE Signal Processing Letters, vol. 11, no.
2, pp. 115–118, 2004.
[5] V.J. Mathews and Z. Xie, “A stochastic gradient adaptive
filter with gradient adaptive step size,” IEEE Trans. on
Signal Processing, vol. 41, pp. 2075–2087, 1993.

6. CONCLUSION
We have demonstrated a novel method for adjusting the learning rate of frequency-domain adaptive filters based on a gradient adaptive estimation of the current misalignment. The proposed method performs better than a double-talk detector and
previous work using direct estimation of the misalignment. In
addition, the closed-loop gradient adaptive estimation of η (ℓ)
makes the algorithm conceptually simple and means that there
are very few important parameters to be tuned. Although the
proposed algorithm is presented in the context of the MDF

[6] W.-P. Ang and B. Farhang-Boroujeny, “A new class
of gradient adaptive step-size LMS algorithms,” IEEE
Trans. on Signal Processing, vol. 49, no. 4, pp. 805–810,
2001.
[7] ITU-T, P.862: Perceptual evaluation of speech quality
(PESQ): An objective method for end-to-end speech quality assessment of narrow-band telephone networks and
speech codecs, 2001.

