(57)【要約】【課題】エコー消去方法において、エコー成分比率推定性能の劣化を避けながら、アルゴリズム上必ず生じる処理遅延を大幅に短くすることを可能にする。【解決手段】残差信号の短時間スペクトルに含まれるエコー成分の比率を、現時点のフレームの多チャネル再生信号の短時間スペクトルだけでなく、過去のフレームの再生信号の短時間スペクトルも一緒に使用して推定する。具体的には、多チャネル再生信号の現時点のフレームと過去のフレームとを、現時点のフレームの第１チャネル再生信号の短時間スペクトルからなる主成分および主成分との相関が除去されたその他の再生信号の短時間スペクトルからなる副成分に分け、主成分のエコーが残差信号に占める割合を求め、副成分のエコーが主成分との相関を除去した残差信号に占める割合を求め、これら２つの割合から収音信号に占める多チャネル再生信号のエコー成分比率を推定する。【選択図】図７
2006121589.tif 000002
【特許請求の範囲】
【請求項１】
複数チャネルの再生信号からエコー信号を予測し、少なくとも１チャネルの収音信号の各収音信号から予測したエコー信号を差し引くことでエコー成分を消去する方法において、
残差信号の短時間スペクトルに含まれるエコー成分の比率を、現時点のフレームの再生信号の短時間スペクトルおよび少なくとも１つの過去のフレームの再生信号の短時間スペクトルから、周波数ごとに推定し、
上記残差信号の短時間スペクトルに含まれるエコー成分の比率から、適応フィルタ係数を修正すること、
を特徴とするエコー消去方法。
【請求項２】
請求項１記載の方法において、
適応フィルタのタップ長よりも短いフレーム長ごとに残差信号の短時間スペクトルに含まれるエコー成分の比率を推定すること、
を特徴とするエコー消去方法。
【請求項３】
請求項１又は２記載の方法において、
適応フィルタ係数をフレームごとに修正すること、
を特徴とするエコー消去方法。
【請求項４】
請求項１〜３の何れかに記載の方法において、
現時点のフレームの第１チャネル再生信号の短時間スペクトルを主成分とし、
主成分との相関を除去したその他の短時間スペクトルを副成分とし、
主成分のエコーが残差信号の短時間スペクトルに占める割合を求め、
副成分のエコーが主成分との相関を除去した残差信号の短時間スペクトルに占める割合を求め、
上記２つの割合から残差信号の短時間スペクトルに占めるエコー成分比率を推定すること、
を特徴とするエコー消去方法。
【請求項５】
請求項４記載の方法において、
残差信号の短時間スペクトルに占めるエコー成分比率Ｓ（ｆ）を、主成分のエコーが残差信号の短時間スペクトルに占める割合γ１２（ｆ）と副成分のエコーが主成分との相関を除去した残差信号の短時間スペクトルに占める割合γ２２（ｆ）から、
【数１】
2006121589.tif 000024
により求めること、
を特徴とするエコー消去方法。
【請求項６】
請求項５記載の方法において、
主成分との相関が除去された残差信号の短時間スペクトルＥ（１）（ｆ）に含まれるエコー成分Ｅ＾（１）（ｆ）を、｜Ｅ（１）（ｆ）−Ｅ＾（１）（ｆ）｜２を最小とする線形和として求め、
副成分のエコーが主成分との相関を除去した残差信号の短時間スペクトルに占める割合γ２２（ｆ）を、
【数２】
2006121589.tif 000025
により求めること、
を特徴とするエコー消去方法。
【請求項７】
請求項１〜６の何れかに記載の方法において、
再生信号として、受話信号に相関変動処理した受話信号を付加した信号を用い、
適応フィルタ更新に付加信号成分が受話信号成分よりも強調された信号を使用すること、
を特徴とするエコー消去方法。
【請求項８】
複数チャネルの再生信号および少なくとも１チャネルの収音信号を受信する手段と、
エコー信号を予測する適応フィルタと、
予測したエコー信号を各収音部からの信号から差し引く手段と、
残差信号の短時間スペクトルに含まれるエコー成分の比率を、現時点のフレームの再生信号の短時間スペクトルおよび少なくとも1つの過去のフレームの再生信号の短時間スペクトルから、周波数ごとに推定する手段と、
上記残差信号の短時間スペクトルに含まれるエコー成分の比率から、適応フィルタ係数を修正する手段と、
を備えるエコー消去装置。
【請求項９】
請求項９記載の装置において、
適応フィルタのタップ長よりも短いフレーム長ごとに残差信号の短時間スペクトルに含まれるエコー成分の比率を推定する手段、
を備えるエコー消去装置。
【請求項１０】
請求項８又は９記載の装置において、
適応フィルタ係数をフレームごとに修正する手段、
を備えるエコー消去装置。
【請求項１１】
請求項８〜１０の何れかに記載の装置において、
現時点のフレームの第１チャネル再生信号の短時間スペクトル（主成分）のエコーが、残差信号の短時間スペクトルに占める割合を求める手段と、
主成分との相関を除去したその他の短時間スペクトル（副成分）のエコーが、主成分との相関を除去した残差信号の短時間スペクトルに占める割合を求める手段と、
上記２つの割合を求める手段の出力から残差信号の短時間スペクトルに占めるエコー成分比率を推定する手段と、
を備えるエコー消去装置。
【請求項１２】
請求項１１記載の装置において、
残差信号の短時間スペクトルに占めるエコー成分比率Ｓ（ｆ）を推定する手段として、主成分のエコーが残差信号の短時間スペクトルに占める割合γ１２（ｆ）と副成分のエコーが主成分との相関を除去した残差信号の短時間スペクトルに占める割合γ２２（ｆ）から、
【数３】
2006121589.tif 000026
により求める手段、
を備えるエコー消去装置。
【請求項１３】
請求項１２記載の装置において、
副成分のエコーが主成分との相関を除去した残差信号の短時間スペクトルに占める割合γ２２（ｆ）を求める手段として、主成分との相関が除去された残差信号の短時間スペクトルＥ（１）（ｆ）に含まれるエコー成分Ｅ＾（１）（ｆ）を、
｜Ｅ（１）（ｆ）−Ｅ＾（１）（ｆ）｜２
を最小とする線形和として求め、副成分のエコーが主成分との相関を除去した残差信号の短時間スペクトルに占める割合γ２２（ｆ）を、
【数４】
2006121589.tif 000027
により求める手段、
を備えるエコー消去装置。
【請求項１４】
請求項８〜１３の何れかに記載の装置において、
受話信号の相関を変動処理する相関変動処理部と、
付加信号成分が受話信号成分よりも強調された信号を生成する信号変換部と、
を備えるエコー消去装置。
【請求項１５】
請求項１〜７に記載のエコー消去方法をコンピュータにより実行するエコー消去プログラム。
【請求項１６】
請求項１５に記載のエコー消去プログラムを記録したコンピュータ読み取り可能な記録媒体。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、例えば多チャネル音響再生系を有する通信会議システムに適用され、ハウリングの原因及び聴覚上の障害となる音響エコーを消去する多チャネル反響消去方法、その装置、そのプログラム及びその記録媒体に関するものである。  
【背景技術】
【０００２】
近年のディジタルネットワークの大容量化により、複数の人が容易に参加でき、より自然な通話環境を提供できる多チャネル拡声型の通信会議システムが検討されている。このシステムでは、受話音声がスピーカから再生されマイクロホンに収音されて音響エコーが生じ、そのまま送信されると通話の障害や不快感などの問題が生じる。さらにスピーカ再生信号の信号パワーよりも音響エコー信号の信号パワーが大きくなる場合には、音響エコーがハウリングを引き起こして通話を不可能にしてしまう。快適な通話環境を実現するにはスピーカからマイクロホンへの音響的回り込みに対応する信号成分を、マイクロホン収音信号から消去する音響エコーキャンセラが必要となる。
【０００３】
Ｍ（≧２）チャネルの再生系とＮ（≧１）チャネルの収音系とで構成される通信会議システムは、図１に示すような構成により音響エコーの消去を行う。すなわち各受話端子１ｍ（ｍ＝１，…，Ｍ）からの受話信号は各スピーカ２ｍ（ｍ＝１，…，Ｍ）で音響信号として再生され、各Ｍ個の音響エコー経路を経て各マイクロホン３ｎ（ｎ＝１，…，Ｎ）に回り込む。受話側の全Ｍチャネルの受話端子１ｍ（ｍ＝１，…，Ｍ）と、Ｎチャネル送話側の送話端子５ｎ（ｎ＝１，…，Ｎ）それぞれとの間にＭチャネルエコーキャンセル部６ｎ（ｎ＝１，…，Ｎ）を接続して音響エコーを消去する。図１の構成は、図２のＭ入力１出力のエコーキャンセル部をＮ個並列に並べた構成になっている。
【０００４】
Ｍチャネルエコーキャンセル部６ｎの内部では、多チャネル受話信号から予測エコー信号が生成され、マイクロホン収音信号との差ｅ（ｋ）および受話信号に基づいて収音信号と予測エコー信号の差が小さくなるように予測エコー生成用の適応フィルタ係数が逐次更新される。収音信号に非エコー信号が含まれる状況では、適応フィルタの係数がこの分だけ誤って更新されるので、適応フィルタが発散する危険もある。発散を回避するには適応フィルタの更新量を十分小さくする必要がある。実際には不必要に小さい更新量を選択するか、発散しない程度の大きさの更新量ではあるが非エコー信号の妨害による不正確な更新をある程度の確率で許容することになり、収束速度を低下させてしまう。
【０００５】
この非エコー成分に対して頑健に動作する多チャネルエコー消去法が特開２００３−２５０１９３号公報（特許文献１）に示されている。この方法と特開２００２−２２３１８２号公報（特許文献２）で示されている適応フィルタ処理での処理遅延が短いマルチディレイ・フィルタ（以下、ＭＤＦと略す。）法とを組み合わせた場合のＭチャネルエコーキャンセル部６の構成例を図３に示す。この方法では、残差信号すなわち収音信号と予測エコー信号との差信号に占めるエコー成分の比率を周波数成分ごとに推定し、この比率に応じて周波数成分ごとに適応フィルタ係数の更新を制御する。また、ＭＤＦ法は、畳み込み処理がより小さいブロック同士のオーバーラップセーブ処理に分割できることを利用する。適応フィルタのタップ長をＬ、分割数をＤ（ただしＬはＤで割り切れる）、ＬＤ＝Ｌ／Ｄとすると、ＭＤＦ法ではＬＤサンプル毎に畳み込み処理が可能なため、ＬＤサンプル毎に適応信号処理を適用することが可能になる。したがって、ハードウェア等の処理能力にかかわらずアルゴリズム上必ず存在する遅延である処理遅延を１／Ｄにできる。
【０００６】
予測エコー生成部６０では、各再生信号ｘｍ（ｋ）（ｍ＝１，…，Ｍ）はＬＤサンプル毎に長さ２ＬＤの信号ベクトルにフレーム化されて周波数領域に変換されて蓄積され、Ｄ個のフレームから周波数領域処理を経由して予測エコー信号が生成される。減算部６１では、マイクロホン３からの収音信号のフレームと予測エコー信号のフレームとの差がとられる。この残差信号フレームは、ＴＦ変換部６２を経てエコー成分比率推定部６３およびフィルタ係数更新部６４に入力される。エコー成分比率推定部６３は、周波数成分ごとに残差信号に含まれるエコー成分の比率を推定する。フィルタ係数更新部６４は、受話信号と残差信号から周波数成分ごとにフィルタ係数の修正方向を求め、その大きさをエコー成分の比率に応じて制御する。
【０００７】
エコー成分比率推定部６３の構成を図４に示す。相関除去部６３１では、多チャネル再生信号のスペクトルＸ１（ｊ，ｆ），…，ＸＭ（ｊ，ｆ）から互いに相関のない多チャネルのスペクトルＸ１（ｊ，ｆ），Ｘ２（１）（ｊ，ｆ），…，ＸＭ（Ｍ−１）（ｊ，ｆ）を求める。相関除去部６３２では、残差信号のスペクトルＥ（ｊ，ｆ）から第１〜第ｍ−１チャネル再生信号の相関成分を除去したスペクトルＥ（ｍ―１）（ｊ，ｆ）（ｍ＝２，…，Ｍ）を求める。コヒーレンス算出部６３３では、コヒーレンス算出部６３３１で第１チャネル再生信号Ｘ１（ｊ，ｆ）と残差信号Ｅ（ｊ，ｆ）のコヒーレンスγ１ｙ２（ｊ，ｆ）を、コヒーレンス算出部６３３ｍで短時間スペクトルＸｍ（ｍ−１）（ｊ，ｆ）とＥ（ｍ−１）（ｊ，ｆ）（ｍ＝２，…，Ｍ）のコヒーレンスγｍｙ（ｍ−１）２（ｊ，ｆ）を求める。エコー成分比率算出部６３４では、次式によりエコー成分比率γ２（ｊ，ｆ）を求める。
【数５】
2006121589.tif 000003
エコー成分比率推定のフローを図５に示す。
【特許文献１】特開２００３−２５０１９３号公報
【特許文献２】特開２００２−２２３１８２号公報
【発明の開示】
【発明が解決しようとする課題】
【０００８】
上記従来法では、適応フィルタのタップ長Ｌを残響時間にほぼ相当する長さに設定する。残響時間は、通常の部屋では３００ｍｓ以上になる。仮に、この時間にあわせてフレーム長２Ｌを取った場合、信号のバッファリングにより生じる遅延（処理遅延）がＬサンプルになり、通話系として非常に話しづらくなってしまう。そのため上記従来法では、この遅延を短くするために適応フィルタのタップ長Ｌよりもフレーム長２ＬＤを短く設定し、エコー成分比率を推定し、適応フィルタ係数の更新量を制御する。
しかし、スピーカから再生されてマイクロホンに収音されるまでに２ＬＤサンプル以上遅延するエコー成分は、非エコー成分として扱われることが問題となる。つまり、フレーム長２ＬＤを、適応フィルタのタップ長Ｌよりも大幅に短く設定した場合には、非エコー成分とみなされるエコー成分が増大する。その結果、エコー成分比率が小さめに設定されたり、エコー成分の推定値が揺らいだりするために、エコー成分比率の推定性能が劣化し、エコー消去性能が劣化してしまう。
【課題を解決するための手段】
【０００９】
この発明では、残差信号の短時間スペクトルＥ（ｊ，ｆ）に含まれるエコー成分の比率を、現時点の多チャネル再生信号フレームから求めた短時間スペクトルＸ１（ｊ，ｆ），…，ＸＭ（ｊ，ｆ）だけでなく、過去の再生信号フレームから求めた短時間スペクトルも一緒に使用して推定する方法を提案する。
この発明では更に、多チャネル再生信号の現時点のフレームと過去のフレームとを、現時点のフレームの第１チャネル再生信号からなる主成分および主成分との相関が除去されたその他のフレームからなる副成分に分け、主成分のエコーが残差信号に占める割合を求め、副成分のエコーが主成分との相関を除去した残差信号に占める割合を求め、これら２つの割合から収音信号に占める多チャネル再生信号のエコー成分比率を推定する方法を提案する。
【発明の効果】
【００１０】
この方法により、過去の信号フレームをエコー成分比率の推定に取り込むことができ、フレーム長が残響時間よりも大幅に短く設定された場合でもエコー成分比率の推定性能劣化を回避し、エコー消去性能の劣化を防ぐことができる。
【発明を実施するための最良の形態】
【００１１】
以下にこの発明の実施形態を図面を参照しながら説明するが、各図中の対応する部分は同一参照番号を付けて重複説明を省略する。
［第１実施形態］
この発明を図１のようなＭ（≧２）チャネル再生系とＮ（≧１）チャネル収音系からなる場合について説明する。上記従来法と同様に、図２のようなＭ入力１出力のエコーキャンセル部をＮ個並列に並べることで、Ｎチャネルの収音系に対応する。この発明では、図３に内部構成を示している図１のＭチャネルエコーキャンセル部６を、図６に内部構造が示されているＭチャネルエコーキャンセル部７に置き換える。図７には図６中のエコー成分比率推定部７３の構成を示す。
【００１２】
以下では、適応フィルタのタップ長をＬ、分割数をＤ（ただしＬはＤで割り切れる）、ＬＤ＝Ｌ／Ｄ、フレーム長を２ＬＤサンプル、シフト長をＬＤサンプル、フレーム時刻をｊとする。フレーム時刻ｊの信号フレームは、サンプル時刻ｋ＝ｊＬＤ−２ＬＤ＋１〜ｊＬＤの信号サンプルからなる。このときの信号のサンプル時刻ｋとフレーム時刻ｊの関係は図１２のようになる。また、過去の再生信号フレームから求めたスペクトルとして、１フレーム前の短時間スペクトルＸ１（ｊ−１，ｆ），…，ＸＭ（ｊ−１，ｆ）を使用する例を説明する。
【００１３】
図６の予測エコー生成部７０内のＴＦ変換部７０１ｍ（ｍ＝１，…，Ｍ）において、各チャネルの時間領域の再生信号ｘｍ（ｋ）をＬＤサンプル毎に長さ２ＬＤの信号ベクトルにフレーム化し、ＦＦＴを使って短時間スペクトルに変換し、フィルタ処理部７０２に蓄積する。
【数６】
2006121589.tif 000004
【００１４】
次にフィルタ処理部７０２において、次式に示すように蓄積された第ｍチャネル再生信号の短時間スペクトルＸｍ（ｊ，ｆ）と周波数領域の適応フィルタ係数Ｗｍ（ｊ，ｄ，ｆ）を乗算処理し、足し合わせることで、第ｍスピーカからマイクロホンに回り込む周波数領域の予測エコー信号Ｙ＾ｍ（ｊ，ｆ）を得る。
【数７】
2006121589.tif 000005
この各スピーカからの予測エコー信号Ｙ＾ｍ（ｊ，ｆ）を足し合わせて周波数領域の予測エコー信号Ｙ＾（ｊ，ｆ）を次式のように得る。
【数８】
2006121589.tif 000006
【００１５】
ＦＴ変換部７０３において、周波数領域の予測エコーを逆ＦＦＴにより時間領域に変換し、時間領域の予測エコーｙ＾（ｋ）を次のように得る。
【数９】
2006121589.tif 000007
【００１６】
減算部７１では、フレーム化部７１１において時間領域の収音信号ｙ（ｋ）をフレーム化し、減算器７１２において予測エコー信号ｙ＾（ｋ）のフレームを差し引いて時間領域での残差信号ｅ（ｋ）のフレームを求める。
【数１０】
2006121589.tif 000008
【００１７】
ＴＦ変換部７２では、次式のように時間領域の残差信号ｅ（ｋ）のフレームにＬＤ個の０を詰めて周波数領域に変換することで周波数領域の残差信号Ｅ（ｊ，ｆ）を求める。
【数１１】
2006121589.tif 000009
【００１８】
エコー成分比率推定部７３において以下のステップＦ１〜Ｆ７により、周波数領域の多チャネル再生信号Ｘｍ（ｊ，ｆ）と周波数領域の残差信号Ｅ（ｊ，ｆ）から、周波数成分ごとに残差信号に含まれるエコー成分の比率Ｓ（ｊ，ｆ）を求める。図８にエコー成分比率Ｓ（ｊ，ｆ）を推定するためのフローを示す。
ステップＦ１
現時点のフレームから求めた多チャネル再生信号の短時間スペクトルＸ１（ｊ，ｆ），…，ＸＭ（ｊ，ｆ）を図７に示す相関除去部７３１内の現時点の蓄積部７３１ａ１に保存する。
【００１９】
ステップＦ２
相関除去部７３１ｂ１では、例えば次式の方法で多チャネル再生信号の短時間スペクトルＸ２（ｊ，ｆ），…，ＸＭ（ｊ，ｆ）からＸ１（ｊ，ｆ）との相関成分を除去して、スペクトルＸ２（１）（ｊ，ｆ），…，ＸＭ（１）（ｊ，ｆ）を得、多チャネル再生信号スペクトルの副成分の一部とする。
【数１２】
2006121589.tif 000010
ここで、ε［］は、平均をとることを意味し、平均処理の一例としては、
【数１３】
2006121589.tif 000011
のように、１フレーム前の処理結果と０〜１の値をとる平滑化定数βを用いる方法がある。
【００２０】
ステップＦ３
相関除去部７３１ｂ２において、過去の蓄積部７３１ａ２に蓄積された１フレーム前の多チャネル再生信号のスペクトルＸ１（ｊ−１，ｆ），…，ＸＭ（ｊ−１，ｆ）から、Ｘ１（ｊ，ｆ）との相関を次のように除去したスペクトルＸ１（１）（ｊ−１，ｆ），…，ＸＭ（１）（ｊ−１，ｆ）を求め、多チャネル再生信号スペクトルの副成分の一部とする。
【数１４】
2006121589.tif 000012
なお、ｎフレーム前の短時間スペクトルＸ１（ｊ−ｎ，ｆ），…，ＸＭ（ｊ−ｎ，ｆ）をエコー成分比率推定に使用する場合にも、同様の計算により得られた結果を多チャネル再生信号スペクトルの副成分の一部とすればよい。
【００２１】
ステップＦ４
相関除去部７３２では、現時点のフレームの残差信号の短時間スペクトルＥ（ｊ，ｆ）からＸ１（ｊ，ｆ）との相関成分を除去したスペクトルＥ（１）（ｊ，ｆ）を求める。
【数１５】
2006121589.tif 000013
【００２２】
ステップＦ５
コヒーレンス算出部７３３１では、多チャネル再生信号スペクトルの主成分である現時点のフレームの第１チャネル再生信号の短時間スペクトルＸ１（ｊ，ｆ）と現時点の残差信号のスペクトルＥ（ｊ，ｆ）から、次のコヒーレンスを求める。
【数１６】
2006121589.tif 000014
【００２３】
ステップＦ６
副成分エコー比率算出部７３３２では、まず相関除去された残差信号スペクトルＥ（１）（ｊ，ｆ）に含まれるエコー成分Ｅ＾（１）（ｊ，ｆ）を求める。エコー成分Ｅ＾（１）（ｊ，ｆ）は、多チャネル再生信号短時間スペクトルの副成分Ｘ２（１）（ｊ，ｆ），…，ＸＭ（１）（ｊ，ｆ），Ｘ１（１）（ｊ−１，ｆ），…，ＸＭ（１）（ｊ−１，ｆ）の線形和
【数１７】
2006121589.tif 000015
のうちで、相関除去された残差信号スペクトルとの誤差
｜Ｅ（１）（ｊ，ｆ）−Ｅ＾（１）（ｊ，ｆ）｜２
が最小となるスペクトルである。この誤差を最小にするスペクトルは、
【数１８】
2006121589.tif 000016
とし、
【００２４】
【数１９】
2006121589.tif 000017
により求められる。さらに、次式により副成分のエコー比率を周波数成分ごとに求める。
【数２０】
2006121589.tif 000018
【００２５】
ステップＦ７
図７のエコー成分比率算出部７３４において、ステップＦ５、Ｆ６で求めた各比率から、残差信号スペクトルＥ（ｊ，ｆ）に占めるエコー成分の比率を求める。
【数２１】
2006121589.tif 000019
【００２６】
図６のフィルタ更新部７４では、まず周波数成分ごとに残差信号Ｅ（ｊ，ｆ）と再生信号Ｘｍ（ｊ，ｆ）から、周波数領域の適応フィルタ係数の修正量ｄＷｍ（ｊ，ｄ，ｆ）を求める。次に周波数領域の残差信号の各周波数成分に含まれる残留エコー成分の比率Ｓ（ｊ，ｆ）を使用して適応フィルタ係数の修正量ｄＷｍ（ｊ，ｄ，ｆ）の大きさを制御しながら、周波数領域の適応フィルタ係数Ｗｍ（ｊ，ｄ，ｆ）を逐次更新する。具体的には、例えば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.の提案と組み合わせ、以下のように修正量ｄＷｍ（ｊ，ｄ，ｆ）および適応フィルタ係数Ｗｍ（ｊ，ｄ，ｆ）を求めることができる。
【数２２】
2006121589.tif 000020
ただし、ｄ＝１，…，Ｄ、Ｘｍ＊（ｊ，ｆ）はＸｍ（ｊ，ｆ）の複素共役である。
【００２７】
【数２３】
2006121589.tif 000021
ただし、ｍ＝１，…，Ｍ、ｄ＝１，…，Ｄ、ｆ＝１，…，２ＬＤ、μは０〜１の値をとるステップサイズ、δは正則化定数である。また、ｐ（ｊ，ｆ）は次式により計算される。
【数２４】
2006121589.tif 000022
【００２８】
このように適応フィルタ係数Ｗｍ（ｊ，ｄ，ｆ）を処理単位ごとに更新することで、次の処理では更新された適応フィルタ係数が使用できる。
現時点のフレームの処理が終了すると、最後に現時点の蓄積部７３１ａ１に蓄積された再生信号情報は過去の蓄積部７３１ａ２に転送され、蓄積される。
なお、蓄積部７３１ａ内で現時点の蓄積部７３１ａ１と過去の蓄積部７３１ａ２とを特に区別し、上記のように一連の処理の最後に現時点の蓄積部７３１ａ１に蓄積された再生信号情報を過去の蓄積部７３１ａ２に転送するのではなく、１つの蓄積部７３１ａに蓄積された情報の中で最新情報を現時点の情報として処理する方法もある。また、図９に示すように処理に利用する現時点の再生信号のスペクトルを、蓄積部から取り出すのではなく、入力された再生信号のスペクトルを直接利用する方法もある。
【００２９】
［第２実施形態］
この発明は、第１実施形態に対してエコー経路推定速度の向上を図ったものであり、図１０に全体構成を、図１１にＭチャネルエコーキャンセル部の構成を示す。
図１０の相関変動処理部９ｍ（ｍ＝１，…，Ｍ）では、例えば次式のようにチャネル毎に受話信号ｕｍ（ｋ）を乱数で振幅変調して元の信号に付加して相互相関が絶えず変動している再生信号ｘｍ（ｋ）を生成し、各スピーカから再生する。
【数２５】
2006121589.tif 000023
【００３０】
この相関変動処理がない場合、Ｍチャネルの受話信号のチャネル間相関が高いために、第ｍチャネル再生信号と第１チャネルの収音信号の入出力関係から、第ｍチャネルの真のエコー経路特性が必ず求められるとは限らない。しかし、この相関変動処理を受話信号に適用することにより、Ｍチャネルのエコー経路特性が確実に推定できる。
図１１のＭチャネルエコーキャンセル部８では、付加信号成分が受話信号成分よりも強調されたチャネル間相関の小さい信号のベクトルを信号変換部８５により生成し、チャネル間相関の大きい再生信号ベクトルの代わりに使用して適応フィルタを更新する。この更新には、例えば特許文献２で提案されている適応アルゴリズムを使用し、そのステップサイズを残差信号に含まれる残留エコー成分の比率Ｓ（ｊ，ｆ）で制御する。チャネル間相関の小さい修正ベクトルを使用することで、適応フィルタによるエコー経路の収束速度が向上する。
【図面の簡単な説明】
【００３１】
【図１】Ｍ入力Ｎ出力のエコー消去装置の一般的構成を示す図。
【図２】Ｍ入力１出力のエコー消去装置の一般的構成を示す図。
【図３】従来のＭチャネルエコーキャンセル部の構成を示す図。
【図４】従来のエコー成分比率推定部の構成を示す図。
【図５】従来のエコー成分比率推定のフローを示す図。
【図６】第１実施形態のＭチャネルエコーキャンセル部の構成を示す図。
【図７】第１実施形態のエコー成分比率推定部の構成を示す図。
【図８】第１実施形態のエコー成分比率推定のフローを示す図。
【図９】第１実施形態のエコー成分比率推定部の変形例の構成を示す図。
【図１０】第２実施形態のＭ入力１出力チャネルのエコー消去装置の一般的構成を示す図。
【図１１】第２実施形態のＭチャネルエコーキャンセル部の構成を示す図。
【図１２】信号のサンプル時刻ｋとフレーム時刻ｊの関係を示す図。
【図１】
2006121589.tif 000028
【図２】
2006121589.tif 000029
【図３】
2006121589.tif 000030
【図４】
2006121589.tif 000031
【図５】
2006121589.tif 000032
【図６】
2006121589.tif 000033
【図７】
2006121589.tif 000034
【図８】
2006121589.tif 000035
【図９】
2006121589.tif 000036
【図１０】
2006121589.tif 000037
【図１１】
2006121589.tif 000038
【図１２】
2006121589.tif 000039