JP4629544B2 - 騒音対策のシミュレーション方法 - Google Patents
騒音対策のシミュレーション方法 Download PDFInfo
- Publication number
- JP4629544B2 JP4629544B2 JP2005279518A JP2005279518A JP4629544B2 JP 4629544 B2 JP4629544 B2 JP 4629544B2 JP 2005279518 A JP2005279518 A JP 2005279518A JP 2005279518 A JP2005279518 A JP 2005279518A JP 4629544 B2 JP4629544 B2 JP 4629544B2
- Authority
- JP
- Japan
- Prior art keywords
- sound
- vector
- sound source
- microphone
- microphones
- 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
Images
Landscapes
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Circuit For Audible Band Transducer (AREA)
Description
本発明は、騒音源の周辺に防音壁を設ける際に、構築する防音壁の効果を予めシミュレーションする方法に関するもので、特に、騒音源が複数の場合や反射音の影響が大きい場合のシミュレーション方法に関する。
従来、高速道路や工場等から発生する騒音を低減するため、その周辺にはコンクリートブロック、あるいは、コンクリート板などの吸遮音設備が設けられている。
ところで、工場内に設置された変電設備のトランス、あるいは、室外空調機などの騒音源となる機器の位置は予め特定することできるが、実際に観測される騒音は上記機器の配置や騒音が洩れる工場の扉や窓の位置等により異なるため、防音壁を新たに構築する場合には、騒音の到来方向を実際に観測して特定するとともに、騒音源と観測点との間に防音壁を設けた場合の到達音の低減効果を予めシミュレーションしておく必要がある。
音の到来方向を推定する方法としては、複数のマイクロフォンで得られる信号の位相差から音の到来方向を推定する、いわゆる音響学的手法が考案されている。ここで、図１３（ａ）に示すように、２つのマイクロフォンＭ１，Ｍ２を所定の間隔ｄだけ離隔して配置し、これらのマイクロフォンＭ１，Ｍ２により、θs方向から平面波として到来する音波を受音する場合を考える。上記θs方向から到来した音波は、まず第１のマイクロフォンＭ１で受音され、次に、上記音波は第２のマイクロフォンＭ２で受音されるとすると、上記第２のマイクロフォンＭ２の受音信号ｘ２（ｔ）は、上記第１のマイクロフォンＭ１の受音信号ｘ１（ｔ）に対して、音波が距離ξ＝ｄ・sinθs進行するのに要した時間τs＝ξ/ｃ（ｃ；音速）だけ遅れた信号となる。したがって、上記マイクロフォンＭ１，Ｍ２間の時間遅れτs（あるいは、マイクロフォンＭ１，Ｍ２間の位相差δ21）を求めることにより、上記音波の到来方向θsを求めることができる。
しかし、マイクロフォン数が２個では音波の到来方向θsを精度よく測定することは困難であるため、実際には、図１３（ｂ）に示すように、多数のマイクロフォンＭ１〜Ｍｍを等間隔に配置したマイクロフォンアレーを構築し、基準となるマイクロフォンＭ１に対する各マイクロフォンＭｉ（ｉ＝２〜ｍ）の位相差δm1から音波の到来方向θsを求めるようにしている。具体的には、各マイクロフォンＭ１〜Ｍｍの後段にそれぞれ遅延器Ｄ１〜Ｄｍを設けるとともに、上記遅延器Ｄ１〜Ｄｍの出力を加算する加算器Σを設けて遅延和アレーを構成する。これにより、上記各マイクロフォンＭ１〜Ｍｍの出力信号に、上記Ｍ１〜Ｍｍの幾何学的配置から求めた時間差をそれぞれ遅延して与えると、仮定した方向θsからの音波の成分はすべて同期化されるが、仮定した方向以外の方向の成分はキャンセルされて小さくなるので、上記遅延された信号を加算処理することにより、音波の到来方向θsを求めることができる。この時間差による遅延和アレー処理を用いて、ある方向から到来する音の成分を強調して取出す方法は、一般に、ビームフォーマ法（または、ビームフォーカシング法）と呼ばれている（例えば、非特許文献１参照）。
ところで、工場内に設置された変電設備のトランス、あるいは、室外空調機などの騒音源となる機器の位置は予め特定することできるが、実際に観測される騒音は上記機器の配置や騒音が洩れる工場の扉や窓の位置等により異なるため、防音壁を新たに構築する場合には、騒音の到来方向を実際に観測して特定するとともに、騒音源と観測点との間に防音壁を設けた場合の到達音の低減効果を予めシミュレーションしておく必要がある。
音の到来方向を推定する方法としては、複数のマイクロフォンで得られる信号の位相差から音の到来方向を推定する、いわゆる音響学的手法が考案されている。ここで、図１３（ａ）に示すように、２つのマイクロフォンＭ１，Ｍ２を所定の間隔ｄだけ離隔して配置し、これらのマイクロフォンＭ１，Ｍ２により、θs方向から平面波として到来する音波を受音する場合を考える。上記θs方向から到来した音波は、まず第１のマイクロフォンＭ１で受音され、次に、上記音波は第２のマイクロフォンＭ２で受音されるとすると、上記第２のマイクロフォンＭ２の受音信号ｘ２（ｔ）は、上記第１のマイクロフォンＭ１の受音信号ｘ１（ｔ）に対して、音波が距離ξ＝ｄ・sinθs進行するのに要した時間τs＝ξ/ｃ（ｃ；音速）だけ遅れた信号となる。したがって、上記マイクロフォンＭ１，Ｍ２間の時間遅れτs（あるいは、マイクロフォンＭ１，Ｍ２間の位相差δ21）を求めることにより、上記音波の到来方向θsを求めることができる。
しかし、マイクロフォン数が２個では音波の到来方向θsを精度よく測定することは困難であるため、実際には、図１３（ｂ）に示すように、多数のマイクロフォンＭ１〜Ｍｍを等間隔に配置したマイクロフォンアレーを構築し、基準となるマイクロフォンＭ１に対する各マイクロフォンＭｉ（ｉ＝２〜ｍ）の位相差δm1から音波の到来方向θsを求めるようにしている。具体的には、各マイクロフォンＭ１〜Ｍｍの後段にそれぞれ遅延器Ｄ１〜Ｄｍを設けるとともに、上記遅延器Ｄ１〜Ｄｍの出力を加算する加算器Σを設けて遅延和アレーを構成する。これにより、上記各マイクロフォンＭ１〜Ｍｍの出力信号に、上記Ｍ１〜Ｍｍの幾何学的配置から求めた時間差をそれぞれ遅延して与えると、仮定した方向θsからの音波の成分はすべて同期化されるが、仮定した方向以外の方向の成分はキャンセルされて小さくなるので、上記遅延された信号を加算処理することにより、音波の到来方向θsを求めることができる。この時間差による遅延和アレー処理を用いて、ある方向から到来する音の成分を強調して取出す方法は、一般に、ビームフォーマ法（または、ビームフォーカシング法）と呼ばれている（例えば、非特許文献１参照）。
一方、計測点に配置された複数のマイクロフォンの出力信号の位相差からではなく、複数のマイクロフォンから互いに交わる直線状に配置された複数のマイクロフォン対を構成し、対となる２つのマイクロフォンＭａ，Ｍｂ間の位相差（時間遅れＤab）と、他の対となる２つのマイクロフォンＭｃ，Ｍｄ間の位相差（時間遅れＤcd）とから音源の方向を推定する方法が提案されている（例えば、非特許文献２，３、特許文献１参照）。
すなわち、図１４に示すように、４個のマイクロフォンＭ１〜Ｍ４を、互いに直交する２直線上にそれぞれ所定の間隔で配置された２組のマイクロフォン対（Ｍ１，
Ｍ３）及びマイクロフォン対（Ｍ２，Ｍ４）を構成するように配置するとともに、第５のマイクロフォンＭ５を上記マイクロフォンＭ１〜Ｍ４の作る平面上にない位置に配置して、更に４組のマイクロフォン対（Ｍ５, Ｍ１）〜（Ｍ５, Ｍ４）を構成した場合、音の入射方向である水平角θと仰角φとは以下の式（１）及び式（２）で表わせる。
ここで、時間遅れＤijは、マイクロフォンＭiに到達する音圧信号と、このマイクロフォンＭiに対して対となるマイクロフォンＭjに到達する音圧信号との時間差であり、この対となる２つのマイクロフォンＭi及びＭjに入力される信号のクロススペクトルＰij（ｆ）を求め、更に、対象とする上記周波数ｆの位相角情報Ψ（ｒａｄ）を用いて、以下の式（３）により算出される。
これにより、上記マイクロフォンアレーを用いて音源の方向を推定する場合に比較して、少ないマイクロフォン数で音源の方向を正確に推定することができるだけでなく、屋外においても、音源の方向を精度よく推定することができる。
なお、上記式（１），（２）は、マイクロフォン間の距離を半波長とする周波数以下の平面波で成立する。また、対象となる音源位置が測定点とほぼ同一平面上にあり、仰角φを必要としない場合には、２組のマイクロフォン対（Ｍ１，Ｍ３）及び（Ｍ２，Ｍ４）のみで音源の方向である水平角θを推定することができる。
大賀寿郎，山崎芳男，金田豊；音響システムとディジタル処理，コロナ社，１９９５ 上明戸昇，野上英和，山下恭弘，財満健史，大脇雅直，杉山武，和田浩之；音情報と画像を組込んだ音源探査システムの開発，日本建築学会計画系論文集，第553号，pp17-22,2002.3 大脇雅直，財満健史，和田浩之，山下恭弘；画像に音情報を組込んだ音源探査システムの開発，電力土木、No.308，pp100-104,2003.11 特開２００３−１１１１８３号公報
すなわち、図１４に示すように、４個のマイクロフォンＭ１〜Ｍ４を、互いに直交する２直線上にそれぞれ所定の間隔で配置された２組のマイクロフォン対（Ｍ１，
Ｍ３）及びマイクロフォン対（Ｍ２，Ｍ４）を構成するように配置するとともに、第５のマイクロフォンＭ５を上記マイクロフォンＭ１〜Ｍ４の作る平面上にない位置に配置して、更に４組のマイクロフォン対（Ｍ５, Ｍ１）〜（Ｍ５, Ｍ４）を構成した場合、音の入射方向である水平角θと仰角φとは以下の式（１）及び式（２）で表わせる。
なお、上記式（１），（２）は、マイクロフォン間の距離を半波長とする周波数以下の平面波で成立する。また、対象となる音源位置が測定点とほぼ同一平面上にあり、仰角φを必要としない場合には、２組のマイクロフォン対（Ｍ１，Ｍ３）及び（Ｍ２，Ｍ４）のみで音源の方向である水平角θを推定することができる。
大賀寿郎，山崎芳男，金田豊；音響システムとディジタル処理，コロナ社，１９９５ 上明戸昇，野上英和，山下恭弘，財満健史，大脇雅直，杉山武，和田浩之；音情報と画像を組込んだ音源探査システムの開発，日本建築学会計画系論文集，第553号，pp17-22,2002.3 大脇雅直，財満健史，和田浩之，山下恭弘；画像に音情報を組込んだ音源探査システムの開発，電力土木、No.308，pp100-104,2003.11
ところで、複数の騒音源と観測点との間に防音壁を設けた場合の到達音の低減効果をシミュレーションする際には、これら複数の騒音源の音源情報を把握しておく必要があるとともに、シミュレーションの結果を、視覚的に分かり易く表示する必要がある。
また、道路や工場等から発生する騒音は、実際には、複数の方向から到来する音の合成音であることが多い。しかしながら、上記従来の音源位置の推定方法では、複数の音源の周波数が異なる場合には、受音する信号から特定の周波数の音信号を抽出して、それぞれの音源の方向を推定すればよいが、周波数が同一である場合には、互いの音が干渉し合うことによって、正しく目的音の方向を推定できない場合があった。また、音源は１個であるが、例えば、騒音源と観測点との間に硬い路面がある場合など、従来の方法で推定した観測音の音源位置は、本例の音源位置よりも低い位置になってしまうなど、路面や壁などによる反射の影響が大きい場合にも、同一周波数の複数の音源が存在する場合と同様に、音源の方向を正確に測定できなかった。
このように、騒音が複数の方向から到来する音の合成音である場合や、騒音源が１個であることが明らかな場合にも、その反射の影響が大きい場合には、防音壁を回折して観測点に到達する音の大きさが、上記騒音源が１箇所であると仮定した場合とは大きく異なるため、防音壁の遮音効果を正確にシミュレーションすることが困難であった。
また、道路や工場等から発生する騒音は、実際には、複数の方向から到来する音の合成音であることが多い。しかしながら、上記従来の音源位置の推定方法では、複数の音源の周波数が異なる場合には、受音する信号から特定の周波数の音信号を抽出して、それぞれの音源の方向を推定すればよいが、周波数が同一である場合には、互いの音が干渉し合うことによって、正しく目的音の方向を推定できない場合があった。また、音源は１個であるが、例えば、騒音源と観測点との間に硬い路面がある場合など、従来の方法で推定した観測音の音源位置は、本例の音源位置よりも低い位置になってしまうなど、路面や壁などによる反射の影響が大きい場合にも、同一周波数の複数の音源が存在する場合と同様に、音源の方向を正確に測定できなかった。
このように、騒音が複数の方向から到来する音の合成音である場合や、騒音源が１個であることが明らかな場合にも、その反射の影響が大きい場合には、防音壁を回折して観測点に到達する音の大きさが、上記騒音源が１箇所であると仮定した場合とは大きく異なるため、防音壁の遮音効果を正確にシミュレーションすることが困難であった。
本発明は、従来の問題点に鑑みてなされたもので、騒音源が複数個である場合や、反射音の影響が大きい場合でも、反射音を含む複数の騒音源の位置や方向を正確に分離して防音壁の遮音効果を正確にシミュレーションすることのできる騒音対策のシミュレーション方法を提供することを目的とする。
本発明者らは、鋭意検討を重ねた結果、複数の騒音源と観測点との間に設ける防音壁の低減効果をシミュレーションする際には、これら複数の騒音源の音源情報を把握するとともに、シミュレーションの結果を、視覚的に分かり易く表示するようにすれば、上記防音壁の防音効果を効果的にシミュレーションできることを見出し、本発明に至ったものである。
また、騒音源は特定できても、壁や床などによる反射などの影響で観測される騒音があたかも複数の方向から騒音が出ているような状態となり、騒音源の位置情報が有効でない場合には、基準となるマイクロフォンに入力する観測音の振幅をその大きさとしたベクトルは、その大きさが各音源から上記基準マイクロフォンに入力する音の振幅の大きさを表し、その位相角が上記観測音に対する位相差を表わす複素ベクトル（以下、音ベクトルという）の和で表わされることから、各マイクロフォンについても、実際に観測された各マイクロフォンに入力する観測音の上記基準マイクロフォンに対する位相差から、各音源の音ベクトルを求めることができるので、これらの音ベクトルの組から各音源毎に各マイクロフォン間の位相差を求め、この位相差から上記複数の音源の方向を推定するとともに、この推定された音源のうちの1つの音源から各マイクロフォンに到来する音の位相差を求めて、この位相差から新たな音ベクトルの組を求め、この新たな音ベクトルの組が最初に仮定した音ベクトルの組に最も近い音ベクトルの組を各音源からの音ベクトルの組とすれば、騒音源の位置情報が有効でない場合でも、複数の音源の方向を分離して特定できるので、この分離した各音源からの音について防音壁の効果をシミュレーションすれば、騒音源が複数個である場合や反射の影響が大きい場合でも、防音壁の効果を正確にシミュレーションすることができる。
また、騒音源は特定できても、壁や床などによる反射などの影響で観測される騒音があたかも複数の方向から騒音が出ているような状態となり、騒音源の位置情報が有効でない場合には、基準となるマイクロフォンに入力する観測音の振幅をその大きさとしたベクトルは、その大きさが各音源から上記基準マイクロフォンに入力する音の振幅の大きさを表し、その位相角が上記観測音に対する位相差を表わす複素ベクトル（以下、音ベクトルという）の和で表わされることから、各マイクロフォンについても、実際に観測された各マイクロフォンに入力する観測音の上記基準マイクロフォンに対する位相差から、各音源の音ベクトルを求めることができるので、これらの音ベクトルの組から各音源毎に各マイクロフォン間の位相差を求め、この位相差から上記複数の音源の方向を推定するとともに、この推定された音源のうちの1つの音源から各マイクロフォンに到来する音の位相差を求めて、この位相差から新たな音ベクトルの組を求め、この新たな音ベクトルの組が最初に仮定した音ベクトルの組に最も近い音ベクトルの組を各音源からの音ベクトルの組とすれば、騒音源の位置情報が有効でない場合でも、複数の音源の方向を分離して特定できるので、この分離した各音源からの音について防音壁の効果をシミュレーションすれば、騒音源が複数個である場合や反射の影響が大きい場合でも、防音壁の効果を正確にシミュレーションすることができる。
すなわち、本願の請求項１に記載の発明は、複数の騒音源と観測点との間に防音壁を設ける際に、複数のマイクロフォンを用いて上記複数の騒音源から観測点に到達する音を観測して、上記複数のマイクロフォンで観測した観測音から上記複数の騒音源の位置を推定した後、防音壁を設けた場合の上記複数の騒音源から観察点に到達する音の大きさを算出して上記防音壁の防音効果をシミュレーションする騒音対策のシミュレーション方法であって、
上記複数のマイクロフォンで観測した観測音を、
複数のマイクロフォンで観測した観測音の、各マイクロフォン間の位相差をそれぞれ検出する第１のステップと、
基準となるマイクロフォンに入力する観測音に対応する複素ベクトルを、大きさが観測音の大きさを表わしその位相角が０°である基準音ベクトルで表わしたとき、各音源から上記基準マイクロフォンに入力する音に対応する音ベクトルが、それぞれ、その大きさが各音源から到来する音の大きさを表わし、その位相角が上記観測音に対する位相差を表わす複素ベクトルであり、かつ、上記各音源の音ベクトルの和が上記基準音ベクトルになるように上記各音源の音ベクトルを設定する第２のステップと、
上記基準となるマイクロフォン以外の各マイクロフォンについて、上記設定された各音ベクトルと大きさが等しく、かつ、その和が当該マイクロフォンに入力する観測音の音ベクトルになるように、上記各マイクロフォンに入力する各音源の音ベクトルをそれぞれ算出し、全てのマイクロフォンについての各音源からの音ベクトルの組を求める第３のステップと、
上記音ベクトルの組の中から１つの音源方向の音ベクトルを特定し、この特定された音ベクトルについて、各マイクロフォン間の位相差を求める第４のステップと、
上記第４のステップで求めた位相差から上記特定された音源の方向を推定する第５のステップと、
上記推定された特定音源から到来する音の各マイクロフォン間の位相差を算出する第６のステップと、
各マイクロフォンについて、上記第６のステップで求めた位相差から、上記推定された特定音源の音ベクトルを各マイクロフォン毎に算出する第７のステップと、
各マイクロフォンについて、上記第３のステップで算出された音ベクトルと、上記第７のステップで算出された音ベクトルとの差のベクトルを求める第８のステップと、
上記第２のステップで設定された各音ベクトルのそれぞれについて、上記差のベクトルを求め、例えば、上記求められた差のベクトルの絶対値の和が最小となるような各音源の音ベクトルの組を特定するなど、上記求められた差のベクトルの大きさに基づいて、最も確からしい各音源の音ベクトルの組を特定し、上記特定された音ベクトルの組の各音ベクトルと上記基準音ベクトルとの位相角の差から、上記各音源の方向をそれぞれ推定する第９のステップと、
上記推定された各音源の方向の映像を採取して、上記推定された各音源の方向と上記採取された映像とから、上記各音源の位置を推定する第１０のステップとを備え、観測音を複数の騒音源からの到達音に分離する複数音源の分離方法により分離した後、上記防音壁の防音効果をシミュレーションするようにしたことを特徴とする。
上記複数のマイクロフォンで観測した観測音を、
複数のマイクロフォンで観測した観測音の、各マイクロフォン間の位相差をそれぞれ検出する第１のステップと、
基準となるマイクロフォンに入力する観測音に対応する複素ベクトルを、大きさが観測音の大きさを表わしその位相角が０°である基準音ベクトルで表わしたとき、各音源から上記基準マイクロフォンに入力する音に対応する音ベクトルが、それぞれ、その大きさが各音源から到来する音の大きさを表わし、その位相角が上記観測音に対する位相差を表わす複素ベクトルであり、かつ、上記各音源の音ベクトルの和が上記基準音ベクトルになるように上記各音源の音ベクトルを設定する第２のステップと、
上記基準となるマイクロフォン以外の各マイクロフォンについて、上記設定された各音ベクトルと大きさが等しく、かつ、その和が当該マイクロフォンに入力する観測音の音ベクトルになるように、上記各マイクロフォンに入力する各音源の音ベクトルをそれぞれ算出し、全てのマイクロフォンについての各音源からの音ベクトルの組を求める第３のステップと、
上記音ベクトルの組の中から１つの音源方向の音ベクトルを特定し、この特定された音ベクトルについて、各マイクロフォン間の位相差を求める第４のステップと、
上記第４のステップで求めた位相差から上記特定された音源の方向を推定する第５のステップと、
上記推定された特定音源から到来する音の各マイクロフォン間の位相差を算出する第６のステップと、
各マイクロフォンについて、上記第６のステップで求めた位相差から、上記推定された特定音源の音ベクトルを各マイクロフォン毎に算出する第７のステップと、
各マイクロフォンについて、上記第３のステップで算出された音ベクトルと、上記第７のステップで算出された音ベクトルとの差のベクトルを求める第８のステップと、
上記第２のステップで設定された各音ベクトルのそれぞれについて、上記差のベクトルを求め、例えば、上記求められた差のベクトルの絶対値の和が最小となるような各音源の音ベクトルの組を特定するなど、上記求められた差のベクトルの大きさに基づいて、最も確からしい各音源の音ベクトルの組を特定し、上記特定された音ベクトルの組の各音ベクトルと上記基準音ベクトルとの位相角の差から、上記各音源の方向をそれぞれ推定する第９のステップと、
上記推定された各音源の方向の映像を採取して、上記推定された各音源の方向と上記採取された映像とから、上記各音源の位置を推定する第１０のステップとを備え、観測音を複数の騒音源からの到達音に分離する複数音源の分離方法により分離した後、上記防音壁の防音効果をシミュレーションするようにしたことを特徴とする。
請求項２に記載の発明は、請求項１に記載の騒音対策のシミュレーション方法において、上記第１のステップで行う観測音の各マイクロフォン間の位相差の検出と、第５のステップで行う音源方向の推定とを、互いに交わる２つの直線上にそれぞれ所定の間隔で配置された２組のマイクロフォン対を構成するマイクロフォン間の位相差を用いて音源の方向を推定する音源位置推定手段を備えた音源位置推定装置を用いて行うようにしたことを特徴とする。
また、請求項３に記載の発明は、請求項１に記載の騒音対策のシミュレーション方法において、上記第１のステップで行う観測音の各マイクロフォン間の位相差の検出と、第５のステップで行う音源方向の推定とを、互いに交わる２つの直線上にそれぞれ所定の間隔で配置された２組のマイクロフォン対と上記２組のマイクロフォン対の作る平面上にない第５のマイクロフォンとから成るマイクロフォン群と、上記２組のマイクロフォン対を構成するマイクロフォン間の位相差、及び、上記第５のマイクロフォンと上記２組のマイクロフォン対を構成する４個のマイクロフォンのそれぞれとで構成される４組のマイクロフォン対を構成するマイクロフォン間の位相差を用いて音源の方向を推定する音源位置推定手段を備えた音源位置推定装置とを用いて行うようにしたことを特徴とする。
また、請求項３に記載の発明は、請求項１に記載の騒音対策のシミュレーション方法において、上記第１のステップで行う観測音の各マイクロフォン間の位相差の検出と、第５のステップで行う音源方向の推定とを、互いに交わる２つの直線上にそれぞれ所定の間隔で配置された２組のマイクロフォン対と上記２組のマイクロフォン対の作る平面上にない第５のマイクロフォンとから成るマイクロフォン群と、上記２組のマイクロフォン対を構成するマイクロフォン間の位相差、及び、上記第５のマイクロフォンと上記２組のマイクロフォン対を構成する４個のマイクロフォンのそれぞれとで構成される４組のマイクロフォン対を構成するマイクロフォン間の位相差を用いて音源の方向を推定する音源位置推定手段を備えた音源位置推定装置とを用いて行うようにしたことを特徴とする。
本発明によれば、複数の騒音源と観測点との間に設ける防音壁の有無による上記各騒音源からの到達音の大きさをシミュレーションする際に、複数のマイクロフォンを用いて上記複数の騒音源から観測点に到達する音を観測するとともに、以下に示すような、複数音源の分離方法を用いて上記複数のマイクロフォンで観察した観察音を複数の騒音源からの到達音に分離し、その結果を、視覚的に分かり易く表示するようにしたので、上記防音壁の防音効果を効果的にシミュレーションすることができる。
複数音源の分離方法は、その大きさが基準となるマイクロフォンに入力する観測音の振幅に対する各音源からの音の振幅の大きさを表し、その位相角が上記観測音に対する位相差を表わす音ベクトルを各音源毎に想定し、この想定された音ベクトルの組を用いて各音源の方向を推定した後、上記推定された音源の中の特定音源から各マイクロフォンに到達する音の位相差を求め、この位相差から新たな音ベクトルを求めるとともに、この新たな音ベクトルと上記設定した上記特定音源の音ベクトルとの差のベクトルが最も小さくなるような音ベクトルを求め、この音ベクトルを含む音ベクトルの組を与える各音源方向を各音源の推定方向として音源の方向を分離した後、上記複数の音源からの到達音の大きさをそれぞれ算出するもので、これにより、複数のマイクロフォンで観察した観察音を複数の騒音源からの到達音に分離できるので、音源が複数個である場合や、反射音の影響が大きい場合でも、防音壁の効果を正確にシミュレーションすることができる。
このとき、上記請求項１に記載の第１のステップで行う観測音の各マイクロフォン間の位相差の検出と、第５のステップで行う音源方向の推定とを、互いに交わる２つの直線上にそれぞれ所定の間隔で配置された２組のマイクロフォン対を構成するマイクロフォン間の位相差を用いて音源の方向を推定する音源位置推定手段を備えた音源位置推定装置、あるいは、互いに交わる２つの直線上にそれぞれ所定の間隔で配置された２組のマイクロフォン対と上記２組のマイクロフォン対の作る平面上にない第５のマイクロフォンとから成るマイクロフォン群と、上記２組のマイクロフォン対を構成するマイクロフォン間の位相差、及び、上記第５のマイクロフォンと上記２組のマイクロフォン対を構成する４個のマイクロフォンのそれぞれとで構成される４組のマイクロフォン対を構成するマイクロフォン間の位相差を用いて音源の方向を推定する音源位置推定手段を備えた音源位置推定装置とを用いて行うようにすれば、少ないマイクロフォン数で、効率よくかつ正確に音源の方向を分離することができる。
複数音源の分離方法は、その大きさが基準となるマイクロフォンに入力する観測音の振幅に対する各音源からの音の振幅の大きさを表し、その位相角が上記観測音に対する位相差を表わす音ベクトルを各音源毎に想定し、この想定された音ベクトルの組を用いて各音源の方向を推定した後、上記推定された音源の中の特定音源から各マイクロフォンに到達する音の位相差を求め、この位相差から新たな音ベクトルを求めるとともに、この新たな音ベクトルと上記設定した上記特定音源の音ベクトルとの差のベクトルが最も小さくなるような音ベクトルを求め、この音ベクトルを含む音ベクトルの組を与える各音源方向を各音源の推定方向として音源の方向を分離した後、上記複数の音源からの到達音の大きさをそれぞれ算出するもので、これにより、複数のマイクロフォンで観察した観察音を複数の騒音源からの到達音に分離できるので、音源が複数個である場合や、反射音の影響が大きい場合でも、防音壁の効果を正確にシミュレーションすることができる。
このとき、上記請求項１に記載の第１のステップで行う観測音の各マイクロフォン間の位相差の検出と、第５のステップで行う音源方向の推定とを、互いに交わる２つの直線上にそれぞれ所定の間隔で配置された２組のマイクロフォン対を構成するマイクロフォン間の位相差を用いて音源の方向を推定する音源位置推定手段を備えた音源位置推定装置、あるいは、互いに交わる２つの直線上にそれぞれ所定の間隔で配置された２組のマイクロフォン対と上記２組のマイクロフォン対の作る平面上にない第５のマイクロフォンとから成るマイクロフォン群と、上記２組のマイクロフォン対を構成するマイクロフォン間の位相差、及び、上記第５のマイクロフォンと上記２組のマイクロフォン対を構成する４個のマイクロフォンのそれぞれとで構成される４組のマイクロフォン対を構成するマイクロフォン間の位相差を用いて音源の方向を推定する音源位置推定手段を備えた音源位置推定装置とを用いて行うようにすれば、少ないマイクロフォン数で、効率よくかつ正確に音源の方向を分離することができる。
以下、本発明の最良の形態について、図面に基づき説明する。
図１は本発明の最良の形態に用いられる騒音対策シミュレーションシステムの概要を示す図で、Ｍ１〜Ｍ５は図示しない騒音源からの雑音の音圧レベルを測定するための測定用のマイクロフォン、１１は音源位置近傍の映像を採取するためのＣＣＤカメラ（以下、カメラという）、１２はローパスフィルタを備えていて、上記マイクロフォンＭ１〜Ｍ５で採取された音響情報から所定の周波数以下の成分を取り出し増幅する増幅器、１３は上記増幅された音響情報（アナログ信号）をデジタル信号に変換するＡ／Ｄ変換器、１４は上記カメラ１１の映像情報信号（アナログ信号）をデジタル信号に変換するビデオ入出力ユニットである。また、２０は上記各マイクロフォンＭ１〜Ｍ５を所定の位置に配列するためのマイクロフォンフレーム、３０は三脚から成る支持部材３１と、この支持部材３１の上部に配設された回転台３２とから成る測定用基台で、この回転台３２により、上記マイクロフォンフレーム２０を回転でき、上記マイクロフォンＭ１〜Ｍ５を水平面内で回転させることができる。
また、４０は入力手段であるキーボード４１とマイクロフォン数やサンプリング周波数などの測定パラメータを記憶するとともに、音源の位置を推定し、推定した音源位置から伝播する騒音を減衰させる防音壁の設定条件をシミュレーションするシミュレーション部４２と画像表示手段であるディスプレイ４３とを備えた音源位置推定装置で、上記シミュレーション部４２は、上記測定パラメータを記憶するパラメータファイル４２ｍを備えたデータ記憶手段４２ａと、上記Ａ／Ｄ変換されたマイクロフォンＭ１〜Ｍ５からの音響情報を用いて騒音源の方向を推定する音源位置推定手段４２ｂと、上記カメラ１１からの映像に上記推定された音源位置を示す画像を付加した画像を生成して上記ディスプレイ４３に送る画像合成手段４２ｃとを備えるとともに、上記騒音源が複数である可能性のある場合、あるいは、上記騒音の反射音の影響が強く、あたかも２音源から音が到来していると考えられる場合に、上記複数の音源の方向、または、直接音の方向と反射音の方向とを分離するための音源分離手段４２ｄと、防音壁の位置や幅や高さなどの設定条件を設定する防音壁条件設定手段４２ｅと、上記防音壁条件設定手段４２ｅで設定された防音壁の有無による上記音源分離手段４２ｄで分離された複数の音源からの音からの到達音の大きさをシミュレーションする遮音効果シミュレーション手段４２ｆと備えている。これにより、上記騒音源が複数である場合や反射音の影響が強い場合などには、上記騒音源の方向を分離して推定するとともに、上記カメラ１１からの映像に、上記推定された複数の音源位置を示す画像を付加した画像を生成して上記ディスプレイ４３に表示することができる。また、上記画像に加えて、設定した防音壁の位置や防音壁の遮音効果などについても上記ディスプレイ４３の同じ画面上に表示することができる。
図１は本発明の最良の形態に用いられる騒音対策シミュレーションシステムの概要を示す図で、Ｍ１〜Ｍ５は図示しない騒音源からの雑音の音圧レベルを測定するための測定用のマイクロフォン、１１は音源位置近傍の映像を採取するためのＣＣＤカメラ（以下、カメラという）、１２はローパスフィルタを備えていて、上記マイクロフォンＭ１〜Ｍ５で採取された音響情報から所定の周波数以下の成分を取り出し増幅する増幅器、１３は上記増幅された音響情報（アナログ信号）をデジタル信号に変換するＡ／Ｄ変換器、１４は上記カメラ１１の映像情報信号（アナログ信号）をデジタル信号に変換するビデオ入出力ユニットである。また、２０は上記各マイクロフォンＭ１〜Ｍ５を所定の位置に配列するためのマイクロフォンフレーム、３０は三脚から成る支持部材３１と、この支持部材３１の上部に配設された回転台３２とから成る測定用基台で、この回転台３２により、上記マイクロフォンフレーム２０を回転でき、上記マイクロフォンＭ１〜Ｍ５を水平面内で回転させることができる。
また、４０は入力手段であるキーボード４１とマイクロフォン数やサンプリング周波数などの測定パラメータを記憶するとともに、音源の位置を推定し、推定した音源位置から伝播する騒音を減衰させる防音壁の設定条件をシミュレーションするシミュレーション部４２と画像表示手段であるディスプレイ４３とを備えた音源位置推定装置で、上記シミュレーション部４２は、上記測定パラメータを記憶するパラメータファイル４２ｍを備えたデータ記憶手段４２ａと、上記Ａ／Ｄ変換されたマイクロフォンＭ１〜Ｍ５からの音響情報を用いて騒音源の方向を推定する音源位置推定手段４２ｂと、上記カメラ１１からの映像に上記推定された音源位置を示す画像を付加した画像を生成して上記ディスプレイ４３に送る画像合成手段４２ｃとを備えるとともに、上記騒音源が複数である可能性のある場合、あるいは、上記騒音の反射音の影響が強く、あたかも２音源から音が到来していると考えられる場合に、上記複数の音源の方向、または、直接音の方向と反射音の方向とを分離するための音源分離手段４２ｄと、防音壁の位置や幅や高さなどの設定条件を設定する防音壁条件設定手段４２ｅと、上記防音壁条件設定手段４２ｅで設定された防音壁の有無による上記音源分離手段４２ｄで分離された複数の音源からの音からの到達音の大きさをシミュレーションする遮音効果シミュレーション手段４２ｆと備えている。これにより、上記騒音源が複数である場合や反射音の影響が強い場合などには、上記騒音源の方向を分離して推定するとともに、上記カメラ１１からの映像に、上記推定された複数の音源位置を示す画像を付加した画像を生成して上記ディスプレイ４３に表示することができる。また、上記画像に加えて、設定した防音壁の位置や防音壁の遮音効果などについても上記ディスプレイ４３の同じ画面上に表示することができる。
次に、本発明による騒音対策のシミュレーション方法について、図２のフローチャートに基づき説明する。なお、説明を簡単にするため、騒音を発生する音源を同一周波数の音を発生する音源Ａと音源Ｂとの２つの音源とし、基準となるマイクロフォンを、図１４に示した四角錐状に配置された５個のマイクロフォンＭ１〜Ｍ５のうちの、上記四角錐の頂点に配置されたマイクロフォンＭ５とした場合について説明する。
最初に、マイクロフォンＭ１〜Ｍ５の出力から、上記各マイクロフォンＭ１〜Ｍ５入力される観測音の音源が１個であると仮定したときの上記観測音の音源方向を音源位置推定手段４２ｂにて推定する（ステップＳ１０）。
このとき実際の測定においては、音源の位置がマイクロフォンの位置から十分（例えば、マイクロフォン間隔の１０倍以上）離れているので、各マイクロフォンＭ１〜Ｍ５に到達する音を平面波とみなすことが可能である。そこで、本例では、音源位置を求める際に、音源の位置がマイクロフォンＭ１〜Ｍ５の位置から十分離れており、音は平面波として各マイクロフォンＭ１〜Ｍ５に入射すると仮定して音源位置を推定する。
平面波近似においては、マイクロフォンＭｉとマイクロフォンＭj間の時間遅れＤijと音源の位置の水平角θ及び仰角φとは、上述した式（１），（２）で表わせるので、各マイクロフォンＭ１〜Ｍ５の出力信号を周波数分析して、対象となる周波数ｆにおける各マイクロフォンＭｉ，Ｍｊ間の時間遅れＤijを算出することにより、上記水平角θ及び仰角φを求めることができる。以下、上記式（１），（２）を再掲する。
また、上記時間遅れＤijは、２つのマイクロフォン対（Ｍｉ，Ｍｊ）に入力される信号のクロススペクトルＰij（ｆ）を求め、更に、対象とする上記周波数ｆの位相角情報Ψ（ｒａｄ）を用いて、以下の式（３）を用いて算出される。
なお、上記推定された音源方向近傍の映像をカメラ１１により採取すれば、上記推定された音源方向と上記カメラ１１で採取した映像情報とから、音源の位置についても推定することができるとともに、音源位置推定装置４０のディスプレイ４３上に表示された音源方向近傍の映像に、上記推定された音源位置を表示することができる。
また、上記音源の位置は、各周波数毎に算出することができる。
但し、この音源位置は、上述したように、観測音の音源が１個であると仮定した場合の音源位置で、本例では、以下に示すように、音源分離手段４２ｄを用いて上記音源位置を音源Ａの位置と音源Ｂの位置とに分離する。
なお、上記観測音の音源方向あるいは音源位置と後に推定する音源Ａ，Ｂの方向あるいは位置とを比較する必要のない場合には、上記時間遅れＤijを算出するだけでよく、上記水平角θ及び仰角φの算出や上記観測音の音源位置の推定については省略してもよい。
最初に、マイクロフォンＭ１〜Ｍ５の出力から、上記各マイクロフォンＭ１〜Ｍ５入力される観測音の音源が１個であると仮定したときの上記観測音の音源方向を音源位置推定手段４２ｂにて推定する（ステップＳ１０）。
このとき実際の測定においては、音源の位置がマイクロフォンの位置から十分（例えば、マイクロフォン間隔の１０倍以上）離れているので、各マイクロフォンＭ１〜Ｍ５に到達する音を平面波とみなすことが可能である。そこで、本例では、音源位置を求める際に、音源の位置がマイクロフォンＭ１〜Ｍ５の位置から十分離れており、音は平面波として各マイクロフォンＭ１〜Ｍ５に入射すると仮定して音源位置を推定する。
平面波近似においては、マイクロフォンＭｉとマイクロフォンＭj間の時間遅れＤijと音源の位置の水平角θ及び仰角φとは、上述した式（１），（２）で表わせるので、各マイクロフォンＭ１〜Ｍ５の出力信号を周波数分析して、対象となる周波数ｆにおける各マイクロフォンＭｉ，Ｍｊ間の時間遅れＤijを算出することにより、上記水平角θ及び仰角φを求めることができる。以下、上記式（１），（２）を再掲する。
また、上記音源の位置は、各周波数毎に算出することができる。
但し、この音源位置は、上述したように、観測音の音源が１個であると仮定した場合の音源位置で、本例では、以下に示すように、音源分離手段４２ｄを用いて上記音源位置を音源Ａの位置と音源Ｂの位置とに分離する。
なお、上記観測音の音源方向あるいは音源位置と後に推定する音源Ａ，Ｂの方向あるいは位置とを比較する必要のない場合には、上記時間遅れＤijを算出するだけでよく、上記水平角θ及び仰角φの算出や上記観測音の音源位置の推定については省略してもよい。
次に、図３（ａ）に示すように、基準マイクロフォンであるマイクロフォンＭ５に入力する観測音を振幅が１で位相角が０°である複素ベクトル（以下、観測音基準ベクトルという）Ｓ５で表現するとともに、音源ＡからマイクロフォンＭ５に入力される音の複素ベクトルＡ５（以下、音ベクトルＡ５という）の振幅ＰAと位相角δA、及び、音源Ｂから上記マイクロフォンＭ５に入力する音の複素ベクトルＢ５（以下、音ベクトルＢ５という）の振幅ＰBと位相角δBとをそれぞれ設定する（ステップＳ１１）。
上記マイクロフォンＭ５に入力される観測音は、音源ＡからマイクロフォンＭ５に入力される音と音源ＢからマイクロフォンＭ５に入力される音との合成音であるので、上記音ベクトルＳ５は上記音ベクトルＡ５と上記音ベクトルＢ５とのベクトル和となる。したがって、上記音ベクトルＡ５の振幅ＰAと位相角δAとを設定すれば、音源ＢからマイクロフォンＭ５に入力される音の音ベクトルＢ５の振幅ＰBと位相角δBとは、上記振幅ＰAと上記位相角δAとを用いて一義的に求めることができる。（ステップＳ１２）。
本例では、図３（ｂ）に示すように、上記ステップＳ１１において設定される音ベクトルＡ５の振幅ＰAを０．４〜４．０の範囲で、ΔＰA＝０．１ずつ変化させるとともに、位相角δAを０°〜３６０°の範囲で、ΔδA＝１°ずつ変化させて、上記ステップＳ１１から以下に記載するステップＳ１８までを繰り返し、上記振幅ＰAと上記位相角δAとの全ての組について音源Ａ及び音源Ｂの方向を推定する。そして、この推定された音源Ａからの音と音源Ｂからの音との合成音が上記観察音に最も近い振幅ＰAと位相角δAとを有する音ベクトルＡ５を求めて上記音源Ａの方向と音源Ｂの方向とを特定する。
上記マイクロフォンＭ５に入力される観測音は、音源ＡからマイクロフォンＭ５に入力される音と音源ＢからマイクロフォンＭ５に入力される音との合成音であるので、上記音ベクトルＳ５は上記音ベクトルＡ５と上記音ベクトルＢ５とのベクトル和となる。したがって、上記音ベクトルＡ５の振幅ＰAと位相角δAとを設定すれば、音源ＢからマイクロフォンＭ５に入力される音の音ベクトルＢ５の振幅ＰBと位相角δBとは、上記振幅ＰAと上記位相角δAとを用いて一義的に求めることができる。（ステップＳ１２）。
本例では、図３（ｂ）に示すように、上記ステップＳ１１において設定される音ベクトルＡ５の振幅ＰAを０．４〜４．０の範囲で、ΔＰA＝０．１ずつ変化させるとともに、位相角δAを０°〜３６０°の範囲で、ΔδA＝１°ずつ変化させて、上記ステップＳ１１から以下に記載するステップＳ１８までを繰り返し、上記振幅ＰAと上記位相角δAとの全ての組について音源Ａ及び音源Ｂの方向を推定する。そして、この推定された音源Ａからの音と音源Ｂからの音との合成音が上記観察音に最も近い振幅ＰAと位相角δAとを有する音ベクトルＡ５を求めて上記音源Ａの方向と音源Ｂの方向とを特定する。
次に、図４を参照して、各マイクロフォンＭ１〜Ｍ４に入力される観察音の音ベクトルＳ１〜Ｓ４から、上記音源Ａ及び音源Ｂから上記マイクロフォンＭ１〜Ｍ４に入力される音の音ベクトルＡ１〜Ａ４及び音ベクトルＢ１〜Ｂ４を求める方法について説明する。
図４（ａ）に示すように、マイクロフォンＭ１の観測音ベクトルＳ１の振幅は上記観測音基準ベクトルＳ５の振幅と等しく１であり、観測音基準ベクトルＳ５と観測音ベクトルＳ１との角度差Δ51は、上記ステップＳ１０で実際に検出された時間遅れＤ51に等しい。また、上記観測音ベクトルＳ１は上記音源ＡからマイクロフォンＭ１に入力される音の音ベクトルＡ１と、上記音源ＢからマイクロフォンＭ１に入力される音の音ベクトルＢ１との和となっており、かつ、上記音ベクトルＡ１，Ｂ１の振幅ＰA1，ＰB1は、それぞれ上記音ベクトルＡ５，Ｂ５の振幅ＰA，ＰBと等しい。したがって、そのベクトル和が上記観測音ベクトルＳ１となる音ベクトルＡ１，Ｂ１としては、図４（ａ）に示すように、［Ａ11，Ｂ11］，［Ａ12，Ｂ12］の２通りが考えられる。
同様に、マイクロフォンＭ２〜Ｍ４に入力される音源Ａ及び音源Ｂからの音の音ベクトルＡ２〜Ａ４，Ｂ２〜Ｂ４についても、図４（ｂ）〜（ｄ）に示すように、それぞれ２組の音ベクトル［Ａi1，Ｂi1］，［Ａi2，Ｂi2］のいずれかになる（ｉ＝２〜４）。したがって、基準となるマイクロフォンＭ５の音ベクトルＡ５，Ｂ５に対してそれぞれ時間遅れがＤ5iである観測音が入力されるマイクロフォンＭｉに入力される音源Ａ及び音源Ｂからの音の音ベクトルの組［Ａik，Ｂik］の数は、１つの音ベクトルＡ５を設定すると、ｉ＝１〜４，ｋ＝１，２であるので、１６通りあることになる（ステップＳ１３）。
図４（ａ）に示すように、マイクロフォンＭ１の観測音ベクトルＳ１の振幅は上記観測音基準ベクトルＳ５の振幅と等しく１であり、観測音基準ベクトルＳ５と観測音ベクトルＳ１との角度差Δ51は、上記ステップＳ１０で実際に検出された時間遅れＤ51に等しい。また、上記観測音ベクトルＳ１は上記音源ＡからマイクロフォンＭ１に入力される音の音ベクトルＡ１と、上記音源ＢからマイクロフォンＭ１に入力される音の音ベクトルＢ１との和となっており、かつ、上記音ベクトルＡ１，Ｂ１の振幅ＰA1，ＰB1は、それぞれ上記音ベクトルＡ５，Ｂ５の振幅ＰA，ＰBと等しい。したがって、そのベクトル和が上記観測音ベクトルＳ１となる音ベクトルＡ１，Ｂ１としては、図４（ａ）に示すように、［Ａ11，Ｂ11］，［Ａ12，Ｂ12］の２通りが考えられる。
同様に、マイクロフォンＭ２〜Ｍ４に入力される音源Ａ及び音源Ｂからの音の音ベクトルＡ２〜Ａ４，Ｂ２〜Ｂ４についても、図４（ｂ）〜（ｄ）に示すように、それぞれ２組の音ベクトル［Ａi1，Ｂi1］，［Ａi2，Ｂi2］のいずれかになる（ｉ＝２〜４）。したがって、基準となるマイクロフォンＭ５の音ベクトルＡ５，Ｂ５に対してそれぞれ時間遅れがＤ5iである観測音が入力されるマイクロフォンＭｉに入力される音源Ａ及び音源Ｂからの音の音ベクトルの組［Ａik，Ｂik］の数は、１つの音ベクトルＡ５を設定すると、ｉ＝１〜４，ｋ＝１，２であるので、１６通りあることになる（ステップＳ１３）。
音ベクトルの組合わせを求めた後には、以下に示すように、音源Ａの方向と音源Ｂの方向とを上記１６通りの音ベクトルの組合わせにつきそれぞれ推定し、どの音ベクトルの組み合わせが実際の音源Ａの方向に近いかを推定する。なお、以下の計算は、仮定した１つの音ベクトルＡ５に対して想定される１６通りの音ベクトルの組［Ａik，Ｂik］について行うが、表現を単純化するため、上記Ａik，Ｂikを単にＡi，Ｂiと表わし、その位相角をδAi，δBiと表わす。
上記各マイクロフォンＭ１〜Ｍ４の音ベクトルＡi(ｉ＝１〜４)の角度δAiと基準マイクロフォンＭ５の音ベクトルＡ５の位相角δAとの角度差δA5iは、図５（ａ）〜（ｄ）に示すように、音源ＡからマイクロフォンＭｉへ入力した音の位相と基準マイクロフォンＭ５へ入力した音の位相との差である。したがって、上記角度差δA5iと、上記角度差δA5iを用いて算出した角度差δA13及び角度差δA24とをそれぞれ時間遅れＤ13,Ｄ24及び時間遅れＤ51〜Ｄ54として、上記式（１），（２）に代入すれば、音源Ａの方向（水平角θAと仰角φA）を推定することができる。また、音源Ｂについても、同様に、上記音ベクトルＢiの角度δBiからδB5i及びδB13及びδB24を求めることにより、音源Ｂの方向（水平角θBと仰角φB）を推定することができる（ステップＳ１４）。
音源Ａ，Ｂの方向の推定が完了した後には、上記推定した１つの音源Ａから平面波が到来してきたと仮定して、上記各マイクロフォンＭ１〜Ｍ５へ入力する平面波の、上記基準となるマイクロフォンＭ５に対する各マイクロフォンＭｉの時間遅れＤ5i(A)をそれぞれ算出する（ステップＳ１５）。上記時間遅れＤ5i(A)は、上記ステップＳ１２で仮定した基準マイクロフォンＭ５の音ベクトルＡ５を用いて推定した音源Ａの方向から上記マイクロフォンＭ５に入力される音の音ベクトルである新たな音ベクトルａ５と各マイクロフォンＭ１〜Ｍ４の新たな音ベクトルａｉ（ｉ＝１〜４）との位相差に相当する。また、上記音ベクトルａｉの振幅が上記音ベクトルＡ５の振幅ＰAに等しい。
ここで、上記新たな音ベクトルａ５を、上記ステップＳ１２で仮定した音ベクトルＡ５に等しいとして、上記新たな音ベクトルａ１〜ａ４と上記ステップＳ１３で求めた音ベクトルＡ１〜Ａ４とを比較する。すなわち、上記新たな音ベクトルａ１〜ａ４と上記ステップＳ１３で求めた音ベクトルＡ１〜Ａ４とがほぼ一致すれば、上記音ベクトルＡ５の設定が実際の音源Ａからの音ベクトルであると考えられる。
そこで、図６（ａ）〜（ｄ）に示すように、上記推定された音源Ａからの音の音ベクトルａｉと、上記ステップＳ１３において上記音ベクトルＡ５に基づいて求めたベクトルＡｉ（ｉ＝１〜４）との差のベクトルΔAiを求めた後、この差のベクトルΔAiの大きさの和を算出して保存する（ステップＳ１６）。上記新たな音ベクトルａ１〜ａ４と上記ステップＳ１３で求めた音ベクトルＡ１〜Ａ４との差が小さいほど上記設定された音ベクトルＡ５が実際の音源Ａからの音ベクトルに近いことから、上記差のベクトルΔAiの大きさの和を、上記設定された音ベクトルＡ５が実際の音源Ａの音ベクトルに近いかどうかの判定基準とすれば、実際の音源Ａの方向を精度よく推定することができる。
上記各マイクロフォンＭ１〜Ｍ４の音ベクトルＡi(ｉ＝１〜４)の角度δAiと基準マイクロフォンＭ５の音ベクトルＡ５の位相角δAとの角度差δA5iは、図５（ａ）〜（ｄ）に示すように、音源ＡからマイクロフォンＭｉへ入力した音の位相と基準マイクロフォンＭ５へ入力した音の位相との差である。したがって、上記角度差δA5iと、上記角度差δA5iを用いて算出した角度差δA13及び角度差δA24とをそれぞれ時間遅れＤ13,Ｄ24及び時間遅れＤ51〜Ｄ54として、上記式（１），（２）に代入すれば、音源Ａの方向（水平角θAと仰角φA）を推定することができる。また、音源Ｂについても、同様に、上記音ベクトルＢiの角度δBiからδB5i及びδB13及びδB24を求めることにより、音源Ｂの方向（水平角θBと仰角φB）を推定することができる（ステップＳ１４）。
音源Ａ，Ｂの方向の推定が完了した後には、上記推定した１つの音源Ａから平面波が到来してきたと仮定して、上記各マイクロフォンＭ１〜Ｍ５へ入力する平面波の、上記基準となるマイクロフォンＭ５に対する各マイクロフォンＭｉの時間遅れＤ5i(A)をそれぞれ算出する（ステップＳ１５）。上記時間遅れＤ5i(A)は、上記ステップＳ１２で仮定した基準マイクロフォンＭ５の音ベクトルＡ５を用いて推定した音源Ａの方向から上記マイクロフォンＭ５に入力される音の音ベクトルである新たな音ベクトルａ５と各マイクロフォンＭ１〜Ｍ４の新たな音ベクトルａｉ（ｉ＝１〜４）との位相差に相当する。また、上記音ベクトルａｉの振幅が上記音ベクトルＡ５の振幅ＰAに等しい。
ここで、上記新たな音ベクトルａ５を、上記ステップＳ１２で仮定した音ベクトルＡ５に等しいとして、上記新たな音ベクトルａ１〜ａ４と上記ステップＳ１３で求めた音ベクトルＡ１〜Ａ４とを比較する。すなわち、上記新たな音ベクトルａ１〜ａ４と上記ステップＳ１３で求めた音ベクトルＡ１〜Ａ４とがほぼ一致すれば、上記音ベクトルＡ５の設定が実際の音源Ａからの音ベクトルであると考えられる。
そこで、図６（ａ）〜（ｄ）に示すように、上記推定された音源Ａからの音の音ベクトルａｉと、上記ステップＳ１３において上記音ベクトルＡ５に基づいて求めたベクトルＡｉ（ｉ＝１〜４）との差のベクトルΔAiを求めた後、この差のベクトルΔAiの大きさの和を算出して保存する（ステップＳ１６）。上記新たな音ベクトルａ１〜ａ４と上記ステップＳ１３で求めた音ベクトルＡ１〜Ａ４との差が小さいほど上記設定された音ベクトルＡ５が実際の音源Ａからの音ベクトルに近いことから、上記差のベクトルΔAiの大きさの和を、上記設定された音ベクトルＡ５が実際の音源Ａの音ベクトルに近いかどうかの判定基準とすれば、実際の音源Ａの方向を精度よく推定することができる。
次にステップＳ１４に戻って、上記差のベクトルΔAiの大きさの和を算出して保存する操作を１つの音ベクトルＡ５が設定されたときに求められた１６通りの組み合わせ全てについて行う（ステップＳ１７）。
１つの音ベクトルＡ５についての差のベクトルΔAiの大きさの和の算出が終了したら、ステップ１１に戻って、上記音ベクトルＡ５の振幅ＰAまたは位相差δAを変化させて上記ステップＳ１２からステップＳ１７までの操作を繰り返し、上記音ベクトルＡ５の振幅ＰAを０．４〜４．０の範囲でΔＰA＝０．１ずつ変化させ、位相差δAを０°〜３６０°の範囲で、ΔδA＝１°ずつ変化させた全ての場合について、上記差のベクトルΔAiの大きさの和を算出して保存する操作を行う（ステップＳ１８）。
そして、上記差のベクトルΔAiの大きさの和を比較し、上記絶対値の和が最小値をとるような音ベクトルＡ５，及びＢ５を特定し、この特定された音ベクトルＡ５，及びＢ５から得られる音源Ａ及び音源Ｂの方向を、観測音を分離して得られた音源Ａ及び音源Ｂの方向とする（ステップＳ１９）。
すなわち、ステップＳ１４，Ｓ１５において、音ベクトルＡ５，及びＢ５が上記特定された音ベクトルＡ５，及びＢ５である場合の音源Ａの方向（水平角θAと仰角φA）及び音源Ｂの方向（水平角θBと仰角φB）が観測音を分離して得られた音源Ａ及び音源Ｂの方向となる。そして、上記推定された音源Ａの方向及び音源Ｂの方向を含む映像をカメラ１１により採取し、音源位置推定装置４０のディスプレイ４３上に表示された映像上に上記音源Ａ，Ｂの方向をそれぞれ表示することにより、上記音源Ａの位置と音源Ｂの位置を決定する（ステップＳ２０）。
１つの音ベクトルＡ５についての差のベクトルΔAiの大きさの和の算出が終了したら、ステップ１１に戻って、上記音ベクトルＡ５の振幅ＰAまたは位相差δAを変化させて上記ステップＳ１２からステップＳ１７までの操作を繰り返し、上記音ベクトルＡ５の振幅ＰAを０．４〜４．０の範囲でΔＰA＝０．１ずつ変化させ、位相差δAを０°〜３６０°の範囲で、ΔδA＝１°ずつ変化させた全ての場合について、上記差のベクトルΔAiの大きさの和を算出して保存する操作を行う（ステップＳ１８）。
そして、上記差のベクトルΔAiの大きさの和を比較し、上記絶対値の和が最小値をとるような音ベクトルＡ５，及びＢ５を特定し、この特定された音ベクトルＡ５，及びＢ５から得られる音源Ａ及び音源Ｂの方向を、観測音を分離して得られた音源Ａ及び音源Ｂの方向とする（ステップＳ１９）。
すなわち、ステップＳ１４，Ｓ１５において、音ベクトルＡ５，及びＢ５が上記特定された音ベクトルＡ５，及びＢ５である場合の音源Ａの方向（水平角θAと仰角φA）及び音源Ｂの方向（水平角θBと仰角φB）が観測音を分離して得られた音源Ａ及び音源Ｂの方向となる。そして、上記推定された音源Ａの方向及び音源Ｂの方向を含む映像をカメラ１１により採取し、音源位置推定装置４０のディスプレイ４３上に表示された映像上に上記音源Ａ，Ｂの方向をそれぞれ表示することにより、上記音源Ａの位置と音源Ｂの位置を決定する（ステップＳ２０）。
防音壁のシミュレーションを行う場合には、まず、防音壁条件設定手段４２ｅにて、上記ディスプレイ４３上に表示された音源Ａ，Ｂの映像を参考に、防音壁の位置や大きさについての条件を設定し（ステップＳ２１）、遮音効果シミュレーション手段４２ｆにて、上記設定された防音壁について、その遮音効果をシミュレーションする（ステップＳ２２）。この場合、騒音源が２つ（音源Ａ及び音源Ｂ）あるとして、上記各音源Ａ，Ｂからの音に対する防音壁の遮音効果についてシミュレーションする。そして、上記設定された防音壁につき、到達音のレベルが所定のレベル以下になっているかどうかを判定し（ステップＳ２３）、遮音効果が十分であれば、上記条件を新たに構築する遮音壁の設定条件とする。また、遮音効果が十分でなければ、ステップＳ２１にもどって、遮音壁の設定条件を変更し、再度シミュレーションを行う。これにより、音源が観測音の方向に１個だけあるとしてシミュレーションした場合と比較してより正確なシミュレーションを行うことができる。
具体的には、上記防音壁の厚さは、音源からの直接音については２５ｄＢ程度減衰するような厚さとするとともに、上記防音壁を、例えば、工場の敷地内などのような、防音壁が構築可能な位置に、上記観測点と音源とを結んだ線に対して垂直に配置すれば、上記音源からの直接音はほぼ完全に遮音される。一方、上記防音壁の上部及び左右からは、上記防音壁で回折された音が観測点に到達するので、上記回折音の大きさを、例えば、前川チャート等の一般的な回折音の算出方法によりシミュレーションして、上記防音壁の高さと幅とを決定する（例えば、山本貢平，高木興一：前川チャートの数式表示について，騒音制御Vol.12.No.4,pp40-43,1991.8 参照）。
上記回折音は音源と防音壁との位置関係により大きく異なる。例えば、図７（ａ）に示すように、音源が２個（音源Ａ，Ｂ）であり、かつ、上記音源が分離できない場合には、上記観測音の推定位置は、同図の黒丸で示すように、上記音源Ａ，Ｂのほぼ中間となるが、これに対して、音源を２個に分離した場合には、同図の白丸で示すように、音源Ａ，Ｂは防音壁の左右側にそれぞれあると特定されるので、防音壁５１の左右から観測点Ｋに到達するの回折音は音源が１個であるとした場合に比べて大きくなっている。したがって、防音壁５１の防音効果を高めるためには、音源が２個であることを考慮して、上記防音壁５１の幅を上記音源が分離できない場合よりも広げるか、上記音源Ａ及び上記音源Ｂに対する防音壁をそれぞれ別個に設けるなどすれば、有効な騒音対策をすることができる。
また音源が１個であっても、壁や舗装路面などによる反射音の影響が強い場合にも、音源を直接音の音源と反射音の音源とに分けて、シミュレーションする必要がある。例えば、図７（ｂ）に示すように、音源Ｃと観測点Ｔとの間に舗装路面５２があった場合には、観測音の音源が１個であるとして音源の位置を推定すると、上記舗装路面５２による反射の影響で、上記観測音の推定位置は、同図の黒丸で示すように、実際の音源Ｃの位置よりも低い位置になるので、防音壁５３の上部からの回折音の遮音が不十分になる。これに対して、上記のように音源を分離する処理を行うと、同図の白丸で示すように、音源Ｃは、上記音源Ｃの他に反射音の音源Ｃ２を有する２音源と見なすことができるので、防音壁５３を設計する際には、上部からの回折音を十分に遮音することができるような高さの防音壁５３を構築することができる。したがって、上記従来の方法に比較して遮音効果の高い防音壁５３を構築することができる。
また音源が１個であっても、壁や舗装路面などによる反射音の影響が強い場合にも、音源を直接音の音源と反射音の音源とに分けて、シミュレーションする必要がある。例えば、図７（ｂ）に示すように、音源Ｃと観測点Ｔとの間に舗装路面５２があった場合には、観測音の音源が１個であるとして音源の位置を推定すると、上記舗装路面５２による反射の影響で、上記観測音の推定位置は、同図の黒丸で示すように、実際の音源Ｃの位置よりも低い位置になるので、防音壁５３の上部からの回折音の遮音が不十分になる。これに対して、上記のように音源を分離する処理を行うと、同図の白丸で示すように、音源Ｃは、上記音源Ｃの他に反射音の音源Ｃ２を有する２音源と見なすことができるので、防音壁５３を設計する際には、上部からの回折音を十分に遮音することができるような高さの防音壁５３を構築することができる。したがって、上記従来の方法に比較して遮音効果の高い防音壁５３を構築することができる。
ここで、目標とする騒音レベルを得るための防音壁の高さを自動的に計算する方法について、図８のフローチャートに基づき説明する。
最初に、目標とする騒音レベルＬtを設定する（ステップＳ３０）とともに、防音壁の高さの初期値Ｈ0を設定する（ステップＳ３１）。そして、カウンターをｎ＝０に初期化して（ステップＳ３２）、予測点Ｘでの騒音レベルＬｘを計算し（ステップＳ３３）、この計算したＬｘと上記目標とする騒音レベルＬtとを比較する（ステップＳ３４）。
上記ステップＳ３４において、上記Ｌｘが目標騒音レベルＬtよりも低い場合には、上記Ｌｘが防音壁の高さが初期値Ｈ0であるときの値かどうかを判定し（ステップＳ３７）、初期値である場合には、オーバースペックの可能性があるので、ステップＳ３８に進んで、防音壁の高さをｈ（例えば、１０ｃｍ）だけ減らし、上記低くした防音壁について、予測点Ｘでの騒音レベルＬｘを計算する（ステップＳ３９）。そして、上記計算したＬｘと上記目標とする騒音レベルＬtとを比較する（ステップＳ４０）。
一方、上記Ｌｘが目標騒音レベルＬtよりも高い場合、すなわち、目標騒音レベルに達していない場合には、ステップＳ３５にてカウンターを１つ進めた後、ステップＳ３６に進んで、防音壁の高さをｈ（例えば、１０ｃｍ）だけ増やして、防音壁の高さをＨ1＝Ｈ0＋ｈに設定した後、上記ステップＳ３３に戻り、上記高くした防音壁について、予測点Ｘでの騒音レベルＬｘを計算し、この計算したＬｘと上記目標とする騒音レベルＬtとを比較する（ステップＳ３４）。この操作を、上記計算したＬｘが目標騒音レベルＬtよりも低くなるまで繰り返す。そして、上記Ｌｘが目標騒音レベルＬtよりも低くなった場合には、上記Ｌｘが目標騒音レベルＬtに達したので、ステップＳ４３に進み、その結果を表示する。
また、上記ステップＳ３９で計算したＬｘが目標騒音レベルＬtよりも高い場合には、防音壁の高さを最適値よりも低く設定してしまったので、ステップＳ４１に進んで、防音壁の高さをｈだけ増やして前の高さに戻した後、予測点Ｘでの騒音レベルＬｘを再度計算し（ステップＳ４２）てから、ステップＳ４３に進み、その結果を表示する。
一方、上記ステップＳ３９で計算したＬｘが目標騒音レベルＬtよりも低い場合には、上記ステップＳ３８に戻って防音壁の高さを更にｈ（例えば、１０ｃｍ）だけ減らした後、上記低くした防音壁について、予測点Ｘでの騒音レベルＬｘを計算し、この計算したＬｘと上記目標とする騒音レベルＬtとを比較する（ステップＳ４０）という操作を、上記Ｌｘが目標騒音レベルＬtよりも高くなるまで繰り返す。そして、上記Ｌｘが目標騒音レベルＬtよりも高くなった場合には、ステップＳ４１に進んで、防音壁の高さをｈだけ増やして前の高さに戻し、予測点Ｘでの騒音レベルＬｘを再度計算して（ステップＳ４２）から、ステップＳ４３に進み、その結果を表示する。
これにより、防音壁の高さを自動的に計算することができる。
なお、防音壁の幅についても、同様の方法により、自動的に計算できる。
最初に、目標とする騒音レベルＬtを設定する（ステップＳ３０）とともに、防音壁の高さの初期値Ｈ0を設定する（ステップＳ３１）。そして、カウンターをｎ＝０に初期化して（ステップＳ３２）、予測点Ｘでの騒音レベルＬｘを計算し（ステップＳ３３）、この計算したＬｘと上記目標とする騒音レベルＬtとを比較する（ステップＳ３４）。
上記ステップＳ３４において、上記Ｌｘが目標騒音レベルＬtよりも低い場合には、上記Ｌｘが防音壁の高さが初期値Ｈ0であるときの値かどうかを判定し（ステップＳ３７）、初期値である場合には、オーバースペックの可能性があるので、ステップＳ３８に進んで、防音壁の高さをｈ（例えば、１０ｃｍ）だけ減らし、上記低くした防音壁について、予測点Ｘでの騒音レベルＬｘを計算する（ステップＳ３９）。そして、上記計算したＬｘと上記目標とする騒音レベルＬtとを比較する（ステップＳ４０）。
一方、上記Ｌｘが目標騒音レベルＬtよりも高い場合、すなわち、目標騒音レベルに達していない場合には、ステップＳ３５にてカウンターを１つ進めた後、ステップＳ３６に進んで、防音壁の高さをｈ（例えば、１０ｃｍ）だけ増やして、防音壁の高さをＨ1＝Ｈ0＋ｈに設定した後、上記ステップＳ３３に戻り、上記高くした防音壁について、予測点Ｘでの騒音レベルＬｘを計算し、この計算したＬｘと上記目標とする騒音レベルＬtとを比較する（ステップＳ３４）。この操作を、上記計算したＬｘが目標騒音レベルＬtよりも低くなるまで繰り返す。そして、上記Ｌｘが目標騒音レベルＬtよりも低くなった場合には、上記Ｌｘが目標騒音レベルＬtに達したので、ステップＳ４３に進み、その結果を表示する。
また、上記ステップＳ３９で計算したＬｘが目標騒音レベルＬtよりも高い場合には、防音壁の高さを最適値よりも低く設定してしまったので、ステップＳ４１に進んで、防音壁の高さをｈだけ増やして前の高さに戻した後、予測点Ｘでの騒音レベルＬｘを再度計算し（ステップＳ４２）てから、ステップＳ４３に進み、その結果を表示する。
一方、上記ステップＳ３９で計算したＬｘが目標騒音レベルＬtよりも低い場合には、上記ステップＳ３８に戻って防音壁の高さを更にｈ（例えば、１０ｃｍ）だけ減らした後、上記低くした防音壁について、予測点Ｘでの騒音レベルＬｘを計算し、この計算したＬｘと上記目標とする騒音レベルＬtとを比較する（ステップＳ４０）という操作を、上記Ｌｘが目標騒音レベルＬtよりも高くなるまで繰り返す。そして、上記Ｌｘが目標騒音レベルＬtよりも高くなった場合には、ステップＳ４１に進んで、防音壁の高さをｈだけ増やして前の高さに戻し、予測点Ｘでの騒音レベルＬｘを再度計算して（ステップＳ４２）から、ステップＳ４３に進み、その結果を表示する。
これにより、防音壁の高さを自動的に計算することができる。
なお、防音壁の幅についても、同様の方法により、自動的に計算できる。
このように、本最良の形態によれば、騒音源と観測点との間に設ける防音壁の有無による上記複数の騒音源からの到達音の大きさをシミュレーションする際に、まず、各マイクロフォンの観測音を表わす音ベクトルが、その大きさが基準となるマイクロフォンに入力する観測音の振幅に対する各音源からの音の振幅の大きさを表し、その位相角が上記観測音に対する位相差を表わす、音源Ａ及び音源Ｂからの音の音ベクトルの和で表わされるとして、基準となるマイクロフォンＭ５に入力される各音源からの音ベクトルＡ５，Ｂ５を想定し、上記音ベクトルＡ５，Ｂ５と各マイクロフォンＭ１〜Ｍ４に入力される観測音の位相差とから、各マイクロフォン毎に音ベクトルＡｉ，Ｂｉを求めて音ベクトルの組（Ａｉ，Ｂｉ）を算出する操作を行って、この算出された音ベクトルの組から各音源Ａ，Ｂの方向を推定する。そして、上記推定された音源Ａの方向から音が到来するとして、上記到来音の基準となるマイクロフォンＭ５に対する各マイクロフォンＭｉの時間遅れＤ5i(A)をそれぞれ算出し、この時間遅れＤ5i(A)から新たな音ベクトルａｉの組を求め、この新たな音ベクトルａｉの組と上記設定した各音源の音ベクトルＡｉの組との差のベクトルΔAiの大きさの和を求める。この差のベクトルΔAiの大きさの和は、上記音ベクトルＡ５が実際の音源Ａの音ベクトルに近いかどうかの判定基準となるので、上記差のベクトルΔAiの大きさの和が最も小さくなるような音ベクトルＡｉの組を求め、この音ベクトルＡｉの組を与える各音源Ａ，Ｂの方向を各音源Ａ，Ｂの推定方向として音源Ａ，Ｂの方向を分離した後、上記複数の音源Ａ，Ｂからの到達音の大きさをそれぞれ算出して、防音壁の効果をシミュレーションするようにしたので、音源が複数個である場合や、反射音の影響が大きい場合でも、防音壁の効果を正確にシミュレーションすることができる。
なお、上記最良の形態では、音ベクトルＡ５の大きさと位相角とを変化させて音源Ａと音源Ｂの方向を推定するようにしたが、図９に示すように、音ベクトルＡ５の位相角と音ベクトルＢ５の位相角とをそれぞれ変化させて音源Ａと音源Ｂの方向を推定するようにしてもよい。この場合、音ベクトルＡ５，Ｂ５の大きさは、音ベクトルＡ５と音ベクトルＢ５との和が観測音基準ベクトルＳ５になることから一義的に求めることができる。
また、上記例では，音源が２個である場合について説明したが、本発明はこれに限るものではなく、音源が３個以上の場合でも、図１０に示すように、観測音基準ベクトルＳ５に対して各音ベクトルＡ５，Ｂ５，Ｃ５，‥‥を想定することにより、各音源Ａ，Ｂ，Ｃ‥‥からの音の到来の方向を推定することができる。例えば、音源が３個である場合には、音ベクトルＡ５，Ｂ５の大きさと位相角となる４つのパラメータを設定するか、音ベクトルＡ５大きさと位相角と、音ベクトルＢ５，Ｃ５の位相角となる４つのパラメータを設定するなどすれば、音ベクトルの組（Ａｉ，Ｂｉ，Ｃｉ）を求めることができる。その後は、音源が２つの場合と同様に、上記音ベクトルＡ５から音源Ａの方向を推定して新たな音ベクトルの組（ａｉ，ｂｉ，ｃｉ）を求めて、上記音ベクトルの組（Ａｉ，Ｂｉ，Ｃｉ）と比較して、最も確からしい（Ａｉ，Ｂｉ，Ｃｉ）の組、すなわち、実際の音源Ａ，Ｂ，Ｃの音ベクトルに近い（Ａｉ，Ｂｉ，Ｃｉ）の組を求めて、上記音源Ａ，Ｂ，Ｃの方向を推定するようにすればよい。
また、音源Ａ及び音源Ｂの方向特定の判定基準を、上述した差のベクトルΔAiの絶対値の和に代えて、差のベクトルΔAiの積を用いても良い。あるいは、上記差のベクトルΔAiの二乗の和の平方根を用いても良い。
また、音源が１個である場合でも、反射音の影響が強い場合には、音源Ａを直接音を発生する音源とし、音源Ｂを上記音源Ａからの反射音を発生する仮想音源（２次音源）とすれば、反射の起こる箇所の特定や反射音の影響を正確に把握することができる。
また、上記例では，音源が２個である場合について説明したが、本発明はこれに限るものではなく、音源が３個以上の場合でも、図１０に示すように、観測音基準ベクトルＳ５に対して各音ベクトルＡ５，Ｂ５，Ｃ５，‥‥を想定することにより、各音源Ａ，Ｂ，Ｃ‥‥からの音の到来の方向を推定することができる。例えば、音源が３個である場合には、音ベクトルＡ５，Ｂ５の大きさと位相角となる４つのパラメータを設定するか、音ベクトルＡ５大きさと位相角と、音ベクトルＢ５，Ｃ５の位相角となる４つのパラメータを設定するなどすれば、音ベクトルの組（Ａｉ，Ｂｉ，Ｃｉ）を求めることができる。その後は、音源が２つの場合と同様に、上記音ベクトルＡ５から音源Ａの方向を推定して新たな音ベクトルの組（ａｉ，ｂｉ，ｃｉ）を求めて、上記音ベクトルの組（Ａｉ，Ｂｉ，Ｃｉ）と比較して、最も確からしい（Ａｉ，Ｂｉ，Ｃｉ）の組、すなわち、実際の音源Ａ，Ｂ，Ｃの音ベクトルに近い（Ａｉ，Ｂｉ，Ｃｉ）の組を求めて、上記音源Ａ，Ｂ，Ｃの方向を推定するようにすればよい。
また、音源Ａ及び音源Ｂの方向特定の判定基準を、上述した差のベクトルΔAiの絶対値の和に代えて、差のベクトルΔAiの積を用いても良い。あるいは、上記差のベクトルΔAiの二乗の和の平方根を用いても良い。
また、音源が１個である場合でも、反射音の影響が強い場合には、音源Ａを直接音を発生する音源とし、音源Ｂを上記音源Ａからの反射音を発生する仮想音源（２次音源）とすれば、反射の起こる箇所の特定や反射音の影響を正確に把握することができる。
本発明による騒音対策のシミュレーションを行った結果を図１１（ａ）に示す。同図において、騒音源は中央の建物内にあり、正面の扉から騒音が建物の外側に伝播する。観測点の正面側に防音壁がない場合に、実測された騒音の大きさは７３．１ｄＢであった。上記騒音源を上述した音源の分離方法で分離したところ、騒音源は正面の扉が直接音の音源であり、手前の路面が反射面であることがわかった。そこで、音源を正面の扉と手前の路面として、観測点の正面側に防音壁を構築したと仮定として騒音の大きさをシミュレーションしたところ、騒音は６６．１ｄＢに低下することが予想された。一方、騒音源を分離しない場合には、音源位置が実際よりも地面に近いとしてシミュレーションすることになるため、防音壁の上部の反射音の回折は無視され、このため、騒音は６３．２ｄＢ（図示せず）まで低下すると予想される。そこで、図１１（ｂ）に示すように、実際に音源を分離してシミュレーションした結果に基づいて防音壁を構築して騒音を測定したところ、騒音の値は６６．９ｄＢとシミュレーション結果にほぼ近い値となった。これにより、本発明によるシミュレーションが有効であることが確認された。
また、図１２（ａ）に示すように、空調室外機が２台稼動している場合についても、上記２台の空調室外機から成る音源を分離して適切な防音壁の大きさをシミュレーションし、図１２（ｂ）に示すような、上記シミュレーションした防音壁と同じ防音壁を構築して騒音を測定した。防音壁がない場合に、実測された騒音の大きさは７１．１ｄＢであり、騒音源を分離した場合には、騒音は６８．８ｄＢまで低下することが予想される。一方、騒音源を分離しない場合には、空調機の位置が２台の空調機のほほ中間にあるとしてシミュレーションすることになるので、防音壁の左右からの反射音の回折は無視され、このため、騒音は６６．１ｄＢ（図示せず）まで低下すると予想される。実際に測定したときの騒音の値は６９．２ｄＢと、上記シミュレーションによる値（６８．８ｄＢ）とほぼ同じ近い値であり、本発明によるシミュレーションが有効であることが確認された。
また、図１２（ａ）に示すように、空調室外機が２台稼動している場合についても、上記２台の空調室外機から成る音源を分離して適切な防音壁の大きさをシミュレーションし、図１２（ｂ）に示すような、上記シミュレーションした防音壁と同じ防音壁を構築して騒音を測定した。防音壁がない場合に、実測された騒音の大きさは７１．１ｄＢであり、騒音源を分離した場合には、騒音は６８．８ｄＢまで低下することが予想される。一方、騒音源を分離しない場合には、空調機の位置が２台の空調機のほほ中間にあるとしてシミュレーションすることになるので、防音壁の左右からの反射音の回折は無視され、このため、騒音は６６．１ｄＢ（図示せず）まで低下すると予想される。実際に測定したときの騒音の値は６９．２ｄＢと、上記シミュレーションによる値（６８．８ｄＢ）とほぼ同じ近い値であり、本発明によるシミュレーションが有効であることが確認された。
以上説明したように、本発明によれば、防音壁の有無による上記複数の騒音源からの到達音の大きさをシミュレーションする際に、上記複数の騒音源の音源情報を把握するとともに、シミュレーションの結果を、視覚的に分かり易く表示するようしたので、防音壁の防音効果を効果的にシミュレーションすることができる。また、騒音源は特定できても、複数の方向から騒音が出ており、騒音源の位置情報が有効でない場合には、複数のマイクロフォンで観測した観測音のデータを用いて騒音源の位置を分離した後、上記複数の騒音源からの到達音の大きさをそれぞれ算出して上記防音壁の効果をシミュレーションすることにより、防音壁の効果を正確にシミュレーションすることができるようにしたので、効率がよくかつ適切な騒音対策を行うことができる。
Ｍ１〜Ｍ５ マイクロフォン、１１ カメラ、１２ 増幅器、
１３ Ａ／Ｄ変換器、１４ ビデオ入出力ユニット、２０ マイクロフォンフレーム、
３０ 基台、３１ 支持部材、３２ 回転台、４０ 音源位置推定装置、
４１ キーボード、４２ シミュレーション部、４２ａ データ記憶手段、
４２ｂ 音源位置推定手段、４２ｃ 画像合成手段、４２ｄ 音源分離手段、
４２ｅ 防音壁条件設定手段、４２ｆ 遮音効果シミュレーション手段、
４２ｍ パラメータファイル、４３ ディスプレイ。
１３ Ａ／Ｄ変換器、１４ ビデオ入出力ユニット、２０ マイクロフォンフレーム、
３０ 基台、３１ 支持部材、３２ 回転台、４０ 音源位置推定装置、
４１ キーボード、４２ シミュレーション部、４２ａ データ記憶手段、
４２ｂ 音源位置推定手段、４２ｃ 画像合成手段、４２ｄ 音源分離手段、
４２ｅ 防音壁条件設定手段、４２ｆ 遮音効果シミュレーション手段、
４２ｍ パラメータファイル、４３ ディスプレイ。
Claims (3)
- 複数の騒音源と観測点との間に防音壁を設ける際に、複数のマイクロフォンを用いて上記複数の騒音源から観測点に到達する音を観測して、上記複数のマイクロフォンで観測した観測音から上記複数の騒音源の位置を推定した後、防音壁を設けた場合の上記複数の騒音源から観察点に到達する音の大きさを算出して、上記防音壁の防音効果をシミュレーションする騒音対策のシミュレーション方法であって、
上記複数のマイクロフォンで観察した観測音を、
複数のマイクロフォンで観測した観測音の、各マイクロフォン間の位相差をそれぞれ検出する第１のステップと、
基準となるマイクロフォンに入力する観測音に対応する複素ベクトルを、大きさが観測音の大きさを表わしその位相角が０°である基準音ベクトルで表わしたとき、各音源から上記基準マイクロフォンに入力する音に対応する音ベクトルが、それぞれ、その大きさが各音源から到来する音の大きさを表わし、その位相角が上記観測音に対する位相差を表わす複素ベクトルであり、かつ、上記各音源の音ベクトルの和が上記基準音ベクトルになるように上記各音源の音ベクトルを設定する第２のステップと、
上記基準となるマイクロフォン以外の各マイクロフォンについて、上記設定された各音ベクトルと大きさが等しく、かつ、その和が当該マイクロフォンに入力する観測音の音ベクトルになるように、上記各マイクロフォンに入力する各音源の音ベクトルをそれぞれ算出し、全てのマイクロフォンについての各音源からの音ベクトルの組を求める第３のステップと、
上記音ベクトルの組の中から１つの音源方向の音ベクトルを特定し、この特定された音ベクトルについて、各マイクロフォン間の位相差を求める第４のステップと、
上記第４のステップで求めた位相差から上記特定された音源の方向を推定する第５のステップと、
上記推定された特定音源から到来する音の各マイクロフォン間の位相差を算出する第６のステップと、
各マイクロフォンについて、上記第６のステップで求めた位相差から、上記推定された特定音源の音ベクトルを各マイクロフォン毎に算出する第７のステップと、
各マイクロフォンについて、上記第３のステップで算出された音ベクトルと、上記第７のステップで算出された音ベクトルとの差のベクトルを求める第８のステップと、
上記第２のステップにおいて、各音源からの音ベクトルの大きさと位相角とを変更して、その変更された音ベクトルについて、上記第２のステップから上記第８のステップを繰り返して、上記第２のステップで設定された各音ベクトルのそれぞれについて、上記差のベクトルを求め、上記求められた差のベクトルの大きさに基づいて、最も確からしい各音源の音ベクトルの組を特定し、上記特定された音ベクトルの組の各音ベクトルと上記基準音ベクトルとの位相角の差から、上記各音源の方向をそれぞれ推定する第９のステップと、
上記推定された音源方向の映像を採取して、上記推定された各音源の方向と上記採取された映像とから、上記各音源の位置を推定する第１０のステップとを備え、観測音を複数の騒音源からの到達音に分離する複数音源の分離方法により分離した後、上記防音壁の防音効果をシミュレーションするようにしたことを特徴とする騒音対策のシミュレーション方法。 - 上記第１のステップで行う観測音の各マイクロフォン間の位相差の検出と、第５のステップで行う音源方向の推定とを、互いに交わる２つの直線上にそれぞれ所定の間隔で配置された２組のマイクロフォン対を構成するマイクロフォン間の位相差を用いて音源の方向を推定する音源位置推定手段を備えた音源位置推定装置を用いて行うようにしたことを特徴とする請求項１に記載の騒音対策のシミュレーション方法。
- 上記第１のステップで行う観測音の各マイクロフォン間の位相差の検出と、第５のステップで行う音源方向の推定とを、互いに交わる２つの直線上にそれぞれ所定の間隔で配置された２組のマイクロフォン対と上記２組のマイクロフォン対の作る平面上にない第５のマイクロフォンとから成るマイクロフォン群と、上記２組のマイクロフォン対を構成するマイクロフォン間の位相差、及び、上記第５のマイクロフォンと上記２組のマイクロフォン対を構成する４個のマイクロフォンのそれぞれとで構成される４組のマイクロフォン対を構成するマイクロフォン間の位相差を用いて音源の方向を推定する音源位置推定手段を備えた音源位置推定装置とを用いて行うようにしたことを特徴とする請求項１に記載の騒音対策のシミュレーション方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2005279518A JP4629544B2 (ja) | 2005-09-27 | 2005-09-27 | 騒音対策のシミュレーション方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2005279518A JP4629544B2 (ja) | 2005-09-27 | 2005-09-27 | 騒音対策のシミュレーション方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2007093251A JP2007093251A (ja) | 2007-04-12 |
| JP4629544B2 true JP4629544B2 (ja) | 2011-02-09 |
Family
ID=37979163
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2005279518A Active JP4629544B2 (ja) | 2005-09-27 | 2005-09-27 | 騒音対策のシミュレーション方法 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP4629544B2 (ja) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP5253268B2 (ja) * | 2009-03-30 | 2013-07-31 | 中部電力株式会社 | 音源・振動源探査システム |
| JP5693201B2 (ja) * | 2010-12-16 | 2015-04-01 | 中部電力株式会社 | 指定領域からの伝播音の再生方法とその装置 |
| US9800973B1 (en) | 2016-05-10 | 2017-10-24 | X Development Llc | Sound source estimation based on simulated sound sensor array responses |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS63228034A (ja) * | 1987-03-17 | 1988-09-22 | Kajima Corp | 騒音予測チヤ−ト図 |
| JPH04221987A (ja) * | 1990-12-21 | 1992-08-12 | Sekisui House Ltd | 外部騒音シミュレーション装置 |
| JPH06109528A (ja) * | 1992-09-25 | 1994-04-19 | Isuzu Motors Ltd | 音源探索方式 |
| JP2003111183A (ja) * | 2001-09-27 | 2003-04-11 | Chubu Electric Power Co Inc | 音源探査システム |
| JP2003214941A (ja) * | 2002-01-28 | 2003-07-30 | Kanagawa Prefecture | 防音壁による騒音予測方法およびその方法をコンピュータに実行させるプログラム |
-
2005
- 2005-09-27 JP JP2005279518A patent/JP4629544B2/ja active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS63228034A (ja) * | 1987-03-17 | 1988-09-22 | Kajima Corp | 騒音予測チヤ−ト図 |
| JPH04221987A (ja) * | 1990-12-21 | 1992-08-12 | Sekisui House Ltd | 外部騒音シミュレーション装置 |
| JPH06109528A (ja) * | 1992-09-25 | 1994-04-19 | Isuzu Motors Ltd | 音源探索方式 |
| JP2003111183A (ja) * | 2001-09-27 | 2003-04-11 | Chubu Electric Power Co Inc | 音源探査システム |
| JP2003214941A (ja) * | 2002-01-28 | 2003-07-30 | Kanagawa Prefecture | 防音壁による騒音予測方法およびその方法をコンピュータに実行させるプログラム |
Also Published As
| Publication number | Publication date |
|---|---|
| JP2007093251A (ja) | 2007-04-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP5089198B2 (ja) | 音源位置推定システム | |
| EP3516350B1 (en) | An acoustic camera and a method for revealing acoustic emissions from various locations and devices | |
| Vorländer | Computer simulations in room acoustics: Concepts and uncertainties | |
| JP5693201B2 (ja) | 指定領域からの伝播音の再生方法とその装置 | |
| JP5253268B2 (ja) | 音源・振動源探査システム | |
| JP4652191B2 (ja) | 複数音源の分離方法 | |
| Yang et al. | Leak location using blind system identification in water distribution pipelines | |
| JP2004085201A (ja) | 振動源探査システム | |
| JP4629544B2 (ja) | 騒音対策のシミュレーション方法 | |
| Angella et al. | Optimal deployment of cameras for video surveillance systems | |
| Barnard et al. | Evaluation of receiving room diffusivity and the effect on low frequency impact insulation class | |
| Torres et al. | Room acoustics analysis using circular arrays: An experimental study based on sound field plane-wave decomposition | |
| Kölsch et al. | Quantification of air leakage paths: a comparison of airflow and acoustic measurements | |
| Ackermann et al. | A ground truth on room acoustical analysis and perception (GRAP) | |
| JP5390447B2 (ja) | 騒音対策効果量の推計方法 | |
| JP6650245B2 (ja) | インパルス応答生成装置及びプログラム | |
| Bilbao et al. | Directional reverberation time and the image source method for rectangular parallelepipedal rooms | |
| JP5172909B2 (ja) | 反射音情報推定装置、反射音情報推定方法、プログラム | |
| Guidorzi et al. | A newly developed low-cost 3D acoustic positioning system: Description and application in a reverberation room | |
| Heydari et al. | Scalable real-time sound source localization method based on TDOA | |
| JP5247220B2 (ja) | 音再生装置及び音再生装置を用いた音対策のシミュレーション方法 | |
| Weigand et al. | Simulation of multiple sound particle diffraction based on the uncertainty relation-a revolution in noise immission prognosis; part ii: Evaluation by measurements | |
| Aspöck et al. | Room geometry acquisition and processing methods for geometrical acoustics simulation models | |
| JP5285665B2 (ja) | 反射音情報推定装置、反射音情報推定方法、プログラム | |
| Riegel et al. | Sonic boom propagation in urban canyons using a combined ray tracing/radiosity method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20080526 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20100528 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20100706 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20100819 |
|
| 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: 20101109 |
|
| 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: 20101111 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20131119 Year of fee payment: 3 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 4629544 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| S531 | Written request for registration of change of domicile |
Free format text: JAPANESE INTERMEDIATE CODE: R313531 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| S531 | Written request for registration of change of domicile |
Free format text: JAPANESE INTERMEDIATE CODE: R313531 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |