JP2007124678A - 適応フィルタを用いた信号処理装置 - Google Patents

適応フィルタを用いた信号処理装置 Download PDF

Info

Publication number
JP2007124678A
JP2007124678A JP2006319008A JP2006319008A JP2007124678A JP 2007124678 A JP2007124678 A JP 2007124678A JP 2006319008 A JP2006319008 A JP 2006319008A JP 2006319008 A JP2006319008 A JP 2006319008A JP 2007124678 A JP2007124678 A JP 2007124678A
Authority
JP
Japan
Prior art keywords
filter
adaptive
output signal
signal power
threshold value
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP2006319008A
Other languages
English (en)
Other versions
JP4928918B2 (ja
Inventor
Hitoshi Nagata
仁史 永田
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.)
Toshiba Corp
Original Assignee
Toshiba 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 Toshiba Corp filed Critical Toshiba Corp
Priority to JP2006319008A priority Critical patent/JP4928918B2/ja
Publication of JP2007124678A publication Critical patent/JP2007124678A/ja
Application granted granted Critical
Publication of JP4928918B2 publication Critical patent/JP4928918B2/ja
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Landscapes

  • Filters That Use Time-Delay Elements (AREA)

Abstract

【課題】目的信号を歪ませたり、背景雑音に対して適応してしまうことなく、確実な適応動作を可能とした適応フィルタを用いた信号処理装置を提供する。
【解決手段】第1の適応フィルタ1と、第1の適応フィルタ1より適応速度の大きい第2の適応フィルタ2とを有する信号処理装置であって、第1の適応フィルタ1は第2の適応フィルタ2の出力信号に基づいてフィルタの更新が制御される。
【選択図】図12

Description

