JP4504782B2 - エコー消去方法、この方法を実施する装置、プログラムおよびその記録媒体 - Google Patents

Description

この発明は、例えば多チャネル音響再生系を有する通信会議システムに適用され、ハウリングの原因及び聴覚上の障害となる音響エコーを消去する多チャネル反響消去方法、その装置、そのプログラム及びその記録媒体に関するものである。

近年のディジタルネットワークの大容量化により、複数の人が容易に参加でき、より自然な通話環境を提供できる多チャネル拡声型の通信会議システムが検討されている。このシステムでは、受話音声がスピーカから再生されマイクロホンに収音されて音響エコーが生じ、そのまま送信されると通話の障害や不快感などの問題が生じる。さらにスピーカ再生信号の信号パワーよりも音響エコー信号の信号パワーが大きくなる場合には、音響エコーがハウリングを引き起こして通話を不可能にしてしまう。快適な通話環境を実現するにはスピーカからマイクロホンへの音響的回り込みに対応する信号成分を、マイクロホン収音信号から消去する音響エコーキャンセラが必要となる。

Ｍ（≧２）チャネルの再生系とＮ（≧１）チャネルの収音系とで構成される通信会議システムは、図１に示すような構成により音響エコーの消去を行う。すなわち各受話端子１_ｍ（ｍ＝１，…，Ｍ）からの受話信号は各スピーカ２_ｍ（ｍ＝１，…，Ｍ）で音響信号として再生され、各Ｍ個の音響エコー経路を経て各マイクロホン３_ｎ（ｎ＝１，…，Ｎ）に回り込む。受話側の全Ｍチャネルの受話端子１_ｍ（ｍ＝１，…，Ｍ）と、Ｎチャネル送話側の送話端子５_ｎ（ｎ＝１，…，Ｎ）それぞれとの間にＭチャネルエコーキャンセル部６_ｎ（ｎ＝１，…，Ｎ）を接続して音響エコーを消去する。図１の構成は、図２のＭ入力１出力のエコーキャンセル部をＮ個並列に並べた構成になっている。

Ｍチャネルエコーキャンセル部６_ｎの内部では、多チャネル受話信号から予測エコー信号が生成され、マイクロホン収音信号との差ｅ（ｋ）および受話信号に基づいて収音信号と予測エコー信号の差が小さくなるように予測エコー生成用の適応フィルタ係数が逐次更新される。収音信号に非エコー信号が含まれる状況では、適応フィルタの係数がこの分だけ誤って更新されるので、適応フィルタが発散する危険もある。発散を回避するには適応フィルタの更新量を十分小さくする必要がある。実際には不必要に小さい更新量を選択するか、発散しない程度の大きさの更新量ではあるが非エコー信号の妨害による不正確な更新をある程度の確率で許容することになり、収束速度を低下させてしまう。

この非エコー成分に対して頑健に動作する多チャネルエコー消去法が特開２００３−２５０１９３号公報（特許文献１）に示されている。この方法と特開２００２−２２３１８２号公報（特許文献２）で示されている適応フィルタ処理での処理遅延が短いマルチディレイ・フィルタ（以下、ＭＤＦと略す。）法とを組み合わせた場合のＭチャネルエコーキャンセル部６の構成例を図３に示す。この方法では、残差信号すなわち収音信号と予測エコー信号との差信号に占めるエコー成分の比率を周波数成分ごとに推定し、この比率に応じて周波数成分ごとに適応フィルタ係数の更新を制御する。また、ＭＤＦ法は、畳み込み処理がより小さいブロック同士のオーバーラップセーブ処理に分割できることを利用する。適応フィルタのタップ長をＬ、分割数をＤ（ただしＬはＤで割り切れる）、Ｌ_Ｄ＝Ｌ／Ｄとすると、ＭＤＦ法ではＬ_Ｄサンプル毎に畳み込み処理が可能なため、Ｌ_Ｄサンプル毎に適応信号処理を適用することが可能になる。したがって、ハードウェア等の処理能力にかかわらずアルゴリズム上必ず存在する遅延である処理遅延を１／Ｄにできる。

