消音システムにおける騒音検出手段の位置決定方法
【課題】三次元閉空間に配置される騒音発生源から発生する騒音を良好に消音するのに、ＡＮＣ技術を充分に生かせる第２騒音検出手段の位置決定方法を得る。
【解決手段】ＡＮＣ技術を地下室等の三次元閉空間で適用するに、騒音発生源から発生する騒音のレベルが高い高音圧レベル周波数を複数求め、その複数の高音圧レベル周波数の音それぞれについて、三次元閉空間内において波動音響解析を実行し、三次元閉空間の音響状態を推定し、定在波の腹部で、三次元閉空間の音圧レベルが、空間内の平均音圧レベルより高い位置を前記第２騒音検出手段を配置する第２位置とする。
【解決手段】ＡＮＣ技術を地下室等の三次元閉空間で適用するに、騒音発生源から発生する騒音のレベルが高い高音圧レベル周波数を複数求め、その複数の高音圧レベル周波数の音それぞれについて、三次元閉空間内において波動音響解析を実行し、三次元閉空間の音響状態を推定し、定在波の腹部で、三次元閉空間の音圧レベルが、空間内の平均音圧レベルより高い位置を前記第２騒音検出手段を配置する第２位置とする。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、騒音発生源から発生する騒音を、第１位置に配置される第１騒音検出手段で検出するとともに、前記第１位置とは異なった第２位置に配置される第２騒音検出手段により検出し、前記第１騒音検出手段により検出される第１騒音情報及び前記第２騒音検出手段により検出される第２騒音情報に基づいて２次音源から制御音を発生して少なくとも前記第２位置で消音する消音システムにおける、前記２次音源又は第２騒音検出手段の位置決定方法に関し、さらに詳細には、三次元閉空間内で共鳴によって増大することがある低周波騒音を低減する消音技術に関する。
【背景技術】
【０００２】
低周波騒音（１００Ｈｚ以下の騒音）は吸音、遮音の効果が小さく、吸音処理、遮音処理ではその消音が困難である。これに対して、逆位相の音を２次音源（逆位相音発生用音源）から発生させることで、低周波音を効果的に低減させるＡＮＣ(Active Noise Control)と呼ばれる技術がある。この技術はヘッドホン、自動車の車内などの密閉空間、排気煙道などの一次元的な伝搬経路の騒音を低減するのに実用化されている。
【０００３】
三次元空間では騒音発生源と２次音源との位置関係が三次元的関係となり、両者の位置が大きくずれるため騒音低減効果は低下する場合が多いが、低周波騒音に限定すれば効果が期待できる。
【０００４】
特許文献１には、ＡＮＣ技術を三次元空間へ適用する事例が紹介されている。
この特許文献１に開示の技術は、円筒形のボイラーの外壁に備えられる送風機から前方に向けて伝播する低周波騒音をボイラーの外部前方で低減することを目的とする。この目的を達成するために、第１騒音検出手段としてのリファレンスマイクで騒音を検出するとともに、第２騒音検出手段としてのエラーマイクを送風機の前方に配置して、その位置でも騒音を検出し、これら両方の検出情報に基づいて、消音のための制御音を送風機の略鉛直上方又は下方に配置されたスピーカから送出して、エラーマイクの配置位置（本願の第２位置に相当する）で消音する。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】２００５−１３４６３２号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
特許文献１に開示の技術は、送風機の開口前方の三次元開放空間において消音を行うものである。
しかしながら、現実には、ＡＮＣの適用を考えた場合、例えば、地下室等の三次元閉空間に空調室外機の原動部が配置される場合等があり、このような事例の場合、その地下室内における騒音を低減するとともに、通気口等から地下室外に漏れる騒音を低減する必要がある。
この種の三次元閉空間に騒音発生源が配置された場合、騒音発生源から発生する騒音の特性が問題となるとともに、その三次元閉空間の音響特性も大きく影響する。結果、第２騒音検出手段の設置位置が不適切であると、充分な消音効果を得ることができない。
【０００７】
本発明の目的は、三次元閉空間に配置される騒音発生源から発生する騒音を良好に消音するのに、ＡＮＣ技術を充分に生かせる２次音源又は第２騒音検出手段の位置決定方法を得ることにある。
【課題を解決するための手段】
【０００８】
上記目的を達成するための、騒音発生源から発生する騒音を、第１位置に配置される第１騒音検出手段で検出するとともに、前記第１位置とは異なった第２位置に配置される第２騒音検出手段により検出し、前記第１騒音検出手段により検出される第１騒音情報及び前記第２騒音検出手段により検出される第２騒音情報に基づいて２次音源から制御音を発生して少なくとも前記第２位置で消音する消音システムにおける、前記２次音源又は第２騒音検出手段の位置決定方法の特徴構成は、
前記騒音発生源が三次元閉空間内に配置される場合に、
三次元開放空間において前記騒音発生源から発生する騒音について、低周波領域において音圧レベルが高い周波数である高音圧レベル周波数を求める第１ステップと、
前記第１ステップで求められる高音圧レベル周波数の音について、前記騒音発生源を配置した状態で、前記三次元閉空間内において、前記高音圧レベル周波数の音が前記騒音発生源から発生した場合の音響状態を波動音響解析により求める第２ステップとを実行し、
求められる前記三次元閉空間の音響状態において、定在波の腹部で、前記三次元閉空間の音圧レベルが、空間内の平均音圧レベルより高い位置を前記２次音源又は前記第２騒音検出手段を配置する第２位置とすることにある。
【０００９】
本願に係る騒音検出手段の位置決定方法は、特に三次元空間が閉空間である場合に、その威力を発揮する。即ち、三次元空間が閉空間である場合は、騒音発生源から発生する騒音について、三次元閉空間の音響特性とも相まって、その音響レベルが大きく変動する腹部及び殆ど変化しない節部が形成される。この様な腹部及び節部は、当然、三次元的な分布を取る。さらに、このような腹部及び節部は、特定の周波数に関して顕著に現れる。この周波数は、主には、騒音発生源の音圧レベルが高い周波数となる。
そこで、本願に係る位置決定方法では、第１ステップで、三次元開放空間での騒音発生源からの騒音の低周波領域の音圧レベルの高い周波数（高音圧レベル周波数）を求める。
ここで、周波数を低周波数に限定しているのは、消音対象とする音の周波数が、低周波領域の音であるためである。
【００１０】
第２ステップでは、このようにして得られる高音圧レベル周波数について、騒音発生源を配置した状態における三次元閉空間内で高音圧レベル周波数の音が騒音発生源から発生した場合の音響状態を波動音響解析で求める。このようにして問題となる可能性が高い周波数音の三次元閉空間内での腹部と節部の形成状況を具体的に把握する。
【００１１】
求められる三次元閉空間の音響状態において、定在波の腹部で、三次元閉空間の音圧レベルが、空間内の平均音圧レベルより高い位置を、２次音源又は第２騒音検出手段或いはそれらの両方を配置する第２位置とする。
【００１２】
このようにすることで、比較的簡易な解析手法を採用しながら、所謂、ＡＮＣを使用して、これまで容易に的確な騒音制御を実現することが困難であった三次元閉空間での消音を実現できるようになった。
【００１３】
即ち、ＡＮＣ技術の効果を好適にするためには、２次音源の位置、エラーマイクである第２騒音検出手段の位置（第２位置）或いはそれらの両方が好適化される必要がある。開放空間であれば騒音発生源付近の任意の場所に例えば、第２騒音検出手段を配置すればよいが、共鳴が起こるような密閉度の高い三次元閉空間においては、配置場所によってその効果が大幅に変わる。即ち、音波の干渉により音が小さくなる定在波の節部に例えば第２騒音検出手段が配置されると検出信号出力が小さくなり、相対的に暗騒音の影響を大きく受けることになる。結果、信号のＳ／Ｎ比が下がり、２次音源と第２騒音検出手段間の伝達関数の同定、消音動作にノイズが混じることになる。従って、消音効果が低下、制御が発散する可能性が高まる。逆に干渉によって音が強められる定在波の腹部では、Ｓ／Ｎ比が高くなるため、効果的な消音が可能になるのである。
よって、本願にあっては、このような三次元閉空間における腹部、節部の確認を、第１ステップ、第２ステップを経て容易且つ確実に行い、２次音源の位置或いは第２騒音検出手段の位置（第２位置）を好適化することで、三次元閉空間を対象としてＡＮＣを最適な状態で適用できるようになった。
【００１４】
このようにして、２次音源、第２騒音検出手段或いはそれら両方の位置を決定するに、以下のように決定するのが好ましい。
即ち、２次音源を配置する２次音源配置腹部と、第２騒音検出手段を配置する第２騒音検出手段配置腹部とを異なった腹部とすることが好ましい。
本願に係る消音システムでは、第２騒音検出手段により２次音源から発生される音による消音効果を確認することとなる。即ち、第２騒音検出手段がエラーマイクロフォンとして働く。このような騒音の検出構造において、三次元閉空間内に形成される複数の腹部に関して、２次音源と第２騒音検出手段とが同一の腹部内に位置されていると、両者が近接し、しかも同一の腹部内に存在するため、第２騒音検出手段において検出される検出音（消音後の音）に関して、その検出音に２次音源からの音（制御音）の影響が大きく残り、本来の消音効果を正確に代表できない場合がある。そこで、２次音源を配置する２次音源配置腹部と、第２騒音検出手段を配置する第２騒音検出手段配置腹部とを、三次元閉空間内に形成される異なった腹部とすることで、真に制御音による消音効果が作用した後の音を第２騒音検出手段により検出して、効果的な消音を実現できる。
さらに、このような構成を採る場合に、２次音源と第２騒音検出手段とを対向した位置関係とすると、第２騒音検出手段配置腹部で２次音源の効果を確実に確認できる。
【００１５】
さらに、三次元閉空間が平面視方形の空間であり、２次音源と第２騒音検出手段とを方形の長辺に沿う位置とすると、平面視方形の三次元閉空間において２次音源から第２騒音検出手段までの距離を確保して、平面視方形の三次元閉空間に於ける消音を有効に実現できる。
【００１６】
さらに、前記第１ステップにおいて、高音圧レベル周波数を音圧レベルに従って限られた数求め、
前記第２ステップにおいて、前記第１ステップで限られた数の求められる高音圧レベル周波数の音それぞれについて、前記騒音発生源を配置した状態で、前記三次元閉空間内において、前記各高音圧レベル周波数の音が前記騒音発生源から発生した場合の音響状態を波動音響解析により求め、
前記第２ステップで求められた、限られた数の高音圧レベル周波数の音に関する前記音響状態に基づいて、前記低周波数領域での前記限られた数の高音圧レベル周波数を含む前記三次元閉空間の音響状態を推定する第３ステップを実行し、
前記第３ステップで求められる前記三次元閉空間の音響状態において、定在波の腹部で、前記三次元閉空間の音圧レベルが空間内の平均音圧レベルより高い位置を前記２次音源又は第２騒音検出手段を配置する位置とすることが好ましい。
【００１７】
通常の騒音発生源から発生する騒音は、その音圧レベルが高い周波数が複数存在する。一方、対象とする周波数を連続的なものとすると、第２ステップで実行する解析が困難、あるいは解析に長時間を要する。そこで、第１ステップで求める高音圧レベル周波数の数は限られた数とする。通常、この種の周波数の数は、２〜５程度で充分である。
そして、第３ステップでは、得られた音響特性を、解析の対象とした周波数に限って合成する。その合成された結果で、低周波数領域での三次元閉空間の音響状態を代表できるものとする。
そして、このようにして第３ステップで求められた音響状態において、その腹部に２次音源、第２騒音検出手段或いはそれらの両方の位置を定めることができる。
【００１８】
さて、上記の位置決定方法では、騒音発生源と２次音源の位置が決まっている場合に第２騒音検出手段の位置である第２位置を決定する方法について説明したが、高音発生源の位置のみが判明しており、２次音源の位置及び第２騒音検出手段の位置（第２位置）を決定する場合に関しては、以下の手順で解析することが好ましい。
即ち、前記第１ステップ及び前記第２ステップを実行し、
前記第２ステップで求められた音響状態に関し、定在波の腹部を前記２次音源の位置と決定するとともに、当該決定された位置に前記２次音源を配置した状態で前記三次元閉空間内において、前記各高音圧レベル周波数の音が前記騒音発生源から発生した場合の音響状態を波動音響解析により求める２次音源位置決定・解析ステップを実行し、
前記２次音源位置決定・解析ステップで得られた音響状態において、定在波の腹部を前記第２騒音検出手段の位置である前記第２位置とする。
このようにすることで、消音対象としている三次元閉空間に関して、消音対象とする騒音の周波数、三次元閉空間内における騒音発生源の位置が特定される状態で、音響解析を行い、その解析により騒音発生源が存在する状態での定在波の形態を特定できる。そして、この定在波の腹部に２次音源の位置を決定することで、騒音発生源から発生する騒音に対して２次音源から発生する２次音による消音効果を有効に発揮させることができる。
次に、前記三次元閉空間に関して、消音対象とする騒音の周波数、三次元閉空間内における騒音発生源の位置及び２次音源の位置が特定される状態で、再度音響解析を実行し当該三次元閉空間の音響状態を求め、得られる定在波の腹の位置を２騒音検出手段の位置である第２位置とすることで、騒音発生源から発生する騒音に対する２次音源から発生する制御音の影響を的確に検出して、効果的な消音を実現できる。
【００１９】
さらに、前記第２騒音検出手段の位置としては、第３ステップで求められる三次元閉空間の音響状態において、定在波の腹部で、三次元閉空間の音圧レベルが最高音圧レベルとなる位置を前記第２位置とすることが好ましい。
このように最高音圧レベルとなる部位を第２騒音検出手段として選択することで、最も良好な消音効果を得ることができる。
【００２０】
さらに、前記第２ステップの波動音響解析において、基礎方程式をヘルムホルツ方程式とするとともに、離散化された境界積分方程式に従って、前記三次元閉空間における各部の音圧レベルを演算することが好ましい。
波動音響解析においては、定常状態を想定して波動方程式を特定の周波数で変数分離したヘルムホルツ方程式を基礎方程式とすることで、容易に三次元空間の解析を実行できる。
さらに、分割された各要素の境界条件に基づいて離散化された境界積分方程式を解くことで、音場の音圧が演算される。
【００２１】
さらに、前記騒音発生源が回転機器である場合に、前記第１ステップにおいて、当該回転機器から発生する騒音に関して、前記音圧レベルが高い基本周波数と当該基本周波数の倍調周波数（２倍調、３倍調等を含む）を、前記高音圧レベル周波数とすることが好ましい。
回転機器の場合、その音圧レベルが高い周波数が、回転機器の通常運転周波数に依存して特定される。さらに、騒音として大きく影響する周波数は、この基本周波数と倍調周波数であるため、これらに対して消音対策を施せば充分だからである。従って、第１ステップを大幅に簡略化できる。
【図面の簡単な説明】
【００２２】
【図１】第１実施形態が対象とする地下室内に配置された空調室外機の原動部の状態を示す図
【図２】消音システムの構成を示す図
【図３】第１実施形態の第２騒音検出手段であるエラーマイクの配置位置の決定フローを示す図
【図４】第１実施形態の原動部から発生する騒音の低周波数領域での音圧分布を示す図
【図５】境界積分方程式の各変数の幾何学的関係を示す図
【図６】離散化された境界積分方程式の各変数の幾何学的関係を示す図
【図７】第１実施形態の波動音響解析で使用した空調室外機の原動部が配置された地下室を代表する格子モデルを示す図
【図８】第１実施形態の２６Ｈｚでの音圧レベルの高いエリア及び水平断面音圧レベル分布の一例を示す図
【図９】第１実施形態の５２Ｈｚでの音圧レベルの高いエリア及び水平断面音圧レベル分布の一例を示す図
【図１０】第１実施形態の７８Ｈｚでの音圧レベルの高いエリア及び水平断面音圧レベル分布の一例を示す図
【図１１】第１実施形態のエラーマイクの配置位置決定の使用した第３ステップで推定した音圧レベルの高いエリア及び水平断面音圧レベル分布の一例を示す図
【図１２】第２実施形態及び第３実施形態の２次音源及び第２騒音検出手段であるエラーマイクの配置位置の決定フローを示す図
【図１３】第２実施形態の地下室の状態を示す図
【図１４】第２実施形態の地下室における騒音源からの音の音響状態を示す図
【図１５】第２実施形態の地下室における２次音源からの音の音響状態を示す図
【図１６】第２実施形態の地下室における消音効果を示す図
【図１７】第２実施形態の地下室における周波数領域における消音効果を示す図
【図１８】第３実施形態の地下室の状態を示す図
【図１９】第３実施形態の地下室における騒音源からの音の音響状態を示す図
【図２０】第３実施形態の地下室における消音効果を示す図
【図２１】第３実施形態の地下室における周波数領域における消音効果を示す図
【発明を実施するための形態】
【００２３】
本願に係る２次音源、第２騒音検出手段或いはそれらの両方の位置決定方法に関して、以下図面に基づいて説明する。
本願では、第一実施形態、第二実施形態、第三実施形態を紹介する。
ここで、第一実施形態は、三次元閉空間内での騒音源、２次音源の位置が定まっている場合において、第２騒音検出手段の位置を決定する方法に関し、第二、第三実施形態は、三次元閉空間内での騒音源の位置が定まっている場合において、２次音源、第２騒音検出手段の位置を決定する方法に関する。
以下の説明では、先ず、第一実施形態について、１．本願に係る消音システムが配置される配置場所の状況、２．消音システムの構成、３．騒音検出手段の位置決定方法の順に説明し、次に、第二、第三実施形態について、１．本願に係る消音システムが配置される配置場所の状況、３．騒音検出手段の位置決定方法の順に説明する。全ての実施形態において、消音システムの構成は、第一実施形態のものと同一であるため、第二、第三実施形態では説明を省略する。
【００２４】
〔第一実施形態〕
１．本願に係る消音システムが配置される配置場所の状況
図１は、本願に係る消音システムが配置される三次元閉空間の一例としての地下室１の状態を示したものである。地下室１は図面からも判明するように、概略方形の閉空間として構成されている。ここで、本願における三次元閉空間の概念は、例えば、通気口といった、地下室１と外空間２とを繋ぐ音響の一部が漏れる程度の開口部３が設けられている構成も含む概念である。本願では、三次元閉空間内で音響的な腹部及び節部が形成されるのを問題とするため、三次元的な腹部、節部が形成される程度の閉鎖性があれば、本願の対象範囲に含まれる。
図１に示すように、地下室１内には、騒音発生源の一例として空調室外機の原動部４（具体的にはガスエンジン）が配置される。同図に示すように、このガスエンジン４もまた、ほぼ方形をしており、地下室１内のほぼ３５％程度の空間を占有する。従って、図１において、原動部４の左下、右上及び左前側の側部に騒音が問題となる周部空間５が存在する。さらに、原動部４の天井側にも騒音が問題となる天井側空間６が広がっている。この天井側空間６は、地下室１の平断面積分の広がりをもっている。
この原動部４において、その騒音は、図１に騒音発生源７として記載する左下の壁部から発生するものと想定している。この想定は、実際の原動部４の筐体に設けられた通気孔に基づいた合理的な想定である。
【００２５】
２．消音システムの構成
消音システム８は、第１騒音検出手段としてのリファレンスマイク９、第２騒音検出手段としてのエラーマイク１０、２次音源としての一対の消音用スピーカ１１と、リファレンスマイク９及びエラーマイク１０により検出される騒音情報に基づいて、消音用スピーカ１１から消音のための制御音を発生するための音響制御情報を生成する制御情報生成装置１２（制御ボード）とから構成されている。
【００２６】
前記リファレンスマイク９は、図１に示すように、上記騒音発生源７の近く直前方に配置する。従って、このリファレンスマイク９は、空調室外機の原動部４から発生し、地下室１内の空間に放出される騒音を、そのまま検出する。
【００２７】
前記エラーマイク１０は、上記騒音発生源７より壁１３側に配置する。このエラーマイク１０は、基本的に消音システム８による消音が良好に働いているか否かの検証用のマイクである。本願では、このエラーマイク１０の配置位置が、独特の手法により決定される。
【００２８】
一対の消音用スピーカ１１は、制御情報生成装置１２から制御音を発生するための制御情報を受取り制御音を発生する。一対の消音用スピーカ１１は、騒音発生源７の近傍に騒音発生源７をその両脇から挟む位置関係で、騒音発生源７から発生する騒音の伝播方向（図１において左下向き）に向けて位置させる。
【００２９】
制御情報生成装置１２は、リファレンスマイク９、エラーマイク１０及び消音用スピーカ１１と電気的に接続されており、リファレンスマイク９及びエラーマイク１０から、それぞれ検出情報を受取る。一方、一対の消音用スピーカ１１には、この制御情報生成装置１２内に設けられている制御ボード１２ａにより、リファレンスマイク９及びエラーマイク１０の検出情報に基づいて決定される制御音を消音用スピーカ１１から発生するための情報を送る。
【００３０】
制御ボード１２ａの働きに関して説明すると、本願で採用するＦｉｌｔｅｒｅｄ−ＸＬＭＳアルゴリズムを採用する消音システム８は、図２（ａ）の構成を取っており、ｘnは第１位置に配置される第１騒音検出手段としてのリファレンスマイク９により検出されるリファレンス信号であり、ｅnは、第２位置に配置される第２騒音検出手段としてのエラーマイク１０に検出される誤差信号である。
そして、一次経路（この経路の伝達関数をＰ（Ｚ）と記載している）を経て伝達される騒音と、二次経路（この経路の伝達関数をＣ（Ｚ）と記載している）を経て伝達される制御音とが、エラーマイク１０が位置する第２位置に到達して残留する騒音である誤差が０となった場合に、消音が良好に行われたこととなる。
【００３１】
従って、この消音システム８では、図２（ｂ）で示すように、リファレンスマイク９により、騒音発生源７からの騒音波形を検出し、この騒音であるリファレンス信号ｘnを基にエラーマイク１０における誤差信号ｅnである音が小さくなるように、リファレンス信号ｘnとは逆位相の制御信号ｙnを生成する。この制御信号ｙnをアンプ１２ｂで増大させ、消音用スピーカ１１から発生させる。また、地下室１に配置される本消音システム８では高周波の騒音は逆に増大させてしまう可能性があるため、制御ボード１２ａ内にローパスフィルタ１２ｃを備えている。
【００３２】
図２（ａ）（ｂ）を参照して、この消音システム８の原理を説明すると、騒音が伝搬する一次経路と制御音源を経由する二次経路の関係は、以下の数１となる。
【００３３】
【数１】
【００３４】
ここで、事前に二次経路の特性を模擬したフィルタが二次経路を精度よく表しているとき、
【００３５】
【数２】
【００３６】
となり、数１、数２から、
【００３７】
【数３】
【００３８】
とできる。従って、数３が成り立つように騒音制御フィルタ（このフィルタの伝達関数をＨ（Ｚ）と記載している）を最適化すればよい。
騒音制御フィルタを求めるためにＬＭＳアルゴリズムは、数４で示される。
【００３９】
【数４】
【００４０】
ここで、μはステップサイズゲインである。
【００４１】
離散化した状態での適応制御は、図２（ｃ）の形態となる。この図では、騒音制御フィルタを実際に使用するＦＩＲフィルタ（Finite Impulse Responseフィルタ（有限インパルス応答フィルタ））と記載し、二次経路フィルタと係数更新アルゴリズムを適応制御アルゴリズムと記載している。
図２（ｃ）から判明するように、本消音システム８では、リファレンス信号を主入力信号ｘ[n]とし、この主入力信号ｘ[n]にフィルタ係数ｈk[n]をかけた出力信号ｙ[n]を生成し、その出力信号ｙ[n]と主入力信号ｘ[n]との合成信号（具体的には和信号で、この信号が誤差信号ε[n]である）が０となるようなｈk[n]が生成される。
【００４２】
ＦＩＲフィルタは、具体的には、入力信号ｘ〔ｎ−ｋ〕にフィルタ係数ｈk〔ｎ〕をかけた結果ｙ〔ｎ〕を出力するフィルタである。
【００４３】
【数５】
【００４４】
さらに、誤差信号ε[n]を最小にするためのフィルタ係数の作成は、公知のＬＭＳ適応制御アルゴリズムを採用することで実行できるが、このアルゴリズムは以下の数６の形態となる。
【００４５】
【数６】
ここで、μはステップサイズパラメータであり、下記に示す例では０．００００２程度を採用した。さらに、タップ数Ｍは、具体的には１０２４とした。
【００４６】
３．騒音検出手段の位置決定方法
本願に係る第２騒音検出手段の配置位置検定方法では、図３に示すように、第１〜第４ステップを経て、第２騒音検出手段であるエラーマイク１０の配置位置を決定する。
【００４７】
第１ステップ
第１ステップでは、三次元開放空間において騒音発生源から発生する騒音に関して、低周波領域において音圧レベルが高い周波数である高音圧レベル周波数を、音圧レベルに従って複数求める。
図４は、先に説明した空調室外機の原動部（騒音発生源７）４から発生する騒音を、開放空間で検出した結果である。空調室外機などの回転機械では、回転により発生する騒音でピーク周波数が決まり、例に示す機器の騒音は２６Ｈｚ（基本周波数）、５２Ｈｚ（２倍調周波数），７８Ｈｚ（３倍調周波数）にピークを持っていた。これら周波数が、本願にいうところの高音圧レベル周波数である。本願で消音に対象とする低周波領域は１００Ｈｚ以下の領域である。
【００４８】
第２ステップ
第２ステップでは、第１ステップで求められる複数の高音圧レベル周波数の音それぞれについて、騒音発生源を配置した状態における三次元閉空間内での、各高音圧レベル周波数の音が前記三次元空間内で発生した場合の音響状態を波動音響解析で求める。
このような波動音響解析には境界要素法を使用する。境界要素法は、波動音響学的に音場を解析するための数値解析手法であり、一般的には定常状態を想定し、波動方程式を特定の周波数で変数分離したヘルムホルツ方程式を基礎方程式とする。そして、ヘルムホルツ方程式を境界上の境界積分方程式の形態に変換した上で、要素に離散化することで音場を解析する。また、境界要素法は、有限要素法と違い、系に単位力が作用したときの系の応答を示す基本解を必要とするため、基本解が知られていない問題には適用できない。
上記の波動音響解析に使用する波動方程式を数７に、ヘルムホルツ方程式を数８に示した。さらに、このヘルムホルツ方程式に対応する境界積分方程式を数９に示した。ここで、この境界積分方程式で使用する三次元ヘルムホルツ方程式の基本解Ｇは数１０に示す形態となる。図５に、境界積分方程式の各変数の幾何学的関係を示した。
【００４９】
【数７】
【００５０】
【数８】
【００５１】
境界積分方程式
【数９】
Ｖ：法線方向粒子速度、ρ：空気密度、Ｃ：立体角
Ｐa:図５のａにおける音圧、Ｐb:図５のｂにおける音圧、Ｇ（ｒa,b）：図５のａｂ間の距離ｒa,bを変数とする基本解、Ｖb：図５における法線方向粒子速度
【００５２】
【数１０】
【００５３】
そして、数９、数１０に使用して境界上で離散化した境界積分方程式は数１１となる。
【００５４】
境界積分方程式の離散化
【数１１】
Ｃ：ｉにおける立体角を４πで割った値
領域の内部においてＣ＝１、滑らかな境界においてＣ＝１／２、領域の外側でＣ＝０となる。
【００５５】
また、図６に、離散化された境界積分方程式の各変数の幾何学的関係を示した。
数１１に示す方程式に基づいた解析に際しては、ＦＯＲＴＲＡＮで記載されたプログラムを使用した。
【００５６】
上述の離散化された境界積分方程式の適用対象とした格子モデルを、図７に示した。
解析に於ける境界の要素分割は、境界要素法では、差分法や有限要素法のように、空間を要素分割する必要はなく、境界のみを要素分割でよい。それら要素の長さは、波動性を適切に計算するために、波長の１／８以下に分割するのが妥当である。そのため、波長の短い高周波では、細かく要素分割する必要があるため、計算負荷が増大する。本願のように、低周波数領域を対象とする解析では、要素の長さは比較的長い合理的な長さとできる。実際の計算では、解析対象とした周波数領域が１００Ｈｚ以下であるため、単一の境界要素の長さは４００ｍｍとした。
【００５７】
また、解析条件として、対象に合わせて境界条件を設定する必要がある。境界条件としては、剛壁を想定した完全反射面、音源を想定した振動面、吸音面を想定したインピーダンス境界がある。実際の計算では、共鳴が問題になる空間は、音を吸収しにくいコンクリートなどの硬い壁が多いため、完全反射面かわずかな吸音を与えるインピーダンス境界を用いた。
【００５８】
図８、図９、図１０には、それぞれ、上記の２６Ｈｚ、５２Ｈｚ，７８Ｈｚについて、高音圧レベル周波数の音圧レベルが高いエリア（左側の図で、閉領域として囲まれた領域）、水平断面音圧レベル分布の一例（右側の図で、音響特性の節部を濃部で示す）を示した。周波数によっては、騒音発生源の側よりも遠方の方が大きな騒音となる場合があることが判る。また、腹部及び節部は三次元的に分布していることが分かる。
【００５９】
第３ステップ
この第３ステップでは、第２ステップで求められた、全ての高音圧レベル周波数の音に関する音響状態に基づいて、低周波数領域での全ての高音圧レベル周波数を含む三次元空間の音響状態を推定する。２６Ｈｚ、５２Ｈｚ，７８Ｈｚの結果を使用して、低周波数領域（２０〜１００Ｈｚ）の騒音として合成（具体的には、単純な合算）した場合での音圧レベルが高いエリア、水平断面音圧レベル分布の一例を図１１に示した。この図から判明するように、地下室の場合、図１１右側のエリア図で、右奥側に示される側壁１３の下部に近い部位に、第２騒音検出手段を配置すべきエリアが広がっていることが分かる。
【００６０】
第４ステップ
この第４ステップでは、第３ステップで求められる三次元空間の音響状態において、定在波の腹部で、前記三次元空間の音圧レベルが平均音圧レベルより高い位置を第２騒音検出手段の配置位置とする。
これまで説明してきた図１１に示した例の場合、空調室外機の配置側（図１の左下側）で壁面近傍の最下部位置が好ましいことが判明した。
このように、騒音が大きい位置にエラーマイクを配置することでＳ／Ｎ比の高い信号が得られ、効果的な消音が可能となる。また、特に騒音を低減したい場所、エラーマイクが劣化しにくい場所なども考慮してエラーマイクの配置位置を選定することが好ましい。
【００６１】
〔第二実施形態〕
１．本願に係る消音システムが配置される配置場所の状況
図１３は、この実施形態に係る消音システムが配置される三次元閉空間の一例としての地下室１の状態を示したものである。地下室１は図面からも判明するように、概略方形の閉空間として構成されている。
図１３に示すように、地下室１の平面視短辺側の壁面上部には、騒音発生源の一例として換気扇４０が配置される。
３．２次音源及び第２騒音検出手段の位置決定方法
この実施形態に係る２次音源及び第２騒音検出手段の配置位置検定方法では、図１２に示すように、第１ステップ、第２ステップ（騒音源からの音による音響状態の導出）、２次音源位置決定・解析ステップ（２次音源の位置決定、及び、決定された２次音源からの音による音響状態の導出）、第２騒音検出手段位置決定ステップ（第２騒音検出手段の位置決定）を実行するのである。
【００６２】
第１ステップ（＃１０）
第１ステップでは、三次元開放空間において騒音発生源から発生する騒音に関して、低周波領域において音圧レベルが高い周波数である高音圧レベル周波数を、音圧レベルに従って求める。
先に説明した換気扇（騒音発生源）４０では、４１Ｈｚ（基本周波数）にピークを持っていた。この周波数が、本願にいうところの高音圧レベル周波数である。本願で消音に対象とする低周波領域は１００Ｈｚ以下の領域である。
【００６３】
第２ステップ（＃２０）
第２ステップでは、第一実施形態と同様に、第１ステップで求められる高音圧レベル周波数の音について、騒音発生源を配置した状態における三次元閉空間内での、各高音圧レベル周波数の音が三次元空間内で発生した場合の音響状態を波動音響解析で求める。解析手法は、先の形態と同様であるので、説明は省略する。
【００６４】
図１４に、上記の４１Ｈｚについての結果を示した。
左側の図面(ａ)が音圧レベルの高い領域（方形四隅にある囲まれた３次元の領域）を示し、右側の図面（ｂ）が音圧分布で、音圧の高い部分を暗い領域として濃度で示した。
これらの図から、平面視方形の空間では、方形の四隅に腹部が形成され、方形の中央側に十字形に節部が形成されていることが判る。また、腹部及び節部は三次元的に分布していることが分かる。
【００６５】
２次音源位置決定・解析ステップ（＃３０）
この２次音源位置決定・解析ステップでは、第２ステップで求められた、高音圧レベル周波数の音に関する音響状態に基づいて、スピーカ１１（２次音源）の位置を、音響状態で腹部となっている部位に決定する。図１３から判明するように、この例では、２次音源の位置を、平面視方形となる三次元閉空間において、騒音源である換気扇４０が位置される短辺側に沿った、騒音源とは反対側の角とした。この部位は、求められた音響状態において音響の腹部の一つの位置である。
【００６６】
このようにして決定した２次音源を配置した状態における三次元閉空間内での、高音圧レベル周波数の音が三次元空間内で発生した場合の音響状態を波動音響解析で求める。解析手法は、先の形態と同様であるので、説明は省略する。
図１５に、上記の４１Ｈｚについての結果を示した。
左側の図面(ａ)が音圧レベルの高い領域（方形四隅にある囲まれた３次元の領域）を示し、右側の図面（ｂ）が音圧分布で、音圧の高い部分を暗い領域として濃度で示した。
これらの図から、平面視方形の空間では、方形の四隅に腹部が形成され、方形の中央側に十字形に節部が形成されていることが判る。
【００６７】
第２騒音検出手段位置決定ステップ（＃４０）
この第２騒音検出手段位置決定ステップでは、２次音源位置決定・解析ステップで求められる三次元空間の音響状態において、定在波の腹部で、三次元空間の音圧レベルが平均音圧レベルより高い位置をエラーマイク１０（第２騒音検出手段）の配置位置とする。
これまで説明してきた図１３に示した例の場合、２次音源に対向した位置が、好適な位置となった。図１６が、消音効果を示した音圧分布であり、図１７が周波数領域における消音効果を示したものである。目的とした４１Ｈｚで消音効果が得られている。
このように、騒音が大きい位置にエラーマイクを配置することでＳ／Ｎ比の高い信号が得られ、効果的な消音が可能となる。
【００６８】
〔第三実施形態〕
１．本願に係る消音システムが配置される配置場所の状況
図１８は、この実施形態に係る消音システムが配置される三次元閉空間の一例としての地下室１の状態を示したものである。地下室１は図面からも判明するように、概略方形の閉空間として構成されている。
図１８に示すように、地下室１の平面視短辺側の壁面上部には、騒音発生源の一例として一対の換気扇４１が配置される。
３．２次音源及び第２騒音検出手段の位置決定方法
この実施形態に係る２次音源及び第２騒音検出手段の配置位置検定方法でも、図１２に示すように、第１ステップ、第２ステップ（騒音源からの音による音響状態の導出）、２次音源位置決定・解析ステップ（２次音源の位置決定、及び、決定された２次音源からの音による音響状態の導出）、第２騒音検出手段位置決定ステップ（第２騒音検出手段の位置決定）を実行する。
【００６９】
第１ステップ（＃１０）
第１ステップでは、三次元開放空間において騒音発生源から発生する騒音に関して、低周波領域において音圧レベルが高い周波数である高音圧レベル周波数を、音圧レベルに従って求める。
先に説明した換気扇（騒音発生源）４１でも、４１Ｈｚ（基本周波数）にピークを持っていた。これら周波数が、本願にいうところの高音圧レベル周波数である。本願で消音に対象とする低周波領域は１００Ｈｚ以下の領域である。
【００７０】
第２ステップ
第２ステップでは、第一実施形態と同様に、第１ステップで求められる高音圧レベル周波数の音について、騒音発生源を配置した状態における三次元閉空間内での、各高音圧レベル周波数の音が前記三次元空間内で発生した場合の音響状態を波動音響解析で求める。解析手法は、先の形態と同様であるので、説明は省略する。
【００７１】
図１９に、上記の４１Ｈｚについての結果を示した。
左側の図面(ａ)が音圧レベルの高い領域（方形四隅にある囲まれた３次元の領域）を示し、右側の図面（ｂ）が音圧分布で、音圧の高い部分を暗い領域として濃度で示した。
これらの図から、平面視方形の空間では、方形の四隅に腹部が形成され、方形の中央側に十字形に節部が形成されていることが判る。また、腹部及び節部は三次元的に分布していることが分かる。
【００７２】
２次音源位置決定・解析ステップ（＃３０）
この２次音源位置決定・解析ステップでは、第２ステップで求められた、高音圧レベル周波数の音に関する音響状態に基づいて、スピーカ１１（２次音源）の位置を、音響状態で腹部となっている部位に２次音源の位置を決定する。この例でも、２次音源の位置を、平面視方形となる三次元閉空間において、騒音源である一方の換気扇４１が位置される短辺側に沿った騒音源側の角とした。この部位は、求められた音響状態において音響の腹部の一つの位置である。
【００７３】
このようにして決定した２次音源を配置した状態における三次元閉空間内での、高音圧レベル周波数の音が三次元空間内で発生した場合の音響状態を波動音響解析で求める。解析手法は、先の形態と同様であるので、説明は省略する。
図１５と同様な結果が得られた。
【００７４】
第２騒音検出手段位置決定ステップ（＃４０）
この第２騒音検出手段位置決定ステップでは、２次音源位置決定・解析ステップで求められる三次元空間の音響状態において、定在波の腹部で、三次元空間の音圧レベルが平均音圧レベルより高い位置をエラーマイク１０（第２騒音検出手段）の配置位置とする。
これまで説明してきた図１８に示した例の場合、２次音源に対向した位置が、好適な位置となった。図２０が、消音効果を示した音圧分布であり、図２１が周波数領域における消音効果を示したものである。目的とした４１Ｈｚで消音効果が得られている。
このように、騒音が大きい位置にエラーマイクを配置することでＳ／Ｎ比の高い信号が得られ、効果的な消音が可能となる。
【産業上の利用可能性】
【００７５】
三次元閉空間に配置される騒音発生源から発生する騒音を良好に消音するのに、ＡＮＣ技術を充分に生かせる第２騒音検出手段の位置決定方法を得ることができた。
【符号の説明】
【００７６】
１ 地下室（三次元閉空間）
４ 空調室外機の原動部
９ リファレンスマイク（第１騒音検出手段）
１０ エラーマイク（第２騒音検出手段）
１１ 消音用スピーカ（２次音源）
１２ 制御情報生成装置（制御ボード）
【技術分野】
【０００１】
本発明は、騒音発生源から発生する騒音を、第１位置に配置される第１騒音検出手段で検出するとともに、前記第１位置とは異なった第２位置に配置される第２騒音検出手段により検出し、前記第１騒音検出手段により検出される第１騒音情報及び前記第２騒音検出手段により検出される第２騒音情報に基づいて２次音源から制御音を発生して少なくとも前記第２位置で消音する消音システムにおける、前記２次音源又は第２騒音検出手段の位置決定方法に関し、さらに詳細には、三次元閉空間内で共鳴によって増大することがある低周波騒音を低減する消音技術に関する。
【背景技術】
【０００２】
低周波騒音（１００Ｈｚ以下の騒音）は吸音、遮音の効果が小さく、吸音処理、遮音処理ではその消音が困難である。これに対して、逆位相の音を２次音源（逆位相音発生用音源）から発生させることで、低周波音を効果的に低減させるＡＮＣ(Active Noise Control)と呼ばれる技術がある。この技術はヘッドホン、自動車の車内などの密閉空間、排気煙道などの一次元的な伝搬経路の騒音を低減するのに実用化されている。
【０００３】
三次元空間では騒音発生源と２次音源との位置関係が三次元的関係となり、両者の位置が大きくずれるため騒音低減効果は低下する場合が多いが、低周波騒音に限定すれば効果が期待できる。
【０００４】
特許文献１には、ＡＮＣ技術を三次元空間へ適用する事例が紹介されている。
この特許文献１に開示の技術は、円筒形のボイラーの外壁に備えられる送風機から前方に向けて伝播する低周波騒音をボイラーの外部前方で低減することを目的とする。この目的を達成するために、第１騒音検出手段としてのリファレンスマイクで騒音を検出するとともに、第２騒音検出手段としてのエラーマイクを送風機の前方に配置して、その位置でも騒音を検出し、これら両方の検出情報に基づいて、消音のための制御音を送風機の略鉛直上方又は下方に配置されたスピーカから送出して、エラーマイクの配置位置（本願の第２位置に相当する）で消音する。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】２００５−１３４６３２号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
特許文献１に開示の技術は、送風機の開口前方の三次元開放空間において消音を行うものである。
しかしながら、現実には、ＡＮＣの適用を考えた場合、例えば、地下室等の三次元閉空間に空調室外機の原動部が配置される場合等があり、このような事例の場合、その地下室内における騒音を低減するとともに、通気口等から地下室外に漏れる騒音を低減する必要がある。
この種の三次元閉空間に騒音発生源が配置された場合、騒音発生源から発生する騒音の特性が問題となるとともに、その三次元閉空間の音響特性も大きく影響する。結果、第２騒音検出手段の設置位置が不適切であると、充分な消音効果を得ることができない。
【０００７】
本発明の目的は、三次元閉空間に配置される騒音発生源から発生する騒音を良好に消音するのに、ＡＮＣ技術を充分に生かせる２次音源又は第２騒音検出手段の位置決定方法を得ることにある。
【課題を解決するための手段】
【０００８】
上記目的を達成するための、騒音発生源から発生する騒音を、第１位置に配置される第１騒音検出手段で検出するとともに、前記第１位置とは異なった第２位置に配置される第２騒音検出手段により検出し、前記第１騒音検出手段により検出される第１騒音情報及び前記第２騒音検出手段により検出される第２騒音情報に基づいて２次音源から制御音を発生して少なくとも前記第２位置で消音する消音システムにおける、前記２次音源又は第２騒音検出手段の位置決定方法の特徴構成は、
前記騒音発生源が三次元閉空間内に配置される場合に、
三次元開放空間において前記騒音発生源から発生する騒音について、低周波領域において音圧レベルが高い周波数である高音圧レベル周波数を求める第１ステップと、
前記第１ステップで求められる高音圧レベル周波数の音について、前記騒音発生源を配置した状態で、前記三次元閉空間内において、前記高音圧レベル周波数の音が前記騒音発生源から発生した場合の音響状態を波動音響解析により求める第２ステップとを実行し、
求められる前記三次元閉空間の音響状態において、定在波の腹部で、前記三次元閉空間の音圧レベルが、空間内の平均音圧レベルより高い位置を前記２次音源又は前記第２騒音検出手段を配置する第２位置とすることにある。
【０００９】
本願に係る騒音検出手段の位置決定方法は、特に三次元空間が閉空間である場合に、その威力を発揮する。即ち、三次元空間が閉空間である場合は、騒音発生源から発生する騒音について、三次元閉空間の音響特性とも相まって、その音響レベルが大きく変動する腹部及び殆ど変化しない節部が形成される。この様な腹部及び節部は、当然、三次元的な分布を取る。さらに、このような腹部及び節部は、特定の周波数に関して顕著に現れる。この周波数は、主には、騒音発生源の音圧レベルが高い周波数となる。
そこで、本願に係る位置決定方法では、第１ステップで、三次元開放空間での騒音発生源からの騒音の低周波領域の音圧レベルの高い周波数（高音圧レベル周波数）を求める。
ここで、周波数を低周波数に限定しているのは、消音対象とする音の周波数が、低周波領域の音であるためである。
【００１０】
第２ステップでは、このようにして得られる高音圧レベル周波数について、騒音発生源を配置した状態における三次元閉空間内で高音圧レベル周波数の音が騒音発生源から発生した場合の音響状態を波動音響解析で求める。このようにして問題となる可能性が高い周波数音の三次元閉空間内での腹部と節部の形成状況を具体的に把握する。
【００１１】
求められる三次元閉空間の音響状態において、定在波の腹部で、三次元閉空間の音圧レベルが、空間内の平均音圧レベルより高い位置を、２次音源又は第２騒音検出手段或いはそれらの両方を配置する第２位置とする。
【００１２】
このようにすることで、比較的簡易な解析手法を採用しながら、所謂、ＡＮＣを使用して、これまで容易に的確な騒音制御を実現することが困難であった三次元閉空間での消音を実現できるようになった。
【００１３】
即ち、ＡＮＣ技術の効果を好適にするためには、２次音源の位置、エラーマイクである第２騒音検出手段の位置（第２位置）或いはそれらの両方が好適化される必要がある。開放空間であれば騒音発生源付近の任意の場所に例えば、第２騒音検出手段を配置すればよいが、共鳴が起こるような密閉度の高い三次元閉空間においては、配置場所によってその効果が大幅に変わる。即ち、音波の干渉により音が小さくなる定在波の節部に例えば第２騒音検出手段が配置されると検出信号出力が小さくなり、相対的に暗騒音の影響を大きく受けることになる。結果、信号のＳ／Ｎ比が下がり、２次音源と第２騒音検出手段間の伝達関数の同定、消音動作にノイズが混じることになる。従って、消音効果が低下、制御が発散する可能性が高まる。逆に干渉によって音が強められる定在波の腹部では、Ｓ／Ｎ比が高くなるため、効果的な消音が可能になるのである。
よって、本願にあっては、このような三次元閉空間における腹部、節部の確認を、第１ステップ、第２ステップを経て容易且つ確実に行い、２次音源の位置或いは第２騒音検出手段の位置（第２位置）を好適化することで、三次元閉空間を対象としてＡＮＣを最適な状態で適用できるようになった。
【００１４】
このようにして、２次音源、第２騒音検出手段或いはそれら両方の位置を決定するに、以下のように決定するのが好ましい。
即ち、２次音源を配置する２次音源配置腹部と、第２騒音検出手段を配置する第２騒音検出手段配置腹部とを異なった腹部とすることが好ましい。
本願に係る消音システムでは、第２騒音検出手段により２次音源から発生される音による消音効果を確認することとなる。即ち、第２騒音検出手段がエラーマイクロフォンとして働く。このような騒音の検出構造において、三次元閉空間内に形成される複数の腹部に関して、２次音源と第２騒音検出手段とが同一の腹部内に位置されていると、両者が近接し、しかも同一の腹部内に存在するため、第２騒音検出手段において検出される検出音（消音後の音）に関して、その検出音に２次音源からの音（制御音）の影響が大きく残り、本来の消音効果を正確に代表できない場合がある。そこで、２次音源を配置する２次音源配置腹部と、第２騒音検出手段を配置する第２騒音検出手段配置腹部とを、三次元閉空間内に形成される異なった腹部とすることで、真に制御音による消音効果が作用した後の音を第２騒音検出手段により検出して、効果的な消音を実現できる。
さらに、このような構成を採る場合に、２次音源と第２騒音検出手段とを対向した位置関係とすると、第２騒音検出手段配置腹部で２次音源の効果を確実に確認できる。
【００１５】
さらに、三次元閉空間が平面視方形の空間であり、２次音源と第２騒音検出手段とを方形の長辺に沿う位置とすると、平面視方形の三次元閉空間において２次音源から第２騒音検出手段までの距離を確保して、平面視方形の三次元閉空間に於ける消音を有効に実現できる。
【００１６】
さらに、前記第１ステップにおいて、高音圧レベル周波数を音圧レベルに従って限られた数求め、
前記第２ステップにおいて、前記第１ステップで限られた数の求められる高音圧レベル周波数の音それぞれについて、前記騒音発生源を配置した状態で、前記三次元閉空間内において、前記各高音圧レベル周波数の音が前記騒音発生源から発生した場合の音響状態を波動音響解析により求め、
前記第２ステップで求められた、限られた数の高音圧レベル周波数の音に関する前記音響状態に基づいて、前記低周波数領域での前記限られた数の高音圧レベル周波数を含む前記三次元閉空間の音響状態を推定する第３ステップを実行し、
前記第３ステップで求められる前記三次元閉空間の音響状態において、定在波の腹部で、前記三次元閉空間の音圧レベルが空間内の平均音圧レベルより高い位置を前記２次音源又は第２騒音検出手段を配置する位置とすることが好ましい。
【００１７】
通常の騒音発生源から発生する騒音は、その音圧レベルが高い周波数が複数存在する。一方、対象とする周波数を連続的なものとすると、第２ステップで実行する解析が困難、あるいは解析に長時間を要する。そこで、第１ステップで求める高音圧レベル周波数の数は限られた数とする。通常、この種の周波数の数は、２〜５程度で充分である。
そして、第３ステップでは、得られた音響特性を、解析の対象とした周波数に限って合成する。その合成された結果で、低周波数領域での三次元閉空間の音響状態を代表できるものとする。
そして、このようにして第３ステップで求められた音響状態において、その腹部に２次音源、第２騒音検出手段或いはそれらの両方の位置を定めることができる。
【００１８】
さて、上記の位置決定方法では、騒音発生源と２次音源の位置が決まっている場合に第２騒音検出手段の位置である第２位置を決定する方法について説明したが、高音発生源の位置のみが判明しており、２次音源の位置及び第２騒音検出手段の位置（第２位置）を決定する場合に関しては、以下の手順で解析することが好ましい。
即ち、前記第１ステップ及び前記第２ステップを実行し、
前記第２ステップで求められた音響状態に関し、定在波の腹部を前記２次音源の位置と決定するとともに、当該決定された位置に前記２次音源を配置した状態で前記三次元閉空間内において、前記各高音圧レベル周波数の音が前記騒音発生源から発生した場合の音響状態を波動音響解析により求める２次音源位置決定・解析ステップを実行し、
前記２次音源位置決定・解析ステップで得られた音響状態において、定在波の腹部を前記第２騒音検出手段の位置である前記第２位置とする。
このようにすることで、消音対象としている三次元閉空間に関して、消音対象とする騒音の周波数、三次元閉空間内における騒音発生源の位置が特定される状態で、音響解析を行い、その解析により騒音発生源が存在する状態での定在波の形態を特定できる。そして、この定在波の腹部に２次音源の位置を決定することで、騒音発生源から発生する騒音に対して２次音源から発生する２次音による消音効果を有効に発揮させることができる。
次に、前記三次元閉空間に関して、消音対象とする騒音の周波数、三次元閉空間内における騒音発生源の位置及び２次音源の位置が特定される状態で、再度音響解析を実行し当該三次元閉空間の音響状態を求め、得られる定在波の腹の位置を２騒音検出手段の位置である第２位置とすることで、騒音発生源から発生する騒音に対する２次音源から発生する制御音の影響を的確に検出して、効果的な消音を実現できる。
【００１９】
さらに、前記第２騒音検出手段の位置としては、第３ステップで求められる三次元閉空間の音響状態において、定在波の腹部で、三次元閉空間の音圧レベルが最高音圧レベルとなる位置を前記第２位置とすることが好ましい。
このように最高音圧レベルとなる部位を第２騒音検出手段として選択することで、最も良好な消音効果を得ることができる。
【００２０】
さらに、前記第２ステップの波動音響解析において、基礎方程式をヘルムホルツ方程式とするとともに、離散化された境界積分方程式に従って、前記三次元閉空間における各部の音圧レベルを演算することが好ましい。
波動音響解析においては、定常状態を想定して波動方程式を特定の周波数で変数分離したヘルムホルツ方程式を基礎方程式とすることで、容易に三次元空間の解析を実行できる。
さらに、分割された各要素の境界条件に基づいて離散化された境界積分方程式を解くことで、音場の音圧が演算される。
【００２１】
さらに、前記騒音発生源が回転機器である場合に、前記第１ステップにおいて、当該回転機器から発生する騒音に関して、前記音圧レベルが高い基本周波数と当該基本周波数の倍調周波数（２倍調、３倍調等を含む）を、前記高音圧レベル周波数とすることが好ましい。
回転機器の場合、その音圧レベルが高い周波数が、回転機器の通常運転周波数に依存して特定される。さらに、騒音として大きく影響する周波数は、この基本周波数と倍調周波数であるため、これらに対して消音対策を施せば充分だからである。従って、第１ステップを大幅に簡略化できる。
【図面の簡単な説明】
【００２２】
【図１】第１実施形態が対象とする地下室内に配置された空調室外機の原動部の状態を示す図
【図２】消音システムの構成を示す図
【図３】第１実施形態の第２騒音検出手段であるエラーマイクの配置位置の決定フローを示す図
【図４】第１実施形態の原動部から発生する騒音の低周波数領域での音圧分布を示す図
【図５】境界積分方程式の各変数の幾何学的関係を示す図
【図６】離散化された境界積分方程式の各変数の幾何学的関係を示す図
【図７】第１実施形態の波動音響解析で使用した空調室外機の原動部が配置された地下室を代表する格子モデルを示す図
【図８】第１実施形態の２６Ｈｚでの音圧レベルの高いエリア及び水平断面音圧レベル分布の一例を示す図
【図９】第１実施形態の５２Ｈｚでの音圧レベルの高いエリア及び水平断面音圧レベル分布の一例を示す図
【図１０】第１実施形態の７８Ｈｚでの音圧レベルの高いエリア及び水平断面音圧レベル分布の一例を示す図
【図１１】第１実施形態のエラーマイクの配置位置決定の使用した第３ステップで推定した音圧レベルの高いエリア及び水平断面音圧レベル分布の一例を示す図
【図１２】第２実施形態及び第３実施形態の２次音源及び第２騒音検出手段であるエラーマイクの配置位置の決定フローを示す図
【図１３】第２実施形態の地下室の状態を示す図
【図１４】第２実施形態の地下室における騒音源からの音の音響状態を示す図
【図１５】第２実施形態の地下室における２次音源からの音の音響状態を示す図
【図１６】第２実施形態の地下室における消音効果を示す図
【図１７】第２実施形態の地下室における周波数領域における消音効果を示す図
【図１８】第３実施形態の地下室の状態を示す図
【図１９】第３実施形態の地下室における騒音源からの音の音響状態を示す図
【図２０】第３実施形態の地下室における消音効果を示す図
【図２１】第３実施形態の地下室における周波数領域における消音効果を示す図
【発明を実施するための形態】
【００２３】
本願に係る２次音源、第２騒音検出手段或いはそれらの両方の位置決定方法に関して、以下図面に基づいて説明する。
本願では、第一実施形態、第二実施形態、第三実施形態を紹介する。
ここで、第一実施形態は、三次元閉空間内での騒音源、２次音源の位置が定まっている場合において、第２騒音検出手段の位置を決定する方法に関し、第二、第三実施形態は、三次元閉空間内での騒音源の位置が定まっている場合において、２次音源、第２騒音検出手段の位置を決定する方法に関する。
以下の説明では、先ず、第一実施形態について、１．本願に係る消音システムが配置される配置場所の状況、２．消音システムの構成、３．騒音検出手段の位置決定方法の順に説明し、次に、第二、第三実施形態について、１．本願に係る消音システムが配置される配置場所の状況、３．騒音検出手段の位置決定方法の順に説明する。全ての実施形態において、消音システムの構成は、第一実施形態のものと同一であるため、第二、第三実施形態では説明を省略する。
【００２４】
〔第一実施形態〕
１．本願に係る消音システムが配置される配置場所の状況
図１は、本願に係る消音システムが配置される三次元閉空間の一例としての地下室１の状態を示したものである。地下室１は図面からも判明するように、概略方形の閉空間として構成されている。ここで、本願における三次元閉空間の概念は、例えば、通気口といった、地下室１と外空間２とを繋ぐ音響の一部が漏れる程度の開口部３が設けられている構成も含む概念である。本願では、三次元閉空間内で音響的な腹部及び節部が形成されるのを問題とするため、三次元的な腹部、節部が形成される程度の閉鎖性があれば、本願の対象範囲に含まれる。
図１に示すように、地下室１内には、騒音発生源の一例として空調室外機の原動部４（具体的にはガスエンジン）が配置される。同図に示すように、このガスエンジン４もまた、ほぼ方形をしており、地下室１内のほぼ３５％程度の空間を占有する。従って、図１において、原動部４の左下、右上及び左前側の側部に騒音が問題となる周部空間５が存在する。さらに、原動部４の天井側にも騒音が問題となる天井側空間６が広がっている。この天井側空間６は、地下室１の平断面積分の広がりをもっている。
この原動部４において、その騒音は、図１に騒音発生源７として記載する左下の壁部から発生するものと想定している。この想定は、実際の原動部４の筐体に設けられた通気孔に基づいた合理的な想定である。
【００２５】
２．消音システムの構成
消音システム８は、第１騒音検出手段としてのリファレンスマイク９、第２騒音検出手段としてのエラーマイク１０、２次音源としての一対の消音用スピーカ１１と、リファレンスマイク９及びエラーマイク１０により検出される騒音情報に基づいて、消音用スピーカ１１から消音のための制御音を発生するための音響制御情報を生成する制御情報生成装置１２（制御ボード）とから構成されている。
【００２６】
前記リファレンスマイク９は、図１に示すように、上記騒音発生源７の近く直前方に配置する。従って、このリファレンスマイク９は、空調室外機の原動部４から発生し、地下室１内の空間に放出される騒音を、そのまま検出する。
【００２７】
前記エラーマイク１０は、上記騒音発生源７より壁１３側に配置する。このエラーマイク１０は、基本的に消音システム８による消音が良好に働いているか否かの検証用のマイクである。本願では、このエラーマイク１０の配置位置が、独特の手法により決定される。
【００２８】
一対の消音用スピーカ１１は、制御情報生成装置１２から制御音を発生するための制御情報を受取り制御音を発生する。一対の消音用スピーカ１１は、騒音発生源７の近傍に騒音発生源７をその両脇から挟む位置関係で、騒音発生源７から発生する騒音の伝播方向（図１において左下向き）に向けて位置させる。
【００２９】
制御情報生成装置１２は、リファレンスマイク９、エラーマイク１０及び消音用スピーカ１１と電気的に接続されており、リファレンスマイク９及びエラーマイク１０から、それぞれ検出情報を受取る。一方、一対の消音用スピーカ１１には、この制御情報生成装置１２内に設けられている制御ボード１２ａにより、リファレンスマイク９及びエラーマイク１０の検出情報に基づいて決定される制御音を消音用スピーカ１１から発生するための情報を送る。
【００３０】
制御ボード１２ａの働きに関して説明すると、本願で採用するＦｉｌｔｅｒｅｄ−ＸＬＭＳアルゴリズムを採用する消音システム８は、図２（ａ）の構成を取っており、ｘnは第１位置に配置される第１騒音検出手段としてのリファレンスマイク９により検出されるリファレンス信号であり、ｅnは、第２位置に配置される第２騒音検出手段としてのエラーマイク１０に検出される誤差信号である。
そして、一次経路（この経路の伝達関数をＰ（Ｚ）と記載している）を経て伝達される騒音と、二次経路（この経路の伝達関数をＣ（Ｚ）と記載している）を経て伝達される制御音とが、エラーマイク１０が位置する第２位置に到達して残留する騒音である誤差が０となった場合に、消音が良好に行われたこととなる。
【００３１】
従って、この消音システム８では、図２（ｂ）で示すように、リファレンスマイク９により、騒音発生源７からの騒音波形を検出し、この騒音であるリファレンス信号ｘnを基にエラーマイク１０における誤差信号ｅnである音が小さくなるように、リファレンス信号ｘnとは逆位相の制御信号ｙnを生成する。この制御信号ｙnをアンプ１２ｂで増大させ、消音用スピーカ１１から発生させる。また、地下室１に配置される本消音システム８では高周波の騒音は逆に増大させてしまう可能性があるため、制御ボード１２ａ内にローパスフィルタ１２ｃを備えている。
【００３２】
図２（ａ）（ｂ）を参照して、この消音システム８の原理を説明すると、騒音が伝搬する一次経路と制御音源を経由する二次経路の関係は、以下の数１となる。
【００３３】
【数１】
【００３４】
ここで、事前に二次経路の特性を模擬したフィルタが二次経路を精度よく表しているとき、
【００３５】
【数２】
【００３６】
となり、数１、数２から、
【００３７】
【数３】
【００３８】
とできる。従って、数３が成り立つように騒音制御フィルタ（このフィルタの伝達関数をＨ（Ｚ）と記載している）を最適化すればよい。
騒音制御フィルタを求めるためにＬＭＳアルゴリズムは、数４で示される。
【００３９】
【数４】
【００４０】
ここで、μはステップサイズゲインである。
【００４１】
離散化した状態での適応制御は、図２（ｃ）の形態となる。この図では、騒音制御フィルタを実際に使用するＦＩＲフィルタ（Finite Impulse Responseフィルタ（有限インパルス応答フィルタ））と記載し、二次経路フィルタと係数更新アルゴリズムを適応制御アルゴリズムと記載している。
図２（ｃ）から判明するように、本消音システム８では、リファレンス信号を主入力信号ｘ[n]とし、この主入力信号ｘ[n]にフィルタ係数ｈk[n]をかけた出力信号ｙ[n]を生成し、その出力信号ｙ[n]と主入力信号ｘ[n]との合成信号（具体的には和信号で、この信号が誤差信号ε[n]である）が０となるようなｈk[n]が生成される。
【００４２】
ＦＩＲフィルタは、具体的には、入力信号ｘ〔ｎ−ｋ〕にフィルタ係数ｈk〔ｎ〕をかけた結果ｙ〔ｎ〕を出力するフィルタである。
【００４３】
【数５】
【００４４】
さらに、誤差信号ε[n]を最小にするためのフィルタ係数の作成は、公知のＬＭＳ適応制御アルゴリズムを採用することで実行できるが、このアルゴリズムは以下の数６の形態となる。
【００４５】
【数６】
ここで、μはステップサイズパラメータであり、下記に示す例では０．００００２程度を採用した。さらに、タップ数Ｍは、具体的には１０２４とした。
【００４６】
３．騒音検出手段の位置決定方法
本願に係る第２騒音検出手段の配置位置検定方法では、図３に示すように、第１〜第４ステップを経て、第２騒音検出手段であるエラーマイク１０の配置位置を決定する。
【００４７】
第１ステップ
第１ステップでは、三次元開放空間において騒音発生源から発生する騒音に関して、低周波領域において音圧レベルが高い周波数である高音圧レベル周波数を、音圧レベルに従って複数求める。
図４は、先に説明した空調室外機の原動部（騒音発生源７）４から発生する騒音を、開放空間で検出した結果である。空調室外機などの回転機械では、回転により発生する騒音でピーク周波数が決まり、例に示す機器の騒音は２６Ｈｚ（基本周波数）、５２Ｈｚ（２倍調周波数），７８Ｈｚ（３倍調周波数）にピークを持っていた。これら周波数が、本願にいうところの高音圧レベル周波数である。本願で消音に対象とする低周波領域は１００Ｈｚ以下の領域である。
【００４８】
第２ステップ
第２ステップでは、第１ステップで求められる複数の高音圧レベル周波数の音それぞれについて、騒音発生源を配置した状態における三次元閉空間内での、各高音圧レベル周波数の音が前記三次元空間内で発生した場合の音響状態を波動音響解析で求める。
このような波動音響解析には境界要素法を使用する。境界要素法は、波動音響学的に音場を解析するための数値解析手法であり、一般的には定常状態を想定し、波動方程式を特定の周波数で変数分離したヘルムホルツ方程式を基礎方程式とする。そして、ヘルムホルツ方程式を境界上の境界積分方程式の形態に変換した上で、要素に離散化することで音場を解析する。また、境界要素法は、有限要素法と違い、系に単位力が作用したときの系の応答を示す基本解を必要とするため、基本解が知られていない問題には適用できない。
上記の波動音響解析に使用する波動方程式を数７に、ヘルムホルツ方程式を数８に示した。さらに、このヘルムホルツ方程式に対応する境界積分方程式を数９に示した。ここで、この境界積分方程式で使用する三次元ヘルムホルツ方程式の基本解Ｇは数１０に示す形態となる。図５に、境界積分方程式の各変数の幾何学的関係を示した。
【００４９】
【数７】
【００５０】
【数８】
【００５１】
境界積分方程式
【数９】
Ｖ：法線方向粒子速度、ρ：空気密度、Ｃ：立体角
Ｐa:図５のａにおける音圧、Ｐb:図５のｂにおける音圧、Ｇ（ｒa,b）：図５のａｂ間の距離ｒa,bを変数とする基本解、Ｖb：図５における法線方向粒子速度
【００５２】
【数１０】
【００５３】
そして、数９、数１０に使用して境界上で離散化した境界積分方程式は数１１となる。
【００５４】
境界積分方程式の離散化
【数１１】
Ｃ：ｉにおける立体角を４πで割った値
領域の内部においてＣ＝１、滑らかな境界においてＣ＝１／２、領域の外側でＣ＝０となる。
【００５５】
また、図６に、離散化された境界積分方程式の各変数の幾何学的関係を示した。
数１１に示す方程式に基づいた解析に際しては、ＦＯＲＴＲＡＮで記載されたプログラムを使用した。
【００５６】
上述の離散化された境界積分方程式の適用対象とした格子モデルを、図７に示した。
解析に於ける境界の要素分割は、境界要素法では、差分法や有限要素法のように、空間を要素分割する必要はなく、境界のみを要素分割でよい。それら要素の長さは、波動性を適切に計算するために、波長の１／８以下に分割するのが妥当である。そのため、波長の短い高周波では、細かく要素分割する必要があるため、計算負荷が増大する。本願のように、低周波数領域を対象とする解析では、要素の長さは比較的長い合理的な長さとできる。実際の計算では、解析対象とした周波数領域が１００Ｈｚ以下であるため、単一の境界要素の長さは４００ｍｍとした。
【００５７】
また、解析条件として、対象に合わせて境界条件を設定する必要がある。境界条件としては、剛壁を想定した完全反射面、音源を想定した振動面、吸音面を想定したインピーダンス境界がある。実際の計算では、共鳴が問題になる空間は、音を吸収しにくいコンクリートなどの硬い壁が多いため、完全反射面かわずかな吸音を与えるインピーダンス境界を用いた。
【００５８】
図８、図９、図１０には、それぞれ、上記の２６Ｈｚ、５２Ｈｚ，７８Ｈｚについて、高音圧レベル周波数の音圧レベルが高いエリア（左側の図で、閉領域として囲まれた領域）、水平断面音圧レベル分布の一例（右側の図で、音響特性の節部を濃部で示す）を示した。周波数によっては、騒音発生源の側よりも遠方の方が大きな騒音となる場合があることが判る。また、腹部及び節部は三次元的に分布していることが分かる。
【００５９】
第３ステップ
この第３ステップでは、第２ステップで求められた、全ての高音圧レベル周波数の音に関する音響状態に基づいて、低周波数領域での全ての高音圧レベル周波数を含む三次元空間の音響状態を推定する。２６Ｈｚ、５２Ｈｚ，７８Ｈｚの結果を使用して、低周波数領域（２０〜１００Ｈｚ）の騒音として合成（具体的には、単純な合算）した場合での音圧レベルが高いエリア、水平断面音圧レベル分布の一例を図１１に示した。この図から判明するように、地下室の場合、図１１右側のエリア図で、右奥側に示される側壁１３の下部に近い部位に、第２騒音検出手段を配置すべきエリアが広がっていることが分かる。
【００６０】
第４ステップ
この第４ステップでは、第３ステップで求められる三次元空間の音響状態において、定在波の腹部で、前記三次元空間の音圧レベルが平均音圧レベルより高い位置を第２騒音検出手段の配置位置とする。
これまで説明してきた図１１に示した例の場合、空調室外機の配置側（図１の左下側）で壁面近傍の最下部位置が好ましいことが判明した。
このように、騒音が大きい位置にエラーマイクを配置することでＳ／Ｎ比の高い信号が得られ、効果的な消音が可能となる。また、特に騒音を低減したい場所、エラーマイクが劣化しにくい場所なども考慮してエラーマイクの配置位置を選定することが好ましい。
【００６１】
〔第二実施形態〕
１．本願に係る消音システムが配置される配置場所の状況
図１３は、この実施形態に係る消音システムが配置される三次元閉空間の一例としての地下室１の状態を示したものである。地下室１は図面からも判明するように、概略方形の閉空間として構成されている。
図１３に示すように、地下室１の平面視短辺側の壁面上部には、騒音発生源の一例として換気扇４０が配置される。
３．２次音源及び第２騒音検出手段の位置決定方法
この実施形態に係る２次音源及び第２騒音検出手段の配置位置検定方法では、図１２に示すように、第１ステップ、第２ステップ（騒音源からの音による音響状態の導出）、２次音源位置決定・解析ステップ（２次音源の位置決定、及び、決定された２次音源からの音による音響状態の導出）、第２騒音検出手段位置決定ステップ（第２騒音検出手段の位置決定）を実行するのである。
【００６２】
第１ステップ（＃１０）
第１ステップでは、三次元開放空間において騒音発生源から発生する騒音に関して、低周波領域において音圧レベルが高い周波数である高音圧レベル周波数を、音圧レベルに従って求める。
先に説明した換気扇（騒音発生源）４０では、４１Ｈｚ（基本周波数）にピークを持っていた。この周波数が、本願にいうところの高音圧レベル周波数である。本願で消音に対象とする低周波領域は１００Ｈｚ以下の領域である。
【００６３】
第２ステップ（＃２０）
第２ステップでは、第一実施形態と同様に、第１ステップで求められる高音圧レベル周波数の音について、騒音発生源を配置した状態における三次元閉空間内での、各高音圧レベル周波数の音が三次元空間内で発生した場合の音響状態を波動音響解析で求める。解析手法は、先の形態と同様であるので、説明は省略する。
【００６４】
図１４に、上記の４１Ｈｚについての結果を示した。
左側の図面(ａ)が音圧レベルの高い領域（方形四隅にある囲まれた３次元の領域）を示し、右側の図面（ｂ）が音圧分布で、音圧の高い部分を暗い領域として濃度で示した。
これらの図から、平面視方形の空間では、方形の四隅に腹部が形成され、方形の中央側に十字形に節部が形成されていることが判る。また、腹部及び節部は三次元的に分布していることが分かる。
【００６５】
２次音源位置決定・解析ステップ（＃３０）
この２次音源位置決定・解析ステップでは、第２ステップで求められた、高音圧レベル周波数の音に関する音響状態に基づいて、スピーカ１１（２次音源）の位置を、音響状態で腹部となっている部位に決定する。図１３から判明するように、この例では、２次音源の位置を、平面視方形となる三次元閉空間において、騒音源である換気扇４０が位置される短辺側に沿った、騒音源とは反対側の角とした。この部位は、求められた音響状態において音響の腹部の一つの位置である。
【００６６】
このようにして決定した２次音源を配置した状態における三次元閉空間内での、高音圧レベル周波数の音が三次元空間内で発生した場合の音響状態を波動音響解析で求める。解析手法は、先の形態と同様であるので、説明は省略する。
図１５に、上記の４１Ｈｚについての結果を示した。
左側の図面(ａ)が音圧レベルの高い領域（方形四隅にある囲まれた３次元の領域）を示し、右側の図面（ｂ）が音圧分布で、音圧の高い部分を暗い領域として濃度で示した。
これらの図から、平面視方形の空間では、方形の四隅に腹部が形成され、方形の中央側に十字形に節部が形成されていることが判る。
【００６７】
第２騒音検出手段位置決定ステップ（＃４０）
この第２騒音検出手段位置決定ステップでは、２次音源位置決定・解析ステップで求められる三次元空間の音響状態において、定在波の腹部で、三次元空間の音圧レベルが平均音圧レベルより高い位置をエラーマイク１０（第２騒音検出手段）の配置位置とする。
これまで説明してきた図１３に示した例の場合、２次音源に対向した位置が、好適な位置となった。図１６が、消音効果を示した音圧分布であり、図１７が周波数領域における消音効果を示したものである。目的とした４１Ｈｚで消音効果が得られている。
このように、騒音が大きい位置にエラーマイクを配置することでＳ／Ｎ比の高い信号が得られ、効果的な消音が可能となる。
【００６８】
〔第三実施形態〕
１．本願に係る消音システムが配置される配置場所の状況
図１８は、この実施形態に係る消音システムが配置される三次元閉空間の一例としての地下室１の状態を示したものである。地下室１は図面からも判明するように、概略方形の閉空間として構成されている。
図１８に示すように、地下室１の平面視短辺側の壁面上部には、騒音発生源の一例として一対の換気扇４１が配置される。
３．２次音源及び第２騒音検出手段の位置決定方法
この実施形態に係る２次音源及び第２騒音検出手段の配置位置検定方法でも、図１２に示すように、第１ステップ、第２ステップ（騒音源からの音による音響状態の導出）、２次音源位置決定・解析ステップ（２次音源の位置決定、及び、決定された２次音源からの音による音響状態の導出）、第２騒音検出手段位置決定ステップ（第２騒音検出手段の位置決定）を実行する。
【００６９】
第１ステップ（＃１０）
第１ステップでは、三次元開放空間において騒音発生源から発生する騒音に関して、低周波領域において音圧レベルが高い周波数である高音圧レベル周波数を、音圧レベルに従って求める。
先に説明した換気扇（騒音発生源）４１でも、４１Ｈｚ（基本周波数）にピークを持っていた。これら周波数が、本願にいうところの高音圧レベル周波数である。本願で消音に対象とする低周波領域は１００Ｈｚ以下の領域である。
【００７０】
第２ステップ
第２ステップでは、第一実施形態と同様に、第１ステップで求められる高音圧レベル周波数の音について、騒音発生源を配置した状態における三次元閉空間内での、各高音圧レベル周波数の音が前記三次元空間内で発生した場合の音響状態を波動音響解析で求める。解析手法は、先の形態と同様であるので、説明は省略する。
【００７１】
図１９に、上記の４１Ｈｚについての結果を示した。
左側の図面(ａ)が音圧レベルの高い領域（方形四隅にある囲まれた３次元の領域）を示し、右側の図面（ｂ）が音圧分布で、音圧の高い部分を暗い領域として濃度で示した。
これらの図から、平面視方形の空間では、方形の四隅に腹部が形成され、方形の中央側に十字形に節部が形成されていることが判る。また、腹部及び節部は三次元的に分布していることが分かる。
【００７２】
２次音源位置決定・解析ステップ（＃３０）
この２次音源位置決定・解析ステップでは、第２ステップで求められた、高音圧レベル周波数の音に関する音響状態に基づいて、スピーカ１１（２次音源）の位置を、音響状態で腹部となっている部位に２次音源の位置を決定する。この例でも、２次音源の位置を、平面視方形となる三次元閉空間において、騒音源である一方の換気扇４１が位置される短辺側に沿った騒音源側の角とした。この部位は、求められた音響状態において音響の腹部の一つの位置である。
【００７３】
このようにして決定した２次音源を配置した状態における三次元閉空間内での、高音圧レベル周波数の音が三次元空間内で発生した場合の音響状態を波動音響解析で求める。解析手法は、先の形態と同様であるので、説明は省略する。
図１５と同様な結果が得られた。
【００７４】
第２騒音検出手段位置決定ステップ（＃４０）
この第２騒音検出手段位置決定ステップでは、２次音源位置決定・解析ステップで求められる三次元空間の音響状態において、定在波の腹部で、三次元空間の音圧レベルが平均音圧レベルより高い位置をエラーマイク１０（第２騒音検出手段）の配置位置とする。
これまで説明してきた図１８に示した例の場合、２次音源に対向した位置が、好適な位置となった。図２０が、消音効果を示した音圧分布であり、図２１が周波数領域における消音効果を示したものである。目的とした４１Ｈｚで消音効果が得られている。
このように、騒音が大きい位置にエラーマイクを配置することでＳ／Ｎ比の高い信号が得られ、効果的な消音が可能となる。
【産業上の利用可能性】
【００７５】
三次元閉空間に配置される騒音発生源から発生する騒音を良好に消音するのに、ＡＮＣ技術を充分に生かせる第２騒音検出手段の位置決定方法を得ることができた。
【符号の説明】
【００７６】
１ 地下室（三次元閉空間）
４ 空調室外機の原動部
９ リファレンスマイク（第１騒音検出手段）
１０ エラーマイク（第２騒音検出手段）
１１ 消音用スピーカ（２次音源）
１２ 制御情報生成装置（制御ボード）
【特許請求の範囲】
【請求項１】
騒音発生源から発生する騒音を、第１位置に配置される第１騒音検出手段で検出するとともに、前記第１位置とは異なった第２位置に配置される第２騒音検出手段により検出し、前記第１騒音検出手段により検出される第１騒音情報及び前記第２騒音検出手段により検出される第２騒音情報に基づいて２次音源から制御音を発生して少なくとも前記第２位置で消音する消音システムにおける、前記２次音源又は前記第２騒音検出手段の位置決定方法であって、
前記騒音発生源が三次元閉空間内に配置される場合に、
三次元開放空間において前記騒音発生源から発生する騒音について、低周波領域において音圧レベルが高い周波数である高音圧レベル周波数を求める第１ステップと、
前記第１ステップで求められる高音圧レベル周波数の音について、前記騒音発生源を配置した状態で、前記三次元閉空間内において、前記高音圧レベル周波数の音が前記騒音発生源から発生した場合の音響状態を波動音響解析により求める第２ステップとを実行し、
求められた前記三次元閉空間の音響状態において、定在波の腹部で、前記三次元閉空間の音圧レベルが空間内の平均音圧レベルより高い位置を前記２次音源又は第２騒音検出手段を配置する位置とする位置決定方法。
【請求項２】
前記２次音源を配置する２次音源配置腹部と、前記第２騒音検出手段を配置する第２騒音検出手段配置腹部とが異なった腹部であり、
前記三次元閉空間において、前記２次音源と前記第２騒音検出手段とが対向した位置関係とする請求項１記載の位置決定方法。
【請求項３】
前記第１ステップにおいて、前記高音圧レベル周波数を音圧レベルに従って限られた数求め、
前記第２ステップにおいて、前記第１ステップで限られた数の求められる高音圧レベル周波数の音それぞれについて、前記騒音発生源を配置した状態で、前記三次元閉空間内において、前記各高音圧レベル周波数の音が前記騒音発生源から発生した場合の音響状態を波動音響解析により求め、
前記第２ステップで求められた、限られた数の高音圧レベル周波数の音に関する前記音響状態に基づいて、前記低周波数領域での前記限られた数の高音圧レベル周波数を含む前記三次元閉空間の音響状態を推定する第３ステップを実行し、
第３ステップで求められる前記三次元閉空間の音響状態において、定在波の腹部で、前記三次元閉空間の音圧レベルが空間内の平均音圧レベルより高い位置を前記２次音源又は第２騒音検出手段を配置する位置とする請求項１又は２記載の位置決定方法。
【請求項４】
前記第２騒音検出手段の位置として、前記第３ステップで求められる前記三次元閉空間の音響状態において、定在波の腹部で、前記三次元閉空間の音圧レベルが最高音圧レベルとなる位置を前記第２位置とする請求項３記載の位置決定方法。
【請求項５】
前記第１ステップ及び前記第２ステップを実行し、
前記第２ステップで求められた音響状態に関し、定在波の腹部を前記２次音源の位置と決定するとともに、当該決定された位置に前記２次音源を配置した状態で前記三次元閉空間内において、前記各高音圧レベル周波数の音が前記騒音発生源から発生した場合の音響状態を波動音響解析により求める２次音源位置決定・解析ステップを実行し、
前記２次音源位置決定・解析ステップで得られた音響状態において、定在波の腹部を前記第２騒音検出手段の位置である前記第２位置とする請求項２記載の位置決定方法。
【請求項６】
前記第２ステップの波動音響解析において、基礎方程式をヘルムホルツ方程式とするとともに、離散化された境界積分方程式に従って、前記三次元閉空間における各部の音圧レベルを演算する請求項１〜５のいずか一項記載の位置決定方法。
【請求項７】
前記騒音発生源が回転機器である場合に、前記第１ステップにおいて、当該回転機器から発生する騒音に関して、前記音圧レベルが高い基本周波数と当該基本周波数の倍調周波数を、前記高音圧レベル周波数とする請求項１〜６の何れか一項記載の位置決定方法。
【請求項１】
騒音発生源から発生する騒音を、第１位置に配置される第１騒音検出手段で検出するとともに、前記第１位置とは異なった第２位置に配置される第２騒音検出手段により検出し、前記第１騒音検出手段により検出される第１騒音情報及び前記第２騒音検出手段により検出される第２騒音情報に基づいて２次音源から制御音を発生して少なくとも前記第２位置で消音する消音システムにおける、前記２次音源又は前記第２騒音検出手段の位置決定方法であって、
前記騒音発生源が三次元閉空間内に配置される場合に、
三次元開放空間において前記騒音発生源から発生する騒音について、低周波領域において音圧レベルが高い周波数である高音圧レベル周波数を求める第１ステップと、
前記第１ステップで求められる高音圧レベル周波数の音について、前記騒音発生源を配置した状態で、前記三次元閉空間内において、前記高音圧レベル周波数の音が前記騒音発生源から発生した場合の音響状態を波動音響解析により求める第２ステップとを実行し、
求められた前記三次元閉空間の音響状態において、定在波の腹部で、前記三次元閉空間の音圧レベルが空間内の平均音圧レベルより高い位置を前記２次音源又は第２騒音検出手段を配置する位置とする位置決定方法。
【請求項２】
前記２次音源を配置する２次音源配置腹部と、前記第２騒音検出手段を配置する第２騒音検出手段配置腹部とが異なった腹部であり、
前記三次元閉空間において、前記２次音源と前記第２騒音検出手段とが対向した位置関係とする請求項１記載の位置決定方法。
【請求項３】
前記第１ステップにおいて、前記高音圧レベル周波数を音圧レベルに従って限られた数求め、
前記第２ステップにおいて、前記第１ステップで限られた数の求められる高音圧レベル周波数の音それぞれについて、前記騒音発生源を配置した状態で、前記三次元閉空間内において、前記各高音圧レベル周波数の音が前記騒音発生源から発生した場合の音響状態を波動音響解析により求め、
前記第２ステップで求められた、限られた数の高音圧レベル周波数の音に関する前記音響状態に基づいて、前記低周波数領域での前記限られた数の高音圧レベル周波数を含む前記三次元閉空間の音響状態を推定する第３ステップを実行し、
第３ステップで求められる前記三次元閉空間の音響状態において、定在波の腹部で、前記三次元閉空間の音圧レベルが空間内の平均音圧レベルより高い位置を前記２次音源又は第２騒音検出手段を配置する位置とする請求項１又は２記載の位置決定方法。
【請求項４】
前記第２騒音検出手段の位置として、前記第３ステップで求められる前記三次元閉空間の音響状態において、定在波の腹部で、前記三次元閉空間の音圧レベルが最高音圧レベルとなる位置を前記第２位置とする請求項３記載の位置決定方法。
【請求項５】
前記第１ステップ及び前記第２ステップを実行し、
前記第２ステップで求められた音響状態に関し、定在波の腹部を前記２次音源の位置と決定するとともに、当該決定された位置に前記２次音源を配置した状態で前記三次元閉空間内において、前記各高音圧レベル周波数の音が前記騒音発生源から発生した場合の音響状態を波動音響解析により求める２次音源位置決定・解析ステップを実行し、
前記２次音源位置決定・解析ステップで得られた音響状態において、定在波の腹部を前記第２騒音検出手段の位置である前記第２位置とする請求項２記載の位置決定方法。
【請求項６】
前記第２ステップの波動音響解析において、基礎方程式をヘルムホルツ方程式とするとともに、離散化された境界積分方程式に従って、前記三次元閉空間における各部の音圧レベルを演算する請求項１〜５のいずか一項記載の位置決定方法。
【請求項７】
前記騒音発生源が回転機器である場合に、前記第１ステップにおいて、当該回転機器から発生する騒音に関して、前記音圧レベルが高い基本周波数と当該基本周波数の倍調周波数を、前記高音圧レベル周波数とする請求項１〜６の何れか一項記載の位置決定方法。
【図１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図１２】
【図１３】
【図１７】
【図１８】
【図２１】
【図８】
【図９】
【図１０】
【図１１】
【図１４】
【図１５】
【図１６】
【図１９】
【図２０】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図１２】
【図１３】
【図１７】
【図１８】
【図２１】
【図８】
【図９】
【図１０】
【図１１】
【図１４】
【図１５】
【図１６】
【図１９】
【図２０】
【公開番号】特開２０１１−１０７６７３（Ｐ２０１１−１０７６７３Ａ）
【公開日】平成２３年６月２日（２０１１．６．２）
【国際特許分類】
【出願番号】特願２０１０−９１６４８（Ｐ２０１０−９１６４８）
【出願日】平成２２年４月１２日（２０１０．４．１２）
【出願人】（００００００２８４）大阪瓦斯株式会社 (2,453)
【Ｆターム（参考）】
【公開日】平成２３年６月２日（２０１１．６．２）
【国際特許分類】
【出願日】平成２２年４月１２日（２０１０．４．１２）
【出願人】（００００００２８４）大阪瓦斯株式会社 (2,453)
【Ｆターム（参考）】
[ Back to top ]