JP7444750B2 - 騒音制御システムおよび騒音制御方法 - Google Patents
騒音制御システムおよび騒音制御方法 Download PDFInfo
- Publication number
- JP7444750B2 JP7444750B2 JP2020160133A JP2020160133A JP7444750B2 JP 7444750 B2 JP7444750 B2 JP 7444750B2 JP 2020160133 A JP2020160133 A JP 2020160133A JP 2020160133 A JP2020160133 A JP 2020160133A JP 7444750 B2 JP7444750 B2 JP 7444750B2
- Authority
- JP
- Japan
- Prior art keywords
- noise
- control
- room
- sound
- error
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims description 42
- 238000012545 processing Methods 0.000 claims description 52
- 238000012546 transfer Methods 0.000 claims description 50
- 230000003044 adaptive effect Effects 0.000 claims description 46
- 238000002360 preparation method Methods 0.000 claims description 14
- 239000000284 extract Substances 0.000 claims description 6
- 238000000605 extraction Methods 0.000 claims description 5
- 230000007274 generation of a signal involved in cell-cell signaling Effects 0.000 claims description 3
- 238000002474 experimental method Methods 0.000 description 29
- 238000005259 measurement Methods 0.000 description 22
- 238000010586 diagram Methods 0.000 description 20
- 230000008569 process Effects 0.000 description 16
- 230000000694 effects Effects 0.000 description 6
- 230000003321 amplification Effects 0.000 description 5
- 238000003199 nucleic acid amplification method Methods 0.000 description 5
- 238000013461 design Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000001902 propagating effect Effects 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 238000009434 installation Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000001131 transforming effect Effects 0.000 description 2
- 238000004378 air conditioning Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000009530 blood pressure measurement Methods 0.000 description 1
- 238000009413 insulation Methods 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
Images
Landscapes
- Soundproofing, Sound Blocking, And Sound Damping (AREA)
Description
本発明は、騒音制御システムおよび騒音制御方法に関する。
閉空間内の共鳴時の音響エネルギーを最小にするアクティブ制御手法が研究されている（例えば、非特許文献１参照）。非特許文献１に記載される方法は、剛壁で囲まれた１室内のモード周波数（共振周波数）に着目している。具体的には、モード周波数が放射されている室内において、制御音源と制御点を適切に配置し、室内の全音響エネルギーが最小となるような制御を行うことで、音圧レベルが全体的に低下して騒音を打ち消している。なお、非特許文献１では、既知の純音を騒音源としている。
非特許文献１におけるシステムの構成図を図１７に示す。図１７に示す従来の騒音制御システム９０１は、主に、参照マイクロホン９０２と、誤差マイクロホン９０３と、制御スピーカ９０４と、誤差経路特性モデル９１２と、適応フィルタ９１３と、ＬＭＳ（Least Mean Square）処理部９１４とを備える。シグナルジェネレータ９０９ｂとスピーカ９０９ａとの組合せを騒音源としておりスピーカ９０９ａから放射される騒音は純音である。
騒音制御システム９０１では、Filtered-x LMSアルゴリズムを含むアクティブノイズコントロール技術を用いている。騒音源近傍において騒音を参照マイクロホン９０２で収音し、参照信号は適応フィルタ９１３に入力される。この適応フィルタ９１３は、例えばFIR適応型デジタルフィルタであり、Filtered-x LMSアルゴリズムに基づいてフィルタ係数が更新される。このフィルタ係数を用いて畳み込み和演算のようなフィルタ処理を行うことにより、フィルタ出力信号が計算される。この出力信号は、「－１」を乗算されて制御スピーカ９０４に供給される。制御スピーカ９０４から放射する音波は、騒音源から誤差マイクロホン９０３に伝搬している騒音に干渉し、騒音を打ち消す。
誤差マイクロホン９０３は、室の角に配置されている。騒音と制御スピーカ９０４からの放射音とが干渉した後の音、すなわち騒音と制御スピーカの放射音との誤差をこの誤差マイクロホン９０３で検出している。誤差マイクロホン９０３の出力は、誤差信号としてＬＭＳ処理部９１４に供給される。ＬＭＳ処理部９１４は、Filtered-x LMSアルゴリズムに基づく処理を行い、参照信号と誤差信号に応じて、適応フィルタ９１３と誤差経路特性モデル９１２とを合成した伝達関数と、参照マイクロホン９０２から誤差マイクロホン９０３までの伝達関数とが等しくなるように、適応フィルタ９１３のフィルタ係数を更新する。適応フィルタ９１３と誤差経路特性モデル９１２とを合成した伝達関数と、参照マイクロホン９０２から誤差マイクロホン９０３までの伝達関数が等しくなることで、誤差マイクロホン９０３において制御スピーカ９０４からの放射音を騒音に干渉させることにより打ち消すことができる。
伊勢 史郎、他２名、「境界要素法による室内のアクティブモード制御の解析」、日本音響学会誌51巻1号、一般社団法人日本音響学会、1995年、pp.25-33
非特許文献１に記載される技術では、１つの室内について議論されており、音源室に隣接する室については考えられていないので、隣室が存在している場合に制御できるのは音源室側のモード周波数のみであって隣室側のモード周波数を制御できない。したがって、騒音源から隣室に伝搬する周波数（特に低周波）が隣室のモード周波数（共振周波数）と一致する場合に音圧が増幅し、隣室に居る人間が不快に感じるという課題があった。
例えば、空調等の設備機械室の周囲に静粛性が求められる会議室が配置されることがある。設備機器から発生する低周波音は壁面を透過して隣室まで伝搬しやすく、室形状により決定するモード周波数（共振周波数）に合致すると増幅する特徴を持つ。その為、会議室に居る人は、設備機器から発生する低周波音を不快に感じる。
なお、重い材料を用いて低周波を遮音することが一般的に行われているが、設計の条件などにより遮音するのに必要な重量の部材を設置することが難しい場合があり、その場合には音圧低減効果を十分に得ることができない。
また、騒音源室の隣室に騒音を低減させるための制御装置を設置することも考えられるが、隣室に制御装置を設置することで利用可能なスペースが制限されてしまう。その為、上記のように騒音源室の隣室が会議室として運用されている場合などでは、制御装置を隣室に設置することは望ましくない。
また、騒音源室の隣室に騒音を低減させるための制御装置を設置することも考えられるが、隣室に制御装置を設置することで利用可能なスペースが制限されてしまう。その為、上記のように騒音源室の隣室が会議室として運用されている場合などでは、制御装置を隣室に設置することは望ましくない。
このような観点から、本発明は、騒音源室内での制御によって隣接する室内の騒音を低減することができる騒音制御システムおよび騒音制御方法を提供する。
本発明に係る騒音制御システムは、騒音源を備える騒音源室に隣接する受音室での騒音を低減させる騒音制御システムである。この騒音制御システムは、前記騒音源室内に設置される参照マイクロホン、誤差マイクロホンおよび制御スピーカと、信号処理を行うことで前記騒音源室内の騒音を打ち消す制御音を発生させる制御信号を前記制御スピーカに出力する信号処理装置と、を備える。
前記参照マイクロホンは、前記制御音によって打ち消される前の騒音を検出可能な場所に設置され、当該騒音を参照信号として前記信号処理装置に出力する。
前記制御スピーカは、前記騒音源室内の室隅に設置され、前記制御音を前記騒音源室内に放出する。
前記誤差マイクロホンは、前記騒音源室内の室隅に設置され、前記制御音によって打ち消された後の騒音を誤差信号として前記信号処理装置に出力する。
前記信号処理装置は、前記参照信号および前記誤差信号から特定の周波数の成分を抽出するバンドパスフィルタと、前記バンドパスフィルタによって抽出された前記参照信号および前記誤差信号に基づき、ＬＭＳアルゴリズムによってフィルタ係数が更新される適応フィルタと、を備える。
前記適応フィルタは、前記バンドパスフィルタによって抽出された前記参照信号を用いて前記制御信号を生成する。前記バンドパスフィルタは、前記騒音源室で音が増幅しやすい第一モード周波数の成分および前記受音室で音が増幅しやすい第二モード周波数の成分を抽出するものである。
前記参照マイクロホンは、前記制御音によって打ち消される前の騒音を検出可能な場所に設置され、当該騒音を参照信号として前記信号処理装置に出力する。
前記制御スピーカは、前記騒音源室内の室隅に設置され、前記制御音を前記騒音源室内に放出する。
前記誤差マイクロホンは、前記騒音源室内の室隅に設置され、前記制御音によって打ち消された後の騒音を誤差信号として前記信号処理装置に出力する。
前記信号処理装置は、前記参照信号および前記誤差信号から特定の周波数の成分を抽出するバンドパスフィルタと、前記バンドパスフィルタによって抽出された前記参照信号および前記誤差信号に基づき、ＬＭＳアルゴリズムによってフィルタ係数が更新される適応フィルタと、を備える。
前記適応フィルタは、前記バンドパスフィルタによって抽出された前記参照信号を用いて前記制御信号を生成する。前記バンドパスフィルタは、前記騒音源室で音が増幅しやすい第一モード周波数の成分および前記受音室で音が増幅しやすい第二モード周波数の成分を抽出するものである。
前記信号処理装置は、前記制御スピーカと前記誤差マイクロホンとの間の伝達関数が組み込まれた誤差経路特性モデルをさらに有するのがよい。
前記適応フィルタのフィルタ係数は、前記適応フィルタと前記誤差経路特性モデルとを合成した伝達関数と、前記参照マイクロホンと前記誤差マイクロホンとの間の伝達関数とが等しくなるように更新される。
前記適応フィルタのフィルタ係数は、前記適応フィルタと前記誤差経路特性モデルとを合成した伝達関数と、前記参照マイクロホンと前記誤差マイクロホンとの間の伝達関数とが等しくなるように更新される。
本発明に係る騒音制御システムにおいては、騒音源室および受音室のモード周波数だけ通過するバンドパスフィルタを用いることで参照信号の絞り込みを行う。これにより、騒音源が放射する周波数と騒音源室および受音室のいずれかのモード周波数とが一致した際に生じる受音室での音圧の増幅を騒音源室内で抑えることができる。その為、受音室全体に渡って音圧低減が可能である。また、当該システムを構成する機器は、騒音源室に設置することができ、受音室にこれらの機器を設置しなくてよい。その為、受音室ではスペースの利用に制限を受けないので、スペースを有効に利用することができる。
また、本発明に係る騒音制御方法は、騒音源を備える騒音源室に隣接する受音室での騒音を低減させる騒音制御方法である。
この騒音制御方法は、参照マイクロホン、誤差マイクロホンおよび制御スピーカを前記騒音源室内に設置する準備工程と、信号処理を行うことで前記騒音源室内の騒音を打ち消す制御音を発生させる制御信号を前記制御スピーカに出力する制御工程と、を有する。
前記準備工程では、前記参照マイクロホンを前記制御音によって打ち消される前の騒音を検出可能な場所に設置し、前記制御スピーカおよび前記誤差マイクロホンを前記騒音源室内の室隅に設置する。
前記制御工程には、抽出工程と制御信号生成工程とが含まれる。抽出工程では、前記参照マイクロホンから取得した参照信号および前記誤差マイクロホンから取得した誤差信号から、バンドパスフィルタを用いて前記騒音源室で音が増幅しやすい第一モード周波数の成分および前記受音室で音が増幅しやすい第二モード周波数の成分を抽出する。
制御信号生成工程では、前記バンドパスフィルタによって抽出された前記参照信号および前記誤差信号に基づき、ＬＭＳアルゴリズムによって適応フィルタのフィルタ係数を更新し、また、前記適応フィルタが前記バンドパスフィルタによって抽出された前記参照信号を用いて前記制御信号を生成する。
この騒音制御方法は、参照マイクロホン、誤差マイクロホンおよび制御スピーカを前記騒音源室内に設置する準備工程と、信号処理を行うことで前記騒音源室内の騒音を打ち消す制御音を発生させる制御信号を前記制御スピーカに出力する制御工程と、を有する。
前記準備工程では、前記参照マイクロホンを前記制御音によって打ち消される前の騒音を検出可能な場所に設置し、前記制御スピーカおよび前記誤差マイクロホンを前記騒音源室内の室隅に設置する。
前記制御工程には、抽出工程と制御信号生成工程とが含まれる。抽出工程では、前記参照マイクロホンから取得した参照信号および前記誤差マイクロホンから取得した誤差信号から、バンドパスフィルタを用いて前記騒音源室で音が増幅しやすい第一モード周波数の成分および前記受音室で音が増幅しやすい第二モード周波数の成分を抽出する。
制御信号生成工程では、前記バンドパスフィルタによって抽出された前記参照信号および前記誤差信号に基づき、ＬＭＳアルゴリズムによって適応フィルタのフィルタ係数を更新し、また、前記適応フィルタが前記バンドパスフィルタによって抽出された前記参照信号を用いて前記制御信号を生成する。
前記準備工程には、モード周波数計測工程と、周波数特性計測工程と、バンドパスフィルタ作成工程とが含まれるのがよい。
モード周波数計測工程では、前記第一モード周波数および前記第二モード周波数を計測する。周波数特性計測工程では、前記騒音の中で音が大きい卓越周波数を前記騒音源室内で計測する。バンドパスフィルタ作成工程では、前記第一モード周波数と前記卓越周波数とが一致する周波数を含む第一狭帯域、および前記第二モード周波数と前記卓越周波数とが一致する周波数を含む第二狭帯域を通過させるバンドパスフィルタを作成する。
モード周波数計測工程では、前記第一モード周波数および前記第二モード周波数を計測する。周波数特性計測工程では、前記騒音の中で音が大きい卓越周波数を前記騒音源室内で計測する。バンドパスフィルタ作成工程では、前記第一モード周波数と前記卓越周波数とが一致する周波数を含む第一狭帯域、および前記第二モード周波数と前記卓越周波数とが一致する周波数を含む第二狭帯域を通過させるバンドパスフィルタを作成する。
本発明に係る騒音制御方法においては、騒音源室および受音室のモード周波数だけ通過するバンドパスフィルタを用いることで参照信号の絞り込みを行う。これにより、騒音源が放射する周波数と騒音源室および受音室のいずれかのモード周波数とが一致した際に生じる受音室での音圧の増幅を騒音源室内で抑えることができる。その為、受音室全体に渡って音圧低減が可能である。また、当該方法で使用する機器は、騒音源室に設置することができ、受音室にこれらの機器を設置しなくてよい。その為、受音室ではスペースの利用に制限を受けないので、スペースを有効に利用することができる。
本発明によれば、騒音源室内での制御によって隣接する室内の騒音を低減することができる。
以下、本発明の実施をするための形態を、適宜図面を参照しながら詳細に説明する。各図は、本発明を十分に理解できる程度に、概略的に示してあるに過ぎない。よって、本発明は、図示例のみに限定されるものではない。なお、各図において、共通する構成要素や同様な構成要素については、同一の符号を付し、それらの重複する説明を省略する。
＜実施形態に係る騒音制御システムの構成＞
図１を参照して、実施形態に係る騒音制御システム１について説明する。図１は、実施形態に係る騒音制御システムの構成図である。騒音制御システム１は、アクティブノイズコントロール技術によって騒音を低減するシステムである。騒音制御システム１を用いれば、騒音源を備える室（「騒音源室」と称する）に隣接する室での騒音を低減できる。騒音源に限定はなく、例えば設備機械などである。騒音源室に隣接する室は、騒音源から放射される騒音が伝搬される側（つまり、受け取る側）なので「受音室」と称する。
図１を参照して、実施形態に係る騒音制御システム１について説明する。図１は、実施形態に係る騒音制御システムの構成図である。騒音制御システム１は、アクティブノイズコントロール技術によって騒音を低減するシステムである。騒音制御システム１を用いれば、騒音源を備える室（「騒音源室」と称する）に隣接する室での騒音を低減できる。騒音源に限定はなく、例えば設備機械などである。騒音源室に隣接する室は、騒音源から放射される騒音が伝搬される側（つまり、受け取る側）なので「受音室」と称する。
図１に示すように、騒音制御システム１の制御手法は、Filtered-x LMS(Least Mean Square)アルゴリズムを含むアクティブノイズコントロール技術を用いたものである。本実施形態では、フィードフォワード制御を想定して説明するが、フィードバック制御であってもよい。騒音制御システム１をフィードバック制御により実現する場合、後記する参照マイクロホン２を備えない構成になる（つまり、参照信号を用いない制御となる）。
騒音制御システム１は、主に、参照マイクロホン２と、誤差マイクロホン３と、制御スピーカ４と、信号処理装置５とを備える。参照マイクロホン２、誤差マイクロホン３および制御スピーカ４は騒音源室に設置され、信号処理装置５の設置場所は特に限定されないが騒音源室であるのがよい。
騒音制御システム１は、主に、参照マイクロホン２と、誤差マイクロホン３と、制御スピーカ４と、信号処理装置５とを備える。参照マイクロホン２、誤差マイクロホン３および制御スピーカ４は騒音源室に設置され、信号処理装置５の設置場所は特に限定されないが騒音源室であるのがよい。
信号処理装置５は、例えばデジタルシグナルプロセッサ（ＤＳＰ：Digital Signal Processor）を備え、デジタルシグナルプロセッサによるプログラム実行処理によって各機能が実現される。信号処理装置５は、参照マイクロホン２および誤差マイクロホン３から騒音源室内の騒音を取得し、信号処理を行うことで騒音を打ち消す制御音を発生させる制御信号を制御スピーカ４に出力する。詳細な処理は後述する。
図２を参照して、騒音制御システム１を構成する機器の配置について説明する。図２は、騒音制御システム１を構成する機器の配置例を示した図である。
参照マイクロホン２は、制御音によって打ち消される前の騒音を検出可能な場所に設置され、例えば騒音源９の近傍に配置される。騒音源９は例えば設備機械であり、通常稼動時に発生する騒音の周波数特定はある程度決まっている。参照マイクロホン２は、騒音を参照信号として信号処理装置５に出力する。
参照マイクロホン２は、制御音によって打ち消される前の騒音を検出可能な場所に設置され、例えば騒音源９の近傍に配置される。騒音源９は例えば設備機械であり、通常稼動時に発生する騒音の周波数特定はある程度決まっている。参照マイクロホン２は、騒音を参照信号として信号処理装置５に出力する。
制御スピーカ４の設置場所は、騒音源室Ｒ１の室隅であればよい。図２では、制御スピーカ４は、受音室Ｒ２に接する壁Ｒａ側の室隅に設置されているが、受音室Ｒ２に接する壁Ｒａ側の室隅でなくてもよい。なお、制御スピーカ４の高さ方向の位置は、限定されないが、本実施形態の制御スピーカ４は床近くに設置される。制御スピーカ４には、騒音を打ち消す制御音を発生させる制御信号が入力され、制御音を騒音源室Ｒ１内に放出する。
誤差マイクロホン３の設置場所は、騒音源室Ｒ１の室隅であればよい。図２に示すように、誤差マイクロホン３は、受音室Ｒ２に接する壁Ｒａ側の室隅に設置されるのが望ましい。なお、誤差マイクロホン３の高さ方向の位置は、限定されないが、本実施形態の誤差マイクロホン３は例えば天井近くに設置される。誤差マイクロホン３は、制御音によって打ち消された後の騒音を誤差信号として信号処理装置５に出力する（すなわち、騒音と制御スピーカ４の放射音との誤差を検出している）。
誤差マイクロホン３の設置場所は、騒音源室Ｒ１の室隅であればよい。図２に示すように、誤差マイクロホン３は、受音室Ｒ２に接する壁Ｒａ側の室隅に設置されるのが望ましい。なお、誤差マイクロホン３の高さ方向の位置は、限定されないが、本実施形態の誤差マイクロホン３は例えば天井近くに設置される。誤差マイクロホン３は、制御音によって打ち消された後の騒音を誤差信号として信号処理装置５に出力する（すなわち、騒音と制御スピーカ４の放射音との誤差を検出している）。
なお、制御スピーカ４および誤差マイクロホン３は、受音室Ｒ２に接する壁Ｒａ側の両方の室隅に設置されてもよい（図２の仮想線参照）。壁Ｒａの両方の室隅に制御スピーカ４および誤差マイクロホン３を配置した場合の信号処理装置５の構成は、例えば図３に示すようになる。図３は、実施形態の変形例に係る騒音制御システム１０１の構成図である。図３に示す騒音制御システム１０１では、Filtered-x LMSアルゴリズムをマルチチャンネル化したMEFX-LMS(Multiple Error Filtered x-LMS)アルゴリズムとなる。MEFX-LMSアルゴルズムについては、例えば「特開平７－２０８８３号公報」に記載がある。
図１を参照して、信号処理装置５について説明する。信号処理装置５は、主に、モード周波数抽出用のバンドパスフィルタ１１ａ，１１ｂと、誤差経路特性モデル１２と、適応フィルタ１３と、ＬＭＳ処理部１４と、パワーアンプ１５とを備える。なお、バンドパスフィルタ１１ａ，１１ｂは同じ機能であるので、区別せずに説明する場合に「バンドパスフィルタ１１」と表記する場合がある。
バンドパスフィルタ１１は、騒音源室Ｒ１（図２参照）で音が増幅しやすい第一モード周波数の成分および受音室Ｒ２（図２参照）で音が増幅しやすい第二モード周波数の成分を抽出する（通過させる）ものである。バンドパスフィルタ１１が抽出する周波数は、例えば騒音源室Ｒ１および受音室Ｒ２での計測により決定される。詳細は後述する。なお、バンドパスフィルタ１１が抽出する第一モード周波数や第二モード周波数の数は、複数であってもよい。
バンドパスフィルタ１１ａには、参照マイクロホン２から参照信号が入力される（なお、参照信号はマイクアンプによって増幅される）。バンドパスフィルタ１１ａによって抽出された参照信号は、誤差経路特性モデル１２および適応フィルタ１３に入力される。
バンドパスフィルタ１１ｂには、誤差マイクロホン３から誤差信号が入力される（なお、誤差信号はマイクアンプによって増幅される）。バンドパスフィルタ１１ｂによって抽出された誤差信号は、ＬＭＳ処理部１４に入力されてFiltered-x LMSアルゴリズムに基づいた演算に使用され、適応フィルタ１３のフィルタ係数が更新される。
誤差経路特性モデル１２は、制御スピーカ４と誤差マイクロホン３との間の伝達関数が組み込まれている。誤差経路特性モデル１２から出力された信号は、ＬＭＳ処理部１４に入力される。
バンドパスフィルタ１１ｂには、誤差マイクロホン３から誤差信号が入力される（なお、誤差信号はマイクアンプによって増幅される）。バンドパスフィルタ１１ｂによって抽出された誤差信号は、ＬＭＳ処理部１４に入力されてFiltered-x LMSアルゴリズムに基づいた演算に使用され、適応フィルタ１３のフィルタ係数が更新される。
誤差経路特性モデル１２は、制御スピーカ４と誤差マイクロホン３との間の伝達関数が組み込まれている。誤差経路特性モデル１２から出力された信号は、ＬＭＳ処理部１４に入力される。
適応フィルタ１３は、バンドパスフィルタ１１ａによって抽出された参照信号から制御信号を生成する。適応フィルタ１３は、Filtered-x LMSアルゴリズムに基づいてフィルタ係数が更新される。適応フィルタ１３は、このフィルタ係数を用いて畳み込み和演算のようなフィルタ処理によってフィルタ出力信号を計算する。フィルタ出力信号に「－１」を乗算することで制御信号となり、制御信号は制御スピーカ４に出力される（なお、制御信号はパワーアンプ１５によって増幅される）。制御スピーカ４から放射する音波は、騒音源９から誤差マイクロホン３に伝搬している騒音に干渉し、騒音を打ち消す。
ＬＭＳ処理部１４は、参照信号と誤差信号に応じて、適応フィルタ１３と誤差経路特性モデル１２とを合成した伝達関数と、参照マイクロホン２から誤差マイクロホン３までの伝達関数とが等しくなるように適応フィルタ１３のフィルタ係数を更新する。適応フィルタ１３と誤差経路特性モデル１２とを合成した伝達関数と、参照マイクロホン２から誤差マイクロホン３までの伝達関数とが等しくなることで、誤差マイクロホン３において制御スピーカ４からの放射音を騒音に干渉させて騒音の音圧を低減させる。
なお、本実施形態では、騒音源９の近傍に参照マイクロホン２を設置することで、騒音源９と参照マイクロホン２間の影響を考慮しないことにする。同様に、制御スピーカ４と誤差マイクロホン３とを近傍に設置した場合、制御スピーカ４と誤差マイクロホン３間の影響を考慮しなくてよい場合もある。その場合、信号処理装置５は、誤差経路特性モデル１２を有しない構成にすることができる。
図３を参照して、複数（例示は二つ）の制御スピーカ４および誤差マイクロホン３を用いる場合の騒音制御システム１０１を説明する。騒音制御システム１０１の制御手法は、Filtered-x LMSアルゴリズムをマルチチャンネル化したMEFX-LMSアルゴリズムを含むアクティブノイズコントロール技術を用いたものである。
騒音制御システム１０１は、主に、参照マイクロホン２と、二つの誤差マイクロホン３1，３2と、二つの制御スピーカ４1，４2と、信号処理装置１０５とを備える。参照マイクロホン２、誤差マイクロホン３1，３2および制御スピーカ４1，４2は騒音源室に設置され、信号処理装置１０５の設置場所は特に限定されないが騒音源室であるのがよい。誤差マイクロホン３1，３2および制御スピーカ４1，４2は、例えば受音室Ｒ２に接する壁Ｒａ側の両方の室隅にそれぞれ設置される（図２参照）。
信号処理装置１０５は、例えばデジタルシグナルプロセッサ（ＤＳＰ：Digital Signal Processor）を備え、デジタルシグナルプロセッサによるプログラム実行処理によって各機能が実現される。
信号処理装置１０５は、主に、モード周波数抽出用のバンドパスフィルタ１１ａ，１１ｂと、四つの誤差経路特性モデル１２11，１２21，１２12，１２22と、二つの適応フィルタ１３1，１３2と、二つのＬＭＳ処理部１４1，１４2とを備える。
信号処理装置１０５は、主に、モード周波数抽出用のバンドパスフィルタ１１ａ，１１ｂと、四つの誤差経路特性モデル１２11，１２21，１２12，１２22と、二つの適応フィルタ１３1，１３2と、二つのＬＭＳ処理部１４1，１４2とを備える。
図３中の符号Ｃ11は、第一制御スピーカ４1－第一誤差マイクロホン３1間の伝達関数を表している。符号Ｃ21は、第一制御スピーカ４1－第二誤差マイクロホン３2間の伝達関数を表している。符号Ｃ12は、第二制御スピーカ４2－第一誤差マイクロホン３1間の伝達関数を表している。符号Ｃ22は、第二制御スピーカ４2－第二誤差マイクロホン３2間の伝達関数を表している。
誤差経路特性モデル１２11は、第一制御スピーカ４1－第一誤差マイクロホン３1間の伝達関数Ｃ11の推定モデルである。誤差経路特性モデル１２21は、第一制御スピーカ４1－第二誤差マイクロホン３2間の伝達関数Ｃ21の推定モデルである。誤差経路特性モデル１２12は、第二制御スピーカ４2－第一誤差マイクロホン３1間の伝達関数Ｃ12の推定モデルである。誤差経路特性モデル１２22は、第二制御スピーカ４2－第二誤差マイクロホン３2間の伝達関数Ｃ22の推定モデルである。誤差経路特性モデルは、伝達関数の記号にハット記号を付して表現される場合があり、図３では誤差経路特性モデル１２11，１２21，１２12，１２22のブロックの内部に当該記号を表記している。
バンドパスフィルタ１１ａには、参照マイクロホン２から参照信号が入力される（なお、参照信号はマイクアンプによって増幅される）。バンドパスフィルタ１１ａによって抽出された参照信号は、誤差経路特性モデル１２11，１２21，１２12，１２22および適応フィルタ１３1，１３2に入力される。
バンドパスフィルタ１１ｂには、誤差マイクロホン３1，３2から誤差信号が入力される（なお、誤差信号はマイクアンプによって増幅される）。バンドパスフィルタ１１ｂによって抽出された誤差信号は、ＬＭＳ処理部１４1，１４2に入力されてMEFX-LMSアルゴリズムに基づいた演算に使用され、適応フィルタ１３1，１３2のフィルタ係数が更新される。
バンドパスフィルタ１１ｂには、誤差マイクロホン３1，３2から誤差信号が入力される（なお、誤差信号はマイクアンプによって増幅される）。バンドパスフィルタ１１ｂによって抽出された誤差信号は、ＬＭＳ処理部１４1，１４2に入力されてMEFX-LMSアルゴリズムに基づいた演算に使用され、適応フィルタ１３1，１３2のフィルタ係数が更新される。
＜実施形態に係る騒音制御システムの準備工程＞
図４を参照して（適宜、図１ないし図３を参照）、実施形態に係る騒音制御システム１，１０１の準備工程について説明する。図４は、騒音制御システム１，１０１の準備工程を示すフローチャートである。
図４に示すように、騒音制御システム１，１０１の準備工程は、主に、「Ｓ１：制御スピーカ・誤差マイクロホンの位置決定工程」、「Ｓ２：騒音源室・受音室内のモード周波数の計測工程」、「Ｓ３：騒音源の周波数特性の計測工程」、「Ｓ４：受音室の周波数特性の計測工程」、「Ｓ５：通過周波数の決定工程」、「Ｓ６：バンドパスフィルタの作成工程」、「Ｓ７：制御スピーカ－誤差マイクロホン間の伝達関数（誤差経路特性）の推定工程」を有する。なお、工程の種類によっては工程の順番を適宜入れ替えることが可能である。
図４を参照して（適宜、図１ないし図３を参照）、実施形態に係る騒音制御システム１，１０１の準備工程について説明する。図４は、騒音制御システム１，１０１の準備工程を示すフローチャートである。
図４に示すように、騒音制御システム１，１０１の準備工程は、主に、「Ｓ１：制御スピーカ・誤差マイクロホンの位置決定工程」、「Ｓ２：騒音源室・受音室内のモード周波数の計測工程」、「Ｓ３：騒音源の周波数特性の計測工程」、「Ｓ４：受音室の周波数特性の計測工程」、「Ｓ５：通過周波数の決定工程」、「Ｓ６：バンドパスフィルタの作成工程」、「Ｓ７：制御スピーカ－誤差マイクロホン間の伝達関数（誤差経路特性）の推定工程」を有する。なお、工程の種類によっては工程の順番を適宜入れ替えることが可能である。
（Ｓ１：制御スピーカ・誤差マイクロホンの位置決定工程）
制御スピーカ４および誤差マイクロホン３を設置する場所を決定する。制御スピーカ４および誤差マイクロホン３が設置されるのは、騒音源室Ｒ１（図２参照）内である。制御スピーカ４の設置箇所は、モード周波数を励起できる室隅とする。誤差信号を得る誤差マイクロホン３の設置箇所は、室隅（室隅はモード周波数の節とならない）かつ伝搬する室（受音室Ｒ２）に接する場所が望ましい。受音室Ｒ２に接する室隅であることで、受音室Ｒ２のモード周波数の制御も行うことが可能となる。
制御スピーカ４および誤差マイクロホン３を設置する場所を決定する。制御スピーカ４および誤差マイクロホン３が設置されるのは、騒音源室Ｒ１（図２参照）内である。制御スピーカ４の設置箇所は、モード周波数を励起できる室隅とする。誤差信号を得る誤差マイクロホン３の設置箇所は、室隅（室隅はモード周波数の節とならない）かつ伝搬する室（受音室Ｒ２）に接する場所が望ましい。受音室Ｒ２に接する室隅であることで、受音室Ｒ２のモード周波数の制御も行うことが可能となる。
（Ｓ２：騒音源室・受音室内のモード周波数の計測工程）
続いて、モード周波数成分を抽出するためのバンドパスフィルタ１１を設計するために、事前に騒音源室Ｒ１、受音室Ｒ２内のモード周波数を計測する。モード周波数の計測方法は特に限定されないが、例えば「適応フィルタ」を用いることでモード周波数を推定することができる。適応フィルタを用いてモード周波数を推定する場合、例えば図５に示す伝達関数推定システム５０を用いる。図５は、適応フィルタを用いた伝達関数推定システムの構成図である。
続いて、モード周波数成分を抽出するためのバンドパスフィルタ１１を設計するために、事前に騒音源室Ｒ１、受音室Ｒ２内のモード周波数を計測する。モード周波数の計測方法は特に限定されないが、例えば「適応フィルタ」を用いることでモード周波数を推定することができる。適応フィルタを用いてモード周波数を推定する場合、例えば図５に示す伝達関数推定システム５０を用いる。図５は、適応フィルタを用いた伝達関数推定システムの構成図である。
伝達関数推定システム５０は、主に、スピーカ５１と、マイクロホン５２と、シグナルジェネレータ５３と、適応フィルタ５４と、ＬＭＳ処理部５５とを備える。適応フィルタ５４は、例えばＦＩＲ（Finite Impulse Response）適応型デジタルフィルタである。
騒音源室Ｒ１内のモード周波数の計測について説明する。騒音源室Ｒ１の角にスピーカ５１、対向する角にマイクロホン５２を設置し、シグナルジェネレータ５３にてホワイトノイズ（白色雑音）を発生させる。ホワイトノイズは、スピーカ５１、適応フィルタ５４、ＬＭＳ処理部５５にそれぞれ入力する。スピーカ５１からはホワイトノイズが騒音源室Ｒ１内に放射され、マイクロホン５２まで伝搬する。適応フィルタ５４のフィルタ係数は、ＬＭＳアルゴリズムに基づいて更新される。このフィルタ係数を用いて畳み込み和演算のようなフィルタ処理により、フィルタ出力信号が計算される。この出力信号は、マイクロホン５２で得られた信号と足し合わされる。足し合わされた信号はＬＭＳ処理部５５に入力される。
ＬＭＳ処理部５５には、ホワイトノイズ、および適応フィルタ５４の出力信号とマイクロホン５２で得られた信号とを足し合わせた信号が入力される。ＬＭＳ処理部５５は、これらの信号を用いてスピーカ５１－マイクロホン５２間の伝達関数と等しくなるように適応フィルタ５４のフィルタ係数を更新する。更新して得られたフィルタ係数をフーリエ変換することで、スピーカ５１－マイクロホン５２間の伝達関数が得られる。騒音源室Ｒ１のモード周波数は、伝達関数で極大値を取っている周波数である。これにより、騒音源室Ｒ１内のモード周波数の計測が完了する。
同様の方法で、受音室Ｒ２内のモード周波数を計測する。つまり、受音室Ｒ２の角にスピーカ５１、対向する角にマイクロホン５２を設置し、シグナルジェネレータ５３にてホワイトノイズ（白色雑音）を発生させ、スピーカ５１－マイクロホン５２間の伝達関数と等しくなるように適応フィルタ５４のフィルタ係数を更新する。そして、更新して得られたフィルタ係数をフーリエ変換することでスピーカ５１－マイクロホン５２間の伝達関数が得られ、伝達関数で極大値を取っている周波数を受音室Ｒ２のモード周波数とする。これにより、受音室Ｒ２内のモード周波数の計測が完了する。
（Ｓ３：騒音源の周波数特性の計測工程）
図４に示すように、続いて、騒音源９の通常稼働時の周波数特性の計測を実施する。この計測は、例えば騒音源９の近くにマイクロホンを設置して行う。騒音源９で音が大きくなっている卓越周波数を求める。
図４に示すように、続いて、騒音源９の通常稼働時の周波数特性の計測を実施する。この計測は、例えば騒音源９の近くにマイクロホンを設置して行う。騒音源９で音が大きくなっている卓越周波数を求める。
（Ｓ４：受音室の周波数特性の計測工程）
「騒音源室・受音室内のモード周波数の計測工程」で計測したモード周波数による増幅が原因の騒音が受音室Ｒ２で既に発生している場合、騒音源９を通常稼働させた状態で受音室Ｒ２内の騒音の周波数特性を計測する。
「騒音源室・受音室内のモード周波数の計測工程」で計測したモード周波数による増幅が原因の騒音が受音室Ｒ２で既に発生している場合、騒音源９を通常稼働させた状態で受音室Ｒ２内の騒音の周波数特性を計測する。
（Ｓ５：通過周波数の決定工程）
工程（Ｓ２）、工程（Ｓ３）、工程（Ｓ４）の計測で得られた結果から、バンドパスフィルタ１１の作成に必要な通過周波数を決定する。具体的には、工程（Ｓ４）と工程（Ｓ２）の結果を比較し、工程（Ｓ４）の卓越周波数と工程（Ｓ２）のモード周波数とが一致する周波数を通過周波数とする。また、通過周波数と工程（Ｓ３）の卓越周波数を比較し、騒音源９から発生している騒音であるかを確認する。なお、騒音源９を設置前などの場合、工程（Ｓ２）のモード周波数をバンドパスフィルタ１１の通過周波数としてもよい。
工程（Ｓ２）、工程（Ｓ３）、工程（Ｓ４）の計測で得られた結果から、バンドパスフィルタ１１の作成に必要な通過周波数を決定する。具体的には、工程（Ｓ４）と工程（Ｓ２）の結果を比較し、工程（Ｓ４）の卓越周波数と工程（Ｓ２）のモード周波数とが一致する周波数を通過周波数とする。また、通過周波数と工程（Ｓ３）の卓越周波数を比較し、騒音源９から発生している騒音であるかを確認する。なお、騒音源９を設置前などの場合、工程（Ｓ２）のモード周波数をバンドパスフィルタ１１の通過周波数としてもよい。
（Ｓ６：バンドパスフィルタの作成工程）
工程（Ｓ５）で得られた通過周波数を基にしてバンドパスフィルタ１１を作成する。例えば、通過周波数を中心に所定の帯域幅（例えば、5Hz）を通過するようなチェビシェフフィルタを設計するのがよい。バンドパスフィルタ１１にチェビシェフフィルタを用いるのがよい理由は、フィルタの減衰特性が急峻で狭帯域を抽出しやすいためである。
工程（Ｓ５）で得られた通過周波数を基にしてバンドパスフィルタ１１を作成する。例えば、通過周波数を中心に所定の帯域幅（例えば、5Hz）を通過するようなチェビシェフフィルタを設計するのがよい。バンドパスフィルタ１１にチェビシェフフィルタを用いるのがよい理由は、フィルタの減衰特性が急峻で狭帯域を抽出しやすいためである。
（Ｓ７：制御スピーカ－誤差マイクロホン間の伝達関数（誤差経路特性）の推定工程）
続いて、アクティブノイズコントロールのフィルタ更新アルゴリズムであるFiltered-x LMSにおける誤差経路特性モデル１２や、MEFX-LMSアルゴリズムにおける誤差経路特性モデル１２11，１２21，１２12，１２22を作成するために、制御スピーカ－誤差マイクロホン間の伝達関数を推定する。推定法としては、例えば適応フィルタによる方法がある。適応フィルタを用いた伝達関数の推定法は、公知の技術であるので詳細な説明は省略する。この推定法の考え方は、例えば「特開平９－１７１３８８号公報」に記載されている。
続いて、アクティブノイズコントロールのフィルタ更新アルゴリズムであるFiltered-x LMSにおける誤差経路特性モデル１２や、MEFX-LMSアルゴリズムにおける誤差経路特性モデル１２11，１２21，１２12，１２22を作成するために、制御スピーカ－誤差マイクロホン間の伝達関数を推定する。推定法としては、例えば適応フィルタによる方法がある。適応フィルタを用いた伝達関数の推定法は、公知の技術であるので詳細な説明は省略する。この推定法の考え方は、例えば「特開平９－１７１３８８号公報」に記載されている。
そして、工程（Ｓ６）で求めたバンドパスフィルタ１１および工程（Ｓ７）で求めた伝達関数を用いて騒音制御システム１，１０１の構築する。以上で騒音制御システム１，１０１の準備工程が完了し、構築した騒音制御システム１，１０１を用いて実際に騒音制御を行う（制御工程）。
＜実施形態に係る騒音制御システムの効果＞
以上のように、本実施形態に係る騒音制御システム１，１０１は、騒音源室Ｒ１および受音室Ｒ２のモード周波数だけ通過するバンドパスフィルタ１１を用いることで参照信号の絞り込みを行う。これにより、騒音源９が放射する周波数と騒音源室Ｒ１および受音室Ｒ２のいずれかのモード周波数とが一致した際に生じる受音室Ｒ２での音圧の増幅を騒音源室Ｒ１内で抑えることができる。その為、受音室Ｒ２全体に渡って音圧低減が可能である。また、騒音制御システム１，１０１を構成する機器は、騒音源室Ｒ１に設置することができ、受音室Ｒ２にこれらの機器を設置しなくてよい。その為、受音室Ｒ２ではスペースの利用に制限を受けないので、スペースを有効に利用することができる。
以上のように、本実施形態に係る騒音制御システム１，１０１は、騒音源室Ｒ１および受音室Ｒ２のモード周波数だけ通過するバンドパスフィルタ１１を用いることで参照信号の絞り込みを行う。これにより、騒音源９が放射する周波数と騒音源室Ｒ１および受音室Ｒ２のいずれかのモード周波数とが一致した際に生じる受音室Ｒ２での音圧の増幅を騒音源室Ｒ１内で抑えることができる。その為、受音室Ｒ２全体に渡って音圧低減が可能である。また、騒音制御システム１，１０１を構成する機器は、騒音源室Ｒ１に設置することができ、受音室Ｒ２にこれらの機器を設置しなくてよい。その為、受音室Ｒ２ではスペースの利用に制限を受けないので、スペースを有効に利用することができる。
本実施形態に係る騒音制御システム１，１０１を用いた騒音制御のイメージを図６に示す。図６では、騒音源室Ｒ１として設備機械がある機械室を想定し、受音室Ｒ２として機械室に隣接する小規模会議室を想定している。図６に示すように、制御スピーカ４を用いた機械室内での騒音制御によって、機械室のみならず小規模会議室の音圧を低減させる。
本実施形態に係る騒音制御システム１，１０１の効果を検証するために実験を行ったので説明する。図７を参照して実験を行った環境について説明する。図７は、実験を行った室を説明するための図である。
図７に示すように、隣接する第一室Ｑ１および第二室Ｑ２を用いて実験を行った。第一室Ｑ１および第二室Ｑ２は、概ね矩形状を呈しており、第二室Ｑ２は第一室Ｑ１に比べて大きい。第一室Ｑ１および第二室Ｑ２内には柱等の障害物Ｍ１～Ｍ４が設けられている。その為、第一室Ｑ１および第二室Ｑ２のモード周波数を計算により求めるのは難しい。
図７に示すように、隣接する第一室Ｑ１および第二室Ｑ２を用いて実験を行った。第一室Ｑ１および第二室Ｑ２は、概ね矩形状を呈しており、第二室Ｑ２は第一室Ｑ１に比べて大きい。第一室Ｑ１および第二室Ｑ２内には柱等の障害物Ｍ１～Ｍ４が設けられている。その為、第一室Ｑ１および第二室Ｑ２のモード周波数を計算により求めるのは難しい。
第一室Ｑ１に騒音源９を模したスピーカ９ａを配置し、スピーカ９ａから騒音を第一室Ｑ１内に放出する。そのため、第一室Ｑ１が騒音源室Ｒ１に相当し、第二室Ｑ２が受音室Ｒ２に相当する。第一室Ｑ１の角（紙面右下の障害物Ｍ２の近く）に仮想の騒音源９としてのスピーカ９ａを配置した。なお、図７での細い破線は、音圧測定点（点線の交点）を示すために補助的に記載したものである。
図８を参照して、準備工程のモード周波数計測工程（Ｓ２）における機器の配置について説明する。図８は、実験でのモード周波数計測工程の機器配置を示す図である。図８に示すように、騒音源室Ｒ１としての第一室Ｑ１の角（紙面右下の障害物Ｍ２の近く）にスピーカ５１を配置し、その対向する角にマイクロホン５２を配置して計測を行った。また、受音室Ｒ２としての第二室Ｑ２の角（紙面左上の室隅）にスピーカ５１を配置し、その対向する角にマイクロホン５２を配置して計測を行った。
各室で実測したフィルタ係数をフーリエ変換したものを図９に示す。図９は、実験を行った室の伝達関数を示すグラフである。図９に示すグラフの横軸は周波数[Hz]であり、縦軸は振幅[dB]である。図９では、騒音源室Ｒ１としての第一室Ｑ１のグラフを実線で示し、受音室Ｒ２としての第二室Ｑ２のグラフを破線で示している。各室の70 [Hz]以下の周波数特性の極大値を求めると以下の通りであり、以上が実験で用いた各室のモード周波数である。
・第一室Ｑ１（騒音源室Ｒ１）：30,53,60,65 [Hz]
・第二室Ｑ２（受音室Ｒ２） ：30,40,60,69 [Hz]
・第一室Ｑ１（騒音源室Ｒ１）：30,53,60,65 [Hz]
・第二室Ｑ２（受音室Ｒ２） ：30,40,60,69 [Hz]
前述した通り、実験では仮想の騒音源９としてスピーカ９ａ（図７参照）を用いている。スピーカ９ａに入力する信号は、任意の純音(単一周波数成分をもつ音)とし、卓越周波数を自由に設定可能とした。実験では、スピーカ９ａに40,53 [Hz]の純音を与えた。スピーカ９ａに40 [Hz]の純音を与えた場合、スピーカ９ａから放出される騒音の卓越周波数は当然に40 [Hz]である。また、スピーカ９ａに53 [Hz]の純音を与えた場合、スピーカ９ａから放出される騒音の卓越周波数は当然に53 [Hz]である。40 [Hz]は、第二室Ｑ２のモード周波数の一つであり、53 [Hz]は、第一室Ｑ１のモード周波数の一つである。実験では、スピーカ９ａから40 [Hz]の純音が放射されているとき、第二室Ｑ２のモード周波数と一致することで第二室Ｑ２内の騒音の増幅が確認できた（準備工程の工程（Ｓ４））。
図１０を参照して、実験での通過周波数の決定工程（Ｓ５）を説明する。上述した通り、スピーカ９ａから出る騒音が40 [Hz]の純音の場合、図１０を参照すると40 [Hz]は第二室Ｑ２のモード周波数であり、40 [Hz]を通過周波数として設定することが可能である。また、スピーカ９ａから出る騒音が53 [Hz]の純音の場合、図１０を参照すると53 [Hz]は第一室Ｑ１のモード周波数であり、53 [Hz]を通過周波数として設定することが可能である。
続いて、実験でのバンドパスフィルタの作成工程（Ｓ６）を説明する。実験では、通過周波数が40 [Hz]に対応させて、低域通過帯域周波数が37.5 [Hz]であり、高域通過帯域周波数が 42.5 [Hz]であるチェビシェフバンドパスフィルタを設計した。作成したバンドパスフィルタを図１１に示す。また、通過周波数が53 [Hz]に対応させて、低域通過帯域周波数が 50.5 [Hz]であり、高域通過帯域周波数が 55.5 [Hz]であるチェビシェフバンドパスフィルタを設計した。作成したバンドパスフィルタを図１２に示す。
続いて、実験での伝達関数（誤差経路特性）の推定工程（Ｓ７）を説明する。実験では、図１３に示す機器配置で制御スピーカ－誤差マイクロホン間の伝達関数（誤差経路特性）を推定した。図１３は、実験での制御スピーカ－誤差マイクロホン間の伝達関数（誤差経路特性）の推定工程の機器配置を示す図である。図１３に示すように、実験では、二つの誤差マイクロホン３1，３2および二つの制御スピーカ４1，４2を用いたMEFX-LMSアルゴリズムを含むアクティブノイズコントロール技術を想定している。
実験での伝達関数（誤差経路特性）の推定工程（Ｓ７）により、誤差経路特性モデル１２11，１２21，１２12，１２22が作成される（図３参照）。誤差経路特性モデル１２11は、第一制御スピーカ４1－第一誤差マイクロホン３1間の伝達関数の推定モデルである。誤差経路特性モデル１２21は、第一制御スピーカ４1－第二誤差マイクロホン３2間の伝達関数の推定モデルである。誤差経路特性モデル１２12は、第二制御スピーカ４2－第一誤差マイクロホン３1間の伝達関数の推定モデルである。誤差経路特性モデル１２22は、第二制御スピーカ４2－第二誤差マイクロホン３2間の伝達関数の推定モデルである。
ここまでの準備工程の結果を基にして、実験では、図１４に示す騒音制御システム２０１を作成した。図１４に示す騒音制御システム２０１は、実施形態の騒音制御システム１０１（図３参照）に対応したものであるが、実験では図１４に示すように、シグナルジェネレータ９ｂとスピーカ９ａとの組合せを騒音源としている。つまり、実験では、騒音源としてシグナルジェネレータ９ｂで純音を生成してスピーカ９ａから放射した。これにより、単一周波数だけ抽出するバンドパスフィルタ１１ａ（図３参照）を模擬している。シグナルジェネレータ９ｂで生成する純音は、40 [Hz], 53 [Hz]の２種類である。なお、誤差マイクロホン３1，３2および制御スピーカ４1，４2の配置は図１３で示す配置と同じである。図１４に示す騒音制御システム２０１を動作させることで、誤差マイクロホン３1，３2における音圧が最小となるように適応フィルタ１３1，１３2が更新される。
図１５および図１６を参照して、騒音制御を行った場合における第二室Ｑ２（受音室Ｒ２）での制御効果について説明する。図１５は、第二室Ｑ２（受音室Ｒ２）における制御前後を比較した40 [Hz]の差分音圧分布である。図１６は、第二室Ｑ２（受音室Ｒ２）における制御前後を比較した53 [Hz]の差分音圧分布である。
図１５および図１６において、横軸は測定点の位置であり、縦軸は制御前後での音圧レベルの差分（「音圧レベル制御前後差分」と称する）である。測定点は、図１３中の破線（Ｘ１～Ｘ１０，Ｙ１～Ｙ７）が交差している点である。音圧レベル制御前後差分が負になると音圧が制御により低減していることを示す。
図１５および図１６において、横軸は測定点の位置であり、縦軸は制御前後での音圧レベルの差分（「音圧レベル制御前後差分」と称する）である。測定点は、図１３中の破線（Ｘ１～Ｘ１０，Ｙ１～Ｙ７）が交差している点である。音圧レベル制御前後差分が負になると音圧が制御により低減していることを示す。
図１５に示すように、シグナルジェネレータ９ｂで生成する純音が40 [Hz]の場合では、受音室Ｒ２である第二室Ｑ２全体の音圧が低減しており、制御効果が現れている。また、図１６に示すように、シグナルジェネレータ９ｂで生成する純音が53 [Hz]の場合では、受音室Ｒ２である第二室Ｑ２の広い範囲で音圧が低減しており、制御効果が現れている。なお、53 [Hz]の場合では、音圧レベル制御前後差分が正になる点も存在するが、全体として音圧が低減する傾向にあるといえる（例えば、全測定点の平均値が負になっており、また、正側の体積よりも負側の体積の大きくなっている）。なお、図示は省略しているが、モード周波数以外の周波数では音圧の増幅はない。
このように、各室（騒音源室、受音室）のモード周波数のみを抽出して制御に利用することで、受音室の広い範囲の音圧を低減できるアクティブノイズコントロールシステムを構築できることが実験結果からも証明された。
以上、本発明の実施形態について説明したが、本発明はこれに限定されるものではなく、特許請求の範囲の趣旨を変えない範囲で実施することができる。
１，１０１ 騒音制御システム
２ 参照マイクロホン
３，３1，３2 誤差マイクロホン
４，４1，４2 制御スピーカ
５，１０５ 信号処理装置
９ 騒音源
１１ａ，１１ｂ，１１ バンドパスフィルタ
１２，１２11，１２21，１２12，１２22 誤差経路特性モデル
１３，１３1，１３2 適応フィルタ
１４，１４1，１４2 ＬＭＳ処理部
１５ パワーアンプ
Ｒ１ 騒音源室
Ｒ２ 受音室
Ｒａ 壁
２ 参照マイクロホン
３，３1，３2 誤差マイクロホン
４，４1，４2 制御スピーカ
５，１０５ 信号処理装置
９ 騒音源
１１ａ，１１ｂ，１１ バンドパスフィルタ
１２，１２11，１２21，１２12，１２22 誤差経路特性モデル
１３，１３1，１３2 適応フィルタ
１４，１４1，１４2 ＬＭＳ処理部
１５ パワーアンプ
Ｒ１ 騒音源室
Ｒ２ 受音室
Ｒａ 壁
Claims (4)
- 騒音源を備える騒音源室に隣接する受音室での騒音を低減させる騒音制御システムであって、
前記騒音源室内に設置される参照マイクロホン、誤差マイクロホンおよび制御スピーカと、
信号処理を行うことで前記騒音源室内の騒音を打ち消す制御音を発生させる制御信号を前記制御スピーカに出力する信号処理装置と、を備え、
前記参照マイクロホンは、前記制御音によって打ち消される前の騒音を検出可能な場所に設置され、当該騒音を参照信号として前記信号処理装置に出力し、
前記制御スピーカは、前記騒音源室内の室隅に設置され、前記制御音を前記騒音源室内に放出し、
前記誤差マイクロホンは、前記騒音源室内の室隅に設置され、前記制御音によって打ち消された後の騒音を誤差信号として前記信号処理装置に出力し、
前記信号処理装置は、
前記参照信号および前記誤差信号から特定の周波数の成分を抽出するバンドパスフィルタと、
前記バンドパスフィルタによって抽出された前記参照信号および前記誤差信号に基づき、ＬＭＳアルゴリズムによってフィルタ係数が更新される適応フィルタと、を備え、
前記適応フィルタは、前記バンドパスフィルタによって抽出された前記参照信号を用いて前記制御信号を生成し、
前記バンドパスフィルタは、前記騒音源室で音が増幅しやすい第一モード周波数の成分および前記受音室で音が増幅しやすい第二モード周波数の成分を抽出するものである、
ことを特徴とする騒音制御システム。 - 前記信号処理装置は、前記制御スピーカと前記誤差マイクロホンとの間の伝達関数が組み込まれた誤差経路特性モデルをさらに有し、
前記適応フィルタのフィルタ係数は、前記適応フィルタと前記誤差経路特性モデルとを合成した伝達関数と、前記参照マイクロホンと前記誤差マイクロホンとの間の伝達関数とが等しくなるように更新される、
ことを特徴とする請求項１に記載の騒音制御システム。 - 騒音源を備える騒音源室に隣接する受音室での騒音を低減させる騒音制御方法であって、
参照マイクロホン、誤差マイクロホンおよび制御スピーカを前記騒音源室内に設置する準備工程と、
信号処理を行うことで前記騒音源室内の騒音を打ち消す制御音を発生させる制御信号を前記制御スピーカに出力する制御工程と、を有し、
前記準備工程では、
前記参照マイクロホンを前記制御音によって打ち消される前の騒音を検出可能な場所に設置し、
前記制御スピーカおよび前記誤差マイクロホンを前記騒音源室内の室隅に設置し、
前記制御工程には、
前記参照マイクロホンから取得した参照信号および前記誤差マイクロホンから取得した誤差信号から、バンドパスフィルタを用いて前記騒音源室で音が増幅しやすい第一モード周波数の成分および前記受音室で音が増幅しやすい第二モード周波数の成分を抽出する抽出工程と、
前記バンドパスフィルタによって抽出された前記参照信号および前記誤差信号に基づき、ＬＭＳアルゴリズムによって適応フィルタのフィルタ係数を更新し、また、前記適応フィルタが前記バンドパスフィルタによって抽出された前記参照信号を用いて前記制御信号を生成する制御信号生成工程と、が含まれる、
ことを特徴とする騒音制御方法。 - 前記準備工程には、
前記第一モード周波数および前記第二モード周波数を計測するモード周波数計測工程と、
前記騒音の中で音が大きい卓越周波数を前記騒音源室内で計測する周波数特性計測工程と、
前記第一モード周波数と前記卓越周波数とが一致する周波数を含む第一狭帯域、および前記第二モード周波数と前記卓越周波数とが一致する周波数を含む第二狭帯域を通過させる前記バンドパスフィルタを作成するバンドパスフィルタ作成工程と、が含まれる、
ことを特徴とする請求項３に記載の騒音制御方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2020160133A JP7444750B2 (ja) | 2020-09-24 | 2020-09-24 | 騒音制御システムおよび騒音制御方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2020160133A JP7444750B2 (ja) | 2020-09-24 | 2020-09-24 | 騒音制御システムおよび騒音制御方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2022053341A JP2022053341A (ja) | 2022-04-05 |
| JP7444750B2 true JP7444750B2 (ja) | 2024-03-06 |
Family
ID=80963219
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2020160133A Active JP7444750B2 (ja) | 2020-09-24 | 2020-09-24 | 騒音制御システムおよび騒音制御方法 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP7444750B2 (ja) |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2001147692A (ja) | 1999-11-19 | 2001-05-29 | Takahiro Yamashita | 騒音抑制装置及び音発生装置 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH11259078A (ja) * | 1998-03-16 | 1999-09-24 | Oki Electric Ind Co Ltd | 能動騒音制御装置 |
-
2020
- 2020-09-24 JP JP2020160133A patent/JP7444750B2/ja active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2001147692A (ja) | 1999-11-19 | 2001-05-29 | Takahiro Yamashita | 騒音抑制装置及び音発生装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| JP2022053341A (ja) | 2022-04-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Kuo et al. | Active noise control system for headphone applications | |
| WO2011118218A1 (ja) | スピーカ装置、音制御装置、スピーカ装置が取り付けられた壁 | |
| JP5562860B2 (ja) | 能動騒音制御装置 | |
| JP2009029405A (ja) | 騒音制御装置 | |
| US11190894B2 (en) | Active room compensation in loudspeaker system | |
| US9191768B2 (en) | Diffracted sound reduction device, diffracted sound reduction method, and filter coefficient determination method | |
| Qiu et al. | A new era for applications of active noise control | |
| Hasegawa et al. | Window active noise control system with virtual sensing technique | |
| JP7444750B2 (ja) | 騒音制御システムおよび騒音制御方法 | |
| JPH0511771A (ja) | 騒音制御装置 | |
| Takahashi et al. | Broadband active sound control system for air-conditioning duct noise | |
| Romeu et al. | Active noise control in ducts in presence of standing waves. Its influence on feedback effect | |
| WO2023021735A1 (ja) | センサモジュール、アクティブ制御装置、アクティブ制御方法およびプログラム | |
| Orlis et al. | Simulation of active low frequency noise control in closed spaces using the FDTD method | |
| García et al. | DSP implementation of the FxLMS algorithm for active noise control: Texas instruments TSM320C6713DSK | |
| KR101603697B1 (ko) | 능동 소음 제어를 이용한 층간 소음 저감 장치 및 방법 | |
| Hu et al. | Active noise control for selective cancellation of external disturbances | |
| de Diego et al. | On the performance of a local active noise control system | |
| Sun et al. | Experiments on performances of active-passive hybrid mufflers | |
| Zorzo et al. | Comparison between a digital and an analog active noise control system for headphones | |
| Michalczyk | Parameterization of LMS-based control algorithms for local zones of quiet | |
| Höber et al. | Generating a Position-Adaptive Quiet Zone in Enclosed Spaces | |
| Kojima et al. | Virtual-Sensing for Active Noise Control Using Reflected Waves | |
| JP2002328681A (ja) | 能動的騒音制御装置 | |
| Hanselka et al. | Converting an Active Noise Blocker for a Tilted Window from Feedforward Control into a Feedback System |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20230214 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20231027 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20231031 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20231226 |
|
| 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: 20240220 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20240222 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 7444750 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |