(57)【要約】  ａ）第１の信号とこの第１の信号のエコーを有する第２の信号とを受け取るステップと、ｂ）有限インパルス応答フィルタを用いて第１の信号をフィルタリングすることによりエコー複製信号を生成するステップであって、有限インパルスフィルタはエコー複製信号を生成するためにフィルタ係数ベクトルを用いる当該ステップと、ｃ）第２の信号からエコー複製信号を減ずることによりエコー消去信号を生成するステップと、ｄ）第２の信号内の干渉の大きさを決定するステップと、ｅ）干渉の大きさに依存してステップベクトルを決定するステップであって、干渉の大きさの増大はステップベクトルのサイズを小さくする当該ステップと、ｆ）ステップベクトルによりフィルタ係数ベクトルを更新するステップとを有するエコーキャンセリング方法である。
2007522720.tif 000002
【特許請求の範囲】
【請求項１】
ａ）第１の信号とこの第１の信号のエコーを有する第２の信号とを受け取るステップと、
ｂ）有限インパルス応答フィルタを用いて前記第１の信号をフィルタリングすることによりエコー複製信号を生成するステップであって、前記有限インパルスフィルタは前記エコー複製信号を生成するためにフィルタ係数ベクトルを用いる当該ステップと、
ｃ）前記第２の信号から前記エコー複製信号を減ずることによりエコー消去信号を生成するステップと、
ｄ）前記第２の信号内の干渉の大きさを決定するステップと、
ｅ）前記干渉の大きさに依存してステップベクトルを決定するステップであって、前記干渉の大きさの増大は前記ステップベクトルのサイズを小さくする当該ステップと、
ｆ）前記ステップベクトルにより前記フィルタ係数ベクトルを更新するステップと
を有するエコーキャンセリング方法。
【請求項２】
前記第２の信号が干渉により影響を及ぼされない場合に、前記フィルタ係数を改善する第１のステップベクトルを生成することと、
前記干渉の大きさに依存して第２のステップベクトルを生成することと、
前記第１及び第２のステップベクトルから前記ステップベクトルを選択し、最も小さいサイズを持つ前記ステップベクトルが選択されることと
により前記ステップベクトルが決定される請求項１記載のエコーキャンセリング方法。
【請求項３】
前記第２のステップベクトルは、前記第２の信号の干渉の無い状態で前記第１のステップベクトルに実質的に対応する請求項２記載のエコーキャンセリング方法。
【請求項４】
前記干渉の大きさが大きくなると、前記第２のステップベクトルは連続的に減少し、前記第２のステップベクトルの前記サイズが前記第１のステップベクトルの前記サイズよりも小さくなる請求項３記載のエコーキャンセリング方法。
【請求項５】
前記干渉の大きさは、前記エコー消去信号のレベルを用いて決定される請求項１ないし４のいずれか１項に記載のエコーキャンセリング方法。
【請求項６】
前記ステップベクトルは合計のエコーリターンロスに依存して決定される請求項１ないし５のいずれか１項に記載のエコーキャンセリング方法。
【請求項７】
ダブルトーク状態が存在するか否かを検出するステップと、
前記ダブルトーク状態が存在するか否かに依存して前記合計のエコーリターンロスを異なるように決定するステップと
を有する請求項５又は６記載のエコーキャンセリング方法。
【請求項８】
前記干渉の大きさは前記エコー消去信号の重み付けされたレベルを用いて決定され、前記エコー消去信号の前記重み付けされたレベルは、前記合計のエコーリターンロスを乗じた前記エコー消去信号の前記レベルよりなる請求項６又は７記載のエコーキャンセリング方法。
【請求項９】
第１の信号を受け取り、エコー複製信号を出力する有限インパルス応答フィルタであって、前記エコー複製信号を生成するためにフィルタ係数ベクトルを用いる当該有限インパルス応答フィルタと、
前記第１の信号のエコーを有する第２の信号から前記エコー複製信号を減ずる減算回路と、
ステップベクトルにより前記フィルタ係数ベクトルを繰り返し更新する係数更新回路と、
前記第２の信号内の干渉の大きさを決定する手段と、
前記干渉の大きさに依存して前記ステップベクトルを決定する手段であって、前記干渉の大きさの増大は前記ステップベクトルのサイズを連続的に小さくする当該手段と
を有するエコーキャンセリング装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、エコーキャンセリング方法及びエコーキャンセリング装置に関する。
【背景技術】
【０００２】
加入者線の受信経路と送信経路との間にカップリングが存在すると、電気通信システムにおいて認知可能なエコーが常に存在する。上記カップリングのために、一方の加入者の受信した信号が、他方の加入者に送信し返され、エコーが生じる。すなわち、他方の加入者は自分自身の信号の遅れているバージョンに気付く。通信の認知性（perceptibility）は上記エコーにより強く低下し、特に遅延が大きい場合に低下する。従って、局所的にエコーを自動的に消去（cancel）する方法及び装置を提供することが望ましい。
【０００３】
エコーの発生の異なる態様のために、音響エコーと電気的エコーとを区別することができる。音響エコーは、テレビ会議システム又はビデオ会議システムのためのハンズフリー電話機のような実際の電気通信の端末のスピーカとマイクロフォンとの音響カップリングから生じる。電気的エコーは、一方の加入者の端末、例えばコードレス電話システムの基地局と２ワイヤの電話線との接続を実現する所謂「２／４ワイヤ」ハイブリッド回路において生成される。技術的な理由のために、これらのハイブリッド回路は、送信経路と受信経路との完全な分離を与えることができない。従って、送信される信号の一部が受信経路にリークし、電気的なラインエコーになる。
【０００４】
図１は、最新技術によるラインエコーキャンセラ（ＬＥＣ）を示している。このラインエコーキャンセラＬＥＣは、２／４ワイヤラインインターフェースにおけるハイブリッドのインピーダンスミスマッチにより生じるエコーを消去する。
【０００５】
図２は、最新技術による音響エコーキャンセラ（ＡＥＣ）を示している。スピーカにより生成される音響信号が、音響エコー経路を介してマイクロフォンに反射される。音響エコーキャンセラＡＥＣに課せられた仕事は、マイクロフォンに反射される信号を抑制することである。
【０００６】
上記両方のエコーキャンセラは、適応フィルタと制御回路とを有している。上記適応フィルタは、２つのノードＳＮＤout，ＲＣＶinの間のエコーを整合させなければならない。ＳＮＤinは基準信号を構成しており。この基準信号からエコー信号が生成される。上記基準信号は、エコーキャンセラに与えられる。上記エコーキャンセラから出力される信号はＳＮＤoutと呼ばれる。上述の信号の命名法は、よく用いられる概念「遠端信号」及び「近端信号」との混乱を避けるべきである。
【０００７】
デジタル／アナログコンバータ（Ｄ／Ａコンバータ）は、デジタル信号ＳＮＤoutをエコー経路を通って受信経路に一部反射されるアナログ信号に変換する。アナログ／デジタルコンバータ（Ａ／Ｄコンバータ）は、エコー信号を含む相手側からの受信信号をＲＣＶinと呼ばれるデジタル信号に変換する。上記適応フィルタは、ＲＣＶinに含まれるエコー信号に対応する信号を生成しなければならない。減算器（＋／−）は、信号ＲＣＶinから上記適応フィルタにより生成される信号を減ずる。従って、残留エコー成分のみを伴う信号ＲＣＶoutが生成される。この誤差信号は、上記適応フィルタの性能を改善し、それに応じてフィルタを変更するために用いられる。制御ユニットは、速度の適応プロセスに影響を及ぼすため、更には音響エコーキャンセラの場合には上記適応フィルタに従う減衰スイッチを制御するために、上記エコーキャンセラのノードの全ての信号を追跡し、他の制御信号（図１及び図２の破線）を導き出す。上記適応フィルタの必要とされる高いフィルタの等級は、適度な限度内において上記取り組みを維持するために最大になされなければならない（室内音響の残響時間は、典型的には２０から２００マイクロ秒までの範囲にあることに注意されたい。）ので、これらの追加の対策が音響エコー消去のために必要であり得る。
【０００８】
従来、適応フィルタとして有限インパルス応答フィルタ（ＦＩＲフィルタ）が用いられている。このＦＩＲフィルタは、入力基準信号ベクトル
2007522720.tif 000003
と係数ベクトル
2007522720.tif 000004
とのコンボリューションにより表される。上記入力信号は、ある数の連続する入力値に関してサンプリングされる。これらの値が、上記入力基準信号ベクトル
2007522720.tif 000005
を構成する。
【０００９】
上記ＦＩＲフィルタにより出力される信号は、以下のように計算される。
【数１】
2007522720.tif 000006
Ｙはスカラーを表している。Ｙは、時間ｋにおいて異なるポイントで連続的に出力される。従って、式１は以下の形に書き換えられ得る。
【数２】
2007522720.tif 000007
【００１０】
各サンプル期間ｋにおける上記フィルタ係数ベクトル
2007522720.tif 000008
は、ベクトルのインクリメント
2007522720.tif 000009
により連続的に更新される。一般的な適応方程式は
【数３】
2007522720.tif 000010
である。
【００１１】
ステップベクトル
2007522720.tif 000011
のノルム
2007522720.tif 000012
は、アルゴリズムの安定性及び適応速度を決定する。最もよく用いられる代表的な適応アルゴリズムは、最急降下法である。上記ステップベクトルの方向は、自乗誤差信号の実際の負の勾配に対応する。最適な係数ベクトルは、平均自乗誤差信号の勾配が０である、すなわち平均自乗誤差信号が最小値に達する場合に達せられる。正規化最小平均自乗アルゴリズム（ＮＬＭＳアルゴリズム）は、以下の式により定義される。
【数４】
2007522720.tif 000013
ｅ（ｋ）は、時間ｋのポイントにおけるエコーとフィルタの出力（エコーの複製）との間の実際の誤差を表している。正規化、すなわち
2007522720.tif 000014
で割ることは、上記入力基準信号のパワーからステップベクトル
2007522720.tif 000015
の依存性を取り除く。
【００１２】
減算回路（＋／−）後の誤差信号ｅ（ｋ）は、出力信号ＲＣＶoutに対応する。誤差信号ｅ（k)は、いかなる干渉も無い残留エコー信号に対応する。これは、信号ＲＣＶinが専ら上記エコー信号を含んでいると仮定している。入力信号ＲＣＶinがエコー信号と他の信号との重ね合わせである場合、誤差信号ｅ（ｋ）はＦＩＲフィルタの誤差の正しい複製ではない。この状況の広く知られている表現法は「ダブルトーク」である。
【００１３】
図２の音響エコーキャンセラＡＥＣでは、すなわち、人間がマイクロフォンに話し、スピーカが音を生成していると、エコー信号と他のノイズ信号との重ね合わせが生じる。上記マイクロフォンに入力される信号は、エコー信号と話者の声との重ね合わせである。フィルタ係数ベクトル
2007522720.tif 000016
の適応は、ダブルトーク中抑制されなければならない。上記エコー複製フィルタ（ＦＩＲフィルタ）はダブルトーク状態の間作用するが、フィルタ係数ベクトル
2007522720.tif 000017
はダブルトークの間には更新されない。ダブルトークの間にエコー特性（エコー経路）が変化しない場合、上記ＦＩＲフィルタの出力はエコー信号を正確に表現するものである。上記エコー信号は、エコー信号と干渉している話者との重ね合わせを構成する入力信号から減じられ得る。このケースでは、話者の音声信号は著しく劣化する。
【００１４】
図３は、ダブルトークの間にフィルタの係数ベクトル更新の禁止のための信号を与える従来のエコーキャンセラを示している。上記ＦＩＲフィルタは、図３ではエコー複製フィルタと呼ばれている。上記フィルタ係数ベクトルを更新する手段は、図３ではフィルタ係数更新部と呼ばれている。ＮＬＭＳは、正規化最小平均自乗アルゴリズムを意味している。これは、以前の段落において紹介されたアルゴリズムである。ＳＩは信号の入力を意味し、ＳＯは信号の出力を意味している。出力信号は、ＳＯによりスピーカによって出力されると仮定され得る。マイクロフォンは、入力信号ＲＩを記録すると仮定され得る。入力信号ＲＩは、エコー信号と妨害する干渉信号との重ね合わせである。ＰＷＦは、プリホワイトニングフィルタを意味している。フィルタＰＷＦは単にオプションであり、エコーキャンセラの機能はプリホワイトニングフィルタＰＷＦに関係なく説明され得る。ＩＰＷは、逆（inverse）プリホワイトニングフィルタを意味している。
【００１５】
上記信号入力ＳＩはプリホワイトニングフィルタを介してエコー複製フィルタ（ＦＩＲフィルタ）に与えられ、その作用はここでは無視され得る。上記エコー複製フィルタ（ＦＩＲフィルタ）は、信号ＲＩから減じられたエコーの複製すなわち出力信号を生成する。結果として得られる信号はＲＯである。上記エコーの複製が正しく、妨害する干渉信号がＲＩに重ね合わせられていない場合、ＲＯはゼロに収束すべきである。ダブルトークが検出される毎にフィルタ係数の更新を禁止する制御回路が設けられている。Ｙ（ｋ）は、エコー複製フィルタの時刻ｋにおける出力サンプルである。ＳＩ（ｋ）は、エコー複製フィルタへの基準信号入力サンプルである。上記エコー複製フィルタの出力は、
【数５】
2007522720.tif 000018
に等しい。
【００１６】
定数Ｎは、ＦＩＲフィルタのフィルタ長のそれぞれの係数の数に対応する。図３において減算回路（＋／−）により出力される信号は、誤差信号ＲＯ(k)に等しい。
ＲＯ（ｋ）＝ＲＩ（ｋ）−Ｙ（ｋ）  （６）
【００１７】
ステップベクトル
2007522720.tif 000019
は、以下の式により決定される。
【数６】
2007522720.tif 000020
αは、制御パラメータであり、実現されるアルゴリズムにより予め決定された安定性の範囲内において適応速度を制御するために用いられる。式（５）ないし（７）により与えられるアルゴリズムの場合には、αは０ないし２の範囲になっていなければならない。別個の制御アルゴリズムは制御パラメータαを導き、ｍは通常の時間の指標ｋ（サンプリングレート）に対するサブサンプリングされた時間の指標である。
【００１８】
信号ＳＩ，ＲＩ及びＲＯは、回路ＳＴＬ及びＬＴＬを介して制御回路に導かれ、これらは異なる時間解像度を用いてレベル計算を行う。回路ＳＴＬは、信号ＳＩ，ＲＩ及びＲＯそれぞれの大きさの短時間減衰（leaky）積分を行う。回路ＬＴＬは、これらの入力信号の長時間減衰積分を行い、制御ブロックの内部状態を介して更に制御される。
【００１９】
所定の期間の信号の積分は、この時間窓内に信号を平均することにより行われ得る。上記時間窓に先行する又は続く信号の値が平均に寄与するべきではない場合、フィルタリング又はコンボリューションが、所定の期間は値１を有し、その他の場合は値０を有するフィルタ関数により行われる。しかしながら、上述したボックス関数以外の他の関数も用いられ得る。
【００２０】
アナログの世界では、信号ｘ（ｔ）の大きさの減衰積分に関する演算子は、
【数７】
2007522720.tif 000021
として定義され得る。
【００２１】
αは０に近づくと、減衰積分は、通常の積分になる。e−α・τは、はるか以前の｜x(t)｜の値が積分に著しく寄与しないという効果を持たせる。長時間減衰積分と短時間減衰積分との違いは、単に異なるαの値により決定される。短時間減衰積分は、αに関して大きい値を有する。ｘ（ｔ）の最近の値のみが積分に寄与する。
【００２２】
離散的な世界では、減衰積分は再帰方式（recursion scheme）
Ｌ＿ｘ（ｋ）＝（１−β）・Ｌ＿ｘ（ｋ−１）＋β・｜ｘ（ｋ）｜  （９）
をもたらす。ｘは図３の信号ＳＩ，ＲＩ又はＲＯであり得る。長時間レベルＬＬ＿ＳＩ，ＬＬ＿ＲＩ及びＬＬ＿ＲＯは、それぞれ、信号ＳＩ，ＲＩ及びＲＯの長時間減衰積分を表す。短時間レベルＬ＿ＳＩ，Ｌ＿ＲＩ及びＬ＿ＲＯは、それぞれ、信号ＳＩ，ＲＩ及びＲＯの短時間減衰積分を表す。サンプルｘ（ｋ）の絶対値（大きさレベル）を減衰積分する代わりに、自乗されたサンプルｘ２（ｋ）（パワーレベル）が減衰積分されてもよい。図３に示されているローカルなエコーキャンセラの制御回路は、２つの状況においてダブルトーク状態を検出する。以下の式が成り立つと、第１のダブルトーク状態が検出される。
Ｌ＿ＳＩ（ｋ）＜Ｌ＿ＲＩ（ｋ）・ＤＢＬ＿ＴＨ１  （１０）
【００２３】
解釈は以下の通りである。すなわち、レベルＬ＿ＲＩが、エコー経路に沿って信号レベルＬ＿ＳＩ（ｋ）から生成される最悪のケースのエコーレベルよりも大きいと、第１のダブルトークの状態に関する指標がマークされる。その結果、フィルタの係数の更新が中断される。
【００２４】
第２のダブルトーク状態の検出は、以下の式により行われる。
【数８】
2007522720.tif 000022
【００２５】
解釈は以下の通りである。すなわち、信号ＲＩと信号ＲＯとの間の現在の減衰である短時間レベルの比Ｌ＿ＲＩ（ｋ）／Ｌ＿ＲＯ（ｋ）が、適応フィルタにより実現されるエコーの減衰の信頼性の高い推定であるべき長時間レベルの比ＬＬ＿ＲＩ（ｍ）／ＬＬ＿ＲＯ（ｍ）よりも小さいと、第２のダブルトーク状態に関する指標がマークされる。その結果、フィルタの適応も禁止される。
【００２６】
より長い時刻の間上記第２のダブルトークの状態のみが続き、アクティブな信号ＳＩが存在する場合、真のダブルトーク状態ではなく、エコー経路の変化に関する確率が与えられる。これは、長いダブルトークが原則ではなく例外であることを仮定している。所定の期間（タイムアウト）の後、たとえ上記第２のダブルトーク状態が続いても、フィルタの更新が再開される。代わりに、エコー経路のシステムの変更が仮定される。推定されたエコーの減衰はもはや信頼性が高くないので、長時間レベルＬＬ＿ＲＩ（ｍ）及びＬＬ＿ＲＯ（ｍ）は現在の短時間レベルＬ＿ＲＩ（ｋ）及びＬ＿ＲＯ（ｋ）にリセットされる。
【００２７】
最新技術によるダブルトーク状態を検出するアルゴリズムは、幾つかの欠点を有している。ダブルトーク状態があまりにも遅く検出されると、フィルタ係数ベクトル
2007522720.tif 000023
は誤差信号に基づいてかなりの量の時間変化し、これは全く間違っている。上記フィルタ係数は、更新が阻止されるまで低下する。ダブルトークの場合には、エコーの減衰又はエコーの消去が非常に不良である。また、たとえダブルトーク状態がもはや存在しない場合であっても、フィルタ係数ベクトル
2007522720.tif 000024
の更新が阻止され得る。この状態では、上記フィルタ係数ベクトルの改善が望ましいが、行われない。ダブルトーク状態を検出する上記アルゴリズムは、１００％の信頼性ではなく、従ってフィルタ係数の適応は、音声通信に関する悪い質及び理解のし易さの影響により誤って行われることが多い。
【発明の開示】
【発明が解決しようとする課題】
【００２８】
本発明の目的は、上記最新技術の欠点を克服するエコーキャンセリング方法及び対応するエコーキャンセリング装置を提供することにある。
【課題を解決するための手段】
【００２９】
上記課題は、
ａ）第１の信号とこの第１の信号のエコーを有する第２の信号とを受け取るステップと、
ｂ）有限インパルス応答フィルタを用いて前記第１の信号をフィルタリングすることによりエコー複製信号を生成するステップであって、前記有限インパルスフィルタは前記エコー複製信号を生成するためにフィルタ係数ベクトルを用いる当該ステップと、
ｃ）前記第２の信号から前記エコー複製信号を減ずることによりエコー消去信号を生成するステップと、
ｄ）前記第２の信号内の干渉の大きさを決定するステップと、
ｅ）前記干渉の大きさに依存してステップベクトルを決定するステップであって、前記干渉の大きさの増大は前記ステップベクトルのサイズを小さくする当該ステップと、
ｆ）前記ステップベクトルにより前記フィルタ係数ベクトル
2007522720.tif 000025
を更新するステップと
を有するエコーキャンセリング方法により解決される。
【００３０】
最新技術による方法と本発明による方法との決定的な違いは、ダブルトーク状態が検出された場合に、フィルタ係数ベクトル
2007522720.tif 000026
の適応が中断されないことである。代わりに、干渉の大きさが決定される。上記干渉の大きさは、ダブルトーク状態に関する確率を表すとみなされ得る。ダブルトークが本当らしいと上記干渉の大きさの値が高く、ダブルトークが本当らしくないと上記干渉の大きさの値が低い。上記フィルタ係数ベクトル
2007522720.tif 000027
の適応は連続的に行われるが、ダブルトーク状態の可能性に依存して減速（slow down）される。干渉の大きさの増大は、上記ステップベクトルのサイズを自動的に小さくする。閾値が達せられると、上記適応は停止されない（硬調なスイッチング）が、連続的に減速される（軟調なコントロール）。上記エコー消去信号のレベルの突然の上昇は、例えば干渉の存在に関する指標として用いられ得る。しかしながら、全体として比較のために閾値を与えること及び適応を中断することの代わりに、本発明の方法は、上記干渉の大きさに依存してステップベクトルのサイズを連続的に小さくする。その結果、閾値に達する前に、全体として適応プロセスは増大した干渉の大きさに反応する。あまりにも遅く検出されるダブルトーク状態のためにエコーキャンセリング能力の大きな低下は生じない。
【００３１】
原則的に、上記干渉の大きさが大きくなると上記ステップベクトルのサイズは小さくされ、上記干渉の大きさが小さくなると上記ステップベクトルのサイズは大きくされる。しかしながら、ステップベクトルの増大は、適応の安定性を保障する上限により好ましく制限される。それにより、適応速度は、高い干渉のレベルの間小さくなり、低い干渉のレベルの間大きくなる。係数ベクトル
2007522720.tif 000028
の変化は、ダブルトーク状態の間はごくわずかであるべきであり、エコー消去は大きく劣化しない。
【００３２】
上記ステップベクトルは、第２の信号が干渉によって影響を及ぼされないようにフィルタ係数ベクトルを改善する第１のステップベクトルを生成することと、上記干渉の大きさに依存して第２のステップベクトルを生成することと、上記第１及び第２のステップベクトルから上記ステップベクトルを選択し、その際、最も小さいサイズを持つステップベクトルが選択されることとにより決定されることが好ましい。
【００３３】
上記第１のステップベクトルの生成は、通常のＮＬＭＳアルゴリズムに従って行われ得る。上記第２のステップベクトルの生成は、干渉の大きさを考慮に入れる。上記干渉の大きさが大きくなると、第２のステップベクトルは連続的に小さくなる。従って、ダブルトークに関する可能性が高い場合、第２のステップベクトルは第１のステップベクトルよりも小さい。結果として得られるステップベクトルは、第１及び第２のステップベクトルのうちの小さい方に対応するので、ステップベクトルは増大する干渉の大きさに依存して連続的に小さくなる。逆に、減少する干渉の大きさは、高すぎる適応速度の危険性を伴って第２のステップベクトルを大きくするが、第１及び第２のステップベクトルのうちの小さい方を再度用いることにより、適応の安定性が保障される。
【００３４】
干渉信号が存在しない場合であっても、上記第２のステップベクトルが同様に第１のステップベクトルよりも著しく小さいと、フィルタ係数ベクトルの適応速度は著しく低下する。従って、第２のステップベクトルは、第２の信号内においていかなる干渉も無い状態で第１のステップベクトルにほぼ対応すべきであるように低下することが好ましい。
【００３５】
上記干渉の大きさは、エコー消去信号のレベルに依存して決定され得る。エコー複製信号が第１の信号のエコーに十分に対応する場合、エコー消去信号の高いレベルは第２の信号内の干渉に関する指標である。正しいテンポでの特定のポイントに関するレベルの計算は、ある時間窓における信号サンプルの大きさ（大きさのレベル）又は自乗された信号サンプル（パワーのレベル）の積分からもたらされる。この計算を正しいテンポで連続的に行うことにより、大きさ又は自乗された信号サンプルへのスムージング又はローパスフィルタリングプロセスが与えられる。上記フィルタリングは、例えばＦＩＲ矩形窓関数（通常の短時間積分）又はＩＩＲ指数窓関数（式９に従う減衰長時間積分）を用いて行われ得る。８ｋＨｚのサンプリングレートで、式（９）の係数βは、約４ないし１６マイクロ秒の時定数を意味する２−５ないし２−７の範囲内にあり得る。好ましくは、上記短時間レベルの減衰積分に関する時定数は、エコー経路インパルス応答の時間窓をカバーすべきである。
【００３６】
上記干渉の大きさは、また、全部の（overall）又は合計の（total）エコーリターンロス（ＥＲＬＴ）を用いて決定され、これは、エコー経路のエコーリターンロス（ＥＲＬ）及びエコーキャンセラのエコーリターンロスの増大（ＥＲＬＥ）を有している。高いエコーリターンロスの増大は、十分に機能するエコーキャンセラを示している。干渉の無い状態において、適応アルゴリズムは速く収束し、高いＥＲＬＥ値をもたらすべきである。その場合、その後干渉が存在する際、高いＥＲＬＥ値はステップベクトルのサイズをより小さくすべきである。
【００３７】
上記エコーリターンロスの推定は、ある時間における２つの信号の長時間レベルの比により決定され、例えば、第１の及び第２の信号の長時間レベルの比によるエコー経路のエコーリターンロス、第２の信号とエコー消去信号との長時間レベルの比によるエコーリターンロスの増大、及び第１の信号とエコー消去信号との長時間レベルの比による全部の又は合計のエコーリターンロスである。エコーリターンロスの推定は、ダブルトーク状態が存在するか否かに依存して決定されることが好ましい。エコーの消去が十分に行われる場合、ダブルトークの間、エコー消去信号は主にダブルトークの干渉により構成される。従って、上記比はダブルトークの場合におけるエコーの減衰を表さない。ダブルトークの間に対応する長時間レベルから直接エコーリターンロスを決定することをできる限り妨げることが好ましい。この目的のために、ダブルトークが検出され、それに応じてエコーリターンロスの計算が妨げられる。あるステップに更に進むために、エコーリターンロスの計算の推定が、ダブルトークを装うシステムの変化の場合に、再適応の長いブロッキングを避けるためにより低いエコーの減衰の方向に推定を引き寄せる又はドライブするよう、特別なリーク係数を用いてダブルトークの異なる状態の間に修正される。
【００３８】
上記エコー消去信号の短時間レベルと合計のエコーリターンロスとは、干渉の大きさを決定するために同時に両方用いられ得る。上記干渉の大きさは、合計のエコーリターンロスを乗じたエコー消去信号の短時間レベルに比例することが好ましい。全てのレベルの計算が信号サンプルの大きさに依存して行われると、上記干渉の大きさが更に自乗される。自乗された信号サンプルに依存するパワーのレベルの計算が選択されると、これはスキップされ得る。この精神において、干渉の大きさを自乗することは、第１及び第２のステップベクトルのサイズに関して同じユニットを得るための変形（適応的動作）である。
【００３９】
上記エコー消去信号の突然の立ち上がりは、フィルタ係数の適応を急速に減速させる。この作用は、ここでは実際に達せられるエコーの減衰と共にますます強くなる。これは、フィルタが良好にエコーを減衰させるほどダブルトーク状態の自動化された制御のプロセスがより急速に及びより信頼性高く作用することを意味する。まさにこの場合には、ダブルトークの制御が機能しなくなり、フィルタベクトルが劣化すると、多くの減衰を放つことができる。それにもかかわらず、ダブルトークの制御機能の良好な開始ポイントは、例えば６ｄＢのできる限り小さい全部の又は合計のエコーロスにより既に達せられている。
【発明を実施するための最良の形態】
【００４０】
本発明の好ましい実施の形態が、添付の図面を参照して以下に説明される。
【００４１】
図４のエコーキャンセリング装置は、有限インパルス応答フィルタ（ＦＩＲフィルタ）と、減算回路
2007522720.tif 000029
と、ステップベクトル
2007522720.tif 000030
を決定する手段（ステップベクトル回路）と、干渉測定回路とを有している。上記ＦＩＲフィルタは基準信号ＳＩを受け取り、この基準信号ＳＩはノードＳＯにおいて出力され、エコー経路を介してノードＲＩに反射される。信号ＲＩは、エコー信号と適応プロセスを妨害する干渉信号との重ね合わせを構成する。上記ＦＩＲフィルタ及び減算回路に課された仕事は、信号ＲＩ内のエコー信号を消去するが、干渉信号を無視することである。上記減算回路の出力は、エコー消去された信号ＲＯである。ＦＩＲフィルタは、入力として基準信号ＳＩを有し、信号Ｙを出力する。エコー複製信号Ｙは、式
【数９】
2007522720.tif 000031
により計算される。
【００４２】
干渉信号が存在しない場合、エコー複製信号Ｙは、収束し、エコー信号を表す。エコー消去信号ＲＯは０に近づくべきである。エコー消去信号ＲＯは、上記干渉測定回路に入力され、この干渉測定回路は信号ＳＩを受け取りもする。干渉測定回路は、信号ＳＩをエコー消去信号ＲＯと比較することにより干渉の大きさを計算する。上記干渉の大きさは、エコー消去信号ＲＯがいかに強く干渉信号により影響を及ぼされているかを示す。計算された干渉の大きさは、ステップベクトル
2007522720.tif 000032
を計算するステップベクトル回路に出力される。干渉の大きさが大きくなると、ステップベクトル
2007522720.tif 000033
は徐々に小さくされる。従って、干渉の大きさが高いと、フィルタ係数ベクトルの適応は減速される。ステップベクトル
2007522720.tif 000034
は、係数更新回路に出力される。この係数更新回路は、フィルタの係数ベクトル
2007522720.tif 000035
を計算し、新しいフィルタ係数ベクトル
2007522720.tif 000036
をＦＩＲフィルタに出力する。ＦＩＲフィルタは、上記
新しいフィルタ係数ベクトル
2007522720.tif 000037
を用いてＹに関する次の値を代わりに計算する。新しい係数ベクトル
2007522720.tif 000038
は、
2007522720.tif 000039
に等しい。
【００４３】
図５は、本発明の第２の実施の形態を示している。図５のエコー複製フィルタは、図４のＦＩＲフィルタに対応している。フィルタ係数更新回路は、図４のステップベクトル回路及び係数更新回路に対応し、ひとまとめにされている。図４の干渉レベル回路は、図５のＥＲＬ回路及びノルム選択回路に対応している。図５のエコー複製フィルタは、プリホワイトニングフィルタＰＷＦを介して基準信号ＳＩを受け取る。上記プリホワイトニングフィルタはプリエンファシスフィルタであり、エコーキャンセリング装置の全体の収束の挙動を改善するために、音声信号のより高いスペクトル部分を持ち上げる（lift up）。エコーキャンセリングアルゴリズムが、フラットなスペクトル（例えば白色ノイズ）を伴うトレーニング信号に関してより良好に実行し、あるトレーニング信号のスペクトルギャップを持ち上げることは良いプラクティスであることはよく知られている。上記プリホワイトニングフィルタの出力及びその挙動は、以下の式により説明され得る。
ＰＷＦ（ｋ）＝ｇ・[ｘ（ｋ）＋（α−１）・ｘ（ｋ−１）] ／ａ  （１３） 
【００４４】
定数ａはフィルタの係数であり、ｇはプリホワイトニングフィルタのゲインファクタである。ｘはプリホワイトニングフィルタに入力される信号であり、ＰＷＦはプリホワイトニングフィルタの出力である。上記ＰＷＦフィルタは、図５のフィルタ係数更新プロセスの動作を改善するが、それは必須の要素ではない。上記ＰＷＦフィルタはエコー複製フィルタの入力部において与えられ、エコー信号ＲＩもまたプリホワイトニングフィルタを用いてフィルタリングされる。更に、図５の回路には逆プリホワイトニングフィルタが設けられている。この逆プリホワイトニングフィルタＩＰＷは、以下の式により説明され得る。
ＩＰＷ（ｋ）＝（１−ａ）・ＩＰＷ（ｋ−１）＋ｇ・ｘ（ｋ）／ａ  （１４）
【００４５】
上記エコー消去信号は、エコー消去信号に対するプリホワイトニングフィルタの影響を取り除くために、逆プリホワイトニングフィルタを通って送られる。
【００４６】
図５のエコーキャンセリング装置は、幾つかの短時間レベルの減衰積分フィルタＳＴＬを有している。信号ＳＩ，ＲＯ及びＲＩは、短時間レベルの減衰積分フィルタＳＴＬに送られる。これらのフィルタの挙動は、以下の式により説明される。
【数１０】
2007522720.tif 000040
【００４７】
その他の場合は、上記式においてｘは短時間レベルの減衰積分フィルタへの入力を表しており、従って、ｘは、どのフィルタが関係するかに依存してＳＩ，ＲＯ又はＲＩに等しい。入力信号レベルの立ち上がり又は立下りに依存して、ＳＴＬフィルタの出力を計算するために異なる定数βが用いられる。現在の入力信号がＳＴＬフィルタの以前の出力値よりも高いと、βＬＰが用いられる。音声信号の開始に関して特に重要視するために時定数パラメータβの異なる設定が用いられる。更に、図５は幾つかの長時間減衰積分フィルタＬＴＬを有している。このＬＴＬフィルタの出力は、以下の式によりある間引き率で計算される。
【数１１】
2007522720.tif 000041
【００４８】
ＬＬ＿ＳＩ（ｍ）は、図５において短時間レベルの減衰積分フィルタＳＴＬを介して基準信号ＳＩに接続された長時間レベルの減衰積分フィルタＬＴＬの出力を意味している。ＬＬ＿ＲＩ（ｍ）は、短時間レベルの減衰積分フィルタＳＴＬを介してエコー信号ＲＩに接続された長時間レベルの減衰積分フィルタＬＴＬの出力を意味している。ＬＬ＿ＲＯ（ｍ）は、上記短時間レベルの減衰積分フィルタを介してエコー消去信号ＲＯに接続された上記長時間レベルの減衰積分フィルタの出力を意味している。図５では状態制御回路も設けられており、これはダブルトーク状態を検出するために用いられる。以下の式において、第１のダブルトーク状態はＤＢＬ＿１ｓｔにより表され、第２のダブルトーク状態はＤＢＬ＿２ｎｄにより示される。長時間レベルの減衰積分フィルタＬＴＬの出力は、検出されるダブルトークの状態に依存して変化する。従って、上記状態制御回路は、図５の長時間の減衰積分フィルタに接続されている。以下の式が当てはまると、第１のダブルトーク状態ＤＢＬ＿１ｓｔが検出される。
【数１２】
2007522720.tif 000042
以下の式が当てはまると、第２のダブルトーク状態ＤＢＬ＿２ｎｄが検出される。
【数１３】
2007522720.tif 000043
【００４９】
ダブルトークの状態制御は、図５のフィルタ係数更新部に接続されていない。２つのダブルトークの状態のうちの１つＤＢＬ＿１ｓｔ又はＤＢＬ＿２ｎｄが図５の状態制御回路により検出されると、上記フィルタ係数の更新は停止されない。ＤＢＬ＿ＴＨ１及びＤＢＬ＿ＴＨ２は定数であり、これらは第１及び第２のダブルトーク状態を検出する閾値を決定する。
【００５０】
図５には、収束速度を制御し、適応ブロッキング及び発散（divergence）を監視する他の制御ユニットも示されている。収束又は適応速度制御ブロックは、実際に達せされるエコーの減衰に依存して、適応フィルタアルゴリズムの適応速度に関する制御パラメータα（ｍ）を切り換える。制御パラメータα（ｍ）は、最適なフィルタ係数に近い収束の挙動を改善するためにプログラム可能な閾値において段階的に小さくされる。ダブルトーク状態ではなく、急激なシステムの変化が仮定されるべきであるので、タイムアウトの状態後における非常に長いが、実際には起こりそうにないフィルタ係数の更新の可能な長いブロッキングを伴う第２のダブルトークの指標の持続時間は回避されなければならない。適応フィルタがエコーを更に増幅する発散の検出の後、エコーキャンセラは全てのフィルタ係数をゼロにリセットすることにより再開される。
【００５１】
干渉の大きさは、図５に示されているＥＲＬ回路及びノルム選択回路により計算される。上記ＥＲＬ回路は、エコーリターンロス回路を意味している。上記ノルム選択回路の出力は、以下の式により与えられる。
【数１４】
2007522720.tif 000044
【００５２】
前述の式は、どちらの値がより大きいかに依存して、
2007522720.tif 000045
に等しい第１の値が選択されるか、又は重み付けされ、自乗された干渉の大きさに等しい第２の値が選択されるかのいずれか一方であることを示している。ノルム選択均等化（equation）回路は上記第１の値
2007522720.tif 000046
と前述の式の第２の値とを受け取り、異なる値のうちのどちらがフィルタ係数更新部に対して出力されるかを選択する。変数ｗ，ＥＲＬＴ（ｍ），Ｌ＿ＲＯ（ｋ）を有する前述の式の第２項は、ＥＲＬ回路により計算される。ｍは、標準的な時間の指標ｋ（サンプリングレート）に対するサブサンプリングされた時間の指標である。Ｌ＿ＲＯ（ｋ）はエコー消去信号の短時間レベルであり、ＥＲＬＴ（ｍ）は合計の達成されたエコーの減衰の推定値を表す長時間レベルＬＬ＿ＳＩ（ｍ）とＬＬ＿ＲＯ（ｍ）とのレベルの比であり、ｗは重み付け定数である。
ＥＲＬＴ（ｍ）＝ＬＬ＿ＳＩ（ｍ）／ＬＬ＿ＲＯ（ｍ）  （２２）
【００５３】
以下の式を用いて図５において上記フィルタの係数の更新が行われる。
【数１５】
2007522720.tif 000047
【００５４】
上記ＥＲＬ回路により計算される重み付けされ、自乗された干渉の大きさが
2007522720.tif 000048
よりも小さいと、係数の更新は標準的な正規化最小平均自乗アルゴリズム（ＮＬＭＳ）に従って行われる。その他の場合は、フィルタ係数が更新されるステップサイズ
2007522720.tif 000049
の制御が、ＥＲＬ回路から出力される重み付けされ、自乗された干渉の大きさにより行われる。
【００５５】
シミュレーションにより、ＥＲＬ回路により計算される重み付けされ、自乗された干渉の大きさが干渉の無い状態で基準信号のノルム
2007522720.tif 000050
に主に対応することが示された。従って、低い干渉のレベルの間の適応速度は、ダブルトークの無い状態でＮＬＭＳアルゴリズムの適応速度に対応する。
【００５６】
図５のエコーキャンセリング装置と最新技術のエコーキャンセリング装置との主な違いは、ダブルトーク状態の間の異なる制御に関連している。図５の状態制御回路により実行されるダブルトークの指示は、フィルタ係数の適応を完全には中断しない。それらは、長い項のレベルＬＬ＿ＳＩ（ｍ）及びＬＬ＿ＲＯ（ｍ）のロバスト更新を経由して合計のエコーロスの信頼性の高い推定を保証するので、間接的な作用のみを持つ。上記ＥＲＬ回路の出力は、より強力なやり方でますます多くの干渉信号を有し、これは、適応速度に対して重み付けされたノルムの選択を経由して自動的に作用する。少なくともシステムが実際に達せられる最適値から著しくドリフトする前に、干渉信号の強い上昇レベルがほとんどすぐに適応を中断する。
【００５７】
図６ａ，６ｂ，６ｃは、図３の既知のエコーキャンセリング装置と比較して図５のエコーキャンセリング装置の強化及びより良好な性能を明らかに示すシミュレーションの例を示している。図６ａは、基準信号ＳＩと図５及び図３の妨害している信号に対応する干渉信号とを示している。これらの信号は、図５及び図３のエコーキャンセラの性能を調べるために用いられる。上記性能は性能指標ＳＤＩＦ（ｍ）により決定される。ＳＤＩＦ（ｍ）は以下の式により計算される。
【数１６】
2007522720.tif 000051
ｈ（ｉ）は、エコー経路の正しいフィルタ係数を表している。
2007522720.tif 000052
は、ＦＩＲフィルタにより実際に用いられるフィルタ係数を表している。従って、ＳＤＩＦは、現在のフィルタ係数の自乗誤差の対数表示である。
【００５８】
図６ｂ及び図６ｃの破線は図３の従来のエコーキャンセリング装置の性能ＳＤＩＦを表し、図６ｂ及び図６ｃの実線は本発明により図５のエコーキャンセリング装置の性能ＳＤＩＦを表している。両方のグラフ６ｂ，６ｃは、条件付きの係数で始まる。これは、干渉信号が加えられる前のある量の時間の間に図５及び図３の回路が基準信号を用いてトレーニングされることを意味する。図６ｂでは、図５及び図３の装置におけるプリホワイトニングフィルタ及び逆プリホワイトニングフィルタＰＷＦ及びＩＰＷはディスエーブルにされ、図６ｃでは、図５及び図３の装置におけるプリホワイトニングフィルタ及び逆プリホワイトニングフィルタＰＷＦ及びＩＰＷはイネーブルにされている。両方の図は、本発明によるエコーキャンセリング装置の性能が妨害している干渉信号の開始により強くは劣化されないことを示している。
【図面の簡単な説明】
【００５９】
【図１】ラインエコーキャンセリング装置の従来の実現を示している。
【図２】音響エコーキャンセリング装置を示している。
【図３】最新技術による他のエコーキャンセラを示している。
【図４】本発明によるエコーキャンセリング装置の第１の実施の形態を示している。
【図５】本発明によるエコーキャンセリング装置の第２の実施の形態を示している。
【図６ａ】反射した基準信号およびダブルトーク状態を引き起こす干渉信号を示している。
【図６ｂ】従来のエコーキャンセリング装置の性能と比較して第２の実施の形態によるエコーキャンセリング装置の性能のグラフを示しており、第２の実施の形態のプリホワイトニングフィルタ及び逆プリホワイトニングフィルタがオフにされている。
【図６ｃ】従来のエコーキャンセリング装置の性能と比較して第２の実施の形態によるエコーキャンセリング装置の性能のグラフを示しており、第２の実施の形態のプリホワイトニングフィルタ及び逆プリホワイトニングフィルタがオンにされている。
【図１】
2007522720.tif 000053
【図２】
2007522720.tif 000054
【図３】
2007522720.tif 000055
【図４】
2007522720.tif 000056
【図５】
2007522720.tif 000057
【図６ａ】
2007522720.tif 000058
【図６ｂ】
2007522720.tif 000059
【図６ｃ】
2007522720.tif 000060
【国際調査報告】
2007522720.tif 000061
2007522720.tif 000062
2007522720.tif 000063