JP2792311B2 - 多チャンネルエコー除去方法および装置 - Google Patents
多チャンネルエコー除去方法および装置Info
- Publication number
- JP2792311B2 JP2792311B2 JP4046185A JP4618592A JP2792311B2 JP 2792311 B2 JP2792311 B2 JP 2792311B2 JP 4046185 A JP4046185 A JP 4046185A JP 4618592 A JP4618592 A JP 4618592A JP 2792311 B2 JP2792311 B2 JP 2792311B2
- Authority
- JP
- Japan
- Prior art keywords
- signal
- delay time
- circuit
- signals
- time difference
- 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 - Fee Related
Links
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
Description
【０００１】
【産業上の利用分野】本発明は、複数の受信信号と単数
または複数の送信信号を有するシステムにおける、エコ
ー除去方法および装置に関するものである。
または複数の送信信号を有するシステムにおける、エコ
ー除去方法および装置に関するものである。
【０００２】
【従来の技術】複数の受信信号と単数または複数の送信
信号を有するシステムにおいて、受信信号が空間音響経
路を伝搬することによって生じるエコーを除去する多チ
ャンネルエコー除去方法あるいは装置に関しては、電子
情報通信学会技術研究報告Ｖｏｌ．８４，Ｎｏ．３３
０，ｐｐ．７−４，ＣＳ−８４−１７８（以下文献１）
において、縦続接続型および線形結合型が提案されてい
る。しかしながら、第６回ディジタル信号処理シンポジ
ウム講演論文集ｐｐ．１４４−１４９，Ａ５−３（以下
文献２）において、縦続接続型および線形結合型には、
ハードウエア規模がチャンネル数の自乗に比例するため
ハードウエア規模が大きくなる、受信信号が強い相互相
関を持つ場合には適応フィルタの収束が遅くなる、適応
フィルタの係数が最適値に収束しない場合がある、など
の問題があることが指摘されている。この問題を解決す
るために、１９９１年電子情報通信学会秋期大会講演論
文集第１巻ｐｐ．８３，Ａ−８３（以下文献３）におい
て、１チャンネル１個の適応フィルタで構成される多チ
ャンネルエコー除去装置が提案されている。以下、文献
３にしたがって、この多チャンネルエコー除去装置を受
信信号、送信信号ともに２チャンネルのテレビ会議シス
テムに適用した場合について説明する。
信号を有するシステムにおいて、受信信号が空間音響経
路を伝搬することによって生じるエコーを除去する多チ
ャンネルエコー除去方法あるいは装置に関しては、電子
情報通信学会技術研究報告Ｖｏｌ．８４，Ｎｏ．３３
０，ｐｐ．７−４，ＣＳ−８４−１７８（以下文献１）
において、縦続接続型および線形結合型が提案されてい
る。しかしながら、第６回ディジタル信号処理シンポジ
ウム講演論文集ｐｐ．１４４−１４９，Ａ５−３（以下
文献２）において、縦続接続型および線形結合型には、
ハードウエア規模がチャンネル数の自乗に比例するため
ハードウエア規模が大きくなる、受信信号が強い相互相
関を持つ場合には適応フィルタの収束が遅くなる、適応
フィルタの係数が最適値に収束しない場合がある、など
の問題があることが指摘されている。この問題を解決す
るために、１９９１年電子情報通信学会秋期大会講演論
文集第１巻ｐｐ．８３，Ａ−８３（以下文献３）におい
て、１チャンネル１個の適応フィルタで構成される多チ
ャンネルエコー除去装置が提案されている。以下、文献
３にしたがって、この多チャンネルエコー除去装置を受
信信号、送信信号ともに２チャンネルのテレビ会議シス
テムに適用した場合について説明する。
【０００３】図１５に、２つのテレビ会議室３０，３１
を結ぶ２チャンネルテレビ会議システムのうち、音声部
分のブロック図を示す。ここでは、第１のテレビ会議室
３０における音響エコー除去を例にとって説明する。
を結ぶ２チャンネルテレビ会議システムのうち、音声部
分のブロック図を示す。ここでは、第１のテレビ会議室
３０における音響エコー除去を例にとって説明する。
【０００４】第２のテレビ会議室３１に第２の話者１８
および第３の話者１９がいると仮定する。第２の話者１
８が発する音声２０および第３の話者１９が発する音声
２２は、空間音響経路を経て第３のマイク２４で収録さ
れて、第２のエコー除去装置１２０2 に供給される。第
３のマイク２４で収録された音声は、第１の受信信号１
として、第１のテレビ会議室３０に送られる。同様に、
第２の話者１８が発する音声２１および第３の話者１９
が発する音声２３は、空間音響経路を経て第４のマイク
２５で収録されて、第２のエコー除去装置１２０2 に供
給される。第４のマイク２５で収録された音声は、第２
の受信信号２として、第１のテレビ会議室３０に送られ
る。
および第３の話者１９がいると仮定する。第２の話者１
８が発する音声２０および第３の話者１９が発する音声
２２は、空間音響経路を経て第３のマイク２４で収録さ
れて、第２のエコー除去装置１２０2 に供給される。第
３のマイク２４で収録された音声は、第１の受信信号１
として、第１のテレビ会議室３０に送られる。同様に、
第２の話者１８が発する音声２１および第３の話者１９
が発する音声２３は、空間音響経路を経て第４のマイク
２５で収録されて、第２のエコー除去装置１２０2 に供
給される。第４のマイク２５で収録された音声は、第２
の受信信号２として、第１のテレビ会議室３０に送られ
る。
【０００５】第１のテレビ会議室３０において、第１の
受信信号１が第１のスピーカ３で再生され、空間音響経
路を経て第１のマイク９に至って生じる第１のエコー５
と、第２の受信信号２が第２のスピーカ４で再生され、
空間音響経路を経て第１のマイク９に至って生じる第２
のエコー６と、第１のマイク９に至った第１の話者１１
の発する音声である第１の送信信号１２が加算されて、
第１の混在信号１４となる。同様に、第１の受信信号１
が第１のスピーカ３で再生され、空間音響経路を経て第
２のマイク１０に至って生じる第３のエコー７と、第２
の受信信号２が第２のスピーカ４で再生され、空間音響
経路を経て第２のマイク１０に至って生じる第４のエコ
ー８と、第２のマイク１０に至った話者１１の発する音
声である第２の送信信号１３が加算されて、第２の混在
信号１５となる。第１および第２の混在信号１４，１５
に含まれるエコー５，６，７，８を除去するために、第
１のエコー除去装置１２０1 を使用する。
受信信号１が第１のスピーカ３で再生され、空間音響経
路を経て第１のマイク９に至って生じる第１のエコー５
と、第２の受信信号２が第２のスピーカ４で再生され、
空間音響経路を経て第１のマイク９に至って生じる第２
のエコー６と、第１のマイク９に至った第１の話者１１
の発する音声である第１の送信信号１２が加算されて、
第１の混在信号１４となる。同様に、第１の受信信号１
が第１のスピーカ３で再生され、空間音響経路を経て第
２のマイク１０に至って生じる第３のエコー７と、第２
の受信信号２が第２のスピーカ４で再生され、空間音響
経路を経て第２のマイク１０に至って生じる第４のエコ
ー８と、第２のマイク１０に至った話者１１の発する音
声である第２の送信信号１３が加算されて、第２の混在
信号１５となる。第１および第２の混在信号１４，１５
に含まれるエコー５，６，７，８を除去するために、第
１のエコー除去装置１２０1 を使用する。
【０００６】遅延時間差推定回路１０１は、第１の受信
信号１および第２の受信信号２を入力とし、２つの受信
信号間の遅延時間差１０２を推定し、検出器１２１に供
給する。検出器１２１は、２つの受信信号１，２のどち
らが遅延時間が短いかを検出して、遅延時間が短い受信
信号を指定する受信信号選択情報１０４としてセレクタ
１０７に供給する。セレクタ１０７は、第１および第２
の受信信号１，２を入力とし、２つの受信信号１，２の
うち、受信信号選択情報１０４が指定する受信信号１０
８を選択し、第１の適応フィルタ１２２および第２の適
応フィルタ１２３に供給する。
信号１および第２の受信信号２を入力とし、２つの受信
信号間の遅延時間差１０２を推定し、検出器１２１に供
給する。検出器１２１は、２つの受信信号１，２のどち
らが遅延時間が短いかを検出して、遅延時間が短い受信
信号を指定する受信信号選択情報１０４としてセレクタ
１０７に供給する。セレクタ１０７は、第１および第２
の受信信号１，２を入力とし、２つの受信信号１，２の
うち、受信信号選択情報１０４が指定する受信信号１０
８を選択し、第１の適応フィルタ１２２および第２の適
応フィルタ１２３に供給する。
【０００７】第１の適応フィルタ１２２は、セレクタ１
０７が選択した受信信号１０８を入力として、第１の混
在信号１４に含まれるエコーに対応するエコーレプリカ
１１１を生成し、第１の減算器１１３に供給する。第１
の減算器１１３は、第１の混在信号１４から、第１の適
応フィルタ１２２の出力であるエコーレプリカ１１１を
差し引き、その結果を第１の出力信号１６とする。第１
の適応フィルタ１２２は、第１の出力信号１６を最小と
するように制御される。
０７が選択した受信信号１０８を入力として、第１の混
在信号１４に含まれるエコーに対応するエコーレプリカ
１１１を生成し、第１の減算器１１３に供給する。第１
の減算器１１３は、第１の混在信号１４から、第１の適
応フィルタ１２２の出力であるエコーレプリカ１１１を
差し引き、その結果を第１の出力信号１６とする。第１
の適応フィルタ１２２は、第１の出力信号１６を最小と
するように制御される。
【０００８】同様に、第２の適応フィルタ１２３は、セ
レクタ１０７が選択した受信信号１０８を入力として、
第２の混在信号１５に含まれるエコーに対応するエコー
レプリカ１１２を生成し、第２の減算器１１４に供給す
る。第２の減算器１１４は、第２の混在信号１５から、
第２の適応フィルタ１２３の出力であるエコーレプリカ
１１２を差し引き、その結果を第２の出力信号１７とす
る。第２の適応フィルタ１２３は、第２の出力信号１７
を最小とするように制御される。
レクタ１０７が選択した受信信号１０８を入力として、
第２の混在信号１５に含まれるエコーに対応するエコー
レプリカ１１２を生成し、第２の減算器１１４に供給す
る。第２の減算器１１４は、第２の混在信号１５から、
第２の適応フィルタ１２３の出力であるエコーレプリカ
１１２を差し引き、その結果を第２の出力信号１７とす
る。第２の適応フィルタ１２３は、第２の出力信号１７
を最小とするように制御される。
【０００９】以下、第２のテレビ会議室３１で、第２の
話者１８および第３の話者１９が交互に発言しており、
第１の話者１１は沈黙している場合を仮定して、解析を
行なう。まず、第３のマイク２４の近くで第２の話者１
８だけが発言しており、第１の話者１１および第３の話
者１９は沈黙している時を考える。第２の話者１８の音
声２０，２１のスペクトルをＸ1(z)、第２の話者１８か
ら第３のマイク２４に至る空間音響経路の伝達関数をＧ
11(z) 、第２の話者１８から第４のマイク２５に至る空
間音響経路の伝達関数をＧ12(z) とすると、第１の受信
信号１のスペクトルＲ1(z)および第２の受信信号２のス
ペクトルＲ2(z)は、 Ｒ1(z)＝Ｇ11(z) Ｘ1(z) （１） Ｒ2(z)＝Ｇ12(z) Ｘ1(z) （２） で与えられる。第１のスピーカ３から第１のマイク９に
至る空間音響経路の伝達関数をＨ11(z) 、第２のスピー
カ４から第１のマイク９に至る空間音響経路の伝達関数
をＨ21(z) とすると、第１の混在信号１４に含まれるエ
コーのスペクトルＳ1(z)は、 Ｓ1(z)＝Ｈ11(z) Ｒ1(z)＋Ｈ21(z) Ｒ2(z) （３） で与えられる。（３）式に（１）式および（２）式を代
入してＲ1(z)およびＲ2(z)を消去すると、
話者１８および第３の話者１９が交互に発言しており、
第１の話者１１は沈黙している場合を仮定して、解析を
行なう。まず、第３のマイク２４の近くで第２の話者１
８だけが発言しており、第１の話者１１および第３の話
者１９は沈黙している時を考える。第２の話者１８の音
声２０，２１のスペクトルをＸ1(z)、第２の話者１８か
ら第３のマイク２４に至る空間音響経路の伝達関数をＧ
11(z) 、第２の話者１８から第４のマイク２５に至る空
間音響経路の伝達関数をＧ12(z) とすると、第１の受信
信号１のスペクトルＲ1(z)および第２の受信信号２のス
ペクトルＲ2(z)は、 Ｒ1(z)＝Ｇ11(z) Ｘ1(z) （１） Ｒ2(z)＝Ｇ12(z) Ｘ1(z) （２） で与えられる。第１のスピーカ３から第１のマイク９に
至る空間音響経路の伝達関数をＨ11(z) 、第２のスピー
カ４から第１のマイク９に至る空間音響経路の伝達関数
をＨ21(z) とすると、第１の混在信号１４に含まれるエ
コーのスペクトルＳ1(z)は、 Ｓ1(z)＝Ｈ11(z) Ｒ1(z)＋Ｈ21(z) Ｒ2(z) （３） で与えられる。（３）式に（１）式および（２）式を代
入してＲ1(z)およびＲ2(z)を消去すると、
【数１】 となる。
【００１０】同様に、第１のスピーカ３から第２のマイ
ク１０に至る空間音響経路の伝達関数をＨ12(z) 、第２
のスピーカ４から第２のマイク１０に至る空間音響経路
の伝達関数をＨ22(z) とすると、第２の混在信号１５に
含まれるエコーのスペクトルＳ2(z)は、 Ｓ2(z)＝Ｈ12(z) Ｒ1(z)＋Ｈ22(z) Ｒ2(z) （５） で与えられる。（５）式に（１）式および（２）式を代
入してＲ1(z)およびＲ2(z)を消去すると、
ク１０に至る空間音響経路の伝達関数をＨ12(z) 、第２
のスピーカ４から第２のマイク１０に至る空間音響経路
の伝達関数をＨ22(z) とすると、第２の混在信号１５に
含まれるエコーのスペクトルＳ2(z)は、 Ｓ2(z)＝Ｈ12(z) Ｒ1(z)＋Ｈ22(z) Ｒ2(z) （５） で与えられる。（５）式に（１）式および（２）式を代
入してＲ1(z)およびＲ2(z)を消去すると、
【数２】 となる。
【００１１】第２の話者１８は第３のマイク２４の近く
で発言しているので、第２の話者１８の音声２０が第１
の受信信号１として第１のテレビ会議室３０に送られる
のに要する遅延時間は、第２の話者１８の音声２１が第
２の受信信号２として第１のテレビ会議室３０に送られ
るのに要する遅延時間より短くなる。したがって、セレ
クタ１０７は第１の受信信号１を選択し、第１の適応フ
ィルタ１２２および第２の適応フィルタ１２３に供給す
る。
で発言しているので、第２の話者１８の音声２０が第１
の受信信号１として第１のテレビ会議室３０に送られる
のに要する遅延時間は、第２の話者１８の音声２１が第
２の受信信号２として第１のテレビ会議室３０に送られ
るのに要する遅延時間より短くなる。したがって、セレ
クタ１０７は第１の受信信号１を選択し、第１の適応フ
ィルタ１２２および第２の適応フィルタ１２３に供給す
る。
【００１２】
【数３】 となる。
【００１３】第１のエコー除去装置１２０1 の第１の出
力信号１６のスペクトルをＥ1(z)、第２の出力信号１７
のスペクトルをＥ2(z)とすると、
力信号１６のスペクトルをＥ1(z)、第２の出力信号１７
のスペクトルをＥ2(z)とすると、
【数４】
【数５】 となる。
【００１４】第１の適応フィルタ１２２および第２の適
応フィルタ１２３は、第１の出力信号１６および第２の
出力信号１７を最小とするように制御されるから、第１
の適応フィルタ１２２の伝達関数Ｗ1(z)の最適値
Ｗ1opt.(z)および第２の適応フィルタ１２３の伝達関数
Ｗ2(z)の最適値Ｗ2opt.(z)は、（１１),（１２）式よ
り、
応フィルタ１２３は、第１の出力信号１６および第２の
出力信号１７を最小とするように制御されるから、第１
の適応フィルタ１２２の伝達関数Ｗ1(z)の最適値
Ｗ1opt.(z)および第２の適応フィルタ１２３の伝達関数
Ｗ2(z)の最適値Ｗ2opt.(z)は、（１１),（１２）式よ
り、
【数６】 で与えられる。
【００１５】次に、同様の手順で、第４のマイク２５の
近くで第３の話者１９だけが発言しており、第１の話者
１１および第２の話者１８は沈黙している時の解析を行
なう。第３の話者１９の音声２２，２３のスペクトルを
Ｘ2(z)、第３の話者１９から第３のマイク２４に至る空
間音響経路の伝達関数をＧ21(z) 、第３の話者１９から
第４のマイク２５に至る空間音響経路の伝達関数をＧ22
(z) とすると、第１の受信信号１のスペクトルＲ1(z)お
よび第２の受信信号２のスペクトルＲ2(z)は、 Ｒ1(z)＝Ｇ21(z) Ｘ2(z) （１５） Ｒ2(z)＝Ｇ22(z) Ｘ2(z) （１６） で与えられる。第１の混在信号１４に含まれるエコーの
スペクトルＳ1(z)および第２の混在信号１５に含まれる
エコーのスペクトルＳ2(z)は、（３）式および（５）式
に（１５）式および（１６）式を代入してＲ1(z)および
Ｒ2(z)を消去すると、
近くで第３の話者１９だけが発言しており、第１の話者
１１および第２の話者１８は沈黙している時の解析を行
なう。第３の話者１９の音声２２，２３のスペクトルを
Ｘ2(z)、第３の話者１９から第３のマイク２４に至る空
間音響経路の伝達関数をＧ21(z) 、第３の話者１９から
第４のマイク２５に至る空間音響経路の伝達関数をＧ22
(z) とすると、第１の受信信号１のスペクトルＲ1(z)お
よび第２の受信信号２のスペクトルＲ2(z)は、 Ｒ1(z)＝Ｇ21(z) Ｘ2(z) （１５） Ｒ2(z)＝Ｇ22(z) Ｘ2(z) （１６） で与えられる。第１の混在信号１４に含まれるエコーの
スペクトルＳ1(z)および第２の混在信号１５に含まれる
エコーのスペクトルＳ2(z)は、（３）式および（５）式
に（１５）式および（１６）式を代入してＲ1(z)および
Ｒ2(z)を消去すると、
【数７】 となる。
【００１６】第３の話者１９は第４のマイク２５の近く
で発言しているので、第３の話者１９の音声２３が第２
の受信信号２として第１のテレビ会議室３０に送られる
のに要する遅延時間は、第３の話者１９の音声２２が第
１の受信信号１として第１のテレビ会議室３０に送られ
るのに要する遅延時間より短くなる。したがって、セレ
クタ１０７は第２の受信信号２を選択し、第１の適応フ
ィルタ１２２および第２の適応フィルタ１２３に供給す
る。
で発言しているので、第３の話者１９の音声２３が第２
の受信信号２として第１のテレビ会議室３０に送られる
のに要する遅延時間は、第３の話者１９の音声２２が第
１の受信信号１として第１のテレビ会議室３０に送られ
るのに要する遅延時間より短くなる。したがって、セレ
クタ１０７は第２の受信信号２を選択し、第１の適応フ
ィルタ１２２および第２の適応フィルタ１２３に供給す
る。
【数８】 となる。
【００１７】第１のエコー除去装置１２０1 の第１の出
力信号１６のスペクトルＥ1(z)および第２の出力信号１
７のスペクトルＥ2(z)は、
力信号１６のスペクトルＥ1(z)および第２の出力信号１
７のスペクトルＥ2(z)は、
【数９】 となる。第１の適応フィルタ１２２および第２の適応フ
ィルタ１２３は、第１の出力信号１６および第２の出力
信号１７を最小とするように制御されるから、（２３),
（２４）式より、第１の適応フィルタ１２２の伝達関数
Ｗ1(z)の最適値Ｗ1opt.(z)および第２の適応フィルタ１
２３の伝達関数Ｗ2(z)の最適値Ｗ2opt.(z)は、
ィルタ１２３は、第１の出力信号１６および第２の出力
信号１７を最小とするように制御されるから、（２３),
（２４）式より、第１の適応フィルタ１２２の伝達関数
Ｗ1(z)の最適値Ｗ1opt.(z)および第２の適応フィルタ１
２３の伝達関数Ｗ2(z)の最適値Ｗ2opt.(z)は、
【数１０】 で与えられる。
【００１８】一般に、第２の話者１８および第３の話者
１９から第３のマイク２４および第４のマイク２５に至
る空間音響経路の伝達関数Ｇ11(z),Ｇ12(z),Ｇ21(z),Ｇ
22(z) は互いに異なる。したがって、（１３）式と（２
５）式、（１４）式と（２６）式を比較することによ
り、第１の適応フィルタ１２２の伝達関数の最適値Ｗ1o
pt.(z)および第２の適応フィルタ１２３の伝達関数の最
適値Ｗ2opt.(z)は、第２の話者１８が発言している時と
第３の話者１９が発言している時では互いに異なること
がわかる。適応フィルタの伝達関数は、適応フィルタの
係数で定められるから、第１の適応フィルタ１２２のフ
ィルタ係数の最適値および第２の適応フィルタ１２３の
フィルタ係数の最適値は、第２の話者１８が発言してい
る時と第３の話者１９が発言している時では互いに異な
る。
１９から第３のマイク２４および第４のマイク２５に至
る空間音響経路の伝達関数Ｇ11(z),Ｇ12(z),Ｇ21(z),Ｇ
22(z) は互いに異なる。したがって、（１３）式と（２
５）式、（１４）式と（２６）式を比較することによ
り、第１の適応フィルタ１２２の伝達関数の最適値Ｗ1o
pt.(z)および第２の適応フィルタ１２３の伝達関数の最
適値Ｗ2opt.(z)は、第２の話者１８が発言している時と
第３の話者１９が発言している時では互いに異なること
がわかる。適応フィルタの伝達関数は、適応フィルタの
係数で定められるから、第１の適応フィルタ１２２のフ
ィルタ係数の最適値および第２の適応フィルタ１２３の
フィルタ係数の最適値は、第２の話者１８が発言してい
る時と第３の話者１９が発言している時では互いに異な
る。
【００１９】さらに、Ｇ11(z),Ｇ12(z),Ｇ21(z),Ｇ
22(z) は、話者の位置によって変化するから、話者が発
言しながら移動した時にも、第１の適応フィルタ１２２
のフィルタ係数の最適値および第２の適応フィルタ１２
３のフィルタ係数の最適値は変化する。
22(z) は、話者の位置によって変化するから、話者が発
言しながら移動した時にも、第１の適応フィルタ１２２
のフィルタ係数の最適値および第２の適応フィルタ１２
３のフィルタ係数の最適値は変化する。
【００２０】第１の適応フィルタ１２２のフィルタ係数
および第２の適応フィルタ１２３のフィルタ係数が各々
の最適値から外れると、第１の混在信号１４および第２
の混在信号１５に含まれるエコーとそれらに対応するエ
コーレプリカ１１１，１１２が一致しなくなり、エコー
を完全に消去できなくなる。したがって、話者が移動し
た時や、話者が交替した時は、第１の適応フィルタ１２
２のフィルタ係数および第２の適応フィルタ１２３のフ
ィルタ係数が各々の最適値に再び一致するまでの間、エ
コーを完全に消去できなくなる。
および第２の適応フィルタ１２３のフィルタ係数が各々
の最適値から外れると、第１の混在信号１４および第２
の混在信号１５に含まれるエコーとそれらに対応するエ
コーレプリカ１１１，１１２が一致しなくなり、エコー
を完全に消去できなくなる。したがって、話者が移動し
た時や、話者が交替した時は、第１の適応フィルタ１２
２のフィルタ係数および第２の適応フィルタ１２３のフ
ィルタ係数が各々の最適値に再び一致するまでの間、エ
コーを完全に消去できなくなる。
【００２１】
【発明が解決しようとする課題】以上、図１５を用いて
詳細に述べたように、従来の多チャンネルエコー除去方
法および装置は、一般に、話者が移動した時や話者が交
替した時からフィルタ係数が再び最適値に収束するまで
に長時間を要し、その間エコー除去性能が低下するとい
う問題があった。本発明の目的は、話者の移動や交替に
高速に追従できる多チャンネルエコー除去方法および装
置を提供することにある。
詳細に述べたように、従来の多チャンネルエコー除去方
法および装置は、一般に、話者が移動した時や話者が交
替した時からフィルタ係数が再び最適値に収束するまで
に長時間を要し、その間エコー除去性能が低下するとい
う問題があった。本発明の目的は、話者の移動や交替に
高速に追従できる多チャンネルエコー除去方法および装
置を提供することにある。
【００２２】
【課題を解決するための手段】第１の発明は、複数の受
信信号と単数または複数の送信信号を有するシステム
で、前記受信信号が空間音響経路を伝搬することによっ
て、又は回線の漏話等によって生じるエコーと前記送信
信号が混在する混在信号から前記エコーを除去する際
に、前記複数の受信信号から２信号を選択する複数の組
合せと１対１に対応させて、前記２信号が信号源から前
記システムに送られるまでに要する遅延時間の差で定義
される２信号間遅延時間差を複数個求め、前記複数の遅
延時間差に基づいて、前記複数の受信信号のうち前記遅
延時間が最も短い受信信号を検出し、前記遅延時間が最
も短い受信信号を入力とし、前記混在信号と１対１に対
応する適応フィルタによって生成された前記エコーのレ
プリカを対応する前記混在信号から差し引いたものを出
力信号とし、前記適応フィルタの係数は、前記出力信号
を最小とするように更新する多チャンネルエコー除去方
法において、前記適応フィルタは、複数組のフィルタ係
数を具備し、前記複数の遅延時間差に応じて使用するフ
ィルタ係数の組を切替えることを特徴としている。
信信号と単数または複数の送信信号を有するシステム
で、前記受信信号が空間音響経路を伝搬することによっ
て、又は回線の漏話等によって生じるエコーと前記送信
信号が混在する混在信号から前記エコーを除去する際
に、前記複数の受信信号から２信号を選択する複数の組
合せと１対１に対応させて、前記２信号が信号源から前
記システムに送られるまでに要する遅延時間の差で定義
される２信号間遅延時間差を複数個求め、前記複数の遅
延時間差に基づいて、前記複数の受信信号のうち前記遅
延時間が最も短い受信信号を検出し、前記遅延時間が最
も短い受信信号を入力とし、前記混在信号と１対１に対
応する適応フィルタによって生成された前記エコーのレ
プリカを対応する前記混在信号から差し引いたものを出
力信号とし、前記適応フィルタの係数は、前記出力信号
を最小とするように更新する多チャンネルエコー除去方
法において、前記適応フィルタは、複数組のフィルタ係
数を具備し、前記複数の遅延時間差に応じて使用するフ
ィルタ係数の組を切替えることを特徴としている。
【００２３】第２の発明は、第１の発明において前記フ
ィルタ係数の組を切替える際に、前記複数の遅延時間差
が変化した時刻から前記複数の遅延時間差が予め定めら
れた期間一定値を保つまでの間、前記フィルタ係数の組
を切替えないことを特徴としている。
ィルタ係数の組を切替える際に、前記複数の遅延時間差
が変化した時刻から前記複数の遅延時間差が予め定めら
れた期間一定値を保つまでの間、前記フィルタ係数の組
を切替えないことを特徴としている。
【００２４】第３の発明は、第１の発明または第２の発
明において前記フィルタ係数を更新する際に、前記複数
の遅延時間差が変化した時刻から、前記複数の遅延時間
差が予め定められた期間一定値を保つまでの間、前記適
応フィルタの係数更新を停止することを特徴としてい
る。
明において前記フィルタ係数を更新する際に、前記複数
の遅延時間差が変化した時刻から、前記複数の遅延時間
差が予め定められた期間一定値を保つまでの間、前記適
応フィルタの係数更新を停止することを特徴としてい
る。
【００２５】第４の発明は、第１、第２または第３の発
明において前記複数の遅延時間差を推定する際に、前記
複数の受信信号から任意の２信号を選択し、第１のトラ
ンスバーサル形適応フィルタで前記２信号のうちの第１
の受信信号を入力として第２の受信信号を予測し、第２
のトランスバーサル形適応フィルタで前記第２の受信信
号を入力として前記第１の受信信号を予測し、前記第１
の適応フィルタの係数と前記第２の適応フィルタの係数
のうち最も絶対値が大きいフィルタ係数から前記２信号
間遅延時間差を推定する操作を、前記複数の遅延時間差
を全て求めるまで繰り返すことを特徴としている。
明において前記複数の遅延時間差を推定する際に、前記
複数の受信信号から任意の２信号を選択し、第１のトラ
ンスバーサル形適応フィルタで前記２信号のうちの第１
の受信信号を入力として第２の受信信号を予測し、第２
のトランスバーサル形適応フィルタで前記第２の受信信
号を入力として前記第１の受信信号を予測し、前記第１
の適応フィルタの係数と前記第２の適応フィルタの係数
のうち最も絶対値が大きいフィルタ係数から前記２信号
間遅延時間差を推定する操作を、前記複数の遅延時間差
を全て求めるまで繰り返すことを特徴としている。
【００２６】第５の発明は、第１、第２または第３の発
明において前記複数の遅延時間差を推定する際に、前記
複数の受信信号から任意の２信号を選択し、前記２信号
間の相互相関関数の絶対値を最大とする時間差を求め、
前記時間差を前記２信号間遅延時間差の推定値とする操
作を、前記複数の遅延時間差を全て求めるまで繰り返す
ことを特徴としている。
明において前記複数の遅延時間差を推定する際に、前記
複数の受信信号から任意の２信号を選択し、前記２信号
間の相互相関関数の絶対値を最大とする時間差を求め、
前記時間差を前記２信号間遅延時間差の推定値とする操
作を、前記複数の遅延時間差を全て求めるまで繰り返す
ことを特徴としている。
【００２７】第６の発明は、第５の発明において前記２
信号間の相互相関関数を計算する際に、前記２信号をそ
れぞれ第１及び第２の信号とし、第１の時刻における前
記第１の信号の値と前記第１の時刻と予め定められた時
間差を有する第２の時刻における前記第２の信号の値と
の乗算をし、前記第１の時刻を予め定められた複数の時
刻に変化させて前記乗算を行なって得られた複数の乗算
結果の重みつき移動平均を求め、前記重みつき移動平均
値をもって前記２信号の前記複数の時間差に対する相互
相関関数の値とすることを特徴としている。
信号間の相互相関関数を計算する際に、前記２信号をそ
れぞれ第１及び第２の信号とし、第１の時刻における前
記第１の信号の値と前記第１の時刻と予め定められた時
間差を有する第２の時刻における前記第２の信号の値と
の乗算をし、前記第１の時刻を予め定められた複数の時
刻に変化させて前記乗算を行なって得られた複数の乗算
結果の重みつき移動平均を求め、前記重みつき移動平均
値をもって前記２信号の前記複数の時間差に対する相互
相関関数の値とすることを特徴としている。
【００２８】第７の発明は、第６の発明における前記重
みつき移動平均を再帰形の積分に置き換えたことを特徴
としている。
みつき移動平均を再帰形の積分に置き換えたことを特徴
としている。
【００２９】第８の発明は、複数の受信信号と単数また
は複数の送信信号を有するシステムで、前記受信信号が
空間音響経路を伝搬することによって、又は回線の漏話
等によって生じるエコーと前記送信信号が混在する混在
信号から前記エコーを除去する際に、前記複数の受信信
号全てを入力し、前記複数の遅延時間差を推定し出力す
る遅延時間差推定回路と、前記遅延時間差推定回路の推
定結果に基づいて、前記受信信号のうち前記遅延時間が
最も短いものを検出して受信信号選択情報とし、前記複
数の適応フィルタが使用するフィルタ係数の組を決定し
てフィルタ係数選択信号とし、前記複数の適応フィルタ
が係数更新を行なうかどうかを決定して係数更新制御信
号とする制御回路と、前記受信信号選択情報に基づいて
入力された前記複数の受信信号のうち前記遅延時間が最
も短いものを選択して出力するセレクタと、前記混在信
号と１対１に対応し、前記フィルタ係数選択信号に基づ
いて複数のフィルタ係数の組の中から使用するフィルタ
係数の組を選択し、前記セレクタで選択された受信信号
を入力して前記エコーのレプリカを生成する複数の適応
フィルタと、前記混在信号と１対１に対応し、前記生成
されたエコーのレプリカを前記混在信号から差し引いて
エコー除去装置の出力信号とする複数の減算器とを少く
とも具備し、前記複数の適応フィルタは、前記係数更新
制御信号に基づいて前記フィルタ係数選択信号で指定さ
れたフィルタ係数の組を前記複数の減算器の出力を最小
とするように更新することを特徴としている。
は複数の送信信号を有するシステムで、前記受信信号が
空間音響経路を伝搬することによって、又は回線の漏話
等によって生じるエコーと前記送信信号が混在する混在
信号から前記エコーを除去する際に、前記複数の受信信
号全てを入力し、前記複数の遅延時間差を推定し出力す
る遅延時間差推定回路と、前記遅延時間差推定回路の推
定結果に基づいて、前記受信信号のうち前記遅延時間が
最も短いものを検出して受信信号選択情報とし、前記複
数の適応フィルタが使用するフィルタ係数の組を決定し
てフィルタ係数選択信号とし、前記複数の適応フィルタ
が係数更新を行なうかどうかを決定して係数更新制御信
号とする制御回路と、前記受信信号選択情報に基づいて
入力された前記複数の受信信号のうち前記遅延時間が最
も短いものを選択して出力するセレクタと、前記混在信
号と１対１に対応し、前記フィルタ係数選択信号に基づ
いて複数のフィルタ係数の組の中から使用するフィルタ
係数の組を選択し、前記セレクタで選択された受信信号
を入力して前記エコーのレプリカを生成する複数の適応
フィルタと、前記混在信号と１対１に対応し、前記生成
されたエコーのレプリカを前記混在信号から差し引いて
エコー除去装置の出力信号とする複数の減算器とを少く
とも具備し、前記複数の適応フィルタは、前記係数更新
制御信号に基づいて前記フィルタ係数選択信号で指定さ
れたフィルタ係数の組を前記複数の減算器の出力を最小
とするように更新することを特徴としている。
【００３０】第９の発明は、第８の発明における制御回
路が、前記遅延時間差推定回路の推定結果に基づいて前
記受信信号のうち前記遅延時間が最も短いものを検出し
て受信信号選択情報とする受信信号選択回路と、前記受
信信号選択情報と前記遅延時間差推定回路の推定結果に
基づいて前記複数の適応フィルタが使用するフィルタ係
数の組を選択してフィルタ係数選択信号とするフィルタ
係数選択回路と、前記複数の適応フィルタが常に係数更
新を行なうように、係数更新制御信号を出力する係数更
新制御回路から構成されることを特徴としている。
路が、前記遅延時間差推定回路の推定結果に基づいて前
記受信信号のうち前記遅延時間が最も短いものを検出し
て受信信号選択情報とする受信信号選択回路と、前記受
信信号選択情報と前記遅延時間差推定回路の推定結果に
基づいて前記複数の適応フィルタが使用するフィルタ係
数の組を選択してフィルタ係数選択信号とするフィルタ
係数選択回路と、前記複数の適応フィルタが常に係数更
新を行なうように、係数更新制御信号を出力する係数更
新制御回路から構成されることを特徴としている。
【００３１】第１０の発明は、第８の発明における制御
回路が、前記遅延時間差推定回路の推定結果を遅延させ
る遅延器と、前記遅延器の出力と前記遅延時間差推定回
路の推定結果を比較することによって前記遅延時間差推
定回路の推定結果の変化を検出する比較器と、前記比較
器が最後に前記遅延時間差推定回路の推定結果の変化を
検出した時からの経過時間を算出し、前記経過時間が予
め定められた時間に達したかどうかを判定する経過時間
判定回路と、前記経過時間判定回路が前記経過時間が予
め定められた時間に達したと判定した時に、前記遅延時
間差推定回路の推定結果に基づいて前記受信信号のうち
最も前記遅延時間が短いものを検出して受信信号選択情
報とする受信信号選択回路と、前記経過時間判定回路が
前記経過時間が予め定められた時間に達したと判定した
時に、前記受信信号選択情報と前記遅延時間差推定回路
の推定結果に基づいて前記複数の適応フィルタが使用す
るフィルタ係数の組を選択してフィルタ係数選択信号と
するフィルタ係数選択回路と、前記複数の適応フィルタ
に常に係数更新を行なうように係数更新制御信号を出力
する係数更新制御回路とから構成されることを特徴とし
ている。
回路が、前記遅延時間差推定回路の推定結果を遅延させ
る遅延器と、前記遅延器の出力と前記遅延時間差推定回
路の推定結果を比較することによって前記遅延時間差推
定回路の推定結果の変化を検出する比較器と、前記比較
器が最後に前記遅延時間差推定回路の推定結果の変化を
検出した時からの経過時間を算出し、前記経過時間が予
め定められた時間に達したかどうかを判定する経過時間
判定回路と、前記経過時間判定回路が前記経過時間が予
め定められた時間に達したと判定した時に、前記遅延時
間差推定回路の推定結果に基づいて前記受信信号のうち
最も前記遅延時間が短いものを検出して受信信号選択情
報とする受信信号選択回路と、前記経過時間判定回路が
前記経過時間が予め定められた時間に達したと判定した
時に、前記受信信号選択情報と前記遅延時間差推定回路
の推定結果に基づいて前記複数の適応フィルタが使用す
るフィルタ係数の組を選択してフィルタ係数選択信号と
するフィルタ係数選択回路と、前記複数の適応フィルタ
に常に係数更新を行なうように係数更新制御信号を出力
する係数更新制御回路とから構成されることを特徴とし
ている。
【００３２】第１１の発明は、第８の発明における制御
回路が、前記遅延時間差推定回路の推定結果に基づいて
前記受信信号のうち前記遅延時間が最も短いものを検出
して受信信号選択情報とする受信信号選択回路と、前記
受信信号選択情報と前記遅延時間差推定回路の推定結果
に基づいて前記複数の適応フィルタが使用するフィルタ
係数の組を選択してフィルタ係数選択信号とするフィル
タ係数選択回路と、前記遅延時間差推定回路の推定結果
を遅延させる遅延器と、前記遅延器の出力と前記遅延時
間差推定回路の推定結果を比較することによって前記時
間差推定回路の推定結果の変化を検出する比較器と、前
記比較器が最後に前記遅延時間差推定回路の推定結果の
変化を検出した時からの経過時間を算出し、前記経過時
間が予め定められた時間に達したかどうかを判定する経
過時間判定回路と、前記比較器が前記遅延時間差推定回
路の推定結果の変化を検出した時にフィルタ係数更新を
停止する係数更新制御信号を、前記経過時間判定回路が
前記経過時間が予め定められた時間に達したと判定した
時にフィルタ係数更新を開始する係数更新制御信号を出
力する係数更新制御回路とから構成されることを特徴と
している。
回路が、前記遅延時間差推定回路の推定結果に基づいて
前記受信信号のうち前記遅延時間が最も短いものを検出
して受信信号選択情報とする受信信号選択回路と、前記
受信信号選択情報と前記遅延時間差推定回路の推定結果
に基づいて前記複数の適応フィルタが使用するフィルタ
係数の組を選択してフィルタ係数選択信号とするフィル
タ係数選択回路と、前記遅延時間差推定回路の推定結果
を遅延させる遅延器と、前記遅延器の出力と前記遅延時
間差推定回路の推定結果を比較することによって前記時
間差推定回路の推定結果の変化を検出する比較器と、前
記比較器が最後に前記遅延時間差推定回路の推定結果の
変化を検出した時からの経過時間を算出し、前記経過時
間が予め定められた時間に達したかどうかを判定する経
過時間判定回路と、前記比較器が前記遅延時間差推定回
路の推定結果の変化を検出した時にフィルタ係数更新を
停止する係数更新制御信号を、前記経過時間判定回路が
前記経過時間が予め定められた時間に達したと判定した
時にフィルタ係数更新を開始する係数更新制御信号を出
力する係数更新制御回路とから構成されることを特徴と
している。
【００３３】第１２の発明は、第８の発明における制御
回路が、前記遅延時間差推定回路の推定結果を遅延させ
る遅延器と、前記遅延器の出力と前記遅延時間差推定回
路の推定結果を比較することによって前記遅延時間差推
定回路の推定結果の変化を検出する比較器と、前記比較
器が最後に前記遅延時間差推定回路の推定結果の変化を
検出した時からの経過時間を算出し、前記経過時間が予
め定められた第１の時間に達したかどうかを判定する第
１の経過時間判定回路と、前記第１の経過時間判定回路
が前記経過時間が予め定められた第１の時間に達したと
判定した時に、前記遅延時間差推定回路の推定結果に基
づいて前記受信信号のうち最も前記遅延時間が短いもの
を検出して受信信号選択情報とする受信信号選択回路
と、前記第１の経過時間判定回路が前記経過時間が予め
定められた第１の時間に達したと判定した時に、前記受
信信号選択情報と前記遅延時間差推定回路の推定結果に
基づいて前記複数の適応フィルタが使用するフィルタ係
数の組を選択してフィルタ係数選択信号とするフィルタ
係数選択回路と、前記比較器が最後に前記遅延時間差推
定回路の推定結果の変化を検出した時からの経過時間を
算出し、前記経過時間が予め定められた第２の時間に達
したかどうかを判定する第２の経過時間判定回路と、前
記比較器が前記遅延時間差推定回路の推定結果の変化を
検出した時にフィルタ係数更新を停止する係数更新制御
信号を、前記第２の経過時間判定回路が前記経過時間が
予め定められた第２の時間に達したと判定した時にフィ
ルタ係数更新を開始する係数更新制御信号を出力する係
数更新制御回路とで構成されることを特徴としている。
回路が、前記遅延時間差推定回路の推定結果を遅延させ
る遅延器と、前記遅延器の出力と前記遅延時間差推定回
路の推定結果を比較することによって前記遅延時間差推
定回路の推定結果の変化を検出する比較器と、前記比較
器が最後に前記遅延時間差推定回路の推定結果の変化を
検出した時からの経過時間を算出し、前記経過時間が予
め定められた第１の時間に達したかどうかを判定する第
１の経過時間判定回路と、前記第１の経過時間判定回路
が前記経過時間が予め定められた第１の時間に達したと
判定した時に、前記遅延時間差推定回路の推定結果に基
づいて前記受信信号のうち最も前記遅延時間が短いもの
を検出して受信信号選択情報とする受信信号選択回路
と、前記第１の経過時間判定回路が前記経過時間が予め
定められた第１の時間に達したと判定した時に、前記受
信信号選択情報と前記遅延時間差推定回路の推定結果に
基づいて前記複数の適応フィルタが使用するフィルタ係
数の組を選択してフィルタ係数選択信号とするフィルタ
係数選択回路と、前記比較器が最後に前記遅延時間差推
定回路の推定結果の変化を検出した時からの経過時間を
算出し、前記経過時間が予め定められた第２の時間に達
したかどうかを判定する第２の経過時間判定回路と、前
記比較器が前記遅延時間差推定回路の推定結果の変化を
検出した時にフィルタ係数更新を停止する係数更新制御
信号を、前記第２の経過時間判定回路が前記経過時間が
予め定められた第２の時間に達したと判定した時にフィ
ルタ係数更新を開始する係数更新制御信号を出力する係
数更新制御回路とで構成されることを特徴としている。
【００３４】第１３の発明は、第８、第９、第１０、第
１１または第１２の発明における遅延時間差推定回路
が、前記複数の受信信号を入力とし、前記時間差を推定
する２つの受信信号を指定する制御信号に応答し、前記
指定された２信号間の遅延時間差を推定し、推定結果を
出力する少なくとも１つの２信号間遅延時間差推定回路
と、前記２信号間遅延時間差推定回路の各々の出力に応
答し、前記複数の時間差を求めるように前記制御信号を
前記２信号間遅延時間差推定回路の各々に出力し、前記
複数の遅延時間差推定結果を出力する制御回路とで構成
されることを特徴としている。
１１または第１２の発明における遅延時間差推定回路
が、前記複数の受信信号を入力とし、前記時間差を推定
する２つの受信信号を指定する制御信号に応答し、前記
指定された２信号間の遅延時間差を推定し、推定結果を
出力する少なくとも１つの２信号間遅延時間差推定回路
と、前記２信号間遅延時間差推定回路の各々の出力に応
答し、前記複数の時間差を求めるように前記制御信号を
前記２信号間遅延時間差推定回路の各々に出力し、前記
複数の遅延時間差推定結果を出力する制御回路とで構成
されることを特徴としている。
【００３５】第１４の発明は、第１３の発明における２
信号間遅延時間差推定回路が、前記複数の受信信号を入
力とし、前記制御信号で指定された２信号の一方を第１
の受信信号として出力する第１のセレクタと、前記複数
の受信信号を入力とし、前記制御信号で指定された２信
号の他方を第２の受信信号として出力する第２のセレク
タと、前記第１の受信信号を入力し、前記第２の受信信
号を予測する第１のトランスバーサル形適応フィルタ
と、前記第２の受信信号を入力し、前記第１の受信信号
を予測する第２のトランスバーサル形適応フィルタと、
前記第１の適応フィルタの係数の絶対値を求める複数の
絶対値計算回路からなる第１の絶対値計算回路群と、前
記第２の適応フィルタの係数の絶対値を求める複数の絶
対値計算回路からなる第２の絶対値計算回路群と、前記
第１および第２の絶対値計算回路群の各絶対値計算回路
の出力に基づいて前記第１および第２の受信信号間の遅
延時間差を推定する判定器とから構成されることを特徴
としている。
信号間遅延時間差推定回路が、前記複数の受信信号を入
力とし、前記制御信号で指定された２信号の一方を第１
の受信信号として出力する第１のセレクタと、前記複数
の受信信号を入力とし、前記制御信号で指定された２信
号の他方を第２の受信信号として出力する第２のセレク
タと、前記第１の受信信号を入力し、前記第２の受信信
号を予測する第１のトランスバーサル形適応フィルタ
と、前記第２の受信信号を入力し、前記第１の受信信号
を予測する第２のトランスバーサル形適応フィルタと、
前記第１の適応フィルタの係数の絶対値を求める複数の
絶対値計算回路からなる第１の絶対値計算回路群と、前
記第２の適応フィルタの係数の絶対値を求める複数の絶
対値計算回路からなる第２の絶対値計算回路群と、前記
第１および第２の絶対値計算回路群の各絶対値計算回路
の出力に基づいて前記第１および第２の受信信号間の遅
延時間差を推定する判定器とから構成されることを特徴
としている。
【００３６】第１５の発明は、第１３の発明における２
信号間遅延時間差推定回路が、前記複数の受信信号を入
力とし、前記制御信号で指定された２信号の一方を第１
の受信信号として出力する第１のセレクタと、前記複数
の受信信号を入力とし、前記制御信号で指定された２信
号の他方を第２の受信信号として出力する第２のセレク
タと、前記第１の受信信号と第２の受信信号の相互相関
関数のうち、予め定められた複数の時間差に対応する値
を計算する相互相関関数計算回路と、前記複数の時間差
に対応する相互相関関数の絶対値を求める複数の絶対値
計算回路からなる絶対値計算回路群と、前記絶対値計算
回路群の各絶対値計算回路の出力に基づいて前記相互相
関関数の絶対値を最大とする時間差を前記第１および第
２の受信信号間遅延時間差の推定値として出力する判定
器とから構成されることを特徴としている。
信号間遅延時間差推定回路が、前記複数の受信信号を入
力とし、前記制御信号で指定された２信号の一方を第１
の受信信号として出力する第１のセレクタと、前記複数
の受信信号を入力とし、前記制御信号で指定された２信
号の他方を第２の受信信号として出力する第２のセレク
タと、前記第１の受信信号と第２の受信信号の相互相関
関数のうち、予め定められた複数の時間差に対応する値
を計算する相互相関関数計算回路と、前記複数の時間差
に対応する相互相関関数の絶対値を求める複数の絶対値
計算回路からなる絶対値計算回路群と、前記絶対値計算
回路群の各絶対値計算回路の出力に基づいて前記相互相
関関数の絶対値を最大とする時間差を前記第１および第
２の受信信号間遅延時間差の推定値として出力する判定
器とから構成されることを特徴としている。
【００３７】第１６の発明は、第１５の発明における相
互相関関数計算回路が、前記第１の受信信号を遅延させ
る第１のタップドディレイラインと、前記第２の受信信
号を遅延させる第２のタップドディレイラインと、前記
第１のタップドディレイラインの各タップ出力と前記第
２の受信信号を互いに乗算する複数の乗算器からなる第
１の乗算器群と、前記第１の乗算器群の各乗算器と１対
１に対応し、各乗算器の出力を積分する複数の積分器か
らなる第１の積分器群と、前記第２のタップドディレイ
ラインの各タップ出力と前記第１の受信信号を互いに乗
算する複数の乗算器からなる第２の乗算器群と、前記第
２の乗算器群の各乗算器と１対１に対応し、各乗算器出
力を積分する複数の積分器からなる第２の積分器群と、
前記第１および第２の受信信号を互いに乗算する第１の
乗算器と、前記第１の乗算器の出力を積分する第１の積
分器とから構成され、前記第１および第２の積分器群と
前記第１の積分器の出力を前記第１および第２の受信信
号の相互相関関数として出力することを特徴としてい
る。
互相関関数計算回路が、前記第１の受信信号を遅延させ
る第１のタップドディレイラインと、前記第２の受信信
号を遅延させる第２のタップドディレイラインと、前記
第１のタップドディレイラインの各タップ出力と前記第
２の受信信号を互いに乗算する複数の乗算器からなる第
１の乗算器群と、前記第１の乗算器群の各乗算器と１対
１に対応し、各乗算器の出力を積分する複数の積分器か
らなる第１の積分器群と、前記第２のタップドディレイ
ラインの各タップ出力と前記第１の受信信号を互いに乗
算する複数の乗算器からなる第２の乗算器群と、前記第
２の乗算器群の各乗算器と１対１に対応し、各乗算器出
力を積分する複数の積分器からなる第２の積分器群と、
前記第１および第２の受信信号を互いに乗算する第１の
乗算器と、前記第１の乗算器の出力を積分する第１の積
分器とから構成され、前記第１および第２の積分器群と
前記第１の積分器の出力を前記第１および第２の受信信
号の相互相関関数として出力することを特徴としてい
る。
【００３８】第１７の発明は、第１６の発明における積
分器が、入力信号を入力する入力端と、出力信号を出力
する出力端と、前記入力信号を遅延させるタップドディ
レイラインと、前記タップドディレイラインの各タップ
出力を定数倍する複数の係数乗算器と、前記複数の係数
乗算器の出力の総和を求め、前記総和を積分結果を表す
前記出力信号として出力する加算器とから構成されるこ
とを特徴としている。
分器が、入力信号を入力する入力端と、出力信号を出力
する出力端と、前記入力信号を遅延させるタップドディ
レイラインと、前記タップドディレイラインの各タップ
出力を定数倍する複数の係数乗算器と、前記複数の係数
乗算器の出力の総和を求め、前記総和を積分結果を表す
前記出力信号として出力する加算器とから構成されるこ
とを特徴としている。
【００３９】第１８の発明は、第１６の発明における積
分器が、入力信号を入力する入力端と、出力信号を出力
する出力端と、前記入力信号を遅延させる第１の遅延器
と、１サンプル前の出力信号を格納する第２の遅延器
と、前記第２の遅延器の値と前記入力信号を加え、前記
第１の遅延器の出力を差し引いたものを前記出力信号と
すると共に、前記第２の遅延器に格納する加算器とから
構成されることを特徴としている。
分器が、入力信号を入力する入力端と、出力信号を出力
する出力端と、前記入力信号を遅延させる第１の遅延器
と、１サンプル前の出力信号を格納する第２の遅延器
と、前記第２の遅延器の値と前記入力信号を加え、前記
第１の遅延器の出力を差し引いたものを前記出力信号と
すると共に、前記第２の遅延器に格納する加算器とから
構成されることを特徴としている。
【００４０】第１９の発明は、第１６の発明における積
分器が、入力信号を入力する入力端と、出力信号を出力
する出力端と、前記入力信号を定数倍する第１の係数乗
算器と、前記出力信号を遅延させるタップドディレイラ
インと、前記タップドディレイラインの各タップ出力を
定数倍する複数の係数乗算器と、前記複数の係数乗算器
の出力と前記第１の係数乗算器の出力の総和を求め、前
記総和を前記出力信号とすると共に前記タップドディレ
イラインに格納する加算器とから構成されることを特徴
としている。
分器が、入力信号を入力する入力端と、出力信号を出力
する出力端と、前記入力信号を定数倍する第１の係数乗
算器と、前記出力信号を遅延させるタップドディレイラ
インと、前記タップドディレイラインの各タップ出力を
定数倍する複数の係数乗算器と、前記複数の係数乗算器
の出力と前記第１の係数乗算器の出力の総和を求め、前
記総和を前記出力信号とすると共に前記タップドディレ
イラインに格納する加算器とから構成されることを特徴
としている。
【００４１】
【作用】本発明の多チャンネルエコー除去方法および装
置は、適応フィルタの係数を複数組用意し、話者の位置
に応じてフィルタ係数の組を切替えて使用することによ
って、話者の移動や交替の際にフィルタ係数を再び最適
値に収束させる必要がなくなるので、話者の移動や交替
に高速に追従することができる。
置は、適応フィルタの係数を複数組用意し、話者の位置
に応じてフィルタ係数の組を切替えて使用することによ
って、話者の移動や交替の際にフィルタ係数を再び最適
値に収束させる必要がなくなるので、話者の移動や交替
に高速に追従することができる。
【００４２】
【実施例】受信信号がスピーカから空間音響経路を伝搬
してマイクで収録される音響エコーを除去する音響エコ
ー除去を例にとり、図１から図１４を参照して、本発明
の実施例を詳細に説明する。
してマイクで収録される音響エコーを除去する音響エコ
ー除去を例にとり、図１から図１４を参照して、本発明
の実施例を詳細に説明する。
【００４３】図１に、本発明において受信信号および送
信信号が２チャンネルである場合の実施例を示す。遅延
時間差推定回路１０１は、第１の受信信号１および第２
の受信信号２を入力とし、２つの受信信号１，２の間の
遅延時間差１０２を推定し、制御回路１０３に供給す
る。
信信号が２チャンネルである場合の実施例を示す。遅延
時間差推定回路１０１は、第１の受信信号１および第２
の受信信号２を入力とし、２つの受信信号１，２の間の
遅延時間差１０２を推定し、制御回路１０３に供給す
る。
【００４４】制御回路１０３は、２信号間遅延時間差１
０２を入力とし、受信信号選択情報１０４、フィルタ係
数選択信号１０５、係数更新制御信号１０６を出力す
る。まず、２信号間遅延時間差１０２に基づいて２つの
受信信号１，２のどちらの遅延時間が短いかを検出し、
受信信号選択情報１０４としてセレクタ１０７に供給す
る。次に、２信号間遅延時間差１０２に基づいて、第１
の適応フィルタ１０９および第２の適応フィルタ１１０
が使用するフィルタ係数を選択し、フィルタ係数選択信
号１０５として第１の適応フィルタ１０９および第２の
適応フィルタ１１０に供給する。最後に、フィルタ係数
更新を行なうかどうかを決定して、係数更新制御信号１
０６として第１の適応フィルタ１０９および第２の適応
フィルタ１１０に供給する。
０２を入力とし、受信信号選択情報１０４、フィルタ係
数選択信号１０５、係数更新制御信号１０６を出力す
る。まず、２信号間遅延時間差１０２に基づいて２つの
受信信号１，２のどちらの遅延時間が短いかを検出し、
受信信号選択情報１０４としてセレクタ１０７に供給す
る。次に、２信号間遅延時間差１０２に基づいて、第１
の適応フィルタ１０９および第２の適応フィルタ１１０
が使用するフィルタ係数を選択し、フィルタ係数選択信
号１０５として第１の適応フィルタ１０９および第２の
適応フィルタ１１０に供給する。最後に、フィルタ係数
更新を行なうかどうかを決定して、係数更新制御信号１
０６として第１の適応フィルタ１０９および第２の適応
フィルタ１１０に供給する。
【００４５】セレクタ１０７は、第１の受信信号１およ
び第２の受信信号２を入力とし、２つの受信信号１，２
のうち、受信信号選択情報１０４で指定された受信信号
を選択し、第１の適応フィルタ１０９および第２の適応
フィルタ１１０に供給する。
び第２の受信信号２を入力とし、２つの受信信号１，２
のうち、受信信号選択情報１０４で指定された受信信号
を選択し、第１の適応フィルタ１０９および第２の適応
フィルタ１１０に供給する。
【００４６】第１の適応フィルタ１０９は、セレクタ１
０７が選択した受信信号１０８を入力として、第１の混
在信号１４に含まれるエコーに対応するエコーレプリカ
１１１を生成し、第１の減算器１１３に供給する。第１
の減算器１１３は、第１の混在信号１４から、第１の適
応フィルタ１０９の出力であるエコーレプリカ１１１を
差し引き、その結果を第１の出力信号１６とする。第１
の適応フィルタ１０９は、第１の出力信号１６を最小と
するように制御される。
０７が選択した受信信号１０８を入力として、第１の混
在信号１４に含まれるエコーに対応するエコーレプリカ
１１１を生成し、第１の減算器１１３に供給する。第１
の減算器１１３は、第１の混在信号１４から、第１の適
応フィルタ１０９の出力であるエコーレプリカ１１１を
差し引き、その結果を第１の出力信号１６とする。第１
の適応フィルタ１０９は、第１の出力信号１６を最小と
するように制御される。
【００４７】同様に、第２の適応フィルタ１１０は、セ
レクタ１０７が選択した受信信号１０８を入力として、
第２の混在信号１５に含まれるエコーに対応するエコー
レプリカ１１２を生成し、第２の減算器１１４に供給す
る。第２の減算器１１４は、第２の混在信号１５から、
第２の適応フィルタ１１０の出力であるエコーレプリカ
１１２を差し引き、その結果を第２の出力信号１７とす
る。第２の適応フィルタ１１０は、第２の出力信号１７
を最小とするように制御される。
レクタ１０７が選択した受信信号１０８を入力として、
第２の混在信号１５に含まれるエコーに対応するエコー
レプリカ１１２を生成し、第２の減算器１１４に供給す
る。第２の減算器１１４は、第２の混在信号１５から、
第２の適応フィルタ１１０の出力であるエコーレプリカ
１１２を差し引き、その結果を第２の出力信号１７とす
る。第２の適応フィルタ１１０は、第２の出力信号１７
を最小とするように制御される。
【００４８】第１の適応フィルタ１０９および第２の適
応フィルタ１１０の制御方法としては、Ｐｒｏｃｅｅｄ
ｉｎｇｓ ｏｆ ｔｈｅ ＩＥＥＥ，Ｖｏｌ．６３，Ｎ
ｏ．１２，ｐｐ．１６９２−１７１６，１９７５，ＵＳ
Ａ（以下文献４）にＬＭＳアルゴリズムが、ＩＥＥＥ
ｔｒａｎｓａｃｔｉｏｎｓ ｏｎ ａｕｔｏｍａｔｉｃ
ｃｏｎｔｒｏｌ，Ｖｏｌ．ＡＣ−１２，Ｎｏ．３，１
９６７，ＵＳＡ（以下文献５）に学習同定法が記載され
ている。例えば、ＬＭＳアルゴリズムに基づくトランス
バーサル形適応フィルタを仮定し、適応フィルタの動作
を説明する。
応フィルタ１１０の制御方法としては、Ｐｒｏｃｅｅｄ
ｉｎｇｓ ｏｆ ｔｈｅ ＩＥＥＥ，Ｖｏｌ．６３，Ｎ
ｏ．１２，ｐｐ．１６９２−１７１６，１９７５，ＵＳ
Ａ（以下文献４）にＬＭＳアルゴリズムが、ＩＥＥＥ
ｔｒａｎｓａｃｔｉｏｎｓ ｏｎ ａｕｔｏｍａｔｉｃ
ｃｏｎｔｒｏｌ，Ｖｏｌ．ＡＣ−１２，Ｎｏ．３，１
９６７，ＵＳＡ（以下文献５）に学習同定法が記載され
ている。例えば、ＬＭＳアルゴリズムに基づくトランス
バーサル形適応フィルタを仮定し、適応フィルタの動作
を説明する。
【００４９】第１の適応フィルタ１０９および第２の適
応フィルタ１１０のタップ数をＮ、時刻ｎにおける第１
の受信信号１と第２の受信信号２の間の遅延時間差をｎ
d 、第１の適応フィルタ１０９および第２の適応フィル
タ１１０の入力信号であるセレクタ１０７で選択された
受信信号１０８をｒ(n) 、第１の適応フィルタ１０９へ
の誤差信号であるエコー除去装置１００の第１の出力信
号１６をｅ1(n)、第１の適応フィルタ１０９の複数組の
フィルタ係数のうち、
応フィルタ１１０のタップ数をＮ、時刻ｎにおける第１
の受信信号１と第２の受信信号２の間の遅延時間差をｎ
d 、第１の適応フィルタ１０９および第２の適応フィル
タ１１０の入力信号であるセレクタ１０７で選択された
受信信号１０８をｒ(n) 、第１の適応フィルタ１０９へ
の誤差信号であるエコー除去装置１００の第１の出力信
号１６をｅ1(n)、第１の適応フィルタ１０９の複数組の
フィルタ係数のうち、
【数１１】 で与えられる。
【００５０】図２に、受信信号がＭチャンネルである場
合の遅延時間差推定回路の一具体例を示す。遅延時間差
推定回路１０１は、Ｋ個（Ｋ≧１）の２信号間遅延時間
差推定回路２１０1,…，２１０K および制御回路２０５
から構成され、複数の受信信号２０１1,…，２０１M を
入力とし、複数の受信信号２０１1,…，２０１M から任
意の２信号を選択する複数の組合せに対応する複数の遅
延時間差の組合せからなる遅延時間差情報２０２を出力
する。Ｋ個の２信号間遅延時間差推定回路２１０1,…，
２１０K は全て同じ構成で、かつ、同様の動作をするの
で、以下、個々の２信号間遅延時間差推定回路について
説明する時には、単に添字ｉを省略して、２信号間遅延
時間差推定回路２１０、制御信号２０４、２信号間遅延
時間差２０３と記述する。
合の遅延時間差推定回路の一具体例を示す。遅延時間差
推定回路１０１は、Ｋ個（Ｋ≧１）の２信号間遅延時間
差推定回路２１０1,…，２１０K および制御回路２０５
から構成され、複数の受信信号２０１1,…，２０１M を
入力とし、複数の受信信号２０１1,…，２０１M から任
意の２信号を選択する複数の組合せに対応する複数の遅
延時間差の組合せからなる遅延時間差情報２０２を出力
する。Ｋ個の２信号間遅延時間差推定回路２１０1,…，
２１０K は全て同じ構成で、かつ、同様の動作をするの
で、以下、個々の２信号間遅延時間差推定回路について
説明する時には、単に添字ｉを省略して、２信号間遅延
時間差推定回路２１０、制御信号２０４、２信号間遅延
時間差２０３と記述する。
【００５１】２信号間遅延時間差推定回路２１０は、Ｍ
チャンネルの受信信号２０１1,…，２０１M を入力と
し、これらの受信信号２０１1,…，２０１M のうち、制
御回路２０５から送られる制御信号２０４によって指定
された２信号の間の遅延時間差２０３を推定し、制御装
置２０５に供給する。
チャンネルの受信信号２０１1,…，２０１M を入力と
し、これらの受信信号２０１1,…，２０１M のうち、制
御回路２０５から送られる制御信号２０４によって指定
された２信号の間の遅延時間差２０３を推定し、制御装
置２０５に供給する。
【００５２】制御装置２０５は、Ｋ個の２信号間遅延時
間差推定回路２１０1,…，２１０Kの出力する２信号間
遅延時間差２０３1,…，２０３K を入力とし、２信号間
遅延時間差を推定する受信信号を指定する制御信号２０
４1,…，２０４K を２信号間遅延時間差推定回路２１０
1,…，２１０K の各々に供給する。複数の受信信号２０
１1,…，２０１M から任意の２信号を選択する複数の組
合せに対応する２信号間遅延時間差が全て揃えて、遅延
時間差情報２０２として出力する。
間差推定回路２１０1,…，２１０Kの出力する２信号間
遅延時間差２０３1,…，２０３K を入力とし、２信号間
遅延時間差を推定する受信信号を指定する制御信号２０
４1,…，２０４K を２信号間遅延時間差推定回路２１０
1,…，２１０K の各々に供給する。複数の受信信号２０
１1,…，２０１M から任意の２信号を選択する複数の組
合せに対応する２信号間遅延時間差が全て揃えて、遅延
時間差情報２０２として出力する。
【００５３】遅延時間差推定回路１０１において、２信
号間遅延時間差推定回路２１０1,…，２１０K は、１つ
の回路を繰り返し使用してもよいし、複数の回路を繰り
返し使用してもよい。受信信号がＭチャンネルの時に、
Ｍ（Ｍ−１）／２個の２信号間遅延時間差推定回路を用
意し、Ｍチャンネルの受信信号から２チャンネルを選ぶ
全ての組合せについて、同時に２信号間遅延時間差を推
定すると、各２信号間遅延時間差推定回路を繰り返し使
用する必要がなく、推定に要する時間を短くすることが
できる。Ｍチャンネルの受信信号から任意に選んだ第１
の受信信号と第２の受信信号の間の遅延時間差をｔ12、
第２の受信信号と第３の受信信号の間の遅延時間差をｔ
23とすると、第１の受信信号と第３の受信信号の間の遅
延時間差ｔ13は ｔ13＝ｔ12＋ｔ23 （３１） で求めることができるので、Ｍ−１個の２信号間遅延時
間差に基づいて、任意の２信号間の遅延時間差を求める
ことができる。なお、受信信号が２チャンネルの場合に
は、１つの２信号間遅延時間差推定回路を１回使用する
だけで、全ての遅延時間差を推定できる。
号間遅延時間差推定回路２１０1,…，２１０K は、１つ
の回路を繰り返し使用してもよいし、複数の回路を繰り
返し使用してもよい。受信信号がＭチャンネルの時に、
Ｍ（Ｍ−１）／２個の２信号間遅延時間差推定回路を用
意し、Ｍチャンネルの受信信号から２チャンネルを選ぶ
全ての組合せについて、同時に２信号間遅延時間差を推
定すると、各２信号間遅延時間差推定回路を繰り返し使
用する必要がなく、推定に要する時間を短くすることが
できる。Ｍチャンネルの受信信号から任意に選んだ第１
の受信信号と第２の受信信号の間の遅延時間差をｔ12、
第２の受信信号と第３の受信信号の間の遅延時間差をｔ
23とすると、第１の受信信号と第３の受信信号の間の遅
延時間差ｔ13は ｔ13＝ｔ12＋ｔ23 （３１） で求めることができるので、Ｍ−１個の２信号間遅延時
間差に基づいて、任意の２信号間の遅延時間差を求める
ことができる。なお、受信信号が２チャンネルの場合に
は、１つの２信号間遅延時間差推定回路を１回使用する
だけで、全ての遅延時間差を推定できる。
【００５４】図３に２信号間遅延時間差推定回路の第１
の具体例を示す。この２信号間遅延時間差推定回路２１
０は、複数の受信信号２０１1,…，２０１M を入力と
し、受信信号のうち、制御信号２０４で指定された２信
号間の遅延時間差２０３を推定する。
の具体例を示す。この２信号間遅延時間差推定回路２１
０は、複数の受信信号２０１1,…，２０１M を入力と
し、受信信号のうち、制御信号２０４で指定された２信
号間の遅延時間差２０３を推定する。
【００５５】２信号間遅延時間差の推定は、２信号を互
いに予測する２つの適応フィルタのフィルタ係数を用い
て行なう。遅延時間が短い受信信号を入力として、遅延
時間が長い受信信号を予測した適応フィルタでは、２信
号間遅延時間差と最も近い遅延を入力信号に与えるタッ
プのフィルタ係数絶対値が最大となる。逆に、遅延時間
が長い受信信号を入力として遅延時間が短い受信信号を
予測する場合は、未来の信号の予測は予測精度が低いの
で、フィルタ係数が十分に成長しない。したがって、フ
ィルタ係数は、遅延時間が長い受信信号を予測する場合
に比べて小さな値となる。
いに予測する２つの適応フィルタのフィルタ係数を用い
て行なう。遅延時間が短い受信信号を入力として、遅延
時間が長い受信信号を予測した適応フィルタでは、２信
号間遅延時間差と最も近い遅延を入力信号に与えるタッ
プのフィルタ係数絶対値が最大となる。逆に、遅延時間
が長い受信信号を入力として遅延時間が短い受信信号を
予測する場合は、未来の信号の予測は予測精度が低いの
で、フィルタ係数が十分に成長しない。したがって、フ
ィルタ係数は、遅延時間が長い受信信号を予測する場合
に比べて小さな値となる。
【００５６】この事実に基づき、２つの適応フィルタで
２信号を互いに予測し、適応フィルタの係数の絶対値が
最大値であるフィルタ係数を見つけることによって、第
１の受信信号２１３および第２の受信信号２１４の間の
遅延時間差を推定できる。
２信号を互いに予測し、適応フィルタの係数の絶対値が
最大値であるフィルタ係数を見つけることによって、第
１の受信信号２１３および第２の受信信号２１４の間の
遅延時間差を推定できる。
【００５７】遅延時間差を推定する対象となる受信信号
を指定する制御信号２０４で指定された２信号の一方
を、複数の受信信号２０１1,…，２０１M を入力とする
第１のセレクタ２１１で選択し、第１の受信信号２１３
として第１の適応フィルタ２３１および第２の減算器２
３８に供給する。前記制御信号２０４で指定された２信
号の他方を受信信号２０１1,…，２０１M を入力とする
第２のセレクタ２１２で選択し、第２の受信信号２１４
として第２の適応フィルタ２３２および第２の減算器２
３７に供給する。
を指定する制御信号２０４で指定された２信号の一方
を、複数の受信信号２０１1,…，２０１M を入力とする
第１のセレクタ２１１で選択し、第１の受信信号２１３
として第１の適応フィルタ２３１および第２の減算器２
３８に供給する。前記制御信号２０４で指定された２信
号の他方を受信信号２０１1,…，２０１M を入力とする
第２のセレクタ２１２で選択し、第２の受信信号２１４
として第２の適応フィルタ２３２および第２の減算器２
３７に供給する。
【００５８】第１の受信信号２１３を入力とする第１の
Ｌタップトランスバーサル形適応フィルタ２３１で、第
２の受信信号２１４を予測し、予測された信号を第１の
減算器２３７に供給する。第１の減算器２３７は、第２
の受信信号２１４から第１のトランスバーサル形適応フ
ィルタ２３１の出力を差し引く。第１のトランスバーサ
ル形適応フィルタ２３１は、第１の減算器２３７の出力
を最小とするように制御される。
Ｌタップトランスバーサル形適応フィルタ２３１で、第
２の受信信号２１４を予測し、予測された信号を第１の
減算器２３７に供給する。第１の減算器２３７は、第２
の受信信号２１４から第１のトランスバーサル形適応フ
ィルタ２３１の出力を差し引く。第１のトランスバーサ
ル形適応フィルタ２３１は、第１の減算器２３７の出力
を最小とするように制御される。
【００５９】第２の受信信号２１４を入力とする第２の
Ｌタップトランスバーサル形適応フィルタ２３２で、第
１の受信信号２１３を予測し、予測された信号を第２の
減算器２３８に供給する。第２の減算器２３８は、第１
の受信信号２１３から第２のトランスバーサル形適応フ
ィルタ２３２の出力を差し引く。第２のトランスバーサ
ル形適応フィルタ２３２は、第２の減算器２３８の出力
を最小とするように制御される。
Ｌタップトランスバーサル形適応フィルタ２３２で、第
１の受信信号２１３を予測し、予測された信号を第２の
減算器２３８に供給する。第２の減算器２３８は、第１
の受信信号２１３から第２のトランスバーサル形適応フ
ィルタ２３２の出力を差し引く。第２のトランスバーサ
ル形適応フィルタ２３２は、第２の減算器２３８の出力
を最小とするように制御される。
【００６０】Ｌ個の絶対値計算回路２３５1,…，２３５
L からなる第１の絶対値計算回路群には第１の適応フィ
ルタ２３１の係数２３３1,…，２３３L が、Ｌ個の絶対
値計算回路２３６1,…，２３６L からなる第２の絶対値
計算回路群には第２の適応フィルタ２３２の係数２３４
1,…，２３４L が供給されており、求められた各々の絶
対値は判定器２３９に供給される。判定器２３９は、絶
対値が最大となるフィルタ係数を求め、どのフィルタ係
数の絶対値が最大であるかに基づいて、２信号間遅延時
間差２０３を推定する。
L からなる第１の絶対値計算回路群には第１の適応フィ
ルタ２３１の係数２３３1,…，２３３L が、Ｌ個の絶対
値計算回路２３６1,…，２３６L からなる第２の絶対値
計算回路群には第２の適応フィルタ２３２の係数２３４
1,…，２３４L が供給されており、求められた各々の絶
対値は判定器２３９に供給される。判定器２３９は、絶
対値が最大となるフィルタ係数を求め、どのフィルタ係
数の絶対値が最大であるかに基づいて、２信号間遅延時
間差２０３を推定する。
【００６１】図４に、制御回路１０３の第１の具体例を
示す。受信信号選択回路４０１は、遅延時間差１０２に
基づいて、複数の受信信号のうち前記遅延時間が最も短
いものを検出して受信信号選択情報１０４とし、フィル
タ係数選択回路４０２に供給するとともに、制御回路の
第１の出力とする。
示す。受信信号選択回路４０１は、遅延時間差１０２に
基づいて、複数の受信信号のうち前記遅延時間が最も短
いものを検出して受信信号選択情報１０４とし、フィル
タ係数選択回路４０２に供給するとともに、制御回路の
第１の出力とする。
【００６２】受信信号が２チャンネルの場合は、遅延時
間差１０２は２信号間遅延時間差となるから、この遅延
時間差の符号に基づいて２信号のどちらが遅延時間が短
いかを判定することができる。受信信号が３チャンネル
以上の時は、遅延時間差１０２は複数の２信号間遅延時
間差となる。この場合は、複数の受信信号から２信号を
選び、この２信号に対応する２信号間遅延時間差の符号
に基づいて２信号のどちらの遅延時間が長いかを判定
し、遅延時間が長いと判定された受信信号を判定対象か
ら除外する操作を、遅延時間が最も短い受信信号だけが
残るまで繰り返せばよい。
間差１０２は２信号間遅延時間差となるから、この遅延
時間差の符号に基づいて２信号のどちらが遅延時間が短
いかを判定することができる。受信信号が３チャンネル
以上の時は、遅延時間差１０２は複数の２信号間遅延時
間差となる。この場合は、複数の受信信号から２信号を
選び、この２信号に対応する２信号間遅延時間差の符号
に基づいて２信号のどちらの遅延時間が長いかを判定
し、遅延時間が長いと判定された受信信号を判定対象か
ら除外する操作を、遅延時間が最も短い受信信号だけが
残るまで繰り返せばよい。
【００６３】フィルタ係数選択回路４０２は、前記遅延
時間差１０２および前記フィルタ係数選択情報１０４に
基づいて複数の適応フィルタが使用するフィルタ係数の
組を選択し、制御回路の第２の出力であるフィルタ係数
選択信号１０５とする。
時間差１０２および前記フィルタ係数選択情報１０４に
基づいて複数の適応フィルタが使用するフィルタ係数の
組を選択し、制御回路の第２の出力であるフィルタ係数
選択信号１０５とする。
【００６４】受信信号がＭチャンネルの場合には、以下
の方法でフィルタ係数の組を選択できる。ここでは、Ｍ
個の適応フィルタ全てが同じフィルタ係数の組番号を使
用すると仮定するが、これによって何ら一般性は失われ
ない。前記遅延時間差１０２および前記フィルタ係数選
択情報１０４に基づいて、各受信信号が、遅延時間が最
も短い受信信号からどれだけ遅れているかを求める。こ
の遅れ時間をそれぞれｔ1,…，ｔM とすると、離散時間
処理においては、一般性を失うことなく、ｔ1,…，ｔM
は整数値０，１，…，ｔmax のいずれかを取るものと仮
定できる。０，１，…，ｔmax から重複を許してＭ個を
選ぶ順列は（ｔmax ＋１）M 通りあるから、（ｔmax ＋
１）M 組のフィルタ係数を用意し、
の方法でフィルタ係数の組を選択できる。ここでは、Ｍ
個の適応フィルタ全てが同じフィルタ係数の組番号を使
用すると仮定するが、これによって何ら一般性は失われ
ない。前記遅延時間差１０２および前記フィルタ係数選
択情報１０４に基づいて、各受信信号が、遅延時間が最
も短い受信信号からどれだけ遅れているかを求める。こ
の遅れ時間をそれぞれｔ1,…，ｔM とすると、離散時間
処理においては、一般性を失うことなく、ｔ1,…，ｔM
は整数値０，１，…，ｔmax のいずれかを取るものと仮
定できる。０，１，…，ｔmax から重複を許してＭ個を
選ぶ順列は（ｔmax ＋１）M 通りあるから、（ｔmax ＋
１）M 組のフィルタ係数を用意し、
【数１２】 番目のフィルタ係数の組を使用するようにすればよい。
【００６５】２チャンネルの場合には、前記遅延時間差
１０２は２信号間遅延時間差となるので、以下の方法で
フィルタ係数の組を選択してもよい。離散時間処理にお
いては、一般性を失うことなく、２信号間遅延時間差ｔ
は２ｔmax ＋１個の整数値−ｔmax,…，０，…，ｔmax
のいずれかを取るものと仮定できる。したがって、２ｔ
max ＋１組のフィルタ係数を用意し、ｔ＋ｔmax ＋１番
目のフィルタ係数の組を使用するようにすればよい。
１０２は２信号間遅延時間差となるので、以下の方法で
フィルタ係数の組を選択してもよい。離散時間処理にお
いては、一般性を失うことなく、２信号間遅延時間差ｔ
は２ｔmax ＋１個の整数値−ｔmax,…，０，…，ｔmax
のいずれかを取るものと仮定できる。したがって、２ｔ
max ＋１組のフィルタ係数を用意し、ｔ＋ｔmax ＋１番
目のフィルタ係数の組を使用するようにすればよい。
【００６６】係数更新制御回路４０３は、前記複数の適
応フィルタが常にフィルタ係数を更新するように、制御
回路の第３の出力である係数更新制御信号１０６を出力
する。この係数更新制御回路４０３は、単なる係数レジ
スタ、あるいは、パルス発生器で実現できる。
応フィルタが常にフィルタ係数を更新するように、制御
回路の第３の出力である係数更新制御信号１０６を出力
する。この係数更新制御回路４０３は、単なる係数レジ
スタ、あるいは、パルス発生器で実現できる。
【００６７】図５に、第１の適応フィルタ１０９および
第２の適応フィルタ１１０として用いる適応フィルタの
具体例を示す。この適応フィルタは、セレクタ１０７で
選択された受信信号１０８、フィルタ係数選択信号１０
５、係数更新制御信号１０６、エコー除去装置の出力信
号５０１を入力、エコーレプリカ５０２を出力とし、受
信信号１０８を格納する信号メモリ５０３、フィルタ係
数を格納するＪ個の係数メモリ５０４1,…，５０４J 、
フィルタ係数を選択するセレクタ５０５、畳み込み演算
やフィルタ係数更新などを行なう演算回路５０６から構
成される。
第２の適応フィルタ１１０として用いる適応フィルタの
具体例を示す。この適応フィルタは、セレクタ１０７で
選択された受信信号１０８、フィルタ係数選択信号１０
５、係数更新制御信号１０６、エコー除去装置の出力信
号５０１を入力、エコーレプリカ５０２を出力とし、受
信信号１０８を格納する信号メモリ５０３、フィルタ係
数を格納するＪ個の係数メモリ５０４1,…，５０４J 、
フィルタ係数を選択するセレクタ５０５、畳み込み演算
やフィルタ係数更新などを行なう演算回路５０６から構
成される。
【００６８】エコーレプリカ５０２の計算およびフィル
タ係数更新は、以下の手順で行なう。まず、演算回路５
０６は、前記受信信号１０８を入力し、信号メモリ５０
３に格納する。次に、セレクタ５０５は、複数のフィル
タ係数メモリ５０４1,…，５０４J のうち、フィルタ係
数選択信号１０５で指定されたものを選択し、演算回路
５０６は、信号メモリ５０３に格納された受信信号１０
８とセレクタ５０５が選択した係数メモリの内容との畳
み込み演算を行ない、エコーレプリカ５０１とする。最
後に、演算回路５０６は、セレクタ５０５が選択した係
数メモリの内容を読み出し、エコー除去装置の出力信号
５０１を最小とするように更新し、セレクタ５０５が選
択した係数メモリに書き込む。
タ係数更新は、以下の手順で行なう。まず、演算回路５
０６は、前記受信信号１０８を入力し、信号メモリ５０
３に格納する。次に、セレクタ５０５は、複数のフィル
タ係数メモリ５０４1,…，５０４J のうち、フィルタ係
数選択信号１０５で指定されたものを選択し、演算回路
５０６は、信号メモリ５０３に格納された受信信号１０
８とセレクタ５０５が選択した係数メモリの内容との畳
み込み演算を行ない、エコーレプリカ５０１とする。最
後に、演算回路５０６は、セレクタ５０５が選択した係
数メモリの内容を読み出し、エコー除去装置の出力信号
５０１を最小とするように更新し、セレクタ５０５が選
択した係数メモリに書き込む。
【００６９】図６に、演算回路の一具体例を示す。以下
の説明では、簡単のため、Ｊ個の係数メモリ５０４
1,…，５０４J のうちセレクタ５０５が選択したもの
を、単に係数メモリ５０４と略す。この演算回路は、信
号メモリ５０３のアドレスを指定する第１のアドレスバ
ス５１０、信号メモリ５０３とのデータ転送を行なう第
１のデータバス５１１、係数メモリ５０４のアドレスを
指定する第２のアドレスバス５１２、係数メモリ５０４
とのデータ転送を行なう第２のデータバス５１３、受信
信号１０８を入力する第１の入力端子５１４、係数更新
制御信号１０６を入力する第２の入力端子５１５、エコ
ー除去装置の出力信号５０１を入力する第３の入力端子
５１６、出力信号５０１を出力する出力端子５１７、第
１のアドレスバス５１０に信号メモリ５０３のアドレス
を供給する第１のアドレス生成回路５１８、第２のアド
レスバス５１２に係数メモリ５０４のアドレスを供給す
る第２のアドレス生成回路５１９、乗算器５２０、乗算
器５２０の入力データを選択する２個のセレクタ５２
１，５２２、加減算などを行なう算術論理演算回路であ
るＡＬＵ５２３、ＡＬＵ５２３の入力データを選択する
セレクタ５２４、Ｉ個のレジスタ５２５1,…，５２５I
で構成されるレジスタバンク５２６、制御プログラムを
格納するプログラムメモリ５２７、演算回路全体を制御
する制御回路５２８で構成される。制御プログラムで
は、算術論理演算命令、データ転送命令、分岐命令の他
に、制御プログラム中に記述されたタップ数やステップ
サイズなどの定数を制御回路５２８を通してレジスタバ
ンク５２６中の任意のレジスタや乗算器５２０などに転
送する定数設定命令を使用できる。
の説明では、簡単のため、Ｊ個の係数メモリ５０４
1,…，５０４J のうちセレクタ５０５が選択したもの
を、単に係数メモリ５０４と略す。この演算回路は、信
号メモリ５０３のアドレスを指定する第１のアドレスバ
ス５１０、信号メモリ５０３とのデータ転送を行なう第
１のデータバス５１１、係数メモリ５０４のアドレスを
指定する第２のアドレスバス５１２、係数メモリ５０４
とのデータ転送を行なう第２のデータバス５１３、受信
信号１０８を入力する第１の入力端子５１４、係数更新
制御信号１０６を入力する第２の入力端子５１５、エコ
ー除去装置の出力信号５０１を入力する第３の入力端子
５１６、出力信号５０１を出力する出力端子５１７、第
１のアドレスバス５１０に信号メモリ５０３のアドレス
を供給する第１のアドレス生成回路５１８、第２のアド
レスバス５１２に係数メモリ５０４のアドレスを供給す
る第２のアドレス生成回路５１９、乗算器５２０、乗算
器５２０の入力データを選択する２個のセレクタ５２
１，５２２、加減算などを行なう算術論理演算回路であ
るＡＬＵ５２３、ＡＬＵ５２３の入力データを選択する
セレクタ５２４、Ｉ個のレジスタ５２５1,…，５２５I
で構成されるレジスタバンク５２６、制御プログラムを
格納するプログラムメモリ５２７、演算回路全体を制御
する制御回路５２８で構成される。制御プログラムで
は、算術論理演算命令、データ転送命令、分岐命令の他
に、制御プログラム中に記述されたタップ数やステップ
サイズなどの定数を制御回路５２８を通してレジスタバ
ンク５２６中の任意のレジスタや乗算器５２０などに転
送する定数設定命令を使用できる。
【００７０】前記受信信号１０８を信号メモリ５０３に
格納する操作は、以下の手順で行なう。 １：受信信号１０８を第１の入力端子５１４から読み込
み、第１のレジスタ５２５1 に格納する。 ２：第１のアドレス生成回路５１８が受信信号１０８を
記憶する信号メモリ５０３のアドレスを第１のアドレス
バス５１０に供給し、第１のレジスタ５２５1 に格納さ
れている受信信号１０８を第１のデータバス５１１に供
給し、受信信号１０８を信号メモリ５０３に格納する。 畳み込み演算は、以下の手順で行なう。 １：定数設定命令を用いて、第１のレジスタ５２５1 を
０に初期化する。 ２：定数設定命令を用いて、第２のレジスタ５２５
2 に、適応フィルタのタップ数を格納する。 ３：第１のアドレス生成回路５１８を、信号メモリ５０
３の最初のアドレスを指示するように初期化する。 ４：第２のアドレス生成回路５１９を、係数メモリ５０
４の最初のアドレスを指示するように初期化する。 ５：第１のアドレス生成回路５１８が信号メモリ５０３
のアドレスを第１のアドレスバス５１０に供給し、第１
のデータバス５１１から信号メモリ５０３の値を読み込
み、第１のセレクタ５２１を通して乗算器５２０に供給
する。 ６：第２のアドレス生成回路５１９が係数メモリ５０４
のアドレスを第２のアドレスバス５１２に供給し、第２
のデータバス５１３から係数メモリ５０４の値を読み込
み、第２のセレクタ５２２を通して乗算器５２０に供給
する。 ７：乗算器５２０が、第１のセレクタ５２１および第２
のセレクタ５２２から供給されたデータを乗算する。 ８：乗算結果を、第３のセレクタ５２４を通してＡＬＵ
５２３に供給する。 ９：第１のレジスタ５２５1 の値をＡＬＵ５２３に供給
する。 10：ＡＬＵ５２３が、第３のセレクタ５２４および第１
のレジスタ５２５1 から供給されたデータを加算する。 11：加算結果を第１のレジスタ５２５1 に格納する。 12：第１のアドレス生成回路５１８が信号メモリ５０３
の次のアドレスを指示するように、アドレスを更新す
る。 13：第２のアドレス生成回路５１９が係数メモリ５０４
の次のアドレスを指示するように、アドレスを更新す
る。 14：第２のレジスタ５２５2 の値をＡＬＵ５２３に供給
する。 15：ＡＬＵ５２３は、第２のレジスタ５２５2 の値から
１を減ずる。 16：減算結果を第２のレジスタ５２５2 に格納する。 17：５から１６の操作を、第２のレジスタ５２５2 の値
が０になるまで繰り返す。 18：第１のレジスタ５２５1 の値を出力端子に供給し
て、出力信号５０１とする。
格納する操作は、以下の手順で行なう。 １：受信信号１０８を第１の入力端子５１４から読み込
み、第１のレジスタ５２５1 に格納する。 ２：第１のアドレス生成回路５１８が受信信号１０８を
記憶する信号メモリ５０３のアドレスを第１のアドレス
バス５１０に供給し、第１のレジスタ５２５1 に格納さ
れている受信信号１０８を第１のデータバス５１１に供
給し、受信信号１０８を信号メモリ５０３に格納する。 畳み込み演算は、以下の手順で行なう。 １：定数設定命令を用いて、第１のレジスタ５２５1 を
０に初期化する。 ２：定数設定命令を用いて、第２のレジスタ５２５
2 に、適応フィルタのタップ数を格納する。 ３：第１のアドレス生成回路５１８を、信号メモリ５０
３の最初のアドレスを指示するように初期化する。 ４：第２のアドレス生成回路５１９を、係数メモリ５０
４の最初のアドレスを指示するように初期化する。 ５：第１のアドレス生成回路５１８が信号メモリ５０３
のアドレスを第１のアドレスバス５１０に供給し、第１
のデータバス５１１から信号メモリ５０３の値を読み込
み、第１のセレクタ５２１を通して乗算器５２０に供給
する。 ６：第２のアドレス生成回路５１９が係数メモリ５０４
のアドレスを第２のアドレスバス５１２に供給し、第２
のデータバス５１３から係数メモリ５０４の値を読み込
み、第２のセレクタ５２２を通して乗算器５２０に供給
する。 ７：乗算器５２０が、第１のセレクタ５２１および第２
のセレクタ５２２から供給されたデータを乗算する。 ８：乗算結果を、第３のセレクタ５２４を通してＡＬＵ
５２３に供給する。 ９：第１のレジスタ５２５1 の値をＡＬＵ５２３に供給
する。 10：ＡＬＵ５２３が、第３のセレクタ５２４および第１
のレジスタ５２５1 から供給されたデータを加算する。 11：加算結果を第１のレジスタ５２５1 に格納する。 12：第１のアドレス生成回路５１８が信号メモリ５０３
の次のアドレスを指示するように、アドレスを更新す
る。 13：第２のアドレス生成回路５１９が係数メモリ５０４
の次のアドレスを指示するように、アドレスを更新す
る。 14：第２のレジスタ５２５2 の値をＡＬＵ５２３に供給
する。 15：ＡＬＵ５２３は、第２のレジスタ５２５2 の値から
１を減ずる。 16：減算結果を第２のレジスタ５２５2 に格納する。 17：５から１６の操作を、第２のレジスタ５２５2 の値
が０になるまで繰り返す。 18：第１のレジスタ５２５1 の値を出力端子に供給し
て、出力信号５０１とする。
【００７１】第２の入力端子５１５から読み込んだ係数
更新制御信号１０６が、フィルタ係数を更新するように
指示している時は、以下の手順でフィルタ係数の更新を
行なう。 １：定数設定命令を用いて、第１のレジスタ５２５
1 に、適応フィルタのタップ数を格納する。 ２：第１のアドレス生成回路５１８を、信号メモリ５０
３の最初のアドレスを指示するように初期化する。 ３：第２のアドレス生成回路５１９を、係数メモリ５０
４の最初のアドレスを指示するように初期化する。 ４：第３の入力端子５１６からエコー除去装置の出力信
号５０１を読み込み、第１のセレクタ５２１を通して乗
算器５２０に供給する。 ５：定数設定命令を用いて、プログラムメモリ５２７に
格納されているステップサイズを制御回路５２８および
第２のセレクタ５２２を通して乗算器５２０に供給す
る。 ６：乗算器５２０が、第１のセレクタ５２１および第２
のセレクタ５２２から供給されたデータを乗算する。 ７：乗算結果を第２のレジスタ５２５2 に格納する。 ８：第２のアドレス生成回路５１９が係数メモリのアド
レスを第２のアドレスバス５１２に供給し、第２のデー
タバス５１３から係数メモリ５０４の値を読み込み、第
３のレジスタ５２５3 に格納する。 ９：第１のアドレス生成回路５１８が信号メモリ５０３
のアドレスを第１のアドレスバス５１０に供給し、第１
のデータバス５１１から信号メモリ５０３の値を読み込
み、第１のセレクタ５２１を通して乗算器５２０に供給
する。 10：第２のレジスタ５２５2 の値を第２のセレクタ５２
２を通して乗算器５２０に供給する。 11：乗算器５２０が、第１のセレクタ５２１および第２
のセレクタ５２２から供給されたデータを乗算する。 12：乗算結果を、第３のセレクタ５２４を通してＡＬＵ
５２３に供給する。 13：第３のレジスタ５２５3 の値をＡＬＵ５２３に供給
する。 14：ＡＬＵ５２３が、第３のセレクタ５２４および第３
のレジスタ５２５3 から供給されたデータを加算する。 15：加算結果を第３のレジスタ５２５3 に格納する。 16：第２のアドレス生成回路５１９が係数メモリ５０４
のアドレスを第２のアドレスバス５１２に供給し、前記
第３のレジスタ５２５3 の値を第２のデータバス５１３
に供給し、係数メモリ５０３に格納する。 17：第１のアドレス生成回路５１８が信号メモリ５０３
の次のアドレスを指示するように、アドレスを更新す
る。 18：第２のアドレス生成回路５１９が係数メモリ５０４
の次のアドレスを指示するように、アドレスを更新す
る。 19：第３のレジスタ５２５3 の値をＡＬＵ５２３に供給
する。 20：ＡＬＵ５２３は、第１のレジスタ５２５1 の値から
１を減ずる。 21：減算結果を第１のレジスタ５２５1 に格納する。 22：８から２１の操作を、第１のレジスタ５２５1 の値
が０になるまで繰り返す。
更新制御信号１０６が、フィルタ係数を更新するように
指示している時は、以下の手順でフィルタ係数の更新を
行なう。 １：定数設定命令を用いて、第１のレジスタ５２５
1 に、適応フィルタのタップ数を格納する。 ２：第１のアドレス生成回路５１８を、信号メモリ５０
３の最初のアドレスを指示するように初期化する。 ３：第２のアドレス生成回路５１９を、係数メモリ５０
４の最初のアドレスを指示するように初期化する。 ４：第３の入力端子５１６からエコー除去装置の出力信
号５０１を読み込み、第１のセレクタ５２１を通して乗
算器５２０に供給する。 ５：定数設定命令を用いて、プログラムメモリ５２７に
格納されているステップサイズを制御回路５２８および
第２のセレクタ５２２を通して乗算器５２０に供給す
る。 ６：乗算器５２０が、第１のセレクタ５２１および第２
のセレクタ５２２から供給されたデータを乗算する。 ７：乗算結果を第２のレジスタ５２５2 に格納する。 ８：第２のアドレス生成回路５１９が係数メモリのアド
レスを第２のアドレスバス５１２に供給し、第２のデー
タバス５１３から係数メモリ５０４の値を読み込み、第
３のレジスタ５２５3 に格納する。 ９：第１のアドレス生成回路５１８が信号メモリ５０３
のアドレスを第１のアドレスバス５１０に供給し、第１
のデータバス５１１から信号メモリ５０３の値を読み込
み、第１のセレクタ５２１を通して乗算器５２０に供給
する。 10：第２のレジスタ５２５2 の値を第２のセレクタ５２
２を通して乗算器５２０に供給する。 11：乗算器５２０が、第１のセレクタ５２１および第２
のセレクタ５２２から供給されたデータを乗算する。 12：乗算結果を、第３のセレクタ５２４を通してＡＬＵ
５２３に供給する。 13：第３のレジスタ５２５3 の値をＡＬＵ５２３に供給
する。 14：ＡＬＵ５２３が、第３のセレクタ５２４および第３
のレジスタ５２５3 から供給されたデータを加算する。 15：加算結果を第３のレジスタ５２５3 に格納する。 16：第２のアドレス生成回路５１９が係数メモリ５０４
のアドレスを第２のアドレスバス５１２に供給し、前記
第３のレジスタ５２５3 の値を第２のデータバス５１３
に供給し、係数メモリ５０３に格納する。 17：第１のアドレス生成回路５１８が信号メモリ５０３
の次のアドレスを指示するように、アドレスを更新す
る。 18：第２のアドレス生成回路５１９が係数メモリ５０４
の次のアドレスを指示するように、アドレスを更新す
る。 19：第３のレジスタ５２５3 の値をＡＬＵ５２３に供給
する。 20：ＡＬＵ５２３は、第１のレジスタ５２５1 の値から
１を減ずる。 21：減算結果を第１のレジスタ５２５1 に格納する。 22：８から２１の操作を、第１のレジスタ５２５1 の値
が０になるまで繰り返す。
【００７２】以上、図１から図６を用いて、本発明の実
施例の１つを説明した。図３に示した２信号間遅延時間
差推定回路２１０の第１の具体例では、２つの適応フィ
ルタで２信号を互いに予測し、適応フィルタの係数の絶
対値が最大値であるフィルタ係数を見つけることによっ
て、２信号間の遅延時間差を推定していた。本発明の他
の実施例として、２信号間の相互相関関数の絶対値を最
大とする時間差を見つけることによって、２信号間の遅
延時間差を推定することもできる。
施例の１つを説明した。図３に示した２信号間遅延時間
差推定回路２１０の第１の具体例では、２つの適応フィ
ルタで２信号を互いに予測し、適応フィルタの係数の絶
対値が最大値であるフィルタ係数を見つけることによっ
て、２信号間の遅延時間差を推定していた。本発明の他
の実施例として、２信号間の相互相関関数の絶対値を最
大とする時間差を見つけることによって、２信号間の遅
延時間差を推定することもできる。
【００７３】図７に２信号間遅延時間差推定回路２１０
の第２の具体例を示す。この２信号間遅延時間差推定回
路２１０は、複数の受信信号２０１1,…，２０１M を入
力とし、これらの受信信号のうち、制御信号２０４で指
定された２信号間の遅延時間差を推定する。
の第２の具体例を示す。この２信号間遅延時間差推定回
路２１０は、複数の受信信号２０１1,…，２０１M を入
力とし、これらの受信信号のうち、制御信号２０４で指
定された２信号間の遅延時間差を推定する。
【００７４】まず、２信号の相互相関関数を用いた遅延
時間差の推定原理を説明する。第１の受信信号２１３お
よび第２の受信信号２１４の時刻ｎにおける値をｘ1
(n), ｘ2(n)とすると、時刻ｎにおける時間差ｍに対す
る２信号の相互相関関数Ｒ12(n,m) は、 Ｒ12(n,m) ＝Ｅ［ｘ1(n)ｘ2(n+m)］ （３２） で与えられる。
時間差の推定原理を説明する。第１の受信信号２１３お
よび第２の受信信号２１４の時刻ｎにおける値をｘ1
(n), ｘ2(n)とすると、時刻ｎにおける時間差ｍに対す
る２信号の相互相関関数Ｒ12(n,m) は、 Ｒ12(n,m) ＝Ｅ［ｘ1(n)ｘ2(n+m)］ （３２） で与えられる。
【００７５】定常な受信信号ｘ1(n)と、それをｎd サン
プルだけ遅延させた信号ｘ2(n)＝ｘ１(n-nd ) を仮定す
ると、
プルだけ遅延させた信号ｘ2(n)＝ｘ１(n-nd ) を仮定す
ると、
【数１３】 となる。したがって、ｍ＝ｎd のとき、Ｒ12(n,m) は最
小値となり、絶対値は最大となる。
小値となり、絶対値は最大となる。
【００７６】この事実より、相互相関関数Ｒ12(n,m) の
絶対値が最大となる時間差ｍで第１の受信信号２１３お
よび第２の受信信号２１４の間の遅延時間差２０３を推
定できることがわかる。
絶対値が最大となる時間差ｍで第１の受信信号２１３お
よび第２の受信信号２１４の間の遅延時間差２０３を推
定できることがわかる。
【００７７】遅延時間差を求める２信号の選択は、２つ
のセレクタ２１１，２１２で行なう。遅延時間差を求め
る対象となる受信信号を指定する制御信号２０４で指定
された２信号の一方を、複数の受信信号２０１1,…，２
０１M を入力とする第１のセレクタ２１１で選択し、第
１の受信信号２１３とする。同様に、前記制御信号２０
４で指定された２信号の他方を受信信号２０１1,…，２
０１M を入力とする第２のセレクタ２１２で選択し、第
２の受信信号２１４とする。
のセレクタ２１１，２１２で行なう。遅延時間差を求め
る対象となる受信信号を指定する制御信号２０４で指定
された２信号の一方を、複数の受信信号２０１1,…，２
０１M を入力とする第１のセレクタ２１１で選択し、第
１の受信信号２１３とする。同様に、前記制御信号２０
４で指定された２信号の他方を受信信号２０１1,…，２
０１M を入力とする第２のセレクタ２１２で選択し、第
２の受信信号２１４とする。
【００７８】相互相関関数計算回路２５１は、前記第１
の受信信号２１３と第２の受信信号２１４の相互相関関
数のうち、予め定められたＪ個の時間差に対応する値２
５２1,…，２５２J を計算する。前記相互相関関数の値
２５２1,…，２５２J と１対１に対応するＪ個の絶対値
計算回路２５３1,…，２５３J からなる絶対値計算回路
群で、前記相互相関関数２５２1,…，２５２J の絶対値
２５４1,…，２５４Jを計算する。
の受信信号２１３と第２の受信信号２１４の相互相関関
数のうち、予め定められたＪ個の時間差に対応する値２
５２1,…，２５２J を計算する。前記相互相関関数の値
２５２1,…，２５２J と１対１に対応するＪ個の絶対値
計算回路２５３1,…，２５３J からなる絶対値計算回路
群で、前記相互相関関数２５２1,…，２５２J の絶対値
２５４1,…，２５４Jを計算する。
【００７９】判定器２５５は、前記相互相関関数２５２
1,…，２５２J の絶対値２５４1,…，２５４J を最大と
する時間差を求め、結果を２信号間遅延時間差２０３と
する。
1,…，２５２J の絶対値２５４1,…，２５４J を最大と
する時間差を求め、結果を２信号間遅延時間差２０３と
する。
【００８０】図８に相互相関関数計算回路２５１の実施
例を示す。この相互相関関数計算回路２５１は、時間差
ｍ＝−Ｌ，…,０，…，＋Ｌに対応するＪ＝２Ｌ＋１個
の相互相関関数値２５２1,…，２５２2L+1を計算する。
相互相関関数のうち、時間差ｍの符号が正の部分の計算
回路は、Ｌ個の遅延器で構成される第１のタップドディ
レイライン２１５、Ｌ個の乗算器２１９1,…，２１
９L 、Ｌ個の積分器２２１1,…，２２１L から構成され
る。第１のタップドディレイライン２１５で、第１の受
信信号２１３を１サンプル周期ずつ遅延させ、Ｌ個の乗
算器２１９1,…，２１９L からなる第１の乗算器群で、
第１のタップドディレイライン２１５の各タップ出力２
１７1,…，２１７L と第２の受信信号２１４を互いに乗
算する。乗算器２１９1,…，２１９L の出力は、Ｌ個の
積分器２２１1,…，２２１L からなる第１の積分器群で
積分される。積分結果が、相互相関関数Ｒ12(n,m) のう
ち、ｍ＝１，２，…，Ｌに対応する値となる。
例を示す。この相互相関関数計算回路２５１は、時間差
ｍ＝−Ｌ，…,０，…，＋Ｌに対応するＪ＝２Ｌ＋１個
の相互相関関数値２５２1,…，２５２2L+1を計算する。
相互相関関数のうち、時間差ｍの符号が正の部分の計算
回路は、Ｌ個の遅延器で構成される第１のタップドディ
レイライン２１５、Ｌ個の乗算器２１９1,…，２１
９L 、Ｌ個の積分器２２１1,…，２２１L から構成され
る。第１のタップドディレイライン２１５で、第１の受
信信号２１３を１サンプル周期ずつ遅延させ、Ｌ個の乗
算器２１９1,…，２１９L からなる第１の乗算器群で、
第１のタップドディレイライン２１５の各タップ出力２
１７1,…，２１７L と第２の受信信号２１４を互いに乗
算する。乗算器２１９1,…，２１９L の出力は、Ｌ個の
積分器２２１1,…，２２１L からなる第１の積分器群で
積分される。積分結果が、相互相関関数Ｒ12(n,m) のう
ち、ｍ＝１，２，…，Ｌに対応する値となる。
【００８１】相互相関関数のうち、時間差ｍの符合が負
の部分の計算回路は、Ｌ個の遅延器で構成される第２の
タップドディレイライン２１６、Ｌ個の乗算器２２０1,
…，２２０L 、Ｌ個の積分器２２２1,…，２２２L から
構成される。第２のタップドディレイライン２１６で、
第２の受信信号２１４を遅延させ、Ｌ個の乗算器２２０
1,…，２２０L からなる第２の乗算器群で、第２のタッ
プドディレイライン２１６の各タップ出力２１８1,…，
２１８L と第１の受信信号２１３を互いに乗算する。乗
算器２２０1,…，２２０L の出力は、Ｌ個の積分器２２
２1,…，２２２L からなる第２の積分器群で積分され
る。積分結果が、相互相関関数Ｒ12(n,m)のうち、ｍ＝
−１，−２，…，−Ｌに対応する値となる。
の部分の計算回路は、Ｌ個の遅延器で構成される第２の
タップドディレイライン２１６、Ｌ個の乗算器２２０1,
…，２２０L 、Ｌ個の積分器２２２1,…，２２２L から
構成される。第２のタップドディレイライン２１６で、
第２の受信信号２１４を遅延させ、Ｌ個の乗算器２２０
1,…，２２０L からなる第２の乗算器群で、第２のタッ
プドディレイライン２１６の各タップ出力２１８1,…，
２１８L と第１の受信信号２１３を互いに乗算する。乗
算器２２０1,…，２２０L の出力は、Ｌ個の積分器２２
２1,…，２２２L からなる第２の積分器群で積分され
る。積分結果が、相互相関関数Ｒ12(n,m)のうち、ｍ＝
−１，−２，…，−Ｌに対応する値となる。
【００８２】相互相関関数のうち、時間差ｍが零の部分
の計算回路は、乗算器２２３、積分器２２４から構成さ
れる。乗算器２２３で第１の受信信号２１３および第２
の受信信号２１４を互いに乗算する。積分器２２４で乗
算器２２３の出力を積分したものが、相互相関関数のう
ちのＲ12(n,0) となる。
の計算回路は、乗算器２２３、積分器２２４から構成さ
れる。乗算器２２３で第１の受信信号２１３および第２
の受信信号２１４を互いに乗算する。積分器２２４で乗
算器２２３の出力を積分したものが、相互相関関数のう
ちのＲ12(n,0) となる。
【００８３】図９に、積分器２２１1,…，２２１L ，２
２２1,…，２２２L ，２２６の第１の具体例であるトラ
ンスバーサル型積分器を示す。この積分器は、入力信号
３０１の積分結果を出力信号３０２とし、タップドディ
レイライン３２１、Ｎ個の係数乗算器３２２1,…，３２
２N 、加算器３２３から構成される。入力信号３０１
は、タップドディレイライン３２１に格納される。タッ
プドディレイライン３２１の各タップ出力と１対１に対
応するＮ個の係数乗算器３２２1,…，３２２N で、タッ
プドディレイライン３２１の各タップ出力に各タップ毎
に予め定められた定数を乗ずる。加算器３２３で、Ｎ個
の係数乗算器３２２1,…，３２２N の各々の乗算結果の
総和を求め、積分器の出力３０２とする。
２２1,…，２２２L ，２２６の第１の具体例であるトラ
ンスバーサル型積分器を示す。この積分器は、入力信号
３０１の積分結果を出力信号３０２とし、タップドディ
レイライン３２１、Ｎ個の係数乗算器３２２1,…，３２
２N 、加算器３２３から構成される。入力信号３０１
は、タップドディレイライン３２１に格納される。タッ
プドディレイライン３２１の各タップ出力と１対１に対
応するＮ個の係数乗算器３２２1,…，３２２N で、タッ
プドディレイライン３２１の各タップ出力に各タップ毎
に予め定められた定数を乗ずる。加算器３２３で、Ｎ個
の係数乗算器３２２1,…，３２２N の各々の乗算結果の
総和を求め、積分器の出力３０２とする。
【００８４】図１０に、積分器２２１1,…，２２１L ，
２２２1,…，２２２L ，２２６の第２の具体例である単
純平均計算回路を示す。この積分器は、入力信号３０１
の積分結果を出力信号３０２とし、Ｄサンプルの遅延器
３０３、遅延器３０４、加算器３０５、係数乗算器３０
６から構成される。入力信号３０１は、加算器３０５に
供給されるとともに、遅延器３０３に格納される。遅延
器３０３は、入力信号３０１をＤサンプル遅延させ、加
算器３０５に供給する。加算器３０５で、遅延器３０４
の値と入力信号３０１を加算し、遅延器３０３の出力を
差し引く。加算器３０５の出力は、遅延器３０４に格納
されると同時に、係数乗算器３０６に供給される。係数
乗算器３０６で加算器３０５の出力をａ倍したものが、
積分器の出力３０２となる。ａは任意の正数であるが、
ａ＝１／Ｄとしたときには、積分結果が入力信号の単純
平均となる。
２２２1,…，２２２L ，２２６の第２の具体例である単
純平均計算回路を示す。この積分器は、入力信号３０１
の積分結果を出力信号３０２とし、Ｄサンプルの遅延器
３０３、遅延器３０４、加算器３０５、係数乗算器３０
６から構成される。入力信号３０１は、加算器３０５に
供給されるとともに、遅延器３０３に格納される。遅延
器３０３は、入力信号３０１をＤサンプル遅延させ、加
算器３０５に供給する。加算器３０５で、遅延器３０４
の値と入力信号３０１を加算し、遅延器３０３の出力を
差し引く。加算器３０５の出力は、遅延器３０４に格納
されると同時に、係数乗算器３０６に供給される。係数
乗算器３０６で加算器３０５の出力をａ倍したものが、
積分器の出力３０２となる。ａは任意の正数であるが、
ａ＝１／Ｄとしたときには、積分結果が入力信号の単純
平均となる。
【００８５】本実施例では、積分器の入力信号３０１を
遅延器３０３および加算器３０５に供給し、係数乗算器
３０６で加算器３０５の出力をａ倍し、乗算結果を積分
器の出力信号３０２としているが、係数乗算器３０６の
位置を変えて、積分器の入力信号３０１を係数乗算器３
０６に供給してａ倍したものを遅延器３０３および加算
器３０５に供給し、加算器３０５の出力を積分器の出力
信号３０２としてもよい。
遅延器３０３および加算器３０５に供給し、係数乗算器
３０６で加算器３０５の出力をａ倍し、乗算結果を積分
器の出力信号３０２としているが、係数乗算器３０６の
位置を変えて、積分器の入力信号３０１を係数乗算器３
０６に供給してａ倍したものを遅延器３０３および加算
器３０５に供給し、加算器３０５の出力を積分器の出力
信号３０２としてもよい。
【００８６】図１１に、積分器２２１1,…，２２１L ，
２２２1,…，２２２L ，２２６の第３の具体例である１
次の再帰形積分器を示す。この積分器は、第１の係数乗
算器３１１、加算器３１２、遅延器３１３、第２の係数
乗算器３１４から構成され、入力信号３０１を積分し
て、結果を出力信号３０２とする。第１の係数乗算器３
１１は積分器の入力信号３０１をα倍し、乗算結果を加
算器３１２に供給する。加算器３１２は第１および第２
の係数乗算器３１１，３１４の乗算結果を加算し、加算
結果を積分器の出力信号３０２として出力するととも
に、遅延器３１３に供給する。遅延器３１３は、加算器
３１２の加算結果を１サンプル遅延させたものを第２の
係数乗算器３１４に供給する。第２の係数乗算器３１４
は、遅延器３１３の出力をβ倍し、乗算結果を加算器３
１２に供給する。ここに、αは任意の正数、βは０＜β
＜１なる定数である。β＝１−α（０＜α＜１）とする
と、積分器の出力信号３０２は、入力信号３０１の重み
付き移動平均となる。
２２２1,…，２２２L ，２２６の第３の具体例である１
次の再帰形積分器を示す。この積分器は、第１の係数乗
算器３１１、加算器３１２、遅延器３１３、第２の係数
乗算器３１４から構成され、入力信号３０１を積分し
て、結果を出力信号３０２とする。第１の係数乗算器３
１１は積分器の入力信号３０１をα倍し、乗算結果を加
算器３１２に供給する。加算器３１２は第１および第２
の係数乗算器３１１，３１４の乗算結果を加算し、加算
結果を積分器の出力信号３０２として出力するととも
に、遅延器３１３に供給する。遅延器３１３は、加算器
３１２の加算結果を１サンプル遅延させたものを第２の
係数乗算器３１４に供給する。第２の係数乗算器３１４
は、遅延器３１３の出力をβ倍し、乗算結果を加算器３
１２に供給する。ここに、αは任意の正数、βは０＜β
＜１なる定数である。β＝１−α（０＜α＜１）とする
と、積分器の出力信号３０２は、入力信号３０１の重み
付き移動平均となる。
【００８７】本実施例では、第１の係数乗算器３１１で
積分器の入力信号３０１をβ倍し、乗算結果を加算器３
１２に供給し、加算器３１２の出力を積分器の出力信号
３０２としているが、係数乗算器３１１の位置を変え
て、積分器の入力信号３０１を加算器３１２に供給し、
加算器３１２の出力を第１の係数乗算器３１１でβ倍
し、乗算結果を積分器の出力信号３０２としてもよい。
積分器の入力信号３０１をβ倍し、乗算結果を加算器３
１２に供給し、加算器３１２の出力を積分器の出力信号
３０２としているが、係数乗算器３１１の位置を変え
て、積分器の入力信号３０１を加算器３１２に供給し、
加算器３１２の出力を第１の係数乗算器３１１でβ倍
し、乗算結果を積分器の出力信号３０２としてもよい。
【００８８】図１１では、１次の積分器の例を示した
が、積分器の次数は、任意の次数を使用できる。再帰形
積分器は、小さなハードウエアで、長時間にわたる積分
が可能であり、ノイズの影響を受けにくいという特徴が
ある。
が、積分器の次数は、任意の次数を使用できる。再帰形
積分器は、小さなハードウエアで、長時間にわたる積分
が可能であり、ノイズの影響を受けにくいという特徴が
ある。
【００８９】以上、図７から図１１を用いて、２信号間
遅延時間差推定回路２１０の他の実施例を説明した。図
４に示した制御回路１０３の第１の実施例では、遅延時
間差１０２が変化した時に、直ちに受信信号の選択およ
びフィルタ係数の切替を行なっていた。本発明の他の実
施例として、遅延時間差１０２が予め定められた期間一
定値を保つまでの間、受信信号の選択およびフィルタ係
数の切替を行なわないことによって、ノイズなどの影響
によって遅延時間差の推定を誤った時に受信信号の選択
およびフィルタ係数の切替を行なって、エコー除去性能
が低下することを防ぐこともできる。
遅延時間差推定回路２１０の他の実施例を説明した。図
４に示した制御回路１０３の第１の実施例では、遅延時
間差１０２が変化した時に、直ちに受信信号の選択およ
びフィルタ係数の切替を行なっていた。本発明の他の実
施例として、遅延時間差１０２が予め定められた期間一
定値を保つまでの間、受信信号の選択およびフィルタ係
数の切替を行なわないことによって、ノイズなどの影響
によって遅延時間差の推定を誤った時に受信信号の選択
およびフィルタ係数の切替を行なって、エコー除去性能
が低下することを防ぐこともできる。
【００９０】図１２に制御回路１０３の第２の具体例を
示す。遅延器４１１は、遅延時間差１０２を１サンプル
遅延させる。比較器４１２は、遅延器４１１の出力と遅
延時間差１０２を比較する。経過時間判定回路４１３
は、比較器４１２が最後に遅延器４１１の出力と遅延時
間差１０２が異なると判定した時から予め定められた時
間を経過したかどうかを判定する。経過時間判定回路４
１３は、比較器４１２が最後に遅延器４１１の出力と遅
延時間差１０２が異なると判定した時にカウンタの値を
予め定められた初期値に設定し、そうでないときはカウ
ンタの値から１を減じ、カウンタの値が０になった時
に、予め定められた信号を出力するプリセット機能つき
カウンタを用いて容易に実現できる。
示す。遅延器４１１は、遅延時間差１０２を１サンプル
遅延させる。比較器４１２は、遅延器４１１の出力と遅
延時間差１０２を比較する。経過時間判定回路４１３
は、比較器４１２が最後に遅延器４１１の出力と遅延時
間差１０２が異なると判定した時から予め定められた時
間を経過したかどうかを判定する。経過時間判定回路４
１３は、比較器４１２が最後に遅延器４１１の出力と遅
延時間差１０２が異なると判定した時にカウンタの値を
予め定められた初期値に設定し、そうでないときはカウ
ンタの値から１を減じ、カウンタの値が０になった時
に、予め定められた信号を出力するプリセット機能つき
カウンタを用いて容易に実現できる。
【００９１】受信信号選択回路４１４は、前記経過時間
判定回路４１３が比較器４１２が最後に遅延器４１１の
出力と遅延時間差１０２が異なると判定した時から予め
定められた時間を経過したと判定した時に、遅延時間差
１０２に基づいて、複数の受信信号のうち前記遅延時間
が最も短いものを検出して受信信号選択情報１０４と
し、フィルタ係数選択回路４１５に供給するとともに、
制御回路の第１の出力とする。
判定回路４１３が比較器４１２が最後に遅延器４１１の
出力と遅延時間差１０２が異なると判定した時から予め
定められた時間を経過したと判定した時に、遅延時間差
１０２に基づいて、複数の受信信号のうち前記遅延時間
が最も短いものを検出して受信信号選択情報１０４と
し、フィルタ係数選択回路４１５に供給するとともに、
制御回路の第１の出力とする。
【００９２】フィルタ係数選択回路４１５は、前記経過
時間判定回路４１３が比較器４１２が最後に遅延器４１
１の出力と遅延時間差１０２が異なると判定した時から
予め定められた時間を経過したと判定した時に、前記遅
延時間差１０２および前記フィルタ係数選択情報１０４
に基づいて複数の適応フィルタが使用するフィルタ係数
の組を選択し、制御回路の第２の出力であるフィルタ係
数選択信号１０５とする。係数更新制御回路４０３は、
前記複数の適応フィルタが常にフィルタ係数を更新する
ように、制御回路１０３の第３の出力である係数更新制
御信号１０６を出力する。
時間判定回路４１３が比較器４１２が最後に遅延器４１
１の出力と遅延時間差１０２が異なると判定した時から
予め定められた時間を経過したと判定した時に、前記遅
延時間差１０２および前記フィルタ係数選択情報１０４
に基づいて複数の適応フィルタが使用するフィルタ係数
の組を選択し、制御回路の第２の出力であるフィルタ係
数選択信号１０５とする。係数更新制御回路４０３は、
前記複数の適応フィルタが常にフィルタ係数を更新する
ように、制御回路１０３の第３の出力である係数更新制
御信号１０６を出力する。
【００９３】図４に示した制御回路１０３の第１の具体
例および図１２に示した制御回路１０３の第２の実施例
では、常にフィルタ係数を更新するように係数更新制御
信号１０６を出力していた。本発明の他の実施例とし
て、遅延時間差１０２が予め定められた期間一定値を保
つまでの間、フィルタ係数の更新を行なわないことによ
って、ノイズなどの影響によって遅延時間差の推定を誤
った時や、話者が交替した直後で複数の話者が発する音
声によるエコーが混在している時にフィルタ係数を更新
して、フィルタ係数が誤った値になることを防ぐことも
できる。
例および図１２に示した制御回路１０３の第２の実施例
では、常にフィルタ係数を更新するように係数更新制御
信号１０６を出力していた。本発明の他の実施例とし
て、遅延時間差１０２が予め定められた期間一定値を保
つまでの間、フィルタ係数の更新を行なわないことによ
って、ノイズなどの影響によって遅延時間差の推定を誤
った時や、話者が交替した直後で複数の話者が発する音
声によるエコーが混在している時にフィルタ係数を更新
して、フィルタ係数が誤った値になることを防ぐことも
できる。
【００９４】図１３に制御回路１０３の第３の具体例を
示す。受信信号選択回路４０１およびフィルタ係数選択
回路４０２は、図４に示した第１の実施例と全く等しい
受信信号選択情報１０４およびフィルタ係数選択信号１
０５を出力する。遅延器４１１および比較器４１２は、
図１２に示した第２の実施例と全く同様にして、遅延時
間差１０２の変化を検出する。経過時間判定回路４２１
は、比較器４１２が最後に遅延器４１１の出力と遅延時
間差１０２が異なると判定した時から予め定められた時
間を経過したかどうかを判定する。係数更新制御回路４
２２は、前記比較器４１２が遅延器４１１の出力と遅延
時間差１０２が異なると判定した時にフィルタ係数更新
を停止する係数更新制御信号１０６を、前記経過時間判
定回路４２１が、比較器４１２が最後に遅延器４１１の
出力と遅延時間差１０２が異なると判定した時から予め
定められた時間を経過した時にフィルタ係数更新を開始
する係数更新制御信号１０６を出力する。
示す。受信信号選択回路４０１およびフィルタ係数選択
回路４０２は、図４に示した第１の実施例と全く等しい
受信信号選択情報１０４およびフィルタ係数選択信号１
０５を出力する。遅延器４１１および比較器４１２は、
図１２に示した第２の実施例と全く同様にして、遅延時
間差１０２の変化を検出する。経過時間判定回路４２１
は、比較器４１２が最後に遅延器４１１の出力と遅延時
間差１０２が異なると判定した時から予め定められた時
間を経過したかどうかを判定する。係数更新制御回路４
２２は、前記比較器４１２が遅延器４１１の出力と遅延
時間差１０２が異なると判定した時にフィルタ係数更新
を停止する係数更新制御信号１０６を、前記経過時間判
定回路４２１が、比較器４１２が最後に遅延器４１１の
出力と遅延時間差１０２が異なると判定した時から予め
定められた時間を経過した時にフィルタ係数更新を開始
する係数更新制御信号１０６を出力する。
【００９５】図４に示した制御回路１０３第１の実施例
および図１３に示した制御回路１０３第３の実施例で
は、遅延時間差１０２が変化した時に、直ちに受信信号
の選択およびフィルタ係数の切替を行なっていた。ま
た、図４に示した制御回路１０３第１の実施例および図
１２に示した制御回路１０３第２の実施例では、常にフ
ィルタ係数を更新するように係数更新制御信号１０６を
出力していた。本発明の他の実施例として、遅延時間差
１０２が予め定められた期間一定値を保つまでの間、受
信信号の選択、フィルタ係数の切替え、フィルタ係数の
更新を行なわないことによって、エコー除去性能の低下
およびフィルタ係数の誤収束を防ぐこともできる。
および図１３に示した制御回路１０３第３の実施例で
は、遅延時間差１０２が変化した時に、直ちに受信信号
の選択およびフィルタ係数の切替を行なっていた。ま
た、図４に示した制御回路１０３第１の実施例および図
１２に示した制御回路１０３第２の実施例では、常にフ
ィルタ係数を更新するように係数更新制御信号１０６を
出力していた。本発明の他の実施例として、遅延時間差
１０２が予め定められた期間一定値を保つまでの間、受
信信号の選択、フィルタ係数の切替え、フィルタ係数の
更新を行なわないことによって、エコー除去性能の低下
およびフィルタ係数の誤収束を防ぐこともできる。
【００９６】図１４に制御回路１０３の第４の具体例を
示す。遅延器４１１および比較器４１２は、第２および
第３の実施例と全く同様にして、遅延時間差１０２の変
化を検出する。第１の経過時間判定回路４１３、受信信
号選択回路４１４およびフィルタ係数選択回路４１５
は、遅延時間差１０２および比較器４１２の比較結果に
基づいて、図１２に示した第２の実施例と全く同様にし
て、受信信号選択情報１０４およびフィルタ係数選択信
号１０５を出力する。第２の経過時間判定回路４２１お
よび係数更新制御回路４２２は、比較器４１２の比較結
果に基づいて、図１３に示した第３の実施例と全く同様
にして、係数更新制御信号１０６を出力する。第１の経
過時間判定回路４１３および第２の経過時間判定回路４
２１は、互いに異なる経過時間を設定できる。
示す。遅延器４１１および比較器４１２は、第２および
第３の実施例と全く同様にして、遅延時間差１０２の変
化を検出する。第１の経過時間判定回路４１３、受信信
号選択回路４１４およびフィルタ係数選択回路４１５
は、遅延時間差１０２および比較器４１２の比較結果に
基づいて、図１２に示した第２の実施例と全く同様にし
て、受信信号選択情報１０４およびフィルタ係数選択信
号１０５を出力する。第２の経過時間判定回路４２１お
よび係数更新制御回路４２２は、比較器４１２の比較結
果に基づいて、図１３に示した第３の実施例と全く同様
にして、係数更新制御信号１０６を出力する。第１の経
過時間判定回路４１３および第２の経過時間判定回路４
２１は、互いに異なる経過時間を設定できる。
【００９７】以上の説明において、図１に示す第１およ
び第２の受信信号１，２と第１および第２の混在信号１
４，１５を有する２チャンネルの場合を例にとって説明
した部分も、複数の受信信号と、単数または複数の送信
信号が存在する場合に容易に拡張可能である。また、受
信信号がスピーカから空間音響経路を伝搬して、マイク
で収録される音響エコーを除去する音響エコーを例にと
っているが、音響エコー以外のエコー、例えば回線の漏
話などによるエコーに対しても適用できる。また、第１
及び第２の適応フィルタ１０９，１１０としては、ＬＭ
Ｓアルゴリズムによるトランスバーサル形適応フィルタ
を用いた例を示したが、本発明では、任意の適応フィル
タが使用可能である。
び第２の受信信号１，２と第１および第２の混在信号１
４，１５を有する２チャンネルの場合を例にとって説明
した部分も、複数の受信信号と、単数または複数の送信
信号が存在する場合に容易に拡張可能である。また、受
信信号がスピーカから空間音響経路を伝搬して、マイク
で収録される音響エコーを除去する音響エコーを例にと
っているが、音響エコー以外のエコー、例えば回線の漏
話などによるエコーに対しても適用できる。また、第１
及び第２の適応フィルタ１０９，１１０としては、ＬＭ
Ｓアルゴリズムによるトランスバーサル形適応フィルタ
を用いた例を示したが、本発明では、任意の適応フィル
タが使用可能である。
【００９８】例えば、文献５に記載されている学習同定
法によるトランスバーサル形適応フィルタを用いた場合
には、適応フィルタの出力であるエコーレプリカは（２
７）式および（２８）式と同じであり、フィルタ係数の
更新は、
法によるトランスバーサル形適応フィルタを用いた場合
には、適応フィルタの出力であるエコーレプリカは（２
７）式および（２８）式と同じであり、フィルタ係数の
更新は、
【数１４】 で与えられる。トランスバーサル形適応フィルタの代わ
りに、再帰形を用いてもよい。また、サブバンド型適応
フィルタや変換領域の適応フィルタを用いてもよい。
りに、再帰形を用いてもよい。また、サブバンド型適応
フィルタや変換領域の適応フィルタを用いてもよい。
【００９９】
【発明の効果】本発明の多チャンネルエコー除去方法お
よび装置は、話者の位置に応じて適応フィルタの係数を
切替えることによって話者の移動や交替に高速に追従で
きるので、従来の技術では話者が移動した時や話者が交
替した時からフィルタ係数が再び最適値に収束するまで
の間、エコー除去性能が低下するという問題を解決でき
る。
よび装置は、話者の位置に応じて適応フィルタの係数を
切替えることによって話者の移動や交替に高速に追従で
きるので、従来の技術では話者が移動した時や話者が交
替した時からフィルタ係数が再び最適値に収束するまで
の間、エコー除去性能が低下するという問題を解決でき
る。
【図１】本発明の多チャンネルエコー除去装置のブロッ
ク図である。
ク図である。
【図２】図１における遅延時間差推定回路の第１の具体
例を示すブロック図である。
例を示すブロック図である。
【図３】図２における２信号間遅延時間差推定回路の第
１の具体例を示すブロック図である。
１の具体例を示すブロック図である。
【図４】図１における制御回路の第１の具体例を示すブ
ロック図である。
ロック図である。
【図５】図１における係数切替え可能適応フィルタの具
体例を示すブロック図である。
体例を示すブロック図である。
【図６】図５の係数切替え可能適応フィルタにおける演
算回路のブロック図である。
算回路のブロック図である。
【図７】図１における２信号間遅延時間差推定回路の第
２の具体例のブロック図である。
２の具体例のブロック図である。
【図８】図７における相互相関関数計算回路のブロック
図である。
図である。
【図９】図８における積分器の第１の具体例であるトラ
ンスバーサル形積分器のブロック図である。
ンスバーサル形積分器のブロック図である。
【図１０】図８における積分器の第２の具体例である単
純平均計算回路のブロック図である。
純平均計算回路のブロック図である。
【図１１】図８における積分器の第３の具体例である再
帰形積分器のブロック図である。
帰形積分器のブロック図である。
【図１２】図１における制御回路の第２の具体例のブロ
ック図である。
ック図である。
【図１３】図１における制御回路の第３の具体例のブロ
ック図である。
ック図である。
【図１４】図１における制御回路の第４の具体例のブロ
ック図である。
ック図である。
【図１５】従来の多チャンネルエコー除去装置のブロッ
ク図である。
ク図である。
１，２ 受信信号 ３，４，２６，２７ スピーカ ５，６，７，８ エコー ９，１０，２４，２５ マイク １１，１８，１９ 話者 １２，１３ 送信信号 １４，１５ 混在信号 １６，１７ エコー除去装置の出力信号 ２０，２１，２２，２３ 話者の音声 ３０，３１ テレビ会議室 １００ 本発明の多チャンネルエコー除去装置 １０１ 遅延時間差推定回路 １０２ 遅延時間差 １０３ 制御回路 １０４ 受信信号選択情報 １０５ フィルタ係数選択信号 １０６ 係数更新制御信号 １０７ セレクタ １０８ 選択された受信信号 １０９，１１０ 適応フィルタ １１１，１１２ エコーのレプリカ １１３，１１４ 減算器 １２０ 従来の多チャンネルエコー除去装置 １２１ 検出器 １２２，１２３ 適応フィルタ ２０１1,…，２０１M 受信信号 ２０２ 遅延時間差情報 ２０３1,…，２０１K ２信号間遅延時間差 ２０４1,…，２０４K 制御信号 ２０５ 制御回路 ２１０1,…，２１０K ２信号間遅延時間差推定回路 ２１１，２１２ セレクタ ２１３，２１４ 受信信号 ２１５，２１６ タップドディレイライン ２１７1,…，２１７L,２１８1,…，２１８L タップ
ドディレイラインのタップ出力 ２１９1,…，２１９L,２２０1,…，２２０L,２２３
乗算器 ２２１1,…，２２１L,２２２1,…，２２２L,２２４
積分器 ２１０ ２信号間遅延時間差推定回路 ２３１，２３２ 適応フィルタ ２３３1,…，２３３L,２３４1,…，２３４L 適応フ
ィルタのフィルタ係数 ２３５1,…，２３５L,２３６ 絶対値計算回路 ２３７，２３８ 減算器 ２３９ 判定器 ２１０ ２信号間遅延時間差推定回路 ２５１ 相互相関関数計算回路 ２５２1,…，２５２J 相互相関関数 ２５３1,…，２５３J 絶対値計算回路 ２５４1,…，２５４J 相互相関関数の絶対値 ２５５ 判定器 ３０１ 積分器の入力信号 ３０２ 積分器の出力信号 ３０３ Ｄサンプルの遅延器 ３０４ 遅延器 ３０５ 加算器 ３０６ 係数乗算器 ３１１，３１４ 係数乗算器 ３１２ 加算器 ３１３ 遅延器 ３２１ タップドディレイライン ３２２1,…，３２２N 係数乗算器 ３２３ 加算器 ４０１，４１４ 受信信号選択回路 ４０２，４１５ フィルタ係数選択回路 ４０３，４２２ 係数更新制御回路 ４１１ 遅延器 ４１２ 比較器 ４１３，４２１ 経過時間判定回路 ５０１ エコー除去装置の出力信号 ５０２ 適応フィルタの出力信号 ５０３ 信号メモリ ５０４1,…，５０４J 係数メモリ ５０５ セレクタ ５０６ 演算回路 ５１０，５１２ アドレスバス ５１１，５１３ データバス ５１４，５１５，５１６ 入力端子 ５１７ 出力端子 ５１８，５１９ アドレス生成回路 ５２０ 乗算器 ５２１，５２２，５２４ セレクタ ５２３ ＡＬＵ ５２５1,…，５０４I レジスタ ５２６ レジスタバンク ５２７ プログラムメモリ ５２８ 制御回路
ドディレイラインのタップ出力 ２１９1,…，２１９L,２２０1,…，２２０L,２２３
乗算器 ２２１1,…，２２１L,２２２1,…，２２２L,２２４
積分器 ２１０ ２信号間遅延時間差推定回路 ２３１，２３２ 適応フィルタ ２３３1,…，２３３L,２３４1,…，２３４L 適応フ
ィルタのフィルタ係数 ２３５1,…，２３５L,２３６ 絶対値計算回路 ２３７，２３８ 減算器 ２３９ 判定器 ２１０ ２信号間遅延時間差推定回路 ２５１ 相互相関関数計算回路 ２５２1,…，２５２J 相互相関関数 ２５３1,…，２５３J 絶対値計算回路 ２５４1,…，２５４J 相互相関関数の絶対値 ２５５ 判定器 ３０１ 積分器の入力信号 ３０２ 積分器の出力信号 ３０３ Ｄサンプルの遅延器 ３０４ 遅延器 ３０５ 加算器 ３０６ 係数乗算器 ３１１，３１４ 係数乗算器 ３１２ 加算器 ３１３ 遅延器 ３２１ タップドディレイライン ３２２1,…，３２２N 係数乗算器 ３２３ 加算器 ４０１，４１４ 受信信号選択回路 ４０２，４１５ フィルタ係数選択回路 ４０３，４２２ 係数更新制御回路 ４１１ 遅延器 ４１２ 比較器 ４１３，４２１ 経過時間判定回路 ５０１ エコー除去装置の出力信号 ５０２ 適応フィルタの出力信号 ５０３ 信号メモリ ５０４1,…，５０４J 係数メモリ ５０５ セレクタ ５０６ 演算回路 ５１０，５１２ アドレスバス ５１１，５１３ データバス ５１４，５１５，５１６ 入力端子 ５１７ 出力端子 ５１８，５１９ アドレス生成回路 ５２０ 乗算器 ５２１，５２２，５２４ セレクタ ５２３ ＡＬＵ ５２５1,…，５０４I レジスタ ５２６ レジスタバンク ５２７ プログラムメモリ ５２８ 制御回路
フロントページの続き (56)参考文献 特開 平６−197050（ＪＰ，Ａ) 特開 平４−284732（ＪＰ，Ａ) 1992年電子情報通信学会春季全国大会 講演論文集，〔分冊１〕 （1992−３− 15） ＰＰ．１・158 1991年電子情報通信学会春季全国大会 講演論文集，〔分冊１〕 （1991−３− 15） ＰＰ．１・202 電子情報通信学会技術研究報告ＣＳ84 −178，Ｖｏｌ．84，Ｎｏ．330 （1985 −３−26） ＰＰ．７−14 (58)調査した分野(Int.Cl.6，ＤＢ名) H04B 3/00 - 3/23 H04M 1/58 - 1/60 H03H 15/00 - 21/00 ＩＮＳＰＥＣ（ＤＩＡＬＯＧ) ＪＩＣＳＴファイル（ＪＯＩＳ) ＷＰＩ（ＤＩＡＬＯＧ)
Claims (19)
- 【請求項１】 複数の受信信号と単数または複数の送信
信号を有するシステムで、前記受信信号が空間音響経路
を伝搬することによって、又は回線の漏話等によって生
じるエコーと前記送信信号が混在する混在信号から前記
エコーを除去する際に、前記複数の受信信号から２信号
を選択する複数の組合せと１対１に対応させて、前記２
信号が信号源から前記システムに送られるまでに要する
遅延時間の差で定義される２信号間遅延時間差を複数個
求め、前記複数の遅延時間差に基づいて、前記複数の受
信信号のうち前記遅延時間が最も短い受信信号を検出
し、前記遅延時間が最も短い受信信号を入力とし、前記
混在信号と１対１に対応する適応フィルタによって生成
された前記エコーのレプリカを対応する前記混在信号か
ら差し引いたものを出力信号とし、前記適応フィルタの
係数は、前記出力信号を最小とするように更新する多チ
ャンネルエコー除去方法において、前記適応フィルタ
は、複数組のフィルタ係数を具備し、前記複数の遅延時
間差に応じて使用するフィルタ係数の組を切替えること
を特徴とする多チャンネルエコー除去方法。 - 【請求項２】 前記フィルタ係数の組を切替える際に、
前記複数の遅延時間差が変化した時刻から前記複数の遅
延時間差が予め定められた期間一定値を保つまでの間、
前記フィルタ係数の組を切替えないことを特徴とする請
求項１に記載の多チャンネルエコー除去方法。 - 【請求項３】 前記フィルタ係数の更新において、前記
複数の遅延時間差が変化した時刻から、前記複数の遅延
時間差が予め定められた期間一定値を保つまでの間、前
記適応フィルタの係数更新を停止することを特徴とする
請求項１または２に記載の多チャンネルエコー除去方
法。 - 【請求項４】 前記複数の遅延時間差を推定する際に、
前記複数の受信信号から任意の２信号を選択し、第１の
トランスバーサル形適応フィルタで前記２信号のうちの
第１の受信信号を入力として第２の受信信号を予測し、
第２のトランスバーサル形適応フィルタで前記第２の受
信信号を入力として前記第１の受信信号を予測し、前記
第１の適応フィルタの係数と前記第２の適応フィルタの
係数のうち最も絶対値が大きいフィルタ係数から前記２
信号間遅延時間差を推定する操作を、前記複数の遅延時
間差を全て求めるまで繰り返すことを特徴とする請求項
１、２又は３に記載の多チャンネルエコー除去方法。 - 【請求項５】 前記複数の遅延時間差を推定する際に、
前記複数の受信信号から任意の２信号を選択し、前記２
信号間の相互相関関数の絶対値を最大とする時間差を求
め、前記時間差を前記２信号間遅延時間差の推定値とす
る操作を、前記複数の遅延時間差を全て求めるまで繰り
返すことを特徴とする請求項１、２又は３に記載の多チ
ャンネルエコー除去方法。 - 【請求項６】 前記２信号間の相互相関関数の計算にお
いて、前記２信号をそれぞれ第１及び第２の信号とし、
第１の時刻における前記第１の信号の値と前記第１の時
刻と予め定められた時間差を有する第２の時刻における
前記第２の信号の値との乗算をし、前記第１の時刻を予
め定められた複数の時刻に変化させて前記乗算を行なっ
て得られた複数の乗算結果の重みつき移動平均を求め、
前記重みつき移動平均値をもって前記２信号の前記複数
の時間差に対する相互相関関数の値とすることを特徴と
する請求項５に記載の多チャンネルエコー除去方法。 - 【請求項７】 前記重みつき移動平均を再帰形の積分に
置き換えたことを特徴とする請求項６に記載の多チャン
ネルエコー除去方法。 - 【請求項８】 複数の受信信号と単数または複数の送信
信号を有するシステムで、前記受信信号が空間音響経路
を伝搬することによって、又は回線の漏話等によって生
じるエコーと前記送信信号が混在する混在信号から前記
エコーを除去する際に、前記複数の受信信号全てを入力
し、前記複数の遅延時間差を推定し出力する遅延時間差
推定回路と、前記遅延時間差推定回路の推定結果に基づ
いて、前記受信信号のうち前記遅延時間が最も短いもの
を検出して受信信号選択情報とし、前記複数の適応フィ
ルタが使用するフィルタ係数の組を決定してフィルタ係
数選択信号とし、前記複数の適応フィルタが係数更新を
行なうかどうかを決定して係数更新制御信号とする制御
回路と、前記受信信号選択情報に基づいて入力された前
記複数の受信信号のうち前記遅延時間が最も短いものを
選択して出力するセレクタと、前記混在信号と１対１に
対応し、前記フィルタ係数選択信号に基づいて複数のフ
ィルタ係数の組の中から使用するフィルタ係数の組を選
択し、前記セレクタで選択された受信信号を入力して前
記エコーのレプリカを生成する複数の適応フィルタと、
前記混在信号と１対１に対応し、前記生成されたエコー
のレプリカを前記混在信号から差し引いてエコー除去装
置の出力信号とする複数の減算器とを少くとも具備し、
前記複数の適応フィルタは、前記係数更新制御信号に基
づいて前記フィルタ係数選択信号で指定されたフィルタ
係数の組を前記複数の減算器の出力を最小とするように
更新することを特徴とする多チャンネルエコー除去装
置。 - 【請求項９】 前記制御回路が、前記遅延時間差推定回
路の推定結果に基づいて前記受信信号のうち前記遅延時
間が最も短いものを検出して受信信号選択情報とする受
信信号選択回路と、前記受信信号選択情報と前記遅延時
間差推定回路の推定結果に基づいて前記複数の適応フィ
ルタが使用するフィルタ係数の組を選択してフィルタ係
数選択信号とするフィルタ係数選択回路と、前記複数の
適応フィルタが常に係数更新を行なうように、係数更新
制御信号を出力する係数更新制御回路から構成されるこ
とを特徴とする請求項８に記載の多チャンネルエコー除
去装置。 - 【請求項１０】 前記制御回路が、前記遅延時間差推定
回路の推定結果を遅延させる遅延器と、前記遅延器の出
力と前記遅延時間差推定回路の推定結果を比較すること
によって前記遅延時間差推定回路の推定結果の変化を検
出する比較器と、前記比較器が最後に前記遅延時間差推
定回路の推定結果の変化を検出した時からの経過時間を
算出し、前記経過時間が予め定められた時間に達したか
どうかを判定する経過時間判定回路と、前記経過時間判
定回路が前記経過時間が予め定められた時間に達したと
判定した時に、前記遅延時間差推定回路の推定結果に基
づいて前記受信信号のうち最も前記遅延時間が短いもの
を検出して受信信号選択情報とする受信信号選択回路
と、前記経過時間判定回路が前記経過時間が予め定めら
れた時間に達したと判定した時に、前記受信信号選択情
報と前記遅延時間差推定回路の推定結果に基づいて前記
複数の適応フィルタが使用するフィルタ係数の組を選択
してフィルタ係数選択信号とするフィルタ係数選択回路
と、前記複数の適応フィルタに常に係数更新を行なうよ
うに係数更新制御信号を出力する係数更新制御回路とか
ら構成されることを特徴とする請求項８に記載の多チャ
ンネルエコー除去装置。 - 【請求項１１】 前記制御回路が、前記遅延時間差推定
回路の推定結果に基づいて前記受信信号のうち前記遅延
時間が最も短いものを検出して受信信号選択情報とする
受信信号選択回路と、前記受信信号選択情報と前記遅延
時間差推定回路の推定結果に基づいて前記複数の適応フ
ィルタが使用するフィルタ係数の組を選択してフィルタ
係数選択信号とするフィルタ係数選択回路と、前記遅延
時間差推定回路の推定結果を遅延させる遅延器と、前記
遅延器の出力と前記遅延時間差推定回路の推定結果を比
較することによって前記時間差推定回路の推定結果の変
化を検出する比較器と、前記比較器が最後に前記遅延時
間差推定回路の推定結果の変化を検出した時からの経過
時間を算出し、前記経過時間が予め定められた時間に達
したかどうかを判定する経過時間判定回路と、前記比較
器が前記遅延時間差推定回路の推定結果の変化を検出し
た時にフィルタ係数更新を停止する係数更新制御信号
を、前記経過時間判定回路が前記経過時間が予め定めら
れた時間に達したと判定した時にフィルタ係数更新を開
始する係数更新制御信号を出力する係数更新制御回路と
から構成されることを特徴とする請求項８に記載の多チ
ャンネルエコー除去装置。 - 【請求項１２】 前記制御回路が、前記遅延時間差推定
回路の推定結果を遅延させる遅延器と、前記遅延器の出
力と前記遅延時間差推定回路の推定結果を比較すること
によって前記遅延時間差推定回路の推定結果の変化を検
出する比較器と、前記比較器が最後に前記遅延時間差推
定回路の推定結果の変化を検出した時からの経過時間を
算出し、前記経過時間が予め定められた第１の時間に達
したかどうかを判定する第１の経過時間判定回路と、前
記第１の経過時間判定回路が前記経過時間が予め定めら
れた第１の時間に達したと判定した時に、前記遅延時間
差推定回路の推定結果に基づいて前記受信信号のうち最
も前記遅延時間が短いものを検出して受信信号選択情報
とする受信信号選択回路と、前記第１の経過時間判定回
路が前記経過時間が予め定められた第１の時間に達した
と判定した時に、前記受信信号選択情報と前記遅延時間
差推定回路の推定結果に基づいて前記複数の適応フィル
タが使用するフィルタ係数の組を選択してフィルタ係数
選択信号とするフィルタ係数選択回路と、前記比較器が
最後に前記遅延時間差推定回路の推定結果の変化を検出
した時からの経過時間を算出し、前記経過時間が予め定
められた第２の時間に達したかどうかを判定する第２の
経過時間判定回路と、前記比較器が前記遅延時間差推定
回路の推定結果の変化を検出した時にフィルタ係数更新
を停止する係数更新制御信号を、前記第２の経過時間判
定回路が前記経過時間が予め定められた第２の時間に達
したと判定した時にフィルタ係数更新を開始する係数更
新制御信号を出力する係数更新制御回路とで構成される
ことを特徴とする請求項８に記載の多チャンネルエコー
除去装置。 - 【請求項１３】 前記遅延時間差推定回路が、前記複数
の受信信号を入力とし、前記時間差を推定する２つの受
信信号を推定する制御信号に応答し、前記指定された２
信号間の遅延時間差を推定し、推定結果を出力する少な
くとも１つの２信号間遅延時間差推定回路と、前記２信
号間遅延時間差推定回路の各々の出力に応答し、前記複
数の時間差を求めるように前記制御信号を前記２信号間
遅延時間差推定回路の各々に出力し、前記複数の遅延時
間差推定結果を出力する制御回路とで構成されることを
特徴とする請求項８、９、１０、１１又は１２に記載の
多チャンネルエコー除去装置。 - 【請求項１４】 前記２信号間遅延時間差推定回路が、
前記複数の受信信号を入力とし、前記制御信号で指定さ
れた２信号の一方を第１の受信信号として出力する第１
のセレクタと、前記複数の受信信号を入力とし、前記制
御信号で指定された２信号の他方を第２の受信信号とし
て出力する第２のセレクタと、前記第１の受信信号を入
力し、前記第２の受信信号を予測する第１のトランスバ
ーサル形適応フィルタと、前記第２の受信信号を入力
し、前記第１の受信信号を予測する第２のトランスバー
サル形適応フィルタと、前記第１の適応フィルタの係数
の絶対値を求める複数の絶対値計算回路からなる第１の
絶対値計算回路群と、前記第２の適応フィルタの係数の
絶対値を求める複数の絶対値計算回路からなる第２の絶
対値計算回路群と、前記第１および第２の絶対値計算回
路群の各絶対値計算回路の出力に基づいて前記第１およ
び第２の受信信号間の遅延時間差を推定する判定器とか
ら構成されることを特徴とする請求項１３に記載の多チ
ャンネルエコー除去装置。 - 【請求項１５】 前記２信号間遅延時間差推定回路が、
前記複数の受信信号を入力とし、前記制御信号で指定さ
れた２信号の一方を第１の受信信号として出力する第１
のセレクタと、前記複数の受信信号を入力とし、前記制
御信号で指定された２信号の他方を第２の受信信号とし
て出力する第２のセレクタと、前記第１の受信信号と第
２の受信信号の相互相関関数のうち、予め定められた複
数の時間差に対応する値を計算する相互相関関数計算回
路と、前記複数の時間差に対応する相互相関関数の絶対
値を求める複数の絶対値計算回路からなる絶対値計算回
路群と、前記絶対値計算回路群の各絶対値計算回路の出
力に基づいて前記相互相関関数の絶対値を最大とする時
間差を前記第１および第２の受信信号間遅延時間差の推
定値として出力する判定器とから構成されることを特徴
とする請求項１３に記載の多チャンネルエコー除去装
置。 - 【請求項１６】 前記相互相関関数計算回路が、前記第
１の受信信号を遅延させる第１のタップドディレイライ
ンと、前記第２の受信信号を遅延させる第２のタップド
ディレイラインと、前記第１のタップドディレイライン
の各タップ出力と前記第２の受信信号を互いに乗算する
複数の乗算器からなる第１の乗算器群と、前記第１の乗
算器群の各乗算器と１対１に対応し、各乗算器の出力を
積分する複数の積分器からなる第１の積分器群と、前記
第２のタップドディレイラインの各タップ出力と前記第
１の受信信号を互いに乗算する複数の乗算器からなる第
２の乗算器群と、前記第２の乗算器群の各乗算器と１対
１に対応し、各乗算器出力を積分する複数の積分器から
なる第２の積分器群と、前記第１および第２の受信信号
を互いに乗算する第１の乗算器と、前記第１の乗算器の
出力を積分する第１の積分器とから構成され、前記第１
および第２の積分器群と前記第１の積分器の出力を前記
第１および第２の受信信号の相互相関関数として出力す
ることを特徴とする請求項１５に記載の多チャンネルエ
コー除去装置。 - 【請求項１７】 前記積分器が、入力信号を入力する入
力端と、出力信号を出力する出力端と、前記入力信号を
遅延させるタップドディレイラインと、前記タップドデ
ィレイラインの各タップ出力を定数倍する複数の係数乗
算器と、前記複数の係数乗算器の出力の総和を求め、前
記総和を積分結果を表す前記出力信号として出力する加
算器とから構成されることを特徴とする請求項１６に記
載の多チャンネルエコー除去装置。 - 【請求項１８】 前記積分器が、入力信号を入力する入
力端と、出力信号を出力する出力端と、前記入力信号を
遅延させる第１の遅延器と、１サンプル前の出力信号を
格納する第２の遅延器と、前記第２の遅延器の値と前記
入力信号を加え、前記第１の遅延器の出力を差し引いた
ものを前記出力信号とすると共に、前記第２の遅延器に
格納する加算器とから構成されることを特徴とする請求
項１６に記載の多チャンネルエコー除去装置。 - 【請求項１９】 前記積分器が、入力信号を入力する入
力端と、出力信号を出力する出力端と、前記入力信号を
定数倍する第１の係数乗算器と、前記出力信号を遅延さ
せるタップドディレイラインと、前記タップドディレイ
ラインの各タップ出力を定数倍する複数の係数乗算器
と、前記複数の係数乗算器の出力と前記第１の係数乗算
器の出力の総和を求め、前記総和を前記出力信号とする
と共に前記タップドディレイラインに格納する加算器と
から構成されることを特徴とする請求項１６に記載の多
チャンネルエコー除去装置。
Priority Applications (5)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP4046185A JP2792311B2 (ja) | 1992-01-31 | 1992-01-31 | 多チャンネルエコー除去方法および装置 |
| DE69321635T DE69321635T2 (de) | 1992-01-31 | 1993-02-01 | Mehrkanalsechokompensation mit adaptiven Filtern welche auswählbare Koeffizientvektoren haben |
| CA002088558A CA2088558C (en) | 1992-01-31 | 1993-02-01 | Multi-channel echo cancellation with adaptive filters having selectable coefficient vectors |
| EP93101526A EP0561133B1 (en) | 1992-01-31 | 1993-02-01 | Multi-channel echo cancellation with adaptive filters having selectable coefficient vectors |
| US08/012,125 US5371789A (en) | 1992-01-31 | 1993-02-01 | Multi-channel echo cancellation with adaptive filters having selectable coefficient vectors |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP4046185A JP2792311B2 (ja) | 1992-01-31 | 1992-01-31 | 多チャンネルエコー除去方法および装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH05218906A JPH05218906A (ja) | 1993-08-27 |
| JP2792311B2 true JP2792311B2 (ja) | 1998-09-03 |
Family
ID=12739988
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP4046185A Expired - Fee Related JP2792311B2 (ja) | 1992-01-31 | 1992-01-31 | 多チャンネルエコー除去方法および装置 |
Country Status (5)
| Country | Link |
|---|---|
| US (1) | US5371789A (ja) |
| EP (1) | EP0561133B1 (ja) |
| JP (1) | JP2792311B2 (ja) |
| CA (1) | CA2088558C (ja) |
| DE (1) | DE69321635T2 (ja) |
Families Citing this family (83)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH084243B2 (ja) * | 1993-05-31 | 1996-01-17 | 日本電気株式会社 | 多チャンネルエコー除去方法および装置 |
| JP2576767B2 (ja) * | 1993-09-10 | 1997-01-29 | 日本電気株式会社 | エコー消去方法およびエコーキャンセラ |
| US6594688B2 (en) | 1993-10-01 | 2003-07-15 | Collaboration Properties, Inc. | Dedicated echo canceler for a workstation |
| US5473686A (en) * | 1994-02-01 | 1995-12-05 | Tandy Corporation | Echo cancellation apparatus |
| US5500892A (en) * | 1994-02-14 | 1996-03-19 | Brooktree Corporation | Echo canceller |
| DE4420212A1 (de) * | 1994-06-04 | 1995-12-07 | Deutsche Bundespost Telekom | Übertragungssystem für gleichzeitige Mehrfachübertragung von mehreren Bild- und Tonsignalen |
| JP2749780B2 (ja) * | 1994-09-30 | 1998-05-13 | 株式会社エイ・ティ・アール人間情報通信研究所 | 適応化型相互相関装置 |
| JP3397269B2 (ja) * | 1994-10-26 | 2003-04-14 | 日本電信電話株式会社 | 多チャネル反響消去方法 |
| US5828756A (en) * | 1994-11-22 | 1998-10-27 | Lucent Technologies Inc. | Stereophonic acoustic echo cancellation using non-linear transformations |
| US5548642A (en) * | 1994-12-23 | 1996-08-20 | At&T Corp. | Optimization of adaptive filter tap settings for subband acoustic echo cancelers in teleconferencing |
| US5495473A (en) * | 1994-12-28 | 1996-02-27 | Intel Corporation | Method and apparatus for implementing hardware independent echo cancellation in a full-duplex audio stream in a computer system |
| JPH08204619A (ja) * | 1995-01-26 | 1996-08-09 | Nec Corp | エコーサプレッサ |
| US5664019A (en) * | 1995-02-08 | 1997-09-02 | Interval Research Corporation | Systems for feedback cancellation in an audio interface garment |
| JP3359460B2 (ja) * | 1995-03-20 | 2002-12-24 | 沖電気工業株式会社 | 適応フィルタ及びエコーキャンセラ |
| US5677951A (en) * | 1995-06-19 | 1997-10-14 | Lucent Technologies Inc. | Adaptive filter and method for implementing echo cancellation |
| JP2856113B2 (ja) * | 1995-07-10 | 1999-02-10 | 日本電気株式会社 | エコーキャンセラ |
| US5663955A (en) * | 1995-08-25 | 1997-09-02 | Lucent Technologies Inc. | Echo canceller system with shared coefficient memory |
| US5664011A (en) * | 1995-08-25 | 1997-09-02 | Lucent Technologies Inc. | Echo canceller with adaptive and non-adaptive filters |
| US6002776A (en) * | 1995-09-18 | 1999-12-14 | Interval Research Corporation | Directional acoustic signal processor and method therefor |
| DE69628618T2 (de) * | 1995-09-26 | 2004-05-13 | Nippon Telegraph And Telephone Corp. | Verfahren und Einrichtung zur mehrkanaligen Kompensation eines akustischen Echos |
| AU4996296A (en) * | 1996-02-26 | 1997-09-10 | Interval Research Corporation | Systems for feedback cancellation in an audio interface garment |
| JP2924762B2 (ja) * | 1996-02-28 | 1999-07-26 | 日本電気株式会社 | アダプティブフィルタ及びその適応化方法 |
| JPH09266456A (ja) * | 1996-03-28 | 1997-10-07 | Mitsubishi Electric Corp | エコーキャンセラ |
| FI106326B (fi) * | 1996-09-23 | 2001-01-15 | Tellabs Oy | Kytkentä kaiunpoistavan moninopeusmodeemin sovittamiseksi linjalle sekä kytkentämenetelmä |
| DE19647276A1 (de) * | 1996-11-15 | 1998-05-20 | Alsthom Cge Alcatel | Verfahren und Anordnung zur adaptiven Echokompensation |
| US5953411A (en) * | 1996-12-18 | 1999-09-14 | Intel Corporation | Method and apparatus for maintaining audio sample correlation |
| US6174806B1 (en) | 1997-01-28 | 2001-01-16 | Micron Technology, Inc. | High pressure anneals of integrated circuit structures |
| US5933797A (en) * | 1997-02-28 | 1999-08-03 | Telefonaktiebolaget Lm Ericsson (Publ) | Adaptive dual filter echo cancellation |
| US5999828A (en) * | 1997-03-19 | 1999-12-07 | Qualcomm Incorporated | Multi-user wireless telephone having dual echo cancellers |
| US6430295B1 (en) | 1997-07-11 | 2002-08-06 | Telefonaktiebolaget Lm Ericsson (Publ) | Methods and apparatus for measuring signal level and delay at multiple sensors |
| US6031908A (en) * | 1997-11-14 | 2000-02-29 | Tellabs Operations, Inc. | Echo canceller employing dual-H architecture having variable adaptive gain settings |
| US6505057B1 (en) * | 1998-01-23 | 2003-01-07 | Digisonix Llc | Integrated vehicle voice enhancement system and hands-free cellular telephone system |
| DE69908463T2 (de) * | 1998-03-05 | 2004-05-13 | Nippon Telegraph And Telephone Corp. | Verfahren und Einrichtung zur mehrkanaligen Kompensation eines akustischen Echos |
| DE69942944D1 (de) * | 1998-04-08 | 2010-12-30 | British Telecomm | Telekonferenzsystem |
| JP3319403B2 (ja) * | 1998-09-02 | 2002-09-03 | 日本電気株式会社 | 多チャネルエコー除去方法及び装置 |
| JP2000196507A (ja) * | 1998-12-28 | 2000-07-14 | Nec Corp | 多重化回線用エコ―除去の方法及び装置 |
| DE60029436T2 (de) | 1999-05-25 | 2006-12-28 | British Telecommunications Public Ltd. Co. | Akustischer Echokompensator |
| US6587974B1 (en) * | 1999-11-19 | 2003-07-01 | Tyco Telecommunications (Us) Inc. | Method for fault detection and fast protection switching in transmission systems |
| US6384873B1 (en) * | 1999-12-03 | 2002-05-07 | Thomson Licensing S.A. | Vector magnitude control of a comb filter |
| JP2002009675A (ja) * | 2000-06-26 | 2002-01-11 | Nec Corp | 多重化回線用エコー除去方法及び装置 |
| US20020054685A1 (en) * | 2000-11-09 | 2002-05-09 | Carlos Avendano | System for suppressing acoustic echoes and interferences in multi-channel audio systems |
| WO2002052543A1 (es) * | 2000-12-26 | 2002-07-04 | Lear Automotive (Eeds) Spain, S.L. | Sistema y procedimiento de comunicacion para el habitaculo de un vehiculo automovil |
| US6934655B2 (en) * | 2001-03-16 | 2005-08-23 | Mindspeed Technologies, Inc. | Method and apparatus for transmission line analysis |
| US20020191779A1 (en) * | 2001-03-23 | 2002-12-19 | Tien Pham | System for convolutional echo cancellation by iterative autocorrelation |
| FR2841721B1 (fr) * | 2002-06-28 | 2004-08-20 | France Telecom | Dispositif de traitement d'echo pour systeme de communication de type monovoie ou multivoies |
| US7251213B2 (en) * | 2002-09-17 | 2007-07-31 | At&T Corp. | Method for remote measurement of echo path delay |
| US20040064622A1 (en) * | 2002-09-30 | 2004-04-01 | Smith Winthrop W. | Signal processing resource with sample-by-sample selective characteristics |
| KR100480789B1 (ko) * | 2003-01-17 | 2005-04-06 | 삼성전자주식회사 | 피드백 구조를 이용한 적응적 빔 형성방법 및 장치 |
| US7218730B2 (en) * | 2003-08-13 | 2007-05-15 | International Business Machines Corporation | System and method for analysis and filtering of signals in a telecommunications network |
| JP4192800B2 (ja) * | 2004-02-13 | 2008-12-10 | ソニー株式会社 | 音声集音装置と方法 |
| US20050285935A1 (en) * | 2004-06-29 | 2005-12-29 | Octiv, Inc. | Personal conferencing node |
| US8457614B2 (en) * | 2005-04-07 | 2013-06-04 | Clearone Communications, Inc. | Wireless multi-unit conference phone |
| US8050398B1 (en) | 2007-10-31 | 2011-11-01 | Clearone Communications, Inc. | Adaptive conferencing pod sidetone compensator connecting to a telephonic device having intermittent sidetone |
| 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 |
| US8462193B1 (en) * | 2010-01-08 | 2013-06-11 | Polycom, Inc. | Method and system for processing audio signals |
| GB2493801B (en) | 2011-08-18 | 2014-05-14 | Ibm | Improved audio quality in teleconferencing |
| US9613634B2 (en) * | 2014-06-19 | 2017-04-04 | Yang Gao | Control of acoustic echo canceller adaptive filter for speech enhancement |
| US9565493B2 (en) | 2015-04-30 | 2017-02-07 | Shure Acquisition Holdings, Inc. | Array microphone system and method of assembling the same |
| US9554207B2 (en) | 2015-04-30 | 2017-01-24 | Shure Acquisition Holdings, Inc. | Offset cartridge microphones |
| US9653060B1 (en) * | 2016-02-09 | 2017-05-16 | Amazon Technologies, Inc. | Hybrid reference signal for acoustic echo cancellation |
| US9659555B1 (en) * | 2016-02-09 | 2017-05-23 | Amazon Technologies, Inc. | Multichannel acoustic echo cancellation |
| US10367948B2 (en) | 2017-01-13 | 2019-07-30 | Shure Acquisition Holdings, Inc. | Post-mixing acoustic echo cancellation systems and methods |
| US10200540B1 (en) * | 2017-08-03 | 2019-02-05 | Bose Corporation | Efficient reutilization of acoustic echo canceler channels |
| US10542153B2 (en) | 2017-08-03 | 2020-01-21 | Bose Corporation | Multi-channel residual echo suppression |
| US10594869B2 (en) * | 2017-08-03 | 2020-03-17 | Bose Corporation | Mitigating impact of double talk for residual echo suppressors |
| US10863269B2 (en) | 2017-10-03 | 2020-12-08 | Bose Corporation | Spatial double-talk detector |
| US10090000B1 (en) * | 2017-11-01 | 2018-10-02 | GM Global Technology Operations LLC | Efficient echo cancellation using transfer function estimation |
| US10522167B1 (en) * | 2018-02-13 | 2019-12-31 | Amazon Techonlogies, Inc. | Multichannel noise cancellation using deep neural network masking |
| EP3804356A1 (en) | 2018-06-01 | 2021-04-14 | Shure Acquisition Holdings, Inc. | Pattern-forming microphone array |
| US11297423B2 (en) | 2018-06-15 | 2022-04-05 | Shure Acquisition Holdings, Inc. | Endfire linear array microphone |
| WO2020061353A1 (en) | 2018-09-20 | 2020-03-26 | Shure Acquisition Holdings, Inc. | Adjustable lobe shape for array microphones |
| US11558693B2 (en) | 2019-03-21 | 2023-01-17 | Shure Acquisition Holdings, Inc. | Auto focus, auto focus within regions, and auto placement of beamformed microphone lobes with inhibition and voice activity detection functionality |
| CN113841421A (zh) | 2019-03-21 | 2021-12-24 | 舒尔获得控股公司 | 具有抑制功能的波束形成麦克风瓣的自动对焦、区域内自动对焦、及自动配置 |
| EP3942842A1 (en) | 2019-03-21 | 2022-01-26 | Shure Acquisition Holdings, Inc. | Housings and associated design features for ceiling array microphones |
| US10964305B2 (en) | 2019-05-20 | 2021-03-30 | Bose Corporation | Mitigating impact of double talk for residual echo suppressors |
| WO2020237206A1 (en) | 2019-05-23 | 2020-11-26 | Shure Acquisition Holdings, Inc. | Steerable speaker array, system, and method for the same |
| US11302347B2 (en) | 2019-05-31 | 2022-04-12 | Shure Acquisition Holdings, Inc. | Low latency automixer integrated with voice and noise activity detection |
| EP4018680A1 (en) | 2019-08-23 | 2022-06-29 | Shure Acquisition Holdings, Inc. | Two-dimensional microphone array with improved directivity |
| US11107488B1 (en) * | 2019-10-24 | 2021-08-31 | Amazon Technologies, Inc. | Reduced reference canceller |
| US11552611B2 (en) | 2020-02-07 | 2023-01-10 | Shure Acquisition Holdings, Inc. | System and method for automatic adjustment of reference gain |
| USD944776S1 (en) | 2020-05-05 | 2022-03-01 | Shure Acquisition Holdings, Inc. | Audio device |
| US11706562B2 (en) | 2020-05-29 | 2023-07-18 | Shure Acquisition Holdings, Inc. | Transducer steering and configuration systems and methods using a local positioning system |
| EP4285605A1 (en) | 2021-01-28 | 2023-12-06 | Shure Acquisition Holdings, Inc. | Hybrid audio beamforming system |
Family Cites Families (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4621172A (en) * | 1982-12-22 | 1986-11-04 | Nec Corporation | Fast convergence method and system for echo canceller |
| FR2538975A1 (fr) * | 1982-12-30 | 1984-07-06 | Trt Telecom Radio Electr | Procede utilise dans un dispositif d'annulation d'echo pour la mesure d'un retard d'echo et dispositif de mise en oeuvre de ce procede |
| US4609787A (en) * | 1984-05-21 | 1986-09-02 | Communications Satellite Corporation | Echo canceller with extended frequency range |
| JPS62120733A (ja) * | 1985-11-21 | 1987-06-02 | Nippon Telegr & Teleph Corp <Ntt> | 反響消去装置 |
| JPS62125722A (ja) * | 1985-11-26 | 1987-06-08 | Matsushita Electric Ind Co Ltd | 音響的エコ−消去装置 |
| JPH0650829B2 (ja) * | 1986-09-16 | 1994-06-29 | 日本電気株式会社 | エコ−キヤンセラ方式モデム |
| EP0272714B1 (en) * | 1986-11-17 | 1992-06-24 | Alcatel N.V. | Adjustable echo canceller |
| JP2810052B2 (ja) * | 1987-11-18 | 1998-10-15 | 株式会社日立製作所 | 伝送装置 |
| FR2628918B1 (fr) * | 1988-03-15 | 1990-08-10 | France Etat | Dispositif annuleur d'echo a filtrage en sous-bandes de frequence |
| US4947425A (en) * | 1989-10-27 | 1990-08-07 | At&T Bell Laboratories | Echo measurement arrangement |
-
1992
- 1992-01-31 JP JP4046185A patent/JP2792311B2/ja not_active Expired - Fee Related
-
1993
- 1993-02-01 EP EP93101526A patent/EP0561133B1/en not_active Expired - Lifetime
- 1993-02-01 CA CA002088558A patent/CA2088558C/en not_active Expired - Fee Related
- 1993-02-01 DE DE69321635T patent/DE69321635T2/de not_active Expired - Fee Related
- 1993-02-01 US US08/012,125 patent/US5371789A/en not_active Expired - Lifetime
Non-Patent Citations (3)
| Title |
|---|
| １９９１年電子情報通信学会春季全国大会講演論文集，〔分冊１〕 （１９９１−３−１５） ＰＰ．１・２０２ |
| １９９２年電子情報通信学会春季全国大会講演論文集，〔分冊１〕 （１９９２−３−１５） ＰＰ．１・１５８ |
| 電子情報通信学会技術研究報告ＣＳ８４−１７８，Ｖｏｌ．８４，Ｎｏ．３３０ （１９８５−３−２６） ＰＰ．７−１４ |
Also Published As
| Publication number | Publication date |
|---|---|
| US5371789A (en) | 1994-12-06 |
| JPH05218906A (ja) | 1993-08-27 |
| CA2088558A1 (en) | 1993-08-01 |
| DE69321635D1 (de) | 1998-11-26 |
| EP0561133B1 (en) | 1998-10-21 |
| EP0561133A1 (en) | 1993-09-22 |
| CA2088558C (en) | 1997-12-23 |
| DE69321635T2 (de) | 1999-03-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP2792311B2 (ja) | 多チャンネルエコー除去方法および装置 | |
| JP2792252B2 (ja) | 多チャンネルエコー除去方法および装置 | |
| US5323459A (en) | Multi-channel echo canceler | |
| EP0627825B1 (en) | Multi-channel echo cancelling method and a device thereof | |
| JP2683490B2 (ja) | 適応性雑音除去装置 | |
| JP3152292B2 (ja) | 多チャネルエコー除去方法及び装置 | |
| US5774562A (en) | Method and apparatus for dereverberation | |
| JPH114288A (ja) | エコーキャンセラ装置 | |
| Buchner et al. | Robust extended multidelay filter and double-talk detector for acoustic echo cancellation | |
| JPH08288890A (ja) | 適応型フィルタリングのための方法および装置 | |
| JP2020519070A (ja) | マルチチャネル干渉除去のための装置および方法 | |
| JP3319403B2 (ja) | 多チャネルエコー除去方法及び装置 | |
| JPH07176991A (ja) | 適応フィルタ装置及びその制御方法 | |
| JP2000035788A (ja) | 多重チャネル適応フィルタリング | |
| JP3403473B2 (ja) | ステレオエコーキャンセラ | |
| JP4041770B2 (ja) | 音響エコー消去方法、その装置、プログラム及びその記録媒体 | |
| JP2542625B2 (ja) | 反響消去装置 | |
| KR100272131B1 (ko) | 계층적 구조의 적응반향 제거장치 | |
| JPH11145878A (ja) | 適応フィルタによる未知システム同定の方法及び装置 | |
| JPH04234212A (ja) | アダプティブ・フィルタ | |
| JP6356087B2 (ja) | エコー消去装置、その方法及びプログラム | |
| JP3339612B2 (ja) | 多チャンネル音声通信会議用反響消去方法 | |
| JP2004357053A (ja) | エコーキャンセラ装置およびエコーキャンセラ方法 | |
| JPH09200092A (ja) | 多チャネルエコー除去方法及び装置 | |
| JPH04245810A (ja) | システム特性推定方法及び装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 19980519 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080619 Year of fee payment: 10 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090619 Year of fee payment: 11 |
|
| LAPS | Cancellation because of no payment of annual fees |