予測エコー生成部６０では、各再生信号ｘ_ｍ（ｋ）（ｍ＝１，…，Ｍ）はＬ_Ｄサンプル毎に長さ２Ｌ_Ｄの信号ベクトルにフレーム化されて周波数領域に変換されて蓄積され、Ｄ個のフレームから周波数領域処理を経由して予測エコー信号が生成される。減算部６１では、マイクロホン３からの収音信号のフレームと予測エコー信号のフレームとの差がとられる。この残差信号フレームは、ＴＦ変換部６２を経てエコー成分比率推定部６３およびフィルタ係数更新部６４に入力される。エコー成分比率推定部６３は、周波数成分ごとに残差信号に含まれるエコー成分の比率を推定する。フィルタ係数更新部６４は、受話信号と残差信号から周波数成分ごとにフィルタ係数の修正方向を求め、その大きさをエコー成分の比率に応じて制御する。

エコー成分比率推定部６３の構成を図４に示す。相関除去部６３１では、多チャネル再生信号のスペクトルＸ_１（ｊ，ｆ），…，Ｘ_Ｍ（ｊ，ｆ）から互いに相関のない多チャネルのスペクトルＸ_１（ｊ，ｆ），Ｘ_２（１）（ｊ，ｆ），…，Ｘ_{Ｍ（Ｍ−１）}（ｊ，ｆ）を求める。相関除去部６３２では、残差信号のスペクトルＥ（ｊ，ｆ）から第１〜第ｍ−１チャネル再生信号の相関成分を除去したスペクトルＥ_{（ｍ―１）}（ｊ，ｆ）（ｍ＝２，…，Ｍ）を求める。コヒーレンス算出部６３３では、コヒーレンス算出部６３３_１で第１チャネル再生信号Ｘ_１（ｊ，ｆ）と残差信号Ｅ（ｊ，ｆ）のコヒーレンスγ_１ｙ ^２（ｊ，ｆ）を、コヒーレンス算出部６３３_ｍで短時間スペクトルＸ_{ｍ（ｍ−１）}（ｊ，ｆ）とＥ_{（ｍ−１）}（ｊ，ｆ）（ｍ＝２，…，Ｍ）のコヒーレンスγ_{ｍｙ（ｍ−１）} ^２（ｊ，ｆ）を求める。エコー成分比率算出部６３４では、次式によりエコー成分比率γ^２（ｊ，ｆ）を求める。

エコー成分比率推定のフローを図５に示す。
特開２００３−２５０１９３号公報 特開２００２−２２３１８２号公報

上記従来法では、適応フィルタのタップ長Ｌを残響時間にほぼ相当する長さに設定する。残響時間は、通常の部屋では３００ｍｓ以上になる。仮に、この時間にあわせてフレーム長２Ｌを取った場合、信号のバッファリングにより生じる遅延（処理遅延）がＬサンプルになり、通話系として非常に話しづらくなってしまう。そのため上記従来法では、この遅延を短くするために適応フィルタのタップ長Ｌよりもフレーム長２Ｌ_Ｄを短く設定し、エコー成分比率を推定し、適応フィルタ係数の更新量を制御する。
しかし、スピーカから再生されてマイクロホンに収音されるまでに２Ｌ_Ｄサンプル以上遅延するエコー成分は、非エコー成分として扱われることが問題となる。つまり、フレーム長２Ｌ_Ｄを、適応フィルタのタップ長Ｌよりも大幅に短く設定した場合には、非エコー成分とみなされるエコー成分が増大する。その結果、エコー成分比率が小さめに設定されたり、エコー成分の推定値が揺らいだりするために、エコー成分比率の推定性能が劣化し、エコー消去性能が劣化してしまう。

