DE19525381B4 - Method for controlling the step size of an adaptive-filter echo canceller - Google Patents
Method for controlling the step size of an adaptive-filter echo canceller Download PDFInfo
- Publication number
- DE19525381B4 DE19525381B4 DE1995125381 DE19525381A DE19525381B4 DE 19525381 B4 DE19525381 B4 DE 19525381B4 DE 1995125381 DE1995125381 DE 1995125381 DE 19525381 A DE19525381 A DE 19525381A DE 19525381 B4 DE19525381 B4 DE 19525381B4
- Authority
- DE
- Germany
- Prior art keywords
- input signal
- power
- step size
- vector
- cross
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04M—TELEPHONIC COMMUNICATION
- H04M9/00—Arrangements for interconnection not involving centralised switching
- H04M9/08—Two-way loud-speaking telephone systems with means for conditioning the signal, e.g. for suppressing echoes for one or both directions of traffic
- H04M9/082—Two-way loud-speaking telephone systems with means for conditioning the signal, e.g. for suppressing echoes for one or both directions of traffic using echo cancellers
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B3/00—Line transmission systems
- H04B3/02—Details
- H04B3/20—Reducing echo effects or singing; Opening or closing transmitting path; Conditioning for transmission in one direction or the other
- H04B3/23—Reducing echo effects or singing; Opening or closing transmitting path; Conditioning for transmission in one direction or the other using a replica of transmitted signal in the time domain, e.g. echo cancellers
Landscapes
- Engineering & Computer Science (AREA)
- Signal Processing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Cable Transmission Systems, Equalization Of Radio And Reduction Of Echo (AREA)
Abstract
Verfahren zur Steuerung der Schrittweite eines Echokompensators mit adaptivem Filter, dessen Koeffizienten durch den NLMS-Algorithmus automatisch angepaßt werden und bei dem die optimale Schrittweite aus dem Quotienten der Schätzwerte für die Streuung des Eingangssignals und des gestörten Fehlers aus der Beziehung abgeleitet wird, dadurch gekennzeichnet, daß der Faktor g ^(k) aus der quadrierten Norm des Vektors der Kreuzkorrelationskoeffizienten und der quadrierten Leistung des Eingangssignals über die Beziehung ermittelt wird, wobei in jedem Abtastschritt sowohl die Kreuzkorrelationskoeffizienten zwischen Eingangs- und Fehlersignal als auch ein Schätzwert für die Leistung des Eingangssignals ermittelt werden, wobei der Schätzwert für die Leistung des Eingangssignals aus der Norm des Teils des Eingangssignalvektors abgeleitet wird, der auch der Ermittlung der Kreuzkorrelationsfunktionen zugrunde liegt.Method for controlling the step size of an adaptive filter echo canceller whose coefficients are automatically adjusted by the NLMS algorithm and in which the optimum step size is the quotient of the estimated values for the scattering of the input signal and the perturbed error from the relationship is derived, characterized in that the factor g ^ (k) from the squared norm of the vector of the cross-correlation coefficients and the squared power of the input signal via the relationship In each sampling step, both the cross correlation coefficients between input and error signal and an estimate of the power of the input signal are determined, the estimated value for the power of the input signal being derived from the norm of the part of the input signal vector, which also determines the Based on cross-correlation functions.
Description
Das erfindungsgemäße Verfahren betrifft die Steuerung der Schrittweite bei einem Echokompensator mit einem FIR-Filter, dessen Koeffizienten durch den NLMS-Algorithmus automatisch angepaßt werden.The inventive method relates to the control of the step size in an echo canceller with a FIR filter, whose coefficients are automatically adjusted by the NLMS algorithm.
Ein Echokompensator der genannten Art ist in dem Zeitschriftenaufsatz S.Yamamoto, S.Kitayama, Trans. of the IECE of Japan, Vol. E 65, No. 1, Jan. 1982, S.1–8 beschrieben. Für die optimale Größe des Schrittweitenfaktors wird der Ausdruck angegeben. σx ist die Streuung des Eingangssignals. Die nicht direkt meßbare Streuung der Störung, σn, wird mittels eines zweiten adaptiv eingestellten Filters geschätzt. Die ebenfalls nicht direkt meßbare Parameterfehlernorm ∥Δh∥ wird aus sogenannten Totzeitkoeffizienten (eine Reihe von Koeffizienten, die bei erfolgtem Abgleich 0 werden müssen) geschätzt.An echo canceller of the type mentioned is described in the journal article S.Yamamoto, S.Kitayama, Trans. Of the IECE of Japan, Vol. 1, Jan. 1982, p.1-8. The expression becomes the optimal size of the step size factor specified. σ x is the dispersion of the input signal. The non-directly measurable scattering of the perturbation, σ n , is estimated by means of a second adaptively adjusted filter. The parameter error standard ∥Δh∥ which is also not directly measurable is estimated from so-called dead time coefficients (a series of coefficients which must become 0 when the adjustment has taken place).
Der Mangel dieser Lösung besteht darin, daß bei Änderungen der Parameter des Echopfades die Parameterfehlernorm falsch geschätzt wird. In der Regel ergeben sich zu kleine Schätzwerte, so daß zu kleine Werte für die Schrittweite erhalten werden. Nachteilig ist ferner, daß ein zweites adaptives Filter benötigt wird.Of the Lack of this solution is that when changes the parameter of the echo path the parameter error standard is misjudged. As a rule, too small estimates result, so too small values for the Increment can be obtained. Another disadvantage is that a second adaptive filter needed becomes.
In
dem Buch R. Frenzel, Fortschr.-Ber. VDI-Reihe 10 Nr. 228, Düsseldorf
1992 ist ein Echokompensator beschrieben, bei dem zur Bestimmung
der optimalen Schrittweite von dem Zusammenhang ausgegangen
wird, den man durch Umformung aus Gl. (1) erhalten kann, wenn man
die Streuung des gestörten
Fehlers einführt:
Auch hier wird die Parameterfehlernorm aus den Totzeitkoeffizienten geschätzt. Zur Verbesserung dieser Schätzung bei Parameteränderungen wird vorgeschlagen, aus den Cepstren der Ausgangssignale von Modell und System eine Größe zu gewinnen, die einerseits der Erkennung von Systemänderungen dient und andererseits eine Größe liefert, mit der die Totzeitkoeffizienten in diesem Fall neu initialisiert werden, um so die Schrittweite zu vergrößern.Also Here, the parameter error standard is estimated from the deadtime coefficients. to Improvement of this estimate with parameter changes is proposed from the cepstres of the output signals of model and system to gain a size which on the one hand serves to detect system changes and on the other hand provides a size with which the deadtime coefficients are reinitialized in this case to increase the step size.
Die praktische Erprobung einer derartigen Schrittweitensteuerung hat gezeigt, daß es bei Sprachsignalen sowohl mit als auch ohne Gegensprechen häufig zu Fehlentscheidungen kommt. Im Ergebnis dessen wird die Schrittweite zu groß eingestellt, so daß die Koeffizienten des adaptiven Filters eine große Streuung besitzen. Außerdem ist der numerische Aufwand für die zum Teil komplexen Berechnungen erheblich.The has practical testing of such step size control shown that it in speech signals both with and without intercommunication too often Wrong decisions comes. As a result, the step size becomes set too high, So that the Coefficients of the adaptive filter have a large dispersion. Besides that is the numerical effort for the sometimes complex calculations considerably.
In
dem Aufsatz T.J.Shan, T. Kailath, Adaptive Algortihms with an Automatic
Gain Control Feature, IEEE Transactions on Circuits and Systems,
Vol. 35, No. 1, Jan. 1988 ist ein Echokompensator mit Schrittweitensteuerung
beschrieben, bei dem die Schrittweite aus
Diese besitzen in der Regel unterschiedliche Vorzeichen. Außerdem basiert der Ausdruck für die Schrittweite nicht auf dem grundlegenden Zusammenhang Gl. (1) und ist damit nicht optimal.These usually have different signs. Also based the expression for the step size is not based on the fundamental Eq. (1) and is therefore not optimal.
In
der Europäischen
Patentschrift
In
der Europäischen
Patentschrift
Diese Schrittweitensteuerung basiert ebenfalls nicht auf dem grundlegenden Zusammenhang von Gleichung (1) und funktioniert damit auch nicht optimal.These Step size control is also not based on the basic Relationship of equation (1) and does not work with it optimal.
Das Ziel der Erfindung ist es, eine Schrittweitensteuerung zu schaffen, die die erwähnten Nachteile nicht aufweist und eine möglichst gute Annäherung der Schrittweite an ihre optimale Größe gemäß Gleichung (1) bzw. Gleichung (2) realisiert. Daraus leitet sich die Aufgabe ab, auf möglichst direktem Wege eine Schätzung für die Parameterfehlernorm zu finden.The The aim of the invention is to provide a step size control, the ones mentioned Disadvantages not and the best possible approximation of Step size to its optimal size according to equation (1) or equation (2) realized. From this the task derives on as possible a direct estimate for the Find parameter error standard.
Nachfolgend wird das erfindungsgemäße Verfahren näher erläutert.following becomes the method according to the invention explained in more detail.
Ausgegangen wird von Gleichung (2). Das Quadrat der Parameterfehlernorm wird durch die Summe aus den quadrierten Differenzen der Koeffizienten von adaptivem Filter und Echopfad gebildet: The starting point is equation (2). The square of the parameter error norm is formed by the sum of the squared differences of the coefficients of adaptive filter and echo path:
Das Quadrat der Parameterfehlernorm kann man bei weißem Rauschen als Eingangssignal auch als Leistung-Übertragungsfaktor g interperetieren.The Square of the parameter error standard can be used as an input signal with white noise also as a power transfer factor interperate.
Gleichung (7) ist direkt nicht auswertbar, da die Koeffizienten des Echopfades ja erst mittels des addaptiven FIR-Filters bestimmt werden sollen. Im Ausgangssignal des Echokompensators sind immer das Störsignal und ein vom Eingangssignal x stammender Teil überlagert. Gerade dieser Teil muß bestimmt werden, wobei das Störsignal unterdrückt werden soll.equation (7) is directly not evaluable, since the coefficients of the echo path yes only by means of the addaptive FIR filter to be determined. The output signal of the echo canceller always contains the interfering signal and a part originating from the input signal x superimposed. Just this part must be determined be, with the interfering signal repressed shall be.
Diese Aufgabe kann durch eine Korrelationsanalyse gelöst werden.These Task can be solved by a correlation analysis.
Bei
einem System mit den Koeffizienten Δh(i) mit weißem Rauschen als Eingangssignal
x gilt für
die Kreuzkorrelationskoeffizienten zwischen Ausgang und Eingang
Dabei wurde angenommen, daß die Störung nicht mit dem Eingangssignal korreliert ist. Für das gesuchte Quadrat der Parameterfehlernorm ergibt sich daher It was assumed that the disturbance is not correlated with the input signal. For the sought square of the parameter error standard is therefore
Setzt man dieses Ergebnis in Gleichung (2) ein, dann erhält man für die Schrittweite Substituting this result in equation (2), one obtains for the step size
In der Praxis muß man die Korrelationsfunktionen durch Zeitmittelwerte als Schätzwerte ersetzen. Je größer die Mittelungszeit Tm gewählt wird, desto besser ist die Störunterdrückung: In practice one must replace the correlation functions by time averages as estimates. The greater the averaging time T m is chosen, the better the interference suppression:
Die Zahl der berechneten Korrelationsfunktionen kann kleiner sein, als die Zahl der Koeffizienten des adaptiven Filter. Falls man die Eingangsleistung durch die quadrierte Vektornorm des Eingangsvektors annähert, wird man in diesem Fall σ 2 / x in Gleichung (9) zweckmäßig aus der Vektornom des verkürzten Eingangsvektors berechnen.The Number of calculated correlation functions may be less than the number of coefficients of the adaptive filter. If you have the input power is approximated by the squared vector norm of the input vector one in this case σ 2 / x in Equation (9) expedient from the vector cathedral of the shortened Calculate input vector.
Die
geschilderten Operationen werden anhand von
Der Eingangsvektor befindet sich in einem Schieberegister SR, welches schon Teil des NLMS-Algorithmus ist. In den Korrelatoren Kn, von denen nur drei gezeigt sind, werden die Kreuzkorrelationsfunktionen gemäß Gleichung (11) gebildet. Sie werden anschließend in den Quadrierstufen Qn quadriert und in dem Summator S addiert.The input vector is located in a shift register SR, which is already part of the NLMS algorithm. In the correlators K n , of which only three are shown, the cross-correlation functions are formed according to equation (11). They are then squared in the squaring stages Q n and added in the summer S.
Der Leistungsschätzer LSx1 bildet einen Schätzwert für die Leistung des Eingangssignals aus den Elementen des eventuell verkürzten Eingangsvektors. Sein Quadrat wird anschließend in der Quadrierstufe Qx1 gebildet. Der Dividierer D1 bildet den Quotienten der Ausgangsgrößen des Summators und der Quadrierstufe Qx1. Damit ist der zweite Teil von Gleichung (10) entsprechend dem Erfindungsgedanken ermittelt.The power estimator LS x1 provides an estimate of the power of the input signal from the elements of the possibly shortened input vector. Its square is then formed in the squaring stage Q x1 . The divider D 1 forms the quotient of the outputs of the summer and the squaring stage Q x1 . Thus, the second part of equation (10) is determined according to the concept of the invention.
Zur vollständigen Berechnung von αopt geben in an sich bekannter Weise zwei weitere Leistungsschätzer LSx und LSe an ihren Ausgängen Schätzwerte für die Leistungen des Eingangs- bzw. des Fehlersignals ab. Ihr Quotient wird durch D0 gebildet und anschließend in dem Multiplizierer M mit dem vorher berechneten Teil multipliziert. Damit ist αopt berechnet.For the complete calculation of α opt , two further power estimators LS x and LS e output estimated values for the powers of the input signal and the error signal, respectively, in their known manner. Its quotient is formed by D 0 and then multiplied in the multiplier M with the previously calculated part. This calculates α opt .
Durch die Streuung der Korrelationskoeffizienten erhält man selbst bei abgeglichenem Modell aus Gleichung (9) einen 'Grundwert' für ∥Δh ^∥2. Diesen kann man näherungsweise aus einer Anzahl von Totzeitkorrelationskoeffizienten ermitteln, die theoretisch 0 sein müssen. Dazu führt man bei dem Fehlersignal e(k) eine Totzeit NTT ein und berechnet die zusätzlichen Korrelationskoeffizienten Due to the scattering of the correlation coefficients, even with a balanced model from equation (9) one obtains a 'basic value' for ∥Δh ^ ∥ 2 . This can be approximated from a number of dead-time correlation coefficients, which theoretically must be 0. For this purpose, a dead time N T T is introduced at the error signal e (k) and the additional correlation coefficients are calculated
Aus der Streuung dieser Kreuzkorrelationskoeffizienten kann auf die Streuung der anderen geschlossen werden. Subtrahiert man diesen Wert von der quadrierten Norm des ersten Vektors, dann erhält man eine verbesserte Schätzung für den Rauschübertragungsfaktor: From the scatter of these cross-correlation coefficients, one can conclude that the others are scattered. Subtracting this value from the squared norm of the first vector gives an improved estimate of the noise transfer factor:
Wegen der relativ großen Schwankungen von g ^(k) ist es zweckmäßig, einen geglätteten Wert zu verwenden.Because of the relatively large one Fluctuations of g ^ (k) it is expedient, a smoothed value to use.
Dieser
geglättete
Wert wird nach der Beziehung
- σσ
- Streuungscattering
- h ^(i)h ^ (i)
- Modellkoeffizientenmodel coefficients
- h(i)Hi)
- Koeffizienten des Echopfadescoefficients of the echo path
- NN
- Anzahl der Modellkoeffizientennumber the model coefficients
- RR
- Korrelationskoeffizientencorrelation coefficients
- xx
- Eingangssignalinput
- ee
- Fehlersignalerror signal
- NT N T
- Totzeitdead
- TT
- Zeit zwischen zwei AbtastungenTime between two scans
- ϛ(k)ς (k)
- Korrelationsfaktorcorrelation factor
Claims (5)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
DE1995125381 DE19525381B4 (en) | 1995-07-12 | 1995-07-12 | Method for controlling the step size of an adaptive-filter echo canceller |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
DE1995125381 DE19525381B4 (en) | 1995-07-12 | 1995-07-12 | Method for controlling the step size of an adaptive-filter echo canceller |
Publications (2)
Publication Number | Publication Date |
---|---|
DE19525381A1 DE19525381A1 (en) | 1997-01-16 |
DE19525381B4 true DE19525381B4 (en) | 2006-01-05 |
Family
ID=7766639
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
DE1995125381 Expired - Fee Related DE19525381B4 (en) | 1995-07-12 | 1995-07-12 | Method for controlling the step size of an adaptive-filter echo canceller |
Country Status (1)
Country | Link |
---|---|
DE (1) | DE19525381B4 (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE19831320A1 (en) * | 1998-07-13 | 2000-01-27 | Ericsson Telefon Ab L M | Digital adaptive filter for communications system, e.g. hands free communications in vehicles, has power estimation unit recursively smoothing increasing and decreasing input power asymmetrically |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4977591A (en) * | 1989-11-17 | 1990-12-11 | Nynex Corporation | Dual mode LMS nonlinear data echo canceller |
US5289539A (en) * | 1990-09-12 | 1994-02-22 | Nec Corporation | Echo canceller using impulse response estimating method |
DE4305256A1 (en) * | 1993-02-20 | 1994-08-25 | Sel Alcatel Ag | Process for improving the acoustic attenuation of electroacoustic systems |
-
1995
- 1995-07-12 DE DE1995125381 patent/DE19525381B4/en not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4977591A (en) * | 1989-11-17 | 1990-12-11 | Nynex Corporation | Dual mode LMS nonlinear data echo canceller |
US5289539A (en) * | 1990-09-12 | 1994-02-22 | Nec Corporation | Echo canceller using impulse response estimating method |
DE4305256A1 (en) * | 1993-02-20 | 1994-08-25 | Sel Alcatel Ag | Process for improving the acoustic attenuation of electroacoustic systems |
Non-Patent Citations (9)
Title |
---|
Yasukawa, H.: Acoustic Echo Cancelling. In: Electronic Letters, 16th July 1992, Vol. 28, No. 15, S. 1403,1404 |
Huhn, Tilmann; Jentschel, Hans-Joachim: Kombination von Geräuschreduktion und Echokompensation beim Freisprechen. In: Nachrichtentech., Elektron., Berlin 43, 1993, H.6, S.274-280 * |
Huhn, Tilmann;Jentschel, Hans-Joachim: Kombination von Geräuschreduktion und Echokompensation beim Freisprechen. In: Nachrichtentech., Elektron., Berlin 43, 1993, H.6, S.274-280 |
Murano, Kazuo, et.al.: Echo Cancellation and Applications. In: IEEE Communications Magazine, Jan. 1990, S. 49-55 * |
Proakis, John G.: Digital Communications,2nd ed., New York, McGraw-Hill Book Company, 1989, S. 683 bis 691, ISBN 0-07-100269-3 * |
R. Frenzel, Fortschr.-Ber. VDI-Reihe 10, Nr. 228, Düsseldorf 1992 |
S. Yamamoto, S. Kitayama: "An Adaptive Echo Can- celler with Variable Step Gain Method". In: Trans. of the IECE of Japan, Vol. E 65, Nr. 1, Jan. 1982, S. 1-8 |
Till, Reinhard: Adaptive Sprecherecho- Kompensation in Modems für die Duplex- Datenübertragung im Fernsprechnetz. In: Frequenz 37, 1983, H.6, S. 145-154 * |
Wehrmann, R., Poltmann, R., u.a.: Signalverarbeitungsverfahren zur Verbesserung der Sprachkommunikation über Freisprecheinrichtungen. In: Der Fernmelde Ingenieur, Teil I: H. 10/94, S. 1-30, Teil II: H. 1/95, S. 4-36 * |
Also Published As
Publication number | Publication date |
---|---|
DE19525381A1 (en) | 1997-01-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
DE69023697T2 (en) | Coefficient generator for digital filters. | |
DE2914945C2 (en) | ||
DE69114999T2 (en) | Echo canceller with improved intercom detector. | |
EP0371567B1 (en) | Echo canceller | |
DE3887459T2 (en) | Digital signal receiving circuit with baud sampling phase control by a power of sampled signals. | |
DE60101986T2 (en) | Channel estimation sequence and method for estimating a transmission channel using such a sequence for channel estimation | |
DE2207141C3 (en) | Circuit arrangement for suppressing unwanted speech signals by means of a predictive filter | |
DE1272978B (en) | Circuit arrangement for correcting the distortions of digital communication signals caused by transmission media of limited frequency bandwidth | |
DE3405010A1 (en) | DEVICE FOR GENERATING A DELAY ESTIMATE FOR AN ECHO CANCELING DEVICE | |
DE3113394C2 (en) | Method and device for adaptive echo cancellation in a two-wire full duplex transmission | |
DE19806015C2 (en) | Process for improving acoustic attenuation in hands-free systems | |
DE19525381B4 (en) | Method for controlling the step size of an adaptive-filter echo canceller | |
DE1787007A1 (en) | VARIABLE EQUALIZER | |
EP1051822B1 (en) | Decision feedback equalizer with noise prediction | |
DE69306893T2 (en) | Method for determining the transmittance of a filter circuit for converting the impulse response of a filter into an impulse response with a minimal phase and filter for carrying out the method | |
DE19525382B4 (en) | Method for controlling the step size of an adaptive-filter echo canceller | |
EP0402519B1 (en) | Method and arrangement for treating the dynamic range of an adaptive recursive network for processing time-discrete signals | |
DE10250361A1 (en) | Device and method for processing pilot symbols for channel estimation using adaptive low-pass filtering | |
EP0092035B1 (en) | Method and circuit arrangement for the effectuation of the method of adaptive echo-cancelling in terminals for duplex data transmission over two-wire lines | |
DE69633454T2 (en) | ADAPTIVE ECHOUNTERPRESSION PROCESS | |
EP0231487A2 (en) | Car radio with nonrecursive digital filter-shaped equalizer | |
DE69734888T2 (en) | Circuit configuration for adapting a multi-speed modem to a line and corresponding adaptation method | |
EP0992151B1 (en) | Method and device for suppressing echo in a hands free device such as a telephone | |
EP0881814A1 (en) | Method for determining the step size Alpha for adjusting the convergence speed in the NLMS algorithm | |
DE2063271C (en) | Adaptive compensation echo canceller |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
OM8 | Search report available as to paragraph 43 lit. 1 sentence 1 patent law | ||
8110 | Request for examination paragraph 44 | ||
8339 | Ceased/non-payment of the annual fee |