JP3381112B2 - エコー除去装置 - Google Patents
エコー除去装置Info
- Publication number
- JP3381112B2 JP3381112B2 JP05002595A JP5002595A JP3381112B2 JP 3381112 B2 JP3381112 B2 JP 3381112B2 JP 05002595 A JP05002595 A JP 05002595A JP 5002595 A JP5002595 A JP 5002595A JP 3381112 B2 JP3381112 B2 JP 3381112B2
- Authority
- JP
- Japan
- Prior art keywords
- echo
- filter
- signal
- value
- tap coefficient
- 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
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/015—Reducing echo effects
-
- 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
- H04R—LOUDSPEAKERS, MICROPHONES, GRAMOPHONE PICK-UPS OR LIKE ACOUSTIC ELECTROMECHANICAL TRANSDUCERS; DEAF-AID SETS; PUBLIC ADDRESS SYSTEMS
- H04R3/00—Circuits for transducers, loudspeakers or microphones
- H04R3/02—Circuits for transducers, loudspeakers or microphones for preventing acoustic reaction, i.e. acoustic oscillatory feedback
Landscapes
- Engineering & Computer Science (AREA)
- Signal Processing (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Otolaryngology (AREA)
- Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Computer Networks & Wireless Communication (AREA)
- Cable Transmission Systems, Equalization Of Radio And Reduction Of Echo (AREA)
- Telephone Function (AREA)
- Filters That Use Time-Delay Elements (AREA)
Description
【０００１】
【産業上の利用分野】本発明は、エコー除去装置に関
し、特に、携帯電話等の小型音声通信端末のスピーカか
らマイクへの音声の廻り込みによるエコーを軽減するた
めのエコー除去装置に関する。
し、特に、携帯電話等の小型音声通信端末のスピーカか
らマイクへの音声の廻り込みによるエコーを軽減するた
めのエコー除去装置に関する。
【０００２】
【従来の技術】携帯電話等のように、音声通信端末の小
型化に伴い、受話用のスピーカから送話用のマイクへの
音声の廻り込みによるエコーの影響が無視できなくなっ
てきている。このような送受話側での音声の廻り込みに
よるエコーを除去するために、例えば図５に示すような
エコー除去装置あるいはエコーキャンセラが知られてい
る。
型化に伴い、受話用のスピーカから送話用のマイクへの
音声の廻り込みによるエコーの影響が無視できなくなっ
てきている。このような送受話側での音声の廻り込みに
よるエコーを除去するために、例えば図５に示すような
エコー除去装置あるいはエコーキャンセラが知られてい
る。
【０００３】この図５において、端子１１には通信相手
から伝送されて受信されスピーカ１２に送られるスピー
カ出力信号ｘ(k) が入力されている。ここでｋは、離散
信号のサンプル番号あるいは時間位置を示す。また、マ
イクロホン１３により収音されて電気信号に変換された
マイク入力信号ｙ(k) が、エコー除去のための減算器１
４に供給され、フィルタ回路１５からの擬似エコー信号
が減算された後、エコー除去信号あるいは残留エコー信
号ｅ(k) として端子１６より取り出される。スピーカ１
２とマイクロホン１３とは、携帯電話の送受話器として
近接して配置されている。
から伝送されて受信されスピーカ１２に送られるスピー
カ出力信号ｘ(k) が入力されている。ここでｋは、離散
信号のサンプル番号あるいは時間位置を示す。また、マ
イクロホン１３により収音されて電気信号に変換された
マイク入力信号ｙ(k) が、エコー除去のための減算器１
４に供給され、フィルタ回路１５からの擬似エコー信号
が減算された後、エコー除去信号あるいは残留エコー信
号ｅ(k) として端子１６より取り出される。スピーカ１
２とマイクロホン１３とは、携帯電話の送受話器として
近接して配置されている。
【０００４】フィルタ回路１５には例えばいわゆるＦＩ
Ｒ（有限インパルス応答）フィルタが用いられており、
フィルタ回路１５のフィルタ係数あるいはタップ係数
は、タップ係数推定回路１７により制御されるようにな
っている。タップ係数推定回路１７は、入力信号である
スピーカ出力信号ｘ(k) や、出力信号である残留エコー
信号ｅ(k) に基づいて、適当な適応アルゴリズムにより
フィルタのタップ係数を逐次更新する。この適応アルゴ
リズムとしては、いわゆる学習同定法あるいは正規化Ｌ
ＭＳ（Least Mean Square ）法が広く用いられている。
Ｒ（有限インパルス応答）フィルタが用いられており、
フィルタ回路１５のフィルタ係数あるいはタップ係数
は、タップ係数推定回路１７により制御されるようにな
っている。タップ係数推定回路１７は、入力信号である
スピーカ出力信号ｘ(k) や、出力信号である残留エコー
信号ｅ(k) に基づいて、適当な適応アルゴリズムにより
フィルタのタップ係数を逐次更新する。この適応アルゴ
リズムとしては、いわゆる学習同定法あるいは正規化Ｌ
ＭＳ（Least Mean Square ）法が広く用いられている。
【０００５】ここで、フィルタのタップ数をＬとし、上
記スピーカ出力信号をフィルタのタップ入力信号とし
て、 Ｘ(k)＝（x(k),x(k-1),…,x(k-L+1)） のベクトルで表し、フィルタのタップ係数を
記スピーカ出力信号をフィルタのタップ入力信号とし
て、 Ｘ(k)＝（x(k),x(k-1),…,x(k-L+1)） のベクトルで表し、フィルタのタップ係数を
【０００６】
【数２】
【０００７】のベクトルで表すとき、ＦＩＲフィルタの
フィルタ回路１５により得られる擬似エコー信号は、上
記タップ入力にタップ係数が乗じられて、
フィルタ回路１５により得られる擬似エコー信号は、上
記タップ入力にタップ係数が乗じられて、
【０００８】
【数３】
【０００９】により求められる。この擬似エコー信号が
上記マイク入力信号ｙ(k) より減算されることにより、
残留エコー信号ｅ(k) が、
上記マイク入力信号ｙ(k) より減算されることにより、
残留エコー信号ｅ(k) が、
【００１０】
【数４】
【００１１】として得られ、これがエコー除去出力信号
あるいはエコーキャンセラ出力信号となる。タップ係数
推定回路１７においては、これらの残留エコー信号ｅ
(k) 及びタップ入力信号Ｘ(k) に基づいて、上記タップ
係数が推定され更新される。
あるいはエコーキャンセラ出力信号となる。タップ係数
推定回路１７においては、これらの残留エコー信号ｅ
(k) 及びタップ入力信号Ｘ(k) に基づいて、上記タップ
係数が推定され更新される。
【００１２】このタップ係数推定回路１７として、上記
学習同定法のアルゴリズムを用いる場合の具体例を図６
に示す。
学習同定法のアルゴリズムを用いる場合の具体例を図６
に示す。
【００１３】この図６において、タップ係数推定回路１
７は、２乗和計算回路１８とタップ係数更新回路１９と
から成り、２乗和計算回路１８は、上記スピーカ出力信
号のタップ入力であるＸ(k)＝（x(k),x(k-1),…,x(k-L+
1)）より、その２乗和
７は、２乗和計算回路１８とタップ係数更新回路１９と
から成り、２乗和計算回路１８は、上記スピーカ出力信
号のタップ入力であるＸ(k)＝（x(k),x(k-1),…,x(k-L+
1)）より、その２乗和
【００１４】
【数５】
【００１５】を計算する。この２乗和は入力信号電力に
相当するものである。タップ係数更新回路１９は、求め
られた２乗和と、上記残留エコー信号ｅ(k) 及びタップ
入力信号Ｘ(k) をもとに、次の（３）式によりタップ係
数を更新する。
相当するものである。タップ係数更新回路１９は、求め
られた２乗和と、上記残留エコー信号ｅ(k) 及びタップ
入力信号Ｘ(k) をもとに、次の（３）式によりタップ係
数を更新する。
【００１６】
【数６】
【００１７】この（３）式中のμは定数であり、０＜μ
＜２の値をとる。この（３）式は、成分毎には、
＜２の値をとる。この（３）式は、成分毎には、
【００１８】
【数７】
【００１９】と表される。
【００２０】
【発明が解決しようとする課題】ところで、フィルタに
対するタップ入力信号Ｘ(k) は音声信号であり、図７に
示すように一般にピッチ周期構造を持つ周期信号であ
る。例えばサンプリング周波数が８ｋHzの場合、上記ピ
ッチ周期は男声で５０〜９０サンプル、女声で２０〜３
５サンプルとなる。このため、フィルタのタップ数Ｌが
上記ピッチ周期のサンプル数に比べて充分に大きくない
場合には、上記（２）式によって、図７の区間ａやｂに
示す範囲での２乗和がとられることになり、求められた
入力信号電力は、図８に示すように大きく変動すること
になる。
対するタップ入力信号Ｘ(k) は音声信号であり、図７に
示すように一般にピッチ周期構造を持つ周期信号であ
る。例えばサンプリング周波数が８ｋHzの場合、上記ピ
ッチ周期は男声で５０〜９０サンプル、女声で２０〜３
５サンプルとなる。このため、フィルタのタップ数Ｌが
上記ピッチ周期のサンプル数に比べて充分に大きくない
場合には、上記（２）式によって、図７の区間ａやｂに
示す範囲での２乗和がとられることになり、求められた
入力信号電力は、図８に示すように大きく変動すること
になる。
【００２１】このように、２乗和あるいは入力信号電力
が変動すると、上記（３）式によって更新されるタップ
係数が変動することになり、安定したエコー除去特性が
得られないことになる。
が変動すると、上記（３）式によって更新されるタップ
係数が変動することになり、安定したエコー除去特性が
得られないことになる。
【００２２】これは、特に、上述したような携帯電話な
どのスピーカ１２からマイクロホン１３への廻り込みに
よるエコー除去を対象とする比較的小型のエコー除去装
置の場合、例えばフィルタのタップ数が５０以下の場合
においては、大きな問題となる。
どのスピーカ１２からマイクロホン１３への廻り込みに
よるエコー除去を対象とする比較的小型のエコー除去装
置の場合、例えばフィルタのタップ数が５０以下の場合
においては、大きな問題となる。
【００２３】本発明は、上述した実情に鑑みてなされた
ものであり、上記入力信号電力の変動を抑えて、安定し
たエコー除去特性あるいはエコー消去特性が得られるよ
うなエコー除去装置の提供を目的とするものである。
ものであり、上記入力信号電力の変動を抑えて、安定し
たエコー除去特性あるいはエコー消去特性が得られるよ
うなエコー除去装置の提供を目的とするものである。
【００２４】
【課題を解決するための手段】本発明は、発音手段から
の音が近接して配置された収音手段に廻り込むことによ
り生じるエコーを除去するエコー除去装置において、上
記発音手段に供給される発音信号に基づいて上記収音手
段に廻り込むエコー成分を推定した擬似エコー信号を出
力するフィルタ手段と、上記収音手段からの収音信号よ
り上記フィルタ手段からの擬似エコー信号を減算して出
力する減算手段と、上記発音信号を平滑化した平滑値を
用いて上記フィルタ手段のフィルタ係数を推定して更新
するフィルタ係数推定手段とを有し、上記フィルタ係数
推定手段は、上記フィルタ手段のタップ数Ｌより大きな
数Ｎについて過去Ｎ個のサンプルの上記発音信号の電力
のＬ個分に相当する電力を上記平滑値として用いること
により、上述の課題を解決する。
の音が近接して配置された収音手段に廻り込むことによ
り生じるエコーを除去するエコー除去装置において、上
記発音手段に供給される発音信号に基づいて上記収音手
段に廻り込むエコー成分を推定した擬似エコー信号を出
力するフィルタ手段と、上記収音手段からの収音信号よ
り上記フィルタ手段からの擬似エコー信号を減算して出
力する減算手段と、上記発音信号を平滑化した平滑値を
用いて上記フィルタ手段のフィルタ係数を推定して更新
するフィルタ係数推定手段とを有し、上記フィルタ係数
推定手段は、上記フィルタ手段のタップ数Ｌより大きな
数Ｎについて過去Ｎ個のサンプルの上記発音信号の電力
のＬ個分に相当する電力を上記平滑値として用いること
により、上述の課題を解決する。
【００２５】ここで、上記フィルタ係数推定手段は、学
習同定法によりフィルタのタップ係数の推定を行い、こ
の学習同定法によるタップ係数推定に用いられる式中の
２乗和の値として、上記発音信号の過去Ｎ個のサンプル
の２乗和の内のＬ個分相当の２乗和の値を上記平滑値と
して用いることが挙げられる。
習同定法によりフィルタのタップ係数の推定を行い、こ
の学習同定法によるタップ係数推定に用いられる式中の
２乗和の値として、上記発音信号の過去Ｎ個のサンプル
の２乗和の内のＬ個分相当の２乗和の値を上記平滑値と
して用いることが挙げられる。
【００２６】
【作用】フィルタ推定の際に平滑値を用いることによ
り、フィルタ係数あるいはタップ係数の変動が抑えら
れ、安定したエコー除去特性あるいはエコー消去特性が
得られる。
り、フィルタ係数あるいはタップ係数の変動が抑えら
れ、安定したエコー除去特性あるいはエコー消去特性が
得られる。
【００２７】
【実施例】以下、本発明に係る好ましい実施例につい
て、図面を参照しながら説明する。
て、図面を参照しながら説明する。
【００２８】図１は、本発明に係るエコー除去装置の一
実施例の概略構成を示し、端子１１には通信相手から伝
送されて受信され発音手段であるスピーカ１２に送られ
る発音信号としてのスピーカ出力信号ｘ(k) が入力され
ている。また、発音手段であるスピーカ１２に近接して
配置された収音手段であるマイクロホン１３により収音
されて電気信号に変換された収音信号としてのマイク入
力信号ｙ(k) が、エコー除去のための減算器１４に供給
され、フィルタ回路１５からの擬似エコー信号が減算さ
れた後、エコー除去信号あるいは残留エコー信号ｅ(k)
として端子１６より取り出される。
実施例の概略構成を示し、端子１１には通信相手から伝
送されて受信され発音手段であるスピーカ１２に送られ
る発音信号としてのスピーカ出力信号ｘ(k) が入力され
ている。また、発音手段であるスピーカ１２に近接して
配置された収音手段であるマイクロホン１３により収音
されて電気信号に変換された収音信号としてのマイク入
力信号ｙ(k) が、エコー除去のための減算器１４に供給
され、フィルタ回路１５からの擬似エコー信号が減算さ
れた後、エコー除去信号あるいは残留エコー信号ｅ(k)
として端子１６より取り出される。
【００２９】このエコー除去装置は例えば小型の携帯電
話等の音声通信端末に用いられるものであり、スピーカ
１２とマイクロホン１３とは、例えば携帯電話の送受話
器として近接して配置されている。
話等の音声通信端末に用いられるものであり、スピーカ
１２とマイクロホン１３とは、例えば携帯電話の送受話
器として近接して配置されている。
【００３０】フィルタ回路１５には例えばいわゆるＦＩ
Ｒ（有限インパルス応答）フィルタが用いられており、
フィルタ回路１５のフィルタ係数あるいはタップ係数
は、フィルタ係数推定手段としてのタップ係数推定回路
２０により制御されるようになっている。タップ係数推
定回路２０は、エコー除去装置への入力信号であるスピ
ーカ出力信号ｘ(k) や、出力信号である残留エコー信号
ｅ(k) に基づいて、適当な適応アルゴリズムによりフィ
ルタのタップ係数を逐次更新するものであり、この適応
アルゴリズムとしては、いわゆる学習同定法あるいは正
規化ＬＭＳ（Least Mean Square ）法が広く用いられて
いる。
Ｒ（有限インパルス応答）フィルタが用いられており、
フィルタ回路１５のフィルタ係数あるいはタップ係数
は、フィルタ係数推定手段としてのタップ係数推定回路
２０により制御されるようになっている。タップ係数推
定回路２０は、エコー除去装置への入力信号であるスピ
ーカ出力信号ｘ(k) や、出力信号である残留エコー信号
ｅ(k) に基づいて、適当な適応アルゴリズムによりフィ
ルタのタップ係数を逐次更新するものであり、この適応
アルゴリズムとしては、いわゆる学習同定法あるいは正
規化ＬＭＳ（Least Mean Square ）法が広く用いられて
いる。
【００３１】タップ係数推定回路２０は、電力平滑値計
算回路２１及びタップ係数更新回路２２を有して構成さ
れている。電力平滑値計算回路２１では、従来のタップ
数Ｌ個のサンプルの２乗和の代わりに、上記発音信号を
平滑化した平滑値、すなわちエコー除去装置への入力信
号である上記スピーカ出力信号ｘ(k) の平滑された電力
値である入力信号電力平滑値Ｐｘ(k)を用いている。こ
の電力平滑値計算回路２１からの平滑値あるいは入力信
号電力平滑値Ｐｘ(k)がタップ係数更新回路２２に送ら
れる。タップ係数更新回路２２では、上記スピーカ出力
信号ｘ(k) のタップ数Ｌ個分のベクトルであるＸ(k)＝
（x(k),x(k-1),…,x(k-L+1)）と、上記残留エコー信号
ｅ(k) と、電力平滑値計算回路２１からの入力信号電力
平滑値Ｐｘ(k)とに基づいて、タップ係数の更新が行わ
れる。
算回路２１及びタップ係数更新回路２２を有して構成さ
れている。電力平滑値計算回路２１では、従来のタップ
数Ｌ個のサンプルの２乗和の代わりに、上記発音信号を
平滑化した平滑値、すなわちエコー除去装置への入力信
号である上記スピーカ出力信号ｘ(k) の平滑された電力
値である入力信号電力平滑値Ｐｘ(k)を用いている。こ
の電力平滑値計算回路２１からの平滑値あるいは入力信
号電力平滑値Ｐｘ(k)がタップ係数更新回路２２に送ら
れる。タップ係数更新回路２２では、上記スピーカ出力
信号ｘ(k) のタップ数Ｌ個分のベクトルであるＸ(k)＝
（x(k),x(k-1),…,x(k-L+1)）と、上記残留エコー信号
ｅ(k) と、電力平滑値計算回路２１からの入力信号電力
平滑値Ｐｘ(k)とに基づいて、タップ係数の更新が行わ
れる。
【００３２】すなわち、例えば学習同定法のアルゴリズ
ムにおいて、前述した（３）式、すなわち、
ムにおいて、前述した（３）式、すなわち、
【００３３】
【数８】
【００３４】の従来の係数更新の式中の分母の２乗和あ
るいは入力信号電力は、
るいは入力信号電力は、
【００３５】
【数９】
【００３６】であればよい。そこで、上記（３）式の分
母として、入力信号であるスピーカ出力信号ｘ(k) の系
列のタップ数Ｌ個分の２乗和を用いる代わりに、何らか
の方法で平滑化された入力信号電力値を用いることによ
り、入力信号レベルの変動を除去することができる。
母として、入力信号であるスピーカ出力信号ｘ(k) の系
列のタップ数Ｌ個分の２乗和を用いる代わりに、何らか
の方法で平滑化された入力信号電力値を用いることによ
り、入力信号レベルの変動を除去することができる。
【００３７】この平滑化された入力信号電力の具体例を
図２に示す。この図２において、曲線ａは上記（３）式
の分母の２乗和により計算される入力信号電力を、また
曲線ｂは上記平滑化された入力信号電力である入力信号
電力平滑値Ｐｘ(k)をそれぞれ示している。この図２か
らも明らかなように、変動の大きな曲線ａの２乗和が、
平滑することで曲線ｂに示すように変動が除去されてい
ることが分かる。
図２に示す。この図２において、曲線ａは上記（３）式
の分母の２乗和により計算される入力信号電力を、また
曲線ｂは上記平滑化された入力信号電力である入力信号
電力平滑値Ｐｘ(k)をそれぞれ示している。この図２か
らも明らかなように、変動の大きな曲線ａの２乗和が、
平滑することで曲線ｂに示すように変動が除去されてい
ることが分かる。
【００３８】このような入力信号電力平滑値Ｐｘ(k)を
用いて、タップ係数更新回路２２におけるタップ係数更
新は、
用いて、タップ係数更新回路２２におけるタップ係数更
新は、
【００３９】
【数１０】
【００４０】の式により行う。この（５）式により計算
されたタップ係数は、変動が除去されたものとなってい
る。
されたタップ係数は、変動が除去されたものとなってい
る。
【００４１】このように、タップ数Ｌが少なくて入力信
号電力である２乗和が変動し、安定したエコー消去特性
が得られない場合でも、入力信号電力を平滑し、この平
滑化された入力信号電力あるいは入力信号電力平滑値Ｐ
ｘ(k)を用いることにより、タップ係数推定の際の変動
による悪影響を防止し、安定したエコー除去特性あるい
はエコー消去特性を実現できる。また、このときのタッ
プ係数推定に要する計算量も従来の学習同定法に要する
計算量と同程度で済む。
号電力である２乗和が変動し、安定したエコー消去特性
が得られない場合でも、入力信号電力を平滑し、この平
滑化された入力信号電力あるいは入力信号電力平滑値Ｐ
ｘ(k)を用いることにより、タップ係数推定の際の変動
による悪影響を防止し、安定したエコー除去特性あるい
はエコー消去特性を実現できる。また、このときのタッ
プ係数推定に要する計算量も従来の学習同定法に要する
計算量と同程度で済む。
【００４２】次に、図１の実施例のタップ係数推定回路
２０内の電力平滑値計算回路２１の具体的ないくつかの
構成例について説明する。
２０内の電力平滑値計算回路２１の具体的ないくつかの
構成例について説明する。
【００４３】先ず、第１の具体例としては、エコー除去
装置への入力信号である上記発音信号としてのスピーカ
出力信号ｘ(k) の系列の電力を計算する際に、過去Ｌ個
分のサンプルの代わりに、より多くの個数Ｎ（Ｎ＞Ｌ）
の過去の入力信号のサンプルの２乗和を求めＬ個分の平
均を計算することにより、入力信号電力の変動を除去す
るものが挙げられる。
装置への入力信号である上記発音信号としてのスピーカ
出力信号ｘ(k) の系列の電力を計算する際に、過去Ｌ個
分のサンプルの代わりに、より多くの個数Ｎ（Ｎ＞Ｌ）
の過去の入力信号のサンプルの２乗和を求めＬ個分の平
均を計算することにより、入力信号電力の変動を除去す
るものが挙げられる。
【００４４】すなわち、このときの入力信号電力平滑値
Ｐｘ(k)は、
Ｐｘ(k)は、
【００４５】
【数１１】
【００４６】により求められる。これは、入力信号系列
の過去のＮ個分のサンプルの２乗和についてのＬサンプ
ル分の平均値あるいは正規化値を求めていることにな
る。この（６）式により求められた入力信号電力平滑値
Ｐｘ(k)を、上記（５）式に代入することにより、現在
のタップ係数から次のタップ係数を推定し更新する。
の過去のＮ個分のサンプルの２乗和についてのＬサンプ
ル分の平均値あるいは正規化値を求めていることにな
る。この（６）式により求められた入力信号電力平滑値
Ｐｘ(k)を、上記（５）式に代入することにより、現在
のタップ係数から次のタップ係数を推定し更新する。
【００４７】次に、第２の具体例として、入力信号電力
平滑値Ｐｘ(k)を指数平滑により求める例について、図
３を参照しながら説明する。
平滑値Ｐｘ(k)を指数平滑により求める例について、図
３を参照しながら説明する。
【００４８】この図３に示す第２の具体例では、上記入
力信号電力平滑値Ｐｘ(k)を、
力信号電力平滑値Ｐｘ(k)を、
【００４９】
【数１２】
【００５０】の式により求めている。ここでλは、０＜
λ＜１の定数であり、この（７）式により計算される入
力信号電力平滑値Ｐｘ(k)は、過去の入力信号サンプル
の２乗値を時定数１／（１−λ）にて指数平滑化したも
のとなる。実際に、上記（７）式は、次の漸化式により
表現できる。
λ＜１の定数であり、この（７）式により計算される入
力信号電力平滑値Ｐｘ(k)は、過去の入力信号サンプル
の２乗値を時定数１／（１−λ）にて指数平滑化したも
のとなる。実際に、上記（７）式は、次の漸化式により
表現できる。
【００５１】
【数１３】
【００５２】図３は、この（８）式を実現するようなタ
ップ係数推定回路２０の構成例を示すものである。
ップ係数推定回路２０の構成例を示すものである。
【００５３】この図３において、タップ係数推定回路２
０の電力平滑値計算回路２１には、エコー除去装置への
入力信号である上記スピーカ出力信号ｘ(k) が供給され
ており、２乗値計算回路２３により２乗値｜ｘ(k)｜２
が計算されて、加算器２４に送られ、加算器２４からの
出力は、１サンプル遅延素子２５及び係数乗算器２６を
介して加算器２４に帰還されている。ここで、加算器２
４からは、入力信号電力平滑値Ｐｘ(k)が出力される。
すなわち、１サンプル遅延素子２５からは１サンプル遅
れた入力信号電力平滑値Ｐｘ(k-1)が出力され、係数乗
算器２６に送られてλが乗算されてλＰｘ(k-1)とな
り、これが加算器２４で｜ｘ(k)｜２と加算されること
により、上記（８）式の入力信号電力平滑値Ｐｘ(k)が
求められる。この加算器２４からの出力である入力信号
電力平滑値Ｐｘ(k)は、タップ係数更新回路２２に送ら
れる。
０の電力平滑値計算回路２１には、エコー除去装置への
入力信号である上記スピーカ出力信号ｘ(k) が供給され
ており、２乗値計算回路２３により２乗値｜ｘ(k)｜２
が計算されて、加算器２４に送られ、加算器２４からの
出力は、１サンプル遅延素子２５及び係数乗算器２６を
介して加算器２４に帰還されている。ここで、加算器２
４からは、入力信号電力平滑値Ｐｘ(k)が出力される。
すなわち、１サンプル遅延素子２５からは１サンプル遅
れた入力信号電力平滑値Ｐｘ(k-1)が出力され、係数乗
算器２６に送られてλが乗算されてλＰｘ(k-1)とな
り、これが加算器２４で｜ｘ(k)｜２と加算されること
により、上記（８）式の入力信号電力平滑値Ｐｘ(k)が
求められる。この加算器２４からの出力である入力信号
電力平滑値Ｐｘ(k)は、タップ係数更新回路２２に送ら
れる。
【００５４】タップ係数更新回路２２では、このような
電力平滑値計算回路２１からの入力信号電力平滑値Ｐｘ
(k)を用いて上記（５）式の計算を行うことにより、フ
ィルタ回路１５のタップ係数を更新する。
電力平滑値計算回路２１からの入力信号電力平滑値Ｐｘ
(k)を用いて上記（５）式の計算を行うことにより、フ
ィルタ回路１５のタップ係数を更新する。
【００５５】なお、上記（５）式のタップ係数ベクトル
の各成分は、
の各成分は、
【００５６】
【数１４】
【００５７】となる。
【００５８】ここで図４は、このような第２の具体例の
アルゴリズムのタップ係数推定を行ったエコー除去装置
のエコー消去特性（曲線ａ）と、従来のアルゴリズムを
用いた場合のエコー消去特性（曲線ｂ）とを比較して示
すものである。この図４の例では、タップ数Ｌが４０の
適応フィルタを用いた場合を示している。
アルゴリズムのタップ係数推定を行ったエコー除去装置
のエコー消去特性（曲線ａ）と、従来のアルゴリズムを
用いた場合のエコー消去特性（曲線ｂ）とを比較して示
すものである。この図４の例では、タップ数Ｌが４０の
適応フィルタを用いた場合を示している。
【００５９】この図４の縦軸のエコー消去特性ERLEは、
【００６０】
【数１５】
【００６１】により定義されるものであり、大きいほど
良好な特性を示すものである。なお、この（１０）式中
のＥ［］は、［］内の期待値を返す関数である。
良好な特性を示すものである。なお、この（１０）式中
のＥ［］は、［］内の期待値を返す関数である。
【００６２】この図４からも明らかなように、本発明の
実施例によるエコー除去装置の方が従来のものよりも安
定した良好なエコー消去特性が得られていことが分か
る。
実施例によるエコー除去装置の方が従来のものよりも安
定した良好なエコー消去特性が得られていことが分か
る。
【００６３】さらに、このような実施例のタップ係数推
定に要する計算量は、従来の学習同定法のアルゴリズム
の場合の計算量とほぼ同程度で済む。
定に要する計算量は、従来の学習同定法のアルゴリズム
の場合の計算量とほぼ同程度で済む。
【００６４】従って、適応フィルタのタップ数が数十程
度と小さい場合、特に小型の携帯電話のエコーキャンセ
ラの場合に用いて有用である。
度と小さい場合、特に小型の携帯電話のエコーキャンセ
ラの場合に用いて有用である。
【００６５】なお、本発明は上述した実施例のみに限定
されるものではなく、例えば、携帯電話の他にも種々の
音声通信端末に適用でき、発音手段や収音手段もスピー
カやマイクロホンに限定されない。また、フィルタ手段
はＦＩＲフィルタ以外にも種々の適応フィルタを使用で
き、フィルタ係数の推定も学習同定法に限定されず、種
々の適応アルゴリズムを使用できる。
されるものではなく、例えば、携帯電話の他にも種々の
音声通信端末に適用でき、発音手段や収音手段もスピー
カやマイクロホンに限定されない。また、フィルタ手段
はＦＩＲフィルタ以外にも種々の適応フィルタを使用で
き、フィルタ係数の推定も学習同定法に限定されず、種
々の適応アルゴリズムを使用できる。
【００６６】
【発明の効果】本発明に係るエコー除去装置によれば、
発音手段に供給される発音信号に基づいてフィルタ手段
により収音手段に廻り込むエコー成分を推定して擬似エ
コー信号を出力し、収音手段からの収音信号より擬似エ
コー信号を減算する際に、フィルタ手段のフィルタ係数
推定を、発音信号を平滑化した平滑値を用いて行うと共
に、フィルタ手段のタップ数Ｌより大きな数Ｎについて
過去Ｎ個のサンプルの上記発音信号の電力のＬ個分に相
当する電力を上記平滑値として用いることにより、フィ
ルタ係数あるいはタップ係数の変動が抑えられ、安定し
たエコー除去特性あるいはエコー消去特性を有するエコ
ー除去装置を提供できる。
発音手段に供給される発音信号に基づいてフィルタ手段
により収音手段に廻り込むエコー成分を推定して擬似エ
コー信号を出力し、収音手段からの収音信号より擬似エ
コー信号を減算する際に、フィルタ手段のフィルタ係数
推定を、発音信号を平滑化した平滑値を用いて行うと共
に、フィルタ手段のタップ数Ｌより大きな数Ｎについて
過去Ｎ個のサンプルの上記発音信号の電力のＬ個分に相
当する電力を上記平滑値として用いることにより、フィ
ルタ係数あるいはタップ係数の変動が抑えられ、安定し
たエコー除去特性あるいはエコー消去特性を有するエコ
ー除去装置を提供できる。
【００６７】しかも、フィルタ係数あるいはタップ係数
の推定に要する計算量は、従来の学習同定法の計算量と
ほぼ同程度で済む。
の推定に要する計算量は、従来の学習同定法の計算量と
ほぼ同程度で済む。
【図１】本発明の一実施例となるエコー除去装置の概略
構成を示すブロック図である。
構成を示すブロック図である。
【図２】入力信号電力の平滑化を説明するための図であ
る。
る。
【図３】本発明の実施例のタップ係数推定回路の具体例
を示すブロック回路図である。
を示すブロック回路図である。
【図４】エコー消去特性を示すグラフである。
【図５】一般的なエコー除去装置の概略構成を示すブロ
ック図である。
ック図である。
【図６】図５のエコー除去装置に用いられるタップ係数
推定回路の具体例を示すブロック図である。
推定回路の具体例を示すブロック図である。
【図７】入力信号のピッチ周期と入力信号電力の計算区
間との関係を示す図である。
間との関係を示す図である。
【図８】フィルタのタップ数が少ない場合の入力信号電
力の変動を示す図である。
力の変動を示す図である。
１２ スピーカ、 １３ マイクロホン、 １４ 減算
器、 １５ フィルタ、 ２０ タップ係数推定回路、
２１ 電力平滑値計算回路 、 ２２ タップ係数更
新回路
器、 １５ フィルタ、 ２０ タップ係数推定回路、
２１ 電力平滑値計算回路 、 ２２ タップ係数更
新回路
─────────────────────────────────────────────────────
フロントページの続き
(56)参考文献 特開 平４−119713（ＪＰ，Ａ)
特開 平６−338828（ＪＰ，Ａ)
特開 平６−97711（ＪＰ，Ａ)
特開 昭51−64815（ＪＰ，Ａ)
Ｋｅｎｓａｋｕ Ｆｕｊｉｉ ａｎｄ
Ｊｕｒｏ Ｏｈｇａ，“Ａ Ｆａｓｔ
Ａｄａｐｔｉｖｅ Ａｌｇｏｒｉｔｈ
ｍ Ｓｕｉｔａｂｌｅ ｆｏｒ Ａｃｏ
ｕｓｔｉｃ Ｅｃｈｏ Ｃａｎｃｅｌｌ
ｅｒ”，ＩＥＩＣＥ Ｔｒａｎｓａｃｔ
ｉｏｎｓ ｏｎ Ｆｕｎｄａｍｅｎｔａ
ｌｓ ｏｆ Ｅｌｅｃｔｒｏｎｉｃｓ,
Ｃｏｍｍｕｎｉｃａｔｉｏｎｓ ａｎｄ
Ｃｏｍｐｕｔｅｒ Ｓｃｉｅｎｃｅ
ｓ，日本，1992年，Ｖｏｌ．Ｅ75−Ａ,
Ｎｏ．11，ｐｐ．1509−1515
Ｋｅｎｓａｋｕ Ｆｕｊｉｉ ａｎｄ
Ｊｕｒｏ Ｏｈｇａ，“Ａ ＦＡＳＴ
ＡＤＡＰＴＩＶＥ ＡＬＧＯＲＩＴＨ
Ｍ ＳＵＩＴＡＢＬＥ ＦＯＲ ＡＣＯ
ＵＳＴＩＣ ＥＣＨＯ ＣＡＮＣＥＬＬ
ＥＲ”，Ｄｅｎｓｈｉ Ｔｏｋｙｏ，日
本，ＩＥＥＥ Ｔｏｋｙｏ，1995年 ２
月，Ｎｏ．33，ｐｐ．244−249
(58)調査した分野(Int.Cl.7，ＤＢ名)
H04B 1/76
H04B 3/00
H04B 7/005 - 7/015
ＪＩＣＳＴファイル（ＪＯＩＳ)
ＩＮＳＰＥＣ（ＤＩＡＬＯＧ)
Claims (3)
- 【請求項１】 発音手段からの音が近接して配置された
収音手段に廻り込むことにより生じるエコーを除去する
エコー除去装置において、 上記発音手段に供給される発音信号に基づいて上記収音
手段に廻り込むエコー成分を推定した擬似エコー信号を
出力するフィルタ手段と、 上記収音手段からの収音信号より上記フィルタ手段から
の擬似エコー信号を減算して出力する減算手段と、 上記発音信号を平滑化した平滑値を用いて上記フィルタ
手段のフィルタ係数を推定して更新するフィルタ係数推
定手段とを有し、 上記フィルタ係数推定手段は、上記フィルタ手段のタッ
プ数Ｌより大きな数Ｎについて過去Ｎ個のサンプルの上
記発音信号の電力のＬ個分に相当する電力を上記平滑値
として用いる ことを特徴とするエコー除去装置。 - 【請求項２】 上記フィルタ係数推定手段は、学習同定
法によりフィルタのタップ係数の推定を行い、この学習
同定法によるタップ係数推定に用いられる式中の２乗和
の値として、上記発音信号の過去Ｎ個のサンプルの２乗
和の内のＬ個分相当の２乗和の値を上記平滑値として用
いることを特徴とする請求項１記載のエコー除去装置。 - 【請求項３】 上記フィルタ係数推定手段は、上記発音
信号の現在のサンプルをｘ(k) とするとき、現在の平滑
値であるＰ ｘ (k)を、 【数１】 により求める ことを特徴とする請求項２記載のエコー除
去装置。
Priority Applications (8)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP05002595A JP3381112B2 (ja) | 1995-03-09 | 1995-03-09 | エコー除去装置 |
| EP96301443A EP0731593B1 (en) | 1995-03-09 | 1996-03-04 | Echo removing apparatus |
| DE1996631708 DE69631708T2 (de) | 1995-03-09 | 1996-03-04 | Einrichtung zur Echounterdrückung |
| US08/611,199 US5737409A (en) | 1995-03-09 | 1996-03-05 | Echo removing apparatus |
| TW085102752A TW286467B (ja) | 1995-03-09 | 1996-03-06 | |
| MYPI96000831A MY117594A (en) | 1995-03-09 | 1996-03-06 | Echo removing apparatus |
| KR1019960006086A KR100415525B1 (ko) | 1995-03-09 | 1996-03-08 | 에코제거장치 |
| CN96106065A CN1118179C (zh) | 1995-03-09 | 1996-03-09 | 回声排除装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP05002595A JP3381112B2 (ja) | 1995-03-09 | 1995-03-09 | エコー除去装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH08251081A JPH08251081A (ja) | 1996-09-27 |
| JP3381112B2 true JP3381112B2 (ja) | 2003-02-24 |
Family
ID=12847468
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP05002595A Expired - Fee Related JP3381112B2 (ja) | 1995-03-09 | 1995-03-09 | エコー除去装置 |
Country Status (8)
| Country | Link |
|---|---|
| US (1) | US5737409A (ja) |
| EP (1) | EP0731593B1 (ja) |
| JP (1) | JP3381112B2 (ja) |
| KR (1) | KR100415525B1 (ja) |
| CN (1) | CN1118179C (ja) |
| DE (1) | DE69631708T2 (ja) |
| MY (1) | MY117594A (ja) |
| TW (1) | TW286467B (ja) |
Families Citing this family (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6256394B1 (en) * | 1996-01-23 | 2001-07-03 | U.S. Philips Corporation | Transmission system for correlated signals |
| CA2249283C (en) * | 1996-04-03 | 2003-06-03 | British Telecommunications Public Limited Company | Acoustic feedback correction |
| JPH10257583A (ja) * | 1997-03-06 | 1998-09-25 | Asahi Chem Ind Co Ltd | 音声処理装置およびその音声処理方法 |
| DE19831320A1 (de) * | 1998-07-13 | 2000-01-27 | Ericsson Telefon Ab L M | Digitales adaptives Filter und akustischer Echokompensator |
| JP2000252881A (ja) * | 1999-02-25 | 2000-09-14 | Mitsubishi Electric Corp | ダブルトーク検知装置並びにエコーキャンセラ装置およびエコーサプレッサー装置 |
| GB2350969B (en) | 1999-06-07 | 2003-11-05 | Ericsson Telefon Ab L M | Loudspeaker volume range control |
| KR100399902B1 (ko) * | 2001-05-02 | 2003-09-29 | 주식회사 하이닉스반도체 | 고속 유선 통신용 액티브 하이브리드 회로 |
| US7778408B2 (en) * | 2004-12-30 | 2010-08-17 | Texas Instruments Incorporated | Method and apparatus for acoustic echo cancellation utilizing dual filters |
| JP4765117B2 (ja) | 2006-03-02 | 2011-09-07 | オンセミコンダクター・トレーディング・リミテッド | エコー防止回路、フィルタ係数設定方法、及びプログラム |
| CN101453244A (zh) * | 2007-11-28 | 2009-06-10 | 博通集成电路（上海）有限公司 | 线回波消除器 |
| CN106534600A (zh) * | 2016-11-24 | 2017-03-22 | 浪潮（苏州）金融技术服务有限公司 | 一种回声消除装置、方法及系统 |
Family Cites Families (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR880001981B1 (ko) * | 1983-01-19 | 1988-10-08 | 에스티씨 피엘씨 | 반향제거기 |
| US4591669A (en) * | 1984-09-26 | 1986-05-27 | At&T Bell Laboratories | Adaptive filter update gain normalization |
| DE3840433A1 (de) * | 1988-12-01 | 1990-06-07 | Philips Patentverwaltung | Echokompensator |
| DE69131443T2 (de) * | 1990-11-01 | 1999-11-18 | Nec Corp., Tokio/Tokyo | Interferenzunterdrücker mit Steuerung der Gewichtungsfaktorenanpassung mit zum Signalleistungspegel umgekehrt proportionaler Stufenhöhe |
| KR950005942B1 (ko) * | 1992-05-28 | 1995-06-07 | 삼성전자주식회사 | 디지탈 적응필터의 동작영역 판단장치 |
| JP2654894B2 (ja) * | 1992-09-30 | 1997-09-17 | 日本電信電話株式会社 | 反響消去装置およびその方法 |
| JP2538176B2 (ja) * | 1993-05-28 | 1996-09-25 | 松下電器産業株式会社 | エコ―制御装置 |
| US5477534A (en) * | 1993-07-30 | 1995-12-19 | Kyocera Corporation | Acoustic echo canceller |
| US5572621A (en) * | 1993-09-21 | 1996-11-05 | U.S. Philips Corporation | Speech signal processing device with continuous monitoring of signal-to-noise ratio |
| JP2947093B2 (ja) * | 1994-11-02 | 1999-09-13 | 日本電気株式会社 | 適応フィルタによるシステム同定の方法および装置 |
-
1995
- 1995-03-09 JP JP05002595A patent/JP3381112B2/ja not_active Expired - Fee Related
-
1996
- 1996-03-04 EP EP96301443A patent/EP0731593B1/en not_active Expired - Lifetime
- 1996-03-04 DE DE1996631708 patent/DE69631708T2/de not_active Expired - Lifetime
- 1996-03-05 US US08/611,199 patent/US5737409A/en not_active Expired - Lifetime
- 1996-03-06 TW TW085102752A patent/TW286467B/zh not_active IP Right Cessation
- 1996-03-06 MY MYPI96000831A patent/MY117594A/en unknown
- 1996-03-08 KR KR1019960006086A patent/KR100415525B1/ko not_active IP Right Cessation
- 1996-03-09 CN CN96106065A patent/CN1118179C/zh not_active Expired - Lifetime
Non-Patent Citations (2)
| Title |
|---|
| Ｋｅｎｓａｋｕ Ｆｕｊｉｉ ａｎｄ Ｊｕｒｏ Ｏｈｇａ，"Ａ ＦＡＳＴ ＡＤＡＰＴＩＶＥ ＡＬＧＯＲＩＴＨＭ ＳＵＩＴＡＢＬＥ ＦＯＲ ＡＣＯＵＳＴＩＣ ＥＣＨＯ ＣＡＮＣＥＬＬＥＲ"，Ｄｅｎｓｈｉ Ｔｏｋｙｏ，日本，ＩＥＥＥ Ｔｏｋｙｏ，１９９５年 ２月，Ｎｏ．３３，ｐｐ．２４４−２４９ |
| Ｋｅｎｓａｋｕ Ｆｕｊｉｉ ａｎｄ Ｊｕｒｏ Ｏｈｇａ，"Ａ Ｆａｓｔ Ａｄａｐｔｉｖｅ Ａｌｇｏｒｉｔｈｍ Ｓｕｉｔａｂｌｅ ｆｏｒ Ａｃｏｕｓｔｉｃ Ｅｃｈｏ Ｃａｎｃｅｌｌｅｒ"，ＩＥＩＣＥ Ｔｒａｎｓａｃｔｉｏｎｓ ｏｎ Ｆｕｎｄａｍｅｎｔａｌｓ ｏｆ Ｅｌｅｃｔｒｏｎｉｃｓ，Ｃｏｍｍｕｎｉｃａｔｉｏｎｓ ａｎｄ Ｃｏｍｐｕｔｅｒ Ｓｃｉｅｎｃｅｓ，日本，１９９２年，Ｖｏｌ．Ｅ７５−Ａ，Ｎｏ．１１，ｐｐ．１５０９−１５１５ |
Also Published As
| Publication number | Publication date |
|---|---|
| EP0731593B1 (en) | 2004-03-03 |
| US5737409A (en) | 1998-04-07 |
| DE69631708D1 (de) | 2004-04-08 |
| JPH08251081A (ja) | 1996-09-27 |
| CN1118179C (zh) | 2003-08-13 |
| EP0731593A2 (en) | 1996-09-11 |
| TW286467B (ja) | 1996-09-21 |
| EP0731593A3 (en) | 2001-08-22 |
| KR100415525B1 (ko) | 2004-04-28 |
| MY117594A (en) | 2004-07-31 |
| KR960036375A (ko) | 1996-10-28 |
| CN1157522A (zh) | 1997-08-20 |
| DE69631708T2 (de) | 2005-01-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| KR100974370B1 (ko) | 에코 억압 방법 및 장치 | |
| US6917688B2 (en) | Adaptive noise cancelling microphone system | |
| US6738480B1 (en) | Method and device for cancelling stereophonic echo with frequency domain filtering | |
| US20110261950A1 (en) | Adaptive filter and echo canceller having this adaptive filter | |
| WO1999003091A1 (en) | Methods and apparatus for measuring signal level and delay at multiple sensors | |
| JP4101317B2 (ja) | 適応識別とそれに関する適応エコーキャンセラのための方法と装置 | |
| CN107134281A (zh) | 一种自适应回声消除中自适应滤波器系数更新方法 | |
| JP3381112B2 (ja) | エコー除去装置 | |
| JP2000502861A (ja) | 収束測定装置および収束測定方法 | |
| CN112863532A (zh) | 回音抑制装置、回音抑制方法以及存储介质 | |
| EP0732838B1 (en) | Acoustic echo cancellor | |
| CN107071196B (zh) | 一种自适应回声消除方法 | |
| US6694020B1 (en) | Frequency domain stereophonic acoustic echo canceller utilizing non-linear transformations | |
| JP4396449B2 (ja) | 残響除去方法及びその装置 | |
| US8036376B2 (en) | Echo prevention circuit having adaptive digital filter | |
| US8934621B2 (en) | Echo prevention circuit having signal subtracter feature | |
| JP2002223182A (ja) | 反響消去方法、その装置、そのプログラム及びその記録媒体 | |
| KR100272131B1 (ko) | 계층적 구조의 적응반향 제거장치 | |
| JPH09261135A (ja) | 音響エコー消去装置 | |
| JP3616341B2 (ja) | 多チャネルエコーキャンセル方法、その装置、そのプログラム及び記録媒体 | |
| JP4903843B2 (ja) | 適応フィルタ及びこれを有するエコーキャンセラ | |
| JP2000252891A (ja) | 信号処理装置 | |
| JP4964267B2 (ja) | 適応フィルタ及びこれを有するエコーキャンセラ | |
| JP3403655B2 (ja) | サブバンド適応フィルタを用いた未知システムの同定方法および装置 | |
| JP2001148644A (ja) | エコーキャンセラ |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20021112 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20071220 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20081220 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20091220 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20091220 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20101220 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20101220 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20111220 Year of fee payment: 9 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20111220 Year of fee payment: 9 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20121220 Year of fee payment: 10 |
|
| LAPS | Cancellation because of no payment of annual fees |