JPH08181639A

JPH08181639A - 多チャネル音声通信会議用反響消去方法

Info

Publication number: JPH08181639A
Application number: JP7050002A
Authority: JP
Inventors: Suehiro Shimauchi; 末廣島内; Shoji Makino; 昭二牧野; Junji Kojima; 順治小島
Original assignee: Nippon Telegraph and Telephone Corp
Current assignee: Nippon Telegraph and Telephone Corp
Priority date: 1994-10-26
Filing date: 1995-03-09
Publication date: 1996-07-12
Anticipated expiration: 2018-04-14
Also published as: EP0709999B1; CA2161358A1; US5661813A; CA2161358C; DE69531136D1; EP0709999A3; JP3397269B2; DE69531136T2; EP0709999A2

Abstract

(57)【要約】【目的】反響路の推定を高速かつ正しく行う。【構成】この発明の第１の方法は、各チャネル間の過
去の受話信号の相互相関に対する現在の受話信号間の相
互相関の変動成分を抽出し、この変動成分を修正ベクト
ルとし、この修正ベクトルを用いて各反響路のインパル
ス応答の推定を逐次修正する。この発明の第２の方法は
受話信号の相互相関を、能動的に変動させる機能を付加
し、この相互相関の変動が付加された受話信号を疑似反
響路ベクトルの修正ベルトルの導出に利用する。

Description

【発明の詳細な説明】

【０００１】

【産業上の利用分野】この発明は、多チャネル再生系を
有する通信会議システムにおいて、ハウリングの原因お
よび聴覚上の障害となる室内反響信号を消去する多チャ
ネル音声通信会議用反響消去方法に関するものである。

【０００２】

【従来の技術】同時通話性能に優れ反響感の少ない拡声
通話システムを提供するために、反響消去装置がある。
まず、１チャネル用の反響消去装置について、その反響
消去方法および装置構成を図３を参照して説明する。拡
声通話において、受話端子１１に相手の発話等で得られ
る受話信号は、そのままスピーカ（再生器）１２から再
生される場合と、スピーカ１２へ送る前に、受話信号の
振幅やパワー等の大きさに応じて自動的に利得を調節す
る受話信号加工部１３が挿入され、受話信号に何らかの
加工が施された後に、スピーカ１２から再生される場合
とがある。このため、この明細書で受話信号ｘ₁(ｋ）と
は、相手からの受話信号そのものとは限らず、受話信号
加工部１３が挿入されたときは、加工された後の受話信
号を指すものとする。ｋは離散時間を表す。反響消去装
置１４は、受話信号ｘ₁(ｋ）がスピーカ１２から反響路
１５を経て、マイクロホン（集音器）１６に集音されて
得られる反響信号ｙ ₁(ｋ）を消去する。ここで、反響信
号ｙ₁(ｋ）は、時刻ｋにおける反響路１５のインパルス
応答をｈ₁₁（ｋ，ｎ）として、ｙ₁(ｋ）＝Σｈ₁₁（ｋ，ｎ）ｘ₁(ｋ−ｎ）（１）Σ はｎ＝０からＬ−１までのような畳み込み演算で得られるものとモデル化でき
る。Ｌはタップ数で、反響路１５の残響時間に対応させ
て、あらかじめ設定しておく定数である。まず、受話信
号蓄積・ベクトル生成部１７において、受話信号ｘ
₁(ｋ）をＬ−１時刻過去のものまで蓄積しておく。蓄積
された信号は、受話信号ベクトルｘ₁(ｋ）、すなわちｘ₁(ｋ）＝［ｘ₁(ｋ),ｘ₁(ｋ−１),…, ｘ₁(ｋ−Ｌ＋１)]^T （２）として出力される。但し、＊^Tはベクトルの転置を表
す。擬似反響信号生成部１８では、式（２）の受話信号
ベクトルｘ₁(ｋ）と、反響路推定部１９から得られる
擬似反響ベクトルｈ＾₁₁（ｋ）との内積演算ｙ＾₁(ｋ）＝ｈ＾₁₁ ^T（ｋ）ｘ₁(ｋ）（３）を行ない、その結果として、擬似反響信号ｙ＾₁(ｋ）を
生成する。この内積演算は、式（１）のような畳み込み
演算と等価である。反響路推定部１９では、擬似反響信
号生成部１８で用いる擬似反響路ベクトルｈ＾
₁₁（ｋ）を生成する。この反響路推定に用いる最も一般
的なアルゴリズムは、ＮＬＭＳアルゴリズム（学習同定
法）である。ＮＬＭＳアルゴリズムでは、時刻ｋにおけ
る受話信号ベクトルｘ₁(ｋ）と、残留反響信号ｅ
₁(ｋ）、すなわち差回路２１でマイクロホン１６の出力
ｙ₁(ｋ）から擬似反響信号ｙ＾₁(ｋ）を差し引いた。

【０００３】ｅ₁(ｋ）＝ｙ₁(ｋ）−ｙ＾₁(ｋ）（４）とから、時刻ｋ＋１において用いる擬似反響路ベクトル
ｈ＾₁₁（ｋ＋１）を次式のように求める。ｈ＾₁₁(k+1) ＝ｈ＾₁₁(k) ＋αｅ₁(k)ｘ₁(k)／（ｘ₁ ^T(k) ｘ₁(k)）（５）但し、αはステップサイズパラメータと呼ばれ０＜α＜
２の範囲で適応動作の調整に用いる。以上のような処理
を繰り返すことにより、反響路推定部１９では、次第に
擬似反響路ベクトルｈ＾₁₁（ｋ）を、真の反響路１５
のインパルス応答ｈ₁₁（ｋ，ｎ）の時系列を各要素とし
て持つ反響路ベクトルｈ₁₁（ｋ）、すなわちｈ₁₁(k) ＝［ｈ₁₁(k,0),ｈ₁₁(k,1),…, ｈ₁₁(k,L-1)]^T （６）と一致させることが可能となり、その結果式（４）の残
留反響信号ｅ₁ （ｋ）を小さくすることができる。

【０００４】一般にＮ（＞２）チャネルの再生系とＭ
（＞１）チャネルの集音系とで構成される通信会議シス
テムの場合の反響の消去は、図４に示すような構成によ
り行われる。すなわち再生側の全Ｎチャネルと集音側の
各１チャネルとの間にＮ入力１出力時系列信号を処理す
るＮチャネル反響消去装置２２₁,２２₂,…，２２_Mをそ
れぞれ接続した反響消去システム２３として実現され
る。この場合システム全体でＮ×Ｍ個の反響路１５
_nm（１＜ｎ＜Ｎ，１＜ｍ＜Ｍ）が存在する。このシステ
ムの構成単位である再生側の全Ｎチャネルと集音側の各
１チャネルとの間に接続されるＮチャネル反響消去装置
２２₁,２２₂,…，２２_Mについては、図３に示した反響
消去装置１４の構成を拡張して、図５に示すように構成
される。これは例えば電子情報通信学会論文誌 '８６／
１０Ｖｏｌ．Ｊ６９−ＡＮｏ．１０「多チャネル適
応ディジタルフィルタ」に詳しく述べられている。ここ
で、集音側が第ｍ集音チャネル（１＜ｍ＜Ｍ）に接続さ
れているＮチャネル反響消去装置２２_mを考える。第ｍ
チャネルの集音器１６_mで集音される反響信号は、各再
生チャネルの受話信号がそれぞれの反響路１５_1m〜１５
_Nmを経て集音側で全て加算されることにより得られるた
めに、反響路推定をどの再生チャネルについても、統一
的に同じ１つの残留反響信号ｅ_m（ｋ）のみを評価して
行なうための工夫が必要となる。まず、各再生チャネル
の受話信号について、受話信号蓄積・ベクトル生成部
（１７₁,１７₂,…，１７_N）により、受話信号ベクトルｘ₁(ｋ）＝［ｘ₁(ｋ),ｘ₁(ｋ−１),…, ｘ₁(ｋ−Ｌ₁ ＋１)]^T （７）ｘ₂(ｋ）＝［ｘ₂(ｋ),ｘ₂(ｋ−１),…, ｘ₂(ｋ−Ｌ₂ ＋１)]^T （８）：ｘ_N(k）＝［ｘ_N(k),ｘ_N(k−1), …, ｘ_N(k−Ｌ_N＋１)]^T （９）を生成する。但し、Ｌ₁,Ｌ₂,…，Ｌ_Nはタップ数で、各
反響路１５_1m，１５_2m，…，１５_Nmの残響時間に対応さ
せて、あらかじめ設定する定数である。これらのベクト
ルをベクトル結合部２４によって、ｘ(k) ＝［ｘ₁ ^T(k),ｘ₂ ^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（ｋ）を生成す
る。このように、各チャネル毎のベクトルを結合して１
つのベクトルとして扱うことにより、基本的な処理の流
れは、図３に示した１チャネル反響消去装置と同様とな
る。

【０００５】

【発明が解決しようとする課題】Ｎチャネルの再生系と
Ｍチャネルの集音系とで構成される通信会議システムに
用いられる従来の反響消去システムの欠点のうち、この
発明が解決しようとするものを、具体例に沿って説明す
る。図６に示すように、Ａ地点とＢ地点との間におい
て、２チャネルで集音・再生を行なうステレオ音声会議
装置に従来の反響消去システムを適用した場合、Ｂ地点
側の反響路が全く変動しない場合においても、Ａ地点側
で発話者が移動あるいは交替を行なう度にその発話に対
するＢ地点からの反響音が増大するという問題がある。
これは、Ｂ地点側の反響消去システムにおいて反響路推
定が正しく行なわれていないために起こる問題である。

【０００６】そこで、この問題を説明するために、Ｂ地
点における反響消去システムの構成単位である２チャネ
ル反響消去装置のうちの１つである第１集音チャネルに
接続された装置２２ｂ₁の動作に注目する。以下では、
各再生チャネルの動作を明確化するために、従来技術の
項で用いた結合ベクトルの各再生チャネルについての要
素が明示された形で表記する。結合ベクトルを用いずに
各再生チャネル毎のベクトルで表記するＢ地点側におけ
る２チャネル受話信号ベクトルをｘ₁(ｋ）とｘ
₂(ｋ）とする。各再生チャネルに対する真の反響路１５
₁₁，１５₂₁の反響路ベクトルをｈ₁₁（ｋ),ｈ
₂₁（ｋ) とすると、これら反響路１５₁₁，１５₂₁を経て
集音される反響信号ｙ₁(ｋ）は、ｙ₁(ｋ）＝ｈ₁₁ ^T（ｋ) ｘ₁(ｋ）＋ｈ₂₁ ^T（ｋ) ｘ₂(ｋ）（14）と表される。一方、反響消去装置内で生成される擬似反
響信号ｙ＾₁(ｋ）は、装置内で生成される擬似反響路
ｈ＾₁₁（ｋ) とｈ＾₂₁（ｋ) とを用いて、ｙ＾₁(k)＝ｈ＾₁₁ ^T(k) ｘ₁(ｋ）＋ｈ＾₂₁ ^T(k) ｘ₂(ｋ）（15）と得られる。Ａ地点から一人の話者が発話した場合に
は、この受話信号ベクトルｘ₁(ｋ）とｘ₂(ｋ）との
間には、非常に強い相互相関がある。ｘ₁(ｋ）とｘ
₂(ｋ）との間に、一定な強い相互相関がある場合には、ｙ＾₁(k)＝ｙ₁(ｋ）（16）の解としての結合ベクトル［ｈ＾₁₁ ^T(k),ｈ＾₂₁ ^T
（ｋ)]は無数に存在し、ｘ₁(ｋ）とｘ₂(ｋ）との相
互相関に固有な部分空間Ｈ_xを形成する。このため、
ＮＬＭＳアルゴリズムのような一般的な逐次誤差最小化
アルゴリズムを用いた場合に、［ｈ＾₁₁ ^T(k),ｈ＾
₂₁ ^T(k)]は初期値から部分空間Ｈ_xまでの距離が最小
となる点に収束し、一般に真値［ｈ₁₁ ^T(k),ｈ₂₁ ^T
(k)]には収束しない。

【０００７】説明を簡単にするために、一定なスカラー
値ｐ₁,ｐ₂ と原信号ベクトルｓ(k) により、受話信号
ベクトルｘ₁(k), ｘ₂(k)が、ｘ₁(ｋ）＝ｐ₁ ｓ(k) ，ｘ₂(ｋ）＝ｐ₂ ｓ(k) （17）と表される場合を考える。［ｈ＾₁₁ ^T(k),ｈ＾₂₁ ^T
(k)]が存在し得る部分空間Ｈ_xは、 p₁ｈ＾₁₁(k)+p₂ｈ＾₂₁(k) ＝p₁ｈ₁₁(k)+p₂ｈ₂₁(k) （18）を満たす図７Ａ中の直線のように見なすことができ、初
期値０から適応を開始した場合、収束点［ｈ＾_11p ^T
(k),ｈ＾_21p ^T(k)]は、ｈ＾_11p(k)=p₁ ²(ｈ₁₁(k)+p₂ｈ₂₁(k)/p₁)/(p₁ ²+p₂ ²) ≠ｈ₁₁(k) （19）ｈ＾_21p(k)=p₂ ²(p₁ｈ₁₁(k)/p₂＋ｈ₂₁(k))/(p₁ ²+p₂ ²)≠ｈ₂₁(k) （20）のように得られる。このため、スカラー値ｐ₁ とｐ₂ と
の間の比が変動した時点で、式（１８）は満足されなく
なるため、反響を消去できなくなり、反響が増大する。

【０００８】以上の例からも分かるように、一般にＮチ
ャネルの再生系とＭチャネルの集音系とで構成される通
信会議システムに、従来の反響消去装置を適用した場合
は、各チャネルの受話信号間に相互相関がある場合に、
反響路の推定が正しく行なわれないため、受話信号間の
相互相関が変動する度に反響が増大する問題が生じる。

【０００９】この発明の目的は、従来の多チャネル音声
通信会議用反響消去システムにおいて、以上に挙げたよ
うな、問題を解決する多チャネル音声通信会議用の反響
消去方法を提供することにある。

【００１０】

【課題を解決するための手段】この発明による反響消去
方法には第１の方法と、第２の方法とがある。第１の方
法としては、残留反響信号と各チャネルの受話信号とか
ら修正ベクトルを求め、この修正ベクトルを用いて、各
反響路のインパルス応答の推定を逐次修正するが、この
発明では各チャネル間の過去の受話信号の相互相関に対
する現在の受話信号間の相互相関の変動成分を抽出し、
この変動成分を修正ベクトルとする。

【００１１】第２の方法としては、受話信号間の相互相
関を、聴感上の違和感のない範囲で能動的に変動させる
機能を付加して、各再生器で音響信号を再生し、この相
互相関の変動が付加された受話信号を疑似反響路ベクト
ルの修正ベクトルの導出に利用する。以下に第１の方法
を具体的に説明する。このためまず、多チャネル受話信
号間の相互相関が変動した場合における、多チャネル反
響消去装置の反響路推定部の動作について考察する。

【００１２】簡単のため、式（１７）で表される２チャ
ネルステレオ受話信号ｘ₁ （ｋ）とｘ₂ （ｋ）との
相互相関が変動した場合として、式（１７）のｐ₁ ，ｐ
₂ がそれらとは比関係の異なるｑ₁ ，ｑ₂ に変化した例
を考える。まず、〔ｈ＾₁₁ ^T（ｋ），ｈ＾
₂₁ ^T（ｋ）〕は、式（１９），（２０）の〔ｈ＾_11p
^T（ｋ），ｈ＾_21p ^T（ｋ）〕に収束する。次に、ｐ
₁ とｐ₂ とがｑ₁ とｑ₂ に変化したとき、〔ｈ＾_11p
^T（ｋ），ｈ＾_21p ^T（ｋ）〕を“初期値”として、
そこから最短距離にある〔ｈ＾_11q ^T（ｋ），ｈ＾
_21q ^T（ｋ）〕に収束する。この動作は、図７Ｂに示す
ように幾何学的に解釈される。ここで、収束点（ｈ＾
_11q，ｈ＾_21q）は収束点（ｈ＾_11p，ｈ＾_21p）から
直線ｈ＾₁₁ｑ₁ ＋ｈ＾₂₁ｑ₂ ＝ｈ₁₁ｑ₁ ＋ｈ₂₁ｑ₂ にお
ろした垂線の交点になる。従って図７Ｂにおけるフィル
タ係数誤差ベクトルｅ_p，ｅ_qのノルムは、一般に‖ｅ
_p‖＞‖ｅ_q‖の関係にあることが自明である。このこ
とは、一般に、受話信号間の相互相関の変化の度に、フ
ィルタ係数誤差ベクトルのノルムが小さくなることを意
味している。つまり、多チャネル受話信号間の相互相関
の変動は、反響消去装置内において反響路推定が誤って
いる場合に、反響の増大をひき起こすが、無限回の変動
が繰り返された末には、反響路推定部における真の反響
路の推定を可能とする有効な情報であると捉えることが
できる。

【００１３】そこで、この発明では前述したように、各
チャネルの受話信号間の相互相関の変動の情報を抽出し
て利用する。従来装置の構成においては、各チャネル間
の受話信号ベクトルを結合して扱う。このため、各チャ
ネルの受話信号間の相互相関は、受話信号結合ベクトル
と過去の受話信号結合ベクトルとの相関関係である自己
相関の中に現れる。よって、各チャネルの受話信号につ
いて、相互相関の変動していない成分を除去して、相互
相関の変動した成分を抽出するためには、受話信号結合
ベクトルの自己相関を小さくする、つまり受話信号結合
ベクトルの自己相関を除去する手法を用いればよい。こ
のような手法として、１チャネル反響消去装置におい
て、射影アルゴリズム、ＲＬＳアルゴリズム、あるいは
線形予測を用いる方法等が提案されている。これらのう
ち、線形予測による手法は、音声信号を線形予測して、
その線形予測係数を用いて音声信号の自己相関を除去す
るものである。ＲＬＳアルゴリズムは、射影アルゴリズ
ムの特殊な場合とみなすことができる。そこで、この発
明では、最も現実的な手法の例として、射影アルゴリズ
ムを用いた場合における受話信号結合ベクトルの自己相
関除去手法を説明する。

【００１４】射影アルゴリズムは、音声のように自己相
関のある信号に対する収束速度を改善する手法として提
案されている。射影アルゴリズムの詳細は、文献「尾
関，梅田，“アフィン部分空間への直交射影を用いた適
応フィルタアルゴリズムとその諸性質”，電子情報通信
学会誌（Ａ），Ｊ６７−Ａ，ｐｐ．１２６−１３２，
（昭５９−２）」に記載されている。

【００１５】一般に、ｐ次の射影アルゴリズムでは、過
去のｐ個の入力信号ベクトルｘ（ｋ），ｘ（ｋ−
１），…，ｘ（ｋ−ｐ＋１）に対して正しい出力ｙ
（ｋ），ｙ（ｋ−１），…，ｙ（ｋ−ｐ＋１）を得るよ
うに擬似反響路ベクトルｈ＾（ｋ）を修正する。すな
わち、ｘ^T（ｋ）ｈ＾（ｋ＋１）＝ｙ（ｋ）ｘ^T（ｋ−１）ｈ＾（ｋ＋１）＝ｙ（ｋ−１）（21）・・・ｘ^T （ｋ−ｐ＋１）ｈ＾（ｋ＋１）＝ｙ（ｋ−ｐ＋１）を満足するｈ＾（ｋ＋１）を求める。但し、ｘ（ｋ）＝〔ｘ（ｋ），ｘ（ｋ−１），…，ｘ（ｋ−Ｌ＋１）〕^T （22）である。未知数の数（タップ数）Ｌより方程式の数ｐが
少ないとき、式（２１）の連立方程式の解ｈ＾（ｋ＋
１）は不定となる。そこで、修正の大きさ‖ｈ＾（ｋ
＋１）−ｈ＾（ｋ）‖が最小となるように修正を行な
う。ｐ次の射影アルゴリズムは、ｈ＾(k+1) ＝ｈ＾(k)+α [Ｘ^T(k)]⁺ｅ(k) ＝ｈ＾(k)+αＸ(k)[Ｘ^T(k) Ｘ(k)]^-1ｅ(k) ＝ｈ＾(k)+αＸ(k) β(k) ＝ｈ＾(k)+α [β₁ ｘ(k)+β₂ ｘ(k-1）＋… ＋β_pｘ(k-p＋1)] （23）但し、Ｘ(k) ＝ [ｘ(k),ｘ(k-1),…, ｘ(k-p+1)] （24）ｅ(k) ＝ [e(k),(1-α)e(k-1),…,(1-α) ^p-1e(k-p+1)]^T （25）ｅ(k) ＝ y(k) − y＾(k) （26）ｙ＾(k) ＝ｈ＾^T(k) ｘ(k) （27） β(k) ＝〔β₁ ，β₂ ，…，β_p〕^T （28）⁺ ：一般化逆行列^-1 ：逆行列と表される。β（ｋ）はｐ元連立一次方程式〔Ｘ^T（ｋ）Ｘ（ｋ）〕β（ｋ）＝ｅ（ｋ）（29）の解である。式（２３）の右辺第２項は修正ベクトルで
あり、修正ベクトルにより擬似反響路ベクトルが逐次更
新される。射影アルゴリズムには演算量の低減をはかっ
た高速算法がいくつか提案されており、それらの詳細は
特願平６−１０２４１８、特願平５−１８４７４２に記
載されている。また、負の時刻の入出力を零として、ｐ
を無限大とした場合が、ＲＬＳアルゴリズムに対応す
る。

【００１６】以上より、ｈ＾（ｋ）の更新式を従来装
置で用いられているＮＬＭＳアルゴリズムによる場合と
射影アルゴリズムによる場合とで比較すると、ＮＬＭＳ
アルゴリズムによる式（１２）では、ｈ＾（ｋ）は入
力信号ベクトルｘ（ｋ）の方向に更新される。一方、
射影アルゴリズムによる式（２３）では、ｈ＾（ｋ）
は式（２３）の第４式右辺第２項以降をｖ(k) ＝β₁ ｘ(k) ＋β₂ ｘ(k-1) ＋…＋β_pｘ(k-p+1) （30）とおいて、このｖ（ｋ）の方向、すなわち、入力信号
ベクトルｘ（ｋ）から過去の入力信号ベクトルｘ
（ｋ−１），…，ｘ（ｋ−ｐ＋１）の全てとの相関を
除去したベクトルの方向に更新される。つまり入力信号
中の過去の信号と同様のものをなるべく除去した信号と
ｖ（ｋ）がなるようにβ₁ 〜β_pを決める、式（３
０）のベクトルｖ（ｋ）は、ｘ（ｋ）が相互相関の
高い多チャネル信号ベクトルの結合ベクトルである場
合、各チャネル間の相互相関の変わらない成分は除去さ
れて、相互相関の変動成分が強調されたベクトルとな
る。

【００１７】次にこの発明の第２の方法を具体的に説明
する。各チャネルの受話信号間の相互相関を能動的に変
動させる機能を付加して、各再生器で音響信号を再生
し、また、このようにして受話信号間の相互相関の変動
が付加された受話信号を疑似反響路ベクトルの修正ベク
トルの導出に利用する。各チャネルの受話信号間の相互
相関に変動付加する方法としては、フィルタ処理によ
る実現方法と、信号乗算処理による実現方法と、信
号加算処理による実現方法と、ピッチシフト処理によ
る実現方法とがある。以下、順に〜各実現方法を説
明する。

【００１８】フィルタ処理による実現方法各チャネルの受話信号ｘ₁(ｋ），ｘ₂(ｋ），…，ｘ
_N（ｋ）を、異なる時変特性をもつ時変フィルタに印
加、すなわち、これら時変フィルタのインパルス応答ｆ
₁(ｋ），ｆ₂(ｋ），…，ｆ_N（ｋ），畳み込み演算を＊
で表すとして、ｘ₁(ｋ）＝ｆ₁(ｋ）＊ｘ₁(ｋ）ｘ₂(ｋ）＝ｆ₂(ｋ）＊ｘ₂(ｋ）ｘ￣_N（ｋ）＝ｆ_N（ｋ）＊ｘ_N（ｋ） (31) を満足するｘ￣₁(ｋ），ｘ￣₂(ｋ），…，ｘ￣_N（ｋ）
に変換することにより、各チャネルの受話信号間の相互
相関への変動の付加を実現する。

【００１９】信号乗算処理による実現方法各チャネルの受話信号ｘ₁(ｋ），ｘ₂(ｋ），…，ｘ
_N（ｋ）を、異なる関数ｇ ₁(ｋ），ｇ₂(ｋ），…，ｇ_N
（ｋ）と乗算、すなわち、ｘ￣₁(ｋ）＝ｇ₁(ｋ）・ｘ₁(ｋ）ｘ￣₂(ｋ）＝ｇ₂(ｋ）・ｘ₂(ｋ）・・・ (32) ｘ￣_N(k）＝ｇ_N（ｋ）・ｘ_N（ｋ）を満足するｘ￣₁(ｋ），ｘ￣₂(ｋ），…，ｘ￣_N（ｋ）
に変換することにより、各チャネルの受話信号間の相互
相関への変動の付加を実現する。

【００２０】信号加算処理による実現方法各チャネルの受話信号ｘ₁(ｋ），ｘ₂(ｋ），…，ｘ
_N（ｋ）を、異なる関数ｎ ₁(ｋ），ｎ₂(ｋ），…，ｎ_N
（ｋ）と加算、すなわち、ｘ￣₁(ｋ）＝ｘ₁(ｋ）＋ｎ₁(ｋ）ｘ￣₂(ｋ）＝ｘ₂(ｋ）＋ｎ₂(ｋ）・・・ (33) ｘ￣_N(k）＝ｘ_N（ｋ）＋ｎ_N（ｋ）を満足するｘ￣₁(ｋ），ｘ￣₂(ｋ），…，ｘ￣_N（ｋ）
に変換することにより、各チャネルの受話信号間の相互
相関への変動の付加を実現する。

【００２１】ピッチシフト処理による実現方法各チャネルの受話信号ｘ₁(ｋ），ｘ₂(ｋ），…，ｘ
_N（ｋ）を、各受話信号の周波数特性について、各チャ
ネルごとに異なる時変の周波数伸縮処理、すなわち、ピ
ッチシフト処理することにより得られるｘ￣₁(ｋ），ｘ
￣₂(ｋ），…，ｘ￣ _N（ｋ）に変換することにより、各
チャネルの受話信号間の相互相関への変動の付加を実現
する。

【００２２】

【実施例】図１にこの発明による第１の方法の実施例を
示す。図５と対応する部分に同一符号を付けてある。こ
の発明では各反響消去装置２２₁〜２２_mに相互相関変
動成分抽出部３１を設ける。相互相関変動成分抽出部３
１に各チャネルの受話信号と残留反響信号とを入力し、
修正ベクトルを生成する。

【００２３】この実施例では受話信号結合ベクトルの自
己相関の除去を射影アルゴリズムを用いて行なうことに
より受話信号の相互相関の変動を検出して、修正ベクト
ルを得る。つまり、相互相関変動成分抽出部３１では、
受話信号結合ベクトルｘ（ｋ）と残留反響信号ｅ
_m（ｋ）の情報を蓄積し、式（２９）に対応する［Ｘ^T（Ｋ）Ｘ（ｋ）］β（ｋ）＝ｅ_m（ｋ）（34）を解いてβ（ｋ）を求め、式（３０）より自己相関除去
ベクトル、つまり修正ベクトルｖ（ｋ）を生成する。
修正部３２では、現在の推定反響路ベクトルｈ＾
_m（ｋ）に対し、修正ベクトルｖ（ｋ）を用いて式
（２３），（３０）より、ｈ＾_m（ｋ＋１）＝ｈ＾_m（ｋ）＋αｖ（ｋ）（35）を演算して擬似反響路結合ベクトルｈ＾_m（ｋ）を更
新し、これを擬似反響信号生成部１８ｍへ出力する。

【００２４】相互相関変動成分抽出部３１により多チャ
ネル受話信号ベクトル間の相互相関の変動成分が抽出さ
れ、前記各受話信号ベクトル間の相互相関の変動成分が
強調された新たな多チャネル信号ベクトルが得られ、こ
の多チャネル信号ベクトルを修正部３２に与えるように
したため、真の反響路の推定速度を速めることができ
る。

【００２５】図９にこの発明の第２の方法の実施例を示
す。図４と対応する部分に同一符号を付けてある。この
発明では、相互相関変動付加部２４が設けられる。相互
相関変動付加部２４で、各チャネルの受話信号ｘ
₁(ｋ），ｘ₂(ｋ），…，ｘ_N（ｋ）は、それらの相互相
関の変動を能動的に付加されて、ｘ￣₁(ｋ），ｘ￣
₂(ｋ），…，ｘ￣_N（ｋ）に変換され、これらｘ￣
₁(ｋ），ｘ￣₂(ｋ），…，ｘ￣_N（ｋ）を、スピーカ１
２₁，１２₂，…，１２_Nから再生するとともに、反響
消去システム２３への入力とする。

【００２６】相互相関変動付加部２４の実現について
は、変換された受話信号ｘ￣₁(ｋ），ｘ￣₂(ｋ），…，
ｘ￣_N（ｋ）が音響信号として再生されるときに、聴覚
上の品質が損なわれないように注意する必要がある。そ
こで、以下に相互相関変動付加部２４の実現において前
節で述べた方法の具体例を挙げる。フィルタ処理による実現の具体例各チャネルの受話信号ｘ₁(ｋ），ｘ₂(ｋ），…，ｘ
_N（ｋ）を、異なる時変特性をもつ時変フィルタのイン
パルス応答ｆ₁(ｋ），ｆ₂(ｋ），…，ｆ_N（ｋ）を用い
て、式（３１）を満足するｘ￣₁(ｋ），ｘ￣₂(ｋ），
…，ｘ￣_N（ｋ）に変換する場合の、時変フィルタの構
成例を示す。人の聴覚による音声知覚においては、主と
して振幅情報が重要であり、位相情報は余り重要はでな
いと言われており、このことは、例えば文献「古井，
“ディジタル音声処理”，東海大学出版会」で述べられ
ている。そこで、入力信号の周波数に対して振幅特性が
平坦な特性を保つように、時変フィルタを実現すれば、
聴覚的にフィルタ処理による影響が小さくなると考えら
れる。

【００２７】この目的に適した構成として、全域通過伝
達関数をもつＩＩＲフィルタがある。このフィルタの伝
達関数Ｆ（ｚ）は、次数をＫとして、Ｆ(z) ＝Ａ｛ｚ^-k−ａ₁ｚ^-(K-1)＋…(-1)^K-1ａ_k-1ｚ^-1＋(-1)^Kａ_K｝／｛１−ａ₁ｚ^-1＋…(-1)^K-1ａ_K-1ｚ^-(K-1)＋(-1)^Kａ_Kｚ^-K｝ (36) と一般に表すことができる。標本化間隔をＴ，角周波数
をωとすると、振幅特性

【００２８】

【数１】式は、フィルタ係数ａ₁，ａ₂，…，ａ_Mがどのように
定められても分母と分子で相殺し合い、ωによらない一
定利得Ａをとる。一方、位相特性はフィルタ係数ａ₁，
ａ₂，…，ａ_Mにより異なる。このａ₁，ａ₂，…，ａ
_Mを時間とともに変化させた場合には、時変の位相特性
を有することになるが、振幅特性への影響はない。

【００２９】さて、ここでは構成が簡単で実現の容易な
２次の全域通過フィルタを考える。２次の全域通過フィ
ルタの構成は、例えば文献「辻井，“ディジタル信号処
理の基礎”，電子情報通信学会」に詳しく述べられてい
る。式（３６）について、簡単のため利得Ａ＝１とする
と、２次の全域通過伝達関数Ｆ（ｚ）は、Ｆ(z) ＝（ｚ^-2−ａ₁ｚ^-1＋ａ₂）／（１−ａ₁ｚ^-1＋ａ₂ｚ^-2） (37) となる。この伝達関数をもとに、ラチス型で構成する
と、図１０に示す通りとなり、式（３７）の係数はａ₁＝−ｒ（１＋ｒ₂） (38) ａ₂＝ｒ₂ (39) と置き換えられＦ(z) ＝｛ｚ^-2＋ｒ₁(1+ｒ₂)ｚ^-1＋ｒ₂｝／｛１＋ｒ₁(1+ｒ₂)ｚ^-1＋ｒ₂ｚ^-2｝ (40) と表される。本発明において、このようなラチス型で構
成する利点は、ｒ₁,ｒ₂が独立に意味をもつパラメータ
として扱えることにある。式（４０）のＦ（ｚ）の群遅
延特性がピークを持つ場合、ピーク周波数はほぼｒ₁の
みに依存して決まり、ピークの急峻さはほぼｒ₂のみに
依存して決まる。つまり、フィルタ係数を時間とともに
変化させる場合、すべての係数を変化させなくとも、例
えば、群遅延特性のピーク周波数に対応したｒ₁を再生
チャネルごとに異なる変化をさせるだけで、効果的な時
変特性を得ることができ、構成が簡単になる。これらの
時変特性をチャネルごとに異なるものとすれば、もとの
多チャネル受話信号に対して、それらの間の相互相関を
変動させる効果がある。

【００３０】信号乗算処理による実現の具体例各チャネルの受話信号ｘ₁(ｋ），ｘ₂(ｋ），…，ｘ
_N（ｋ）を、異なる関数ｇ ₁(ｋ），ｇ₂(ｋ），…，ｇ_N
（ｋ）と乗算して、式（３２）を満足するｘ￣₁(ｋ），
ｘ￣₂(ｋ），…，ｘ￣_N（ｋ）に変換する場合の、関数
ｇ₁(ｋ），ｇ₂(ｋ),…，ｇ_N（ｋ）の例を示す。

【００３１】式（３２）において、もし関数ｇ₁(ｋ），
ｇ₂(ｋ），…，ｇ_N（ｋ）が正負の値をとる場合には、
その変動の大きさによらず受話信号の品質は大幅に劣化
すると考えられる。そこで、関数ｇ₁(ｋ），ｇ₂(ｋ），
…，ｇ_N（ｋ）は、常に正、あるいは常に負の値のみを
とり、また、変動の幅も制限されていることが望まし
い。

【００３２】そこで、関数ｇ₁(ｋ），ｇ₂(ｋ），…，ｇ
_N（ｋ）が以下のような形である場合を考える。最大振
幅が１に正規化された異なる関数ｗ₁(ｋ），ｗ₂(ｋ），
…，ｗ_N（ｋ）と、スカラー値ｋ₁，ｋ₂，…，ｋ
_N（０＜‖ｋ₁‖，‖ｋ₂‖，…，‖ｋ_N‖＜１），お
よび利得係数Ａ₁，Ａ₂，…，Ａ_Nを用いて、ｇ
₁(ｋ），ｇ₂(ｋ），…，ｇ_N（ｋ）を、ｇ₁(ｋ）＝Ａ₁(１＋ｋ₁ω₁(ｋ））ｇ₂(ｋ）＝Ａ₂(１＋ｋ₂ω₂(ｋ））・・・ (41) ｇ_N(k）＝Ａ_N（１＋ｋ_Nω_N（ｋ））と表す。また任意の時刻ｋにおいて、−１＜ｋ₁ω
₁(ｋ），ｋ₂ω₂(ｋ），…，ｋ_Nｗ_N（ｋ）＜１が成り
立つから、ｇ₁(ｋ），ｇ₂(ｋ），…，ｇ_N（ｋ）の値の
符号はそれぞれ一定である。聴覚的なバランスを考えれ
ば、利得係数Ａ₁，Ａ ₂，…，Ａ_Nの値はすべて等しい
か、あるいは少なくとも同符号の値であることが望まし
く、この場合、ｇ₁(ｋ），ｇ₂(ｋ），…，ｇ_N（ｋ）は
すべて同符号の値を取る。また、スカラー値ｋ₁，
ｋ₂，…，ｋ_Nの値により、処理後の信号の歪み量を調
節することができる。

【００３３】また、周波数シフトという観点から、式
（４１）とは別に、ｇ₁(ｋ）＝Ａ₁cos(ω₁(ｋ）・ｋ＋φ₁）ｇ₂(ｋ）＝Ａ₂cos(ω₂(ｋ）・ｋ＋φ₂）・・・ (42) ｇ_N(k）＝Ａ_Ncos(ω_N(k) ・ｋ＋φ_N）として、各チャネルごとに異なる時変特性をもち、ある
変動幅の中で変化する関数ω₁(ｋ），ω₂(ｋ），…，ω
_N(k) によって、受話信号ｘ₁(ｋ），ｘ₂(ｋ),…，ｘ_N
（ｋ）の周波数特性をシフトさせることにより、もとの
多チャネル受話信号に対して、相互相関の変動付加を実
現することもできる。

【００３４】信号加算処理による実現の具体例各チャネルの受話信号ｘ₁(ｋ），ｘ₂(ｋ),…，ｘ
_N（ｋ）を、異なる関数ｎ₁（ｋ），ｎ₂(ｋ),…，ｎ_N
（ｋ）と加算して、式（３３）を満足するｘ￣₁(ｋ），
ｘ￣₂(ｋ），…，ｘ￣_N（ｋ）に変換する場合の、関数
ｎ₁(ｋ），ｎ₂(ｋ），…，ｎ_N（ｋ）の例を示す。

【００３５】関数ｎ₁(ｋ），ｎ₂(ｋ），…，ｎ_N（ｋ）
は、聴覚的な影響を考慮すると、各チャネルの受話信号
ｘ₁(ｋ），ｘ₂(ｋ),…，ｘ_N（ｋ）と比べて、振幅値が
適度に小さいことが望ましい。そこで、例えば、各チャ
ネルの受話信号ｘ₁(ｋ），ｘ ₂(ｋ),…，ｘ_N（ｋ）の振
幅の定格値をＡ_rとして、最大振幅がこのＡ_rとなるよ
うに規格化された異なる関数ｎ_r1（ｋ），ｎ_r2（ｋ），
…，ｎ_2N（ｋ）に対して、スカラー値λ₁，λ₂，…，
λ_Nを掛けることにより、ｎ₁(ｋ），ｎ₂(ｋ），…，ｎ
_N（ｋ）を、ｎ₁(ｋ）＝λ₁ｎ_r1（ｋ）ｎ₂(ｋ）＝λ₂ｎ_r2（ｋ）・・・ (43) ｎ_N(k) ＝λ_Nｎ_rN（ｋ）と表す。スカラー値λ₁,λ₂，…，λ_Nの値により、受
話信号ｘ₁(ｋ),ｘ₂(ｋ),…，ｘ_N（ｋ）に加える信号の
大きさを定格値Ａ_rからどの程度小さくするのかを決定
できる。

【００３６】ピッチシフト処理による実現の具体例各チャネルの受話信号ｘ₁(ｋ），ｘ₂(ｋ),…，ｘ
_N（ｋ）を、各受話信号の周波数特性について、各チャ
ネルごとに異なる時変の周波数伸縮処理、すなわち、ピ
ッチシフト処理することにより得られるｘ￣₁(ｋ），ｘ
￣₂(ｋ），…，ｘ￣_N（ｋ）に変換するための実施例を
示す。

【００３７】ピッチシフト処理は、周波数伸縮ととも
に、時間伸縮をともない、結局、時間領域で処理するこ
とが可能であるが、時間軸を縮小させれば、もとの信号
が占めていた占有時間と等しくするための補間処理が必
要となり、逆に時間軸を伸長させれば、もとの信号より
も占有時間が長くなるため間引き処理が必要となる。こ
れらの補間・間引き処理は、無音区間を検出するなどし
て、その無音区間を引き延ばしたり、削除することによ
り実現できる。

【００３８】

【発明の効果】この発明方法は、２地点間において、多
チャネルの伝送系と、多チャネルの集音・再生系を持つ
端末装置とを用いることにより、双方の音響空間情報の
伝達を可能とする臨場感の高い音声通信を実現する際に
反響を消去するために適用することができる。

【００３９】簡単のため、図６に示したような、２地点
間を２チャネルの伝送系で結び、２チャネルの集音・再
生系を持つ端末装置を用いて、ステレオ音声通信を行な
う場合を考える。例えば、Ａ地点から複数の異なる話者
が一人ずつそれぞれの席から発話した場合に，集音され
るステレオ音声信号には、それぞれの発話者と２本のマ
イクロホンの位置関係に依存した相互相関がある。それ
らの音声信号がＢ地点側において再生されるとき、従来
の反響消去装置を用いた場合では、反響路推定がステレ
オ信号の相互相関に依存して誤るため、話者の交替の度
にＡ地点側に大きな反響が戻ってしまう。ところが、こ
の発明の反響消去方法を用いた場合、話者の交替による
ステレオ再生信号の相互相関の変動を抽出して、反響路
推定の誤りを修正するため、話者の交替の度に起きてい
た反響の増大を低減できる。また、人間の発話音声を複
数のマイクロホンによって集音するとき、固定した位置
から発話した場合でも、集音された多チャネル音声信号
の相互相関は微妙に変動している。この発明方法では、
このような信号間の相互相関の微妙な変動をも有効に利
用できる。

【００４０】図２Ａは、体や頭を固定して発話した音声
をステレオ集音した信号を受話信号として用いた計算機
シミュレーション結果であり、従来方法とこの発明方法
とにおいて、反響路推定部で生成された擬似反響路結合
ベクトルと真の反響路結合ベクトルとの誤差ベクトルの
大きさを比較したものである。ここで、この発明方法の
相互相関変動成分抽出には、２次の射影アルゴリズムを
用いている。この結果より、この発明方法が信号間の相
互相関の微妙な変動をも捉え、反響路推定の誤りを修正
できることが分かる。この発明ではチャネル間の微妙な
変動を過去のデータを利用することにより、有効にとら
えることにより、行なわれていると言える。

【００４１】図２Ｂに図２Ａと同様なステレオ集音の場
合の擬似反響路ベクトルの収束の様子を射影アルゴリズ
ムの次数をパラメータとして示す。タップ数を５００＋
５００、ステップサイズα＝０．５とした場合である。
この図からもこの発明が有効であることが理解される。
ステレオ再生装置を有する端末において、図６に示した
ように相手側で２チャネル集音されたステレオ信号を受
信してそのまま再生する方法の他に、多地点間通信等で
は、対地毎の受話信号を受信側で任意に音像定位処理す
ることにより、受話環境を快適にする方法が考えられて
いる。この発明方法は、このような多地点間通信用端末
にも適用可能である。図８は、４地点間通信の構成を示
したものである。各地点において、集音は１チャネル
（モノラル）である。いま、Ｄ地点を例にとって説明す
る。Ｄ地点において、Ａ、Ｂ、Ｃ地点からの受話信号に
ついて、それぞれ右、中央、左に定位するように音像定
位処理を受話信号に処理部１３ｄで施すことにより、新
たな２チャネルステレオ受話信号を生成して２チャネル
のステレオ再生を行なうものとする。この場合において
も、この発明方法を用いることにより、先の場合と同様
に、異なる対地の話者へ話者交替するとき等に、ステレ
オ受話信号の相互相関の変動を抽出して、反響路推定の
誤りを修正するため、従来方法を用いた場合よりも反響
の増大を低減できる。

【００４２】ところで図１１は、一人の話者の発話する
音声を１本のマイクロホンにより集音した後、計算機上
で２系統の異なる時不変フィルタにそれぞれ通過させる
ことにより得られた２チャネル信号を受話信号として用
いた。つまり、受話信号間の相互相関が一定の場合であ
り、この場合は、受話信号間の相互相関の変動抽出・利
用をしようにも、受話信号間の相互相関の変動そのもの
がないため、従来法と比較しても、反響路推定は停滞し
て改善されない。

【００４３】そこで、上記方法に本発明の第２の方法を
適用すれば、受話信号間の相互相関は常に変動すること
になり、反響路推定が停滞することがなくなる。このと
きの構成は、図９と同じであるが、Ｎチャネル反響消去
装置２２₁，２２₂，…，２２_Mの構成は図１のように
なる。図１２に示す計算機シミュレーション結果は、図
１１で用いたのと同じ受話信号、すなわち、チャネル間
の相互相関が一定なステレオ音声信号に対して、相互相
関の変動を付加した場合と、付加しなかった場合とにお
ける受話信号の相互相関の変動成分を抽出・利用する方
法の真の反響路と疑似反響路との誤差を時間経過ととも
に描いたものである。相互相関の変動を付加しなかった
場合のデータは、図１１からの再掲であり、真の反響路
への収束は停滞している。一方、相互相関の変動を付加
した場合には、これが停滞することなく、図２Ａと同様
な良好な結果が得られている。また、このときのステレ
オ（２チャネル）信号への相互相関の変動の付加は、式
（３２）に示す信号乗算処理により実現し、ｇ₁(ｋ），
ｇ₂(ｋ）は、ｇ₁(ｋ）＝１＋0.０５・ω（ｋ） (44) ｇ₂(ｋ）＝１−0.０５・ω（ｋ） (45) とし、但し、ω（ｋ）は、最大振幅が１の白色雑音信号
である。

【００４４】すなわち、一般に多チャネル再生系と１チ
ャネル以上の集音系で構成される通信会議システムに、
従来の反響消去方法を適用した場合、各チャネルの受話
信号間に相互相関がある場合に、疑似反響路の推定が正
しく行われないため、受話信号間に相互相関が変動する
度に反響が増大する問題が生じていた。この発明の第１
の方法においては、受話信号間の相互相関の変動成分を
抽出・利用することにより、反響路推定の誤りを修正す
るため、上記のような問題を改善する効果がある。

【００４５】また、本発明の第２の方法においては、受
話信号間の相互相関に変動を付加する機能を有するた
め、反響路推定において誤った解に完全に停滞すること
がなく、真の反響路との誤差を小さくする方向へ反響路
推定を続けるので、上記のような問題を改善する効果が
ある。

【図面の簡単な説明】

【図１】この発明による第１の方法を適用した多チャネ
ル反響消去装置の構成を示すブロック図。

【図２】この発明方法による擬似反響路ベクトルの収束
の様子を示す図。

【図３】１チャネル反響消去装置を示すブロック図。

【図４】多チャネル反響消去装置システムを示すブロッ
ク図。

【図５】従来の多チャネル反響消去装置を示すブロック
図。

【図６】ステレオ音声通信系を示すブロック図。

【図７】Ａは受話信号間に相互相関がある場合の反響路
推定動作を示す図。Ｂは前記相互相関が変動した場合の
反響路推定動作を示す図である。

【図８】４地点間通信において音像定位機能をもたせる
音声通信会議方式を示すブロック図。

【図９】この発明による第２の方法を適用した多チャネ
ル反響消去装置の構成を示すブロック図。

【図１０】この発明の第２の方法を実現するディジタル
フィルタの一例を示す図。

【図１１】相互相関の変動を抽出する方法と、従来法と
の真の反響路と疑似反響路との誤差ベクトルの大きさを
比較した図。

【図１２】相互相関が一定なステレオ信号を受話したと
きに、相関変動抽出法に、この発明の第２の方法を適用
した場合と、しなかった場合を比較した図。

Claims

【特許請求の範囲】

【請求項１】各チャネルの受話信号をチャネルごとの
再生器で音響信号に再生し、これら各音響信号が、上記各再生器から集音器に到る各
反響路を経由して集音され、その集音器からの反響信号から擬似反響信号を差し引
き、その残りである残留反響信号と、上記各チャネルの受話
信号とから修正ベクトルを求め、この修正ベクトルを用いて上記各反響路のインパルス応
答の推定を逐次修正し、その修正されたインパルス応答をもつ各擬似反響路を生
成し、これら各擬似反響路に上記受話信号の対応するものをそ
れぞれ印加して上記擬似反響信号を生成する多チャネル
音声通信会議用反響消去方法において、上記各チャネルの受話信号の相互相関の変動成分を抽出
し、この変動成分を上記修正ベクトルとすることを特徴とす
る多チャネル音声通信会議用反響消去方法。
【請求項２】上記各チャネルの受話信号の系列よりな
る受話信号ベクトルを結合させた受話信号結合ベクトル
を求め、この受話信号結合ベクトルの過去のものに対す
る現在のそれとの相関の変動を検出して上記相互相関の
変動成分とすることを特徴とする請求項１記載の多チャ
ネル音声通信会議用反響消去方法。
【請求項３】各チャネルの受話信号をチャネルごとの
再生器で音響信号に再生し、これら各音響信号が、上記
各再生器から集音器に到る各反響路を経由して集音さ
れ、その集音器から疑似反響信号を差し引き、その残り
である残留反響信号と、上記各チャネルの受話信号とか
ら修正ベクトルを求め、この修正ベクトルを用いて上記
各反響路のインパルス応答の推定を逐次修正し、その修
正されたインパルス応答をもつ各疑似反響路を生成し、
これら各疑似反響路に上記受話信号の対応するものをそ
れぞれ印加して、上記疑似反響信号を生成する多チャネ
ル音声通信会議用反響消去方法において、上記各チャネルの受話信号間の相互相関に変動成分を付
加して、上記再生器で音響信号を再生し、相互相関の変
動成分が付加された受話信号を上記修正ベクトルの導出
に利用することを特徴とする多チャネル音声通信会議用
反響消去方法。
【請求項４】再生チャネル数をＮとして、上記各チャ
ネルの受話信号それぞれを、離散時間ｋの関数として、
ｘ₁(ｋ），ｘ₂(ｋ），…，ｘ_N（ｋ）と表したとき、こ
れらの受話信号を各チャネルに対して異なる時変特性を
もつ時変フィルタに印加、これら時変フィルタのインパ
ルス応答ｆ₁(ｋ），ｆ₂(ｋ），…，ｆ _N（ｋ），畳み込
み演算を＊で表すとして、ｘ￣₁(ｋ）＝ｆ₁(ｋ）＊ｘ₁(ｋ）ｘ￣₂(ｋ）＝ｆ₂(ｋ）＊ｘ₂(ｋ）・・・ｘ￣_N（ｋ）＝ｆ_N（ｋ）＊ｘ_N（ｋ）を満足するｘ￣₁(ｋ），ｘ￣₂(ｋ），…，ｘ￣_N（ｋ）
に変換することにより、各チャネルの受話信号間の相互
相関への変動の付加を実現することを特徴とする請求項
３記載の多チャネル音声通信会議用反響消去方法。
【請求項５】再生チャネル数をＮとして、上記各チャ
ネルの受話信号それぞれを、離散時間ｋの関数として、
ｘ₁(ｋ），ｘ₂(ｋ），…，ｘ_N（ｋ）と表したとき、こ
れらの受話信号を各チャネルに対して異なる関数ｇ
₁(ｋ），ｇ₂(ｋ），…，ｇ_N（ｋ）と乗算し、ｘ￣₁(ｋ）＝ｇ₁(ｋ）・ｘ₁(ｋ）ｘ￣₂(ｋ）＝ｇ₂(ｋ）・ｘ₂(ｋ）・・・ｘ￣_N（ｋ）＝ｇ_N（ｋ）・ｘ_N（ｋ）を満足するｘ￣₁(ｋ），ｘ￣₂(ｋ），…，ｘ￣_N（ｋ）
に変換することにより、上述の各チャネルの受話信号間
の相互相関への変動の付加を実現することを特徴とする
請求項３記載の多チャネル音声通信会議用反響消去方
法。
【請求項６】再生チャネル数をＮとして、上記各チャ
ネルの受話信号それぞれを、離散時間ｋの関数として、
ｘ₁(ｋ），ｘ₂(ｋ），…，ｘ_N（ｋ）と表したとき、こ
れらの受話信号を各チャネルに対して異なる関数ｎ
₁(ｋ），ｎ₂(ｋ),…，ｎ_N（ｋ）と加算し、ｘ￣₁(ｋ）＝ｘ₁(ｋ）＋ｎ₁(ｋ）ｘ￣₂(ｋ）＝ｘ₂(ｋ）＋ｎ₂(ｋ）・・・ｘ￣_N（ｋ）＝ｘ_N（ｋ）＋ｎ_N（ｋ）を満足するｘ￣₁(ｋ），ｘ￣₂(ｋ），…，ｘ￣_N（ｋ）
に変換することにより、上述の各チャネルの受話信号間
の相互相関への変動の付加を実現することを特徴とする
請求項３記載の多チャネル音声通信会議用反響消去方
法。
【請求項７】再生チャネル数をＮとして、上記各チャ
ネルの受話信号それぞれを、離散時間ｋの関数として、
ｘ₁(ｋ），ｘ₂(ｋ），…，ｘ_N（ｋ）と表したとき、こ
れらの受話信号を各受話信号の周波数特性について、各
チャネルごとに異なる時変の周波数伸縮処理することに
より得られるｘ￣₁(ｋ），ｘ￣₂(ｋ),…，ｘ￣_N（ｋ）
に変換することにより、各チャネルの受話信号間の相互
相関への変動の付加を実現することを特徴とする請求項
３記載の多チャネル音声通信会議用反響消去方法。