この発明では、残差信号の短時間スペクトルＥ（ｊ，ｆ）に含まれるエコー成分の比率を、現時点の多チャネル再生信号フレームから求めた短時間スペクトルＸ_１（ｊ，ｆ），…，Ｘ_Ｍ（ｊ，ｆ）だけでなく、過去の再生信号フレームから求めた短時間スペクトルも一緒に使用して推定する方法を提案する。
この発明では更に、多チャネル再生信号の現時点のフレームと過去のフレームとを、現時点のフレームの第１チャネル再生信号からなる主成分および主成分との相関が除去されたその他のフレームからなる副成分に分け、主成分のエコーが残差信号に占める割合を求め、副成分のエコーが主成分との相関を除去した残差信号に占める割合を求め、これら２つの割合から収音信号に占める多チャネル再生信号のエコー成分比率を推定する方法を提案する。

この方法により、過去の信号フレームをエコー成分比率の推定に取り込むことができ、フレーム長が残響時間よりも大幅に短く設定された場合でもエコー成分比率の推定性能劣化を回避し、エコー消去性能の劣化を防ぐことができる。

以下にこの発明の実施形態を図面を参照しながら説明するが、各図中の対応する部分は同一参照番号を付けて重複説明を省略する。
［第１実施形態］
この発明を図１のようなＭ（≧２）チャネル再生系とＮ（≧１）チャネル収音系からなる場合について説明する。上記従来法と同様に、図２のようなＭ入力１出力のエコーキャンセル部をＮ個並列に並べることで、Ｎチャネルの収音系に対応する。この発明では、図３に内部構成を示している図１のＭチャネルエコーキャンセル部６を、図６に内部構造が示されているＭチャネルエコーキャンセル部７に置き換える。図７には図６中のエコー成分比率推定部７３の構成を示す。

以下では、適応フィルタのタップ長をＬ、分割数をＤ（ただしＬはＤで割り切れる）、Ｌ_Ｄ＝Ｌ／Ｄ、フレーム長を２Ｌ_Ｄサンプル、シフト長をＬ_Ｄサンプル、フレーム時刻をｊとする。フレーム時刻ｊの信号フレームは、サンプル時刻ｋ＝ｊＬ_Ｄ−２Ｌ_Ｄ＋１〜ｊＬ_Ｄの信号サンプルからなる。このときの信号のサンプル時刻ｋとフレーム時刻ｊの関係は図１２のようになる。また、過去の再生信号フレームから求めたスペクトルとして、１フレーム前の短時間スペクトルＸ_１（ｊ−１，ｆ），…，Ｘ_Ｍ（ｊ−１，ｆ）を使用する例を説明する。

図６の予測エコー生成部７０内のＴＦ変換部７０１_ｍ（ｍ＝１，…，Ｍ）において、各チャネルの時間領域の再生信号ｘ_ｍ（ｋ）をＬ_Ｄサンプル毎に長さ２Ｌ_Ｄの信号ベクトルにフレーム化し、ＦＦＴを使って短時間スペクトルに変換し、フィルタ処理部７０２に蓄積する。

次にフィルタ処理部７０２において、次式に示すように蓄積された第ｍチャネル再生信号の短時間スペクトルＸ_ｍ（ｊ，ｆ）と周波数領域の適応フィルタ係数Ｗ_ｍ（ｊ，ｄ，ｆ）を乗算処理し、足し合わせることで、第ｍスピーカからマイクロホンに回り込む周波数領域の予測エコー信号Ｙ＾_ｍ（ｊ，ｆ）を得る。

この各スピーカからの予測エコー信号Ｙ＾_ｍ（ｊ，ｆ）を足し合わせて周波数領域の予測エコー信号Ｙ＾（ｊ，ｆ）を次式のように得る。

ＦＴ変換部７０３において、周波数領域の予測エコーを逆ＦＦＴにより時間領域に変換し、時間領域の予測エコーｙ＾（ｋ）を次のように得る。

減算部７１では、フレーム化部７１１において時間領域の収音信号ｙ（ｋ）をフレーム化し、減算器７１２において予測エコー信号ｙ＾（ｋ）のフレームを差し引いて時間領域での残差信号ｅ（ｋ）のフレームを求める。

ＴＦ変換部７２では、次式のように時間領域の残差信号ｅ（ｋ）のフレームにＬ_Ｄ個の０を詰めて周波数領域に変換することで周波数領域の残差信号Ｅ（ｊ，ｆ）を求める。

エコー成分比率推定部７３において以下のステップＦ１〜Ｆ７により、周波数領域の多チャネル再生信号Ｘ_ｍ（ｊ，ｆ）と周波数領域の残差信号Ｅ（ｊ，ｆ）から、周波数成分ごとに残差信号に含まれるエコー成分の比率Ｓ（ｊ，ｆ）を求める。図８にエコー成分比率Ｓ（ｊ，ｆ）を推定するためのフローを示す。
ステップＦ１
現時点のフレームから求めた多チャネル再生信号の短時間スペクトルＸ_１（ｊ，ｆ），…，Ｘ_Ｍ（ｊ，ｆ）を図７に示す相関除去部７３１内の現時点の蓄積部７３１ａ１に保存する。

ステップＦ２
相関除去部７３１ｂ１では、例えば次式の方法で多チャネル再生信号の短時間スペクトルＸ_２（ｊ，ｆ），…，Ｘ_Ｍ（ｊ，ｆ）からＸ_１（ｊ，ｆ）との相関成分を除去して、スペクトルＸ_２（１）（ｊ，ｆ），…，Ｘ_Ｍ（１）（ｊ，ｆ）を得、多チャネル再生信号スペクトルの副成分の一部とする。

ここで、ε［］は、平均をとることを意味し、平均処理の一例としては、

のように、１フレーム前の処理結果と０〜１の値をとる平滑化定数βを用いる方法がある。

ステップＦ３
相関除去部７３１ｂ２において、過去の蓄積部７３１ａ２に蓄積された１フレーム前の多チャネル再生信号のスペクトルＸ_１（ｊ−１，ｆ），…，Ｘ_Ｍ（ｊ−１，ｆ）から、Ｘ_１（ｊ，ｆ）との相関を次のように除去したスペクトルＸ_１（１）（ｊ−１，ｆ），…，Ｘ_Ｍ（１）（ｊ−１，ｆ）を求め、多チャネル再生信号スペクトルの副成分の一部とする。

なお、ｎフレーム前の短時間スペクトルＸ_１（ｊ−ｎ，ｆ），…，Ｘ_Ｍ（ｊ−ｎ，ｆ）をエコー成分比率推定に使用する場合にも、同様の計算により得られた結果を多チャネル再生信号スペクトルの副成分の一部とすればよい。

ステップＦ４
相関除去部７３２では、現時点のフレームの残差信号の短時間スペクトルＥ（ｊ，ｆ）からＸ_１（ｊ，ｆ）との相関成分を除去したスペクトルＥ_（１）（ｊ，ｆ）を求める。

ステップＦ５
コヒーレンス算出部７３３_１では、多チャネル再生信号スペクトルの主成分である現時点のフレームの第１チャネル再生信号の短時間スペクトルＸ_１（ｊ，ｆ）と現時点の残差信号のスペクトルＥ（ｊ，ｆ）から、次のコヒーレンスを求める。

ステップＦ６
副成分エコー比率算出部７３３_２では、まず相関除去された残差信号スペクトルＥ_（１）（ｊ，ｆ）に含まれるエコー成分Ｅ＾_（１）（ｊ，ｆ）を求める。エコー成分Ｅ＾_（１）（ｊ，ｆ）は、多チャネル再生信号短時間スペクトルの副成分Ｘ_２（１）（ｊ，ｆ），…，Ｘ_Ｍ（１）（ｊ，ｆ），Ｘ_１（１）（ｊ−１，ｆ），…，Ｘ_Ｍ（１）（ｊ−１，ｆ）の線形和

のうちで、相関除去された残差信号スペクトルとの誤差
｜Ｅ_（１）（ｊ，ｆ）−Ｅ＾_（１）（ｊ，ｆ）｜^２
が最小となるスペクトルである。この誤差を最小にするスペクトルは、

とし、

により求められる。さらに、次式により副成分のエコー比率を周波数成分ごとに求める。

ステップＦ７
図７のエコー成分比率算出部７３４において、ステップＦ５、Ｆ６で求めた各比率から、残差信号スペクトルＥ（ｊ，ｆ）に占めるエコー成分の比率を求める。

図６のフィルタ更新部７４では、まず周波数成分ごとに残差信号Ｅ（ｊ，ｆ）と再生信号Ｘ_ｍ（ｊ，ｆ）から、周波数領域の適応フィルタ係数の修正量ｄＷ_ｍ（ｊ，ｄ，ｆ）を求める。次に周波数領域の残差信号の各周波数成分に含まれる残留エコー成分の比率Ｓ（ｊ，ｆ）を使用して適応フィルタ係数の修正量ｄＷ_ｍ（ｊ，ｄ，ｆ）の大きさを制御しながら、周波数領域の適応フィルタ係数Ｗ_ｍ（ｊ，ｄ，ｆ）を逐次更新する。具体的には、例えばJia-sien Soo, Khee K. Pang, “Multidelay Block Frequency Domain Adaptive Filter,” IEEE Transactions on Acoustics, Speech, and Signal Processing, Vol.38, No.2, pp.373-376, 1990.の提案と組み合わせ、以下のように修正量ｄＷ_ｍ（ｊ，ｄ，ｆ）および適応フィルタ係数Ｗ_ｍ（ｊ，ｄ，ｆ）を求めることができる。

ただし、ｄ＝１，…，Ｄ、Ｘ_ｍ ^＊（ｊ，ｆ）はＸ_ｍ（ｊ，ｆ）の複素共役である。

ここでＶ _ｍ （ｊ＋1，ｄ，ｆ）は、既に計算済みの適応フィルタ係数Ｗ _ｍ （ｊ，ｄ，ｆ）及び修正量ｄＷ _ｍ （ｊ，ｄ，ｆ）から以下の式で計算される。

ただし、ｍ＝１，…，Ｍ、ｄ＝１，…，Ｄ、ｆ＝１，…，２Ｌ_Ｄ、μは０〜１の値をとるステップサイズ、δは正則化定数である。また、ｐ（ｊ，ｆ）は次式により計算される。

このように適応フィルタ係数Ｗ_ｍ（ｊ，ｄ，ｆ）を処理単位ごとに更新することで、次の処理では更新された適応フィルタ係数が使用できる。
現時点のフレームの処理が終了すると、最後に現時点の蓄積部７３１ａ１に蓄積された再生信号情報は過去の蓄積部７３１ａ２に転送され、蓄積される。
なお、蓄積部７３１ａ内で現時点の蓄積部７３１ａ１と過去の蓄積部７３１ａ２とを特に区別し、上記のように一連の処理の最後に現時点の蓄積部７３１ａ１に蓄積された再生信号情報を過去の蓄積部７３１ａ２に転送するのではなく、１つの蓄積部７３１ａに蓄積された情報の中で最新情報を現時点の情報として処理する方法もある。また、図９に示すように処理に利用する現時点の再生信号のスペクトルを、蓄積部から取り出すのではなく、入力された再生信号のスペクトルを直接利用する方法もある。

［第２実施形態］
この発明は、第１実施形態に対してエコー経路推定速度の向上を図ったものであり、図１０に全体構成を、図１１にＭチャネルエコーキャンセル部の構成を示す。
図１０の相関変動処理部９_ｍ（ｍ＝１，…，Ｍ）では、例えば次式のようにチャネル毎に受話信号ｕ_ｍ（ｋ）を乱数で振幅変調して元の信号に付加して相互相関が絶えず変動している再生信号ｘ_ｍ（ｋ）を生成し、各スピーカから再生する。

この相関変動処理がない場合、Ｍチャネルの受話信号のチャネル間相関が高いために、第ｍチャネル再生信号と第１チャネルの収音信号の入出力関係から、第ｍチャネルの真のエコー経路特性が必ず求められるとは限らない。しかし、この相関変動処理を受話信号に適用することにより、Ｍチャネルのエコー経路特性が確実に推定できる。
図１１のＭチャネルエコーキャンセル部８では、付加信号成分が受話信号成分よりも強調されたチャネル間相関の小さい信号のベクトルを信号変換部８５により生成し、チャネル間相関の大きい再生信号ベクトルの代わりに使用して適応フィルタを更新する。この更新には、例えば特許文献２で提案されている適応アルゴリズムを使用し、そのステップサイズを残差信号に含まれる残留エコー成分の比率Ｓ（ｊ，ｆ）で制御する。チャネル間相関の小さい修正ベクトルを使用することで、適応フィルタによるエコー経路の収束速度が向上する。

Ｍ入力Ｎ出力のエコー消去装置の一般的構成を示す図。 Ｍ入力１出力のエコー消去装置の一般的構成を示す図。 従来のＭチャネルエコーキャンセル部の構成を示す図。 従来のエコー成分比率推定部の構成を示す図。 従来のエコー成分比率推定のフローを示す図。 第１実施形態のＭチャネルエコーキャンセル部の構成を示す図。 第１実施形態のエコー成分比率推定部の構成を示す図。 第１実施形態のエコー成分比率推定のフローを示す図。 第１実施形態のエコー成分比率推定部の変形例の構成を示す図。 第２実施形態のＭ入力１出力チャネルのエコー消去装置の一般的構成を示す図。 第２実施形態のＭチャネルエコーキャンセル部の構成を示す図。 信号のサンプル時刻ｋとフレーム時刻ｊの関係を示す図。

Claims (14)

1. 複数チャネル（Ｍチャネル）の再生信号から適応フィルタを用いてエコー信号を予測し、少なくとも１チャネルの収音信号の各収音信号から予測したエコー信号を差し引くことでエコー成分を消去する方法において、
   現時点のフレームの第１チャネル再生信号の短時間スペクトルを主成分とし、
   現時点のフレームの第２から第Ｍチャネルまでの再生信号、及び、少なくとも１フレーム過去の第１から第Ｍチャネルまでの再生信号について、各々の短時間スペクトルから主成分とした短時間スペクトルとの相関を除去して、副成分を構成する複数の短時間スペクトルを求め、
   主成分のエコーが残差信号の短時間スペクトルに占める割合を求め、
   副成分のエコーが主成分との相関を除去した残差信号の短時間スペクトルに占める割合を求め、
   上記２つの割合から残差信号の短時間スペクトルに占めるエコー成分比率を周波数ごとに推定し、
   周波数ごとに推定された上記エコー成分比率が大きいほど、適応フィルタ係数を修正する量が大きくなるように、対応する周波数の適応フィルタ係数の値を更新する、
   ことを特徴とするエコー消去方法。
2. 請求項１記載の方法において、
   残差信号の短時間スペクトルに占めるエコー成分比率Ｓ（ｆ）を、主成分のエコーが残差信号の短時間スペクトルに占める割合γ_１ ^２（ｆ）と副成分のエコーが主成分との相関を除去した残差信号の短時間スペクトルに占める割合γ_２ ^２（ｆ）から、
   

   

   により求めること、
   を特徴とするエコー消去方法。
3. 請求項２記載の方法において、
   主成分との相関が除去された残差信号の短時間スペクトルＥ_（１）（ｆ）に含まれるエコー成分Ｅ＾_（１）（ｆ）を、｜Ｅ_（１）（ｆ）−Ｅ＾_（１）（ｆ）｜^２を最小とする線形和として求め、
   副成分のエコーが主成分との相関を除去した残差信号の短時間スペクトルに占める割合γ_２ ^２（ｆ）を、
   

   

   により求めること、
   を特徴とするエコー消去方法。
