JP3145269B2 - エコーキャンセラ制御方法 - Google Patents
エコーキャンセラ制御方法Info
- Publication number
- JP3145269B2 JP3145269B2 JP06000595A JP6000595A JP3145269B2 JP 3145269 B2 JP3145269 B2 JP 3145269B2 JP 06000595 A JP06000595 A JP 06000595A JP 6000595 A JP6000595 A JP 6000595A JP 3145269 B2 JP3145269 B2 JP 3145269B2
- Authority
- JP
- Japan
- Prior art keywords
- echo
- value
- equation
- erle
- step gain
- 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
Landscapes
- Filters That Use Time-Delay Elements (AREA)
- Cable Transmission Systems, Equalization Of Radio And Reduction Of Echo (AREA)
Description
【０００１】
【産業上の利用分野】本発明は、電話回線や電子会議シ
ステム等で発生するエコーの排除機能を持つエコーキャ
ンセラの制御方法に関し、特に、エコーを迅速且つ高精
度に除去できるようにエコーキャンセラの特性を制御す
るものである。
ステム等で発生するエコーの排除機能を持つエコーキャ
ンセラの制御方法に関し、特に、エコーを迅速且つ高精
度に除去できるようにエコーキャンセラの特性を制御す
るものである。
【０００２】
【従来の技術】電話回線の２線式回線と４線式回線との
接続部では、インピーダンス不整合のために受信信号の
一部が送信側に漏れる、いわゆるエコーが発生する。同
様のエコーは電子会議システム等でも発生する。
接続部では、インピーダンス不整合のために受信信号の
一部が送信側に漏れる、いわゆるエコーが発生する。同
様のエコーは電子会議システム等でも発生する。
【０００３】電話回線の場合では、図５に示すように、
ディジタル受信信号ｘj はＰＣＭコーデック２でアナロ
グ信号に変換された後、２線／４線変換回路８を経て電
話器９から音声出力され、一方、電話器９から入力され
たアナログ音声信号は、２線／４線変換回路８を経てＰ
ＣＭコーデック２に入力し、ディジタル信号に変換され
て送信出力端10から相手方に送信される。このとき、Ｐ
ＣＭコーデック２でアナログ変換された受信信号の一部
は、２線／４線変換回路８から送信側に漏れ、ＰＣＭコ
ーデック２に再入力し、送信信号に混ざって相手方に送
り返される。
ディジタル受信信号ｘj はＰＣＭコーデック２でアナロ
グ信号に変換された後、２線／４線変換回路８を経て電
話器９から音声出力され、一方、電話器９から入力され
たアナログ音声信号は、２線／４線変換回路８を経てＰ
ＣＭコーデック２に入力し、ディジタル信号に変換され
て送信出力端10から相手方に送信される。このとき、Ｐ
ＣＭコーデック２でアナログ変換された受信信号の一部
は、２線／４線変換回路８から送信側に漏れ、ＰＣＭコ
ーデック２に再入力し、送信信号に混ざって相手方に送
り返される。
【０００４】これがエコーであり、このエコーを抑圧す
るために、エコーキャンセラ３が受信入力端１と送信出
力端10との間に配置される。従来のエコーキャンセラ３
は、受信入力端１から入力した受信入力ｘj を記憶する
シフトレジスタ４と、ＰＣＭコーデック２のアナログ出
力からアナログ入力に至るエコー経路のインパルス応答
の推定値Ｈj を記憶するメモリ５と、受信入力ｘj とエ
コーの推定値Ｈj とを畳込み演算して擬似エコーｙ'jを
発生する畳込み演算器６と、送信入力から擬似エコーを
差引いた値ｅj を送信出力端８に出力する減算器７とを
備えている。
るために、エコーキャンセラ３が受信入力端１と送信出
力端10との間に配置される。従来のエコーキャンセラ３
は、受信入力端１から入力した受信入力ｘj を記憶する
シフトレジスタ４と、ＰＣＭコーデック２のアナログ出
力からアナログ入力に至るエコー経路のインパルス応答
の推定値Ｈj を記憶するメモリ５と、受信入力ｘj とエ
コーの推定値Ｈj とを畳込み演算して擬似エコーｙ'jを
発生する畳込み演算器６と、送信入力から擬似エコーを
差引いた値ｅj を送信出力端８に出力する減算器７とを
備えている。
【０００５】このエコーキャンセラ３は、受信入力端１
から受信入力ｘj が入力すると、それをシフトレジスタ
４にベクトルＸj として取り込む。受信入力ｘj は、同
時にＰＣＭコーデック２に入力し、ＰＣＭコーデック２
からアナログ信号に変換されて出力され、その一部がエ
コーｙ（ｔ）となって再びＰＣＭコーデック２に入力す
る。
から受信入力ｘj が入力すると、それをシフトレジスタ
４にベクトルＸj として取り込む。受信入力ｘj は、同
時にＰＣＭコーデック２に入力し、ＰＣＭコーデック２
からアナログ信号に変換されて出力され、その一部がエ
コーｙ（ｔ）となって再びＰＣＭコーデック２に入力す
る。
【０００６】ＰＣＭコーデック２には、その他に周囲騒
音等をも含む送信信号ｎ（ｔ）が入力する。従って、Ｐ
ＣＭコーデック２のディジタル出力信号はｓj＝ｙj＋ｎ
j となり、送信信号ｎj は、エコー信号ｙj に対する外
乱として影響する。
音等をも含む送信信号ｎ（ｔ）が入力する。従って、Ｐ
ＣＭコーデック２のディジタル出力信号はｓj＝ｙj＋ｎ
j となり、送信信号ｎj は、エコー信号ｙj に対する外
乱として影響する。
【０００７】一方、エコーキャンセラ３内部のメモリ５
には、ＰＣＭコーデック２のアナログ出力端からアナロ
グ入力端に至るエコー経路の推定インパルス応答がベク
トルＨj として記憶されており、畳込み演算器６は、こ
のＨj と受信入力列Ｘj との畳込み演算を行なって擬似
エコーｙ'jを発生する。
には、ＰＣＭコーデック２のアナログ出力端からアナロ
グ入力端に至るエコー経路の推定インパルス応答がベク
トルＨj として記憶されており、畳込み演算器６は、こ
のＨj と受信入力列Ｘj との畳込み演算を行なって擬似
エコーｙ'jを発生する。
【０００８】減算器７は、送信入力ｓj から擬似エコー
ｙ'jを差引いた残留エコーｅj ＝ｓj −ｙ'jを算出し、
これを送信出力端10から送信する。その結果、エコーの
戻りが抑制される。
ｙ'jを差引いた残留エコーｅj ＝ｓj −ｙ'jを算出し、
これを送信出力端10から送信する。その結果、エコーの
戻りが抑制される。
【０００９】次いで、エコーキャンセラ３は、エコー抑
圧精度を高めるために、式（１）（２）に従って推定値
Ｈj を補正する。このアルゴリズムは学習同定法と呼ば
れ、エコーキャンセラでは最もよく用いられている。
圧精度を高めるために、式（１）（２）に従って推定値
Ｈj を補正する。このアルゴリズムは学習同定法と呼ば
れ、エコーキャンセラでは最もよく用いられている。
【００１０】 エコー消去演算； ｅj ＝ｓj−Ｘj T・Ｈj （１） 推定値補正演算； Ｈj+1＝Ｈj＋α（ｅj／‖Ｘj‖2）Ｘj （２） ここで、Ｘj＝（xj，xj-1,‥,xj-N+1）T ：受信入力ベ
クトル Ｈj＝（h1 j，h2 j, ‥,hN j ）T ：推定インパルス応答
ベクトル であり、 Ｘj T・Ｈj＝Σxj-m・hm j （畳込み演算、Σはm=0からN-
1まで加算）、 ‖Ｘj‖2＝Σxj-m 2 （Σはm=0からN-1まで加算） である。また、α（０＜α＜２）は、ステップゲインと
呼ばれる補正演算の補正幅を与えるパラメータであり、
従来のエコーキャンセラでは外乱に対する推定精度の劣
化を考慮してα＝0.125としている。
クトル Ｈj＝（h1 j，h2 j, ‥,hN j ）T ：推定インパルス応答
ベクトル であり、 Ｘj T・Ｈj＝Σxj-m・hm j （畳込み演算、Σはm=0からN-
1まで加算）、 ‖Ｘj‖2＝Σxj-m 2 （Σはm=0からN-1まで加算） である。また、α（０＜α＜２）は、ステップゲインと
呼ばれる補正演算の補正幅を与えるパラメータであり、
従来のエコーキャンセラでは外乱に対する推定精度の劣
化を考慮してα＝0.125としている。
【００１１】この式（１）によるエコーの消去と、式
（２）による推定値の補正とを交互に繰返すことによ
り、外乱ｎj のレベルがエコーｙj のレベルより低い期
間では、推定値Ｈj が真のインパルス応答に漸次収束
し、エコー抑圧効果が増加する。
（２）による推定値の補正とを交互に繰返すことによ
り、外乱ｎj のレベルがエコーｙj のレベルより低い期
間では、推定値Ｈj が真のインパルス応答に漸次収束
し、エコー抑圧効果が増加する。
【００１２】エコーキャンセラの収束特性は、受信入力
ｘj 及び外乱ｎj を定常白色ガウス雑音とした場合に、
式（３）のように理論的に求められる。 <ｅ0j+1 2>＝｛１−α（２−α）／Ｎ｝<ｅ0j 2>＋（α2／Ｎ）<ｎj 2> （３） ここで、ｅ0jは時刻ｊにおけるＨj の推定誤差による残
留エコー成分であり、真のインパルス応答をＨとすると
ｅ0j＝Ｘj T・（Ｈ−Ｈj）である。また、<ｅ0j 2>、<ｎj
2 >は各々の信号の２乗の集合平均値を表している。Ｎ
は推定インパルス応答長（タップ数）である。
ｘj 及び外乱ｎj を定常白色ガウス雑音とした場合に、
式（３）のように理論的に求められる。 <ｅ0j+1 2>＝｛１−α（２−α）／Ｎ｝<ｅ0j 2>＋（α2／Ｎ）<ｎj 2> （３） ここで、ｅ0jは時刻ｊにおけるＨj の推定誤差による残
留エコー成分であり、真のインパルス応答をＨとすると
ｅ0j＝Ｘj T・（Ｈ−Ｈj）である。また、<ｅ0j 2>、<ｎj
2 >は各々の信号の２乗の集合平均値を表している。Ｎ
は推定インパルス応答長（タップ数）である。
【００１３】エコーキャンセラは、外乱<ｎj 2>の存在下
で残留エコーの２乗値<ｅ0j 2 >をより速く減少させるよ
うに動作することが望まれる。式（３）より、<ｅ0j 2 >
は１サンプル当たり平均｛１−α（２−α）／Ｎ｝倍で
減少し、最終的には<ｅ0∞ 2>＝｛α／（２−α）｝<ｎ
2 >に収束することが分かる。ただし、<ｎj 2>＝<ｎj+1 2
>＝<ｎ2 >＝一定とする。
で残留エコーの２乗値<ｅ0j 2 >をより速く減少させるよ
うに動作することが望まれる。式（３）より、<ｅ0j 2 >
は１サンプル当たり平均｛１−α（２−α）／Ｎ｝倍で
減少し、最終的には<ｅ0∞ 2>＝｛α／（２−α）｝<ｎ
2 >に収束することが分かる。ただし、<ｎj 2>＝<ｎj+1 2
>＝<ｎ2 >＝一定とする。
【００１４】このエコーキャンセラの性能は打消量ＥＲ
ＬＥ＝（<ｙj 2>／<ｅ0j 2 >）で表される。式（３）の両
辺を<ｙ2 >で割ることにより、式（４）が得られる。
ＬＥ＝（<ｙj 2>／<ｅ0j 2 >）で表される。式（３）の両
辺を<ｙ2 >で割ることにより、式（４）が得られる。
【００１５】 1／ERLEj+1＝｛１−α（２−α）／Ｎ｝(1／ERLEj)＋(α2／Ｎ)(１／SNR） （４） ここで、ＳＮＲ（＝<ｎj 2>／<ｙj 2>）はエコー信号とエ
コー経路で加わる雑音及び非線形歪（外乱）とのレベル
比であり、ここではＳＮＲ一定としている。
コー経路で加わる雑音及び非線形歪（外乱）とのレベル
比であり、ここではＳＮＲ一定としている。
【００１６】図６には、この式（４）によって打消量Ｅ
ＲＬＥが時間とともに増加する従来のエコーキャンセラ
の収束特性を示している。（１）は回線用エコーキャン
セラ（ＳＮＲ＝１００００（４０ｄＢ））の特性であ
り、（２）は音響用エコーキャンセラ（ＳＮＲ＝３２
（１５ｄＢ））の特性である。
ＲＬＥが時間とともに増加する従来のエコーキャンセラ
の収束特性を示している。（１）は回線用エコーキャン
セラ（ＳＮＲ＝１００００（４０ｄＢ））の特性であ
り、（２）は音響用エコーキャンセラ（ＳＮＲ＝３２
（１５ｄＢ））の特性である。
【００１７】
【発明が解決しようとする課題】このエコーキャンセラ
の特性は、ステップゲインαを大きく取ると、ＥＲＬＥ
の立上りは速くなるが定常状態における誤差が増大し、
αを小さく取ると、逆の傾向が現れる。従来のエコーキ
ャンセラでは、収束速度と推定精度との兼ね合いによ
り、αを固定値に設定しているため、収束開始時に求め
られる大きな収束速度と、定常時に求められる高い推定
精度とを共に満足する収束特性を持つことができなかっ
た。
の特性は、ステップゲインαを大きく取ると、ＥＲＬＥ
の立上りは速くなるが定常状態における誤差が増大し、
αを小さく取ると、逆の傾向が現れる。従来のエコーキ
ャンセラでは、収束速度と推定精度との兼ね合いによ
り、αを固定値に設定しているため、収束開始時に求め
られる大きな収束速度と、定常時に求められる高い推定
精度とを共に満足する収束特性を持つことができなかっ
た。
【００１８】本発明は、こうした従来の問題点を解決す
るものであり、収束開始時には大きな収束速度で収束
し、定常時には高い推定精度を維持することができるよ
うにエコーキャンセラを制御する方法を提供することを
目的としている。
るものであり、収束開始時には大きな収束速度で収束
し、定常時には高い推定精度を維持することができるよ
うにエコーキャンセラを制御する方法を提供することを
目的としている。
【００１９】
【課題を解決するための手段】そこで、本発明では、受
信入力列とエコー経路のインパルス応答の推定値とを畳
込み演算して疑似エコーを発生し、これをエコー信号か
ら差し引くことによってエコーを抑圧し、このときに得
られる残留エコーの２乗値を零に近づけるように前記推
定値を逐次補正するエコーキャンセラの制御方法におい
て、この推定値の逐次補正におけるステップゲインαj
を式 αj+1＝（Ｎ−αj）αj／（Ｎ−αj 2） （但し、Ｎは推定インパルス応答長（タップ数））によ
って設定している。
信入力列とエコー経路のインパルス応答の推定値とを畳
込み演算して疑似エコーを発生し、これをエコー信号か
ら差し引くことによってエコーを抑圧し、このときに得
られる残留エコーの２乗値を零に近づけるように前記推
定値を逐次補正するエコーキャンセラの制御方法におい
て、この推定値の逐次補正におけるステップゲインαj
を式 αj+1＝（Ｎ−αj）αj／（Ｎ−αj 2） （但し、Ｎは推定インパルス応答長（タップ数））によ
って設定している。
【００２０】また、この推定値の逐次補正におけるステ
ップゲインαj を式 αj+1＝｛１−αj（１−αj）／Ｎ｝αj によって設定している。
ップゲインαj を式 αj+1＝｛１−αj（１−αj）／Ｎ｝αj によって設定している。
【００２１】また、この推定値の逐次補正におけるステ
ップゲインαj を式 αj＝ERLEj -1／(ERLEj -1＋SNR-1) （但し、ERLEjは時刻ｊにおける打消量、SNRはエコー信
号とエコー信号経路で加わる外乱とのレベル比）によっ
て設定している。
ップゲインαj を式 αj＝ERLEj -1／(ERLEj -1＋SNR-1) （但し、ERLEjは時刻ｊにおける打消量、SNRはエコー信
号とエコー信号経路で加わる外乱とのレベル比）によっ
て設定している。
【００２２】また、ステップゲインαj の値が式 αj＝ERLEj -1／(ERLEj -1＋SNR-1) の関係から大きくずれたとき、ステップゲインαj の値
をこの式によって再設定し、その後、ステップゲインα
j の値を式 αj+1＝（Ｎ−αj）αj／（Ｎ−αj 2） によって逐次設定している。
をこの式によって再設定し、その後、ステップゲインα
j の値を式 αj+1＝（Ｎ−αj）αj／（Ｎ−αj 2） によって逐次設定している。
【００２３】また、ステップゲインαj の値が式 αj＝ERLEj -1／(ERLEj -1＋SNR-1) の関係から大きくずれたとき、ステップゲインαj の値
をこの式によって再設定し、その後、ステップゲインα
j の値を式 αj+1＝｛１−αj（１−αj）／Ｎ｝αj によって逐次設定している。
をこの式によって再設定し、その後、ステップゲインα
j の値を式 αj+1＝｛１−αj（１−αj）／Ｎ｝αj によって逐次設定している。
【００２４】
【作用】式（４）から、打消量ＥＲＬＥを最大にする最
適ステップゲインが αj ＝ERLEj -1／(ERLEj -1＋SNR-1) （５） として求められ、この式を変形整理することにより αj+1＝（Ｎ−αj）αj／（Ｎ−αj 2） （６） が得られる。この式（６）に従って推定値の逐次補正に
おけるステップゲインαj を制御することにより、αj
は、最適値を取りながら収束の進行に伴って減少する。
その結果、エコーキャンセラは収束開始時の収束速度が
大きく、定常時の推定精度が高い収束特性を得ることが
できる。
適ステップゲインが αj ＝ERLEj -1／(ERLEj -1＋SNR-1) （５） として求められ、この式を変形整理することにより αj+1＝（Ｎ−αj）αj／（Ｎ−αj 2） （６） が得られる。この式（６）に従って推定値の逐次補正に
おけるステップゲインαj を制御することにより、αj
は、最適値を取りながら収束の進行に伴って減少する。
その結果、エコーキャンセラは収束開始時の収束速度が
大きく、定常時の推定精度が高い収束特性を得ることが
できる。
【００２５】また、式（６）では、分母に変数αj を含
むため、割り算が必須となる。式（６）を近似して得た αj+1＝｛１−αj（１−αj）／Ｎ｝αj （７） は、１／Ｎが定数であるため、演算に当たって割り算が
不要になり、ステップゲインαj の逐次計算が容易にな
る。
むため、割り算が必須となる。式（６）を近似して得た αj+1＝｛１−αj（１−αj）／Ｎ｝αj （７） は、１／Ｎが定数であるため、演算に当たって割り算が
不要になり、ステップゲインαj の逐次計算が容易にな
る。
【００２６】また、打消量ＥＲＬＥを適当な時間間隔で
計算し、それを基に式（５）からステップゲインαj を
求めることにより、最適ステップゲインαj を収束の進
行状況に合わせて減少させることもできる。
計算し、それを基に式（５）からステップゲインαj を
求めることにより、最適ステップゲインαj を収束の進
行状況に合わせて減少させることもできる。
【００２７】また、式（５）と式（６）、または式
（５）と式（７）とを用いるエコーキャンセラ制御方法
では、収束開始時や通話回線の切替えでエコーパスが変
化したときなど、推定値のずれの大きい段階では、その
時点での打消量ＥＲＬＥを求めて式（５）によりαj を
制御し、その後の収束進行時には式（６）または式
（７）による逐次制御を行なう。
（５）と式（７）とを用いるエコーキャンセラ制御方法
では、収束開始時や通話回線の切替えでエコーパスが変
化したときなど、推定値のずれの大きい段階では、その
時点での打消量ＥＲＬＥを求めて式（５）によりαj を
制御し、その後の収束進行時には式（６）または式
（７）による逐次制御を行なう。
【００２８】
（第１実施例）本発明の制御方法を実施するエコーキャ
ンセラは、機能ブロックとして、図１に示すように、推
定値Ｈj の補正におけるステップゲインαを制御するα
制御部11を備えている。その他の構成は従来の装置（図
５）と変わりがない。
ンセラは、機能ブロックとして、図１に示すように、推
定値Ｈj の補正におけるステップゲインαを制御するα
制御部11を備えている。その他の構成は従来の装置（図
５）と変わりがない。
【００２９】このエコーキャンセラ３では、受信入力端
１から受信入力ｘj が入力すると、シフトレジスタ４が
ベクトルＸj として取り込み、畳込み演算器６が、メモ
リ５に記憶しているエコー経路の推定インパルス応答ベ
クトルＨj と受信入力列Ｘjとの畳込み演算を行なって
擬似エコーｙ'jを発生し、減算器７が、送信入力ｓjか
ら擬似エコーｙ'jを差引いた残留エコーｅj ＝ｓj −
ｙ'jを送信出力端10に出力する。
１から受信入力ｘj が入力すると、シフトレジスタ４が
ベクトルＸj として取り込み、畳込み演算器６が、メモ
リ５に記憶しているエコー経路の推定インパルス応答ベ
クトルＨj と受信入力列Ｘjとの畳込み演算を行なって
擬似エコーｙ'jを発生し、減算器７が、送信入力ｓjか
ら擬似エコーｙ'jを差引いた残留エコーｅj ＝ｓj −
ｙ'jを送信出力端10に出力する。
【００３０】そして、送信出力ｓj のレベルを計算し、
残留エコーｅj のレベルを計算し、α制御部11で、後述
する式（６）によりステップゲインαj の値を計算し、
それらの結果を用いて次式（１）（２’）によってメモ
リ５に記憶した推定値Ｈj を補正する。エコーキャンセ
ラはこの手順を順次繰返す。
残留エコーｅj のレベルを計算し、α制御部11で、後述
する式（６）によりステップゲインαj の値を計算し、
それらの結果を用いて次式（１）（２’）によってメモ
リ５に記憶した推定値Ｈj を補正する。エコーキャンセ
ラはこの手順を順次繰返す。
【００３１】推定値Ｈj の補正に用いる式（１）
（２’）は、 エコー消去演算； ｅj ＝ｓj−Ｘj T・Ｈj （１） 推定値補正演算； Ｈj+1＝Ｈj＋αj（ｅj／‖Ｘj‖2）Ｘj （２'） ここで、Ｘj＝（xj，xj-1,‥,xj-N+1）T ：受信入力ベ
クトル Ｈj＝（h1 j，h2 j,‥,hN j ）T ：推定インパルス応答
ベクトル であり、 Ｘj T・Ｈj＝Σxj-m・hm j （畳込み演算、Σはm=0からN-
1まで加算）、 ‖Ｘj‖2＝Σxj-m 2 （Σはm=0からN-1まで加算） である。この式は、推定値補正演算（２’）のステップ
ゲインαが、固定値では無く、時刻ｊによって変化する
値αj である点で、従来のエコーキャンセラの用いてい
る式と異なっている。
（２’）は、 エコー消去演算； ｅj ＝ｓj−Ｘj T・Ｈj （１） 推定値補正演算； Ｈj+1＝Ｈj＋αj（ｅj／‖Ｘj‖2）Ｘj （２'） ここで、Ｘj＝（xj，xj-1,‥,xj-N+1）T ：受信入力ベ
クトル Ｈj＝（h1 j，h2 j,‥,hN j ）T ：推定インパルス応答
ベクトル であり、 Ｘj T・Ｈj＝Σxj-m・hm j （畳込み演算、Σはm=0からN-
1まで加算）、 ‖Ｘj‖2＝Σxj-m 2 （Σはm=0からN-1まで加算） である。この式は、推定値補正演算（２’）のステップ
ゲインαが、固定値では無く、時刻ｊによって変化する
値αj である点で、従来のエコーキャンセラの用いてい
る式と異なっている。
【００３２】第１実施例の制御方法では、このαj を、
初期値をα0 ＝１として、式（６）により逐次求める。
初期値をα0 ＝１として、式（６）により逐次求める。
【００３３】 αj+1＝（Ｎ−αj）αj／（Ｎ−αj 2） （６）
【００３４】この式を導出した根拠について説明する。
エコーキャンセラの収束特性を表す式（４）は、次のよ
うに変形することができる。
エコーキャンセラの収束特性を表す式（４）は、次のよ
うに変形することができる。
【００３５】１／ERLEj+1＝{(ERLEj -1＋SNR-1)／Ｎ}[α
−{ERLEj -1／(ERLEj -1＋SNR-1)}]2＋ERLEj -1−ERLEj -1
／Ｎ(ERLEj -1＋SNR-1) この式から、ＥＲＬＥj+1 を最大にする最適ステップゲ
インを次式（５）によって求めることができる。 α＝ERLEj -1／(ERLEj -1＋SNR-1) （５） この式（５）におけるαとＥＲＬＥとの関係を図２に示
している。このα−ＥＲＬＥ特性はタップ数Ｎには依存
せず、ＳＮＲが大きい程、ＥＲＬＥに対応するα値が高
い値を示す。
−{ERLEj -1／(ERLEj -1＋SNR-1)}]2＋ERLEj -1−ERLEj -1
／Ｎ(ERLEj -1＋SNR-1) この式から、ＥＲＬＥj+1 を最大にする最適ステップゲ
インを次式（５）によって求めることができる。 α＝ERLEj -1／(ERLEj -1＋SNR-1) （５） この式（５）におけるαとＥＲＬＥとの関係を図２に示
している。このα−ＥＲＬＥ特性はタップ数Ｎには依存
せず、ＳＮＲが大きい程、ＥＲＬＥに対応するα値が高
い値を示す。
【００３６】この式（５）から、αj+1＝ERLEj+1 -1／(E
RLEj+1 -1＋SNR-1)、αj＝ERLEj -1／(ERLEj -1＋SNR-1)と
置いてERLEj+1、ERLEjについて解き、式（４）に代入整
理することにより前記式（６）を得ることができる。式
（６）は外乱レベルが一定であれば、ＳＮＲに依存せず
に成り立つ式である。
RLEj+1 -1＋SNR-1)、αj＝ERLEj -1／(ERLEj -1＋SNR-1)と
置いてERLEj+1、ERLEjについて解き、式（４）に代入整
理することにより前記式（６）を得ることができる。式
（６）は外乱レベルが一定であれば、ＳＮＲに依存せず
に成り立つ式である。
【００３７】第１実施例のエコーキャンセラでは、推定
値Ｈj の補正に際して、α制御部11がαj を式（６）に
よって逐次的に求め、そのαj を用いて、前記式（１）
（２'）により推定値Ｈj の補正が実行される。ステッ
プゲインαj は常に最適値を取りながら、式（６）の関
係に従って時間とともに減少する。そのため、収束初期
の段階では大きな収束速度が有し、定常状態では高い推
定精度を保持する。
値Ｈj の補正に際して、α制御部11がαj を式（６）に
よって逐次的に求め、そのαj を用いて、前記式（１）
（２'）により推定値Ｈj の補正が実行される。ステッ
プゲインαj は常に最適値を取りながら、式（６）の関
係に従って時間とともに減少する。そのため、収束初期
の段階では大きな収束速度が有し、定常状態では高い推
定精度を保持する。
【００３８】図３には、αj を式（６）によって制御す
るときの回線用エコーキャンセラ（ＳＮＲ＝１００００
（４０ｄＢ））及び音響用エコーキャンセラ（ＳＮＲ＝
３２（１５ｄＢ））における最適αの時間変化を示して
いる。時間が経過し、収束が進むに連れてαj が小さく
なり、より高い推定精度が得られることが分かる。ま
た、タップ数Ｎが大きいほど収束速度は遅くなるが、タ
ップ数Ｎが大きいときにはαj の減少する速度も遅くな
り、長い時間に渡って大きな収束速度が得られるように
調整されることが分かる。
るときの回線用エコーキャンセラ（ＳＮＲ＝１００００
（４０ｄＢ））及び音響用エコーキャンセラ（ＳＮＲ＝
３２（１５ｄＢ））における最適αの時間変化を示して
いる。時間が経過し、収束が進むに連れてαj が小さく
なり、より高い推定精度が得られることが分かる。ま
た、タップ数Ｎが大きいほど収束速度は遅くなるが、タ
ップ数Ｎが大きいときにはαj の減少する速度も遅くな
り、長い時間に渡って大きな収束速度が得られるように
調整されることが分かる。
【００３９】図４は式（６）のαj を式（４）に代入
し、逐次計算させて描いた収束特性である。従来例の図
６と比べると、収束開始時は約６倍の速さ（α＝０．７
６固定の時とほぼ同等）で収束し、定常状態でもＥＲＬ
Ｅは単調増加を続け、全体として滑らかな収束特性が得
られる。
し、逐次計算させて描いた収束特性である。従来例の図
６と比べると、収束開始時は約６倍の速さ（α＝０．７
６固定の時とほぼ同等）で収束し、定常状態でもＥＲＬ
Ｅは単調増加を続け、全体として滑らかな収束特性が得
られる。
【００４０】（第２実施例）前記式（６）は分母に変数
αj を含むため除算の実行が必須であるが、第２実施例
の制御方法では、この点を改善し、簡単な演算で変数α
j を求めることができるように構成している。この制御
方法を実施するエコーキャンセラの機能ブロックは第１
実施例の場合（図１）と同じであり、ただ、α制御部11
が、推定値Ｈj の補正のステップゲインαj を求めるた
めに用いる式だけを異にしている。
αj を含むため除算の実行が必須であるが、第２実施例
の制御方法では、この点を改善し、簡単な演算で変数α
j を求めることができるように構成している。この制御
方法を実施するエコーキャンセラの機能ブロックは第１
実施例の場合（図１）と同じであり、ただ、α制御部11
が、推定値Ｈj の補正のステップゲインαj を求めるた
めに用いる式だけを異にしている。
【００４１】この制御方法では、推定値補正演算
（２’）のステップゲインαj を、初期値をα0 ＝１と
して、次式（７）により逐次求める。 αj+1＝｛１−αj（１−αj）／Ｎ｝αj （７）
（２’）のステップゲインαj を、初期値をα0 ＝１と
して、次式（７）により逐次求める。 αj+1＝｛１−αj（１−αj）／Ｎ｝αj （７）
【００４２】この式を導出した根拠について説明する。
図３のｔ−α特性を観察すると、ある最適α値に到達す
るまでの時間が（１／Ｎ）に比例するように見て取れ
る。これは収束速度が（１／Ｎ）にほぼ比例することか
ら妥当な推定であると言える。そこで漸化式（６）を、
“αj の差分が（１／Ｎ）に比例する”と仮定して近似
すると、以下の展開となり、式（７）を得ることができ
る。
図３のｔ−α特性を観察すると、ある最適α値に到達す
るまでの時間が（１／Ｎ）に比例するように見て取れ
る。これは収束速度が（１／Ｎ）にほぼ比例することか
ら妥当な推定であると言える。そこで漸化式（６）を、
“αj の差分が（１／Ｎ）に比例する”と仮定して近似
すると、以下の展開となり、式（７）を得ることができ
る。
【００４３】 αj+1−αj＝（Ｎ−αj）αj／（Ｎ−αj 2）−αj ＝（αj 2−αj）／（Ｎ−αj 2）＊αj ≒（αj 2−αj）／Ｎ＊αj ∵１／Ｎに比例，Ｎ>>αj 2 ∴ αj+1＝｛１−αj（１−αj）／Ｎ｝αj （７）
【００４４】これは非常によい近似であり、回線用エコ
ーキャンセラ（ＳＮＲ＝１００００（４０ｄＢ））及び
音響用エコーキャンセラ（ＳＮＲ＝３２（１５ｄＢ））
における最適αの時間変化は図３と同じようになる。ま
た、式（７）のαj を式（４）に代入して逐次計算させ
て描いた収束特性も図４と同じようになる。
ーキャンセラ（ＳＮＲ＝１００００（４０ｄＢ））及び
音響用エコーキャンセラ（ＳＮＲ＝３２（１５ｄＢ））
における最適αの時間変化は図３と同じようになる。ま
た、式（７）のαj を式（４）に代入して逐次計算させ
て描いた収束特性も図４と同じようになる。
【００４５】この式（７）では（１／Ｎ）が定数なの
で、式（７）を用いてステップゲインαj を逐次的に求
めるα制御部11では除算を行なう必要がなく、演算を容
易に行なうことができる。
で、式（７）を用いてステップゲインαj を逐次的に求
めるα制御部11では除算を行なう必要がなく、演算を容
易に行なうことができる。
【００４６】（第３実施例）第３実施例の制御方法で
は、変数αj を求めるに当たって、前述した式（５）を
用いている。この制御方法を実施するエコーキャンセラ
の機能ブロックは第１実施例の場合（図１）と同じであ
り、ただ、α制御部11が、推定値Ｈj の補正のステップ
ゲインαj を求めるために用いる式だけを異にしてい
る。
は、変数αj を求めるに当たって、前述した式（５）を
用いている。この制御方法を実施するエコーキャンセラ
の機能ブロックは第１実施例の場合（図１）と同じであ
り、ただ、α制御部11が、推定値Ｈj の補正のステップ
ゲインαj を求めるために用いる式だけを異にしてい
る。
【００４７】このエコーキャンセラでは、受信入力ｘj
の入力に応じて、畳込み演算器６が擬似エコーｙ'jを発
生し、減算器７が、送信入力ｓj から擬似エコーｙ'jを
差引いた残留エコーｅj ＝ｓj −ｙ'jを送信出力端10に
出力し、その後、送信出力ｓj のレベルを計算し、残留
エコーｅj のレベルを計算し、また、α制御部11で、打
消量ＥＲＬＥを計算してから、次式 αj＝ERLEj -1／(ERLEj -1＋SNR-1) （５） を用いてαj を計算し、それらの結果を用いて前記式
（１）（２’）によってメモリ５に記憶した推定値Ｈj
を補正する。エコーキャンセラはこの手順を順次繰返
す。
の入力に応じて、畳込み演算器６が擬似エコーｙ'jを発
生し、減算器７が、送信入力ｓj から擬似エコーｙ'jを
差引いた残留エコーｅj ＝ｓj −ｙ'jを送信出力端10に
出力し、その後、送信出力ｓj のレベルを計算し、残留
エコーｅj のレベルを計算し、また、α制御部11で、打
消量ＥＲＬＥを計算してから、次式 αj＝ERLEj -1／(ERLEj -1＋SNR-1) （５） を用いてαj を計算し、それらの結果を用いて前記式
（１）（２’）によってメモリ５に記憶した推定値Ｈj
を補正する。エコーキャンセラはこの手順を順次繰返
す。
【００４８】なお、打消量ＥＲＬＥは、送信信号ｎ
（ｔ）の小さいシングルトーク状態の時にその値を求め
る。
（ｔ）の小さいシングルトーク状態の時にその値を求め
る。
【００４９】このときのα−ＥＲＬＥ特性を図２に示し
ている。この特性はタップ数Ｎには依存しない。また、
ＳＮＲが大きい程、ＥＲＬＥに対応するαが高い。その
ため、ＳＮＲが大きい場合には、大きなα値が設定さ
れ、高い収束速度が得られようになる。
ている。この特性はタップ数Ｎには依存しない。また、
ＳＮＲが大きい程、ＥＲＬＥに対応するαが高い。その
ため、ＳＮＲが大きい場合には、大きなα値が設定さ
れ、高い収束速度が得られようになる。
【００５０】（第４実施例）第４実施例の制御方法で
は、第１実施例の制御方法と第３実施例の制御方法とを
組合せてエコーキャンセラの推定値の補正を行なってい
る。
は、第１実施例の制御方法と第３実施例の制御方法とを
組合せてエコーキャンセラの推定値の補正を行なってい
る。
【００５１】この制御方法を実施するエコーキャンセラ
の機能ブロックは、α制御部11で行なわれる制御の内容
を除いては、第１実施例の場合（図１）と変わりがな
い。
の機能ブロックは、α制御部11で行なわれる制御の内容
を除いては、第１実施例の場合（図１）と変わりがな
い。
【００５２】α制御部11は、収束開始時または通話中の
回線切替えによりエコーパスが変化した時などに、その
時点での打消量ＥＲＬＥj を求め、推定値Ｈj の補正に
用いるステップゲインαj の値を式（５）から計算す
る。その後の収束進行時には式（６）によって最適なα
j の値を逐次計算する。
回線切替えによりエコーパスが変化した時などに、その
時点での打消量ＥＲＬＥj を求め、推定値Ｈj の補正に
用いるステップゲインαj の値を式（５）から計算す
る。その後の収束進行時には式（６）によって最適なα
j の値を逐次計算する。
【００５３】この場合、α制御部11自身が、残留エコー
の急増などを基にエコーパスの変化を検出するようにし
ても良いし、エコーパス変化の検出結果を他所から受け
て、α制御部11がステップゲインαj の制御を切替えて
も良い。
の急増などを基にエコーパスの変化を検出するようにし
ても良いし、エコーパス変化の検出結果を他所から受け
て、α制御部11がステップゲインαj の制御を切替えて
も良い。
【００５４】この制御方法では、収束速度が大きく、従
って、ステップゲインαj の値が式（５） αj＝ERLEj -1／(ERLEj -1＋SNR-1) の関係から大きく外れる段階では、式（５）によってα
j を再設定し、収束速度が落ち着きを見せた段階ではα
j を式（６）によって逐次制御する。そのため、通話の
開始時や通話途中でエコーの急増する事態が発生したと
きにも、エコーを素早く収束させることが可能になる。
って、ステップゲインαj の値が式（５） αj＝ERLEj -1／(ERLEj -1＋SNR-1) の関係から大きく外れる段階では、式（５）によってα
j を再設定し、収束速度が落ち着きを見せた段階ではα
j を式（６）によって逐次制御する。そのため、通話の
開始時や通話途中でエコーの急増する事態が発生したと
きにも、エコーを素早く収束させることが可能になる。
【００５５】（第５実施例）第５実施例の制御方法で
は、第２実施例の制御方法と第３実施例の制御方法とを
組合せてエコーキャンセラの推定値の補正を行なってい
る。
は、第２実施例の制御方法と第３実施例の制御方法とを
組合せてエコーキャンセラの推定値の補正を行なってい
る。
【００５６】この制御方法を実施するエコーキャンセラ
の機能ブロックは、α制御部11で行なわれる制御の内容
を除いては、第１実施例の場合（図１）と変わりがな
い。
の機能ブロックは、α制御部11で行なわれる制御の内容
を除いては、第１実施例の場合（図１）と変わりがな
い。
【００５７】α制御部11は、収束開始時または通話中の
回線切替えなどによるエコーパスの変化が検出された時
には、その時点での打消量ＥＲＬＥj を求め、推定値Ｈ
j の補正に用いるステップゲインαj の値を式（５）か
ら計算する。その後の収束進行時には式（７）によって
最適なαj の値を逐次計算する。
回線切替えなどによるエコーパスの変化が検出された時
には、その時点での打消量ＥＲＬＥj を求め、推定値Ｈ
j の補正に用いるステップゲインαj の値を式（５）か
ら計算する。その後の収束進行時には式（７）によって
最適なαj の値を逐次計算する。
【００５８】この制御方法でも、第４実施例と同様、通
話の開始時や通話途中でエコーの急増する事態が発生し
たときに、エコーを素早く収束させることができる。
話の開始時や通話途中でエコーの急増する事態が発生し
たときに、エコーを素早く収束させることができる。
【００５９】
【発明の効果】以上の実施例の説明から明らかなよう
に、本発明のエコーキャンセラ制御方法では、ステップ
ゲインαj を最適値に制御しながら、その値の大きさを
収束の進行に合わせて減少させているため、エコーキャ
ンセラは、大きい収束速度と高い推定精度との両方の特
性を持つことができ、エコーの素早い減衰と、高い精度
でのエコーの抑制が可能になる。従って、利用者は、常
にエコーの少ない快適な通話を行なうことができる。
に、本発明のエコーキャンセラ制御方法では、ステップ
ゲインαj を最適値に制御しながら、その値の大きさを
収束の進行に合わせて減少させているため、エコーキャ
ンセラは、大きい収束速度と高い推定精度との両方の特
性を持つことができ、エコーの素早い減衰と、高い精度
でのエコーの抑制が可能になる。従って、利用者は、常
にエコーの少ない快適な通話を行なうことができる。
【図１】本発明の実施例におけるエコーキャンセラの機
能ブロック図、
能ブロック図、
【図２】第３実施例のエコーキャンセラにおける打消量
ＥＲＬＥ対ステップゲインαの特性図、
ＥＲＬＥ対ステップゲインαの特性図、
【図３】第１及び第２実施例のエコーキャンセラにおけ
る時間対ステップゲインαの特性図、
る時間対ステップゲインαの特性図、
【図４】第１及び第２実施例のエコーキャンセラにおけ
る時間対打消量ＥＲＬＥの特性図、
る時間対打消量ＥＲＬＥの特性図、
【図５】従来のエコーキャンセラの機能ブロック図、
【図６】従来のエコーキャンセラにおける時間対打消量
ＥＲＬＥの特性図である。
ＥＲＬＥの特性図である。
１ 受信入力端 ２ ＰＣＭコーデック ３ エコーキャンセラ ４ 受信入力用シフトレジスタ ５ 推定インパルス応答用メモリ ６ 畳込み演算器 ７ 減算器 ８ ２線／４線変換回路 ９ 電話機 10 送信出力端 11 α制御部
───────────────────────────────────────────────────── フロントページの続き (56)参考文献 特開 平６−296119（ＪＰ，Ａ) 特開 平４−284030（ＪＰ，Ａ) 特開 平４−81133（ＪＰ，Ａ) 特開 平４−4619（ＪＰ，Ａ) 特開 平３−106231（ＪＰ，Ａ) 特開 平１−220530（ＪＰ，Ａ) 特開 平１−125030（ＪＰ，Ａ) 特開 昭63−234732（ＪＰ，Ａ) 特開 昭63−228817（ＪＰ，Ａ) 特開 昭60−240226（ＪＰ，Ａ) (58)調査した分野(Int.Cl.7，ＤＢ名) H04B 3/23 H03H 21/00
Claims (5)
- 【請求項１】 受信入力列とエコー経路のインパルス応
答の推定値とを畳込み演算して疑似エコーを発生し、こ
れをエコー信号から差し引くことによってエコーを抑圧
し、このときに得られる残留エコーの２乗値を零に近づ
けるように前記推定値を逐次補正するエコーキャンセラ
の制御方法において、 前記推定値の逐次補正におけるステップゲインαj を式 αj+1＝（Ｎ−αj）αj／（Ｎ−αj 2） （但し、Ｎは推定インパルス応答長）によって設定する
ことを特徴とするエコーキャンセラ制御方法。 - 【請求項２】 受信入力列とエコー経路のインパルス応
答の推定値とを畳込み演算して疑似エコーを発生し、こ
れをエコー信号から差し引くことによってエコーを抑圧
し、このときに得られる残留エコーの２乗値を零に近づ
けるように前記推定値を逐次補正するエコーキャンセラ
の制御方法において、 前記推定値の逐次補正におけるステップゲインαj を式 αj+1＝｛１−αj（１−αj）／Ｎ｝αj （但し、Ｎは推定インパルス応答長）によって設定する
ことを特徴とするエコーキャンセラ制御方法。 - 【請求項３】 受信入力列とエコー経路のインパルス応
答の推定値とを畳込み演算して疑似エコーを発生し、こ
れをエコー信号から差し引くことによってエコーを抑圧
し、このときに得られる残留エコーの２乗値を零に近づ
けるように前記推定値を逐次補正するエコーキャンセラ
の制御方法において、 前記推定値の逐次補正におけるステップゲインαj を式 αj＝ERLEj -1／(ERLEj -1＋SNR-1) （但し、ERLEjは時刻ｊにおける打消量、SNRはエコー信
号とエコー信号経路で加わる外乱とのレベル比）によっ
て設定することを特徴とするエコーキャンセラ制御方
法。 - 【請求項４】 前記ステップゲインαj の値が前記式 αj＝ERLEj -1／(ERLEj -1＋SNR-1) の関係から大きくずれたとき、ステップゲインαj の値
を前記式によって再設定し、その後、ステップゲインα
j の値を式 αj+1＝（Ｎ−αj）αj／（Ｎ−αj 2） （但し、Ｎは推定インパルス応答長）によって逐次設定
することを特徴とする請求項３に記載のエコーキャンセ
ラ制御方法。 - 【請求項５】 前記ステップゲインαj の値が前記式 αj＝ERLEj -1／(ERLEj -1＋SNR-1) の関係から大きくずれたとき、ステップゲインαj の値
を前記式によって再設定し、その後、ステップゲインα
j の値を式 αj+1＝｛１−αj（１−αj）／Ｎ｝αj （但し、Ｎは推定インパルス応答長）によって逐次設定
することを特徴とする請求項３に記載のエコーキャンセ
ラ制御方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP06000595A JP3145269B2 (ja) | 1995-02-24 | 1995-02-24 | エコーキャンセラ制御方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP06000595A JP3145269B2 (ja) | 1995-02-24 | 1995-02-24 | エコーキャンセラ制御方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH08237174A JPH08237174A (ja) | 1996-09-13 |
| JP3145269B2 true JP3145269B2 (ja) | 2001-03-12 |
Family
ID=13129552
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP06000595A Expired - Fee Related JP3145269B2 (ja) | 1995-02-24 | 1995-02-24 | エコーキャンセラ制御方法 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3145269B2 (ja) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6580793B1 (en) * | 1999-08-31 | 2003-06-17 | Lucent Technologies Inc. | Method and apparatus for echo cancellation with self-deactivation |
| JP2002009675A (ja) * | 2000-06-26 | 2002-01-11 | Nec Corp | 多重化回線用エコー除去方法及び装置 |
| JP5061853B2 (ja) | 2007-11-06 | 2012-10-31 | 沖電気工業株式会社 | エコーキャンセラ及びエコーキャンセルプログラム |
| JP5061976B2 (ja) | 2008-03-12 | 2012-10-31 | 沖電気工業株式会社 | エコーキャンセラ、エコーキャンセル方法及びプログラム |
| JP5223576B2 (ja) | 2008-10-02 | 2013-06-26 | 沖電気工業株式会社 | エコーキャンセラ、エコーキャンセル方法及びプログラム |
| JP5493817B2 (ja) | 2009-12-17 | 2014-05-14 | 沖電気工業株式会社 | エコーキャンセラ |
-
1995
- 1995-02-24 JP JP06000595A patent/JP3145269B2/ja not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| JPH08237174A (ja) | 1996-09-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP2836277B2 (ja) | エコーキャンセル装置 | |
| JP2518433B2 (ja) | ダブルト―ク検出回路 | |
| US7856097B2 (en) | Echo canceling apparatus, telephone set using the same, and echo canceling method | |
| CA2075843C (en) | Echo canceller | |
| WO1995031052A1 (fr) | Compensateur d'echos et procede d'estimation de la trajectoire de l'echo | |
| US6185299B1 (en) | Adaptive echo cancellation device in a voice communication system | |
| JP3145269B2 (ja) | エコーキャンセラ制御方法 | |
| JP5293952B2 (ja) | 信号処理方法、信号処理装置、並びに信号処理プログラム | |
| US6580794B1 (en) | Acoustic echo canceler with a peak impulse response detector | |
| JP3244416B2 (ja) | エコーキャンセラ | |
| JP3002374B2 (ja) | エコーキャンセラと併用する音声スイッチの制御方法 | |
| JP3452341B2 (ja) | エコーキャンセラ | |
| JP2551869B2 (ja) | 反響消去装置 | |
| JP2907706B2 (ja) | エコーキャンセラと併用する音声スイッチの制御方式 | |
| JP3147457B2 (ja) | エコーキャンセラ制御法 | |
| JP2841952B2 (ja) | エコー消去装置 | |
| JPH0946277A (ja) | エコーキャンセラ装置 | |
| JP3121969B2 (ja) | 音響反響除去装置 | |
| JP2516260B2 (ja) | エコ―キャンセラ | |
| JP3152815B2 (ja) | 音響反響除去装置 | |
| JPH05206898A (ja) | エコー消去装置 | |
| JPH07303063A (ja) | エコーキャンセラおよびエコーパス推定方法 | |
| JP2693183B2 (ja) | エコーキャンセラ | |
| JP3121988B2 (ja) | 音響反響除去装置 | |
| JPH05252079A (ja) | エコーキャンセラ |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| LAPS | Cancellation because of no payment of annual fees |