JP3654470B2 - サブバンド多チャネル音声通信会議用反響消去方法 - Google Patents
サブバンド多チャネル音声通信会議用反響消去方法 Download PDFInfo
- Publication number
- JP3654470B2 JP3654470B2 JP24352496A JP24352496A JP3654470B2 JP 3654470 B2 JP3654470 B2 JP 3654470B2 JP 24352496 A JP24352496 A JP 24352496A JP 24352496 A JP24352496 A JP 24352496A JP 3654470 B2 JP3654470 B2 JP 3654470B2
- Authority
- JP
- Japan
- Prior art keywords
- echo
- signal
- channel
- subband
- signals
- 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.)
- Expired - Lifetime
Links
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04M—TELEPHONIC COMMUNICATION
- H04M9/00—Arrangements for interconnection not involving centralised switching
- H04M9/08—Two-way loud-speaking telephone systems with means for conditioning the signal, e.g. for suppressing echoes for one or both directions of traffic
- H04M9/082—Two-way loud-speaking telephone systems with means for conditioning the signal, e.g. for suppressing echoes for one or both directions of traffic using echo cancellers
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04M—TELEPHONIC COMMUNICATION
- H04M1/00—Substation equipment, e.g. for use by subscribers
- H04M1/60—Substation equipment, e.g. for use by subscribers including speech amplifiers
Landscapes
- Engineering & Computer Science (AREA)
- Signal Processing (AREA)
- Cable Transmission Systems, Equalization Of Radio And Reduction Of Echo (AREA)
- Telephonic Communication Services (AREA)
Description
【産業上の利用分野】
この発明は、多チャネル再生系を有する通信会議システムにおいて、ハウリングの原因および聴覚上の障害となる室内反響信号を消去する多チャネル音声通信会議用反響消去方法に関し、特に分割された周波数帯域毎に擬似反響路の推定インパルス応答を射影アルゴリズム等により修正するサブバンド多チャネル音声通信会議用反響消去方法に関するものである。
【0002】
【従来の技術】
1チャネル用反響消去
同時通話性能に優れ反響感の少ない拡声通話システムを提供するために、反響消去装置がある。
まず、1チャネル用の反響消去装置について、その反響消去方法および装置構成を図1を参照して説明する。拡声通話において、受話端子11に相手の発話等で得られる受話信号は、そのままスピーカ(再生器)12から再生される場合と、スピーカ12へ送る前に、受話信号の振幅やパワー等の大きさに応じて自動的に利得を調節する受話信号加工部13が挿入され、受話信号に何らかの加工が施された後に、スピーカ12から再生される場合とがある。このため、この明細書で受話信号x1(k)とは、相手からの受話信号そのものとは限らず、受話信号加工部13が挿入されたときは、加工された後の受話信号を指すものとする。kは離散時間を表す。反響消去装置14は、受話信号x1(k)がスピーカ12から反響路15を経て、マイクロホン(集音器)16に集音されて得られる反響信号y1(k)を消去する。ここで、反響信号y1(k)は、時刻kにおける反響路15のインパルス応答をh11(k, l)として、
のような畳み込み演算で得られるものとモデル化できる。Lはタップ数で、反響路15の残響時間に対応させて、あらかじめ設定しておく定数である。まず、受話信号蓄積・ベクトル生成部17において、受話信号x1(k)をL-1時刻過去のものまで蓄積しておく。蓄積された信号は、受話信号ベクトルx 1(k)、すなわち
x 1(k)=[x1(k), x1(k-1),…, x1(k−L+1)]T (2)
として出力される。但し、*Tはベクトルの転置を表す。擬似反響路18では、式(2)の受話信号ベクトルx 1(k)と、反響路推定部19から得られる擬似反響
路ベクトルh^11(k)との内積演算
y^1(k)=h^11 T(k)x 1(k) (3)
を行ない、その結果として、擬似反響信号y^1(k)を生成する。この内積演算は、式(1)のような畳み込み演算と等価である。反響路推定部19では、擬似反響路18で用いる擬似反響路ベクトルh^11(k)を生成する。
【0003】
ここでスピーカ12からマイクロホン16に至る反響路15のインパルス応答h11(k, l)は例えば人や物の移動による音場の変化により変動するので、この経時変動に追従して擬似反響路ベクトルh^11(k)を変化させる必要がある。こ
の構成例において反響消去装置14は適応形FIRフィルタを用いて構成する。この反響路推定に用いる最も一般的なアルゴリズムは、NLMSアルゴリズム(学習同定法)である。NLMSアルゴリズムでは、時刻kにおける受話信号ベクトルx 1(k)と、残留反響信号e1(k)、すなわち減算器21でマイクロホン16
の出力y1(k)から擬似反響信号y^1(k)を差し引いた
e1(k)=y1(k)−y^1(k) (4)
とから、時刻k+1において用いる擬似反響路ベクトルh^11(k+1)を次式の
ように求める。
但し、μはステップサイズパラメータと呼ばれ0<μ<2の範囲で適応動作の調整に用いる。以上のような処理を繰り返すことにより、反響路推定部19では、次第に擬似反響路ベクトルh^11(k)を、真の反響路15のインパルス応答h1
1(k, l)の時系列を各要素として持つ反響路ベクトルh 11(k)、すなわち
h 11(k)=[h11(k,0),h11(k,1),…, h11(k,L-1)]T (6)
と一致させることが可能となり、その結果式(4)の残留反響信号e1(k)を小さくすることができる。
【0004】
現在、この反響路推定に用いる最も有効なアルゴリズムは、射影アルゴリズムあるいはES射影アルゴリズムである。射影アルゴリズムは、アルゴリズム内部で入力信号の自己相関を取り除くことにより、音声のように相関のある信号に対する収束速度を改善するという考え方に基づいている。自己相関成分を除去することは時間領域において信号の白色化を行っていることを意味している。射影アルゴリズムの詳細は、文献「尾関、梅田、“アフィン部分空間への直交射影を用いた適応フィルタアルゴリズムとその諸性質”、電子情報通信学会誌(A), J67−A, pp. 126−132, (昭59−2)」に記載されている。
【0005】
一般に、p次の射影アルゴリズムでは、過去のp個の入力信号ベクトルx(k
), x(k-1), …, x(k-p+1)に対して正しい出力y(k), y(k-1), …, y(k-
p+1)を得るように擬似反響路ベクトルh^(k)を修正する。すなわち、
x T(k)h^(k+1)=y(k)
x T(k-1)h^(k+1)=y(k-1) (7)
・
・
・
x T(k-p+1)h^(k+1)=y(k-p+1)
を満足するh^(k+1)を求める。但し、
x(k)=〔x(k), x(k-1), …, x(k−L+1)〕T (8)
である。未知数の数(タップ数)Lより方程式の数pが少ないとき、式(7)の連立方程式の解h^(k+1)は不定となる。そこで、修正の大きさ‖h^(k+1)−
h^(k)‖が最小となるように修正を行なう。p次の射影アルゴリズムは、
但し、
+ :一般化逆行列
-1 :逆行列
と表される。β(k)はp元連立一次方程式
〔X T(k)X(k)〕β(k)=e(k) (15)
の解である。逆行列演算における不安定性を回避するために小さな正の定数δを用いて
〔X T(k)X(k)+δI〕β(k)=e(k) (15)'
としてもよい。ただしIは単位行列である。式(9)の右辺第2項は修正ベクト
ルであり、修正ベクトルにより擬似反響路ベクトルが逐次更新される。式(9)中のX(k)β(k)は入力信号の自己相関を取り除く処理を表しており、従ってこ
れは時間領域内での信号の白色化処理を意味している。即ち、射影アルゴリズムは時間領域における入力信号の白色化によりインパルス応答の修正速度を高めているといえる。射影アルゴリズムには演算量の低減をはかった高速算法がいくつか提案されており、それらの詳細は特開平7−312535号公報、特開平7−92980号公報に記載されている。また、負の時刻の入出力を零として、pを無限大とした場合が、RLSアルゴリズムに対応する。
【0006】
ES射影アルゴリズムは、射影アルゴリズムに反響路の変動特性に着目したESアルゴリズムを組み合わせたもので、射影アルゴリズムよりさらに収束速度の大きな反響消去装置を提供することができる。p次のES射影アルゴリズムは、
ただし、
A= diag[α1, α2, …,αL]:ステップサイズ行列
αi= α0λi-1(i=1,2,…,L)
λ:インパルス応答変動量の減衰率(0<λ<1)
μ:第2のステップサイズ(スカラ量)
と表される。β(k)はp元連立一次方程式
〔X T(k)AX(k)〕β(k)=e(k) (17)
の解である。逆行列演算における不安定性を回避するために小さな正の定数δを用いて
〔X T(k)AX(k)+δI〕β(k)=e(k) (17)'
としてもよい。ただしIは単位行列である。
【0007】
擬似反響路18がディジタルFIRフィルタで構成される場合には, そのフィルタ係数h^11(k)は室内反響路15のインパルス応答h 11(k)を直接模擬し
たものとなっている。従って、反響路15の変動に応じて必要なフィルタ係数修正の大きさは、反響路インパルス応答h 11(k)の変動量と一致する。そこで、
フィルタ係数修正動作における修正幅を表わすステップサイズ行列Aはインパ
ルス応答の経時変動特性で重み付けられている。一般に室内音場におけるインパルス応答変動量は減衰率λを用いた指数関数として表わされる。ステップサイズ行列Aの対角成分αl(l=1,2,…,L)は図2Aに示すようにlの増加に伴ってα0
からインパルス応答の指数減衰特性と同じ傾きで指数減衰し0に漸近する。このアルゴリズムは、人や物の移動により室内反響路のインパルス応答が変動する場合に、その変動量(インパルス応答の差分)はインパルス応答と同じ減衰率で指数減衰するという音響学的知見を利用したものである。変動の大きいインパルス応答の初期の係数は大きなステップで、変動の小さいインパルス応答の後期の係数は小さなステップで修正することにより、収束速度の大きな反響消去装置を提供することができる。
【0008】
反響消去装置を複数のDSP(Digital Signal Processor)チップで構成する場合には、図2Bに示すようにステップサイズαlの指数減衰曲線を階段状に近似し、各チップ毎に一定のαlを設定する。これにより、従来の射影アルゴリズムとほぼ同等の演算量と記憶容量でES射影アルゴリズムを実現することができる。ES射影アルゴリズムの詳細は特開平5-244043、文献S.Makino and Y.Kaneda,“Exponentially weighted step-size projection algorithm for acoustic echo cancellers”, Trans. IEICE Japan, vol.E75-A, no.11, pp.1500-1508, Nov.1992 に記載されている。
【0009】
以上より、h^(k)の更新式を従来装置で用いられているNLMSアルゴリ
ズムによる場合と(ES)射影アルゴリズムによる場合とで比較すると、NLMSアルゴリズムによる式(5)では、h^(k)は入力信号ベクトルx(k)の方向に
更新される。一方、(ES)射影アルゴリズムによる式(9), (16)では、h^(k
)は式(9), (16)の第4式右辺第2項以降を
v(k)=β1 x(k)+β2 x(k-1)+…+βp x(k-p+1) (18)
とおいて、このv(k)の方向、すなわち、入力信号ベクトルx(k)から過去の
入力信号ベクトルx(k-1), …, x(k-p+1)の全てとの相関を除去したベクト
ルの方向に更新される。つまり修正ベクトルv(k)が入力信号中の過去の信号
と同様のものをなるべく除去した信号となるようにβ1〜βpを決める。
【0010】
また、従来において、信号を複数の周波数帯域に分割し、その各周波数帯域について、NLMSアルゴリズム等により、反響路15の変動にもとづく擬似反響路18のフィルタ係数の逐次修正を行い、これら各周波数帯域ごとの残差を合成して出力するサブバンド法が知られている。これは例えば、米国特許第5,272,695、文献S.Gay and R.Mammone,“Fast converging subband acoustic echo cancellation using RAP on the WER DSP16A”, Proc.ICASSP90, pp.1141-1144, Apr.1990などに示されている。このサブバンド法により、周波数領域で信号の平坦化、いわゆる白色化が行われ、反響路の変動時の擬似反響路フィルタ係数の推定における収束速度が速くなる。
【0011】
通信会議システム用反響消去
さて、一般にI(>2)チャネルの再生系とJ(>1)チャネルの集音系とで構成される通信会議システムの場合の反響の消去は、図3に示すような構成により行われる。すなわち再生側の全Iチャネルと集音側の各1チャネルとの間にI入力1出力時系列信号を処理するIチャネル反響消去装置221,222,…, 22Jをそれぞれ接続した反響消去システム23として実現される。この場合システム全体でI×J個の反響路15ij(1<i<I, 1<j<J)が存在する。このシステムの構成単位である再生側の全Iチャネルと集音側の各1チャネルとの間に接続されるIチャネル反響消去装置221,222,…, 22Jについては、図1に示した反響消去装置14の構成を拡張して、図4に示すように構成される。これは例えば電子情報通信学会論文誌 '86/10 Vol. J69−A No. 10「多チャネル適応ディジタルフィルタ」に詳しく述べられている。ここで、集音側が第j集音チャネル(1<j<J)に接続されているIチャネル反響消去装置22jを考える。第jチャネルの集音器16jで集音される反響信号は、各再生チャネルの受話信号がそれぞれの反響路151j〜15Ijを経て集音側で全て加算されることにより得られるために、反響路推定をどの再生チャネルについても、統一的に同じ1つの残留反響信号ej(k)のみを評価して行なうための工夫が必要となる。まず、各再生チャネルの受話信号について、受話信号蓄積・ベクトル生成部(171,172,…, 17I)により、受話信号ベクトル
x 1(k)=[x1(k), x1(k-1),…, x1(k−L1+1)]T (19)
x 2(k)=[x2(k), x2(k-1),…, x2(k−L2+1)]T (20)
:
x I(k)=[xI(k), xI(k-1), …, xI(k−LI+1)]T (21)
を生成する。但し、L1,L2,…, LIはタップ数で、各反響路151j, 152j, …, 15Ijの残響時間に対応させて、あらかじめ設定する定数である。これらのベクトルをベクトル結合部24によって、
x(k)=[x 1 T(k), x 2 T(k),…, x I T(k)]T (22)
と結合する。また、反響路推定部19jにおいても、各再生チャネルと第j集音チャネルとの間のI個の反響路を模擬するための、各擬似反響路ベクトルh^
1j(k),h^2j(k),…, h^Ij(k)を結合して
として扱う。擬似反響路結合ベクトルh^j(k)の更新は、NLMSアルゴリズ
ムを用いた場合、
のように行なわれる。擬似反響路18jでは、内積演算
y^j(k)=h^j T(k)x(k) (25)
により、第j集音チャネルで集音された反響信号yj(k)に対する擬似反響信号y^j(k)を生成する。このように、各チャネル毎のベクトルを結合して1つのベクトルとして扱うことにより、基本的な処理の流れは、図1に示した1チャネル反響消去装置と同様となる。
【0012】
【発明が解決しようとする課題】
Iチャネルの再生系とJチャネルの集音系とで構成される通信会議システムに用いられる従来の反響消去システムの欠点のうち、この発明が解決しようとするものを、具体例に沿って説明する。
図5に示すように、A地点とB地点との間において、2チャネルで集音・再生を行なうステレオ音声会議装置に従来の反響消去システムを適用した場合、B地点側の反響路が全く変動しない場合においても、A地点側で発話者が移動あるいは交替を行なう度にその発話に対するB地点からの反響音が増大するという問題がある。これは、B地点側の反響消去システムにおいて反響路推定が正しく行なわれていないために起こる問題である。
【0013】
そこで、この問題を説明するために、B地点における反響消去システムの構成単位である2チャネル反響消去装置のうちの1つである第1集音チャネルに接続された装置22b1の動作に注目する。以下では、各再生チャネルの動作を明確化するために、従来技術の項で用いた結合ベクトルの各再生チャネルについての要素が明示された形で表記する。結合ベクトルを用いずに各再生チャネル毎のベクトルで表記するB地点側における2チャネル受話信号ベクトルをx 1(k)と
x 2(k)とする。各再生チャネルに対する真の反響路1511, 1521の反響路ベ
クトルをh 11(k),h 21(k)とすると、これら反響路1511, 1521を経て集
音される反響信号y1(k)は、
y1(k)=h 11 T(k)x 1(k)+h 21 T(k)x 2(k) (26)
と表される。一方、反響消去装置内で生成される擬似反響信号y^1(k)は、装置内で生成される擬似反響路h^11(k)とh^21(k)とを用いて、
y^1(k)=h^11 T(k)x 1(k)+h^21 T(k)x 2(k) (27)
と得られる。A地点から一人の話者が発話した場合には、この受話信号ベクトルx 1(k)とx 2(k)との間には、非常に強い相互相関がある。x 1(k)とx 2(
k)との間に、一定な強い相互相関がある場合には、
y^1(k)=y1(k) (28)
の解としての結合ベクトル[h^11 T(k),h^21 T(k)]は無数に存在し、x
1(k)とx 2(k)との相互相関に固有な部分空間H xを形成する。このため、N
LMSアルゴリズムのような一般的な逐次誤差最小化アルゴリズムを用いた場合に、[h^11 T(k),h^21 T(k)]は初期値から部分空間H xまでの距離が最
小となる点に収束し、一般に真値[h 11 T(k),h 21 T(k)]には収束しない。
【0014】
説明を簡単にするために、一定なスカラー値p1,p2と原信号ベクトルs(k)
により、受話信号ベクトルx 1(k), x 2(k)が、
x 1(k)=p1 s(k), x 2(k)=p2 s(k) (29)
と表される場合を考える。[h^11 T(k),h^21 T(k)]が存在し得る部分空間
H xは、
p1 h^11(k)+p2 h^21(k)=p1 h 11(k)+p2 h 21(k) (30)
を満たす図6A中の直線のように見なすことができ、初期値0から適応を開始した場合、収束点[h^11p T(k),h^21p T(k)]は、
のように得られる。このため、スカラー値p1とp2との間の比が変動した時点で、式(30)は満足されなくなるため、反響を消去できなくなり、反響が増大する。
【0015】
以上の例からも分かるように、一般にIチャネルの再生系とJチャネルの集音系とで構成される通信会議システムに、従来の反響消去装置を適用した場合は、各チャネルの受話信号間に相互相関がある場合に、反響路の推定が正しく行なわれないため、受話信号間の相互相関が変動する度に反響が増大する問題が生じる。
【0016】
この問題に対して、射影アルゴリズム等を用いて各チャネルの受話信号間の相互相関の変動の情報を抽出して利用することにより、真値への収束速度を改善できることが特願平7−50002で提案されている。しかし、全周波数帯域におけるインパルス応答に対して射影アルゴリズム等により修正を行なうこの提案の方法法では、まだまだ収束速度がおそいという問題があった。また、収束速度を向上させるためには、非常に高い射影次数が必要であり、演算量が増加するという問題がある。
【0017】
さて、前述のように射影アルゴリズム及びES射影アルゴリズムは入力信号の自己相関を取り除くことにより収束速度を上げているものである。一方、信号を複数の周波数帯域に分割して処理するサブバンド法は入力信号を周波数領域で平坦化(白色化)することにより収束速度を上げているものである。ここで、サブバンド構成にすることにより、相互相関に対してどのような影響があるかは、全く知られていない。むしろ、ステレオエコーキャンセラ(ステレオ反響消去装置)にサブバンド法を適用しても前述したステレオエコーキャンセラの問題は何ら解決されず、無駄だと思われてきた。そのため、これまでずっと、フルバンド法が用いられてきた。
この発明の目的は、従来の多チャネル音声通信会議用反響消去システムにおいて、以上に挙げたような問題を解決する多チャネル音声通信会議用の反響消去方法を提供することにある。
【0018】
【課題を解決するための手段】
この発明の方法を説明するに先立ち、この発明の一部の考え方のもとになる多チャネル受話信号間の相互相関が変動した場合における、多チャネル反響消去装置の反響路推定部の動作について考察する。
簡単のため、式(29)で表される2チャネルステレオ受話信号x 1(k)とx 2(
k)との相互相関が変動した場合として、式(29)のp1, p2がそれらとは比関係の異なるq1, q2に変化した例を考える。まず、〔h^11 T(k), h^21 T(k)〕
は、式(31), (32)の〔h^11p T(k), h^21p T(k)〕に収束する。次に、p1
とp2とがq1とq2に変化したとき、〔h^11p T(k), h^21p T(k)〕を“初
期値”として、そこから最短距離にある〔h^11q T(k), h^21q T(k)〕に収
束する。この動作は、図6Bに示すように幾何学的に解釈される。ここで、収束点(h^11q, h^21q)は収束点(h^11p, h^21p)から直線h^11q
1+h^21q2=h 11q1+h 21q2におろした垂線の交点になる。従って図
6Bにおけるフィルタ係数誤差ベクトルep, eqのノルムは、一般に‖ep‖>‖eq‖の関係にあることが自明である。このことは、一般に、受話信号間の相互相関の変化の度に、フィルタ係数誤差ベクトルのノルムが小さくなることを意味している。つまり、多チャネル受話信号間の相互相関の変動は、反響消去装置内において反響路推定が誤っている場合に、反響の増大をひき起こすが、無限回の変動が繰り返された末には、反響路推定部における真の反響路の推定を可能とする有効な情報であると捉えることができる。
【0019】
この発明による多チャネル反響消去方法には第1の方法と、第2の方法と、第3の方法とがある。
第1の方法としては、多チャネル反響消去方法において、サブバンド法を用いる。サブバンド法では間引きにより各サンプル間の時間間隔が広がっているため、相互相関の変動の情報が強調される。もちろん1チャネル反響消去方法にサブバンド法を用いた場合と同様に、入力信号の周波数領域での平坦化(白色化)も同時に行なわれる。これにより真値推定の収束速度が向上する。
【0020】
第2の方法としては、第1の方法において残留反響信号と各チャネルの受話信号とから修正ベクトルを求め、この修正ベクトルを用いて、各反響路のインパルス応答の推定を逐次修正するが、この発明では各チャネル間の過去の受話信号の相互相関に対する現在の受話信号間の相互相関の変動成分を抽出し、この変動成分を修正ベクトルとする。一例として、各チャネル間の相互相関の変動を強調する方法として(ES)射影アルゴリズムを用いる。このように、サブバンド多チャネルエコーキャンセラに(ES)射影アルゴリズムを用いると、サブバンド法では間引きによりタップ数が減っているため、サブバンドの相互相関の変動の強調にさらに(ES)射影アルゴリズムの相互相関の変動の強調効果が相乗され、予期せぬ効果を生み、少ない射影次数で非常に速い収束速度が得られる。さらに、射影次数を各周波数帯域で最適に設定することにより、少ない演算量で最大の効果を得ることができる。
【0021】
第3の方法としては、第1の方法において受話信号間の相互相関を、聴感上の違和感のない範囲で能動的に変動させる機能を付加して、各再生器で音響信号を再生し、この相互相関の変動が付加された受話信号を擬似反響路ベクトルの修正ベクトルの導出に利用する。さらに, 付加する機能を各周波数帯域で異なるものに最適に設定することにより、聴感上の劣化を最小限に押えることができる。
【0022】
【発明の実施の形態】
図7は前記第1の方法および第2の方法を用いたこの発明の反響消去方法を適用した反響消去装置の機能構成例を示し、図3と対応する部分には同一符号を付けてある。受話信号xi(k')はスピーカ12iで音声に再生されると共に周波数帯域分割部51i内で周波数帯域別のN個の実数信号xin(k)(n=0,1,…,N-1)に分割される。一方、スピーカ12iからの再生音声が反響路15ijを経てマイクロホン16jにより反響信号yj(k')としてピックアップされる。この反響信号yj(k')は周波数帯域分割部52j内で周波数帯域分割部51iと同一の分割特性で周波数帯域別のN個の実数信号yjn(k)に分割される。以降、jについてはすべて同様の処理となるので、簡単のためjは省略して記述する。
【0023】
周波数帯域分割部51iで分割されたそれぞれの周波数帯域には受話信号蓄積・ベクトル生成部17in、ベクトル結合部24n、擬似反響路18n、反響路推定部19nからなるIチャネル反響消去装置22nが設けられてあり、それぞれの擬似反響路18nは擬似反響信号y^n(k)を生成し、減算器21nでこれら擬似反響信号y^n(k)を対応する周波数帯域の反響信号yn(k)から差し引くことにより反響信号yn(k)は消去される。ここで擬似反響路18nは反響路15の経時変動に追従する必要があり、残差en(k)=yn(k)−y^n(k)が0に近づくように、射影アルゴリズムまたはES射影アルゴリズムを用いた推定部19nによって擬似反響路ベクトルh^n(k)が逐次推定され、擬似反響路18nの修正が行なわれることによっ
て、常に最適な反響消去が維持される。各周波数帯域の誤差信号(残差信号)en(k)は周波数帯域合成部53jで全周波数帯域の誤差信号e(k')に合成され、送信端4に出力される。
【0024】
図8は周波数帯域分割部51の内部の機能構成を示しており、A/D 変換器5により受話信号x(t)がサンプル値化され、帯域通過フィルタ20nにより帯域制限される。または、帯域通過フィルタ20nを用いる代わりに、文献Crochiere and Rabiner,“Multirate Digital Signal Processing”, Englewood Cliffs, NJ:Prentice-Hall, pp.52-56, 1983 に記載されているSSB(Single Side Band)法を用いても、帯域制限された実数信号xn(k')を得ることができる。帯域制限された実数信号xn(k')は間引き率Rで間引かれ実数信号xn(k)とされる。以下、分割された各周波数帯域のこの信号xn(k)をサブバンド受話信号と呼ぶ。反響信号y(k')に対する周波数帯域分割部52も図8の周波数帯域分割部51と同様に構成されており、これによって周波数帯域毎に分割された反響信号yn(k)をサブバンド反響信号と呼ぶ。
【0025】
サブバンド法では周波数帯域幅が狭まっており、間引きにより各サンプル間の時間間隔が広がっているため、相互相関の変動の情報が強調される。このため、真値への収束速度が早まる。
さらに、第2の方法では前述したように、各チャネルの受話信号間の相互相関の変動の情報を抽出して利用する。図9に各周波数帯域における反響消去装置22jを示す。図4と対応する部分に同一符号を付けてある。この発明では各反響消去装置221〜22jに相互相関変動成分抽出部31を設ける。相互相関変動成分抽出部31では、受話信号結合ベクトルx(k)と残留反響信号ej(k)の情報
を蓄積し、自己相関除去ベクトル、つまり修正ベクトルv(k)を生成する。こ
のベクトルv(k)は、x(k)が相互相関の高い多チャネル信号ベクトルの結合
ベクトルである場合、各チャネル間の相互相関の変わらない成分は除去されて、相互相関の変動成分が強調されたベクトルとなる。修正部32では、現在の推定反響路ベクトルh^j(k)に対し、修正ベクトルv(k)を用いて擬似反響路結
合ベクトルh^j(k)を更新し、これを擬似反響路18jへ出力する。
【0026】
従来装置の構成においては、各チャネル間の受話信号ベクトルを結合して扱う。このため、各チャネルの受話信号間の相互相関は、受話信号結合ベクトルと過去の受話信号結合ベクトルとの相関関係である自己相関の中に現れる。よって、各チャネルの受話信号について、相互相関の変動していない成分を除去して、相互相関の変動した成分を抽出するためには、受話信号結合ベクトルの自己相関を小さくする、つまり受話信号結合ベクトルの自己相関を除去する手法を用いればよい。このような手法として、1チャネル反響消去装置において、射影アルゴリズム、RLSアルゴリズム、あるいは線形予測を用いる方法等が提案されている。これらのうち、線形予測による手法は、音声信号を線形予測して、その線形予測係数を用いて音声信号の自己相関を除去するものである。RLSアルゴリズムは、射影アルゴリズムの特殊な場合とみなすことができる。そこで、この発明では、最も現実的な手法の例として、射影アルゴリズムを用いた場合における受話信号結合ベクトルの自己相関除去手法を説明する。
【0027】
図10は第n番目の周波数帯域の推定部19nの内部の機能構成例としてpn次のES射影アルゴリズム(A=Iの時には射影アルゴリズム)を用いた場合を
示す。
受話信号記憶部31はサブバンド受話信号xn(k)が与えられる毎にそれぞれ連続するLn個ずつのサブバンド受話信号からなるpn個のサブバンド信号ベクトル
を生成し、更にこれらベクトルからサブバンド受話信号行列
X n(k)= [x n(k),x n(k-1),…,x n(k-pn+1)]
を生成する。Lnはn番目の周波数帯域の擬似反響路18nを構成するFIR フィルタのタップ数である。なおこれらの式中のx n(k),x n(k-1),…はそれぞれ式
(22)で与えられるベクトル結合したベクトルである。つまり受話信号記憶部31には1乃至Iチャネルのサブバンド受話信号x n(k)がベクトル結合した
X n(k)が入力される。
【0028】
ただし、n番目の周波数帯域ではpn次の射影アルゴリズム又はES射影アルゴリズムを実行するものとする。ステップサイズ行列記憶部32には第1のステップサイズ行列A n=diag[αn1,αn2,…,αnLn] が記憶される。ステップサイズ
行列A nは対応する周波数帯域におけるインパルス応答の変動特性で重み付け
られている。
【0029】
一般的な室内を対象とする場合には、第n番目の周波数帯域におけるインパルス応答変動量は減衰率γnを用いた指数関数として表わされる。自己相関演算部33では第1のステップサイズ行列A nで重み付けたサブバンド受話信号行列
X n(k)の自己相関行列X n(k)T A n X n(k)が演算される。この演算された
自己相関行列と、残差記憶部34からの残差en(k)とがβ n(k)演算部35に供
給されて次式のpn元連立一次方程式
[X n(k)T A n X n(k)]β n(k)= e n(k) (33)
を解くことにより、定数β n(k)を求める。逆行列演算における不安定性を回避
するために、小さな正の定数δnを用いて式(33)の代わりに次式
[X n(k)T A n X n(k)+δn I]β n(k)= e n(k) (33)'
を使ってもよい。ただし、Iは単位行列である。ステップサイズ記憶部36か
らの第2のステップサイズμnと、A n, X n(k), β n(k)とが修正ベクトル
演算部37に供給されて修正ベクトル
μn A n X n(k)β n(k) (34)
が演算され、その修正ベクトルは加算器38へ供給されてタップ係数記憶部39からのLn個の要素からなる擬似反響路ベクトルh^n(k)に加算されて修正さ
れた擬似反響路ベクトルh^n(k+1)が得られる。演算結果h^n(k+1)は擬似
反響路18nへ出力されると同時に、タップ係数記憶部39の値を更新する。
【0030】
以上の操作をサブバンド受話信号xn(k)が与えられる毎に繰り返すことにより、擬似反響路18nは次に示す(35)式に従って逐次修正され、擬似反響路ベクトルh^n(k)は真の擬似反響路15の第n周波数帯域におけるインパルス応答
h n(k)に近づいてゆく。
ここで特筆すべきことは、相互相関の変動は各周波数帯域によって異なり、一般に、低域では小さく高域では大きいという点である。この点を考慮すれば、相互相関の変動の小さい周波数帯域では射影次数を大きくしてして収束速度を早め、相互相関の変動の大きい周波数帯域では射影次数を小さくすることにより、少ない演算量で全体としての収束速度を均一にすることができる。
【0031】
もう一つ特筆すべきことは、各帯域のタップ数Lnは間引きにより少なくなっているという点である。(ES)射影アルゴリズムでは、次数pと収束特性の関係はタップ数Lにより影響を受ける。すなわち、タップ数Lが大きい場合には、完全に収束速度が飽和するためには次数pを大きくする必要があり、タップ数Lが小さい場合には、小さな次数pで完全に収束速度が飽和する。例えば、タップ数Lが1000程度の場合には、完全に収束速度が飽和するためにはpが50〜100次程度必要であるのに対して、帯域分割と間引きによりタップ数Lnが64程度に小さくなった場合には、16〜32次程度で十分となる。
【0032】
このように、サブバンド多チャネルエコーキャンセラにさらに(ES)射影アルゴリズムを用いると、サブバンド法では間引きによりタップ数が減っているため、サブバンドの相互相関の変動の強調にさらに(ES)射影アルゴリズムの相互相関の変動の強調効果が相乗され、予期せぬ効果を生み、少ない射影次数で非常に速い収束速度が得られる。さらに、射影次数pを各周波数帯域で最適に設定することにより、少ない演算量で最大の効果を得ることができる。
【0033】
図11に示す一般的な音声の周波数スペクトルの例では、受話信号を複数の周波数帯域F0,F1,…,FN-1に分割した場合には、各帯域でそのスペクトル形状が異なる。その結果、各帯域の収束速度が飽和するために必要となる(ES)射影アルゴリズムの次数も異なる。一例を挙げれば、図11中の帯域Fg では、スペクトルが複雑に変化しているため、必要な(ES)射影アルゴリズムの次数は大きく、帯域Fnではスペクトルが平坦に近いため、必要な(ES)射影アルゴリズムの次数は小さい。
【0034】
従って、この発明では各帯域においてその帯域に適した、つまり収束速度が最も速く、かつ次数が最も少ないような次数を決める。このため、図7中に示す次数決定制御部54を設ける。その内部の機能構成例を図12に示す。それぞれの周波数帯域における(ES)射影アルゴリズムの次数pnは、次数決定制御部54で次のように決定する。
【0035】
方法1:次数pnを変えながらynとenとから収束の様子を観測する。この場合の次数決定制御部54は図12に示すように、それぞれの周波数帯域毎に設けられた次数決定部540〜54N-1から構成されている。各周波数帯域において、例えば次数pn=1と決め、サブバンド受話信号xn(k)が入力される毎に得られるサブバンド反響信号ynと残差信号enの比(echo return loss enhancement:ERLE)rpn=20log10(yn/en)である反響消去量を求めることにより例えば図13に模式的に示すような反響消去量ERLEの収束曲線が得られる。同様にpnを順次大にした場合についても同様のERLE収束曲線が得られる。そこで、方法1では、時間k=0 からpn次の(ES)射影アルゴリズムによる反響消去処理を開始し、次数決定部54nにおいて予め決めた時点k=kKでのERLEの値rsを求めることを単調増加する一連の次数pn=ps,s=1,2,…, 例えばps+1=ps+1,p1=1、について順次繰り返し実行し、各次数pn=psについて反響消去量rsを求める毎にその前の次数pn=ps-1で求めた反響消去量rs-1からの増加の、次数の増加(ps-ps-1)に対する割合R=(rs-rs-1)/(ps-ps-1)を求め、この値が予め決めた閾値Rthより最初に小さくなる次数pn=psを決定する。
【0036】
図14は第1の次数決定方法による具体的な次数決定手順の例を示す。次数決定部54nはステップS1で整数パラメータs, kをそれぞれ1と0に初期設定し、ステップS2で次数pn=psとし、対応する周波数帯域の推定部19nに与える。ステップS3でサブバンド受話信号xn(k)が入力されると、ステップS4で受話信号ベクトルx n(k)を擬似反響路18nに入力して擬似反響信号y^n(k)=h
^n(k)T x n(k)を得て、その擬似反響信号y^n(k)の、サブバンド反響信号yn(k)
に対する推定誤差en(k)=yn(k)-y^n(k)を減算器21nで求める。更に、推定部19nにおいて時点kを最新時点とするpn+Ln個のサブバンド受話信号xn(k),xn(k-1),…,xn(k-pn-Ln+2)とpn個の推定誤差en(k),en(k-1),…,en(k-pn+1)とを使ってpn=ps次の(ES)射影アルゴリズムにより擬似反響路ベクトルh^n(k+1)を推定
する。更に、推定したh^n(k+1)を対応する擬似反響路18nに設定する。即
ち、ステップS4ではサブバンド受話信号xn(k)が入力される毎に、図7の反響消去装置における第n番目の周波数帯域に係わる部分によりpn=ps次の(ES)射影アルゴリズムに基づく反響消去処理が行われる。
【0037】
ステップS5で時間kが予め決めた時間kKに達したか判定し、達していなければステップS6でkを1歩進してステップS3に戻り、ステップS3, S4で同様のpn次(ES)射影アルゴリズムによる反響消去処理を実行する。時間kがkKに達するまで反響消去処理を繰り返し実行し、k=kKとなった時点でステップS7において次数決定部54nはその時のサブバンド反響信号ys=yn(kK)と推定誤差es=en(kK)から反響消去量rs=20log10(ys/es)を求め、更に、次数増加(ps-ps-1)に対する反響消去量の増加の割合(即ち収束速度の飽和の程度)R=(rs-rs-1)/(ps-ps-1)を求める。ただしs=1の場合R=∞とする。次にステップS8で収束速度の飽和の程度Rが予め決めた値Rthより小さくなったか判定し、Rthより小となっていない(収束速度が飽和していない)場合はステップS10でkを0にリセットし、sを1歩進してステップS2に戻り、次数pn=psについて再び前述と同様の処理ステップS3〜S9を実行する。ステップS8でR<Rthを満足すれば、収束速度は十分飽和したと判定し、ステップS10でその時の次数psをn番目の周波数帯域における(ES)射影アルゴリズムの次数pnとして決定し、対応する推定部19nに出力する。
【0038】
方法2:あるいは、図13において、一定消去量rKに達する時間が次数pnと共に短くなることを利用して、反響消去量の飽和値に対し、予め決めた割合の閾値rK、例えば30dBを決める。単調増加する次数pn=ps, s=1,2,… についてそれぞれ時点k=0からpn=ps次の(ES)射影アルゴリズムによる反響消去動作を開始し、サブバンド受話信号xn(k)が入力される毎に図12の次数決定部54nで反響消去量rs=20log10(ys/es)を求めることを繰り返し、rs>rKとなる時点ksを得る。前の次数ps-1で得た時点ks-1との差(ks-1-ks)の、次数の増加(ps-ps-1)に対する割合R=(ks-1-ks)/(ps-ps-1)が予め決めた閾値Rthより小さくなる次数pn=psを決定する。この場合の次数決定処理手順を図15に示す。
【0039】
図14の場合と同様にステップS1で整数パラメータs, kをそれぞれ1と0に初期設定し、ステップS2で次数pnをpsに決め、ステップS3でサブバンド受話信号xn(k)が入力される毎にステップS4でps次の(ES)射影アルゴリズムによる反響消去処理を行う。この第2の方法ではステップS5で次数決定部54nはサブバンド反響信号ys=yn(k)と推定誤差es=en(k)から反響消去量rs=20log10(ys/es)を求め、ステップS6で反響消去量rsが予め決めた値rKより大となったか判定し、rs>rKとなっていなければステップS7でkを1歩進してステップS3に戻り、以下ステップS4, S5, S6で同様の処理を繰り返す。rs>rKとなるとステップS8で現在の時点kを次数pn=psで消去量rsがrKに達した時点ksと判定し、次数pn=ps-1の場合にrKに達した時点ks-1からの時間差(ks-1-ks)の、次数変化(ps-ps-1)に対する割合、即ち収束速度の飽和の程度をR=(ks-1-ks)/(ps-ps-1)として求める。ただし、s=1の時はR=∞とする。次にステップS9でRが予め決めた閾値Rthより小さくなったか判定し、なっていなければステップS10でkを0にリセットし、sを1歩進してステップS2に戻り、ステップS2〜S9の処理を繰り返す。ステップS9でR<Rthと判定された場合は、ステップS11でその時点での次数pn=psをn番目の周波数帯域における(ES)射影アルゴリズムの次数pnと決定し、対応する推定部19nに設定する。
【0040】
上述の第1及び第2のいずれの方法の場合も、n番目の周波数帯域において、単調増加させる(ES)射影アルゴリズムの次数pn=ps, s=1,2,…として、例えばps+1=ps+d,dは1以上の一定の整数、としてもよいし、ps+1=ps+sdとしてもよいし、あるいはps+1=cps、cは2以上の一定整数、としてもよい。その他、漸次増加する関数であればどのような単調増加関数であってもよい。
【0041】
方法3:pn次の(ES)射影アルゴリズムは、入力信号xnを(pn-1)次の白色化FIRフィルタ(即ちタップ数pnの線形予測フィルタ)に通すことと等価である。そこで、図16に示すように次数決定制御部54の各次数決定部54n(n=0,1,…,N-1)を(pn-1)次の白色化FIRフィルタ54Fnと予測分析部54Pnで構成する。予測分析部54Pnは白色化FIR フィルタ54Fnに入力される一連のLn個の各サブバンド受話信号xn(k-j),j=0,1,…,Ln-1,に対し、Ln個の予測誤差εn(k-j),j=0,1,…,Ln-1, の2乗和
Σεn 2(k-j)= Σ{xn(k-j)-Σanixn(k-j-i)}2 (40)
ただし、左辺及び右辺第1のΣはj=0〜Lnについて、右辺第2のΣはi=1〜pn-1についての総和とする、
が最小となる予測係数ani,i=1,2,…,pn-1を求める(一般に線形予測分析と呼ばれ、周知の技術である)。これによって得られた予測係数ani をフィルタ54Fnに設定して出力に得られるLn個の予測誤差εn(k-j), j=0,1,…,Ln-1をベクトル
ε n(k)=[εn(k),εn(k-1),…,εn(k-Ln+1)] (41)
で表すと、予測誤差ベクトルε n(k)は次式
ε n(k)=X n(k)a n(k) (42)
と表すことができる。ただし
a n(k)=[1,-an1,-an2,…,-anpn-1]T (43)
式(42)はサブバンド受話信号行列X n(k)に対する白色化処理を表しており、こ
の白色化処理により得られた予測誤差ベクトルε n(k)の共分散行列
Q n(k)=ε n(k)ε n(k)T (44)
はLn個の固有値λn0,λn1,…,λnLn-1を有している。それらの固有値のうちの最大値λmaxと最小値λminの比Cs=λmax/λminはサブバンド信号行列X n(k)に対
する白色化の程度を表している。この比Csの値が小さいほど白色化の程度は大きい。完全な白色化が達成されると比Cs=1となる。そこで、射影次数pnをpn=ps, s=1,2,… と順次増加させたときの比Cs, s=1,2,…を求め、比Csが予め決めた閾値Cthより小さくなったとき、あるいはCsの変化量ΔCs=Cs-1-Csが予め決めた閾値ΔCthより小さくなったときの最初の次数psを第n番目の周波数帯域における射影の次数pnと決める。
【0042】
図17は第3の方法を使って第n番目の周波数帯域における(ES)射影アルゴリズムの次数pnを決める手順を示す。ステップS1で整数パラメータsを1とし、ステップS2で次数pnをpsに設定し、ステップS3で受話信号行列X n(k)を
構成するのに必要な数のサブバンド受話信号xn(k),xn(k-1),…,xn(k-pn-Ln+2)を入力する。次にステップS4で入力されたサブバンド受話信号に対し式(40)による線形予測分析を行って予測係数ani,i=1,2,…,pn-1を求め、ステップS5で式(42)によりサブバンド受話信号に対し白色化処理を行って推定誤差ベクトルε n
(k)を求める。次にステップS6で式(44)で表される推定誤差ベクトルε n(k)
の共分散行列Q n(k)を求め、更にその共分散行列のLn個の固有値λ0,λ1,…,
λLn-1を求める。ステップS7で固有値の最大値と最小値の比Cs=λmax/λminを計算し、ステップS8で比Csが閾値Cthより小となったか判定し、小となっていなければステップS9でsを1歩進してステップS2に戻り、再びステップS2〜S8による同様の処理を繰り返す。ステップS8で比Csが閾値Cthより小と判定されると、ステップS10で次数psを第n番目の周波数帯域における(ES)射影アルゴリズムの次数pnと決定し、出力する。
【0043】
上述のステップS8では比Csを閾値Cthと比較した場合を示したが、前述のようにステップS7で前回の次数ps-1で得た比Cs-1と今回の次数psで得た比Csとの差分ΔCs=Cs-1-Csを求め、ステップS8でこの差分ΔCsを予め決めた閾値ΔCthと比較し、ΔCsがΔCthより小になったらステップS10でその時の次数pn=psを出力してもよい。次数pn=psを単調に増加させる関数としては前述の第1及び第2の次数決定方法で説明したと同様な単調増加関数を使うことができる。
【0044】
例えば工場出荷時に、方法1〜3のいずれかを用いていろいろな音声について(ES)射影アルゴリズムの次数pnを決定して推定部19nに予め設定しておいてもよい。なお、方法1〜3のいずれかを用いていろいろな音声、いろいろな分割数N、タップ数Lnについて(ES)射影アルゴリズムの次数pnを決定してその標準的な値を予めROM化しておき、使用者が、その反響消去装置の所望の帯域分割数Nと対応して前記ROMから推定部19nにpnを設定するようにしてもよい。
【0045】
DSPなどを用いて反響消去装置を構成する場合には、実時間処理の制約上、(ES)射影アルゴリズムの次数を完全に収束速度が飽和するレベルまで上げることができない場合も多い。このような場合には与えられた演算時間の中で装置全体としての性能が最高となるように各帯域の次数pnを決定する。
以上のように各周波数帯域における(ES)射影アルゴリズムの次数を、複数の周波数帯域のそれぞれに対してその装置に許される最も好ましい次数に設定することにより、(ES)射影アルゴリズムの次数を上げることによる相互相関の変動の強調と入力信号の白色化の効果を十分に引き出すことができ、収束速度の大きな反響消去装置を得ることができる。
【0046】
図18にこの発明の他の実施例を適用した機能構成を示す。この実施例は図7に示した実施例とは周波数帯域分割方法が異なっており、信号Xn(k), Yn(k), En(k)および擬似反響路ベクトルH^n(k)が複素数である点が異なる。同様な周波
数帯域分割方法を用いた反響消去装置は前記S.Gay and R.Mammone の文献に示されている。受話信号x(k')は周波数帯域分割部61iで周波数帯域別のN個のサブバンド複素数信号Xin(k)(n=0,1,…,N-1)に分割される。同様に反響信号yj(k')は周波数帯域分割部62jで周波数帯域別のN個のサブバンド複素数信号Yjn(k)に分割される。以降、jについてはすべて同様の処理となるので、簡単のためjは省略して記述する。
【0047】
それぞれの周波数帯域には受話信号蓄積ベクトル生成部117n、ベクトル結合部124n、擬似反響路65nが設けられてあり、擬似反響路65nからの擬似反響信号Y^n(k)をサブバンド反響信号Yn(k)から減算器66nで差し引くことによりサブバンド反響信号Yn(k)は消去される。
擬似反響路(複素FIRフィルタ)65nは擬似反響路15の経時変動に追従する必要があり、残差En(k)=Yn(k)−Y^n(k)が0に近づくように、複素射影アルゴリズム又は複素ES射影アルゴリズムを用いた推定部67nによって逐次推定され、擬似反響路65nの修正が行なわれることによって、常に最適な反響消去が維持される。
各帯域の誤差信号En(k)は周波数帯域合成部63jで全帯域信号e(k')に合成される。この分割・合成過程はN点FFTを用いて効率的に実行できる。
【0048】
図19に周波数帯域分割部61の内部を示しており、A/D変換器5により受話信号x(t)がサンプル値化され、その受話信号x(k')に乗算器68nでWN -k'n=exp{-j2πk'n/N}が乗算され、通過周波数帯域幅-π/N〜π/Nの低域通過フィルタ71nを用いて帯域制限され、N個の周波数帯域に分割される。周波数帯域制限された信号Xn(k')は間引き率Rで間引かれ、サブバンド信号Xn(k)となる。全分割周波数帯域のサブバンド信号Xo(k)〜XN-1(k)は短時間スペクトルに対応する。N=16の場合の周波数帯域を図20に示す。16帯域信号のうち0と8が実数、残りは複素数である。帯域8に対して対称なもの(例えば7と9、等)は複素共役の関係にあり、全部で9つ(2実数、7複素数)の周波数帯域信号があれば全周波数帯域信号を合成することができる。
【0049】
図21は第n番目の周波数帯域の推定部67nの内部の一例としてpn次の複素ES射影アルゴリズム(A=Iの時には複素射影アルゴリズム)を用いた構成
例を示したものである。
サブバンド受話信号Xn(k)は受話信号記憶部75でサブバンド受話信号行列
X n(k)とされる。ここで受話信号記憶部75には1乃至Iチャネルの各サブバ
ンド受話信号x n(k)が図9中のベクトル結合部24と同様の処理でベクトル結
合されたX n(k)が入力される。ステップサイズ行列記憶部76には第1のステ
ップサイズ行列A nが記憶される。ステップサイズ行列A nは対応する周波数
帯域におけるインパルス応答の変動特性で重み付けられている。一般的な室内を対象とする場合には、第n番目の周波数帯域におけるインパルス応答変動量は減衰率γnを用いた指数関数として表わされる。自己相関演算部77では第1のステップサイズ行列A nで重み付けた受話信号行列X n(k)の自己相関行列X n
(k)T A n X n *(k)が演算される。ただし*は複素共役を表す。この演算された
自己相関行列と、残差記憶部78からの残差E n(k)は、β n(k)演算部79に
供給されてpn元連立一次方程式
[X n(k)T A n X n *(k)]β n(k)=E n(k) (45)
を解くことにより、定数β n(k)を求める。逆行列演算における不安定性を回避
するために、小さな正の定数δnを用いて式(45)の代わりに次式
[X n(k)T A n X n *(k)+δn I]β n(k)=E n(k) (45)'
を使ってもよい。ただし、Iは単位行列である。
【0050】
A n, X n(k), β n(k)とステップサイズ記憶部81からのμnとが修正情
報生成部82に供給されて次式
μn A n X n *(k)β n(k) (46)
が演算され、その出力は加算器83へ供給されてタップ係数記憶部84からの
H^n(k)に加算されてH^n(k+1)が得られる。演算結果H^n(k+1)は擬似反
響路65nへ出力されると同時に、タップ係数記憶部84の値を更新する。
【0051】
以上の操作により、擬似反響路65nは次の式(47)に従って逐次修正され、擬似反響路ベクトルH^n(k)は真の反響路のインパルス応答H n(k)に近づいて
ゆく。
*:複素共役ここで図7の実施例と同様に、各帯域のタップ数は間引きにより少なくなっている。また、各帯域で相互相関の変動量やスペクトルが異なる。その結果、各帯域の収束速度が飽和するために必要となる複素(ES)射影アルゴリズムの最小の次数も異なる。各周波数帯域における複素(ES)射影アルゴリズムの次数pnの決定方法は図7の実施例で示した3つの方法と同様な方法が可能である。
【0052】
図22は図14中の第1又は第2の次数決定方法を適用する場合の次数決定制御部64の機能構成例を示したもので, それぞれの周波数帯域に対応して次数決定部640〜64N-1が設けられる。各次数決定部64nはサブバンド複素受話信号Xn(k)、サブバンド複素反響信号Yn(k)、複素誤差信号En(k)が与えられ、複素(ES)射影アルゴリズムの次数pnを図12の場合と同様に次のように決定する。
【0053】
第1及び第2の方法では、次数pnをpn=ps, s=1,2,…と変えながら、Yn=YsとEn=Esとから求めた反響消去量rs=20log10(Ys/Es)の収束の様子を観測する。次数psを大きくするにつれ収束速度が大となり、図13で示したと同様に収束速度が飽和値に近づく。収束速度が充分飽和に近づいたと判定された時の次数pn=psを出力して推定部67nに設定する。即ち、図13で説明したと同様に第1の方法では各次数毎に反響消去処理開始後、予め決めた時点kKでの反響消去量rs=20log10(Ys/Es)を求め、前回の次数ps-1で求めた反響消去量rs-1からの変化の割合 Rs=(rs-rs-1)/(ps-ps-1)が閾値Rthより小さくなったときの次数psを複素(ES)射影アルゴリズムの次数pnと決める。第2の方法では各次数毎に反響消去処理開始後、反響消去量 rs=20log10(Ys/Es)が予め決めた値rKに達した時点までの時間ksを求め、前回の次数ps-1について求めた時間ks-1からの短縮率(ks-1-ks)/(ps-ps-1)が閾値Rthより小さくなったときの次数psを複素(ES)射影アルゴリズムの次数pnと決める。
【0054】
図23は図14において第3の次数決定方法を適用する場合の次数決定制御部64の構成例を示し、図16の場合と同様に各帯域の次数決定部64n(n=0,1,…,N-1)は(pn-1)次の白色化FIR フィルタ64Fnと予測分析部64Pnとから構成されている。サブバンド複素信号Xnを(pn-1)次の白色化 FIRフィルタ64Fnに通し、予測分析部64Pnにより線形予測分析を行って予測係数ani,i=1,2,…,pn-1を決め、その時の推定誤差ベクトルε n(k)をX n(k)a n(k)から計算する。更に、推定
誤差ベクトルの共分散行列Q n(k)=ε n(k)ε n(k)Tを求め、そのLn個の固有
値を求め、それらの固有値の最大値と最小値の比Cs=λmax/λminを求める。この比Csが次数psの増加につれ小さくなり、予め決めた値Cthより小となったときの次数psを複素(ES)射影アルゴリズムの次数pnと決める。
【0055】
図7の実施例と同様に、いろいろな音声およびいろいろな分割数N、タップ数Lnについて次数pnを決定し、その標準的な値を予めROM化しておき、例えば電源立ち上げの際にこのROMから推定部6nに与えるようにしてもよい。また、DSPなどを用いて反響消去装置を構成する場合には、実時間処理の制約上、複素(ES)射影アルゴリズムの次数を完全に収束速度が飽和するレベルまで上げることができない場合も多い。このような場合には与えられた演算時間の中で装置全体としての性能が最高となるように各帯域の次数pnを決定する。
【0056】
以上のように各周波数帯域における複素(ES)射影アルゴリズムの次数を、複数の周波数帯域のそれぞれに対して最適とされる次数に設定することにより、複素(ES)射影アルゴリズムの次数を上げることによる相互相関の変動量の強調と入力信号の白色化の効果を十分に引き出すことができ、収束速度の大きな反響消去装置を得ることができる。
【0057】
また、室内の残響時間は一般に、低い周波数帯域では長く、高い周波数帯域では短い点から、擬似反響路65nのタップ数Lnを低い周波数帯域に対しては大に、低い周波数帯域に対しては小にし、従って射影次数pnは低い周波数帯域に対しては大きく、高い周波数帯域に対しては小さくするのが好ましい。例えば音声信号では老若男女によらず一般に図11に示したように低い周波数対により多くのエネルギーが集中していることを利用して、低い周波数帯域ではタップ数Lnを大きく、従って射影の次数pnも十分大きな値とするが、高い周波数帯域ではタップ数Lnを小さくし、かつ射影の次数pnも小さくすることが可能である。更に、人間の聴覚感度は一般に低い周波数帯域では高く、高い周波数帯域では低いことからも、低い周波数帯域に対してはタップ数Lnを大きく、従って射影次数pnも大きく、低い周波数帯域に対してはタップ数Lnを小さく、従って射影次数pnも小さくしてよい。この様に、高い周波数帯域のタップ数Lnを小さくすることにより、あるいは更に射影次数pnを小さくすることによって全体としての演算量を削減することができる。音声の場合、例えば周波数帯域を32帯域に分割し、周波数の低い方から1番目と2番目の周波数帯域の射影次数をそれぞれ16とし、3及び4番目の周波数帯域の射影次数を8とし、5〜8番目の周波数帯域の射影次数を4とし、9〜16番目の周波数帯域の射影次数を2とし、それより高い周波数帯域の次数を全て1にする。
【0058】
なお、上述において周波数帯域の分割数Nは32〜64程度が現実的と考えられ、分割数が多過ぎると遅延が大となる。また分割方法は例えば図11に示すように非等分割にする場合に限らず、等分割に分けてもよい。各帯域の最適次数pnの決定は方法1の場合は誤りが生じない。
図24にこの発明の第3の方法の実施例を示す。図7と対応する部分に同一符号を付けてある。この発明では、相互相関変動付加部24が設けられる。相互相関変動付加部24で、各チャネルの受話信号x1(k), x2(k), …, xI(k)は、それらの相互相関の変動を能動的に付加されて、x ̄1(k), x ̄2(k), …, x ̄I(k)に変換され、これらx ̄1(k), x ̄2(k), …, x ̄I(k)を、スピーカ121, 122, …, 12Iから再生するとともに、反響消去システム23への入力とし、擬似反響路ベクトルの修正ベクトルの導出に利用する。
【0059】
相互相関変動付加部24の実現については、変換された受話信号x ̄1(k), x ̄2(k), …, x ̄I(k)が音響信号として再生されるときに、聴覚上の品質が損なわれないように注意する必要がある。ここで注目すべきことは、人間の聴覚特性は各周波数帯域で異なり、各実現方法による聴覚上の品質劣化は各周波数帯域で異なるという点である。
【0060】
各チャネルの受話信号間の相互相関に変動付加する方法としては、▲1▼フィルタ処理による実現方法と、▲2▼信号乗算処理による実現方法と、▲3▼信号加算処理による実現方法と、▲4▼ピッチシフト処理による実現方法とがある。以下、順に▲1▼〜▲4▼各実現方法の具体例を説明する。
▲1▼ フィルタ処理による実現方法
各チャネルの受話信号x1(k), x2(k), …, xI(k)を、異なる時変特性をもつ時変フィルタに印加、すなわち、これら時変フィルタのインパルス応答f1(k), f2(k), …, fI(k), 畳み込み演算を*で表すとして、
x ̄1(k)=f1(k)*x1(k)
x ̄2(k)=f2(k)*x2(k)
x ̄I(k)=fI(k)*xI(k) (48)
を満足するx ̄1(k), x ̄2(k), …, x ̄I(k)に変換することにより、各チャネルの受話信号間の相互相関への変動の付加を実現する。
【0061】
人の聴覚による音声知覚においては、主として振幅情報が重要であり、位相情報は余り重要はでないと言われており、このことは、例えば文献「古井, “ディジタル音声処理”, 東海大学出版会」で述べられている。そこで、入力信号の周波数に対して振幅特性が平坦な特性を保つように、時変フィルタを実現すれば、聴覚的にフィルタ処理による影響が小さくなると考えられる。
【0062】
この目的に適した構成として、全域通過伝達関数をもつIIRフィルタがある。このフィルタの伝達関数F(z)は、次数をKとして、
と一般に表すことができる。標本化間隔をT、角周波数をωとすると、振幅特性
|F(exp(jωT))|
式は、フィルタ係数a1, a2, …, aKがどのように定められても分母と分子で相殺し合い、ωによらない一定利得Aをとる。一方、位相特性はフィルタ係数a1, a2, …, aKにより異なる。このa1, a2, …, aKを時間とともに変化させた場合には、時変の位相特性を有することになるが、振幅特性への影響はない。
【0063】
さて、ここでは構成が簡単で実現の容易な2次の全域通過フィルタを考える。2次の全域通過フィルタの構成は、例えば文献「辻井, “ディジタル信号処理の基礎”, 電子情報通信学会」に詳しく述べられている。式(49)について、簡単のため利得A=1とすると、2次の全域通過伝達関数F(z)は、
F(z)=(z-2−a1z-1+a2)/(1−a1z-1+a2z-2) (50)
となる。この伝達関数をもとに、ラチス型で構成すると、式(50)の係数は
a1=−r(1+r2) (51)
a2=r2 (52)
と置き換えられ
と表される。この発明において、このようなラチス型で構成する利点は、r1,r2が独立に意味をもつパラメータとして扱えることにある。式(53)のF(z)の群遅延特性がピークを持つ場合、ピーク周波数はほぼr1のみに依存して決まり、ピークの急峻さはほぼr2のみに依存して決まる。つまり、フィルタ係数を時間とともに変化させる場合、すべての係数を変化させなくとも、例えば、群遅延特性のピーク周波数に対応したr1を再生チャネルごとに異なる変化をさせるだけで、効果的な時変特性を得ることができ、構成が簡単になる。これらの時変特性をチャネルごとに異なるものとすれば、もとの多チャネル受話信号に対して、それらの間の相互相関を変動させる効果がある。
【0064】
▲2▼ 信号乗算処理による実現方法
各チャネルの受話信号x1(k), x2(k), …, xI(k)を、異なる関数g1(k), g2(k), …, gI(k)と乗算、すなわち、
x ̄1(k)=g1(k)・x1(k)
x ̄2(k)=g2(k)・x2(k)
・
・
・ (54)
x ̄I(k)=gI(k)・xI(k)
を満足するx ̄1(k), x ̄2(k), …, x ̄I(k)に変換することにより、各チャネルの受話信号間の相互相関への変動の付加を実現する。
【0065】
式(54)において、もし関数g1(k), g2(k), …, gI(k)が正負の値をとる場合には、その変動の大きさによらず受話信号の品質は大幅に劣化すると考えられる。そこで、関数g1(k), g2(k), …, gI(k)は、常に正、あるいは常に負の値のみをとり、また、変動の幅も制限されていることが望ましい。
そこで、関数g1(k), g2(k), …, gI(k)が以下のような形である場合を考える。最大振幅が1に正規化された異なる関数ω1(k), ω2(k), …, ωI(k)と、スカラー値k1, k2, …, kI(0<‖k1‖, ‖k2‖, …, ‖kI‖<1), および利得係数A1, A2, …, AIを用いて、g1(k), g2(k), …, gI(k)を、
g1(k)=A1(1+k1ω1(k))
g2(k)=A2(1+k2ω2(k))
・
・
・ (55)
gI(k)=AI(1+kIωI(k))
と表す。また任意の時刻kにおいて、−1<k1ω1(k), k2ω2(k), …, kIωI(k)<1が成り立つから、g1(k), g2(k), …, gI(k)の値の符号はそれぞれ一定である。聴覚的なバランスを考えれば、利得係数A1, A2, …, AIの値はすべて等しいか、あるいは少なくとも同符号の値であることが望ましく、この場合、g1(k), g2(k), …, gI(k)はすべて同符号の値を取る。また、スカラー値k1, k2, …, kIの値により、処理後の信号の歪み量を調節することができる。
【0066】
また、周波数シフトという観点から、式(55)とは別に、
g1(k)=A1cos(ω1(k)・k+φ1)
g2(k)=A2cos(ω2(k)・k+φ2)
・
・
・ (56)
gI(k)=AIcos(ωI(k)・k+φI)
として、各チャネルごとに異なる時変特性をもち、ある変動幅の中で変化する関数ω1(k), ω2(k), …, ωI(k)によって、受話信号x1(k), x2(k),…, xI(k)の周波数特性をシフトさせることにより、もとの多チャネル受話信号に対して、相互相関の変動付加を実現することもできる。
【0067】
▲3▼ 信号加算処理による実現方法
各チャネルの受話信号x1(k), x2(k), …, xI(k)を、異なる関数n1(k), n2(k), …, nI(k)と加算、すなわち、
x ̄1(k)=x1(k)+n1(k)
x ̄2(k)=x2(k)+n2(k)
・
・
・ (57)
x ̄I(k)=xI(k)+nI(k)
を満足するx ̄1(k), x ̄2(k), …, x ̄I(k)に変換することにより、各チャネルの受話信号間の相互相関への変動の付加を実現する。
【0068】
関数n1(k), n2(k), …, nI(k)は、聴覚的な影響を考慮すると、各チャネルの受話信号x1(k), x2(k),…, xI(k)と比べて、振幅値が適度に小さいことが望ましい。そこで、例えば、各チャネルの受話信号x1(k), x2(k),…, xI(k)の振幅の定格値をAr として、最大振幅がこのAr となるように規格化された異なる関数nr1(k), nr2(k), …, n2I(k)に対して、スカラー値λ1, λ2, …, λIを掛けることにより、n1(k), n2(k), …, nI(k)を、
n1(k)=λ1nr1(k)
n2(k)=λ2nr2(k)
・
・
・ (58)
nI(k)=λInrI(k)
と表す。スカラー値λ1,λ2, …, λIの値により、受話信号x1(k), x2(k),…,xI(k)に加える信号の大きさを定格値Arからどの程度小さくするのかを決定できる。
【0069】
▲4▼ ピッチシフト処理による実現方法
各チャネルの受話信号x1(k), x2(k), …, xI(k)を、各受話信号の周波数特性について、各チャネルごとに異なる時変の周波数伸縮処理、すなわち、ピッチシフト処理することにより得られるx ̄1(k), x ̄2(k), …, x ̄I(k)に変換することにより、各チャネルの受話信号間の相互相関への変動の付加を実現する。
【0070】
ピッチシフト処理は、周波数伸縮とともに、時間伸縮をともない、結局、時間領域で処理することが可能であるが、時間軸を縮小させれば、もとの信号が占めていた占有時間と等しくするための補間処理が必要となり、逆に時間軸を伸長させれば、もとの信号よりも占有時間が長くなるため間引き処理が必要となる。これらの補間・間引き処理は、無音区間を検出するなどして、その無音区間を引き延ばしたり、削除することにより実現できる。
【0071】
【発明の効果】
この発明方法は、2地点間において、多チャネルの伝送系と、多チャネルの集音・再生系を持つ端末装置とを用いることにより、双方の音響空間情報の伝達を可能とする臨場感の高い音声通信を実現する際に反響を消去するために適用することができる。
【0072】
簡単のため、図5に示したような、2地点間を2チャネルの伝送系で結び、2チャネルの集音・再生系を持つ端末装置を用いて、ステレオ音声通信を行なう場合を考える。例えば、A地点から複数の異なる話者が一人ずつそれぞれの席から発話した場合に、集音されるステレオ音声信号には、それぞれの発話者と2本のマイクロホンの位置関係に依存した相互相関がある。それらの音声信号がB地点側において再生されるとき、従来の反響消去装置を用いた場合では、反響路推定がステレオ信号の相互相関に依存して誤るため、話者の交替の度にA地点側に大きな反響が戻ってしまう。ところが、この発明の反響消去方法を用いた場合、話者の交替によるステレオ再生信号の相互相関の変動を抽出して、反響路推定の誤りを修正するため、話者の交替の度に起きていた反響の増大を低減できる。また、人間の発話音声を複数のマイクロホンによって集音するとき、固定した位置から発話した場合でも、集音された多チャネル音声信号の相互相関は微妙に変動している。この発明方法では、このような信号間の相互相関の微妙な変動をも有効に利用できる。
【0073】
図7の実施例において、周波数帯域分割部51、52内の帯域通過フィルタ200〜20N-1(分割部52の構成は51の構成と同様であり図示してない)による帯域制限の代わりに、前述のCrochiere and Rabinerの文献に示されているSSB法を使って帯域制限された実信号xn(k)、yn(k)を得て反響消去処理を行った場合の収束特性の計算機シミュレーション結果を図25に示す。計算機シミュレーションには実測したインパルス応答(512タップ、サンプリング周波数16kHz)を使用した。体や頭を固定して発話した音声をステレオ集音した信号を受話信号として用いた。反響信号にはS/N比=35dBとなるように近端雑音を加えた。
【0074】
図25は、反響路推定部で生成された擬似反響路結合ベクトルと真の反響路結合ベクトルとの誤差ベクトルの大きさを帯域分割数Nをパラメータとして比較したものである。ここで、この発明方法の相互相関変動成分抽出には、p=1次の射影アルゴリズム、すなわち、従来のNLMSアルゴリズムを用いている。帯域分割数N=1,4,8,32である。間引き率R=N/4である。各帯域のタップ数はN=1(フルバンド)のとき512+512とし、サブバンドでは間引き率に応じて少なく設定した。定常消去量が等しくなるように、第2のステップサイズμnおよび小さな正の定数δnを調節した。この結果より、この発明方法が信号間の相互相関の微妙な変動をも捉え、反響路推定の誤りを修正できること、さらに帯域分割数Nが大きくなるほど収束速度が速くなることが分かる。これは従来は予期し難い効果である。この発明ではチャネル間の微妙な変動を帯域分割手法を用いて有効にとらえることにより、収束速度を早めていると言える。
【0075】
図26に図25と同様なステレオ集音の場合の擬似反響路ベクトルの収束の様子を帯域分割数N=32の場合について射影アルゴリズムの次数をパラメータとして示す。帯域分割数N=32、間引き率R=8である。各帯域のタップ数は64+64とした場合である。射影アルゴリズムの次数は、pn=1,2,4,8,32である。ここでは、定常消去量が等しくなるように、第2のステップサイズμnおよび小さな正の定数δnを調節した。射影アルゴリズムの次数pnを増加させるに従って、収束速度が増加することがわかる。しかも従来における周波数分割処理法で得られる収束速度の高速化よりも、可成り高速に収束している。この図から、帯域分割数を行いさらに射影アルゴリズムを用いることが有効であることが理解される。このように、この発明を用いれば、小さな射影次数で速い収束速度が得られる。拡声通話系では人の移動などによる反響路の変動が多く、これに迅速に適応できることは大きな利点となる。
【0076】
更に、周波数帯域に応じて適切な射影次数を決めることにより、一層少ない演算量で高速の収束が可能となる。
ステレオ再生装置を有する端末において、図5に示したように相手側で2チャネル集音されたステレオ信号を受信してそのまま再生する方法の他に、多地点間通信等では、対地毎の受話信号を受信側で任意に音像定位処理することにより、受話環境を快適にする方法が考えられている。この発明方法は、このような多地点間通信用端末にも適用可能である。図27は、4地点間通信の構成を示したものである。各地点において、集音は1チャネル(モノラル)である。いま、D地点を例にとって説明する。D地点において、A、B、C地点からの受話信号について、それぞれ右、中央、左に定位するように音像定位処理を受話信号に処理部13dで施すことにより、新たな2チャネルステレオ受話信号を生成して2チャネルのステレオ再生を行なうものとする。この方法では、受話信号間の相互相関が一定であり、この場合は、受話信号間の相互相関の変動抽出・利用をしようにも、受話信号間の相互相関の変動そのものがないため、従来法と比較しても、反響路推定は停滞して改善されない。そこで、上記方法にこの発明の第3の方法を適用すれば、受話信号間の相互相関は常に変動することになり、先の場合と同様に、ステレオ受話信号の相互相関の変動を抽出して、反響路推定の誤りを修正するため、反響路推定が停滞することがなく、従来方法を用いた場合よりも反響の増大を低減できる。
【0077】
すなわち、一般に多チャネル再生系と1チャネル以上の集音系で構成される通信会議システムに、従来の反響消去方法を適用した場合、各チャネルの受話信号間に相互相関がある場合に、擬似反響路の推定が正しく行われないため、受話信号間の相互相関が変動する度に反響が増大する問題が生じていた。この発明の第1の方法においては、サブバンド法を用いて受話信号間の相互相関の変動成分を抽出・利用することにより、反響路推定の誤りを修正するため、上記のような問題を改善する効果がある。
【0078】
また、この発明の第2の方法においては、射影アルゴリズム等を用いて受話信号間の相互相関の変動成分を抽出・利用することにより、反響路推定の誤りを修正するため、上記のような問題を改善する効果があり、更に、周波数帯域に応じて適切な次数を決めることにより、一層少ない演算量で高速の収束が可能となる。
【0079】
また、この発明の第3の方法においては、受話信号間の相互相関に変動を付加する機能を有するため、反響路推定において誤った解に完全に停滞することがなく、真の反響路との誤差を小さくする方向へ反響路推定を続けるので、上記のような問題を改善する効果があり、更に、周波数帯域に応じて適切な相互相関変動付加機能を決めることにより、聴感上違和感を与えることなく高速の収束が可能となる。
【図面の簡単な説明】
【図1】従来の1チャネル反響消去装置の機能構成を示すブロック図。
【図2】Aはステップサイズ行列Aの対角成分αlを示す説明図、Bはステップサイズ行列Aの対角成分αlを階段状に近似する一例を示す説明図。
【図3】従来の多チャネル反響消去装置システムを示すブロック図。
【図4】従来の多チャネル反響消去装置の機能構成を示すブロック図。
【図5】ステレオ音声通信系を示すブロック図。
【図6】Aは受話信号間に相互相関がある場合の反響路推定動作を示す図。Bは前記相互相関が変動した場合の反響路推定動作を示す図である。
【図7】この発明の実施例の機能構成例を示すブロック図。
【図8】図7中の周波数帯域分割部51の内部の機能構成例を示すブロック図。
【図9】この発明による第1の方法を適用した多チャネル反響消去装置の構成を示すブロック図。
【図10】図7中の第n番目の周波数帯域の推定部19nの内部の機能構成例を示すブロック図。
【図11】音声信号の周波数スペクトルの一例を示す説明図。
【図12】図7中の第1又は第2の次数決定方法を適用する場合の次数決定制御部54の機能構成例を示すブロック図。
【図13】射影次数を増加した場合の反響消去量の収束速度の変化を模式的に示す図。
【図14】第1の次数決定方法による次数決定手順を示すフロー図。
【図15】第2の次数決定方法による次数決定手順を示すフロー図。
【図16】第3の次数決定方法を実施する場合の次数決定制御部54の構成例を示すブロック図。
【図17】第3の次数決定方法による次数決定手順を示すフロー図。
【図18】この発明のもう1つの実施例の機能構成例を示すブロック図。
【図19】図18中の周波数帯域分割部61の内部の機能構成例を示すブロック図。
【図20】周波数帯域信号の一例を示す説明図。
【図21】図18中の第n番目の周波数帯域の推定部67nの内部の機能構成例を示すブロック図。
【図22】図14中の第1又は第2の次数決定方法を適用する場合の次数決定制御部64の構成例を示すブロック図。
【図23】図14中の第3の次数決定方法を適用する場合の次数決定制御部64の構成例を示すブロック図。
【図24】この発明による第3の方法を適用した多チャネル反響消去装置の構成を示すブロック図。
【図25】この発明方法による第1の方法による擬似反響路ベクトルの収束の様子を示す図。
【図26】この発明方法による第2の方法による擬似反響路ベクトルの収束の様子を示す図。
【図27】4地点間通信において音像定位機能をもたせる音声通信会議方式を示すブロック図。
Claims (15)
- 各チャネルの受話信号をチャネルごとの再生器で音響信号に再生し、
これら各音響信号が、上記各再生器から集音器に到る各反響路を経由して集音され、
その集音器からの反響信号から擬似反響信号を差し引き、
その残りである残留反響信号と、上記各チャネルの受話信号とから修正ベクトルを求め、
この修正ベクトルを用いて上記各反響路のインパルス応答の推定値を逐次修正し、
その修正されたインパルス応答をもつ各擬似反響路を生成し、
これら各擬似反響路に上記受話信号の対応するものをそれぞれ印加して上記擬似反響信号を生成する多チャネル音声通信会議用反響消去方法において、
(a) 上記受話信号と上記反響信号とをそれぞれN個の周波数帯域に分割してN個のサブバンド受話信号とN個のサブバンド反響信号を生成するステップと、Nは2以上の整数であり、
(b) 上記反響路の上記N個の周波数帯域におけるそれぞれのインパルスレスポンスを模擬するそれぞれ予め決めたタップ数のフィルタ係数が与えられるディジタルフィルタによりそれぞれ構成されたN個の擬似反響路に上記N個のサブバンド受話信号をそれぞれ通してN個の擬似反響信号を生成するステップと、
(c) 上記N個のサブバンド反響信号から対応する上記N個の擬似反響信号をそれぞれ減算することによりN個の周波数帯域の反響消去誤差信号を生成するステップと、
(d) それぞれの上記N個の反響消去誤差信号と、上記N個のサブバンド受話信号の対応するものとから、対応する上記反響消去誤差信号を最小とするように上記ディジタルフィルタに与える上記フィルタ係数をそれぞれ逐次的に修正するステップと、
(e) 上記N個の周波数帯域の上記反響消去誤差信号を合成して上記反響信号が抑圧された全周波数帯域の送信信号を生成するステップ、
とを含み、
上記各チャネルの受話信号の系列よりなる受話信号ベクトルを結合させた受話信号結合ベクトルを求め、この受話信号結合ベクトルの過去のものに対する現在のそれとの相関の変動を検出して上記各チャネルの受話信号の相互相関の変動成分を抽出し、この変動成分を上記修正ベクトルとする。 - 請求項1のサブバンド多チャネル音声通信会議用反響消去方法において、
上記各チャネルの上記受話信号結合ベクトルの過去のものに対する現在のそれとの相関の変動を検出する検出方法を上記N個の周波数帯域で最適に設定する。 - 請求項2のサブバンド多チャネル音声通信会議用反響消去方法において、
上記各チャネルの上記受話信号結合ベクトルの過去のものに対する現在のそれとの相関の変動を検出する検出方法が射影アルゴリズムまたはES射影アルゴリズムであり、射影次数を上記N個の周波数帯域で最適に設定する。 - 請求項3のサブバンド多チャネル音声通信会議用反響消去方法において、
上記N個の周波数帯域の上記射影アルゴリズム又はES射影アルゴリズムのそれぞれの次数は上記受話信号に対しそれぞれの周波数帯域において反響消去量の収束速度がほぼ飽和する最小の値にそれぞれ決められている。 - 請求項3のサブバンド多チャネル音声通信会議用反響消去方法において、
上記N個の周波数帯域の上記射影アルゴリズム又はES射影アルゴリズムのそれぞれの次数は、上記受話信号を線形予測フィルタにより白色化したときの推定誤差ベクトルの白色化がほぼ飽和する最小の値にそれぞれ決められている。 - 請求項3のサブバンド多チャネル音声通信会議用反響消去方法において,
上記N個の周波数帯域のそれぞれの擬似反響路を構成する上記ディジタルフィルタのタップ数は、所望の受話信号の周波数領域におけるエネルギー分布、室内残響特性、及び人間の聴覚特性の少なくとも1つに基づいて予め決められている。 - 請求項3のサブバンド多チャネル音声通信会議用反響消去方法において、
上記N個の周波数帯域の低い周波数帯域に対応する上記ディジタルフィルタのタップ数は、高い周波数帯域の対応するディジタルフィルタのタップ数より大とされている。 - 請求項3のサブバンド多チャネル音声通信会議用反響消去方法において、
上記低い周波数帯域における上記射影アルゴリズム又はES射影アルゴリズムの次数は、上記高い周波数帯域における上記射影アルゴリズム又はES射影アルゴリズムの次数より大とされている。 - 各チャネルの受話信号をチャネルごとの再生器で音響信号に再生し、
これら各音響信号が、上記各再生器から集音器に到る各反響路を経由して集音され、
その集音器からの反響信号から擬似反響信号を差し引き、
その残りである残留反響信号と、上記各チャネルの受話信号とから修正ベクトルを求め、
この修正ベクトルを用いて上記各反響路のインパルス応答の推定値を逐次修正し、
その修正されたインパルス応答をもつ各擬似反響路を生成し、
これら各擬似反響路に上記受話信号の対応するものをそれぞれ印加して上記擬似反響信号を生成する多チャネル音声通信会議用反響消去方法において、
(a) 上記受話信号と上記反響信号とをそれぞれN個の周波数帯域に分割してN個のサブバンド受話信号とN個のサブバンド反響信号を生成するステップと、Nは2以上の整数であり、
(b) 上記反響路の上記N個の周波数帯域におけるそれぞれのインパルスレスポンスを模擬するそれぞれ予め決めたタップ数のフィルタ係数が与えられるディジタルフィルタによりそれぞれ構成されたN個の擬似反響路に上記N個のサブバンド受話信号をそれぞれ通してN個の擬似反響信号を生成するステップと、
(c) 上記N個のサブバンド反響信号から対応する上記N個の擬似反響信号をそれぞれ減算することによりN個の周波数帯域の反響消去誤差信号を生成するステップと、
(d) それぞれの上記N個の反響消去誤差信号と、上記N個のサブバンド受話信号の対応するものとから、対応する上記反響消去誤差信号を最小とするように上記ディジタルフィルタに与える上記フィルタ係数をそれぞれ逐次的に修正するステップと、
(e) 上記N個の周波数帯域の上記反響消去誤差信号を合成して上記反響信号が抑圧された全周波数帯域の送信信号を生成するステップ、
とを含み、
上記各チャネルの受話信号間の相互相関に変動成分を付加して、上記再生器で音響信号を再生し、相互相関の変動成分が付加された受話信号を上記修正ベクトルの導出に利用する。 - 請求項9のサブバンド多チャネル音声通信会議用反響消去方法において、
再生チャネル数をIとして、上記各チャネルの受話信号それぞれを、離散時間kの関数として、x1(k), x2(k), …, xI(k)と表したとき、これらの受話信号を各チャネルに対して異なる時変特性をもつ時変フィルタに印加、これら時変フィルタのインパルス応答f1(k), f2(k), …, fI(k), 畳み込み演算を*で表すとして、
x ̄1(k)=f1(k)*x1(k)
x ̄2(k)=f2(k)*x2(k)
・
・
・
x ̄I(k)=fI(k)*xI(k)
を満足するx ̄1(k), x ̄2(k), …, x ̄I(k)に変換することにより、各チャネルの受話信号間の相互相関への変動の付加を実現する。 - 請求項9のサブバンド多チャネル音声通信会議用反響消去方法において,
再生チャネル数をIとして、上記各チャネルの受話信号それぞれを、離散時間kの関数として、x1(k), x2(k), …, xI(k)と表したとき、これらの受話信号を各チャネルに対して異なる関数g1(k), g2(k), …, gI(k)と乗算し、
x ̄1(k)=g1(k)・x1(k)
x ̄2(k)=g2(k)・x2(k)
・
・
・
x ̄I(k)=gI(k)・xI(k)
を満足するx ̄1(k), x ̄2(k), …, x ̄I(k)に変換することにより、上述の各チャネルの受話信号間の相互相関への変動の付加を実現する。 - 請求項9のサブバンド多チャネル音声通信会議用反響消去方法において,
再生チャネル数をIとして、上記各チャネルの受話信号それぞれを、離散時間kの関数として、x1(k), x2(k), …, xI(k)と表したとき、これらの受話信号を各チャネルに対して異なる関数n1(k), n2(k),…, nI(k)と加算し、
x ̄1(k)=x1(k)+n1(k)
x ̄2(k)=x2(k)+n2(k)
・
・
・
x ̄I(k)=xI(k)+nI(k)
を満足するx ̄1(k), x ̄2(k), …, x ̄I(k)に変換することにより、上述の各チャネルの受話信号間の相互相関への変動の付加を実現する。 - 請求項9のサブバンド多チャネル音声通信会議用反響消去方法において,
再生チャネル数をIとして、上記各チャネルの受話信号それぞれを、離散時間kの関数として、x1(k), x2(k), …, xI(k)と表したとき、これらの受話信号を各受話信号の周波数特性について、各チャネルごとに異なる時変の周波数伸縮処理することにより得られるx ̄1(k), x ̄2(k),…, x ̄I(k)に変換することにより、各チャネルの受話信号間の相互相関への変動の付加を実現する。 - 請求項9〜13のサブバンド多チャネル音声通信会議用反響消去方法において、
各チャネルの受話信号間の相互相関への変動の付加方法を、聴覚上の品質劣化が少なくなるように上記N個の周波数帯域で最適に設定する。 - 請求項1〜14のサブバンド多チャネル音声通信会議用反響消去方法において、
上記N個のサブバンド受話信号と上記N個のサブバンド反響信号を生成するステップは、上記受話信号及び上記反響信号をそれぞれ上記N個の周波数帯域に分割してからそれぞれ予め決めた間引き率で信号を間引くステップを含む。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP24352496A JP3654470B2 (ja) | 1996-09-13 | 1996-09-13 | サブバンド多チャネル音声通信会議用反響消去方法 |
US08/927,961 US6246760B1 (en) | 1996-09-13 | 1997-09-11 | Subband echo cancellation method for multichannel audio teleconference and echo canceller using the same |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP24352496A JP3654470B2 (ja) | 1996-09-13 | 1996-09-13 | サブバンド多チャネル音声通信会議用反響消去方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
JPH1093680A JPH1093680A (ja) | 1998-04-10 |
JP3654470B2 true JP3654470B2 (ja) | 2005-06-02 |
Family
ID=17105193
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP24352496A Expired - Lifetime JP3654470B2 (ja) | 1996-09-13 | 1996-09-13 | サブバンド多チャネル音声通信会議用反響消去方法 |
Country Status (2)
Country | Link |
---|---|
US (1) | US6246760B1 (ja) |
JP (1) | JP3654470B2 (ja) |
Families Citing this family (67)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3710606B2 (ja) * | 1997-09-04 | 2005-10-26 | 松下電器産業株式会社 | エコーキャンセラ |
DE19806015C2 (de) * | 1998-02-13 | 1999-12-23 | Siemens Ag | Verfahren zur Verbesserung der akustischen Rückhördämpfung in Freisprecheinrichtungen |
DE19814971A1 (de) * | 1998-04-03 | 1999-10-07 | Daimlerchrysler Aerospace Ag | Verfahren zur Störbefreiung eines Mikrophonsignals |
AU3430199A (en) * | 1998-04-08 | 1999-11-01 | British Telecommunications Public Limited Company | Echo cancellation |
JP3781902B2 (ja) * | 1998-07-01 | 2006-06-07 | 株式会社リコー | 音像定位制御装置および音像定位制御方式 |
US7072465B1 (en) * | 1999-01-06 | 2006-07-04 | Agere Systems, Inc. | Adaptive filter for network echo cancellation |
US6744886B1 (en) * | 1999-01-06 | 2004-06-01 | Lucent Technologies Inc. | Adaptive filter utilizing proportional affine projection algorithm |
JP2003500937A (ja) * | 1999-05-25 | 2003-01-07 | ブリティッシュ・テレコミュニケーションズ・パブリック・リミテッド・カンパニー | 音響エコー消去 |
WO2001001571A1 (en) * | 1999-06-24 | 2001-01-04 | Koninklijke Philips Electronics N.V. | Acoustic echo and noise cancellation |
US6694020B1 (en) * | 1999-09-14 | 2004-02-17 | Agere Systems, Inc. | Frequency domain stereophonic acoustic echo canceller utilizing non-linear transformations |
FI114258B (fi) * | 2000-06-09 | 2004-09-15 | Nokia Corp | Menetelmä häiriön vaikutuksen vähentämiseksi ja vastaanotin |
US6700978B2 (en) | 2000-09-08 | 2004-03-02 | Intel Corporation | Method and apparatus for fast converging affine projection based echo canceller |
US6909782B2 (en) * | 2000-09-08 | 2005-06-21 | Intel Corporation | Fast converging affine projection based echo canceller for sparse multi-path channels |
KR100355033B1 (ko) * | 2000-12-30 | 2002-10-19 | 주식회사 실트로닉 테크놀로지 | 선형예측 분석을 이용한 워터마크 삽입/추출 장치 및 그방법 |
US6766021B2 (en) * | 2001-03-13 | 2004-07-20 | Adaptive Digital Technologies | Echo canceller |
US20020191779A1 (en) * | 2001-03-23 | 2002-12-19 | Tien Pham | System for convolutional echo cancellation by iterative autocorrelation |
JP3608525B2 (ja) | 2001-05-09 | 2005-01-12 | ヤマハ株式会社 | 2チャンネルエコーキャンセル用フィルタのインパルス応答設定方法および2チャンネルエコーキャンセラ並びに双方向2チャンネル音声伝送装置 |
GB2379369B (en) * | 2001-08-29 | 2005-06-29 | Zarlink Semiconductor Inc | Subband echo location and double-talk detection in communication systems |
FR2841721B1 (fr) * | 2002-06-28 | 2004-08-20 | France Telecom | Dispositif de traitement d'echo pour systeme de communication de type monovoie ou multivoies |
CA2399159A1 (en) * | 2002-08-16 | 2004-02-16 | Dspfactory Ltd. | Convergence improvement for oversampled subband adaptive filters |
NO318401B1 (no) * | 2003-03-10 | 2005-03-14 | Tandberg Telecom As | Et audio-ekkokanselleringssystem og en fremgangsmate for a tilveiebringe et ekkodempet utgangssignal fra et ekkotillagt signal |
US20040252652A1 (en) * | 2003-06-10 | 2004-12-16 | Alexander Berestesky | Cross correlation, bulk delay estimation, and echo cancellation |
DE10351793B4 (de) * | 2003-11-06 | 2006-01-12 | Herbert Buchner | Adaptive Filtervorrichtung und Verfahren zum Verarbeiten eines akustischen Eingangssignals |
JP4663630B2 (ja) * | 2004-03-22 | 2011-04-06 | ティーオーエー株式会社 | マルチチャンネルシステム同定装置 |
JP4630956B2 (ja) * | 2004-03-30 | 2011-02-09 | 学校法人早稲田大学 | ハウリング周波数成分強調方法およびその装置、ハウリング検出方法およびその装置、ハウリング抑圧方法およびその装置、ピーク周波数成分強調方法およびその装置 |
US7352858B2 (en) * | 2004-06-30 | 2008-04-01 | Microsoft Corporation | Multi-channel echo cancellation with round robin regularization |
US7903137B2 (en) * | 2004-10-15 | 2011-03-08 | Lifesize Communications, Inc. | Videoconferencing echo cancellers |
US8457614B2 (en) | 2005-04-07 | 2013-06-04 | Clearone Communications, Inc. | Wireless multi-unit conference phone |
JP4581114B2 (ja) * | 2005-05-16 | 2010-11-17 | 株式会社国際電気通信基礎技術研究所 | 適応型ビーム形成器 |
US20070036342A1 (en) * | 2005-08-05 | 2007-02-15 | Boillot Marc A | Method and system for operation of a voice activity detector |
US7426270B2 (en) * | 2005-08-10 | 2008-09-16 | Clarity Technologies, Inc. | Method and system for clear signal capture |
KR101277041B1 (ko) | 2005-09-01 | 2013-06-24 | 파나소닉 주식회사 | 멀티 채널 음향 신호 처리 장치 및 방법 |
US8189765B2 (en) * | 2006-07-06 | 2012-05-29 | Panasonic Corporation | Multichannel echo canceller |
US8199927B1 (en) | 2007-10-31 | 2012-06-12 | ClearOnce Communications, Inc. | Conferencing system implementing echo cancellation and push-to-talk microphone detection using two-stage frequency filter |
US8050398B1 (en) | 2007-10-31 | 2011-11-01 | Clearone Communications, Inc. | Adaptive conferencing pod sidetone compensator connecting to a telephonic device having intermittent sidetone |
CN101878637B (zh) * | 2007-11-29 | 2013-09-18 | 艾利森电话股份有限公司 | 用于对语音信号进行回声消除的方法和配置 |
US8170226B2 (en) * | 2008-06-20 | 2012-05-01 | Microsoft Corporation | Acoustic echo cancellation and adaptive filters |
JP5284475B2 (ja) * | 2008-08-25 | 2013-09-11 | ドルビー ラボラトリーズ ライセンシング コーポレイション | 前白色化を伴うlmsアルゴリズムによって適応させられる適応フィルタの更新済みフィルタ係数を決定する方法 |
US8605890B2 (en) * | 2008-09-22 | 2013-12-10 | Microsoft Corporation | Multichannel acoustic echo cancellation |
JP5113794B2 (ja) * | 2009-04-02 | 2013-01-09 | 日本電信電話株式会社 | 適応マイクロホンアレイ残響抑圧装置、適応マイクロホンアレイ残響抑圧方法及びプログラム |
JP5161157B2 (ja) * | 2009-06-15 | 2013-03-13 | 日本電信電話株式会社 | 周波数領域エコー除去装置、周波数領域エコー除去方法、プログラム |
US8724798B2 (en) | 2009-11-20 | 2014-05-13 | Adobe Systems Incorporated | System and method for acoustic echo cancellation using spectral decomposition |
GB2493801B (en) | 2011-08-18 | 2014-05-14 | Ibm | Improved audio quality in teleconferencing |
JP5698110B2 (ja) * | 2011-11-30 | 2015-04-08 | 日本電信電話株式会社 | マルチチャネルエコー消去方法、マルチチャネルエコー消去装置、およびプログラム |
CN102457632B (zh) * | 2011-12-29 | 2014-07-30 | 歌尔声学股份有限公司 | 一种多受话端回声消除方法 |
US9173025B2 (en) | 2012-02-08 | 2015-10-27 | Dolby Laboratories Licensing Corporation | Combined suppression of noise, echo, and out-of-location signals |
US8712076B2 (en) | 2012-02-08 | 2014-04-29 | Dolby Laboratories Licensing Corporation | Post-processing including median filtering of noise suppression gains |
US20140016794A1 (en) * | 2012-07-13 | 2014-01-16 | Conexant Systems, Inc. | Echo cancellation system and method with multiple microphones and multiple speakers |
CN102876969B (zh) * | 2012-07-31 | 2015-03-04 | 宝山钢铁股份有限公司 | 一种超高强度高韧性耐磨钢板及其制造方法 |
CN103888630A (zh) | 2012-12-20 | 2014-06-25 | 杜比实验室特许公司 | 用于控制声学回声消除的方法和音频处理装置 |
US20150371655A1 (en) * | 2014-06-19 | 2015-12-24 | Yang Gao | Acoustic Echo Preprocessing for Speech Enhancement |
EP3275208B1 (en) | 2015-03-25 | 2019-12-25 | Dolby Laboratories Licensing Corporation | Sub-band mixing of multiple microphones |
US9712866B2 (en) | 2015-04-16 | 2017-07-18 | Comigo Ltd. | Cancelling TV audio disturbance by set-top boxes in conferences |
US9947337B1 (en) * | 2017-03-21 | 2018-04-17 | Omnivision Technologies, Inc. | Echo cancellation system and method with reduced residual echo |
US10013995B1 (en) * | 2017-05-10 | 2018-07-03 | Cirrus Logic, Inc. | Combined reference signal for acoustic echo cancellation |
GB201709849D0 (en) * | 2017-06-20 | 2017-08-02 | Nokia Technologies Oy | Processing audio signals |
US10200540B1 (en) * | 2017-08-03 | 2019-02-05 | Bose Corporation | Efficient reutilization of acoustic echo canceler channels |
CN111149370B (zh) | 2017-09-29 | 2021-10-01 | 杜比实验室特许公司 | 会议***中的啸叫检测 |
CN111213359B (zh) | 2017-10-04 | 2021-09-10 | 主动音频有限公司 | 回声消除器和用于回声消除器的方法 |
US10566008B2 (en) * | 2018-03-02 | 2020-02-18 | Cirrus Logic, Inc. | Method and apparatus for acoustic echo suppression |
JP6894402B2 (ja) * | 2018-05-23 | 2021-06-30 | 国立大学法人岩手大学 | システム同定装置及び方法及びプログラム及び記憶媒体 |
US11417351B2 (en) | 2018-06-26 | 2022-08-16 | Google Llc | Multi-channel echo cancellation with scenario memory |
CN110782911A (zh) * | 2018-07-30 | 2020-02-11 | 阿里巴巴集团控股有限公司 | 音频信号处理方法、装置、设备和存储介质 |
US11600287B2 (en) * | 2019-03-28 | 2023-03-07 | Bose Corporation | Cancellation of vehicle active sound management signals for handsfree systems |
CN111294474B (zh) * | 2020-02-13 | 2021-04-16 | 杭州国芯科技股份有限公司 | 一种双端通话检测方法 |
EP4205108A1 (en) * | 2020-09-02 | 2023-07-05 | Huawei Technologies Co., Ltd. | Acoustic processing device for multichannel nonlinear acoustic echo cancellation |
WO2022048736A1 (en) * | 2020-09-02 | 2022-03-10 | Huawei Technologies Co., Ltd. | Acoustic processing device for mimo acoustic echo cancellation |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH04123621A (ja) * | 1990-09-14 | 1992-04-23 | Nippon Telegr & Teleph Corp <Ntt> | 反響消去装置 |
JP2508574B2 (ja) * | 1992-11-10 | 1996-06-19 | 日本電気株式会社 | 多チャンネルエコ―除去装置 |
JP3397269B2 (ja) * | 1994-10-26 | 2003-04-14 | 日本電信電話株式会社 | 多チャネル反響消去方法 |
FR2729024A1 (fr) * | 1994-12-30 | 1996-07-05 | Matra Communication | Annuleur d'echo acoustique avec filtrage en sous-bandes |
CA2174366C (en) * | 1995-04-20 | 2001-03-27 | Shoji Makino | Subband echo cancellation method using projection algorithm |
DE69634027T2 (de) * | 1995-08-14 | 2005-12-22 | Nippon Telegraph And Telephone Corp. | Akustischer Teilband-Echokompensator |
DE69628618T2 (de) * | 1995-09-26 | 2004-05-13 | Nippon Telegraph And Telephone Corp. | Verfahren und Einrichtung zur mehrkanaligen Kompensation eines akustischen Echos |
US5721782A (en) * | 1996-03-25 | 1998-02-24 | Motorola, Inc. | Partitioned echo canceler utilizing decimation echo location |
-
1996
- 1996-09-13 JP JP24352496A patent/JP3654470B2/ja not_active Expired - Lifetime
-
1997
- 1997-09-11 US US08/927,961 patent/US6246760B1/en not_active Expired - Lifetime
Also Published As
Publication number | Publication date |
---|---|
US6246760B1 (en) | 2001-06-12 |
JPH1093680A (ja) | 1998-04-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP3654470B2 (ja) | サブバンド多チャネル音声通信会議用反響消去方法 | |
JP3397269B2 (ja) | 多チャネル反響消去方法 | |
CN111292759B (zh) | 一种基于神经网络的立体声回声消除方法及*** | |
Shimauchi et al. | Stereo projection echo canceller with true echo path estimation | |
Gilloire et al. | Using auditory properties to improve the behaviour of stereophonic acoustic echo cancellers | |
RU2495506C2 (ru) | Устройство и способ расчета параметров управления фильтра эхоподавления и устройство и способ расчета величины задержки | |
EP0739102B1 (en) | Subband echo cancellation method using projection algorithm | |
KR100825214B1 (ko) | 비대칭 다중 채널 필터 | |
JP3099870B2 (ja) | 音響エコーキャンセラ | |
US6744887B1 (en) | Acoustic echo processing system | |
JPH04123621A (ja) | 反響消去装置 | |
Schobben | Real-time adaptive concepts in acoustics: Blind signal separation and multichannel echo cancellation | |
CN111213359A9 (zh) | 回声消除器和用于回声消除器的方法 | |
Benesty et al. | A hybrid mono/stereo acoustic echo canceler | |
JP3756828B2 (ja) | 反響消去方法、この方法を実施する装置、プログラムおよびその記録媒体 | |
JP3673727B2 (ja) | 反響消去方法、その装置、そのプログラム及びその記録媒体 | |
JP3402427B2 (ja) | 多チャネル反響消去方法及び装置 | |
JP3616341B2 (ja) | 多チャネルエコーキャンセル方法、その装置、そのプログラム及び記録媒体 | |
Surin et al. | An adaptive noise decorrelation technique for stereophonic acoustic echo cancellation | |
JP3707572B2 (ja) | サブバンド反響消去方法 | |
JP3339612B2 (ja) | 多チャンネル音声通信会議用反響消去方法 | |
Romoli et al. | An interactive optimization procedure for stereophonic acoustic echo cancellation systems | |
Chang et al. | Distributed Kalman Filtering for Speech Dereverberation and Noise Reduction in Acoustic Sensor Networks | |
Bendoumia | New two-microphone simplified sub-band forward algorithm based on separated variable step-sizes for acoustic noise reduction | |
Nezamdoust et al. | Frequency-Domain Functional Links For Nonlinear Feedback Cancellation In Hearing Aids |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20050201 |
|
RD01 | Notification of change of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7426 Effective date: 20050224 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20050224 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090311 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090311 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100311 Year of fee payment: 5 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110311 Year of fee payment: 6 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110311 Year of fee payment: 6 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120311 Year of fee payment: 7 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130311 Year of fee payment: 8 |
|
S531 | Written request for registration of change of domicile |
Free format text: JAPANESE INTERMEDIATE CODE: R313531 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
EXPY | Cancellation because of completion of term |