本発明は、適応フィルタを用いた信号処理装置に係り、特に音声認識装置やテレビ会議装置などで音声入力のためにマイクロホンアレイを用いて妨害雑音を抑圧し目的信号の音声を抽出するマイクロホンアレイ処理装置に適した信号処理装置に関する。
配列された複数のマイクロホンからなるマイクロホンアレイを用い、このマイクロホンアレイの出力信号を処理することにより、雑音を抑圧して目的音源からの目的信号のみを抽出するマイクロホンアレイ処理装置は、音声認識装置やテレビ会議装置などにおける音声入力を目的として開発が進められている。中でも、少ない数のマイクロホンで大きな雑音抑圧効果が得られる適応ビームフォーマを利用したマイクロホンアレイ処理装置については、例えば文献1:電子情報通信学会編「音響システムとデジタル処理」、文献2:Heykin:“Adaptive Filter Theory (Plentice Hall)”に述べられているように、一般化サイドローブキャンセラ、フロスト型ビームフォーマおよび参照信号法など、種々の方法が知られている。
適応ビームフォーマ処理は、基本的には妨害雑音の到来方向に死角を有する指向性ビームを形成したフィルタにより妨害雑音を抑圧する処理である。この適応ビームフォーマ処理では、目的信号が入力した際に目的信号の一部が除去されて歪んだり、目的信号と妨害雑音が同時に入力した際に目的信号が歪んだりするという問題を避けるため、目的信号が入力した際に適応速度を小さくすることが行われている。一般には、適応ビームフォーマの出力信号パワーが大きい場合に、目的信号が入力したものとみなして適応を停止する処理が行われる。ところが、出力信号パワーは目的信号が入力した場合だけでなく、妨害雑音などの状況が変化した場合にも増大するため、このように出力信号パワーが大きいときに完全に適応を停止することには問題がある。
一方、文献3: Julie E.Greenburg et.al.:“Evaluation of an adaptive beamforming method for hearing aids ”,pp.1662-1676,Jarnul of Acous.Soc.of Am.91(3),1992 に開示されている補聴器用のマイクロホンアレイでは、適応フィルタの適応アルゴリズムとして一般的なNLMS(Normalized Least Mean Square)法によるフィルタ更新演算の式(1)を改良し、式(2)のように参照信号パワーと出力信号パワーの和に基づいて適応速度を正規化することにより、出力信号パワーの増大時に適応を完全に停止せず、目的信号の歪みを抑えるようにしている。
n+1 =W−2μe/P (1)
n+1 =W−2μe/(P+P) (2)
ここで、nは時刻、Wはフィルタ係数、μはステップサイズ、e(=d−W )は誤差信号、dは希望応答、Xは参照信号xの最新のサンプルから過去のサンプルをフィルタのタップ点数分だけ並べてベクトルとした参照信号ベクトル、Pは参照信号パワー、Pは出力信号パワーである。参照信号xは、一般に妨害雑音と目的信号の一部および背景雑音が重畳した信号となっている。式(1)(2)のいずれの場合も、適応フィルタの適応速度は右辺第2項の大きさで決まる。
図14に、従来の適応フィルタの構成を示す。この適応フィルタは、参照信号xにフィルタ演算を施すフィルタ部101と、希望応答dからフィルタ部101の出力信号yを減算する減算器102と、フィルタ部101のフィルタ係数について式(1)または(2)に示した更新演算を行うフィルタ更新部103からなっている。
通常のNLMS法や、文献3に記載の改良されたNLMS法においては、妨害雑音も目的信号も存在しない状況では、参照信号パワーPが非常に小さくなることから、式(1)(2)より理解されるように適応速度が非常に大きくなるため、適応フィルタが背景雑音に適応してしまうという問題がある。背景雑音に対して適応が進んだ場合、フィルタ係数は本来除去したい信号を除去するのに適したフィルタ係数から遠ざかってしまう。すなわち、妨害雑音が例えば音声のような非定常な信号である場合には、妨害雑音パワーの小さい部分でフィルタ係数が本来の値から遠ざかり、妨害雑音パワーの大きい部分でフィルタ係数が収束するという動作の繰り返しとなり、この繰り返しの過程でフィルタ係数の収束までの間、高いレベルの妨害雑音が出力されてしまう。
この問題に対し、例えばエコーキャンセラでは参照信号のレベルに基づいて適応を制御している。すなわち、参照信号レベルがしきい値以下の場合には適応を停止し、参照信号レベルがしきい値より大きい場合にのみ適応を行うことによって、雑音の影響を低減している。しかし、適応ビームフォーマでは、上述したように妨害雑音と目的信号の一部および背景雑音が重畳したものが参照信号となっているため、単に参照信号レベルが大きいときだけ適応するようにした場合、しきい値レベル以下の妨害雑音があるときはこれを除去できないことになる。
上述したように、適応フィルタにおいて目的信号の入力時に目的信号が歪むという現象を抑えるために、出力信号パワーが大きいときに適応を停止する処理を行う従来の方法では、目的信号が入力した場合だけでなく状況が変化した場合にも出力パワーが増大して適応を停止してしまうため、正しい適応動作を行うことができないという問題がある。
また、従来のNLMS法や改善されたNLMS法においては、妨害雑音も目的信号もなく参照信号パワーが小さくなる状況での適応速度が非常に大きくなるため、適応フィルタが背景雑音に対して適応してしまうという問題があり、この問題を避けるためにエコーキャンセラのように参照信号のレベルがしきい値以下の場合には適応を停止し、しきい値より大きい場合にのみ適応を行う方法では、適応ビームフォーマにおいてしきい値レベル以下の妨害雑音があるときはこれを除去できないという問題が発生する。
本発明は、上述した従来技術の問題点を解決するためになされたもので、目的信号を歪ませたり、背景雑音に対して適応してしまうことなく、確実な適応動作を可能とした適応フィルタを用いた信号処理装置を提供することを目的とする。
上記の課題を解決するため、本発明に係る適応フィルタは、出力信号パワーが予め定めたしきい値以下のときには適応フィルタの更新速度を決める係数であるステップサイズが非零のほぼ一定値であり、出力信号パワーがしきい値を越えた場合にはステップサイズが出力信号パワーの増加に対して減少するように適応モードが制御されることを基本的な特徴とし、目的信号の歪み防止と背景雑音への適応防止の両立を可能としている。
すなわち、適応フィルタに目的信号が入力した場合は、出力信号パワーがしきい値より大きくなり、ステップサイズが出力信号パワーの増加に対して減少することで適応速度が小さくなるため、目的信号が歪んでしまうことはない。一方、適応フィルタに妨害雑音のみ入力した場合と背景雑音のみが入力した場合は、出力信号パワーがしきい値以下となり、ステップサイズ固定の適応フィルタとして動作する。この場合、ステップサイズが非零のほぼ一定値であるために、適応速度は入力レベルに比例することになり、レベルの高い妨害雑音に対しては速く適応できると共に、レベルの低い背景雑音に適応してしまうことが防止される。
より具体的には、本発明に係る適応フィルタは、少なくとも1チャネルの参照信号にフィルタ演算を行って1チャネルの信号を出力するフィルタ手段と、希望応答からフィルタ手段の出力信号を差し引いて適応フィルタの出力信号を得る減算手段と、参照信号と減算手段の出力信号に基づいてフィルタ手段を更新するフィルタ更新手段とを有し、フィルタ更新手段において減算手段の出力信号パワーを予め定められたしきい値と比較し、出力信号パワーがしきい値以下のときはフィルタ手段の更新速度を決める係数を非零のほぼ一定値とし、出力信号パワーがしきい値を越えたときは該係数を出力信号パワーの増加に対して減少させることを特徴とする。
本発明に係る他の適応フィルタは、複数チャネルの入力信号を処理して1チャネルの信号を出力するフィルタ手段と、このフィルタ手段の特性変化範囲に関して予め定めた拘束条件の下で該フィルタ手段を更新するフィルタ更新手段とを有し、フィルタ更新手段においてフィルタ手段の出力信号パワーを予め定められたしきい値と比較し、出力信号パワーがしきい値以下のときはフィルタ手段の更新速度を決める係数を非零のほぼ一定値とし、出力信号パワーがしきい値を越えたときは該係数を出力信号パワーの増加に対して減少させることを特徴とする。
本発明に係る適応フィルタにおけるフィルタ更新手段は、より具体的には適応フィルタの出力信号パワーを予め定められたしきい値と比較し、出力信号パワーがしきい値以下のときはフィルタ手段の更新速度を決める係数を非零のほぼ一定値とし、出力信号パワーがしきい値を越えたときは該係数を出力信号パワーの増加に対して減少させる制御を行う制御手段と、この制御手段により制御された係数を含む更新演算式に従ってフィルタ手段の更新演算を行う更新演算手段とにより構成される。
本発明に係る信号処理装置は、上述した適応フィルタの前に複数チャネルの信号から目的信号に関する成分を除去するブロッキングフィルタを配置し、このブロッキングフィルタの出力信号を適応フィルタに参照信号として入力する。
本発明に係る他の信号処理装置は、第1の適応フィルタと、この第1の適応フィルタより適応速度の大きい第2の適応フィルタとを具備し、第1の適応フィルタは第2の適応フィルタの出力信号に基づいてフィルタの更新が制御されることを特徴とする。このようにすると、収束速度がさらに改善される。
ここで、第1の適応フィルタにおいては、その出力信号パワーを予め定められたしきい値と比較し、出力信号パワーがしきい値以下のときは適応速度がほぼ一定値となり、出力信号パワーがしきい値を越えたときは第1および第2の適応フィルタの出力信号パワーに応じて適応速度が変化するように、フィルタの更新が制御される。さらに具体的には、第1の適応フィルタはその出力信号パワーを予め定められたしきい値と比較し、出力信号パワーがしきい値以下のときはフィルタの更新速度を決める係数を非零のほぼ一定値とし、出力信号パワーがしきい値を越えたときは該係数を第1および第2の適応フィルタの出力信号パワーの比の増加に対して減少させるようにする。
本発明によれば目的信号の歪み防止と雑音への適応防止の両立を可能とした適応フィルタを用いた信号処理装置を提供することができる。
以下、図面を参照して本発明の実施の形態を説明する。
(第1の実施形態)
図1を参照して、本発明による適応フィルタを適応ノイズキャンセラに適用した場合の実施形態について説明する。この適応フィルタは、参照信号xに対しフィルタ演算を行って出力信号yを得る遅延線タップフィルタにより構成されるフィルタ部11と、希望応答dからフィルタ部11の出力信号yを減算して適応フィルタの出力信号である誤差信号eを得る減算器12と、参照信号xと誤差信号eに基づいてフィルタ部11のフィルタ特性(フィルタ係数)を更新するフィルタ更新部13からなる。以下、フィルタ部11のフィルタ係数を更新する処理を単にフィルタ更新という。
ここで、フィルタ更新部13が図14に示した従来の適応フィルタと異なっている。すなわち、このフィルタ更新部13は、誤差信号eのパワーによって適応モードを切り換え制御する適応モード制御部14と、この適応モード制御部14によって制御される適応モードに従って実際にフィルタ更新演算を行うフィルタ更新演算部15とにより構成される。フィルタ更新演算の詳細については、後述する。
本実施形態では、適応フィルタの適応アルゴリズムとして最も一般的なLMS(Least Mean Squre)法を例として説明するが、適応速度を制御できるアルゴリズムであれば、射影LMS法でもRLS(Recursive Least Squre )法でも、これら以外のアルゴリズムでもよい。
適応ノイズキャンセラでは、例えば図2に示すようにマイクロホンM1を目的信号源の近傍に置き、このマイクロホンM1の出力信号を希望応答dとし、もう一つのマイクロホンM2を妨害雑音源の近傍に置き、このマイクロホンM2の出力信号を参照信号xとして使用する。この場合、マイクロホンM1には目的信号と背景雑音が入力し、マイクロホンM2には妨害雑音と背景雑音が入力する。マイクロホンM2には目的信号がなるべく入力しないようにマイクロホン配置などの設定を行うが、完全に目的信号が入力しないようにすることは困難である。すなわち、マイクロホンM2には実際には妨害雑音と背景雑音のほか、目的信号の成分も一部入力する。
一般に、このような適応ノイズキャンセラでは参照信号xと相関のある成分を希望応答dの中から除去するため、マイクロホンM1,M2間で相関の大きい妨害雑音が除去される。実際には、上述のように参照信号x中に目的信号成分が混入した場合、妨害雑音だけでなく目的信号が除去されて目的信号が歪むことがあるため、目的信号が入力した際に適応を遅くする必要がある。
LMS法による適応フィルタのフィルタ更新演算は、従来では先の式(1)または(2)に基づいて行われるが、本実施形態では適応フィルタの出力信号である誤差信号eのパワー(以下、出力信号パワーという)Pに応じて適応モードを切り換え、次式に基づいてフィルタ更新演算を行う。
=d−W (3)
≦Pのとき
n+1 =W−2μe (4)
>Pのとき
n+1 =W−2μe/P (5)
ここで、Pは予め定められた適応モード切換えのしきい値、Wはフィルタ部11のn回の更新後のフィルタ係数、μはステップサイズ、eは適応フィルタの出力信号(誤差信号)、dは希望応答、Xは参照信号xの最新のサンプルから過去のサンプルをフィルタ部11のタップ点数分だけ並べてベクトルとした参照信号ベクトルである。予め設定される各種定数の値は、例えばフィルタ部11のタップ数を100、ステップサイズμを0.001、しきい値Pを平均出力信号パワーの0.001倍とするが、状況に応じて実験的に決めることが望ましい。
次に、本実施形態におけるフィルタ更新動作について図3を用いて説明する。図3は、フィルタ更新演算式(2)(4)(5)における更新量を表す右辺第2項において、誤差信号eと参照信号Xの積eに掛かる係数(フィルタ係数の更新速度を決める係数)の大きさの変化を図示したもので、この係数の大きさをステップサイズμと区別して便宜上、実効ステップサイズと呼ぶことにする。図3(a)は、従来の改良されたNLMS法におけるフィルタ更新演算式である式(2)において、しきい値P以下で適応停止を行った場合の実効ステップサイズを参照信号パワーPと出力信号パワーPの和に対して示したものであり、(b)は本実施形態に基づくフィルタ更新演算式(4)(5)の場合の実効ステップサイズを出力信号パワーPに対して示している。
NLMS法では、従来の技術でも述べたように、妨害雑音も目的信号もなく参照信号パワーP、または参照信号パワーPと出力信号パワーPの和が低い場合には、背景雑音に対して適応が進み、フィルタ係数が本来除去したい妨害雑音を除去するための値から遠ざかってしまうことを防ぐために、PまたはP+Pがしきい値P以下になったときに適応を停止している。この場合、図3(a)に示すように実効ステップサイズは0になるので、PまたはP+Pがしきい値P以下となるような妨害雑音に対しては、全く適応できないことになる。
これに対し、本実施形態では出力信号パワーPがしきい値P以下の場合は式(4)によりフィルタ更新を行い、出力信号パワーPがしきい値Pより大きい場合は式(5)によりフィルタ更新を行うように適応モードを切り換えている。式(4)はステップサイズ固定のLMS法に基づくフィルタ更新演算式、式(5)はNLMS法に基づくフィルタ更新演算式であり、実効ステップサイズは出力信号パワーPに対して図3(b)に示したように変化する。
すなわち、本実施形態ではP≦Pの場合は実効ステップサイズは2μ×Pなる非零の一定値であり、P>Pの場合は実効ステップサイズは2μ*P/PとなってPに反比例する。なお、P≦Pの場合の実効ステップサイズは厳密に非零の一定値である必要はなく、ほぼ一定値であればよい。
従って、出力信号パワーPがしきい値P以下の場合の式(4)によるフィルタ更新でも、完全に適応を停止しないので、妨害雑音に適応することができ、しかも式(4)はステップサイズ固定のフィルタ更新であるため、レベルの低い背景雑音に適応してしまうことは少ない。一方、出力信号パワーPがしきい値Pより大きい場合の式(5)によるフィルタ更新では、図3(b)において出力信号パワーPがしきい値P以上の部分に示されるように、実効ステップサイズが出力信号パワーPに反比例し適応が遅くなってくるので、目的信号の歪みを低減できる。
次に、図4に示すフローチャートを参照して上述した適応フィルタ処理の流れを説明する。
まず、初期設定でステップサイズμ、適応モード切り換えのしきい値P、フィルタ長Kを設定し、フィルタ係数を例えば全て0、フィルタ更新回数nを0とおく(ステップS11)。
次に、参照信号xの複数のサンプルから参照信号ベクトルXを構成した後、
式(3)により適応フィルタ出力信号(誤差信号)eを計算し、この出力信号eのパワーPを計算する。出力信号パワーPの計算は、例えば信号eの最新の値からフィルタ部11のタップ数分までの過去の値の平均値を用いるようにする(ステップS12)。
次に、ステップS13において出力信号パワーPとしきい値Pを比較し、P≦Pの場合はステップS14に進み、P>Pの場合はステップS15に進む。ステップS14では、式(4)に従ってフィルタ更新を行い、ステップS16に進む。ステップS15では、式(5)に従ってフィルタ更新を行い、ステップS16に進む。ステップS16では、nをインクリメントし、ステップS12に戻る。そして、以上の処理を入力データの数だけ繰り返す。
このように本実施形態では、出力信号パワーPに応じて適応フィルタの適応モードをステップサイズ固定のLMS法とNLMS法との間で切り換え制御することにより、目的信号の歪み低減と背景雑音への適応防止の両立を図ることができる。
本実施形態では、P>Pの場合の実効ステップサイズを式(5)に示されるようにμ*P/Pとして出力信号パワーPに反比例するようにしたが、必ずしも正確に反比例するようにする必要はなく、出力信号パワーPの増加に対して減少させるようにすればよい。また、このようにP>Pの場合の実効ステップサイズを出力信号パワーPの増加のみに依存して減少させる代りに、従来技術におけるフィルタ特性の更新演算式である式(2)に示されるように、参照信号パワーPと出力信号パワーPとの和の増加に対して減少させるようにしてもよい。
(第2の実施形態)
本実施形態では、本発明による適応フィルタを多チャネルに拡張した場合につしいて説明する。
この適応フィルタは、図5に示すように複数チャネル(Mチャネル)の参照信号xに単位フィルタ部21−1〜21−Mでフィルタ演算を施した後、それらの出力信号の和yを加算器26で求めて出力するフィルタ部21と、希望応答dからフィルタ部21の出力信号yを減算して適応フィルタの出力信号である誤差信号eを得る減算器22と、参照信号xと誤差信号eに基づいてフィルタ部21の特性(フィルタ係数)を更新するフィルタ更新部23からなる。
フィルタ更新部23は、誤差信号eのパワーに応じて適応モードを切り換え制御する適応モード制御部24と、この適応モード制御部24によって制御される適応モードに従って実際にフィルタ更新演算を行うフィルタ更新演算部25とにより構成される。
次に、本実施形態の適応フィルタの動作を説明する。チャネル数をM、1チャネル当たりの単位フィルタ部21−1〜21−Mの各々のタップ数をKとする。適応フィルタの適応アルゴリズムとしては、最も一般的なLMS法を例として説明するが、適応速度を制御できるアルゴリズムならば射影LMS法でもRLS法でも何でもよいことは第1の実施形態と同様である。
多チャネルのフィルタ部21は、図6に示すようなタップ付き遅延線とその各タップ出力に係数を乗じる係数乗算器とからなる遅延線タップフィルタをM個並べて単位フィルタ部21−1〜21−Mが構成される。この多チャネルフィルタ部21のフィルタ特性は、mチャネル目の遅延線タップフィルタのk番目の係数wm,k を並べたベクトル
=(w1,1 ,w1,2 ,…,w1,K
…,wM−1,1 ,…,wM−1,K (6)
により表し、参照信号xは図5においてmチャネル目の参照信号の最新のサンプルxm,n からフィルタタップ数分過去までのサンプルを並べた次式により表すことにする。
=(x1,n ,x1,n−1 ,…,x1,n−K+1
…,xM,n ,xM,n−1 ,…,xM,n−K+1 (7)
フィルタ更新演算は上のW,Xを用いると、式(3)(4)(5)と全く同様に、次式に基づいて行われる。
=d−W (8)
≦Pのとき
n+1 =W−2μe (9)
>Pのとき
n+1 =W−2μe/P (10)
ここで、Pは予め定められた適応モード切換えのしきい値、eは適応フィルタの出力信号(誤差信号)、Pは出力信号パワー、μはステップサイズである。予め設定される各種定数の値は、例えば単位フィルタ部11−1〜11−Mの各々のタップ数を100、ステップサイズμを0.001、しきい値Pを平均出力信号パワーの0.001倍とするが、状況に応じて実験的に決めるのが望ましい。
本実施形態においても、第1の実施形態の場合と同様に、適応フィルタの出力信号パワーPにより適応モードをステップサイズ固定のLMSとNLMSとの間で切り換えており、式(9)のステップサイズ固定のLMS法では雑音への適応を防止し、式(10)のNLMS法では目的信号の歪みを低減している。
本実施形態の適応フィルタの処理の流れは第1の実施形態と同様であるので、改めて述べない。
このように本実施形態では、多チャネル入力の適応フィルタにおいて出力信号パワーPに応じて適応フィルタの適応モードをステップサイズ固定のLMS法とNLMS法との間で切り換え制御することにより、目的信号の歪み低減と背景雑音への適応防止の両立を図ることができる。
本実施形態では、P>Pの場合の実効ステップサイズを式(10)に示されるようにμ*P/Pとして出力信号パワーPに反比例するようにしたが、必ずしも正確に反比例するようにする必要はなく、出力信号パワーPの増加に対して減少させるようにすればよい。また、このようにP>Pの場合の実効ステップサイズを出力信号パワーPの増加のみに依存して減少させる代りに、従来技術におけるフィルタ特性の更新演算式である式(2)に示されるように、参照信号パワーPと出力信号パワーPとの和の増加に対して減少させるようにしてもよい。
(第3の実施形態)
次に、第1または第2の実施形態で述べた適応フィルタを、妨害雑音を抑圧して特定の方向から到来する目的信号のみを抽出する適応ビームフォーマに適用した実施形態について図7を用いて説明する。
図7は、適応ビームフォーマの中でも一般的な一般化サイドローブキャンセル(GSC)と呼ばれる構成を示しており、複数チャネル(Mチャネル)の入力信号sから目的信号を除去してM−1チャネルの信号を出力するブロッキングフィルタ41と、第1または第2の実施形態で説明した適応フィルタ42と、Mチャネルの入力信号sの加算を行う加算器43からなる。
M=2の場合はM−1=1であるので、適応フィルタ41として第1の実施形態で説明した1チャネル入力の適応フィルタを用い、M>2の場合は適応フィルタ41として第2の実施形態で説明した多チャネル入力の適応フィルタを用いることになる。ここでは、M>2として説明するが、M=2の場合は単にチャネル数を減らして同様の処理を行うだけである。
GSCは、目的信号を除去するブロッキングフィルタの構成により種々のものがあるが、ここでは最も簡単なGriffith-Jim型GSCを例にとって説明する。本発明の適応フィルタはブロッキングフィルタの方式に依存していないので、このGriffith-Jim型以外のGSCでも適用可能である。
Griffith-Jim型GSCのブロッキングフィルタは、図8に示すようにMチャネルの入力信号の隣り合ったチャネル間の差分をとるフィルタであり、全てのチャネルに目的信号が同位相で入力した場合、目的信号が消去された信号がM−1チャネル分出力される。mチャネル目のブロッキングフィルタの出力信号は、時刻をnとし、mチャネル目の入力信号をsm,n として、
m,n =sm+1,n −sm,n
(m=1,…,M−1) (11)
である。
このブロッキングフィルタ41の出力信号を第2の実施形態で述べた適応フィルタ42に参照信号xとして入力してフィルタ更新を行う場合、更新演算式は、適応フィルタ42のフィルタ係数をW=(w1,1 ,w1,2 ,…,w1,K ,…,wM−1,1 ,…,wM−1,K 、参照信号ベクトルをX=(x1,n ,x1,n−1 ,…,x1,n−K+1 ,…,xM,n ,xM,n−1 ,…,xM,n−K+1 として、
=ΣxM,n (12)
=y−W (13)
≦Pのとき
n+1 =W−2μe (14)
>Pのとき
n+1 =W−2μe/P (15)
と表わされる。
このような適応ビームフォーマを例えば前述したマイクロホンアレイ処理装置に適用して、音声を目的信号として抽出する場合、音声の伝播揺らぎや音源が完全な点音源でないことなどの原因により、ブロッキングフィルタ41のみでは目的信号を完全に除去することは困難である。
これに対し、第1または第2の実施形態で説明した適応フィルタをブロッキングフィルタ41の後段に適応フィルタ42として配置すれば、この適応フィルタは前述したように目的信号が入力した場合の目的信号の歪み低減と、妨害雑音と目的信号がない場合の適応フィルタの背景雑音への適応防止が可能であるため、より性能を向上できる。
すなわち、本実施形態では適応フィルタ42において出力信号パワーPがしきい値P以下の場合は式(14)によりフィルタ更新を行い、出力信号Pがしきい値Pより大きい場合は式(15)によりフィルタ更新を行うように適応モードを切り換える。式(14)はステップサイズ固定のLMS法なので、パワーの低い背景雑音への適応が少なく、また完全に適応を停止しないので、しきい値P以下の妨害雑音に対しては適応ができる。
一方、出力信号パワーが大きい場合の式(15)による更新では、ステップサイズが出力信号パワーに反比例するので目的信号の歪みを低減できる。
次に、図9に示すフローチャートを参照して本実施形態による適応ビームフォーマの処理の流れを説明する。
まず、初期設定でステップサイズμ、適応モード切り換えのしきい値P、チャネル数M、フィルタ長Kを設定し、フィルタの係数を例えば全て0、フィルタ更新の回数nを0とおく(ステップS21)。
次に、入力チャネル間の差分信号xm,n を式(11)により計算する(ステップS22)。
次に、入力信号の和yを式(12)により計算する(ステップS23)。
次に、差分信号xm,n から参照信号ベクトルXを構成し、誤差信号eを式(13)により計算する(ステップS24)。
次に、誤差信号eのパワーPを計算して、ステップS25でPとしきい値Pを比較し、P≦Pの場合はステップS26に進み、P>Pの場合はステップS27に進む。
ステップS26では、式(12)に従ってフィルタを更新し、ステップS28に進む。ステップS27では、式(13)に従ってフィルタを更新し、ステップS28に進む。次にステップS28では、nをインクリメントし、ステップS22に戻る。そして、以上の処理を入力データの数だけ繰り返す。
このように本実施形態では適応ビームフォーマ、特にGSCに本発明の適応フィルタを適用することにより、音声などの非定常な信号が入力する場合の雑音抑圧性能を大幅に向上することが可能となる。
本実施形態では、P>Pの場合の実効ステップサイズを式(15)に示されるように、μ*P/Pとして出力信号パワーPに反比例させたが、必ずしも正確に反比例するようにする必要はなく、出力信号パワーPの増加に対して減少させるようにすればよい。また、このようにP>Pの場合の実効ステップサイズを出力信号パワーPの増加のみに依存して減少させる代りに、従来技術におけるフィルタ特性の更新演算式である式(2)に示されるように、参照信号パワーPと出力信号パワーPとの和の増加に対して減少させるようにしてもよい。
(第4の実施形態)
次に、拘束条件の下でパワーを最小化するフロスト型ビームフォーマに本発明の適応フィルタを適用した場合の実施形態について述べる。
本実施形態のフロスト型ビームフォーマは図10に示す構成であり、複数チャネル(Mチャネル)の入力信号xに単位フィルタ部51−1〜51−Mでフィルタ演算を施した後、それらの出力信号の和yを加算器56で求めて出力するフィルタ部51と、出力信号yに基づいてフィルタ部51の特性(フィルタ係数)を更新するフィルタ更新部53からなる。フィルタ部51は、例えば図6に示した多チャネルの遅延線タップフィルタによって構成される。チャネル数はM、タップ数はKとする。フィルタ更新部53は、適応モード制御部54とフィルタ更新演算部55および拘束条件設定部57からなる。
フィルタ更新部53では、処理に先立って拘束条件設定部55においてフィルタ部51の特性変化範囲を拘束するための拘束条件を設定し、この拘束条件設定部57で設定された条件の下で、入力信号xに基づいてフィルタ更新演算部55でフィルタ部51の更新演算を行う。このフィルタ更新演算の際、出力信号yのパワーに基づき、適応モード制御部54により適応速度の制御を行う。
ここで、フロスト型ビームフォーマの処理に関して説明する。フロスト型ビームフォーマは、例えば、文献4:O.L.Frost,III,“An Algorithm for Linearly Constrained Adaptive Array Processing ”,Proceeding of The IEEE,Vol.60,No.8,pp.926-935(1972) に詳述されているように、LMSアルゴリズムにより次式に従ってフィルタ更新を行う。
=W (16)
n+1 =P[W−μy]+F (17)
ここでnは時刻、W=(w1,1 ,w1,2 ,…,w1,K ,…,wM−1,1 ,…,wM−1,K はフィルタ係数、X=(x1,n ,x1,n−1 ,…,x1,n−K+1 ,…,xM,n ,xM,n−1 ,…,xM,n−K+1 は参照信号ベクトル、Pは拘束条件により定まる部分空間への射影行列、Fは該部分空間から拘束条件を満たす空間への平行移動ベクトル、μはステップサイズ、yはフィルタ出力である。
射影行列Pと平行移動ベクトルFは、
P=I−A(AA)−1 (18)
F=A(AA)−1G (19)
により、適応処理の前に求めておく。ここで、AとGはフィルタの拘束条件を定める行列であり、例えば、拘束条件の数をJとし、AはA=[a,a,…,a,…,a]、GはJ/2+1番めのみ値が1で、他の要素は0のベクトルとする。aは例えば、
i=K*(m−1)+(K/2+1)−(J/2+1)+j (20)
(m=1,2,…,M)
で計算される番号iの要素の値が1で、他は0のベクトルを用いる。なお、上式において割り算は切り下げとする。
本実施形態においては、このフロスト型ビームフォーマの更新演算式(16)(17)を変更し、次式に基づいてフィルタ更新を行う。
=W (21)
≦Pのとき
n+1 =P[W−μy]+F (22)
>Pのとき
n+1 =P[W−μy/P]+F (23)
ここで、Pは出力信号yのパワー、Pは予め定めたしきい値である。先の実施形態と同様、出力信号パワーPに応じて適応速度を切り換え制御するようにしており、出力信号パワーPがしきい値以下の場合はステップサイズ固定のLMS法、大きければNLMS法に切り換える。
これによりGSCの場合と同様に、目的信号が入力した場合には式(21)のNLMS法により適応速度を遅くして目的信号の歪みを抑え、目的信号も妨害雑音も入力していない場合には式(21)によるステップサイズ固定のLMS法により、背景雑音への適応を少なくしている。
次に、図11を参照して本実施形態の処理の流れを説明する。
まず、初期設定でステップサイズμ、適応モード切り換えしきい値P、チャネル数M、フィルタ長Kを設定し、フィルタの値を例えば全て0、フィルタ更新の回数nを0とおく(ステップS31)。
次に、拘束条件を定める行列A,Gを構成し、式(18)(19)により行列P,Fを求める(ステップS32)。
次に、参照信号ベクトルXを構成し、フィルタ出力yを式(21)により計算する(ステップS33)。
次に、フィルタ出力信号yのパワーPを計算して、ステップS34でPとしきい値Pとを比較し、P≦Pの場合はステップS35に進み、P>Pの場合はステップS36に進む。
ステップS35では、式(22)に従ってフィルタを更新し、ステップS37に進む。ステップS36では、式(23)に従ってフィルタを更新し、ステップS37に進む。次にステップS37では、nをインクリメントし、ステップS33に戻る。そして、以上の処理を入力データの数だけ繰り返す。
以上に述べたように、フロスト型ビームフォーマに本発明の適応フィルタを適用することにより、音声などの非定常な信号を入力する場合の雑音抑圧性能を大幅に向上することが可能となった。
(第5の実施形態)
次に、第5の実施形態として、本発明による適応フィルタを第1の適応フィルタとし、これと適応速度がより大きい第2の適応フィルタを組み合わせて、第2の適応フィルタの出力信号を利用して第1の適応フィルタの適応速度を制御することにより、第1の適応フィルタの性能をさらに向上させた信号処理装置の実施形態について説明する。これは第1〜第4の実施形態で説明した本発明の適応フィルタの収束速度をさらに改善するものであり、特に定常的な信号が妨害雑音となっている場合のノイズキャンセラやビームフォーマにおいて、目的信号の歪みの低減と収束速度向上を両立することができる。
図12に、本実施形態による信号処理装置の構成を示す。この信号処理装置は第1の適応フィルタ1と第2の適応フィルタ2により構成され、両者に同じ入力信号が与えられる。
第1の適応フィルタ1は第1〜第4の実施形態のいずれかの適応フィルタと同じ構成であり、参照信号xに対してフィルタ演算を行って出力信号y1を得るフィルタ部11−1と、希望応答dからフィルタ部11−1の出力信号を減算して誤差信号e1を得る減算器12−1と、参照信号xと誤差信号e1に基づいてフィルタ部11−1のフィルタ係数を更新する適応モード制御部14−1およびフィルタ更新演算部15−1とからなる。この第1の適応フィルタ1の出力信号が信号処理装置の出力となる。
第2の適応フィルタ2は、一般的なNLMSアルゴリズムによる適応フィルタの構成であり、参照信号xに対してフィルタ演算を行って出力信号y2を得るフィルタ部11−2と、希望応答dからフィルタ部11−2の出力信号y2を減算して誤差信号e2を得る減算器12−2と、参照信号xと誤差信号eに基づいてフィルタ更新の演算を行うフィルタ更新演算部15−2とからなる。
ここで、第2の適応フィルタ2は第1の適応フィルタ1よりも適応速度が高く設定されており、第2の適応フィルタ2の出力信号パワーに基づき、適応モード制御部14−1において第1の適応フィルタ1の適応速度を制御するようにする。適応速度は、例えばステップサイズの値を大きくしたり、フィルタのタップ数を少なくすることにより大きい値に設定できる。全体の処理は、時刻をn、第1の適応フィルタ1のフィルタ係数をW1、第2の適応フィルタ2のフィルタ係数をW2とおき、次式により行う。
e1=d−W1 (24)
e2=d−W2 (25)
W2n+1 =W2−2μe2/P (26)
P1≦Pのとき
W1n+1 =W1−2μe1 (27)
P1>Pのとき
W1n+1 =W1−2μe1
/P1*(P1/P2)α (28)
ここで、e1,e2はそれぞれ適応フィルタ1,2の出力信号、Pは参照信号パワー、dは希望応答、μ,μはそれぞれ適応フィルタ1,2のステップサイズである。μの値はP*μよりも大きい値とし、例えば、P*μ=0.000003、μ=0.2とする。式(28)におけるαは、適応速度の加速係数であり、例えば、1.5とする。
式(26)の適応フィルタ2の処理は一般的なNLMS法におけるフィルタ更新演算式であり、ステップサイズμを適応フィルタ1のステップサイズμよりも大きくしているため適応が速い。従って、第2の適応フィルタ2では目的信号の歪みが大きいが、素早く妨害雑音を除去することができる。
一方、第1の適応フィルタ1は逆に適応速度を比較的遅くしているため、妨害雑音が入力した場合、初期収束する前の時点で除去されずに出力された妨害雑音のパワーによって適応速度が遅くなる。
そこで、式(28)に示すように、第1の適応フィルタ1と第2の適応フィルタ2の出力信号パワーの比P1/P2に応じて、第1の適応フィルタ1のフィルタ更新の際の変化を大きくするようにする。妨害雑音が入力した際、収束が速い第2の適応フィルタ2の出力信号パワーP2は、第1の適応フィルタ1の出力信号パワーP1よりも小さいので、パワー比P1/P2は1より大きくなり、第1の低既往フィルタ1の適応が速くなる。一方、目的信号が入力した場合、パワー比P1/P2は1に近い値となるため、第1の適応フィルタ1の適応速度は遅いままであり、入力信号の歪み低減は保持される。
次に、図13を参照して本実施形態の処理の流れを説明する。
まず、初期設定でステップサイズμ,μ、適応モード切り換えのしきい値P、適応速度の加速係数α、フィルタ長Kを設定し、フィルタ係数の値を例えば全て0、時刻nを0とおく(ステップS41)。
次に、参照信号xから参照信号ベクトルXを構成して、誤差信号e1,e2を式(24)(25)により計算し、さらに誤差信号e1,e2のパワーP1,P2を計算する(ステップS42)。
次に、参照信号ベクトルXから参照信号パワーPを求め、式(26)に従い、適応フィルタ2の更新を行う(ステップS43)。
次に、ステップS44において誤差信号e1のパワーP1としきい値Pを比較し、P1≦Pの場合はステップS45に進み、P1>Pの場合はステップS46に進む。
ステップS45では、式(27)に従ってフィルタを更新し、ステップS47に進む。ステップS46では、式(28)に従ってフィルタを更新し、ステップS47に進む。ステップS47では、nをインクリメントし、ステップS42に戻る。そして、以上の処理を入力データの数だけ繰り返す。
このように本実施形態では、出力信号パワーに基づいて適応速度を切換え制御するようにした本発明の適応フィルタを第1の適応フィルタ1として用い、これより適応速度を速くした第2の適応フィルタ2の出力信号に応じて第1の適応フィルタ1の適応速度を変化させることにより、目的信号の歪みを低減したまま収束速度を向上させることが可能となる。
これにより、適応ノイズキャンセラや適応ビームフォーマにおいて、特に定常的な妨害雑音が入力した場合の妨害雑音抑圧性能を大幅に向上させることができる。本実施形態における適応ビームフォーマは、第3の実施形態で述べたGSC型、第4の実施形態で述べたフロスト型いずれにも適用可能である。
なお、本実施形態では参照信号が1チャネルの場合について述べたが、第2の実施形態と同様に複数チャネルの参照信号を扱うように構成することも可能である。
本発明の第1の実施形態に係る適応フィルタの構成を示すブロック図 適応ノイズキャンセラの原理を説明するための図 本発明と従来のフィルタ更新動作を比較して説明するための図 第1の実施形態における適応フィルタ処理の流れを示すフローチャート 本発明の第2の実施形態に係る多チャネル化した適応フィルタの構成を示すブロック図 第2の実施形態における多チャネルフィルタ部の説明図 本発明の第3の実施形態に係る信号処理装置である適応ビームフォーマの構成を示すブロック図 第3の実施形態で用いるブロッキングフィルタの一例の構成を示す図 第3の実施形態における処理の流れを示すフローチャート 本発明の第4の実施形態に係る信号処理装置であるフロスト型ビームフォーマの構成を示すブロック図 第4の実施形態における処理の流れを示すフローチャート 本発明の第5の実施形態に係る信号処理装置の構成を示すブロック図 第5の実施形態における処理の流れを示すフローチャート 従来の適応フィルタの構成を示すブロック図
符号の説明
1…第1の適応フィルタ
2…第2の適応フィルタ
11,11−1,11−2…フィルタ部
12,12−1,12−2…減算器
13…フィルタ更新部
14,14−1…適応モード制御部
15,15−1,15−2…フィルタ更新演算部
21…フィルタ部
21−1〜21−M…単位フィルタ部
22…減算器
23…フィルタ更新部
24…適応モード制御部
25…フィルタ更新演算部
26…加算器
41…ブロッキングフィルタ
42…適応フィルタ
51…フィルタ部
51−1〜51−M…単位フィルタ部
52…減算器
53…フィルタ更新部
54…適応モード制御部
55…フィルタ更新演算部
56…加算器
57…拘束条件設定部

Claims (5)

  1. 第1の適応フィルタと、この第1の適応フィルタより適応速度の大きい第2の適応フィルタとを具備し、第1の適応フィルタは第2の適応フィルタの出力信号に基づいてフィルタの更新が制御されることを特徴とする信号処理装置。
  2. 前記第1の適応フィルタは、その出力信号パワーを予め定められたしきい値と比較し、出力信号パワーがしきい値以下のときは適応速度が非零の一定値となり、出力信号パワーがしきい値を越えたときは該適応速度が第1および第2の適応フィルタの出力信号パワーに応じて変化するように、フィルタの更新が制御されることを特徴とする請求項1に記載の信号処理装置。
  3. 前記第1の適応フィルタは、参照信号に対してフィルタ演算を行って出力信号を得るフィルタ部と、希望応答から前記フィルタ部の出力信号を減算して誤差信号を得る減算器と、前記参照信号と前記誤差信号に基づいてフィルタ部のフィルタ係数を更新するフィルタ更新演算部と、前記フィルタ部の出力信号パワーを予め定められたしきい値と比較し、前記出力信号パワーがしきい値以下のときは適応速度が非零の一定値となり、前記出力信号パワーがしきい値を越えたときは該適応速度が前記出力信号パワーおよび第2の適応フィルタの出力信号パワーに応じて変化するように前記フィルタ係数の更新を制御する適応モード制御部とを有することを特徴とする請求項2記載の信号処理装置。
  4. 前記第1の適応フィルタは、その出力信号パワーを予め定められたしきい値と比較し、出力信号パワーがしきい値以下のときはフィルタの更新速度を決める係数を非零の一定値とし、出力信号パワーがしきい値を越えたときは該係数を第1および第2の適応フィルタの出力信号パワーの比の増加に対して減少させることを特徴とする請求項1に記載の信号処理装置。
  5. 前記第1の適応フィルタは、参照信号に対してフィルタ演算を行って出力信号を得るフィルタ部と、希望応答から前記フィルタ部の出力信号を減算して誤差信号を得る減算器と、前記参照信号と前記誤差信号に基づいてフィルタ部のフィルタ係数を更新するフィルタ更新演算部と、前記フィルタ部の出力信号パワーを予め定められたしきい値と比較し、前記出力信号パワーがしきい値以下のときは適応速度が非零の一定値となり、出力信号パワーがしきい値を越えたときは該適応速度が前記出力信号パワーおよび第2の適応フィルタの出力信号パワーに応じて変化するように前記フィルタ係数の更新を制御する適応モード制御部とを有することを特徴とする請求項4記載の信号処理装置。
JP2006319008A 2006-11-27 2006-11-27 適応フィルタを用いた信号処理装置 Expired - Fee Related JP4928918B2 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2006319008A JP4928918B2 (ja) 2006-11-27 2006-11-27 適応フィルタを用いた信号処理装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2006319008A JP4928918B2 (ja) 2006-11-27 2006-11-27 適応フィルタを用いた信号処理装置

Related Parent Applications (1)

Application Number Title Priority Date Filing Date
JP17420297A Division JP4068182B2 (ja) 1997-06-30 1997-06-30 適応フィルタ

Publications (2)

Publication Number Publication Date
JP2007124678A true JP2007124678A (ja) 2007-05-17
JP4928918B2 JP4928918B2 (ja) 2012-05-09

Family

ID=38147922

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2006319008A Expired - Fee Related JP4928918B2 (ja) 2006-11-27 2006-11-27 適応フィルタを用いた信号処理装置

Country Status (1)

Country Link
JP (1) JP4928918B2 (ja)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2013114797A1 (ja) 2012-02-02 2013-08-08 旭化成エレクトロニクス株式会社 物理量計測装置及び物理量計測方法
JP2020502562A (ja) * 2016-11-23 2020-01-23 テレフオンアクチーボラゲット エルエム エリクソン(パブル) 相関分離フィルタの適応制御のための方法および装置

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH04271622A (ja) * 1991-02-27 1992-09-28 Nec Corp エコーキャンセラ
JPH0918291A (ja) * 1995-06-30 1997-01-17 Nec Corp 雑音消去方法及び雑音消去装置

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH04271622A (ja) * 1991-02-27 1992-09-28 Nec Corp エコーキャンセラ
JPH0918291A (ja) * 1995-06-30 1997-01-17 Nec Corp 雑音消去方法及び雑音消去装置

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2013114797A1 (ja) 2012-02-02 2013-08-08 旭化成エレクトロニクス株式会社 物理量計測装置及び物理量計測方法
CN104081161A (zh) * 2012-02-02 2014-10-01 旭化成微电子株式会社 物理量测量装置以及物理量测量方法
US9740661B2 (en) 2012-02-02 2017-08-22 Asahi Kasei Microdevices Corporation Physical quantity measuring apparatus and physical quantity measuring method
JP2020502562A (ja) * 2016-11-23 2020-01-23 テレフオンアクチーボラゲット エルエム エリクソン(パブル) 相関分離フィルタの適応制御のための方法および装置
JP2021101242A (ja) * 2016-11-23 2021-07-08 テレフオンアクチーボラゲット エルエム エリクソン(パブル) 相関分離フィルタの適応制御のための方法および装置
US11501785B2 (en) 2016-11-23 2022-11-15 Telefonaktiebolaget Lm Ericsson (Publ) Method and apparatus for adaptive control of decorrelation filters
JP7201721B2 (ja) 2016-11-23 2023-01-10 テレフオンアクチーボラゲット エルエム エリクソン(パブル) 相関分離フィルタの適応制御のための方法および装置
US11942098B2 (en) 2016-11-23 2024-03-26 Telefonaktiebolaget Lm Ericsson (Publ) Method and apparatus for adaptive control of decorrelation filters

Also Published As

Publication number Publication date
JP4928918B2 (ja) 2012-05-09

Similar Documents

Publication Publication Date Title
JP5451876B2 (ja) 音響マルチチャネル相殺
US7773743B2 (en) Integration of a microphone array with acoustic echo cancellation and residual echo suppression
JP4957810B2 (ja) 音処理装置、音処理方法及び音処理プログラム
CN107483761B (zh) 一种回波抑制方法及装置
US20130322655A1 (en) Method and device for microphone selection
JP5662232B2 (ja) エコー消去装置、その方法及びプログラム
WO2007123048A1 (ja) 適応アレイ制御装置、方法、プログラム、及びこれを利用した適応アレイ処理装置、方法、プログラム
JP4928918B2 (ja) 適応フィルタを用いた信号処理装置
JP4068182B2 (ja) 適応フィルタ
JP4879195B2 (ja) 騒音低減装置
JP2005513874A (ja) スペクトルエコーのテイル予測器を有するエコーキャンセラ
JP2008005094A (ja) エコー抑圧方法、装置、エコー抑圧プログラム、記録媒体
JP3982823B2 (ja) 音声処理装置、音声処理方法および音声処理プログラム
KR100545832B1 (ko) 간섭신호에 강인한 음향 반향 제거장치
JP2004349796A (ja) 音響エコー消去方法、その装置、プログラム及びその記録媒体
JPH0667693A (ja) 雑音低減装置
JP3609611B2 (ja) エコー除去方法およびエコーキャンセラ
JP2015201787A (ja) エコー消去装置、その方法及びプログラム
JP6356087B2 (ja) エコー消去装置、その方法及びプログラム
JP4324036B2 (ja) エコーキャンセル装置
US20230137830A1 (en) Wideband adaptation of echo path changes in an acoustic echo canceller
JP2006014359A (ja) エコー抑制処理システム
JP2000165300A (ja) サブバンド適応フィルタを用いた未知システムの同定方法および装置
JPH08250981A (ja) フィルタ係数の推定装置
JP2001077730A (ja) 適応フィルタの係数推定装置

Legal Events

Date Code Title Description
A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20091215

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20100215

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20101207

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20110131

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20110802

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20110930

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20120117

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20120213

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20150217

Year of fee payment: 3

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20150217

Year of fee payment: 3

LAPS Cancellation because of no payment of annual fees