JP2001251224A - 反響消去方法および反響消去装置 - Google Patents

反響消去方法および反響消去装置

Info

Publication number
JP2001251224A
JP2001251224A JP2000061265A JP2000061265A JP2001251224A JP 2001251224 A JP2001251224 A JP 2001251224A JP 2000061265 A JP2000061265 A JP 2000061265A JP 2000061265 A JP2000061265 A JP 2000061265A JP 2001251224 A JP2001251224 A JP 2001251224A
Authority
JP
Japan
Prior art keywords
signal
echo
vector
echo path
band
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.)
Pending
Application number
JP2000061265A
Other languages
English (en)
Inventor
Akira Emura
暁 江村
Suehiro Shimauchi
末廣 島内
Shigeaki Aoki
茂明 青木
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nippon Telegraph and Telephone Corp
Original Assignee
Nippon Telegraph and Telephone Corp
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Nippon Telegraph and Telephone Corp filed Critical Nippon Telegraph and Telephone Corp
Priority to JP2000061265A priority Critical patent/JP2001251224A/ja
Publication of JP2001251224A publication Critical patent/JP2001251224A/ja
Pending legal-status Critical Current

Links

Landscapes

  • Telephone Function (AREA)
  • Cable Transmission Systems, Equalization Of Radio And Reduction Of Echo (AREA)

Abstract

(57)【要約】 【課題】 反響路インパルス応答の周波数に依存する性
質に着目して適応フィルタによる真の反響路インパルス
応答への収束特性を向上させる。 【解決手段】 受話信号x(k)は再生器2により音響信号
に再生され、収音器3により収音された信号y(k)は減算
器7により擬似反響路6からの擬似反響信号u(k)が
減算されて残留信号e(k)が得られる。受話信号ベクトル
生成部5aにて受話信号x(k)から受話信号ベクトルX(k)が
生成され、修正方向算出部5bにて正の対称行列Bを受話
信号ベクトルX(k)に乗算させて修正方向BX(k)が算出さ
れる。修正情報生成部5cで受話信号ベクトルX(k),修正
方向BX(k),誤差e(k)から反響路減衰特性の周波数によ
る相違を反映する修正ベクトルを算出し、その修正ベク
トルとタップ係数記憶部5dの前回の擬似反響路推定値W
(k)とを加算器5eにて加算して次時刻の擬似反響路推
定値W(k+1) を得る。

Description

【発明の詳細な説明】
【0001】
【発明の属する技術分野】本発明は、テレビ会議等の拡
声通話系において受話音が送話音に混合することによっ
て生ずる反響、会話を妨げハウリングの原因となる反響
を消去する反響消去方法および反響消去装置に関する。
【0002】
【従来の技術】図10は従来の反響消去方法の一例を示
すブロック図である。受話信号x(k)を受ける受話入力端
1からスピーカ2に至る受話系と、マイクロホン3から
送話出力端4に至る送話系とからなる拡声通話系におい
て、受話信号x(k)が擬似反響路6へ供給され、その擬似
反響路6からの擬似反響信号u(k)を、サンプル値化され
た反響信号y(k)から減算器7で差し引くことにより反響
信号y(k)は消去される。また、残差
【0003】
【数1】e(k)=y(k)−u(k) が0に近づくように推定回路5によって、擬似反響路6の
修正が逐次行われることによって、擬似反響路6は反響
路の経時変動に追従する。
【0004】学習同定法では、この擬似反響路6のイン
パルス応答W(k)は残差e(k)と受話信号からなるベクトル
X(k)をもちいて、次の(式1)のように逐次修正され
る。
【0005】
【数2】
【0006】これにより、擬似反響路のインパルス応答
W(k)は、真の反響路インパルス応答に近づいていくが、
1回の修正では、擬似反響路のインパルス応答W(k)は、
受話信号ベクトルX(k)の方向にしか修正されない。
【0007】真の反響路インパルス応答では、スピーカ
から直接マイクロホンに収音される直接波のパワーは大
きく、壁等を反射してマイクロホンに収音される間接波
のパワーは反射を繰り返すほどパワーが減衰している。
実際の会議室にスピーカとマイクロホンを設置して反
響路のインパルス応答とその変動を測定すると、インパ
ルス応答の振幅包絡は一般に指数減衰することが知られ
ている。
【0008】擬似反響路インパルス応答の真の反響路イ
ンパルス応答への収束特性を改善するために、反響路の
インパルス応答およびその変動特性が指数減衰特性を持
つことに着目して擬似反響路インパルス応答のタップ、
すなわち各係数毎にステップゲインを設定するESアルゴ
リズムが提案されている(S. Makino, Y. Kaneda, N.Ko
izumi:“ Exponentially Weighted Step size NLMS Ada
ptive Filter Basedon the Statistics of a Room Impu
lse Response”, IEEE Transactions on SPEECG AND AU
DIO PROCESSING, Vol. 1, No. 1, pp. 101-108, Januar
y 1993)。
【0009】この反響消去方法は、インパルス応答W(k)
の中で直接波に対応する部分を大きく、間接波に対応す
る部分をパワーの減衰に応じて小さく修正するように、
受信信号ベクトルに定数行列Aを乗算して修正方向AX(k)
を求め、インパルス応答W(k)を次の(式2)のように逐
次修正する。
【0010】
【数3】
【0011】演算量と記憶量をほとんど増加させること
なく、この反響消去方法をデジタル・シグナルプロセッ
サに実装するため、反響路インパルス応答の指数減衰特
性を階段状に簡略化することも提案されている。
【0012】いずれの場合にも、従来の反響消去方法で
は、反響路インパルス応答の減衰特性は周波数によら
ず、ほぼ均一と仮定している。
【0013】
【発明が解決しようとする課題】しかしながら、実際の
音声通信会議の場では、反響路インパルス応答の減衰特
性が周波数によらず均一とみなせるとは限らない。例え
ば、オーディオ帯域のような広帯域で反響消去を行うと
き、10kHz以上では壁面による吸音効果が大きくなり、
反響路減衰特性は周波数によらず均一とはみなせなくな
る。
【0014】また、電話帯域でも、スピーカ―マイクロ
ホン間の音響結合量を低減するために、音声通信会議に
は指向性のあるマイクロホンを用いる場合が多い。この
ような音響系の反響路インパルス応答を時間―周波数領
域で分析すると、反響路の減衰特性は周波数により均一
でない。
【0015】本発明は、このような反響路インパルス応
答の周波数に依存する性質に着目することで、適応フィ
ルタによる真の反響路インパルス応答への収束特性を向
上させた反響消去方法および反響消去装置を提供するこ
とを目的とする。
【0016】
【課題を解決するための手段】上記目的を達成するた
め、請求項1の反響消去方法の発明は、受話信号を再生
器で音響信号に再生するステップと、再生された前記音
響信号を反響路を経由して収音器で収音するステップ
と、収音された収音信号から擬似反響信号を減算するス
テップと、反響路減衰特性の周波数による相違を反映す
る行列を前記受話信号ベクトルに乗算して修正方向を算
出するステップと、前記受話信号ベクトル,前記修正方
向,前記減算の残りである残留信号とから修正ベクトル
を算出するステップと、前記修正ベクトルを用いて前記
反響路のインパルス応答の推定を逐次修正するステップ
と、前記修正されたインパルス応答をもつ擬似反響路を
生成するステップと、前記擬似反響路に前記受話信号を
印加することで前記擬似反響信号を生成するステップと
を有することを特徴とする。
【0017】ここで、周波数毎に異なる反響路減衰特性
を反映する前記行列がKブロック(Kは2以上の正の整
数)に分割されており、前記受話信号ベクトルから前記
修正方向の各係数を求めるFIRフィルタ係数が前記各
ブロック毎に設定されていることを特徴とすることがで
きる。
【0018】また、前記受話信号ベクトルから前記修正
方向を算出するための前記行列の乗算が、シフトバッフ
ァと、該シフトバッファ間のFIRフィルタ、もしくは
前記受話信号ベクトルの一部から前記シフトバッファへ
のFIRフィルタにより遂行されることを特徴とするこ
とができる。
【0019】また、前記擬似反響路に印加される前記受
話信号を帯域通過フィルタ通過後にあらかじめ決められ
た間引き率で間引かれたサブバンド受話信号とし、 前
記収音器からの収音信号を帯域通過フィルタ通過後にあ
らかじめ決められた間引き率で間引かれたサブバンド収
音信号とし、前記修正ベクトルが該サブバンド受話信号
と該サブバンド収音信号から生成され、反響消去後の信
号をあらかじめ決めれらた補間率により補間して帯域通
過フィルタを通過させることを特徴とすることができ
る。
【0020】また、M個(Mは2以上の正の整数)の収
音器により収音された信号を対応するそれぞれM個のフ
ィルタにより処理するステップと、前記M個のフィルタ
により処理された信号を加算するステップとを有し、該
加算して得られた信号を前記反響路の出力である収音信
号として取り扱うことを特徴とすることができる。
【0021】上記目的を達成するため、請求項6の反響
消去装置の発明は、受話信号を音響信号に再生する再生
器と、再生された前記音響信号を反響路を経由して収音
する収音器と、収音された収音信号から擬似反響信号を
減算する手段と、反響路減衰特性の周波数による相違を
反映する行列を前記受話信号ベクトルに乗算して修正方
向を算出する手段と、前記受話信号ベクトル,前記修正
方向,前記減算の残りである残留信号とから修正ベクト
ルを算出する手段と、前記修正ベクトルを用いて前記反
響路のインパルス応答の推定を逐次修正する手段と、前
記修正されたインパルス応答をもつ擬似反響路を生成す
る手段と、前記擬似反響路に前記受話信号を印加するこ
とで前記擬似反響信号を生成する手段とを有することを
特徴とする。
【0022】(作用)本発明では、反響路減衰特性の周
波数による相違を擬似反響路インパルス応答の逐次修正
に反映させるようにしているので、擬似反響路インパル
ス応答の反響路インパルス応答への収束速度の向上が得
られる。具体的には、図2のような行列Bを入力信号ベ
クトルに乗算することで、周波数ごとの減衰特性を調整
し推定する反響路インパルス応答の減衰特性を反映する
修正ベクトルを生成する。擬似反響路インパルス応答の
逐次修正は次の(式3)のようになる。
【0023】
【数4】
【0024】
【発明の実施の形態】以下、図面を参照して、本発明の
実施の形態を詳細に説明する。
【0025】(発明の原理)まず、最初に、本発明の原
理を説明する。学習同定法(NLMSアルゴリズム)で
は、擬似反響路インパルス応答の修正方向dW(k)は、
【0026】
【数5】 という拘束条件のもとで、評価関数
【0027】
【数6】 を最小にするベクトルとして、Lagrange(ラグランジ
ュ)の未定乗数法により導かれる。 そして、
【0028】
【数7】 を擬似反響路インパルス応答の更新式とする。ただし、
μ、W(k)、X(k)は(式1)で定義した通りであ
る。また、これら(式4)、(式5)、(式6)の詳細
については、例えば、S. Haykin,『Adaptive Filter T
heory』,Prentice-Hall,pp. 352―356,1991に記載さ
れている。しかし、この評価関数V1(k)を用いると、疑
似反響路インパルス応答における各係数が、振幅包絡が
指数減衰しているインパルス応答内での位置に依存して
残差に影響することが考慮されていない。
【0029】そこで、次の評価関数
【0030】
【数8】
【0031】
【数9】 という修正方向が導びかれ、(C+CT-1の設定によっ
て、反響路の減衰特性を修正方向に反映させることが可
能となる。ここで、
【0032】
【数10】 とおけば、上記の(式6)および(式8)から、
【0033】
【数11】 となり、上記の(式3)の擬似反響路インパルス応答の
更新式が導かれる。行列Bとして、固有値が全て正の対
称行列を選び、入力信号ベクトルに乗算させて修正方向
を求めるとき、行列
【0034】
【数12】 の固有値は全て正になり、(式7)のV2(k)が評価関数
の性質を満たす。このため、反響路インパルス応答の推
定に問題は生じない。行列Bの一例として、図2に示す
ような行列が考えられる。この行列は5つのL×Lブロッ
ク行列D(p)からなり、各行列D(p)は同図の下部に示すよ
うに、1×L行列、(L−2)×L行列、1×L行列の3つ
の行列に分割できる。行列Bの大きさは、5L×5Lにな
る。
【0035】修正方向BXの中で信号x(n)に対応する
成分bnは、行列D(p)の成分と受話信号の乗算により、
【0036】
【数13】 になる。この乗算は受信信号を伝達関数
【0037】
【数14】
【0038】のデジタルフィルタに通過させて1サンプ
ルずらすことと等価である。サンプリング周波数がFの
とき、0〜1の値をとる正規化周波数fr=2f/Fを
もちいると、このフィルタの周波数特性H(fr)は
【0039】
【数15】H(fr)=p+(1−2p)e-j π fr+pe-j2 π fr になる。例えばp=0.1,0.2では、図3に示すような周波
数特性|H(fr)|を持つ。
【0040】従って、例えば0<p1<p2<p3<p4
0.25のように、ブロック行列D(p)の設定をブロック
毎に変えことにより、後部ほど高域成分が低域成分より
も減衰している修正方向を得ることが可能となる。
【0041】また逆に、0.25>p1>p2>p3>p4
>0のように、ブロック行列D(p)の設定をブロック毎に
変えことにより、後部ほど低域成分が高域成分よりも減
衰している修正方向を得ることが可能となる。
【0042】(第1の実施形態)図1に本発明の第1の
実施形態の構成を示す。受話信号x(k)は再生器(スピー
カ)2により音響信号に再生される。収音器(マイクロ
ホン)3により収音された信号y(k)は、減算器7により
擬似反響路(FIRフィルタ;有限インパルスフィル
タ)6からの擬似反響信号u(k)が減算されて残留信
号e(k)が得られる。
【0043】
【数16】 推定部5の受話信号ベクトル生成部5aにおいて、受話信
号x(k)から次式で表わされる受話信号ベクトルX(k)が生
成され、
【0044】
【数17】
【0045】修正方向算出部5bにおいて、固有値が全て
正の対称行列Bを受話信号ベクトルX(k)に乗算させて修
正方向BX(k)が生成される。推定部5の修正情報生成部5
cでは、受話信号ベクトルX(k),修正方向BX(k),誤差
(残留信号)e(k)から、修正ベクトル
【0046】
【数18】 ただし、μ:ステップゲイン(0<μ<2)を計算し、
この計算結果とタップ係数記憶部5dに記憶されている前
回の擬似反響路推定値W(k)を加算器5eにおいて加算
して、次時刻の擬似反響路推定値W(k+1) を得る((式
9)を参照)。
【0047】推定部5で得られた次時刻の擬似反響路推
定値W(k+1) は擬似反響路6とタップ係数記憶部5dへ供
給され、擬似反響路6においてこの擬似反響路推定値W
(k+1)と受話信号ベクトルX(k+1)とから擬似反響信号u
(k+1)が生成される。
【0048】(第2の実施形態)図2に、本発明の第2
の実施形態における、ブロックに分割された非対角成分
をもつ行列の一例を挙げる。図2の行列は、非対角成分
をもつ5つのブロック行列D(p)からなり、
【0049】
【数19】 のとき、行列D(p)はさらに行方向に3つの行列に分割さ
れる。図2において、行列Bは全部で13のブロックに分
割される。このとき、受話信号ベクトルX(k)に行列D(p)
の中央の行列を乗算することは、受話信号x(k)を伝達関
【0050】
【数20】 のデジタルフィルタに通過させることと等価になってい
る。
【0051】(第3の実施形態)図4に、本発明の第3
の実施形態における、非対角成分を持つ行列Bによる受
話信号ベクトルX(k)から修正方向BX(k)を算出する演算
がシフトバッファにより実行される、修正方向算出法の
一例を示す。この図4の構成は、本発明の第2の実施形
態の修正方向BX(k)をシフトバッファにより構成したと
きの図である。
【0052】図4中のフィルタF1,F2,F3は、行列D(p)の
一部
【0053】
【数21】 をそれぞれ対応する受話信号ベクトルX(k)の一部に乗算
させることに対応している。
【0054】シフトバッファによる修正方向の算出は、
次のステップで行われる。 (a)入力信号の受話信号ベクトルX(k−1)を1つシフ
トする。 (b)第1〜第5のシフトバッファの値を1つシフトす
る。 (c)入力信号の受話信号x(k)を入力信号の受話信号ベ
クトルの先頭と第1シフトバッファの先頭に入れる。 (d)第2,3,4,5のシフトバッファの先頭に、対応す
る入力信号ベクトルのサンプルをFIRフィルタF2 により
処理した結果を入れる。 (e)第1〜第5のシフトバッファの値を修正方向BX(k)
に転送する。 (f)FIRフィルタF1, F3による処理結果の1サンプル
を、修正方向BX(k)中の対応するサンプルに入れる。
【0055】このように、シフト処理を用いることで、
行列の計算をすることなく、少ない演算量で修正方向BX
(k)を得ることができる。
【0056】(第4の実施形態)図5に、本発明の第4
の実施形態における、周波数帯域型の反響消去方法の一
例を示す。
【0057】推定部5および擬似反響路6へ供給される
サブバンド受話信号x(k)は、受話信号を帯域通過フィル
タ10に通過させてから、間引き部11においてあらかじめ
決められた間引き率で間引くことにより得られる。ま
た、収音器3からの信号を、帯域通過フィルタ12に通過
させ、あらかじめ決められた間引き率で間引き部13にお
いて間引くことで、サブバンド収音信号y(k)が得られ
る。
【0058】減算器7においてサブバンド収音信号y(k)
から擬似反響信号u(k)を減算して得られた反響消去後の
信号e(k)は、あらかじめ決めれらた補間率により補間部
14において補間され、帯域通過フィルタ15を通過するこ
とで、間引き前の周波数帯域の信号に戻す。
【0059】この第4の実施形態の反響消去方法に従っ
て、残響時間250msの部屋においてサンプリング周波数4
0kHzで実測した室内インパルス応答の10〜15kHz帯域の
伝達特性を用いて数値シミュレーションを行った。この
伝達特性は、室内インパルス応答を10〜15kHz帯域通過
フィルタで処理し、4分の1に間引いて、1500タップで
打ち切って得られたFIRフィルタでシミュレートした。
この反響路インパルス応答とスペクトログラムを図6と
図7に示す。収音器3のSN比を−35dBに設定した。適応
フィルタのタップ数は1250タップとした。
【0060】サブバンド受話信号から修正方向を算出す
る行列Bは、 250×250の行列
【0061】
【数22】 とr=0.6をもちいて次のように設定された1250×1
250の行列である。
【0062】
【数23】 受話信号として白色雑音を入力したときの、学習同定
法、ESアルゴリズムと本発明による実験結果の残留反響
レベルのグラフを図8に示す。残留反響信号のレベルが
収音信号と比較して−20dBに達した時間により収束速度
を評価すると、図8のグラフでは、本発明による収束速
度は学習同定法の約50%程度であり、ESアルゴリズムの
収束速度の約90%となっている。
【0063】(第5の実施形態)図9に、本発明の第5
の実施形態における反響消去方法を示す。反響を含む信
号y(k)は、M個の収音器31〜3Mにより収音された信号の
それぞれを対応するフィルタ201〜20Mに通過させて、加
算器21により加算することで得られる。
【0064】本実施形態では、複数の収音器31〜3M、複
数のフィルタ201〜20Mと加算器21を組み合わせる構成
により、指向性のある収音が可能になる。指向性のある
収音により、スピーカ2―マイクロホン3間の音響結合
量が低減され、テレビ会議が円滑に進行できるようにな
る。収音に指向性を持たせた場合には、反響路インパル
ス応答の減衰特性が周波数に大きく依存するため、反響
路減衰特性の周波数による相違を反映する本発明の反響
消去方法が有効となる。
【0065】(他の実施形態)なお、本発明の目的は、
前述した実施の形態の機能を実現するソフトウエアのプ
ログラムコードを記録した記録媒体(記憶媒体)を、シ
ステムあるいは装置に供給し、そのシステムあるいは装
置のコンピュータ(またはCPUやMPU)が記録媒体
に格納されたプログラムコードを読み出し、実行するこ
とによっても、達成されることは言うまでもない。この
場合、記録媒体から読み出されたプログラムコード自体
が前述した実施の形態の機能を実現することになり、そ
のプログラムコードを記録した記録媒体は本発明を構成
することになる。そのプログラムコードを記録し、また
テーブル等の変数データを記録する記録媒体としては、
例えばフロッピディスク、ハードディスク、光ディス
ク、光磁気ディスク、CD−ROMなどを用いことがで
きる。
【0066】
【発明の効果】以上説明したように、本発明によれば、
収音信号から擬似反響信号を減算し、反響路減衰特性の
周波数による相違を反映する行列を受話信号ベクトルに
乗算して修正方向を算出し、受話信号ベクトル,修正方
向,減算の残りとから修正ベクトルを算出し、この修正
ベクトルを用いて反響路のインパルス応答の推定を逐次
修正するようにしたので、適応フィルタによる真の反響
路インパルス応答への収束特性を向上させることができ
る。
【0067】
【図面の簡単な説明】
【図1】本発明の第1の実施形態の構成を示すブロック
図である。
【0068】
【図2】本発明の第2の実施形態におけるブロックに分
割された行列を示す説明図である。
【0069】
【図3】伝達関数p+(1−2p)z-1+pz-2のデジ
タルフィルタの周波数振幅特性を示すグラフである。
【0070】
【図4】本発明の第3の実施形態におけるシフトバッフ
ァにより構成される修正方向算出の方法を示す概念図で
ある。
【0071】
【図5】本発明の第4の実施形態における周波数帯域型
反響消去方法の構成を示すブロック図である。
【0072】
【図6】実際にスピーカと会議用マイクロホンによって
測定された10〜15kHzの反響路インパルス応答を示す
グラフである。
【0073】
【図7】実際にスピーカと会議用マイクロホンよって測
定された10〜15kHzの反響路インパルス応答のスペク
トログラムである。
【0074】
【図8】本発明の第4の形態において学習同定法とES法
と本発明の実験結果の残留反響レベルを比較したグラフ
である。
【0075】
【図9】本発明の第5実施形態における、複数の収音器
により収音され、それぞれFIRフィルタにより処理され
加算されて得られた信号を反響路の出力とする反響消去
方法の構成を示すブロック図である。
【0076】
【図10】従来の反響消去方法の構成を示すブロック図
である。
【0077】
【符号の説明】
1 受話入力端 2 再生器(スピーカ) 3 収音器(マイクロホン) 31〜3M 収音器 4 送話出力端 5 推定部(推定回路) 5a 受話信号ベクトル生成部 5b 修正方向算出部 5c 修正方情報生成部 5d タップ係数記憶部 5e 加算器 6 擬似反響路(FIRフィルタ) 7 減算器 10、12、15 帯域通過フィルタ 11、13 間引き部 14 補間部 201〜20M フィルタ 21 加算器
───────────────────────────────────────────────────── フロントページの続き (72)発明者 青木 茂明 東京都千代田区大手町二丁目3番1号 日 本電信電話株式会社内 Fターム(参考) 5K027 AA07 BB03 DD10 DD11 DD14 HH03 JJ03 5K046 BB01 HH24 HH52 HH79

Claims (10)

    【特許請求の範囲】
  1. 【請求項1】 受話信号を再生器で音響信号に再生する
    ステップと、 再生された前記音響信号を反響路を経由して収音器で収
    音するステップと、 収音された収音信号から擬似反響信号を減算するステッ
    プと、 反響路減衰特性の周波数による相違を反映する行列を前
    記受話信号ベクトルに乗算して修正方向を算出するステ
    ップと、 前記受話信号ベクトル,前記修正方向,前記減算の残り
    である残留信号とから修正ベクトルを算出するステップ
    と、 前記修正ベクトルを用いて前記反響路のインパルス応答
    の推定を逐次修正するステップと、 前記修正されたインパルス応答をもつ擬似反響路を生成
    するステップと、 前記擬似反響路に前記受話信号を印加することで前記擬
    似反響信号を生成するステップとを有することを特徴と
    する反響消去方法。
  2. 【請求項2】 周波数毎に異なる反響路減衰特性を反映
    する前記行列がKブロック(Kは2以上の正の整数)に
    分割されており、 前記受話信号ベクトルから前記修正方向の各係数を求め
    るFIRフィルタ係数が前記各ブロック毎に設定されて
    いることを特徴とする請求項1に記載の反響消去方法。
  3. 【請求項3】 前記受話信号ベクトルから前記修正方向
    を算出するための前記行列の乗算が、シフトバッファ
    と、該シフトバッファ間のFIRフィルタ、もしくは前
    記受話信号ベクトルの一部から前記シフトバッファへの
    FIRフィルタにより遂行されることを特徴とする請求
    項2に記載の反響消去方法。
  4. 【請求項4】 前記擬似反響路に印加される前記受話信
    号を帯域通過フィルタ通過後にあらかじめ決められた間
    引き率で間引かれたサブバンド受話信号とし、 前記収
    音器からの収音信号を帯域通過フィルタ通過後にあらか
    じめ決められた間引き率で間引かれたサブバンド収音信
    号とし、前記修正ベクトルが該サブバンド受話信号と該
    サブバンド収音信号から生成され、 反響消去後の信号をあらかじめ決めれらた補間率により
    補間して帯域通過フィルタを通過させることを特徴とす
    る請求項1ないし3のいずれかに記載の反響消去方法。
  5. 【請求項5】 M個(Mは2以上の正の整数)の収音器
    により収音された信号を対応するそれぞれM個のフィル
    タにより処理するステップと、 前記M個のフィルタにより処理された信号を加算するス
    テップとを有し、 該加算して得られた信号を前記反響路の出力である収音
    信号として取り扱うことを特徴とする請求項1ないし4
    のいずれかに記載の反響消去方法。
  6. 【請求項6】 受話信号を音響信号に再生する再生器
    と、 再生された前記音響信号を反響路を経由して収音する収
    音器と、 収音された収音信号から擬似反響信号を減算する手段
    と、 反響路減衰特性の周波数による相違を反映する行列を前
    記受話信号ベクトルに乗算して修正方向を算出する手段
    と、 前記受話信号ベクトル,前記修正方向,前記減算の残り
    である残留信号とから修正ベクトルを算出する手段と、 前記修正ベクトルを用いて前記反響路のインパルス応答
    の推定を逐次修正する手段と、 前記修正されたインパルス応答をもつ擬似反響路を生成
    する手段と、 前記擬似反響路に前記受話信号を印加することで前記擬
    似反響信号を生成する手段とを有することを特徴とする
    反響消去装置。
  7. 【請求項7】 周波数ごとに異なる反響路減衰特性を反
    映する前記行列がKブロック(Kは2以上の正の整数)
    に分割されており、 前記受話信号ベクトルから前記修正方向の各係数を求め
    るFIRフィルタ係数が前記各ブロック毎に設定されて
    いることを特徴とする請求項6に記載の反響消去装置。
  8. 【請求項8】 前記受話信号ベクトルから前記修正方向
    を算出するための前記行列の乗算演算が、シフトバッフ
    ァと、該シフトバッファ間のFIRフィルタ、もしくは
    前記受話信号ベクトルの一部から前記シフトバッファへ
    のFIRフィルタにより遂行されることを特徴とする請
    求項7に記載の反響消去装置。
  9. 【請求項9】 前記擬似反響路に印加される前記受話信
    号を帯域通過フィルタ通過後にあらかじめ決められた間
    引き率で間引かれたサブバンド受話信号とし、 前記収
    音器からの収音信号を帯域通過フィルタ通過後にあらか
    じめ決められた間引き率で間引かれたサブバンド収音信
    号とし、前記修正ベクトルが該サブバンド受話信号と該
    サブバンド収音信号から生成され、 反響消去後の信号をあらかじめ決めれらた補間率により
    補間して帯域通過フィルタを通過させることを特徴とす
    る請求項6ないし8のいずれかに記載の反響消去装置。
  10. 【請求項10】 M個(Mは2以上の正の整数)の収音
    器により収音された信号をそれぞれ処理するM個のフィ
    ルタと、 前記M個のフィルタにより処理された信号を加算する加
    算器とを有し、 該加算器で加算して得られた信号を前記反響路の出力で
    ある収音信号として取り扱うことを特徴とする請求項6
    ないし9のいずれかに記載の反響消去装置。
JP2000061265A 2000-03-06 2000-03-06 反響消去方法および反響消去装置 Pending JP2001251224A (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2000061265A JP2001251224A (ja) 2000-03-06 2000-03-06 反響消去方法および反響消去装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2000061265A JP2001251224A (ja) 2000-03-06 2000-03-06 反響消去方法および反響消去装置

Publications (1)

Publication Number Publication Date
JP2001251224A true JP2001251224A (ja) 2001-09-14

Family

ID=18581426

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2000061265A Pending JP2001251224A (ja) 2000-03-06 2000-03-06 反響消去方法および反響消去装置

Country Status (1)

Country Link
JP (1) JP2001251224A (ja)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005086377A (ja) * 2003-09-05 2005-03-31 Japan Radio Co Ltd 遅延波キャンセラ
WO2005091677A1 (ja) * 2004-03-22 2005-09-29 Toa Corporation マルチチャンネルシステム同定装置
JP2009271045A (ja) * 2008-05-07 2009-11-19 Sensitive Object 物体の表面上の衝撃位置を決定する方法
WO2018230062A1 (ja) * 2017-06-12 2018-12-20 株式会社オーディオテクニカ 音声信号処理装置と音声信号処理方法と音声信号処理プログラム

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005086377A (ja) * 2003-09-05 2005-03-31 Japan Radio Co Ltd 遅延波キャンセラ
JP4579518B2 (ja) * 2003-09-05 2010-11-10 日本無線株式会社 遅延波キャンセラ
WO2005091677A1 (ja) * 2004-03-22 2005-09-29 Toa Corporation マルチチャンネルシステム同定装置
JP2009271045A (ja) * 2008-05-07 2009-11-19 Sensitive Object 物体の表面上の衝撃位置を決定する方法
WO2018230062A1 (ja) * 2017-06-12 2018-12-20 株式会社オーディオテクニカ 音声信号処理装置と音声信号処理方法と音声信号処理プログラム
JPWO2018230062A1 (ja) * 2017-06-12 2020-05-21 株式会社オーディオテクニカ 音声信号処理装置と音声信号処理方法と音声信号処理プログラム
US11227618B2 (en) 2017-06-12 2022-01-18 Audio-Technica Corporation Sound signal processing device, sound signal processing method and sound signal processing program
JP7122756B2 (ja) 2017-06-12 2022-08-22 株式会社オーディオテクニカ 音声信号処理装置と音声信号処理方法と音声信号処理プログラム

Similar Documents

Publication Publication Date Title
JP5075042B2 (ja) エコー消去装置、エコー消去方法、そのプログラム、記録媒体
US7620172B2 (en) Method and system for eliminating noises and echo in voice signals
JP4567655B2 (ja) 音声信号中の暗騒音を抑制する方法及び装置、並びに反響消去を伴う対応する装置
US9866792B2 (en) Display apparatus and echo cancellation method thereof
JP3359460B2 (ja) 適応フィルタ及びエコーキャンセラ
JP2006174456A (ja) 受信オーディオを制限するためのシステム
CN110956975B (zh) 回声消除方法及装置
JPH04120812A (ja) Fir形エコーキャンセラ
KR20010043833A (ko) 스펙트럼 종속 지수 이득 함수 평균화를 이용한 스펙트럼공제에 의한 신호 잡음 저감
JP3507020B2 (ja) 反響抑圧方法、反響抑圧装置及び反響抑圧プログラム記憶媒体
CN112863532A (zh) 回音抑制装置、回音抑制方法以及存储介质
CN112509595A (zh) 音频数据处理方法、系统及存储介质
CN109215672B (zh) 一种声音信息的处理方法、装置及设备
WO2019098178A1 (ja) 音声コミュニケーション装置、音声コミュニケーション方法、プログラム
JP2003188776A (ja) 音響エコー消去方法、装置及び音響エコー消去プログラム
JP3403473B2 (ja) ステレオエコーキャンセラ
JP4119328B2 (ja) 収音方法、その装置、そのプログラム、およびその記録媒体。
JP2005051744A (ja) 通話装置
JP3864914B2 (ja) エコー抑圧装置
JP2001251224A (ja) 反響消去方法および反響消去装置
JP2003250193A (ja) 反響消去方法、この方法を実施する装置、プログラムおよびその記録媒体
CN115579016A (zh) 一种声学回声消除的方法和系统
JP4041770B2 (ja) 音響エコー消去方法、その装置、プログラム及びその記録媒体
JP3724705B2 (ja) 反響抑圧方法・反響抑圧装置
JP2008060715A (ja) エコー消去装置、その方法、そのプログラム、およびその記録媒体