JP2002261659A - 多チャネルエコーキャンセル方法、その装置、そのプログラム及び記録媒体 - Google Patents
多チャネルエコーキャンセル方法、その装置、そのプログラム及び記録媒体Info
- Publication number
- JP2002261659A JP2002261659A JP2001051387A JP2001051387A JP2002261659A JP 2002261659 A JP2002261659 A JP 2002261659A JP 2001051387 A JP2001051387 A JP 2001051387A JP 2001051387 A JP2001051387 A JP 2001051387A JP 2002261659 A JP2002261659 A JP 2002261659A
- Authority
- JP
- Japan
- Prior art keywords
- channel
- signal
- echo
- impulse response
- step size
- 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.)
- Granted
Links
Landscapes
- Cable Transmission Systems, Equalization Of Radio And Reduction Of Echo (AREA)
- Telephone Function (AREA)
Abstract
(57)【要約】
【課題】 収束速度を速くする。
【解決手段】 スピーカ１２1 ，１２2 とマイクロホン
１６1 間のインパルス応答ｈ1 （ｋ），ｈ2 （ｋ）
の各要素の２乗和を求め（３２1 ，３２2 ）、その平方
根を求め（３３1 ，３３2 ）、μ1 ＝２ａ／（ａ＋
ｂ），μ2 ＝２ｂ／（ａ＋ｂ）を計算し（３５1 ，３５
2 ）、受話信号ｘ1 （ｋ），ｘ2 （ｋ）の一定時間の２
乗和ＸT （ｋ）Ｘ（ｋ）を求め（２７）、推定部１
９1 ，１９2 でそれぞれｈ＾1 （ｋ＋１）＝ｈ＾1
（ｋ）＋（μ1 ｅ（ｋ）／（ＸT （ｋ）Ｘ
（ｋ）））Ｘ1 （ｋ），ｈ＾2 （ｋ＋１）＝ｈ＾
2 （ｋ）＋（μ2ｅ（ｋ）／（ＸT （ｋ）Ｘ
（ｋ）））Ｘ2 （ｋ）を計算して疑似反響路１８1 ，
１８2 を構成し、これにｘ1 （ｋ），ｘ2 （ｋ）を通
し、疑似反響信号ｙ＾1 （ｋ），ｙ＾2 （ｋ）を得、マ
イクロホン１６の出力ｙ（ｋ）からｙ＾1 （ｋ）とｙ＾
2 （ｋ）を引き、ｅ（ｋ）を得る。
１６1 間のインパルス応答ｈ1 （ｋ），ｈ2 （ｋ）
の各要素の２乗和を求め（３２1 ，３２2 ）、その平方
根を求め（３３1 ，３３2 ）、μ1 ＝２ａ／（ａ＋
ｂ），μ2 ＝２ｂ／（ａ＋ｂ）を計算し（３５1 ，３５
2 ）、受話信号ｘ1 （ｋ），ｘ2 （ｋ）の一定時間の２
乗和ＸT （ｋ）Ｘ（ｋ）を求め（２７）、推定部１
９1 ，１９2 でそれぞれｈ＾1 （ｋ＋１）＝ｈ＾1
（ｋ）＋（μ1 ｅ（ｋ）／（ＸT （ｋ）Ｘ
（ｋ）））Ｘ1 （ｋ），ｈ＾2 （ｋ＋１）＝ｈ＾
2 （ｋ）＋（μ2ｅ（ｋ）／（ＸT （ｋ）Ｘ
（ｋ）））Ｘ2 （ｋ）を計算して疑似反響路１８1 ，
１８2 を構成し、これにｘ1 （ｋ），ｘ2 （ｋ）を通
し、疑似反響信号ｙ＾1 （ｋ），ｙ＾2 （ｋ）を得、マ
イクロホン１６の出力ｙ（ｋ）からｙ＾1 （ｋ）とｙ＾
2 （ｋ）を引き、ｅ（ｋ）を得る。
Description
【０００１】
【発明の属する技術分野】この発明は、多チャネル再生
系を有する例えば通信会議システムにおいて、ハウリン
グの原因および聴覚上の障害となる室内反響信号を消去
する多チャネル反響消去方法、その装置、そのプログラ
ム及びその記録媒体に関するものである。
系を有する例えば通信会議システムにおいて、ハウリン
グの原因および聴覚上の障害となる室内反響信号を消去
する多チャネル反響消去方法、その装置、そのプログラ
ム及びその記録媒体に関するものである。
【０００２】
【従来の技術】同時通話性能に優れ反響感の少ない拡声
通話システムを提供するために、反響消去装置がある。
まず、１チャネル用の反響消去装置について、その反響
消去方法および装置構成を図４を参照して説明する。拡
声通話において、受話端子１１に相手の発話等で得られ
る受話信号は、そのままスピーカ１２から再生される場
合と、スピーカ１２へ送る前に、受話信号の振幅やパワ
ー等の大きさに応じて自動的に利得を調節するなど、受
話信号に何らかの加工が施された後に、スピーカ１２か
ら再生される場合とがある。このため、この明細書で受
話信号ｘ1(ｋ）とは、相手からの受話信号そのものとは
限らず、受話信号に対し加工が施された場合は、その加
工された後の受話信号を指すものとする。ｋは離散時間
を表す。反響消去装置１４は、受話信号ｘ1(ｋ）がスピ
ーカ１２から反響路１５を経て、マイクロホン１６に集
音されて得られる反響信号ｙ1(ｋ）を消去する。ここ
で、反響信号ｙ1(ｋ）は、時刻ｋにおける反響路１５の
インパルス応答をｈ11（ｋ，ｎ）として、 ｙ1(ｋ）＝Σｈ11（ｋ，ｎ）ｘ1(ｋ−ｎ） （１）Σ はｎ＝０からＬ−１までのような畳み込み演算で得
られるものとモデル化できる。Ｌはタップ数で、反響路
１５の残響時間に対応させて、あらかじめ設定しておく
定数である。まず、受話信号蓄積・ベクトル生成部１７
において、受話信号ｘ1(ｋ）をＬ−１時刻過去のものま
で蓄積しておく。蓄積された信号は、受話信号ベクトル
ｘ1(ｋ）、すなわち ｘ1(ｋ）＝［ｘ1(ｋ),ｘ1(ｋ−１),…, ｘ1(ｋ−Ｌ＋１)]T （２） として出力される。但し、＊T はベクトルの転置を表
す。疑似反響信号生成部１８では、式（２）の受話信号
ベクトルｘ1(ｋ）と、反響路推定部１９から得られる
疑似反響ベクトルｈ＾11（ｋ）との内積演算 ｙ＾1(ｋ）＝ｈ＾11 T （ｋ）ｘ1(ｋ） （３） を行ない、その結果として、疑似反響信号ｙ＾1(ｋ）を
生成する。この内積演算は、式（１）のような畳み込み
演算と等価である。反響路推定部１９では、疑似反響信
号生成部１８で用いる疑似反響路ベクトルｈ＾
11（ｋ）を生成する。この反響路推定に用いる最も一般
的なアルゴリズムは、ＮＬＭＳアルゴリズム（学習同定
法）である。ＮＬＭＳアルゴリズムでは、時刻ｋにおけ
る受話信号ベクトルｘ1(ｋ）と、残留反響信号ｅ
1(ｋ）、すなわち差回路２１でマイクロホン１６の出力
ｙ1(ｋ）から疑似反響信号ｙ＾1(ｋ）を差し引いた信号 ｅ1(ｋ）＝ｙ1(ｋ）−ｙ＾1(ｋ） （４） とから、時刻ｋ＋１において用いる疑似反響路ベクトル
ｈ＾11（ｋ＋１）を次式のように求める。ｈ ＾11(k+1) ＝ｈ＾11(k) ＋μｅ1(k)ｘ1(k)／（ｘ1 T (k) ｘ1(k)） （５） 但し、μはステップサイズパラメータと呼ばれ０＜μ＜
２の範囲で適応動作の調整に用いる。以上のような処理
を繰り返すことにより、反響路推定部１９では、次第に
疑似反響路ベクトルｈ＾11（ｋ）を、真の反響路１５
のインパルス応答ｈ11（ｋ，ｎ）の時系列を各要素とし
て持つ反響路ベクトルｈ11（ｋ）、すなわち ｈ11(k) ＝［ｈ11(k,0),ｈ11(k,1),…, ｈ11(k,L-1)]T （６） と一致させることが可能となり、その結果式（４）の残
留反響信号ｅ1 （ｋ）を小さくすることができる。
通話システムを提供するために、反響消去装置がある。
まず、１チャネル用の反響消去装置について、その反響
消去方法および装置構成を図４を参照して説明する。拡
声通話において、受話端子１１に相手の発話等で得られ
る受話信号は、そのままスピーカ１２から再生される場
合と、スピーカ１２へ送る前に、受話信号の振幅やパワ
ー等の大きさに応じて自動的に利得を調節するなど、受
話信号に何らかの加工が施された後に、スピーカ１２か
ら再生される場合とがある。このため、この明細書で受
話信号ｘ1(ｋ）とは、相手からの受話信号そのものとは
限らず、受話信号に対し加工が施された場合は、その加
工された後の受話信号を指すものとする。ｋは離散時間
を表す。反響消去装置１４は、受話信号ｘ1(ｋ）がスピ
ーカ１２から反響路１５を経て、マイクロホン１６に集
音されて得られる反響信号ｙ1(ｋ）を消去する。ここ
で、反響信号ｙ1(ｋ）は、時刻ｋにおける反響路１５の
インパルス応答をｈ11（ｋ，ｎ）として、 ｙ1(ｋ）＝Σｈ11（ｋ，ｎ）ｘ1(ｋ−ｎ） （１）Σ はｎ＝０からＬ−１までのような畳み込み演算で得
られるものとモデル化できる。Ｌはタップ数で、反響路
１５の残響時間に対応させて、あらかじめ設定しておく
定数である。まず、受話信号蓄積・ベクトル生成部１７
において、受話信号ｘ1(ｋ）をＬ−１時刻過去のものま
で蓄積しておく。蓄積された信号は、受話信号ベクトル
ｘ1(ｋ）、すなわち ｘ1(ｋ）＝［ｘ1(ｋ),ｘ1(ｋ−１),…, ｘ1(ｋ−Ｌ＋１)]T （２） として出力される。但し、＊T はベクトルの転置を表
す。疑似反響信号生成部１８では、式（２）の受話信号
ベクトルｘ1(ｋ）と、反響路推定部１９から得られる
疑似反響ベクトルｈ＾11（ｋ）との内積演算 ｙ＾1(ｋ）＝ｈ＾11 T （ｋ）ｘ1(ｋ） （３） を行ない、その結果として、疑似反響信号ｙ＾1(ｋ）を
生成する。この内積演算は、式（１）のような畳み込み
演算と等価である。反響路推定部１９では、疑似反響信
号生成部１８で用いる疑似反響路ベクトルｈ＾
11（ｋ）を生成する。この反響路推定に用いる最も一般
的なアルゴリズムは、ＮＬＭＳアルゴリズム（学習同定
法）である。ＮＬＭＳアルゴリズムでは、時刻ｋにおけ
る受話信号ベクトルｘ1(ｋ）と、残留反響信号ｅ
1(ｋ）、すなわち差回路２１でマイクロホン１６の出力
ｙ1(ｋ）から疑似反響信号ｙ＾1(ｋ）を差し引いた信号 ｅ1(ｋ）＝ｙ1(ｋ）−ｙ＾1(ｋ） （４） とから、時刻ｋ＋１において用いる疑似反響路ベクトル
ｈ＾11（ｋ＋１）を次式のように求める。ｈ ＾11(k+1) ＝ｈ＾11(k) ＋μｅ1(k)ｘ1(k)／（ｘ1 T (k) ｘ1(k)） （５） 但し、μはステップサイズパラメータと呼ばれ０＜μ＜
２の範囲で適応動作の調整に用いる。以上のような処理
を繰り返すことにより、反響路推定部１９では、次第に
疑似反響路ベクトルｈ＾11（ｋ）を、真の反響路１５
のインパルス応答ｈ11（ｋ，ｎ）の時系列を各要素とし
て持つ反響路ベクトルｈ11（ｋ）、すなわち ｈ11(k) ＝［ｈ11(k,0),ｈ11(k,1),…, ｈ11(k,L-1)]T （６） と一致させることが可能となり、その結果式（４）の残
留反響信号ｅ1 （ｋ）を小さくすることができる。
【０００３】一般にＮ（＞２）チャネルの再生系とＭ
（＞１）チャネルの集音系とで構成される通信会議シス
テムの場合の反響の消去は、図５に示すような構成によ
り行われる。すなわち再生側の全Ｎチャネルと集音側の
各１チャネルとの間にＮ入力１出力時系列信号を処理す
るＮチャネル反響消去装置２２1,２２2,…，２２M をそ
れぞれ接続した反響消去システム２３として実現され
る。この場合システム全体でＮ×Ｍ個の反響路１５
nm（１＜ｎ＜Ｎ，１＜ｍ＜Ｍ）が存在する。このシステ
ムの構成単位である再生側の全Ｎチャネルと集音側の各
１チャネルとの間に接続されるＮチャネル反響消去装置
２２1,２２2,…，２２M については、図４に示した反響
消去装置１４の構成を拡張して、図６に示すように構成
される。これは例えば電子情報通信学会論文誌 '８６／
１０ Ｖｏｌ．Ｊ６９−Ａ Ｎｏ．１０「多チャネル適
応ディジタルフィルタ」に詳しく述べられている。ここ
で、集音側が第ｍ集音チャネル（１＜ｍ＜Ｍ）に接続さ
れているＮチャネル反響消去装置２２m を考える。第ｍ
チャネルのマイクロホン１６m で集音される反響信号
は、各再生チャネルの受話信号がそれぞれの反響路１５
1m〜１５Nmを経て集音側で全て加算されることにより得
られるために、反響路推定をどの再生チャネルについて
も、統一的に同じ１つの残留反響信号ｅm （ｋ）のみを
評価して行なうための工夫が必要となる。まず、各再生
チャネルの受話信号について、受話信号蓄積・ベクトル
生成部（１７1,１７2,…，１７N ）により、受話信号ベ
クトル ｘ1(ｋ）＝［ｘ1(ｋ),ｘ1(ｋ−１),…, ｘ1(ｋ−Ｌ1 ＋１)]T （７） ｘ2(ｋ）＝［ｘ2(ｋ),ｘ2(ｋ−１),…, ｘ2(ｋ−Ｌ2 ＋１)]T （８） ： ｘN (k）＝［ｘN (k),ｘN (k−1), …, ｘN (k−ＬN ＋１)]T （９） を生成する。但し、Ｌ1,Ｌ2,…，ＬN はタップ数で、各
反響路１５1m，１５2m，…，１５Nmの残響時間に対応さ
せて、あらかじめ設定する定数である。これらのベクト
ルをベクトル結合部２４によって、 ｘ(k) ＝［ｘ1 T (k),ｘ2 T (k),…，ｘN T (k)]T （10） と結合する。また、反響路推定部１９m においても、各
再生チャネルと第ｍ集音チャネルとの間のＮ個の反響路
を模疑するための、各疑似反響路ベクトルｈ＾ 1m(k),
ｈ＾2m(k),…，ｈ＾Nm(k) を結合して ｈ＾m (k) ＝［ｈ＾1m T (k),ｈ＾2m T (k),…，ｈ＾Nm T (k)]T （11） として扱う。疑似反響路結合ベクトルｈ＾m (k) の更
新は、ＮＬＭＳアルゴリズムを用いた場合、 ｈ＾m (k+1)=ｈ＾m (k)+μｅm (k) ｘ(k) ／（ｘT (k) ｘ(k)) （12） のように行なわれる。疑似反響信号生成部１８m では、
内積演算 ｙ＾m （ｋ）＝ｈ＾m T (k) ｘ（ｋ） （13） により、第ｍ集音チャネルで集音された反響信号ｙ
m （ｋ）に対する疑似反響信号ｙ＾m （ｋ）を生成す
る。このように、各チャネル毎のベクトルを結合して１
つのベクトルとして扱うことにより、基本的な処理の流
れは、図４に示した１チャネル反響消去装置と同様とな
る。
（＞１）チャネルの集音系とで構成される通信会議シス
テムの場合の反響の消去は、図５に示すような構成によ
り行われる。すなわち再生側の全Ｎチャネルと集音側の
各１チャネルとの間にＮ入力１出力時系列信号を処理す
るＮチャネル反響消去装置２２1,２２2,…，２２M をそ
れぞれ接続した反響消去システム２３として実現され
る。この場合システム全体でＮ×Ｍ個の反響路１５
nm（１＜ｎ＜Ｎ，１＜ｍ＜Ｍ）が存在する。このシステ
ムの構成単位である再生側の全Ｎチャネルと集音側の各
１チャネルとの間に接続されるＮチャネル反響消去装置
２２1,２２2,…，２２M については、図４に示した反響
消去装置１４の構成を拡張して、図６に示すように構成
される。これは例えば電子情報通信学会論文誌 '８６／
１０ Ｖｏｌ．Ｊ６９−Ａ Ｎｏ．１０「多チャネル適
応ディジタルフィルタ」に詳しく述べられている。ここ
で、集音側が第ｍ集音チャネル（１＜ｍ＜Ｍ）に接続さ
れているＮチャネル反響消去装置２２m を考える。第ｍ
チャネルのマイクロホン１６m で集音される反響信号
は、各再生チャネルの受話信号がそれぞれの反響路１５
1m〜１５Nmを経て集音側で全て加算されることにより得
られるために、反響路推定をどの再生チャネルについて
も、統一的に同じ１つの残留反響信号ｅm （ｋ）のみを
評価して行なうための工夫が必要となる。まず、各再生
チャネルの受話信号について、受話信号蓄積・ベクトル
生成部（１７1,１７2,…，１７N ）により、受話信号ベ
クトル ｘ1(ｋ）＝［ｘ1(ｋ),ｘ1(ｋ−１),…, ｘ1(ｋ−Ｌ1 ＋１)]T （７） ｘ2(ｋ）＝［ｘ2(ｋ),ｘ2(ｋ−１),…, ｘ2(ｋ−Ｌ2 ＋１)]T （８） ： ｘN (k）＝［ｘN (k),ｘN (k−1), …, ｘN (k−ＬN ＋１)]T （９） を生成する。但し、Ｌ1,Ｌ2,…，ＬN はタップ数で、各
反響路１５1m，１５2m，…，１５Nmの残響時間に対応さ
せて、あらかじめ設定する定数である。これらのベクト
ルをベクトル結合部２４によって、 ｘ(k) ＝［ｘ1 T (k),ｘ2 T (k),…，ｘN T (k)]T （10） と結合する。また、反響路推定部１９m においても、各
再生チャネルと第ｍ集音チャネルとの間のＮ個の反響路
を模疑するための、各疑似反響路ベクトルｈ＾ 1m(k),
ｈ＾2m(k),…，ｈ＾Nm(k) を結合して ｈ＾m (k) ＝［ｈ＾1m T (k),ｈ＾2m T (k),…，ｈ＾Nm T (k)]T （11） として扱う。疑似反響路結合ベクトルｈ＾m (k) の更
新は、ＮＬＭＳアルゴリズムを用いた場合、 ｈ＾m (k+1)=ｈ＾m (k)+μｅm (k) ｘ(k) ／（ｘT (k) ｘ(k)) （12） のように行なわれる。疑似反響信号生成部１８m では、
内積演算 ｙ＾m （ｋ）＝ｈ＾m T (k) ｘ（ｋ） （13） により、第ｍ集音チャネルで集音された反響信号ｙ
m （ｋ）に対する疑似反響信号ｙ＾m （ｋ）を生成す
る。このように、各チャネル毎のベクトルを結合して１
つのベクトルとして扱うことにより、基本的な処理の流
れは、図４に示した１チャネル反響消去装置と同様とな
る。
【０００４】
【発明が解決しようとする課題】従来の多チャネルエコ
ー消去方法はその疑似反響路が真の反響路にほぼ収束す
るのが遅い。この点を改善すべく研究した結果、各チャ
ネルの受話信号の再生信号ごとの反響路についての疑似
反響路において、その反響路の逐次推定における更新ス
テップサイズμが同一値である点に問題があると考える
に至った。説明を簡単にするため、図７に示すように受
話チャネルが２つ、集音チャネルが１つの場合を例とす
る。各受話信号ｘ1 （ｋ），ｘ2 （ｋ）はそれぞれチャ
ネル疑似反響路１８1 ，１８2 へ印加され、チャネル疑
似反響路１８1 ，１８2 よりの各チャネル疑似反響信号
ｙ＾1 （ｋ），ｙ＾2 （ｋ）は合成回路２６で合成され
て総合疑似反響信号ｙ＾（ｋ）とされて差回路２１へ供
給され、マイクロホン１６1 よりの集音信号ｙ（ｋ）か
らｙ＾（ｋ）が引算され、残留反響信号ｅ（ｋ）が求め
られる。残留反響信号ｅ（ｋ）と、受話信号ｘ
1 （ｋ），ｘ2 （ｋ）がチャネル反響路推定部１９1 ，
１９2 へ供給され、チャネル反響路推定部１９ 1 ，１９
2 はそれぞれ、受話信号ｘ1 （ｋ），ｘ2 （ｋ）の各再
生信号がマイクロホン１６1 に至る反響路１５1 ，１５
2 のインパルス応答ｈ1 （ｋ），ｈ 2 （ｋ）を逐次
推定して、これら推定したインパルス応答ｈ＾
1 （ｋ），ｈ＾2 （ｋ）をチャネル疑似反響路１
８1 ，１８2 にそれぞれ設定する。
ー消去方法はその疑似反響路が真の反響路にほぼ収束す
るのが遅い。この点を改善すべく研究した結果、各チャ
ネルの受話信号の再生信号ごとの反響路についての疑似
反響路において、その反響路の逐次推定における更新ス
テップサイズμが同一値である点に問題があると考える
に至った。説明を簡単にするため、図７に示すように受
話チャネルが２つ、集音チャネルが１つの場合を例とす
る。各受話信号ｘ1 （ｋ），ｘ2 （ｋ）はそれぞれチャ
ネル疑似反響路１８1 ，１８2 へ印加され、チャネル疑
似反響路１８1 ，１８2 よりの各チャネル疑似反響信号
ｙ＾1 （ｋ），ｙ＾2 （ｋ）は合成回路２６で合成され
て総合疑似反響信号ｙ＾（ｋ）とされて差回路２１へ供
給され、マイクロホン１６1 よりの集音信号ｙ（ｋ）か
らｙ＾（ｋ）が引算され、残留反響信号ｅ（ｋ）が求め
られる。残留反響信号ｅ（ｋ）と、受話信号ｘ
1 （ｋ），ｘ2 （ｋ）がチャネル反響路推定部１９1 ，
１９2 へ供給され、チャネル反響路推定部１９ 1 ，１９
2 はそれぞれ、受話信号ｘ1 （ｋ），ｘ2 （ｋ）の各再
生信号がマイクロホン１６1 に至る反響路１５1 ，１５
2 のインパルス応答ｈ1 （ｋ），ｈ 2 （ｋ）を逐次
推定して、これら推定したインパルス応答ｈ＾
1 （ｋ），ｈ＾2 （ｋ）をチャネル疑似反響路１
８1 ，１８2 にそれぞれ設定する。
【０００５】チャネル反響路推定部１９1 ，１９2 にお
ける各推定チャネルインパルス応答ｈ＾1 （ｋ），
ｈ＾2 （ｋ）は式（１２）から以下のようになる。式
（１２）は受話チャネルが２つであるから式（１４）と
なる。
ける各推定チャネルインパルス応答ｈ＾1 （ｋ），
ｈ＾2 （ｋ）は式（１２）から以下のようになる。式
（１２）は受話チャネルが２つであるから式（１４）と
なる。
【０００６】
【数１】
【０００７】反響路１５1 ，１５2 の残響時間により決
るタップ数Ｌ1 とＬ2 を等しくＬとした。このように
ＸT （ｋ）Ｘ（ｋ）はスカラーであるから、式（１
４）は行列の各行ごとに次のように書ける。 ｈ＾1 （ｋ＋１）＝ｈ＾1 （ｋ）＋（μｅ（ｋ）／（ＸT （ｋ）Ｘ（ ｋ）））Ｘ1 （ｋ） （15） ｈ＾2 （ｋ＋１）＝ｈ＾2 （ｋ）＋（μｅ（ｋ）／（ＸT （ｋ）Ｘ（ ｋ）））Ｘ2 （ｋ） （16） つまりチャネル反響路推定部１９1 ，１９2 はそれぞれ
式（１５），（１６）によりチャネル反響路１５1 ，１
５2 のインパルス応答ｈ1 （ｋ），ｈ2 （ｋ）を逐
次推定することになる。Ｘ1 （ｋ）／（ＸT （ｋ）
Ｘ（ｋ）），Ｘ2 （ｋ）／（ＸT （ｋ）Ｘ
（ｋ））はそれぞれ全チャネルの受話信号のパワーで規
格化された受話信号（ベクトル）と言える。Ｘ
T （ｋ）Ｘ（ｋ）はパワー計算部２７で計算されて各
推定部１９1 ，１９2 へ供給される。また図６中には図
５中の受話信号蓄積ベクトル生成部は必要とするが省略
してある。式（１５），（１６）より、その更新ステッ
プサイズは共にμと等しいものとなっている。
るタップ数Ｌ1 とＬ2 を等しくＬとした。このように
ＸT （ｋ）Ｘ（ｋ）はスカラーであるから、式（１
４）は行列の各行ごとに次のように書ける。 ｈ＾1 （ｋ＋１）＝ｈ＾1 （ｋ）＋（μｅ（ｋ）／（ＸT （ｋ）Ｘ（ ｋ）））Ｘ1 （ｋ） （15） ｈ＾2 （ｋ＋１）＝ｈ＾2 （ｋ）＋（μｅ（ｋ）／（ＸT （ｋ）Ｘ（ ｋ）））Ｘ2 （ｋ） （16） つまりチャネル反響路推定部１９1 ，１９2 はそれぞれ
式（１５），（１６）によりチャネル反響路１５1 ，１
５2 のインパルス応答ｈ1 （ｋ），ｈ2 （ｋ）を逐
次推定することになる。Ｘ1 （ｋ）／（ＸT （ｋ）
Ｘ（ｋ）），Ｘ2 （ｋ）／（ＸT （ｋ）Ｘ
（ｋ））はそれぞれ全チャネルの受話信号のパワーで規
格化された受話信号（ベクトル）と言える。Ｘ
T （ｋ）Ｘ（ｋ）はパワー計算部２７で計算されて各
推定部１９1 ，１９2 へ供給される。また図６中には図
５中の受話信号蓄積ベクトル生成部は必要とするが省略
してある。式（１５），（１６）より、その更新ステッ
プサイズは共にμと等しいものとなっている。
【０００８】しかし、例えば図７に示すように、スピー
カ１２1 ，１２2 とマイクロホン１６1 ，１６2 を互い
に対向させて、方形の各コーナに位置させた場合、マイ
クロホン１６1 に対し正面のスピーカ１２1 からの反響
路１５1 のインパルス応答ｈ 1 （ｋ）の方が、対角位置
のスピーカ１２2 からの反響路１５2 のインパルス応答
ｈ2 （ｋ）よりエネルギーが大きい。このことは、ス
ピーカとマイクロホンの方形コーナの配置に限らず、他
の配置においても反響路によりインパルス応答のエネル
ギーに差があり、かつ３チャネル以上の多チャネルエコ
ー消去方法についても同様のことが言える。
カ１２1 ，１２2 とマイクロホン１６1 ，１６2 を互い
に対向させて、方形の各コーナに位置させた場合、マイ
クロホン１６1 に対し正面のスピーカ１２1 からの反響
路１５1 のインパルス応答ｈ 1 （ｋ）の方が、対角位置
のスピーカ１２2 からの反響路１５2 のインパルス応答
ｈ2 （ｋ）よりエネルギーが大きい。このことは、ス
ピーカとマイクロホンの方形コーナの配置に限らず、他
の配置においても反響路によりインパルス応答のエネル
ギーに差があり、かつ３チャネル以上の多チャネルエコ
ー消去方法についても同様のことが言える。
【０００９】そこでこの発明の目的はこのインパルス応
答のエネルギーの差を考慮して従来よりも収束が速い多
チャネルエコーキャンセル方法、その装置、そのプログ
ラム及びその記録媒体を提供することにある。
答のエネルギーの差を考慮して従来よりも収束が速い多
チャネルエコーキャンセル方法、その装置、そのプログ
ラム及びその記録媒体を提供することにある。
【００１０】
【課題を解決するための手段】この発明によれば各チャ
ネルの受話信号をチャネル疑似反響路に印加してチャネ
ル疑似反響信号をそれぞれ生成し、これらチャネル疑似
反響信号を加算して総合疑似反響信号を生成し、この総
合疑似反響信号を集音信号から差し引いて残留反響信号
を生成し、特に各チャネルの受話信号の再生信号と集音
信号との各インパルス応答と対応する信号の大きさ、つ
まり振幅又はパワーから各チャネルの更新ステップサイ
ズを生成し、各チャネルごとにその更新ステップサイズ
と、各受話信号と、上記残留反響信号とからそのチャネ
ルの反響路のインパルス応答を逐次推定し、その推定し
たインパルス応答によりそのチャネルの疑似反響路を生
成する。
ネルの受話信号をチャネル疑似反響路に印加してチャネ
ル疑似反響信号をそれぞれ生成し、これらチャネル疑似
反響信号を加算して総合疑似反響信号を生成し、この総
合疑似反響信号を集音信号から差し引いて残留反響信号
を生成し、特に各チャネルの受話信号の再生信号と集音
信号との各インパルス応答と対応する信号の大きさ、つ
まり振幅又はパワーから各チャネルの更新ステップサイ
ズを生成し、各チャネルごとにその更新ステップサイズ
と、各受話信号と、上記残留反響信号とからそのチャネ
ルの反響路のインパルス応答を逐次推定し、その推定し
たインパルス応答によりそのチャネルの疑似反響路を生
成する。
【００１１】
【発明の実施の形態】図１にこの発明の実施形態を示
す。説明を簡単にするため、受話２チャネル、集音１チ
ャネルを例とし、図７と対応する部分に同一参照符号を
付けてある。この実施形態では、各チャネルの反響路１
５1 ，１５2 のインパルス応答ｈ1 （ｋ），ｈ
2 （ｋ）と対応する信号をステップサイズ生成部３１に
入力して、各チャネルごとのステップサイズμ1 ，μ2
をそのインパルス応答に対応する信号の大きさに応じて
決定する。インパルス応答ｈ1 （ｋ），ｈ2 （ｋ）
と対応する信号としては受話信号ｘ1 （ｋ），ｘ
2 （ｋ）の各再生信号と集音信号ｙ（ｋ）との間の各イ
ンパルス応答、つまり反響路１５1 ，１５2 の各インパ
ルス応答ｈ1 （ｋ），ｈ2 （ｋ）それ自体、あるい
はスピーカ１２1 ，１２2 、マイクロホン１６1 の配置
に応じて経験的に推測される大略の値、またはチャネル
疑似反響路１８1 ，１８2 に設定される推定インパルス
応答などが用いられる。インパルス応答ｈ1 （ｋ），
ｈ2 （ｋ）を用いる場合は、予め、各１つのスピーカ
のみからインパルス音を放射して、マイクロホン１６の
出力信号を測定して求める。
す。説明を簡単にするため、受話２チャネル、集音１チ
ャネルを例とし、図７と対応する部分に同一参照符号を
付けてある。この実施形態では、各チャネルの反響路１
５1 ，１５2 のインパルス応答ｈ1 （ｋ），ｈ
2 （ｋ）と対応する信号をステップサイズ生成部３１に
入力して、各チャネルごとのステップサイズμ1 ，μ2
をそのインパルス応答に対応する信号の大きさに応じて
決定する。インパルス応答ｈ1 （ｋ），ｈ2 （ｋ）
と対応する信号としては受話信号ｘ1 （ｋ），ｘ
2 （ｋ）の各再生信号と集音信号ｙ（ｋ）との間の各イ
ンパルス応答、つまり反響路１５1 ，１５2 の各インパ
ルス応答ｈ1 （ｋ），ｈ2 （ｋ）それ自体、あるい
はスピーカ１２1 ，１２2 、マイクロホン１６1 の配置
に応じて経験的に推測される大略の値、またはチャネル
疑似反響路１８1 ，１８2 に設定される推定インパルス
応答などが用いられる。インパルス応答ｈ1 （ｋ），
ｈ2 （ｋ）を用いる場合は、予め、各１つのスピーカ
のみからインパルス音を放射して、マイクロホン１６の
出力信号を測定して求める。
【００１２】更新ステップサイズμ1 ，μ2 は各チャネ
ルのインパルス応答ｈ1 （ｋ），ｈ2 （ｋ）と対応
する信号（以下ｈ1 （ｋ），ｈ2 （ｋ）で代表す
る）の大きさ、つまり振幅又はパワーに応じ、大きけれ
ば大きな値とする。特に各チャネルのインパルス応答
ｈ1 （ｋ），ｈ2 （ｋ）のパワーの、これらパワー
の和に対する比に応じて更新ステップサイズμ1 ，μ2
を決定する。この場合、これら各チャネルの更新ステッ
プサイズμ1 2とμ2 2の和が２となるようにする。例えば
次のようにして求める。いまｈ1 （ｋ），ｈ
2 （ｋ）がそれぞれ１要素のみとして、これらをｈ平
面上に書くと図２に示すように、式（１４）の更新ベク
トルｈ＾（ｋ）は、ｈ1 （ｋ）とｈ2 （ｋ）のパ
ワーが等しい場合は、横軸に対し４５度のベクトル４１
上にあることになる。ベクトルｈ1 （ｋ），ｈ
2 （ｋ）の各長さをａ＝√（ｈ1 T （ｋ）ｈ
1 （ｋ）），ｂ＝√（ｈ2 T （ｋ）ｈ2 （ｋ））と
すると、更新ベクトルｈ＾（ｋ）は図２中のベクトル
４２上のように設定することで、その疑似反響路が真の
反響路にほぼ収束するのを速くすることができる。この
更新ベクトルｈ＾（ｋ）に応じてステップサイズ
μ1 ，μ2 を、√（ａ2 ）と√（ｂ2 ）の比と等しくす
る、つまり μ1 ：μ2 ＝√（ａ2 ）：√（ｂ2 ） とすればよい。従って下記の連立方程式が得られる。
ルのインパルス応答ｈ1 （ｋ），ｈ2 （ｋ）と対応
する信号（以下ｈ1 （ｋ），ｈ2 （ｋ）で代表す
る）の大きさ、つまり振幅又はパワーに応じ、大きけれ
ば大きな値とする。特に各チャネルのインパルス応答
ｈ1 （ｋ），ｈ2 （ｋ）のパワーの、これらパワー
の和に対する比に応じて更新ステップサイズμ1 ，μ2
を決定する。この場合、これら各チャネルの更新ステッ
プサイズμ1 2とμ2 2の和が２となるようにする。例えば
次のようにして求める。いまｈ1 （ｋ），ｈ
2 （ｋ）がそれぞれ１要素のみとして、これらをｈ平
面上に書くと図２に示すように、式（１４）の更新ベク
トルｈ＾（ｋ）は、ｈ1 （ｋ）とｈ2 （ｋ）のパ
ワーが等しい場合は、横軸に対し４５度のベクトル４１
上にあることになる。ベクトルｈ1 （ｋ），ｈ
2 （ｋ）の各長さをａ＝√（ｈ1 T （ｋ）ｈ
1 （ｋ）），ｂ＝√（ｈ2 T （ｋ）ｈ2 （ｋ））と
すると、更新ベクトルｈ＾（ｋ）は図２中のベクトル
４２上のように設定することで、その疑似反響路が真の
反響路にほぼ収束するのを速くすることができる。この
更新ベクトルｈ＾（ｋ）に応じてステップサイズ
μ1 ，μ2 を、√（ａ2 ）と√（ｂ2 ）の比と等しくす
る、つまり μ1 ：μ2 ＝√（ａ2 ）：√（ｂ2 ） とすればよい。従って下記の連立方程式が得られる。
【００１３】 μ1 √（ｂ2 ）＝μ2 √（ａ2 ） μ1 2＋μ2 2＝２ この連立方程式を解くと以下のようになる。 μ1 ＝２√（ａ2 ）／（√（ａ2 ）＋√（ｂ2 ）） μ2 ＝２√（ｂ2 ）／（√（ａ2 ）＋√（ｂ2 ）） つまり μ1 ＝２√（ｈ1 T （ｋ）ｈ1 （ｋ））／（√（ｈ1 T （ｋ）ｈ1 （ｋ））＋√（ｈ2 T （ｋ）ｈ2 （ｋ））） （17） μ2 ＝２√（ｈ2 T （ｋ）ｈ2 （ｋ））／√（ｈ1 T （ｋ）ｈ1 （ ｋ））＋√（ｈ2 T （ｋ）ｈ2 （ｋ））） （18） をステップサイズ生成部３１で計算すればよい。
【００１４】従ってインパルス応答ｈ1 （ｋ），ｈ
2 （ｋ）がそれぞれ２乗和回路３２ 1 ，３２2 に入力さ
れて、各ベクトルの要素の２乗和 ｈ1 T （ｋ）ｈ1 （ｋ）＝Σi=0 L-1 ｈ1 2（ｋ，ｉ） ｈ2 T （ｋ）ｈ2 （ｋ）＝Σi=0 L-1 ｈ2 2（ｋ，ｉ） がそれぞれ計算され、これら２乗和は開平回路３３1 ，
３３2 で平方根がとられ、ａ＝√（ｈ1 T （ｋ）ｈ
1 （ｋ），ｂ＝√（ｈ2 T （ｋ）ｈ2 （ｋ））が求
められる。これらａとｂが加算回路３４で加算され、式
（１７），（１８）の分母が求められる。乗除算回路３
５1 で加算回路３４の加算値でａを除算した値に２を乗
算して式（１７）の計算結果μ1 が求められる。同様に
乗除算回路３５2 で加算回路３４の加算値でｂを除算し
た値に２を乗算して式（１８）の計算結果μ2 が求めら
れる。
2 （ｋ）がそれぞれ２乗和回路３２ 1 ，３２2 に入力さ
れて、各ベクトルの要素の２乗和 ｈ1 T （ｋ）ｈ1 （ｋ）＝Σi=0 L-1 ｈ1 2（ｋ，ｉ） ｈ2 T （ｋ）ｈ2 （ｋ）＝Σi=0 L-1 ｈ2 2（ｋ，ｉ） がそれぞれ計算され、これら２乗和は開平回路３３1 ，
３３2 で平方根がとられ、ａ＝√（ｈ1 T （ｋ）ｈ
1 （ｋ），ｂ＝√（ｈ2 T （ｋ）ｈ2 （ｋ））が求
められる。これらａとｂが加算回路３４で加算され、式
（１７），（１８）の分母が求められる。乗除算回路３
５1 で加算回路３４の加算値でａを除算した値に２を乗
算して式（１７）の計算結果μ1 が求められる。同様に
乗除算回路３５2 で加算回路３４の加算値でｂを除算し
た値に２を乗算して式（１８）の計算結果μ2 が求めら
れる。
【００１５】これら更新ステップサイズμ1 ，μ2 がそ
れぞれ対応するチャネル反響推定部１９1 ，１９2 へ供
給される。チャネル反響路推定部１９1 ，１９2 には図
６に示した場合と同様にそれぞれ受話信号ｘ1 （ｋ），
ｘ2 （ｋ）が供給されると共にパワー計算部２７よりの
所定時間の全チャネルの受話信号パワーの和と残留反響
信号ｅ（ｋ）が共通に供給される。チャネル反響路推定
部１９1 ，１９2 ではそれぞれ以下の式によりチャネル
反響路のインパルス応答の推定ｈ＾1 （ｋ＋１），
ｈ＾2 （ｋ＋１）が逐次行われる。
れぞれ対応するチャネル反響推定部１９1 ，１９2 へ供
給される。チャネル反響路推定部１９1 ，１９2 には図
６に示した場合と同様にそれぞれ受話信号ｘ1 （ｋ），
ｘ2 （ｋ）が供給されると共にパワー計算部２７よりの
所定時間の全チャネルの受話信号パワーの和と残留反響
信号ｅ（ｋ）が共通に供給される。チャネル反響路推定
部１９1 ，１９2 ではそれぞれ以下の式によりチャネル
反響路のインパルス応答の推定ｈ＾1 （ｋ＋１），
ｈ＾2 （ｋ＋１）が逐次行われる。
【００１６】 ｈ＾1 （ｋ＋１）＝ｈ＾1 （ｋ）＋（μ1 ｅ（ｋ）／（ＸT （ｋ）Ｘ （ｋ）））Ｘ1 （ｋ） （19） ｈ＾2 （ｋ＋２）＝ｈ＾2 （ｋ）＋（μ2 ｅ（ｋ）／（ＸT （ｋ）Ｘ （ｋ）））Ｘ1 （ｋ） （20） これら推定されたインパルス応答ｈ＾1 （ｋ＋１），
ｈ＾2 （ｋ＋１）がそれぞれチャネル疑似反響路１８
1 ，１８2 に設定される。その他は図６に示した装置と
同様な処理を行う。
ｈ＾2 （ｋ＋１）がそれぞれチャネル疑似反響路１８
1 ，１８2 に設定される。その他は図６に示した装置と
同様な処理を行う。
【００１７】このように更新ステップサイズμ1 ，μ2
が反響路１５1 ，１５2 のインパルス応答ｈ
1 （ｋ），ｈ2 （ｋ）のエネルギーの相対比に応じ
て、大きいエネルギーのチャネルに対応する更新ステッ
プサイズが大きくされ、それだけ全体としての疑似反響
路１８1 ，１８2 の収束が速くなる。ステップサイズ生
成部３１で各時刻ｋごとにμ1 （ｋ），μ2 （ｋ）を計
算することにより、反響路１５1 ，１５2 が変動して
ｈ1 （ｋ），ｈ2 （ｋ）が変動してもこれに追従し
て、適切なμ1 （ｋ），μ2 （ｋ）が設定されることに
なる。また最初にステップサイズ生成部３１に、インパ
ルス応答ｈ1 （ｋ），ｈ2 （ｋ）の経験的推定値を
入力しても、その後、チャネル疑似反響路１８1，１８
2 の各インパルス応答ｈ＾1 （ｋ），ｈ＾2 （ｋ）
をステップサイズ生成部３１に入力することにより、チ
ャネル疑似反響路１８1 ，１８2 が反響路１５1 ，１５
2 に収束するに従って、更新ステップサイズμ1 ，μ2
が適切な値μ1 （ｋ），μ2 （ｋ）に更新される。ある
いはμ1 ，μ2 の初期値はμ1 ＝μ 2 とし従来において
使用していた適当な値を設定し、その後は、ｈ＾
1 （ｋ），ｈ＾2 （ｋ）をステップサイズ生成部３１
に入力してμ1 ，μ2 を更新してもよい。
が反響路１５1 ，１５2 のインパルス応答ｈ
1 （ｋ），ｈ2 （ｋ）のエネルギーの相対比に応じ
て、大きいエネルギーのチャネルに対応する更新ステッ
プサイズが大きくされ、それだけ全体としての疑似反響
路１８1 ，１８2 の収束が速くなる。ステップサイズ生
成部３１で各時刻ｋごとにμ1 （ｋ），μ2 （ｋ）を計
算することにより、反響路１５1 ，１５2 が変動して
ｈ1 （ｋ），ｈ2 （ｋ）が変動してもこれに追従し
て、適切なμ1 （ｋ），μ2 （ｋ）が設定されることに
なる。また最初にステップサイズ生成部３１に、インパ
ルス応答ｈ1 （ｋ），ｈ2 （ｋ）の経験的推定値を
入力しても、その後、チャネル疑似反響路１８1，１８
2 の各インパルス応答ｈ＾1 （ｋ），ｈ＾2 （ｋ）
をステップサイズ生成部３１に入力することにより、チ
ャネル疑似反響路１８1 ，１８2 が反響路１５1 ，１５
2 に収束するに従って、更新ステップサイズμ1 ，μ2
が適切な値μ1 （ｋ），μ2 （ｋ）に更新される。ある
いはμ1 ，μ2 の初期値はμ1 ＝μ 2 とし従来において
使用していた適当な値を設定し、その後は、ｈ＾
1 （ｋ），ｈ＾2 （ｋ）をステップサイズ生成部３１
に入力してμ1 ，μ2 を更新してもよい。
【００１８】この発明は受話チャネルが３チャネル以上
の場合にも適用できる。例えばＮチャネルの場合、ステ
ップサイズ生成部３１で各チャネルｎ（ｎ＝１，２，・・
・，Ｎ）に対し、次の計算をすればよい。 μｎ＝Ｎ√（ｈn T（ｋ）ｈn（ｋ））／√（ｈ1 T
（ｋ）ｈ1（ｋ）＋(ｈ2 T（ｋ）ｈ2（ｋ）・・・＋
ｈN T（ｋ）ｈN（ｋ）） この式から理解されるように、各チャネルのインパルス
応答の大きさ（振幅又はパワー）の、チャネル間の相対
的大きさ比により、インパルス応答の大きさが大きいチ
ャネルの疑似反響路の収束速度が大きくなるようにステ
ップサイズを求めればよい。また上記実施形態ではチャ
ネル反響路推定部において、全チャネルの受話信号のパ
ワーの和を求めたものを用いたが、必ずしもその必要は
ない。例えば電子情報通信学会論文誌’８６／１０ Ｖ
ｏｌ．Ｊ６９−Ａ Ｎｏ．１０「多チャネル適応ディジ
タルフィルタ」で提案しているように、チャネルごとに
その受話信号をその電力で正規化したものを用いてもよ
い。
の場合にも適用できる。例えばＮチャネルの場合、ステ
ップサイズ生成部３１で各チャネルｎ（ｎ＝１，２，・・
・，Ｎ）に対し、次の計算をすればよい。 μｎ＝Ｎ√（ｈn T（ｋ）ｈn（ｋ））／√（ｈ1 T
（ｋ）ｈ1（ｋ）＋(ｈ2 T（ｋ）ｈ2（ｋ）・・・＋
ｈN T（ｋ）ｈN（ｋ）） この式から理解されるように、各チャネルのインパルス
応答の大きさ（振幅又はパワー）の、チャネル間の相対
的大きさ比により、インパルス応答の大きさが大きいチ
ャネルの疑似反響路の収束速度が大きくなるようにステ
ップサイズを求めればよい。また上記実施形態ではチャ
ネル反響路推定部において、全チャネルの受話信号のパ
ワーの和を求めたものを用いたが、必ずしもその必要は
ない。例えば電子情報通信学会論文誌’８６／１０ Ｖ
ｏｌ．Ｊ６９−Ａ Ｎｏ．１０「多チャネル適応ディジ
タルフィルタ」で提案しているように、チャネルごとに
その受話信号をその電力で正規化したものを用いてもよ
い。
【００１９】またチャネル反響路推定のアルゴリズムと
しては学習同定法に限らず、射影法その他のアルゴリズ
ムを用いてもよい。更にこの多チャネルエコーキャンセ
ル装置を、コンピュータによりプログラムを実行させて
機能させることもできる。この場合はこの多チャネルエ
コーキャンセルプログラムを、ＣＤ−ＲＯＭ、フロッピ
ー（登録商標）ディスク、磁気ディスクなどからコンピ
ュータのプログラムメモリにインストールし、あるいは
通信回線を通じてダウンロードして、そのプログラムを
実行することになる。
しては学習同定法に限らず、射影法その他のアルゴリズ
ムを用いてもよい。更にこの多チャネルエコーキャンセ
ル装置を、コンピュータによりプログラムを実行させて
機能させることもできる。この場合はこの多チャネルエ
コーキャンセルプログラムを、ＣＤ−ＲＯＭ、フロッピ
ー（登録商標）ディスク、磁気ディスクなどからコンピ
ュータのプログラムメモリにインストールし、あるいは
通信回線を通じてダウンロードして、そのプログラムを
実行することになる。
【００２０】
【発明の効果】図１に示した実施形態において、マイク
ロホン１６1 と正面スピーカ１２1 間の距離が約１．５
ｍ、スピーカ１２1 と１２2 間の距離が約１．２ｍの場
合、インパルス応答ｈ1 （ｋ）とｈ2 （ｋ）のパワ
ー比は壁の反射などにより変わり、１：０．７５〜０．
８５程度となるが、１：０．８とした。この場合の計算
機シミュレーションの結果を図３に示す。図３中の実線
はこの発明方法による疑似反響路の収束状態を示し、破
線は従来法による収束状態を示す。反響路を変化させて
から係数誤差（インパルス応答推定値誤差）が−１０ｄ
Ｂになる時間が、従来法では８．３秒掛ったが、この発
明方法によれば５．９秒となり、従来法より約１．５倍
の収束速度が得られ、この発明方法が優れていることが
確認できた。
ロホン１６1 と正面スピーカ１２1 間の距離が約１．５
ｍ、スピーカ１２1 と１２2 間の距離が約１．２ｍの場
合、インパルス応答ｈ1 （ｋ）とｈ2 （ｋ）のパワ
ー比は壁の反射などにより変わり、１：０．７５〜０．
８５程度となるが、１：０．８とした。この場合の計算
機シミュレーションの結果を図３に示す。図３中の実線
はこの発明方法による疑似反響路の収束状態を示し、破
線は従来法による収束状態を示す。反響路を変化させて
から係数誤差（インパルス応答推定値誤差）が−１０ｄ
Ｂになる時間が、従来法では８．３秒掛ったが、この発
明方法によれば５．９秒となり、従来法より約１．５倍
の収束速度が得られ、この発明方法が優れていることが
確認できた。
【００２１】以上述べたようにこの発明によれば、各チ
ャネルの反響路のインパルス応答に応じてそのチャネル
の反響路推定における更新ステップサイズを決定するこ
とにより、多チャネルエコーキャンセルの収束速度を向
上させることができる。
ャネルの反響路のインパルス応答に応じてそのチャネル
の反響路推定における更新ステップサイズを決定するこ
とにより、多チャネルエコーキャンセルの収束速度を向
上させることができる。
【図１】この発明の実施形態の機能構成を示す図。
【図２】インパルス応答に応じてステップサイズを決め
るための説明図。
るための説明図。
【図３】この発明の効果を確認するための、収束状態の
計算機シミュレーションを示す図。
計算機シミュレーションを示す図。
【図４】従来の１チャネルエコーキャンセル装置の機能
構成を示す図。
構成を示す図。
【図５】従来の多チャネルエコーキャンセル装置の構成
図。
図。
【図６】従来の多チャネルエコーキャンセル装置の具体
的機能構成を示す図。
的機能構成を示す図。
【図７】従来の２チャネルエコーキャンセル装置の機能
構成を示す図。
構成を示す図。
Claims (6)
- 【請求項１】 Ｎチャネル（Ｎは２以上の整数）の受話
信号をそれぞれ対応するチャネル疑似反響路に印加して
チャネル疑似反響信号を生成し、 これらチャネル疑似反響信号を加算して総合疑似反響信
号を生成し、 その総合疑似反響信号を、上記各受話信号の再生信号に
対する集音信号から差し引いて残留反響信号を生成し、 上記各チャネルの受話信号の再生信号と上記集音信号と
の間の各インパルス応答の大きさから各チャネルの更新
ステップサイズを算出し、 上記各チャネルの更新ステップサイズと、各チャネルの
受話信号と、上記残留反響信号とからそのチャネルの上
記再生信号の反響路のインパルス応答を逐次推定し、 その逐次推定したインパルス応答により対応するチャネ
ルの疑似反響路を生成することを特徴とする多チャネル
エコーキャンセル方法。 - 【請求項２】 上記更新ステップサイズの算出は、各チ
ャネルのインパルス応答対応パワーの、チャネル間の相
対比を求め、この比に基づきそのチャネルの上記更新ス
テップサイズを決定することを特徴とする請求項１記載
の多チャネルエコーキャンセル方法。 - 【請求項３】 Ｎチャネル（Ｎは２以上の整数）の受話
信号がそれぞれ印加されて、チャネル疑似反響信号を生
成するＮ個のチャネル疑似反響路と、 上記Ｎ個のチャネル疑似反響信号を加算して総合疑似信
号を生成する合成回路と、 上記Ｎ個の受話信号の再生信号に対する集音信号から上
記総合疑似信号を差し引いて残留反響信号を生成する減
算回路と、 上記Ｎ個の受話信号の再生信号と上記集音信号との間の
各インパルス応答を入力して、その大きさに応じたＮ個
の更新ステップサイズを求めるステップサイズ生成部
と、 上記各チャネルごとに上記更新ステップサイズと、各チ
ャネルの受話信号と、上記残留反響信号とが入力され、
そのチャネルの上記再生信号の反響路のインパルス応答
を逐次推定し、その推定したインパルス応答を対応する
上記チャネル疑似反響路に設定するＮ個の反響路推定部
と、を具備する多チャネルエコーキャンセル装置。 - 【請求項４】 上記ステップサイズ生成部は各インパル
ス応答対応のパワーの、チャネル間の相対比に基づき各
対応チャネルの更新ステップサイズを決定して出力する
回路であることを特徴とする請求項３記載の多チャネル
エコーキャンセル装置。 - 【請求項５】 請求項１又は２記載の多チャネルエコー
キャンセル方法をコンピュータに実行させるプログラ
ム。 - 【請求項６】 請求項５記載のプログラムを記録したコ
ンピュータ読み出し可能な記録媒体。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2001051387A JP3616341B2 (ja) | 2001-02-27 | 2001-02-27 | 多チャネルエコーキャンセル方法、その装置、そのプログラム及び記録媒体 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2001051387A JP3616341B2 (ja) | 2001-02-27 | 2001-02-27 | 多チャネルエコーキャンセル方法、その装置、そのプログラム及び記録媒体 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2002261659A true JP2002261659A (ja) | 2002-09-13 |
| JP3616341B2 JP3616341B2 (ja) | 2005-02-02 |
Family
ID=18912179
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2001051387A Expired - Fee Related JP3616341B2 (ja) | 2001-02-27 | 2001-02-27 | 多チャネルエコーキャンセル方法、その装置、そのプログラム及び記録媒体 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3616341B2 (ja) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2007067854A (ja) * | 2005-08-31 | 2007-03-15 | Nippon Telegr & Teleph Corp <Ntt> | 反響消去方法、反響消去装置、プログラム、記録媒体 |
| JP2007067855A (ja) * | 2005-08-31 | 2007-03-15 | Nippon Telegr & Teleph Corp <Ntt> | 反響消去方法、反響消去装置、プログラム、記録媒体 |
| JP2009284096A (ja) * | 2008-05-20 | 2009-12-03 | Nippon Telegr & Teleph Corp <Ntt> | ステレオ音響エコーキャンセル方法、ステレオ音響エコーキャンセル装置、ステレオ音響エコーキャンセルプログラム、その記録媒体 |
| JP2011114696A (ja) * | 2009-11-27 | 2011-06-09 | Nippon Telegr & Teleph Corp <Ntt> | エコーキャンセル方法、エコーキャンセル装置、エコーキャンセルプログラム |
| JP2015028639A (ja) * | 2008-11-20 | 2015-02-12 | アップル インコーポレイテッド | 音声信号補償を用いるアクティブノイズコントロールのためのシステム |
-
2001
- 2001-02-27 JP JP2001051387A patent/JP3616341B2/ja not_active Expired - Fee Related
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2007067854A (ja) * | 2005-08-31 | 2007-03-15 | Nippon Telegr & Teleph Corp <Ntt> | 反響消去方法、反響消去装置、プログラム、記録媒体 |
| JP2007067855A (ja) * | 2005-08-31 | 2007-03-15 | Nippon Telegr & Teleph Corp <Ntt> | 反響消去方法、反響消去装置、プログラム、記録媒体 |
| JP4504891B2 (ja) * | 2005-08-31 | 2010-07-14 | 日本電信電話株式会社 | 反響消去方法、反響消去装置、プログラム、記録媒体 |
| JP4504892B2 (ja) * | 2005-08-31 | 2010-07-14 | 日本電信電話株式会社 | 反響消去方法、反響消去装置、プログラム、記録媒体 |
| JP2009284096A (ja) * | 2008-05-20 | 2009-12-03 | Nippon Telegr & Teleph Corp <Ntt> | ステレオ音響エコーキャンセル方法、ステレオ音響エコーキャンセル装置、ステレオ音響エコーキャンセルプログラム、その記録媒体 |
| JP2015028639A (ja) * | 2008-11-20 | 2015-02-12 | アップル インコーポレイテッド | 音声信号補償を用いるアクティブノイズコントロールのためのシステム |
| JP2011114696A (ja) * | 2009-11-27 | 2011-06-09 | Nippon Telegr & Teleph Corp <Ntt> | エコーキャンセル方法、エコーキャンセル装置、エコーキャンセルプログラム |
Also Published As
| Publication number | Publication date |
|---|---|
| JP3616341B2 (ja) | 2005-02-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103428385B (zh) | 用于处理音频信号的方法及用于处理音频信号的电路布置 | |
| CN104158990B (zh) | 用于处理音频信号的方法和音频接收电路 | |
| US8189765B2 (en) | Multichannel echo canceller | |
| JP4654777B2 (ja) | 音響エコーキャンセル装置 | |
| EP2749016B1 (en) | Processing audio signals | |
| US5774562A (en) | Method and apparatus for dereverberation | |
| JPH08181639A (ja) | 多チャネル音声通信会議用反響消去方法 | |
| US9667803B2 (en) | Nonlinear acoustic echo cancellation based on transducer impedance | |
| EP3829151B1 (en) | Echo suppression device, echo suppression method, and echo suppression program | |
| JP2004349806A (ja) | 多チャネル音響エコー消去方法、その装置、そのプログラム及びその記録媒体 | |
| US6381272B1 (en) | Multi-channel adaptive filtering | |
| JP2003309493A (ja) | 反響低減方法、反響低減装置、反響低減プログラム | |
| JP3616341B2 (ja) | 多チャネルエコーキャンセル方法、その装置、そのプログラム及び記録媒体 | |
| JP3673727B2 (ja) | 反響消去方法、その装置、そのプログラム及びその記録媒体 | |
| JP2003250193A (ja) | 反響消去方法、この方法を実施する装置、プログラムおよびその記録媒体 | |
| JP4041770B2 (ja) | 音響エコー消去方法、その装置、プログラム及びその記録媒体 | |
| JP3628267B2 (ja) | 多チャネル反響消去方法、その装置、そのプログラム及びその記録媒体 | |
| KR102045953B1 (ko) | 칼만필터 기반의 다채널 입출력 음향학적 반향 제거 방법 | |
| JP2002062900A (ja) | 収音装置及び受信装置 | |
| JP3583998B2 (ja) | 多チャネルエコー消去方法、その装置及びプログラム記録媒体 | |
| JP6075783B2 (ja) | エコー消去装置、エコー消去方法及びプログラム | |
| JP4663630B2 (ja) | マルチチャンネルシステム同定装置 | |
| JP6356087B2 (ja) | エコー消去装置、その方法及びプログラム | |
| JP3339612B2 (ja) | 多チャンネル音声通信会議用反響消去方法 | |
| JP3700838B2 (ja) | 多チャネル反響消去方法・多チャネル反響消去装置・多チャネル反響消去プログラム |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A977 | Report on retrieval |
Effective date: 20040601 Free format text: JAPANESE INTERMEDIATE CODE: A971007 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20040713 |
|
| A521 | Written amendment |
Effective date: 20040910 Free format text: JAPANESE INTERMEDIATE CODE: A523 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Effective date: 20041012 Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
| A61 | First payment of annual fees (during grant procedure) |
Effective date: 20041104 Free format text: JAPANESE INTERMEDIATE CODE: A61 |
|
| FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20071112 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20081112 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20091112 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20101112 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20101112 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20111112 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (prs date is renewal date of database) |
Year of fee payment: 7 Free format text: PAYMENT UNTIL: 20111112 |
|
| FPAY | Renewal fee payment (prs date is renewal date of database) |
Year of fee payment: 8 Free format text: PAYMENT UNTIL: 20121112 |
|
| FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20121112 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (prs date is renewal date of database) |
Year of fee payment: 9 Free format text: PAYMENT UNTIL: 20131112 |
|
| LAPS | Cancellation because of no payment of annual fees |