4. 請求項１〜３の何れかに記載の方法において、
   適応フィルタのタップ長よりも短いフレーム長ごとに残差信号の短時間スペクトルに含まれるエコー成分の比率を推定すること、
   を特徴とするエコー消去方法。
5. 請求項１〜４の何れかに記載の方法において、
   適応フィルタ係数をフレームごとに修正すること、
   を特徴とするエコー消去方法。
6. 請求項１〜５の何れかに記載の方法において、
   再生信号として、受話信号に相関変動処理した受話信号を付加した信号を用い、
   適応フィルタ更新に付加信号成分が受話信号成分よりも強調された信号を使用すること、
   を特徴とするエコー消去方法。
7. 複数チャネルの再生信号および少なくとも１チャネルの収音信号を受信する手段と、
   エコー信号を予測する適応フィルタと、
   予測したエコー信号を各収音部からの信号から差し引く手段と、
   現時点のフレームの第１チャネル再生信号の短時間スペクトル（主成分）のエコーが、残差信号の短時間スペクトルに占める割合を求める手段と、
   現時点のフレームの第２から第Ｍチャネルまでの再生信号、及び、少なくとも１フレーム過去の第１から第Ｍチャネルまでの再生信号について、各々の短時間スペクトルから主成分とした短時間スペクトルとの相関を除去して、副成分を構成する複数の短時間スペクトルを求める手段と、
   副成分のエコーが主成分との相関を除去した残差信号の短時間スペクトルに占める割合を求める手段と、
   上記２つの割合から残差信号の短時間スペクトルに占めるエコー成分比率を周波数ごとに推定する手段と、
   周波数ごとに推定された上記エコー成分比率が大きいほど、適応フィルタ係数を修正する量が大きくなるように、対応する周波数の適応フィルタ係数の値を更新する手段と、
   を備えるエコー消去装置。
8. 請求項７記載の装置において、
   残差信号の短時間スペクトルに占めるエコー成分比率Ｓ（ｆ）を推定する手段として、主成分のエコーが残差信号の短時間スペクトルに占める割合γ_１ ^２（ｆ）と副成分のエコーが主成分との相関を除去した残差信号の短時間スペクトルに占める割合γ_２ ^２（ｆ）から、
   

   

   により求める手段、
   を備えるエコー消去装置。
9. 請求項８記載の装置において、
   副成分のエコーが主成分との相関を除去した残差信号の短時間スペクトルに占める割合γ_２ ^２（ｆ）を求める手段として、主成分との相関が除去された残差信号の短時間スペクトルＥ_（１）（ｆ）に含まれるエコー成分Ｅ＾_（１）（ｆ）を、
   ｜Ｅ_（１）（ｆ）−Ｅ＾_（１）（ｆ）｜^２
   を最小とする線形和として求め、副成分のエコーが主成分との相関を除去した残差信号の短時間スペクトルに占める割合γ_２ ^２（ｆ）を、
   

   

   により求める手段、
   を備えるエコー消去装置。
10. 請求項７〜９の何れかに記載の装置において、
    適応フィルタのタップ長よりも短いフレーム長ごとに残差信号の短時間スペクトルに含まれるエコー成分の比率を推定する手段、
    を備えるエコー消去装置。
11. 請求項７〜１０の何れかに記載の装置において、
    適応フィルタ係数をフレームごとに修正する手段、
    を備えるエコー消去装置。
12. 請求項７〜１１の何れかに記載の装置において、
    受話信号の相関を変動処理する相関変動処理部と、
    付加信号成分が受話信号成分よりも強調された信号を生成する信号変換部と、
    を備えるエコー消去装置。
13. 請求項１〜６の何れかに記載のエコー消去方法をコンピュータにより実行するエコー消去プログラム。
14. 請求項１３に記載のエコー消去プログラムを記録したコンピュータ読み取り可能な記録媒体。

Images