JP3393738B2 - 適応信号処理方法 - Google Patents
適応信号処理方法Info
- Publication number
- JP3393738B2 JP3393738B2 JP16025995A JP16025995A JP3393738B2 JP 3393738 B2 JP3393738 B2 JP 3393738B2 JP 16025995 A JP16025995 A JP 16025995A JP 16025995 A JP16025995 A JP 16025995A JP 3393738 B2 JP3393738 B2 JP 3393738B2
- Authority
- JP
- Japan
- Prior art keywords
- adaptive
- signal
- stage
- time
- processing method
- 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
Description
【０００１】
【産業上の利用分野】本発明は、１個または複数の信号
から評価点における信号を合成するように適応ディジタ
ルフィルタの特性を定め、これにより音響系または振動
系のインパルス応答を解析するようにした適応信号処理
方法。
から評価点における信号を合成するように適応ディジタ
ルフィルタの特性を定め、これにより音響系または振動
系のインパルス応答を解析するようにした適応信号処理
方法。
【０００２】
【従来の技術】従来より、音波の空間伝達特性の推定、
機器等からの騒音振動の測定やその低減のために、適応
信号処理を用いた方法が数多く提案されている。
機器等からの騒音振動の測定やその低減のために、適応
信号処理を用いた方法が数多く提案されている。
【０００３】そのうち、音源または振動源の影響の度合
を調べる寄与診断方法等の従来例として、特開平５−２
６７２２号が挙げられる。
を調べる寄与診断方法等の従来例として、特開平５−２
６７２２号が挙げられる。
【０００４】これは 機器等の周辺の複数点で検出され
る騒音、もしくは振動の各信号をそれぞれ適応フィルタ
に入力し、その適応フィルタの出力信号の和と、騒音も
しくは振動を評価する評価点で検出される騒音もしくは
振動の信号との残差が最小になるように、前記適応フィ
ルタの係数を更新し、一定値に収束した時の各適応フィ
ルタの出力信号を用いて、評価点における騒音もしくは
振動に対して音源もしくは振動源が寄与している度合を
解析するものである。
る騒音、もしくは振動の各信号をそれぞれ適応フィルタ
に入力し、その適応フィルタの出力信号の和と、騒音も
しくは振動を評価する評価点で検出される騒音もしくは
振動の信号との残差が最小になるように、前記適応フィ
ルタの係数を更新し、一定値に収束した時の各適応フィ
ルタの出力信号を用いて、評価点における騒音もしくは
振動に対して音源もしくは振動源が寄与している度合を
解析するものである。
【０００５】また、音波の空間伝達特性の推定では、信
号発生源からの信号を適応フィルタに入力し、前記適応
フィルタの出力信号と、前記信号発生源からの信号を音
波として空間に放射し、設置した測定点で検出される騒
音情報の信号との差の残差信号が最小になるように、適
応フィルタの係数を更新し、そのフィルタ係数を収束さ
せ、収束した適応フィルタ係数を用いて測定点までの空
間の伝達特性を推定していた。
号発生源からの信号を適応フィルタに入力し、前記適応
フィルタの出力信号と、前記信号発生源からの信号を音
波として空間に放射し、設置した測定点で検出される騒
音情報の信号との差の残差信号が最小になるように、適
応フィルタの係数を更新し、そのフィルタ係数を収束さ
せ、収束した適応フィルタ係数を用いて測定点までの空
間の伝達特性を推定していた。
【０００６】
【発明が解決しようとする課題】上記した従来の適応信
号処理方法では、高残響場等の残響時間の長い環境での
場合、インパルス応答が精度よく推定できないため、適
応フィルタの出力信号の和と、騒音もしくは振動の評価
点で検出される騒音もしくは振動の信号との残差が適応
し切れずに、多く残ってしまうという問題がある。
号処理方法では、高残響場等の残響時間の長い環境での
場合、インパルス応答が精度よく推定できないため、適
応フィルタの出力信号の和と、騒音もしくは振動の評価
点で検出される騒音もしくは振動の信号との残差が適応
し切れずに、多く残ってしまうという問題がある。
【０００７】さらに、直接波と反射波の時間間隔が長い
場合には、反射波まで表現しようとすると、適応フィル
タのタップ数を多くとる必要があった。
場合には、反射波まで表現しようとすると、適応フィル
タのタップ数を多くとる必要があった。
【０００８】本発明の目的は、高残響場等での騒音・振
動の測定や音源・振動源の寄与診断等を正確に行え、構
成の簡単な適応信号処理方法を提供することにある。
動の測定や音源・振動源の寄与診断等を正確に行え、構
成の簡単な適応信号処理方法を提供することにある。
【０００９】
【課題を解決するための手段】本発明は、第１の１個又
は複数個の時系列信号を１個又は複数個の１段目適応フ
ィルタの各々に入力し、該適応フィルタの出力信号の和
と前記第１の時系列信号とは別の第２の時系列信号との
差を表す第１の残差信号が最小になるように、前記１段
目適応フィルタの係数を更新するとともに、前記第１の
時系列信号を１個又は複数個の２段目適応フィルタの各
々に入力し、該適応フィルタの出力信号の和と前記第１
の残差信号との差を表す第２の残差信号が最小になるよ
うに前記２段目適応フィルタを更新することを特徴とす
る適応信号処理方法を開示する。
は複数個の時系列信号を１個又は複数個の１段目適応フ
ィルタの各々に入力し、該適応フィルタの出力信号の和
と前記第１の時系列信号とは別の第２の時系列信号との
差を表す第１の残差信号が最小になるように、前記１段
目適応フィルタの係数を更新するとともに、前記第１の
時系列信号を１個又は複数個の２段目適応フィルタの各
々に入力し、該適応フィルタの出力信号の和と前記第１
の残差信号との差を表す第２の残差信号が最小になるよ
うに前記２段目適応フィルタを更新することを特徴とす
る適応信号処理方法を開示する。
【００１０】また本発明は、前記１段目及び２段目の適
応フィルタに入力される前記第１の時系列信号を、その
遅延時間が可変な遅延手段により遅延できるようにした
ことを特徴とする適応信号処理方法を開示する。
応フィルタに入力される前記第１の時系列信号を、その
遅延時間が可変な遅延手段により遅延できるようにした
ことを特徴とする適応信号処理方法を開示する。
【００１１】
【作用】適応フィルタを２段設け、さらに、遅延手段に
より１段目と２段目適応フィルタとの間に時間差を設
け、遅延時間を調節することにより、１段目と２段目適
応フィルタをつなぎ合わせた系の特性でより正確に対象
とする系を近似でき、また、直接波と反射波の間隔がか
なり長い系の場合には、直接波に１段目適応フィルタ
を、反射波に２段目適応フィルタを用いれば、より少な
いタップ数でその系を近似できる。
より１段目と２段目適応フィルタとの間に時間差を設
け、遅延時間を調節することにより、１段目と２段目適
応フィルタをつなぎ合わせた系の特性でより正確に対象
とする系を近似でき、また、直接波と反射波の間隔がか
なり長い系の場合には、直接波に１段目適応フィルタ
を、反射波に２段目適応フィルタを用いれば、より少な
いタップ数でその系を近似できる。
【００１２】又、上記の方法を用いることで、測定信号
の評価点への寄与率や、与えられた騒音または振動伝達
経路の伝達特性をより容易にかつ正確に検出または推定
できる。
の評価点への寄与率や、与えられた騒音または振動伝達
経路の伝達特性をより容易にかつ正確に検出または推定
できる。
【００１３】又、上記推定した伝達特性を利用して評価
点に於ける騒音または振動を効果的に抑圧することが出
来、更に上記検出した寄与率を用いて評価点に生じる騒
音または振動を模擬的に生成することが出来る。
点に於ける騒音または振動を効果的に抑圧することが出
来、更に上記検出した寄与率を用いて評価点に生じる騒
音または振動を模擬的に生成することが出来る。
【００１４】
【実施例】以下、本発明の実施例を図面に基づいて説明
する。図２は、本発明を寄与率解析システムに適用した
時の全体構成を示す図である。同図に示した寄与率解析
システムは、機器１から発生された騒音もしくは振動の
情報を得るマイクロフォン等のセンサ２ａ、２ｂ、２ｃ
と、機器１より一定距離離れた評価点（人間の耳元付
近）に設置したマイクロフォン３と、センサ２ａ、２
ｂ、２ｃで検出された信号がマイクロフォン３に対して
どのくらい寄与しているかを解析する寄与率解析装置４
と、表示部５から構成されている。
する。図２は、本発明を寄与率解析システムに適用した
時の全体構成を示す図である。同図に示した寄与率解析
システムは、機器１から発生された騒音もしくは振動の
情報を得るマイクロフォン等のセンサ２ａ、２ｂ、２ｃ
と、機器１より一定距離離れた評価点（人間の耳元付
近）に設置したマイクロフォン３と、センサ２ａ、２
ｂ、２ｃで検出された信号がマイクロフォン３に対して
どのくらい寄与しているかを解析する寄与率解析装置４
と、表示部５から構成されている。
【００１５】図１は上記の寄与率解析装置４の一実施例
を示すブロック図で、前記センサ２ａ、２ｂ、２ｃで得
られた騒音もしくは振動の信号およびマイクロフォン３
で得られた騒音信号は、それぞれローパスフィルタ４０
１〜４０４で帯域制限を受け、さらにＡ／Ｄ変換器４１
１〜４１４でディジタル信号に変換される。センサ２
ａ、２ｂ、２ｃに基づくディジタル信号は、それぞれ１
段目の適応ディジタルフィルタ（ＡＤＦ）４３１〜４３
３に入力され、ここでフィルタ係数との畳み込み演算が
施されて出力信号ｙ１、ｙ２、ｙ３が生成される。この
とき、必要に応じて遅延回路４２１〜４２３によりディ
ジタル信号を遅延する。次に適応ディジタルフィルタ４
３１〜４３３の出力ｙ１〜ｙ３は加算回路４５１で加算
されてＳ１とされ、これと評価点に設置されたマイクロ
フォンからの騒音信号をディジタル化したディジタル信
号ｄとの差信号ｅ１が減算回路４６で算出される。適応
アルゴリズム回路４４１〜４４３は、上記差信号ｅ１が
最小となるように適応ディジタルフィルタ４３１〜４３
３の係数を更新していく。
を示すブロック図で、前記センサ２ａ、２ｂ、２ｃで得
られた騒音もしくは振動の信号およびマイクロフォン３
で得られた騒音信号は、それぞれローパスフィルタ４０
１〜４０４で帯域制限を受け、さらにＡ／Ｄ変換器４１
１〜４１４でディジタル信号に変換される。センサ２
ａ、２ｂ、２ｃに基づくディジタル信号は、それぞれ１
段目の適応ディジタルフィルタ（ＡＤＦ）４３１〜４３
３に入力され、ここでフィルタ係数との畳み込み演算が
施されて出力信号ｙ１、ｙ２、ｙ３が生成される。この
とき、必要に応じて遅延回路４２１〜４２３によりディ
ジタル信号を遅延する。次に適応ディジタルフィルタ４
３１〜４３３の出力ｙ１〜ｙ３は加算回路４５１で加算
されてＳ１とされ、これと評価点に設置されたマイクロ
フォンからの騒音信号をディジタル化したディジタル信
号ｄとの差信号ｅ１が減算回路４６で算出される。適応
アルゴリズム回路４４１〜４４３は、上記差信号ｅ１が
最小となるように適応ディジタルフィルタ４３１〜４３
３の係数を更新していく。
【００１６】さらに、前記センサ２ａ、２ｂ、２ｃから
得られたディジタル信号は、遅延回路４２４〜４２６で
必要に応じて一定時間遅延された後、それぞれ２段目の
適応ディジタルフィルタ（ＡＤＦ）４３４〜４３６に入
力され、ここでフィルタ係数との畳み込み演算が施され
て出力信号ｙ４、ｙ５、ｙ６が生成される。ここで遅延
回路４２４〜４２６の遅延時間は、遅延設定回路４７１
により１段目の遅延回路４２１〜４２３の遅延時間に連
動して設定される。次に適応ディジタルフィルタ４３４
〜４３６の出力ｙ４〜ｙ６は加算回路４５２で加算され
てｓ２とされ、これと１段目の適応ディジタルフィルタ
４３１〜４３３で処理しきれなかった前記の差信号ｅ１
との差信号ｅ２が減算回路４６２で算出される。適応ア
ルゴリズム回路４４４〜４４６は、上記差信号ｅ２が最
小となるように、適応ディジタルフィルタ４３４〜４３
６の係数を更新していく。なお、残差判定回路７２は、
１段目の適応ディジタルフィルタに基づく差信号ｅ１の
大きさに対するしきい値を設定し、差信号ｅ１がそのし
きい値より大きければ、２段目の適応ディジタルフィル
タを用い、小さければ用いないようにするものである。
得られたディジタル信号は、遅延回路４２４〜４２６で
必要に応じて一定時間遅延された後、それぞれ２段目の
適応ディジタルフィルタ（ＡＤＦ）４３４〜４３６に入
力され、ここでフィルタ係数との畳み込み演算が施され
て出力信号ｙ４、ｙ５、ｙ６が生成される。ここで遅延
回路４２４〜４２６の遅延時間は、遅延設定回路４７１
により１段目の遅延回路４２１〜４２３の遅延時間に連
動して設定される。次に適応ディジタルフィルタ４３４
〜４３６の出力ｙ４〜ｙ６は加算回路４５２で加算され
てｓ２とされ、これと１段目の適応ディジタルフィルタ
４３１〜４３３で処理しきれなかった前記の差信号ｅ１
との差信号ｅ２が減算回路４６２で算出される。適応ア
ルゴリズム回路４４４〜４４６は、上記差信号ｅ２が最
小となるように、適応ディジタルフィルタ４３４〜４３
６の係数を更新していく。なお、残差判定回路７２は、
１段目の適応ディジタルフィルタに基づく差信号ｅ１の
大きさに対するしきい値を設定し、差信号ｅ１がそのし
きい値より大きければ、２段目の適応ディジタルフィル
タを用い、小さければ用いないようにするものである。
【００１７】以上の動作をくり返し、１段目および２段
目の適応ディジタルフィルタ系の係数が一定値に収束し
たとき、各適応ディジタルフィルタ４３１〜４３６の出
力ｙ１〜ｙ６から、評価点に対する各センサ設置位置か
らの騒音もしくは振動の寄与率を解析する。また、解析
された寄与率、各適応ディジタルフィルタの係数等は表
示部５に出力し表示する。
目の適応ディジタルフィルタ系の係数が一定値に収束し
たとき、各適応ディジタルフィルタ４３１〜４３６の出
力ｙ１〜ｙ６から、評価点に対する各センサ設置位置か
らの騒音もしくは振動の寄与率を解析する。また、解析
された寄与率、各適応ディジタルフィルタの係数等は表
示部５に出力し表示する。
【００１８】寄与率の解析方法としては、前記センサ２
ａ、２ｂ、２ｃに対応する適応ディジタルフィルタ出力
信号のパワーを、個々の適応ディジタルフィルタ４３１
〜４３６の出力信号ｙ１〜ｙ６及び差信号ｅ２それぞれ
のパワーの和で除した値を、前記各センサに関する寄与
率とする。例えば、前記センサ２ａの寄与率は、次式で
表わされる。
ａ、２ｂ、２ｃに対応する適応ディジタルフィルタ出力
信号のパワーを、個々の適応ディジタルフィルタ４３１
〜４３６の出力信号ｙ１〜ｙ６及び差信号ｅ２それぞれ
のパワーの和で除した値を、前記各センサに関する寄与
率とする。例えば、前記センサ２ａの寄与率は、次式で
表わされる。
【数１】
【００１９】本実施例によると、１段目の適応ディジタ
ルフィルタ４３１〜４３６へ入力する手前に、遅延回路
４２１〜４２３を設け、これにより遅延時間を調節する
ことができる。従って直接波の到達に時間のかかるよう
な系の場合には、前記遅延回路にて時間遅れ分だけ遅延
させることにより、適応ディジタルフィルタの同一タッ
プ数で表現できるインパルス応答が長くできる。また、
２段目の適応ディジタルフィルタ４３４〜４３６へ入力
する手前に、遅延回路４２４〜４２６を設けたことによ
り、１段目の適応ディジタルフィルタと組み合わせて、
様々な音環境でのインパルス応答を表現できる。なお、
これら遅延回路の遅延時間の設定は、適応ディジタルフ
ィルタの遅延時間を種々変化させて適応ディジタルフィ
ルタの係数調整を行い、その誤差がより小さくなるよう
な遅延時間とすることにより行う。
ルフィルタ４３１〜４３６へ入力する手前に、遅延回路
４２１〜４２３を設け、これにより遅延時間を調節する
ことができる。従って直接波の到達に時間のかかるよう
な系の場合には、前記遅延回路にて時間遅れ分だけ遅延
させることにより、適応ディジタルフィルタの同一タッ
プ数で表現できるインパルス応答が長くできる。また、
２段目の適応ディジタルフィルタ４３４〜４３６へ入力
する手前に、遅延回路４２４〜４２６を設けたことによ
り、１段目の適応ディジタルフィルタと組み合わせて、
様々な音環境でのインパルス応答を表現できる。なお、
これら遅延回路の遅延時間の設定は、適応ディジタルフ
ィルタの遅延時間を種々変化させて適応ディジタルフィ
ルタの係数調整を行い、その誤差がより小さくなるよう
な遅延時間とすることにより行う。
【００２０】図３はそのインパルス応答の１例を示して
いる。この図に示すような直接波と反射波の時間間隔が
長いような系では、遅延回路により、直接波は１段目適
応フィルタ、反射波は２段目適応フィルタと分けてイン
パルス応答を求めることができ、タップ数を有効に使え
る。また、１段目、２段目の適応ディジタルフィルタの
タップ数を、それぞれ可変設定可能な構造としておけ
ば、遅延回路と組合せて適応ディジタルフィルタの調整
がより正確且つ容易に行えるようになる。さらに、適応
フィルタは２段だけではなく、状況、周辺環境に応じ
て、３段、４段と並列に複数の適応フィルタを設けるこ
とができる。また、差信号ｅ１及びｅ２をＤ／Ａ変換し
て音声として出力することで、適応ディジタルフィルタ
の集束状況を知ることができる。
いる。この図に示すような直接波と反射波の時間間隔が
長いような系では、遅延回路により、直接波は１段目適
応フィルタ、反射波は２段目適応フィルタと分けてイン
パルス応答を求めることができ、タップ数を有効に使え
る。また、１段目、２段目の適応ディジタルフィルタの
タップ数を、それぞれ可変設定可能な構造としておけ
ば、遅延回路と組合せて適応ディジタルフィルタの調整
がより正確且つ容易に行えるようになる。さらに、適応
フィルタは２段だけではなく、状況、周辺環境に応じ
て、３段、４段と並列に複数の適応フィルタを設けるこ
とができる。また、差信号ｅ１及びｅ２をＤ／Ａ変換し
て音声として出力することで、適応ディジタルフィルタ
の集束状況を知ることができる。
【００２１】本実施例は、騒音を発生する一般機器、例
えば、掃除機、空調機、車両等にも適用できる。さら
に、騒音もしくは振動の情報を得るセンサ、及びスピー
カ等はその個数を増やすことで測定点を増やすことがで
きる。また、センサとしては、マイクロフォン、振動セ
ンサだけでなく、時間的に変動する時系列信号を検出で
きるセンサ、例えば、アンテナ、光センサ等であっても
本実施例を適用できることはいうまでもない。
えば、掃除機、空調機、車両等にも適用できる。さら
に、騒音もしくは振動の情報を得るセンサ、及びスピー
カ等はその個数を増やすことで測定点を増やすことがで
きる。また、センサとしては、マイクロフォン、振動セ
ンサだけでなく、時間的に変動する時系列信号を検出で
きるセンサ、例えば、アンテナ、光センサ等であっても
本実施例を適用できることはいうまでもない。
【００２２】図４は、本発明になる能動消音装置の一実
施例を示すブロック図である。本装置は、機器６からの
騒音の情報を得るマイクロフォン等の基準センサ７と、
騒音を低減したい目標位置に前記騒音から伝播する音波
と逆位相同振幅の音波（付加音波）を送るための付加音
源であるスピーカ８と、前記目標位置に設置されて機器
６からの音波及びスピーカ８からの付加音波を検出する
消音エラーセンサとしてのマイクロフォン１０と、スピ
ーカ８からの付加音波の発生を制御するコントローラ９
から成っている。
施例を示すブロック図である。本装置は、機器６からの
騒音の情報を得るマイクロフォン等の基準センサ７と、
騒音を低減したい目標位置に前記騒音から伝播する音波
と逆位相同振幅の音波（付加音波）を送るための付加音
源であるスピーカ８と、前記目標位置に設置されて機器
６からの音波及びスピーカ８からの付加音波を検出する
消音エラーセンサとしてのマイクロフォン１０と、スピ
ーカ８からの付加音波の発生を制御するコントローラ９
から成っている。
【００２３】コントローラ９には、ここの制御に必要
な、スピーカ８への入力信号とマイクロフォン１０の出
力信号との間の推定伝達関数９３が用いられている。図
５は、この推定伝達関数９３の推定を行うための測定シ
ステムを示すブロック図で、本発明の適応信号処理方法
を用いたシステムとなっている。そして図５において、
スピーカ８とマイクロフォン１０は、図４の実際の減音
装置で用いるときと同じ環境下の同じ位置に配置されて
いるものとする。図５に於て、雑音信号発生源５０１か
らの基準信号ｓ２１はＡ／Ｄ変換器５０２でディジタル
形態の基準信号ｓ２２とされ、該基準信号ｓ２２はＤ／
Ａ変換器５０３とパワーアンプ５０４を介してスピーカ
８に与えられ、音波が放射される。また、前記基準信号
ｓ２２は、１段目の適応ディジタルフィルタ（ＡＤＦ）
５０５に入力され、フィルタ係数との畳み込み演算が行
われて信号ｓ２３が生成される。スピーカ８からの放射
音波はマイクロフォン１０で検出され、Ａ／Ｄ変換器５
０９によりディジタル化され、ディジタル信号ｓ２５と
される。この信号ｓ２５と上記信号ｓ２３との差信号ｅ
２１が減算回路５１０で求められ、この差信号ｅ２１が
最小となるように、適応アルゴリズム５０７は適応ディ
ジタルフィルタ５０５の係数を更新していく。このと
き、必要に応じて遅延回路５１１により基準信号ｓ２２
を遅延する。
な、スピーカ８への入力信号とマイクロフォン１０の出
力信号との間の推定伝達関数９３が用いられている。図
５は、この推定伝達関数９３の推定を行うための測定シ
ステムを示すブロック図で、本発明の適応信号処理方法
を用いたシステムとなっている。そして図５において、
スピーカ８とマイクロフォン１０は、図４の実際の減音
装置で用いるときと同じ環境下の同じ位置に配置されて
いるものとする。図５に於て、雑音信号発生源５０１か
らの基準信号ｓ２１はＡ／Ｄ変換器５０２でディジタル
形態の基準信号ｓ２２とされ、該基準信号ｓ２２はＤ／
Ａ変換器５０３とパワーアンプ５０４を介してスピーカ
８に与えられ、音波が放射される。また、前記基準信号
ｓ２２は、１段目の適応ディジタルフィルタ（ＡＤＦ）
５０５に入力され、フィルタ係数との畳み込み演算が行
われて信号ｓ２３が生成される。スピーカ８からの放射
音波はマイクロフォン１０で検出され、Ａ／Ｄ変換器５
０９によりディジタル化され、ディジタル信号ｓ２５と
される。この信号ｓ２５と上記信号ｓ２３との差信号ｅ
２１が減算回路５１０で求められ、この差信号ｅ２１が
最小となるように、適応アルゴリズム５０７は適応ディ
ジタルフィルタ５０５の係数を更新していく。このと
き、必要に応じて遅延回路５１１により基準信号ｓ２２
を遅延する。
【００２４】さらに、前記基準信号ｓ２２は遅延回路５
１２により必要に応じて一定時間遅延され、２段目の適
応ディジタルフィルタ（ＡＤＦ）５０６に入力され、フ
ィルタ係数との畳み込み演算が行われて信号ｓ２４が生
成される。この出力信号ｓ２４と、１段目の適応ディジ
タルフィルタで処理しきれなかった前記差信号ｅ２１と
の差ｅ２２が減算回路５１３で求められ、この差信号ｅ
２２が最小となるように、適応アルゴリズム５０８は適
応ディジタルフィルタ５０６の係数を更新していく。以
上の動作をくり返して各適応ディジタルフィルタの係数
を収束させ、その収束した値から推定伝達関数９３が求
められる。
１２により必要に応じて一定時間遅延され、２段目の適
応ディジタルフィルタ（ＡＤＦ）５０６に入力され、フ
ィルタ係数との畳み込み演算が行われて信号ｓ２４が生
成される。この出力信号ｓ２４と、１段目の適応ディジ
タルフィルタで処理しきれなかった前記差信号ｅ２１と
の差ｅ２２が減算回路５１３で求められ、この差信号ｅ
２２が最小となるように、適応アルゴリズム５０８は適
応ディジタルフィルタ５０６の係数を更新していく。以
上の動作をくり返して各適応ディジタルフィルタの係数
を収束させ、その収束した値から推定伝達関数９３が求
められる。
【００２５】上記のようにして求められた推定伝達関数
９３を用いた図４の能動消音装置の動作を次ぎに説明す
る。機器１の騒音情報として前記基準センサ７から出力
された検出信号ｓ１１は、前記Ａ／Ｄ変換器９２でディ
ジタル化されて基準信号ｓ１２とされ、さらに適応ディ
ジタルフィルタ９０の畳み込み演算が行われて信号ｓ１
６に変換される。この信号ｓ１６は、Ｄ／Ａ変換器９４
でアナログ化され、パワーアンプで増幅されてスピーカ
８から放射される。前記マイクロフォン１０はこの放射
音と機器１からの騒音の加わった音波を検出し、その出
力はＡ／Ｄ変換器９５でディジタル信号ｓ１４とされ
る。一方、前記基準信号ｓ１２と前記推定伝達関数９３
との畳み込み演算により仮想入力信号ｓ１３が生成さ
れ、この信号ｓ１３と前記信号ｓ１４を用いて、適応ア
ルゴリズム９１により信号ｓ１４の自乗値が最小となる
ように、適応ディジタルフィルタ９０の係数が更新され
る。
９３を用いた図４の能動消音装置の動作を次ぎに説明す
る。機器１の騒音情報として前記基準センサ７から出力
された検出信号ｓ１１は、前記Ａ／Ｄ変換器９２でディ
ジタル化されて基準信号ｓ１２とされ、さらに適応ディ
ジタルフィルタ９０の畳み込み演算が行われて信号ｓ１
６に変換される。この信号ｓ１６は、Ｄ／Ａ変換器９４
でアナログ化され、パワーアンプで増幅されてスピーカ
８から放射される。前記マイクロフォン１０はこの放射
音と機器１からの騒音の加わった音波を検出し、その出
力はＡ／Ｄ変換器９５でディジタル信号ｓ１４とされ
る。一方、前記基準信号ｓ１２と前記推定伝達関数９３
との畳み込み演算により仮想入力信号ｓ１３が生成さ
れ、この信号ｓ１３と前記信号ｓ１４を用いて、適応ア
ルゴリズム９１により信号ｓ１４の自乗値が最小となる
ように、適応ディジタルフィルタ９０の係数が更新され
る。
【００２６】本実施例は、騒音を発生する一般機器、例
えば、空調機、車両等にも適用できる。また上記実施例
において、騒音もしくは振動の情報を得るセンサ、及び
スピーカ等は個数を増やすことができる。さらに、図５
の推定伝達関数を求める際の適応フィルタは２段だけで
はなく、状況、周辺環境に応じて、３段、４段と並列に
複数の適応フィルタを設けることができる。
えば、空調機、車両等にも適用できる。また上記実施例
において、騒音もしくは振動の情報を得るセンサ、及び
スピーカ等は個数を増やすことができる。さらに、図５
の推定伝達関数を求める際の適応フィルタは２段だけで
はなく、状況、周辺環境に応じて、３段、４段と並列に
複数の適応フィルタを設けることができる。
【００２７】図６は、本発明の適応信号処理方法を用い
た音環境シミュレータの実施例を示すブロック図で、騒
音源である機器１１の騒音もしくは振動の情報を得るマ
イクロフォン等のセンサ１２ａ、１２ｂ、１２ｃと、機
器１１より一定距離離れた評価点（人間の耳元付近）に
設置されたマイクロフォン１３と、センサ１２ａ、１２
ｂ、１２ｃで検出された信号がマイクロフォン１３に対
してどのくらい寄与しているかを解析するための寄与率
解析装置１４と、その寄与率解析装置１４で得られた寄
与率等をセンサ１２ａ、１２ｂ、１２ｃで検出した信号
対応に分配し、評価点での騒音を合成するための音合成
手段１５と、合成した音を出力するスピーカ１６ａ、１
６ｂ、及び合成音を再生する部屋１７に設置されたダミ
ーヘッド１８から成っている。ここで音声合成部１５の
構成は図７に示されている。
た音環境シミュレータの実施例を示すブロック図で、騒
音源である機器１１の騒音もしくは振動の情報を得るマ
イクロフォン等のセンサ１２ａ、１２ｂ、１２ｃと、機
器１１より一定距離離れた評価点（人間の耳元付近）に
設置されたマイクロフォン１３と、センサ１２ａ、１２
ｂ、１２ｃで検出された信号がマイクロフォン１３に対
してどのくらい寄与しているかを解析するための寄与率
解析装置１４と、その寄与率解析装置１４で得られた寄
与率等をセンサ１２ａ、１２ｂ、１２ｃで検出した信号
対応に分配し、評価点での騒音を合成するための音合成
手段１５と、合成した音を出力するスピーカ１６ａ、１
６ｂ、及び合成音を再生する部屋１７に設置されたダミ
ーヘッド１８から成っている。ここで音声合成部１５の
構成は図７に示されている。
【００２８】この構成において、寄与率解析装置１４
は、図１の実施例で説明したものであって、センサ１２
ａ〜１２ｃが検出した騒音が、マイクロフォン１３で検
出した評価点の騒音にどれくらい寄与しているかの寄与
率を算出して出力する。この各寄与率は、音声合成手段
１５の寄与率分配回路１５１により各センサ１２ａ〜１
２ｃ対応の寄与率に分配され、ディジタルフィルタ１５
２ａ〜１５２ｃへそれぞれ入力される。また、前記寄与
率解析手段１４の収束時の各適応ディジタルフィルタ係
数がディジタルフィルタ１５２ａ〜１５２ｃにその係数
としてコピーされる。各ディジタルフィルタ１５２ａ〜
１５２ｃに於ては、そのフィルタの伝達関数と入力寄与
率との畳み込み演算が行われ、ディジタル信号ｓａ、ｓ
ｂ、ｓｃが出力される。これらの信号は加算回路１５３
で加算され、合成信号ｓをとして空間補正回路１５４へ
入力される。
は、図１の実施例で説明したものであって、センサ１２
ａ〜１２ｃが検出した騒音が、マイクロフォン１３で検
出した評価点の騒音にどれくらい寄与しているかの寄与
率を算出して出力する。この各寄与率は、音声合成手段
１５の寄与率分配回路１５１により各センサ１２ａ〜１
２ｃ対応の寄与率に分配され、ディジタルフィルタ１５
２ａ〜１５２ｃへそれぞれ入力される。また、前記寄与
率解析手段１４の収束時の各適応ディジタルフィルタ係
数がディジタルフィルタ１５２ａ〜１５２ｃにその係数
としてコピーされる。各ディジタルフィルタ１５２ａ〜
１５２ｃに於ては、そのフィルタの伝達関数と入力寄与
率との畳み込み演算が行われ、ディジタル信号ｓａ、ｓ
ｂ、ｓｃが出力される。これらの信号は加算回路１５３
で加算され、合成信号ｓをとして空間補正回路１５４へ
入力される。
【００２９】記憶媒体１５５には、様々な伝達特性（残
響特性等）を持つ空間での伝達特性の測定値がデータベ
ース化して格納されており、空間補正回路１５４では、
記憶媒体１５５から呼び出した伝達特性と上記合成信号
ｓとの畳み込み演算を行う。さらに、合成音を再生する
部屋１７内でのスピーカ１６ａ、１６ｂから評価する人
間もしくはダミーヘッド１８までの逆伝達特性を表わす
逆フィルタ１５６と上記空間補正回路１５４出力との畳
み込み演算を行い、演算結果をＤ／Ａ変換器でアナログ
信号に変換し、スピーカ１６ａ、１６ｂから出力する。
このようにして様々な特性を有する空間での製品音が再
現できる。
響特性等）を持つ空間での伝達特性の測定値がデータベ
ース化して格納されており、空間補正回路１５４では、
記憶媒体１５５から呼び出した伝達特性と上記合成信号
ｓとの畳み込み演算を行う。さらに、合成音を再生する
部屋１７内でのスピーカ１６ａ、１６ｂから評価する人
間もしくはダミーヘッド１８までの逆伝達特性を表わす
逆フィルタ１５６と上記空間補正回路１５４出力との畳
み込み演算を行い、演算結果をＤ／Ａ変換器でアナログ
信号に変換し、スピーカ１６ａ、１６ｂから出力する。
このようにして様々な特性を有する空間での製品音が再
現できる。
【００３０】本実施例は騒音を発生する一般機器、例え
ば、掃除機、空調機、車両等にも適用出来る。さらに、
騒音もしくは振動の情報を得るセンサ、および、スピー
カ等は個数を増やすことが出来る。さらに、寄与率解析
装置１４に於る適応フィルタは２段だけではなく、状
況、周辺環境に応じて、３段、４段と並列に複数の適応
フィルタを設けるようにしてもよい。
ば、掃除機、空調機、車両等にも適用出来る。さらに、
騒音もしくは振動の情報を得るセンサ、および、スピー
カ等は個数を増やすことが出来る。さらに、寄与率解析
装置１４に於る適応フィルタは２段だけではなく、状
況、周辺環境に応じて、３段、４段と並列に複数の適応
フィルタを設けるようにしてもよい。
【００３１】
【発明の効果】本発明によれば、適応信号処理におい
て、適応フィルタを２段設け、さらに、各適応フィルタ
と遅延時間の組み合わせによって、高残響場等様々な環
境下での音響系のインパルス応答の算出や音場の伝達関
数の推定を高精度に行え、また同様な環境下での消音装
置や音環境シミュレータの実現が可能になる効果があ
る。
て、適応フィルタを２段設け、さらに、各適応フィルタ
と遅延時間の組み合わせによって、高残響場等様々な環
境下での音響系のインパルス応答の算出や音場の伝達関
数の推定を高精度に行え、また同様な環境下での消音装
置や音環境シミュレータの実現が可能になる効果があ
る。
【図１】寄与率解析装置の一実施例を示すブロック図で
ある。
ある。
【図２】寄与率解析システムの全体構成を示すブロック
図である。
図である。
【図３】インパルスレスポンスの一例を示す図である。
【図４】能動消音装置の全体構成を示すブロック図であ
る。
る。
【図５】図４に於る推定伝達関数を求めるシステムの一
実施例を示す図である。
実施例を示す図である。
【図６】本発明の音環境シミュレータの一実施例を示す
全体構成図である。
全体構成図である。
【図７】図６の音声合成手段の詳細を示すブロック図で
ある。
ある。
１、６、１１ 騒音を発生する機器
２ａ、２ｂ、２ｃ、７ センサ
３、１０、１３ マイクロフォン
４、１４ 寄与率解析装置
４２１〜４２６ 遅延回路
４３１〜４３３ １段目の適応ディジタルフィルタ
４３４〜４３６ ２段目の適応ディジタルフィルタ
４５１、４５２ 加算回路
４６１、４６２ 減算回路
４８１ 残差判定回路
５０５、５０６ 適応ディジタルフィルタ
５１０、５１３ 減算回路
５１１、５１２ 遅延回路
８ 消音用スピーカ
９ コントローラ
９０ 適応ディジタルフィルタ
９３ 推定伝達関数
１６ａ、１６ｂ 音再現用スピーカ
１５ 音合成手段
１８ 人間もしくはダミーヘッド
フロントページの続き
(72)発明者 大木 克夫
茨城県土浦市神立町502番地 株式会社
日立製作所機械研究所内
(72)発明者 太田 義注
茨城県土浦市神立町502番地 株式会社
日立製作所機械研究所内
(58)調査した分野(Int.Cl.7，ＤＢ名)
G01H 17/00
G01H 3/00
G10K 11/178
Claims (7)
- 【請求項１】 第１の１個又は複数個の時系列信号を１
個又は複数個の１段目適応フィルタの各々に入力し、該
適応フィルタの出力信号の和と前記第１の時系列信号と
は別の第２の時系列信号との差を表す第１の残差信号が
最小になるように、前記１段目適応フィルタの係数を更
新するとともに、 前記第１の時系列信号を１個又は複数個の２段目適応フ
ィルタの各々に入力し、該適応フィルタの出力信号の和
と前記第１の残差信号との差を表す第２の残差信号が最
小になるように前記２段目適応フィルタを更新すること
を特徴とする適応信号処理方法。 - 【請求項２】 前記第１の時系列信号を入力とする３段
目あるいは必要に応じてそれ以降の適応フィルタを設
け、該各適応フィルタの係数をその出力の和と１段前の
適応フィルタにより生成された残差信号との差が最小に
なるように更新することを特徴とする請求項１記載の適
応信号処理方法。 - 【請求項３】 前記１段目及び２段目の適応フィルタに
入力される前記第１の時系列信号を、その遅延時間が可
変な遅延手段により遅延できるようにしたことを特徴と
する請求項１記載の適応信号処理方法。 - 【請求項４】 前記２段目適応フィルタ側の遅延時間
を、前記１段目適応フィルタ側の遅延時間に連動して設
定するようにしたことを特徴とする請求項３記載の適応
信号処理方法。 - 【請求項５】 前記１段目適応フィルタ側の遅延時間
を、当該遅延時間を変化させたときに得られる前記第１
の残差信号が最小となる値に設定し、 さらに前記２段目適応フィルタ側の遅延時間を、当該遅
延時間を変化させたときに得られる前記第２の残差信号
が最小となる値に設定することを特徴とする請求項３記
載の適応信号処理方法。 - 【請求項６】 前記１段目及び２段目適応フィルタの各
々を、そのタップ数が可変な構造としたことを特徴する
請求項１記載の適応信号処理方法。 - 【請求項７】 前記第１の残差信号に対するしきい値を
設け、前記残差信号が前記しきい値より大きいときだけ
前記２段目の適応フィルタに於る係数の更新処理を行う
ようにしたことを特徴とする請求項１記載の適応信号処
理方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP16025995A JP3393738B2 (ja) | 1995-06-27 | 1995-06-27 | 適応信号処理方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP16025995A JP3393738B2 (ja) | 1995-06-27 | 1995-06-27 | 適応信号処理方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH0915035A JPH0915035A (ja) | 1997-01-17 |
| JP3393738B2 true JP3393738B2 (ja) | 2003-04-07 |
Family
ID=15711146
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP16025995A Expired - Fee Related JP3393738B2 (ja) | 1995-06-27 | 1995-06-27 | 適応信号処理方法 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3393738B2 (ja) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP5253268B2 (ja) * | 2009-03-30 | 2013-07-31 | 中部電力株式会社 | 音源・振動源探査システム |
| AT512755B1 (de) * | 2013-06-18 | 2015-03-15 | Avl List Gmbh | Verfahren zur Erstellung eines dreidimensionalen, akustischen Netzes |
-
1995
- 1995-06-27 JP JP16025995A patent/JP3393738B2/ja not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| JPH0915035A (ja) | 1997-01-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP5654513B2 (ja) | 音識別方法および装置 | |
| US4653102A (en) | Directional microphone system | |
| US8355510B2 (en) | Reduced latency low frequency equalization system | |
| JP5229053B2 (ja) | 信号処理装置、および信号処理方法、並びにプログラム | |
| US5774562A (en) | Method and apparatus for dereverberation | |
| KR102191736B1 (ko) | 인공신경망을 이용한 음성향상방법 및 장치 | |
| JP2006517072A (ja) | マルチチャネル信号を用いて再生部を制御する方法および装置 | |
| CN105103219A (zh) | 降低噪音的方法 | |
| JP3283423B2 (ja) | マイクロホン装置 | |
| JP4130779B2 (ja) | 音場制御システム及び音場制御方法 | |
| JP3393738B2 (ja) | 適応信号処理方法 | |
| JPH0511771A (ja) | 騒音制御装置 | |
| US20080285768A1 (en) | Method and System for Modifying and Audio Signal, and Filter System for Modifying an Electrical Signal | |
| EP3787311A1 (en) | Sound image reproduction device, sound image reproduction method and sound image reproduction program | |
| US9445195B2 (en) | Directivity control method and device | |
| JP3657770B2 (ja) | 適応等化システムにおける遅延時間測定方式 | |
| JP3380571B2 (ja) | 騒音キャンセル装置 | |
| JP2714098B2 (ja) | 音響周波数特性の補正方法 | |
| JP2000099039A (ja) | 拡声音の明瞭度改善方法及び装置 | |
| JP4522509B2 (ja) | オーディオ装置 | |
| JP3556427B2 (ja) | オーディオ装置の制御帯域決定方法 | |
| JP4176938B2 (ja) | 遅延時間設定方式 | |
| JPH11167383A (ja) | 適応等化システムの接続確認方式 | |
| JPH11262081A (ja) | 遅延時間設定方式 | |
| JP3264141B2 (ja) | 音環境シミュレータ |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| LAPS | Cancellation because of no payment of annual fees |