JP3692402B2

JP3692402B2 - Sound source search method and apparatus

Info

Publication number: JP3692402B2
Application number: JP2002308880A
Authority: JP
Inventors: 達哉石井; 健一郎長井; 秀司生沼; 克巳武田; 輝臣中谷
Original assignee: Japan Aerospace Exploration Agency JAXA
Current assignee: Japan Aerospace Exploration Agency JAXA
Priority date: 2002-10-23
Filing date: 2002-10-23
Publication date: 2005-09-07
Anticipated expiration: 2022-10-23
Also published as: JP2004144579A

Description

【０００１】
【発明の属する技術分野】
この発明は、風洞中に置かれた模型から発生する騒音、任意の軌道を飛行する飛翔体から発生する騒音、その他各種機械装置から発生する騒音等の音源を特定するために用いる音源探査方法及び装置に関する。
【０００２】
【従来の技術】
航空機などの飛翔体、回転翼、気流中に置かれた風洞模型、その他の各種の機械・プラント装置からは様々な周波数の騒音が様々な位置から放出されている。例えば、図４は、従来から知られている航空機の騒音源を例示した模式図であり、（ａ）は航空機の片側半分を示す模式図、（ｂ）はギア騒音源を例に取ったときの等音圧線のイメージ図である。図４（ａ）に示すように、航空機５０の騒音源には、エンジンに関しては、主として前方で発生するファン・圧縮機騒音５１、主として後方で発生するジェット騒音５２がある。主翼に関してはフラップ騒音／フラップ・ジェット干渉騒音５３と翼端渦音５４とがあり、尾翼に関しては排気干渉騒音５５と翼端渦音５６とがあり、胴体に関しては、境界層騒音５７やギヤ騒音５８がある。図４（ｂ）にギヤ騒音５８を例に取って示すように、騒音源５８ａは必ずしも点状の騒音源ではなくある広がりを持つが、騒音源５８ａからは距離ｒが離れるほど音圧が低下していく球状の等圧線５９が描かれる。
【０００３】
従来の音源探査技術では、複数のマイクロホンを規則的に並べてマイクロホンの受信信号を解析することで音の入射方向を求める方式が一般であり、使用されるマイクロホン群はマイクロホンアレイとかフェーズドアレイと呼ばれる。マイクロホンアレイ法には、二つの重要な技術的要素が存在する。一つはマイクロホンが受信したデータの解析法であり、もう一つはマイクロホンの配置法である。
【０００４】
データ解析法については、各マイクロホンの受信信号に重み係数をかけて加重平均する方式が一般的である。基本的な考え方は次のとおりである。
仮想解析面にある音源から発生し自由空間に放射される音は球状に伝播して、各マイクロホン（座標は既知とする）で受信される。この際、音のエネルギーは距離ｒの二乗に反比例し、且つ距離ｒを伝播するのに要する時間遅れτを伴う（伝播時間τは伝播距離ｒと音速ｃから求まる。τ＝ｒ／ｃ）。また、音圧は、音源とマイクロホンとの伝播距離ｒに反比例して低下する。図５に示すように、マイクロホンアレイ６０の表面６１に配設されたマイクロホン６２のうちｉ番目のマイクロホン６２ｉが受信する音Ｐｍ_i （ｔ）（ｔは時間）は、航空機５０の表面６３上の幾つかの想定した音源６４のうちｊ番目の音源６４ｊが発生した音Ｐｓ_j が音源６４ｊからそのマイクロホン６２ｉまでの距離ｒ_i,j を時間遅れτ_i,jを伴って伝播したときの、すべてのｊについての音源６４ｊから伝播した音の合成であり、次の式（１）で表される。この場合、音源６４の音圧がゼロの場合を含む。
Ｐｍ_i （ｔ）＝αΣ_j ｛Ｐｓ_j （ｔ−τ_i,j ）／ｒ_i,j ｝・・・（１）
なお、主流がある場合の音源と見かけマイクロホン６２との伝播距離ｒ’から求まる伝播時間（τ＝ｒ’／ｃ）を考慮すると、音源の特性を入力すれぱ各マイクロホン６２の受信音の特性を求めることができる。
【０００５】
音源の分析を行いたい面（音源面）を多数の仮想音源領域に分割し、伝播距離と伝播時間とを用いてマイクロホンの受信信号を積算すると、マイクロホンの受信信号から音源面上の音圧分布を逆算することができる。音圧分布は音圧等圧線で表されるので、音圧等圧線の極値に音源が存在すると想定される。受信信号を時系列信号のまま積算するのが時間領域法、周波数分析してから積算するのが周波数領域法である。
【０００６】
ｉ番目のマイクロホン６２ｉの信号をｊ番目の焦点（仮想音源位置）に戻すには、音圧を距離分だけ増幅し、且つ位相を伝播時間分戻せばよく、次の式（２）で表される。
Ｐｓ_i,j ＝ｗ_i ×ｒ_i,j ×Ｐｍ_j （ｔ＋τ_i,j ）・・・式（２）
ここで、ｗ_i は重み関数である。
更に、図６に示すように、航空機５０の表面においてｊ番目の焦点６４ｊの音圧について、すべてのマイクロホン６２ｉの受信信号を積算すると、ｊ番目の焦点に仮想した音源６４ｊの強さＰｓ_j が式（３）として求められる。
Ｐｓ_j （ｔ）＝Σｗ_i ×ｒ_i,j ×Ｐｍ_i （ｔ＋τ_i,j ）／Σｗ_i ・・式（３）
音源６４が存在しない領域では、積算される信号間の相関が弱いため、積算値が小さくなり、音源６４が存在する焦点における音圧との差（音圧差）が大きくなる。周波数領域法では、予め、異なるマイクロホン間のクロススペクトルを求めておき、これに距離による振幅と位相の補正を施す。例えば、ｋ番目とｌ番目のマイクロホン受信信号のクロススペクトルをＣ_i,j とし、ｊ番目の焦点とｋ，ｌマイクロホンとの距離をそれぞれ、ｒ_j,k ，ｒ_j,l 、音速をｃ、周波数をｆとする。とのとき、マイクロホンｋ，ｌでのクロススペクトルを焦点ｊに換算したときのクロススペクトルＣｓ_j,k,l は、以下のように位相と振幅について補正される。

【０００７】
列車や自動車のような音を発生しながら移動する移動体の音源位置を特定するために、移動体の走行方向に並べた一次元のマイクロホンアレイに、マイクロホンアレイに対して垂直面に２組のインテンシテイプローブを付加することが提案されている（例えば、特許文献１参照）。マイクロホンアレイの受信信号を前記の加重平均することで移動体走行方向での音源位置が求められる一方、インテンシテイプローブの受信信号によって音の強さのベクトルが求められるので走行方向に垂直方向の音源位置が特定できるとしている。二次元分布を解析することに優位性があるが、インテンシテイプローブとマイクロホンアレイの組み合わせは、計測システムとしては大掛かりなものとなり、機動性に欠ける。
【０００８】
マイクロホンの受信信号から得られる周波数スペクトルをもとに、マイクロホンと音源との間の音波の伝播経路を使って音源の周波数スペクトルを求め、両周波数スペクトルを比較し、音波の伝播経路を最適化することで音源位置を特定する方法が提案されている（例えば、特許文献２参照）。最適化手法には遺伝的アルゴリズムが用いられる。遺伝的アルゴリズムとは最適化手法の一つであり、遺伝子樽造と遺伝子の持つ適用度を使って遺伝子の選択、交配、突然変異を幾世代に渡って繰り返すことで、目的の環境条件に適合する遺伝子即ち最適値を得る手法である。因みに、遺伝的アルゴリズムの適用例は、空調冷凍機器の能動的騒音制御に際しての適応フィルタ係数の探査や、自動車内の振動騒音特性の経年変化に対する適応フィルタ係数の追従などに見受けられる。遺伝的アルゴリズムは最適化法の一種であり、最小二乗法に比べると、第二位以下の極値がある場合にも有効である反面、データ処理や制御などのある程度実時間性を要求される場合には、一般的に計算負荷が多く、現在の計算機技術では実用化は困難である。
【０００９】
一方、マイクロホンの配置法については、図１２に示すように、一次元配置（図１２（ａ））、一次元アレイの十字状、Ｘ字状、Ｖ字状等の交差配置（同（ｂ）〜（ｄ））、インテンシテイプローブとの組み合わせ（同（ｅ））、円・円弧・環状配置（同（ｆ））、球などが考えられてきた。実用上は二次元音源分布を計測する場合が多く、そのためにはマイクロホンの二次元配置、例えば、二次元格子状配置（同（ｇ）、二次元分散アレイ（同（ｈ））で十分である。
【００１０】
以上、幾つかの計測法とマイクロホンの配置法が発展的に提案されてきたが、音源探査用マイクロホンアレイによる計測法、付帯装置、校正法について従来技術については、幾つかの問題点がある。即ち、先ず、マイクロホンの配置は、解析結果に大きな影響を及ぼす。一次元分布のマイクロホンの配置では、音源の二次元分析できない。図１２（ｇ）に示した格子配置のような規則配置されたマイクロホンでは、図１２（ｉ）に示すように、二つのマイクロホン８０ａ，８０ｂの間の垂直二等分線８１上の音源の音に対して、マイクロホン８０ａ，８０ｂが同位相となるので、複数の焦点８２ａ，８２ｂ，８２ｃ・・・が垂直二等分線８１に存在してしまう。図１２（ｊ）に示すように、音源位置８３から格子状に音圧の尾根８４ａ，８４ｂが発生することが知られており、音源周りにノイズ領域８５が広がり、解析上、真音源９０と偽音源９１とが現れる。
【００１１】
また、二次元分散配置の場合、マイクロホン座標を最適化するには真音源音圧レベルのみを参照しており、周波数帯域毎に真音源音圧と偽音源音圧との差及び真音源座標のずれを考慮していない。その結果、偽音源音圧レベルが低くても真音源位置がぼやけたり、その逆が起こりうるため、解析精度に悪影響が現れる。また、最適化は最小二乗法等による各マイクロホンの指向性フィルタ更新が一般的であるが、極値が多数存在する場合には有効ではない。更に、遺伝的アルゴリズムは複数の極値が存在する場合にも最適解を判別する能力に優れるが、演算処理に時間を要し、振動騒音制御など実時間性を求められる用途よりも設計に用いる方が有効である。
【００１２】
これらの騒音源を持つ対象について、低予算で、短時間の間に、しかも適正な騒音対策を行うためには、煩いと感じる周波数帯域の音が何処から（騒音源位置）どの程度（騒音の大きさ）で発生しているかを三次元情報として精度良く突き止める必要がある。このような所謂、音源探査を行うことによって注目する騒音の発生位置、周波数帯域、音圧等が特定されれば、騒音原因判別のために要する模型試験や数値解析の作業負担が軽減され、騒音低減目標をどのくらいにすべきか、騒音発生源部位の改造具合、低減装置付加の可否、使用条件変更の有無、周囲防音等の対策の規模、機械装置の根本的な変更の是非等について、確度の高い情報を得ることができ、これらの対策決定に要する時間を短縮し、より的確で有効な対策を選択することができる。
【００１３】
【特許文献１】
特開平０８−０１５００３号公報（第６コラム、段落［００２３］〜第８コラム、段落［００３４］、図１〜図３）
【特許文献２】
特開平０９−０２１８６３号公報（第１１コラム、段落［００５６］〜第１２コラム、段落［００５８］）
【００１４】
【発明が解決しようとする課題】
そこで、解析面の推定音源領域において逆算して求めた音圧分布から真音源と偽音源とを区別して真音源を求める際に、真音源を効率的に且つ精度良く求めることができる音圧検出素子の最適配置を求める点で解決すべき課題がある。
【００１５】
この発明の目的は、音源から伝播された音圧を検出する複数の音圧検出素子の最適配置を効率良く求めて、真音源を偽音源と区別して精度良く探査することができる音源探査方法及び装置を提供することである。
【００１６】
【課題を解決するための手段】
上記の課題を解決するため、この発明による音源探査方法は、空間内に存在する音源から伝播した音の音圧を二次元平面又は三次元曲面の表面上に分散配置された複数の音圧検出素子で検出し、前記各音圧検出素子で検出した音圧信号に基づいて前記音源について推定された仮想音源の音源音圧分布を求め、前記音源音圧分布から定められる評価関数に基づいて前記音圧検出素子の前記表面上の最適配置を決定することから成っている。
【００１７】
また、この発明による音源探査装置は、二次元平面又は三次元曲面の表面を持ち且つ前記表面上に空間内に存在する音源から伝播した音の音圧を検出する複数の音圧検出素子が分散配置された音圧検出素子配列体、及び前記各音圧検出素子で検出した音圧信号に基づいて前記音源について推定された仮想音源の前記音源音圧分布を求め、前記音源音圧分布から定められる評価関数に基づいて前記音圧検出素子の前記表面上の最適配置を求める演算装置から成っている。
【００１８】
この発明による音源探査方法及び音源探査装置によれば、空間内に存在する音源から伝播した音の音圧は、音源探査装置の二次元平面又は三次元曲面の表面上に分散配置された複数の音圧検出素子で検出される。各音圧検出素子で検出した音圧信号に基づいて、例えば逆算演算によって、音源について空間中に推定された仮想音源の音源音圧分布が求められる。検出された音圧信号に基づいて求めた音源音圧分布からでは、偽音源が現れる等の現象に起因して真音源を一義的に決定することが一般には困難であるが、音源音圧分布から定められる評価関数に基づいて音圧検出素子の前記表面上の最適配置を演算にて求めることができ、そうして求められた最適配置に配置された音圧検出素子で音を検出して仮想音源の音源音圧分布を求めることで、音源探査装置による真音源を簡単に且つ効率よく決定することが可能になる。
【００１９】
この音源探査方法及び音源探査装置において、前記評価関数は、前記音圧分布から求められる前記音源の真音源位置と偽音源位置との音圧差についての非負関数及び前記真音源位置近傍音圧勾配についての非負関数の線形結合又は積であり、前記音圧検出素子の最適配置は、前記評価関数を最大化するときの前記音圧検出素子の配置座標として決定することができる。評価関数をこのように定義すると、音源の真音源位置と偽音源位置との音圧差が大きいほど、また真音源位置近傍音圧勾配が急であればあるほで、各非負関数の関数値が大きくなり、それらの線形結合又は積で定義される評価関数も大きくなる。従って、評価関数が最大化するときに最も真音源と偽音源との区別が付きやすくなり、そのときの音圧検出素子の配置座標を音圧検出素子の最適配置とすることができる。
【００２０】
評価関数を非負関数の線形結合又は積とした音源探査方法及び音源探査装置において、前記評価関数を最大化する計算過程には、前記音圧検出素子の前記配置座標を遺伝子座に並べた多数の個体から群を形成し、前記群内の前記個体の選択、交叉及び突然変異から成る過程を通じた変遷を幾世代に渡って繰り返す遺伝的アルゴリズムを適用し、前記評価関数を前記選択の適応度として用いることができる。個体の選択、交叉及び突然変異から成る過程を通じた変遷を幾世代に渡って繰り返す遺伝的アルゴリズムを適用し且つ評価関数を選択の適応度として用いることによって、評価関数に低い関数値を与える個体（音圧検出素子の配置）は、自然と排除されて、評価関数に高い関数値を与える良質な個体の割合が増加し、最適な音圧検出素子の配置が偽音源に惑わされることなく効率的に得ることができる。
【００２１】
この音源探査方法及び音源探査装置において、前記音圧検出素子で検出したすべての前記音圧信号を収録し、収録した前記音圧信号のうち前記音源音圧分布を求めるのに用いる前記音圧信号を、探査すべき周波数に応じて選択する際に、前記評価関数による前記最適配置の決定の手法を適用することができる。予め、多数の音圧検出素子で検出したすべての音圧信号を収録しておき、解析の際に、探査すべき特に注目する周波数に応じて、収録した前記音圧信号のうちの一部を選択して、音源音圧分布を求める。収録した音圧信号のうち音源音圧分布を求めるのに用いる音圧信号をどれを選択するかは、採用する一部の音圧検出素子をどこに配置するかという問題に置き換えることができるから、評価関数による最適配置の決定の手法を適用することができる。
【００２２】
この音源探査方法において、開始信号発生装置からの開始信号に基づいて、前記音圧検出素子で検出した前記音圧信号を離散化信号として同時収録することを開始することができる。音源探査装置においては、前記音圧検出素子で検出した前記音圧信号を離散化信号として同時に収録する同時収録装置、及び前記同時収録装置による収録を開始させる開始信号を出力する開始信号発生装置を備えることができる。例えば、通過する航空機や車両のような対象の騒音を探査する場合のように、通過に合わせて音圧検出を開始する必要がある場合がある。このような場合には、音よりも速い光や画像手段で対象が通過するのを検出する開始信号発生装置からの開始信号に基づいて、音圧検出素子で検出した音圧信号を離散化信号として同時収録する。
【００２３】
この音源探査方法及び音源探査装置において、前記音圧検出素子は、マイクロホン又は圧力センサとすることができる。音圧を電気信号に変換可能なトランスデューサとすることが好ましい。
【００２４】
【発明の実施の形態】
以下、添付した図面に基づいて、この発明による音源探査方法及び装置の実施例を説明する。図１はこの発明による音源探査方法の概要を説明する説明図であり、（ａ）は音源が存在する音源面と音圧検出素子が配置される面との位置関係を示した解析モデルの一例であり、同（ｂ）は音源解析結果を示す音圧分布、同（ｃ）はその一断面を取った音圧グラフである。
【００２５】
図１に示す説明図によれば、この発明による音源探査方法は、次の（１）〜（７）に示す手順を踏んで行なわれる。
（１）この音源探査方法及び装置の特徴は、各音圧検出素子座標を解析的に設計することにある。解析モデルとして、図１（ａ）に示すように、音源探査装置１の表面（この例では、二次元平板状平面）２の上に、音圧検出素子である複数のマイクロホン３の配置座標の初期値を設定する。また、音源探査対象として航空機や車両、或いは工場のような対象物１０の音源面１１に真音源１２の座標と音圧とを設定する。即ち、真音源１２は、空間内の既知座標に置かれる。また、後の計算における収束判定で用い入る判定値εを設定する。
【００２６】
（２）各マイクロホン３における受信信号を計算することにより、真音源１２からの放射音の計算を行なう。
真音源１２から各マイクロホン３に到達する音は、球状の放射特性を考慮すれば、振幅が真音源１２からマイクロホン３までの距離ｒに反比例しかつ当該距離ｒを音が伝播する時間τだけ時間遅れがある。よって、同じ真音源１２が放射した音は各マイクロホン３にて異なる振幅と異なる時間遅れ（位相という）をもって受信される。この真音源１２とマイクロホン３の音圧信号との関係は、図５に基づいて既に説明した通りである。この関係を逆手にとり、図６に既に示したように、すべてのマイクロホン３の受信信号を、真音源１２があると想定される音源装想定位置（焦点という）に、距離分増幅しかつ時間を進ませて集中積算することで、焦点の音を計算することができる。
【００２７】
（３）音源面１１に仮想音源を配置してその音圧を計算し、音源面１１の音圧等高線を計算することにより、音源探査を行なう。
真音源１２を含む解析面を細分割して各部分を焦点とみなして上記の計算を行い、焦点の最大音圧を求めると、図１（ｂ）に示すように、解析面上に音の等高線１４が得られる。等高線１４は、既知音源（即ち、真音源１２）の位置で最大値をもたらすはずであるが、既知音源位置とは異なり本来音源を置いていない位置にも極値が発生し、偽音源１３が現れる。
【００２８】
（４）二つのパラメータを計算する。
図１（ｃ）に示すように、一つのパラメータは、既知音源位置、即ち真音源１２と偽音源１３との音圧差ξであり、もう一つのパラメータは真音源音圧勾配φである。計測では、真音源１２の位置での音圧（最大値）と偽音源１３の位置での音圧（極大値）の音圧差ξを大きくすること即ち信号対ノイズ比を大きくすることが望ましい。一方で、真音源１２に相当する最大値近傍の音圧等高線の勾配φが急峻であることは音源位置の誤差範囲が小さいことを意味し、位置精度の点で望ましい。
【００２９】
（５）評価関数Ｈ（ξ，φ）を計算する。
音圧差ξ及び真音源音圧勾配φの両方を変数とし、これらの非負関数の線形結合和又は積として表される以下の評価関数（目的関数とも言う）が定義される。
Ｈ（ξ，φ）＝αＦ（ξ）＋βＧ（φ）又は
Ｈ（ξ，φ）＝αＦ（ξ）×Ｇ（φ）
但し、Ｆ（ｘ），Ｇ（ｘ）：正の非負関数、α，β：正数
この発明では、この評価関数を最大化するようにマイクロホン座標の更新と評価関数の計算を繰り返し、評価関数が最大値を取るときのマイクロホン座標を、音源探査用マイクロホンの配置として決定する。
【００３０】
（６）収束判定をする。収束判定は、Ｈ（ξ，φ）の変化がε以下であるか否かで判定する。この条件を満たせば収束したと判定する。
（７）収束判定がＮｏであれば、マイクロホン３の座標を変更して（２）に戻る。Ｙｅｓであれば終了し、マイクロホン３の座標を決定する。
【００３１】
上記の例において、音源探査装置１の表面２は、この例では二次元平面板表面であったが、これに限らず、球面などの三次元曲面板表面とすることもできる。また、音圧検出素子としては、マイクロホン３としたが、これに限らず、圧力センサ等の電気信号として検出信号を出力することができるトランスデューサとすることが好ましい。実用に差し支えない程度の探査精度を得るには、埋め込まれるマイクロホン３の個数については４個以上とすることが好ましい。
【００３２】
評価関数を最適化する計算過程において、遺伝的アルゴリズムを用いることができる。遺伝的アルゴリズムは、群に含まれる個体の遺伝子について三つの過程（選択、交叉、突然変異）を幾世代に渡って繰り返すことにより、適応度の高い遺伝子を持った群（即ち、最適解）を得る手法であり、巡回セールスマン問題、ナップサック問題等の解法として良く知られており、更に、その他、航空機機体形状最適化問題や製造プロセス管理、金融・投資の応用例にも適用されている。
【００３３】
遺伝的アルゴリズムにおいては、群にはＮ個の個体があるとし、各世代で群内の個体数は不変とされる。マイクロホンの配置決定においては、各個体はマイクロホン座標群の情報を持っており、具体的には、遺伝子座にはマイクロホンアレイに用いる数量分のマイクロホン座標が並べられる。こういった個体を数十個から数百個準備して群とする。群内の個体を選択、交叉、突然変異という三つの過程を通じて変遷させ、これを幾世代に渡って繰り返す。このとき、選択の適応度を、真音源音圧と偽音源音圧との音圧差ξ及び音源音圧勾配φを変数として一意に定まる上記の評価関数Ｈとする。
【００３４】
各個体はＭ個（＝マイクロホン個数）の遺伝子座を持ち、各遺伝子座にはマイクロホンの座標情報が書き込まれている。ｉ番目の遺伝子座に書き込まれる座標情報が、例えば、（ｘ，ｙ）＝（５，３）（単位ｍｍ）としたとき、整数表示の場合は、（ＸＹ）は（５３）であるが、二進数表示の場合には（０１０１００１１）となる。これをＭ個のマイクロホンについて繋げて個体を次のようにあらわすことができる。実数表示では、ｊ番目の個体（９１０・…・５３・…・１３）であり、中程の５３がｉ番目のマイクロホン情報に関する記述である。これに対して、二進数表示では、ｊ番目の個体（１００１１０１０・・…０１０１００１１・・…０１１０１）である。なお、初期値は乱数を使って、全個体の全遺伝子座を設定する。
【００３５】
上記の条件下で、各世代について次の三つの過程を繰り返す。
（１）選択
各個体はＭ個のマイクロホンの座標情報を有するので、各個体について前記の解析面の音圧等高線を計算し、評価関数（即ち、適応度）を求めることができ、個体と評価関数は１対１に対応する。個体については、適応度の高い順に並べることができる（ランキングと称される）。初期条件で、適応度の下限を定めておき、下限値を下回った適応度を持つ個体は不適格個体として廃棄する。群の個体数を保つために、廃棄された個数だけ優良な個体が複製される。複製の仕方は例えば、ランキング比率から決定してもよい。例えば、群が１０個の個体を持っており３個が廃棄された場合には、新たな複製を適合となった上位３つを複製する。選択によって適応度の低い、評価関数の悪い個体（マイクロホン配置）は自然と排除され、優良な適応度を持つ個体の割合が増えることになる。選択の過程にはエリート選抜という行為を加えることもできる。これは、予め設定した適応度の上限を超えた個体については、以下の交叉や突然変異の確率を極度に下げるかゼロとする行為である。
【００３６】
（２）交叉
予め設定する交叉確率で選ばれた異なる伺体同士が遺伝子座の一部（その位置と長さを確率的に選択できる）を交換する行為である。交叉によって、群の個体の傾向に変化を与え、複数の極値から最大値を探査することにつながる。交叉の例として、次に示す個体Ａと個体Ｂが中程に記載部分の遺伝子座を交換して、それぞれが個体Ａ’個体Ｂ’となるときは以下のようになる。
個体Ａ個体Ｂ
（０１１０…１１０１０…０１０１）（１０００…１０００１…１１１０）
→
個体Ａ’ 個体Ｂ’
（０１１０…１０００１…０１０１）（１０００…１１０１０…１１１０）
【００３７】
（３）突然変異
予め設定する突然変異確率で選ばれた個体が、遺伝子座の一部又は全てを交換する行為である。突然変異により、群集団の変化が活性化される。遺伝子座の値については乱数を使ったり、二進法であればビット変換するなど様々な手法が用いられる。
例を、以下に示す。
個体Ａ個体Ａ’
（０１１０…１１０１０…０１０１）→（１００１…００１０１…１０１０）
以上の三つの過程を数百〜数千世代繰り返すうちに、優良な個体の割合が増し、最適なマイクロホン配置が決定される。
【００３８】
次に、図２を参照して、音圧検出素子の選択的使用について説明する。図２は、この発明による音源探査方法及び装置において、音圧検出素子の選択的使用の概要を模式的に説明する説明図である。上記の例では、音源探査に使用する音圧検出用素子の数とその座標とを、実際の計測前に検出対象に合わせて最適化したが、図２に示す例では、予め分散配置又は格子配置された多数の音圧検出用素子で音圧データを収録しておき、騒音解析の際に、その解析で用いる周波数に応じてマイクロホンの数と座標とを選択する。即ち、図２（ａ）に示すように、音源探査装置２０は、平面状の表面２２上の既知の座標に、分散して配置（図示の例では格子配置）した多数のマイクロホン２３を備えている。例えば、分散配置された１００個のマイクロホン２３のすべてを用いて同時に騒音（模式的な音波２４）の計測が行われ、計測した騒音データはすべて収録される。
【００３９】
音圧検出素子の選択的使用においては、とにかく多数の点で騒音計測をしておき、騒音解析の段階で解析に用いるマイクロホン２３が選ばれる。マイクロホンの組合せとその組合せのマイクロホン座標を用いて、評価関数を最適化する解析に使用するマイクロホンを周波数毎に選択する。例えば、周波数１００Ｈｚではマイクロホン２３を１０個を選択し、実際の解析では、その選択されたマイクロホン２３が計測した騒音データが採用される。同様に、１０００Ｈｚでは５０個を、２０００Ｈｚでは１００個全てを使うという選択方法で、マイクロホン２３の選択が行なわれる。マイクロホン２３の選択方法は、上記した評価関数を用いる方法が利用され、最適配置を示したマイクロホン２３を選択すればよい。
【００４０】
具体的には、図２（ｂ）に示すように、少数のマイクロホン２３を用いて位置精度は高くはないが、周波数帯ｗ１，ｗ２に応じて騒音源２５，２６の位置が分析して求められる。図２（ｂ）で求めた位置に騒音源２５，２６が存在すると仮定して、上記した音圧差ξと音圧勾配φの関数である評価関数Ｈを用いる方法を適用して、評価関数Ｈが最大となるマイクロホン２３の組合せが選択される。図２（ｃ）に示すように、選択されたマイクロホン２３の組合せを用いて再度、精密な騒音源探査を行い、騒音源２５，２６のそれぞれについて高い位置精度で音源探査することができる。
【００４１】
本発明による遺伝的アルゴリズムを用いたマイクロホンアレイによる音源探査のシミュレーションを実行した。以下に、そのシミュレーションの条件及び結果を示す。図７は、マイクロホンアレイを設計するための数学モデルを示す図である。図７中、７０は４１個のマイクロホン７２が存在する計測面であり、原点７１を基準に座標が定められている。原点７１に１個のマイクロホン７２が配置され、原点７１を中心にそれぞれ８個のマイクロホン７２が配置される５つの同じ領域が円周方向に等間隔で存在するように設定されている。なお、５角形に限らず他の多角形でもよい。図７中、７５は仮想音源７６，７７が存在する解析面であり、計測面７０に平行とし、面間距離ｈは４ｍに設定されている。音源７６の座標は［０．０，０．０，４．０］に設定し、音源７７の座標は［−０．７，−０．７，４．０］に設定された。両音源７５，７７共に、２０００Ｈｚで音圧１の純音を発生するとする。なお、解析面７５は、実際は計測面７０と平行でなくてもよく、複数の解析面を選択することで三次元の音源探査が可能である。
【００４２】
遺伝的アルゴリズムのパラメータは、交叉確率が０．８、突然変異確率が０．０５、エリート確率が０．０５、群個体数が４０、個体当遺伝子長さが１４４である。また、適応度Ｈは、真音源音圧と偽音源音圧との差をＦ、真音源音圧勾配をＧ、ペナルティ関数（マイクロホン最小間隔、マイクロホンアレイ寸法、最大音圧値のずれ）をＰとしたとき、Ｈ＝Ｆ×Ｇ×Ｐで定める。
【００４３】
図８は、最大適応度を持つ個体（即ち、マイクロホン座標の組合せ）を異なる世代ごとに表した図である。図８（ａ）〜（ｄ）は、それぞれ、第１世代、第６世代、第２１世代、及び第５１世代におけるマイクロホン座標を示しており、適応度は、それぞれ、２３．６０、２９．３５、４２．６６、及び４５．４２である。世代が進んで進化するにつれて適応度が増加し、マイクロホンの配置が分散していくのが解る。また、図９は、最大適応度を持つ個体によって推定される音源の位置を、図８と同様の異なる世代（第１、第６、第２１、及び第５１世代）毎に、二次元等高線表示（左側）と三次元表示（右側）とで示した図である。図中の小さい点はマイクロホンを表しており、図８に示すマイクロホンの座標に対応している。世代が進むにつれて真音源音圧と周囲の雑音圧との音圧差が広がると共に、音源位置での音圧勾配が急峻になる様子を見て取ることができる。
【００４４】
図１０は、上記の最適設計で得られたマイクロホンアレイ（図８に示す第５１世代の座標を使用）を用いて音源探査をシミュレーションしたときの音圧分布の結果を示す図であり、図１０（ａ）は二次元等高線表示であって、図１０（ｂ）は（ａ）の三次元表示である。また、図１１は、既存の十字状のマイクロホンアレイを用いて音源探査をシミュレーションしたときの音圧分布の結果を示す図であり、図１１（ａ）は二次元等高線表示であって、図１１（ｂ）は（ａ）の三次元表示である。図１０及び図１１中の点は、マイクロホン位置を示している。使用する音源は、周波数２０００Ｈｚ、音圧１の点音源とし、座標は［０．５，０．５，４］にあるものとする。図１１に示す十字状のマイクロホンアレイでは、音源位置を中心にしてアレイの配置方向に音圧の峰は偽音源として現れる。これに対して、本発明による遺伝的アルゴリズムによって最適設計されたマイクロホンアレイを用いて測定・解析すると、音源位置以外では音圧の増加が抑制されるとともに、音源位置での音圧勾配が急峻になっていることが見て取れる。
【００４５】
図３は、騒音探査の対象物を飛行中又は走行中の航空機としたとき、本発明による音源探査装置を用いた騒音源探査の態様の概要を示す概略図である。音源探査装置３０は、図１（ａ）に示した音源探査装置１と同様に、平板状の表面２上に多数のマイクロホン３を備えている。表面２を与える平板４と、多数のマイクロホン３とは、この音源探査装置３０において、音圧検出素子配列体５を定めている。このシステムでは、騒音源探査の対象物は飛行中又は走行中の航空機４０であり、エンジン４１のジェット騒音を真音源４２と想定している。音源探査装置３０は、航空機４０が探査可能な領域を通過した瞬間をレーザー光線の反射として光学的に捕らえて捕捉信号を出力することができる開始信号発生装置３１を備えている。
【００４６】
音源探査装置３０においては、開始信号発生装置３１からの捕捉信号に応答して、電子計算機３４が行なう制御により、すべてのマイクロホン３からの音圧信号は増幅器３２で増幅され、増幅された音圧信号は同時収録装置３３によって同時に離散化信号として収録される。収録されたデータは、電子計算機３４によって解析されてデータ記憶装置３５に記憶され、またディスプレイのような結果表示装置３６に表示される。同時収録装置３３に入力されたデータは、マイクロホン３からの音圧信号と連動させて解析することで、航空機４０の瞬時位置ごとに騒音源４２を特定することができる。開始信号発生装置３１は、同時収録を開始させる始動装置の役割を果たすが、航空機４０の通過については、レーザー光を用いた捕捉以外にも、光学的にはＣＣＤ画像により、音響学的には超音波反射により、更にその他、航空機４０の通過に伴う圧力変化を電圧に変換することによって検出することができる。電子計算機３４は、各音圧検出素子３で検出した音圧信号に基づいて音源４２について推定された仮想音源の音源音圧分布を求め、音源音圧分布から定められる評価関数に基づいて音圧検出素子３の表面２上の最適配置を求める音源探査装置３０についての演算装置として機能している。
【００４７】
この発明による音源探査を、航空機における騒音源の探査を例に採って説明したが、本発明の適用は、この例に限らず、風洞実験における模型或いは実機の騒音源探査、自動車、トラック等の車両、鉄道車両等の移動音源を持つ物体の騒音源探査、工場、プラント、発電設備、建築物等の大型設備の騒音源探査、軸流或いは遠心式圧縮機、送風機、タービン等の各種回転機械についての騒音源探査、或いはＯＡ機器や家電製品の騒音異常診断にも適用できることは明らかである。
【００４８】
【発明の効果】
この発明による音源探査方法及び装置は、上記のように、空間内に存在する音源から伝播した音の音圧を二次元平面又は三次元曲面の表面上に分散配置された複数の音圧検出素子で検出し、各音圧検出素子で検出した音圧信号に基づいて音源について推定された仮想音源の音源音圧分布を求め、音源音圧分布から定められる評価関数に基づいて音圧検出素子の前記表面上の最適配置を決定しているので、そうして求められた最適配置に配置された音圧検出素子で音を検出して仮想音源の音源音圧分布を求めることで、真音源と偽音源とを区別して、音源探査装置による真音源を簡単に効率よく且つ精度良く決定することができる。このように、騒音源の周波数帯域、騒音源位置、及び騒音の大きさについて三次元情報として確度の高い情報を得ることが可能になるので、騒音原因判別のための作業負担が軽減され、騒音低減目標、騒音発生源部位の改造程度、低減装置付加の可否、使用条件変更の有無、周囲防音等の対策の規模、機械装置の根本的な変更の是非等のより的確で有効な騒音対策を短時間で選択することができる。
【図面の簡単な説明】
【図１】この発明による音源探査方法の概要を説明する説明図である。
【図２】この発明による音源探査方法及び装置において、音圧検出素子の選択的使用の概要を模式的に説明する説明図である。
【図３】騒音探査の対象物を飛行中又は走行中の航空機としたとき、本発明による音源探査装置を用いた騒音源探査の態様の概要を示す概略図である。
【図４】航空機の騒音源を示す模式図である。
【図５】物体の複数の騒音源から伝播した音を一つのマイクロホンが検出する様子を示す概略図である。
【図６】音源探査装置に備わる複数のマイクロホンが検出した音から物体上の騒音源の音の強さを求める様子を示す概略図である。
【図７】この発明による遺伝的アルゴリズムによってマイクロホンアレイを設計するための数学モデルを示す図である。
【図８】図７に示す数学モデルに基づいて、最大適応度を持つマイクロホン座標の組合せを異なる世代ごとに表した図である。
【図９】図８に対応して、最大適応度を持つ個体によって推定される音源の位置を、異なる世代毎に、二次元等高線表示（左側）と三次元表示（右側）とで示した図である。
【図１０】図８に示す第５１世代の座標を持つ最適設計で得られたマイクロホンアレイを用いて音源探査をシミュレーションしたときの音圧分布の結果を示す図であって、（ａ）は二次元等高線表示、（ｂ）は（ａ）の三次元表示である。
【図１１】既存の十字状のマイクロホンアレイを用いて音源探査をシミュレーションしたときの音圧分布の結果を示す図であって、（ａ）は二次元等高線表示、（ｂ）は（ａ）の三次元表示である。
【図１２】従来の音源探査装置の概要を説明する図である。
【符号の説明】
１，２０，３０探査装置２，２２表面
３，２３音圧検出素子（マイクロホン）４平板
５音圧検出素子配列体
１０，４０対象物（航空機）１１騒音面
１２，２５，２６，４２騒音源（真音源）１３偽音源
１４仮想音源の音源音圧分布（等高線）
３１開始信号発生装置３２増幅器
３３同時収録装置３４電子計算機（演算装置）３５データ記憶装置３６結果表示装置
ξ 音圧差 φ 音源音圧勾配[0001]
BACKGROUND OF THE INVENTION
The present invention relates to a sound source exploration method used for specifying a sound source such as noise generated from a model placed in a wind tunnel, noise generated from a flying object flying in an arbitrary orbit, and noise generated from various other mechanical devices, and the like. Relates to the device.
[0002]
[Prior art]
Noises of various frequencies are emitted from various positions from flying objects such as aircraft, rotor blades, wind tunnel models placed in an air stream, and other various machines and plant devices. For example, FIG. 4 is a schematic diagram illustrating a conventionally known aircraft noise source, (a) is a schematic diagram showing one half of an aircraft, and (b) is an example of a gear noise source. It is an image figure of no sound pressure line. As shown in FIG. 4A, the noise sources of the aircraft 50 include fan / compressor noise 51 generated mainly in the front and jet noise 52 generated mainly in the rear in relation to the engine. There are flap noise / flap jet interference noise 53 and wing tip vortex sound 54 for the main wing, exhaust interference noise 55 and wing tip vortex sound 56 for the tail wing, and boundary layer noise 57 and gear noise for the fuselage. There are 58. As shown in FIG. 4B by taking the gear noise 58 as an example, the noise source 58a is not necessarily a point-like noise source but has a certain spread, but the sound pressure decreases from the noise source 58a as the distance r increases. A spherical isobar 59 is drawn.
[0003]
In the conventional sound source exploration technique, a method of obtaining a sound incident direction by regularly arranging a plurality of microphones and analyzing a reception signal of the microphones is generally used, and a group of microphones used is called a microphone array or a phased array. There are two important technical elements in the microphone array method. One is an analysis method of data received by a microphone, and the other is a microphone arrangement method.
[0004]
As a data analysis method, a method of performing weighted averaging by applying a weighting coefficient to a reception signal of each microphone is common. The basic idea is as follows.
Sound generated from a sound source on the virtual analysis surface and radiated to free space propagates in a spherical shape and is received by each microphone (coordinates are known). At this time, the sound energy is inversely proportional to the square of the distance r and is accompanied by a time delay τ required to propagate the distance r (the propagation time τ is obtained from the propagation distance r and the sound speed c, τ = r / c). The sound pressure decreases in inverse proportion to the propagation distance r between the sound source and the microphone. As shown in FIG. 5, the sound Pm received by the i-th microphone 62 i among the microphones 62 disposed on the surface 61 of the microphone array 60. _i (T) (t is time) is the sound Ps generated by the jth sound source 64j out of several assumed sound sources 64 on the surface 63 of the aircraft 50. _j Is the distance r from the sound source 64j to the microphone 62i. _{i, j} The time delay τ _{i, j} Is the synthesis of the sound propagated from the sound source 64j for all j, and is expressed by the following equation (1). In this case, the case where the sound pressure of the sound source 64 is zero is included.
Pm _i (T) = αΣ _j {Ps _j (T-τ _{i, j} ) / R _{i, j} } (1)
In consideration of the propagation time (τ = r ′ / c) determined from the propagation distance r ′ between the sound source and the apparent microphone 62 when there is a mainstream, the characteristics of the received sound of each microphone 62 can be obtained by inputting the sound source characteristics. Can be sought.
[0005]
Dividing the surface (sound source surface) where the sound source is to be analyzed into a number of virtual sound source regions, and integrating the microphone reception signal using the propagation distance and propagation time, the sound pressure distribution on the sound source surface from the microphone reception signal Can be calculated backwards. Since the sound pressure distribution is represented by a sound pressure isobar, it is assumed that a sound source exists at the extreme value of the sound pressure isobar. The time domain method integrates received signals as time series signals, and the frequency domain method integrates after frequency analysis.
[0006]
In order to return the signal of the i-th microphone 62i to the j-th focal point (virtual sound source position), the sound pressure is amplified by the distance and the phase is returned by the propagation time, which is expressed by the following equation (2). The
Ps _{i, j} = W _i Xr _{i, j} × Pm _j (T + τ _{i, j} ) ... Formula (2)
Where w _i Is a weight function.
Further, as shown in FIG. 6, when the reception signals of all microphones 62i are integrated for the sound pressure of the jth focal point 64j on the surface of the aircraft 50, the strength Ps of the sound source 64j virtually assumed at the jth focal point. _j Is obtained as equation (3).
Ps _j (T) = Σw _i Xr _{i, j} × Pm _i (T + τ _{i, j} ) / Σw _i ..Formula (3)
In the region where the sound source 64 does not exist, the correlation between the signals to be integrated is weak, so the integrated value decreases, and the difference (sound pressure difference) from the sound pressure at the focal point where the sound source 64 exists increases. In the frequency domain method, a cross spectrum between different microphones is obtained in advance, and the amplitude and phase are corrected according to the distance. For example, the cross spectrum of the kth and lth microphone received signals is represented by C _{i, j} And the distance between the jth focal point and the k, l microphones, respectively, r _{j, k} , R _{j, l} Let c be the sound speed and f be the frequency. And the cross spectrum Cs when the cross spectrum at the microphones k and l is converted into the focus j. _{j, k, l} Is corrected for phase and amplitude as follows.

[0007]
In order to specify the sound source position of a moving body that generates sounds such as trains and automobiles, two sets of two-dimensional microphone arrays are arranged in a vertical plane with respect to the microphone array. It has been proposed to add an intensity probe (see, for example, Patent Document 1). While the sound source position in the traveling direction of the moving body is obtained by weighted averaging of the received signals of the microphone array, the sound intensity vector is obtained from the received signal of the intensity probe, so that the sound source in the direction perpendicular to the traveling direction is obtained. The position can be specified. Although there is an advantage in analyzing a two-dimensional distribution, the combination of an intensity probe and a microphone array becomes a large-scale measurement system and lacks mobility.
[0008]
Based on the frequency spectrum obtained from the received signal of the microphone, the sound wave frequency path between the microphone and the sound source is used to determine the frequency spectrum of the sound source, both frequency spectra are compared, and the sound wave propagation path is optimized. Thus, a method for specifying the sound source position has been proposed (see, for example, Patent Document 2). A genetic algorithm is used as an optimization method. A genetic algorithm is an optimization method that uses gene barrels and the applicability of genes to repeat gene selection, mating, and mutation over several generations to meet target environmental conditions. This is a technique for obtaining an optimal value, that is, an optimal value. Incidentally, examples of the application of the genetic algorithm can be found in the search for an adaptive filter coefficient in active noise control of an air-conditioning refrigeration equipment, the tracking of the adaptive filter coefficient with respect to secular change of vibration noise characteristics in an automobile, and the like. Genetic algorithm is a kind of optimization method, and it is effective even when there is an extreme value less than or equal to the second place compared to the least square method, but it requires some real-time performance such as data processing and control. In some cases, the calculation load is generally large, and it is difficult to put it to practical use with the current computer technology.
[0009]
On the other hand, as shown in FIG. 12, the microphone arrangement method is a one-dimensional arrangement (FIG. 12A), a cross arrangement of a one-dimensional array such as a cross shape, an X shape, or a V shape (FIG. 12B). (D)), combinations with intensity probes (same (e)), circles, arcs, annular arrangements (same (f)), spheres, and the like have been considered. In practice, a two-dimensional sound source distribution is often measured, and for that purpose, a two-dimensional arrangement of microphones, for example, a two-dimensional lattice arrangement (same (g), two-dimensional distributed array (same (h)) is sufficient. .
[0010]
As described above, several measurement methods and microphone placement methods have been proposed, but there are some problems with the conventional techniques regarding the measurement method, auxiliary device, and calibration method using a microphone array for sound source search. That is, first, the arrangement of the microphones greatly affects the analysis result. A two-dimensional analysis of the sound source cannot be performed with a one-dimensional distribution of microphones. In microphones that are regularly arranged such as the lattice arrangement shown in FIG. 12G, the sound of the sound source on the vertical bisector 81 between the two

microphones

80a and 80b is obtained as shown in FIG. On the other hand, since the

microphones

80a and 80b have the same phase, a plurality of

focal points

82a, 82b, 82c,... As shown in FIG. 12 (j), it is known that

sound pressure ridges

84a and 84b are generated in a grid pattern from the sound source position 83, and a noise region 85 spreads around the sound source. A false sound source 91 appears.
[0011]
In the case of a two-dimensional distributed arrangement, only the true sound source sound pressure level is referred to optimize the microphone coordinates, and the difference between the true sound source sound pressure and the false sound source sound pressure for each frequency band and the true sound source coordinate Does not consider deviation. As a result, even if the false sound source sound pressure level is low, the true sound source position may be blurred or vice versa, which will adversely affect the analysis accuracy. In addition, optimization is generally performed by updating the directional filter of each microphone by the least square method or the like, but is not effective when there are many extreme values. Furthermore, the genetic algorithm is excellent in the ability to discriminate the optimal solution even when there are multiple extreme values, but it takes time for the arithmetic processing and is used for designing rather than applications that require real-time performance such as vibration noise control. Is more effective.
[0012]
In order to carry out appropriate noise countermeasures for a target with these noise sources in a short budget, in a short period of time, from where (sound source position) the sound in the frequency band you feel is annoying (noise level) It is necessary to accurately determine whether or not it is generated as 3D information. If the noise generation position, frequency band, sound pressure, etc. to be noticed are identified by performing so-called sound source exploration, the workload of model tests and numerical analysis required for noise cause determination can be reduced, and noise can be reduced. The accuracy of the reduction target, how much the noise source part is modified, whether or not a reduction device can be added, whether or not usage conditions have been changed, the scale of measures such as ambient soundproofing, and whether or not to radically change the machinery High information can be obtained, the time required to determine these measures can be shortened, and more accurate and effective measures can be selected.
[0013]
[Patent Document 1]
JP-A-08-015003 (6th column, paragraph [0023] to 8th column, paragraph [0034], FIGS. 1 to 3)
[Patent Document 2]
JP 09-021863 A (11th column, paragraph [0056] to 12th column, paragraph [0058])
[0014]
[Problems to be solved by the invention]
Therefore, when the true sound source is obtained by distinguishing the true sound source and the false sound source from the sound pressure distribution obtained by calculating back in the estimated sound source region on the analysis surface, the sound pressure detection can be obtained efficiently and accurately. There is a problem to be solved in terms of finding the optimum arrangement of elements.
[0015]
An object of the present invention is to obtain a sound source search method capable of efficiently searching for an optimum arrangement of a plurality of sound pressure detection elements for detecting sound pressure propagated from a sound source, and to accurately search a true sound source from a false sound source, and Is to provide a device.
[0016]
[Means for Solving the Problems]
In order to solve the above-described problem, a sound source search method according to the present invention includes a plurality of sound pressure detections in which sound pressures of sound propagated from sound sources existing in space are distributed on a two-dimensional plane or a three-dimensional curved surface. The sound source sound pressure distribution of the virtual sound source estimated for the sound source is obtained based on the sound pressure signal detected by the element and detected by each sound pressure detection element, and based on the evaluation function determined from the sound source sound pressure distribution And determining the optimal placement of the sound pressure sensing element on the surface.
[0017]
Also, the sound source exploration device according to the present invention has a two-dimensional plane or a three-dimensional curved surface, and a plurality of sound pressure detection elements for detecting the sound pressure of sound propagated from a sound source existing in space on the surface is dispersed. Obtaining the sound source sound pressure distribution of the virtual sound source estimated for the sound source based on the arranged sound pressure detecting element array and the sound pressure signal detected by each sound pressure detecting element, and determining from the sound source sound pressure distribution And an arithmetic unit that obtains the optimum arrangement of the sound pressure detecting elements on the surface based on the evaluation function.
[0018]
According to the sound source exploration method and sound source exploration device according to the present invention, the sound pressure of the sound propagated from the sound source existing in the space is a plurality of distributed and arranged on the surface of the two-dimensional plane or the three-dimensional curved surface of the sound source exploration device. It is detected by the sound pressure detecting element. Based on the sound pressure signal detected by each sound pressure detection element, the sound source sound pressure distribution of the virtual sound source estimated in the space for the sound source is obtained, for example, by back calculation. From the sound source sound pressure distribution obtained based on the detected sound pressure signal, it is generally difficult to uniquely determine the true sound source due to a phenomenon such as the appearance of a false sound source. Based on the evaluation function determined from the above, the optimal arrangement on the surface of the sound pressure detecting element can be obtained by calculation, and the sound is detected by the sound pressure detecting element arranged in the optimum arrangement thus obtained. By obtaining the sound source sound pressure distribution of the virtual sound source, it becomes possible to easily and efficiently determine the true sound source by the sound source search device.
[0019]
In the sound source search method and the sound source search device, the evaluation function is a non-negative function for a sound pressure difference between the true sound source position and the false sound source position of the sound source obtained from the sound pressure distribution and a sound pressure gradient near the true sound source position. The optimal arrangement of the sound pressure detecting elements can be determined as the arrangement coordinates of the sound pressure detecting elements when the evaluation function is maximized. When the evaluation function is defined in this way, the sound pressure difference between the true sound source position and the false sound source position of the sound source is larger, and the sound pressure gradient near the true sound source position is steeper, and the function value of each non-negative function is The evaluation function defined by the linear combination or product thereof becomes larger. Therefore, when the evaluation function is maximized, the true sound source and the false sound source can be most easily distinguished, and the arrangement coordinates of the sound pressure detection elements at that time can be set as the optimum arrangement of the sound pressure detection elements.
[0020]
In a sound source search method and a sound source search device in which the evaluation function is a linear combination or product of non-negative functions, a calculation process for maximizing the evaluation function includes a number of arrangements of the arrangement coordinates of the sound pressure detection elements arranged at genetic loci. Applying a genetic algorithm that forms a group from individuals and repeats the generation through the process of selection, crossover and mutation of the individuals in the group for generations, and uses the evaluation function as the fitness of the selection Can be used. Individuals who give a low function value to the evaluation function by applying a genetic algorithm that repeats the transition through the process of individual selection, crossover and mutation over generations and using the evaluation function as the fitness of selection ( The arrangement of sound pressure detection elements is naturally eliminated, and the proportion of high-quality individuals that give a high function value to the evaluation function increases, and the optimal arrangement of sound pressure detection elements is efficient without being confused by false sound sources. Can get to.
[0021]
In this sound source exploration method and sound source exploration device, all the sound pressure signals detected by the sound pressure detecting element are recorded, and the sound pressure signal used for obtaining the sound source sound pressure distribution among the recorded sound pressure signals Is selected according to the frequency to be searched, the method of determining the optimum arrangement by the evaluation function can be applied. All sound pressure signals detected by a large number of sound pressure detection elements are recorded in advance, and in the analysis, a part of the recorded sound pressure signals is selected according to the frequency of particular attention to be searched. Select to determine the sound source sound pressure distribution. Which sound pressure signal to use to determine the sound source sound pressure distribution among the recorded sound pressure signals can be replaced with the problem of where to place some of the sound pressure detection elements to be adopted, It is possible to apply a method for determining an optimum arrangement by using an evaluation function.
[0022]
In this sound source search method, it is possible to start simultaneously recording the sound pressure signal detected by the sound pressure detection element as a discretization signal based on the start signal from the start signal generator. In the sound source exploration device, a simultaneous recording device that simultaneously records the sound pressure signal detected by the sound pressure detection element as a discretization signal, and a start signal generation device that outputs a start signal for starting recording by the simultaneous recording device Can be provided. For example, it may be necessary to start sound pressure detection in accordance with the passage, such as when searching for the noise of an object such as an aircraft or a vehicle that passes through. In such a case, the sound pressure signal detected by the sound pressure detecting element based on the start signal from the start signal generating device that detects the passage of the object with light faster than the sound or image means is a discrete signal. As a simultaneous recording.
[0023]
In the sound source search method and the sound source search device, the sound pressure detection element may be a microphone or a pressure sensor. A transducer capable of converting sound pressure into an electric signal is preferable.
[0024]
DETAILED DESCRIPTION OF THE INVENTION
Hereinafter, embodiments of a sound source searching method and apparatus according to the present invention will be described with reference to the accompanying drawings. FIG. 1 is an explanatory diagram for explaining an outline of a sound source search method according to the present invention. FIG. 1A is an example of an analysis model showing a positional relationship between a sound source surface where a sound source exists and a surface where a sound pressure detecting element is arranged. (B) is a sound pressure distribution showing a sound source analysis result, and (c) is a sound pressure graph taking one section thereof.
[0025]
According to the explanatory diagram shown in FIG. 1, the sound source search method according to the present invention is carried out by following the procedures shown in the following (1) to (7).
(1) The sound source search method and apparatus are characterized in that each sound pressure detection element coordinate is designed analytically. As an analysis model, as shown in FIG. 1A, the arrangement coordinates of a plurality of microphones 3 which are sound pressure detection elements are placed on the surface (in this example, a two-dimensional flat plane) 2 of the sound source exploration device 1. Set the initial value. Further, the coordinates and sound pressure of the true sound source 12 are set on the sound source surface 11 of the object 10 such as an aircraft, a vehicle, or a factory as a sound source search target. That is, the true sound source 12 is placed at known coordinates in the space. Further, a determination value ε to be used for convergence determination in later calculation is set.
[0026]
(2) Calculate the sound emitted from the true sound source 12 by calculating the received signal at each microphone 3.
The sound arriving at each microphone 3 from the true sound source 12 takes a time τ in which the amplitude is inversely proportional to the distance r from the true sound source 12 to the microphone 3 and the sound propagates through the distance r, considering spherical radiation characteristics. There is a delay. Therefore, sounds radiated from the same true sound source 12 are received by the microphones 3 with different amplitudes and different time delays (referred to as phases). The relationship between the true sound source 12 and the sound pressure signal of the microphone 3 is as already described with reference to FIG. Taking this relationship in reverse, as already shown in FIG. 6, the reception signals of all microphones 3 are amplified by the distance to the assumed sound source equipment position (referred to as the focal point) where the true sound source 12 is assumed, and the time is increased. The focus sound can be calculated by advancing and accumulating.
[0027]
(3) The sound source is searched by arranging a virtual sound source on the sound source surface 11, calculating the sound pressure, and calculating the sound pressure contour of the sound source surface 11.
When the above calculation is performed by subdividing the analysis surface including the true sound source 12 and regarding each portion as a focal point and obtaining the maximum sound pressure at the focal point, as shown in FIG. A contour line 14 is obtained. The contour line 14 should have the maximum value at the position of the known sound source (that is, the true sound source 12). However, unlike the known sound source position, an extreme value is generated at a position where the sound source is not originally placed, and the false sound source 13 is generated. appear.
[0028]
(4) Calculate two parameters.
As shown in FIG. 1C, one parameter is the known sound source position, that is, the sound pressure difference ξ between the true sound source 12 and the false sound source 13, and the other parameter is the true sound source sound pressure gradient φ. In measurement, it is desirable to increase the sound pressure difference ξ between the sound pressure (maximum value) at the position of the true sound source 12 and the sound pressure (maximum value) at the position of the false sound source 13, that is, to increase the signal-to-noise ratio. On the other hand, the steep slope φ of the sound pressure contour near the maximum value corresponding to the true sound source 12 means that the error range of the sound source position is small, which is desirable in terms of position accuracy.
[0029]
(5) The evaluation function H (ξ, φ) is calculated.
The following evaluation function (also called an objective function) expressed as a linear combination sum or product of these non-negative functions is defined with both the sound pressure difference ξ and the true sound source sound pressure gradient φ as variables.
H (ξ, φ) = αF (ξ) + βG (φ) or
H (ξ, φ) = αF (ξ) × G (φ)
Where F (x), G (x): positive non-negative function, α, β: positive number
In the present invention, the updating of the microphone coordinates and the calculation of the evaluation function are repeated so as to maximize the evaluation function, and the microphone coordinates when the evaluation function takes the maximum value are determined as the arrangement of the sound source exploration microphones.
[0030]
(6) Determine convergence. The convergence determination is made based on whether or not the change in H (ξ, φ) is equal to or less than ε. If this condition is satisfied, it is determined that the convergence has occurred.
(7) If the convergence determination is No, the coordinates of the microphone 3 are changed and the process returns to (2). If Yes, the process ends and the coordinates of the microphone 3 are determined.
[0031]
In the above example, the surface 2 of the sound source exploration device 1 is a two-dimensional flat plate surface in this example, but is not limited thereto, and may be a three-dimensional curved plate surface such as a spherical surface. Further, although the microphone 3 is used as the sound pressure detection element, the present invention is not limited to this, and a transducer capable of outputting a detection signal as an electric signal from a pressure sensor or the like is preferable. In order to obtain a search accuracy that does not interfere with practical use, the number of microphones 3 to be embedded is preferably 4 or more.
[0032]
A genetic algorithm can be used in the calculation process for optimizing the evaluation function. The genetic algorithm repeats three processes (selection, crossover, mutation) for individual genes included in a group over several generations, thereby generating a group (ie, optimal solution) having a highly adaptable gene. This technique is well known as a solution for the traveling salesman problem, knapsack problem, and the like, and is also applied to aircraft aircraft shape optimization problems, manufacturing process management, and financial / investment applications.
[0033]
In the genetic algorithm, there are N individuals in a group, and the number of individuals in the group is unchanged at each generation. In determining the placement of microphones, each individual has information on the microphone coordinate group. Specifically, the number of microphone coordinates used for the microphone array is arranged in the gene locus. Dozens or hundreds of these individuals are prepared to form a group. The individuals in the group are changed through three processes: selection, crossover, and mutation, and this process is repeated for generations. At this time, the fitness of selection is the evaluation function H that is uniquely determined by using the sound pressure difference ξ and the sound source sound pressure gradient φ between the true sound source sound pressure and the false sound source sound pressure as variables.
[0034]
Each individual has M (= number of microphones) gene loci, and microphone coordinate information is written in each gene locus. For example, when the coordinate information written in the i-th gene locus is (x, y) = (5, 3) (unit: mm), in the case of integer display, (XY) is (53). In the case of binary number display, it becomes (01010011). This can be connected to M microphones to represent the individual as follows. In the real number display, it is the j-th individual (910... 53... 13), and the middle 53 is a description regarding the i-th microphone information. On the other hand, in the binary representation, it is the jth individual (10011010... 0101010011... 01101). The initial value is set for all loci of all individuals using random numbers.
[0035]
Under the above conditions, the following three processes are repeated for each generation.
(1) Selection
Since each individual has coordinate information of M microphones, the sound pressure contours of the analysis surface can be calculated for each individual, and an evaluation function (ie, fitness) can be obtained. Corresponding to 1. Individuals can be arranged in descending order of fitness (referred to as ranking). Under the initial conditions, a lower limit of fitness is set, and individuals with fitness below the lower limit are discarded as ineligible individuals. In order to keep the number of individuals in the group, only the number of good individuals that are discarded are duplicated. For example, the duplication method may be determined from the ranking ratio. For example, if a group has 10 individuals and 3 are discarded, the top three that are matched for new replication are replicated. Individuals with low fitness and poor evaluation function (microphone placement) are naturally excluded by selection, and the proportion of individuals with good fitness increases. In the process of selection, an act of elite selection can be added. This is an act of extremely reducing or setting the probability of crossover or mutation below to zero for an individual that exceeds a preset upper limit of fitness.
[0036]
(2) Crossover
This is an act of exchanging a part of loci (the position and length can be selected stochastically) between different members selected with a preset crossover probability. Crossover changes the tendency of individuals in the group and leads to exploring the maximum value from multiple extreme values. As an example of crossover, when the following individuals A and B exchange the loci of the portion described in the middle, and each becomes an individual A ′ individual B ′, the following occurs.
Individual A Individual B
(0110 ... 11010 ... 0101) (1000 ... 10001 ... 1110)
→
Individual A 'Individual B'
(0110 ... 10001 ... 0101) (1000 ... 11010 ... 1110)
[0037]
(3) Mutation
An individual selected with a preset mutation probability is an act of exchanging a part or all of a gene locus. Mutations activate group population changes. Various methods are used for the value of the locus, such as using a random number or converting the bit value if it is binary.
An example is shown below.
Individual A Individual A '
(0110 ... 11010 ... 0101) → (1001 ... 00101 ... 1010)
As the above three processes are repeated for several hundred to several thousand generations, the proportion of excellent individuals increases, and the optimum microphone arrangement is determined.
[0038]
Next, selective use of the sound pressure detection element will be described with reference to FIG. FIG. 2 is an explanatory view schematically illustrating an outline of selective use of a sound pressure detecting element in the sound source searching method and apparatus according to the present invention. In the above example, the number of sound pressure detection elements used for sound source exploration and their coordinates are optimized according to the detection target before actual measurement. However, in the example shown in FIG. Sound pressure data is recorded by a large number of arranged elements for detecting sound pressure, and the number and coordinates of microphones are selected according to the frequency used in the analysis at the time of noise analysis. That is, as shown in FIG. 2A, the sound source exploration device 20 includes a large number of microphones 23 that are dispersedly arranged (lattice arrangement in the illustrated example) at known coordinates on a planar surface 22. Yes. For example, noise (schematic sound wave 24) is simultaneously measured using all 100 microphones 23 that are dispersedly arranged, and all measured noise data is recorded.
[0039]
In the selective use of the sound pressure detecting element, noise is measured at many points anyway, and the microphone 23 used for analysis is selected at the noise analysis stage. Using the combination of microphones and the microphone coordinates of the combination, a microphone to be used for analysis for optimizing the evaluation function is selected for each frequency. For example, ten microphones 23 are selected at a frequency of 100 Hz, and noise data measured by the selected microphones 23 is employed in actual analysis. Similarly, the microphone 23 is selected by a selection method of using 50 at 1000 Hz and 100 at 2000 Hz. As a method for selecting the microphone 23, the method using the above-described evaluation function is used, and the microphone 23 having the optimum arrangement may be selected.
[0040]
Specifically, as shown in FIG. 2B, the position accuracy is not high by using a small number of microphones 23, but the positions of the noise sources 25 and 26 are analyzed and obtained according to the frequency bands w1 and w2. It is done. Assuming that the noise sources 25 and 26 exist at the positions obtained in FIG. 2B, the evaluation function H is applied by applying the above-described method using the evaluation function H that is a function of the sound pressure difference ξ and the sound pressure gradient φ. The combination of the microphones 23 that maximizes is selected. As shown in FIG. 2C, the precise noise source search is performed again using the selected combination of the microphones 23, and the sound source search can be performed with high positional accuracy for each of the noise sources 25 and 26.
[0041]
A simulation of sound source search by a microphone array using a genetic algorithm according to the present invention was performed. The simulation conditions and results are shown below. FIG. 7 is a diagram illustrating a mathematical model for designing a microphone array. In FIG. 7, reference numeral 70 denotes a measurement surface on which 41 microphones 72 exist, and coordinates are determined with reference to the origin 71. One microphone 72 is arranged at the origin 71, and five identical regions where eight microphones 72 are arranged around the origin 71 are set so as to exist at equal intervals in the circumferential direction. The polygon is not limited to a pentagon. In FIG. 7, reference numeral 75 denotes an analysis surface where the

virtual sound sources

76 and 77 exist, which is parallel to the measurement surface 70, and the inter-plane distance h is set to 4 m. The coordinates of the sound source 76 are set to [0.0, 0.0, 4.0], and the coordinates of the sound source 77 are set to [−0.7, −0.7, 4.0]. It is assumed that both

sound sources

75 and 77 generate a pure sound with a sound pressure of 1 at 2000 Hz. The analysis surface 75 may not actually be parallel to the measurement surface 70, and a three-dimensional sound source search can be performed by selecting a plurality of analysis surfaces.
[0042]
The parameters of the genetic algorithm are a crossover probability of 0.8, a mutation probability of 0.05, an elite probability of 0.05, a group population of 40, and an individual gene length of 144. The fitness H is F for the difference between the true sound source sound pressure and the false sound source sound pressure, G for the true sound source sound pressure gradient, and P for the penalty function (minimum microphone interval, microphone array dimensions, maximum sound pressure value deviation). , H = F × G × P.
[0043]
FIG. 8 is a diagram showing individuals having the maximum fitness (that is, combinations of microphone coordinates) for different generations. FIGS. 8A to 8D show microphone coordinates in the first generation, the sixth generation, the 21st generation, and the 51st generation, respectively, and the fitness levels are 23.60 and 29.35, respectively. 42.66 and 45.42. It can be seen that as the generation progresses and evolves, the fitness increases and the arrangement of microphones is dispersed. Further, FIG. 9 shows a two-dimensional contour display of the position of the sound source estimated by the individual having the maximum fitness for each different generation (first, sixth, twenty-first, and fifty-first generations) similar to FIG. It is the figure shown by (left side) and three-dimensional display (right side). A small point in the figure represents a microphone and corresponds to the coordinates of the microphone shown in FIG. As the generation progresses, the sound pressure difference between the true sound source sound pressure and the surrounding noise pressure increases, and the sound pressure gradient at the sound source position becomes steep.
[0044]
FIG. 10 is a diagram showing the result of sound pressure distribution when a sound source exploration was simulated using the microphone array (using the 51st generation coordinates shown in FIG. 8) obtained by the above-described optimum design. (A) is a two-dimensional contour display, and FIG. 10 (b) is the three-dimensional display of (a). FIG. 11 is a diagram showing the result of sound pressure distribution when sound source exploration is simulated using an existing cross-shaped microphone array, and FIG. 11A is a two-dimensional contour line display. (B) is the three-dimensional display of (a). The points in FIGS. 10 and 11 indicate the microphone positions. The sound source to be used is a point sound source having a frequency of 2000 Hz and a sound pressure of 1, and coordinates are in [0.5, 0.5, 4]. In the cross-shaped microphone array shown in FIG. 11, a sound pressure peak appears as a false sound source in the arrangement direction of the array around the sound source position. On the other hand, when measurement and analysis are performed using a microphone array optimally designed by the genetic algorithm according to the present invention, an increase in sound pressure is suppressed at positions other than the sound source position, and a sound pressure gradient at the sound source position is steep. You can see that it is.
[0045]
FIG. 3 is a schematic diagram showing an outline of an aspect of noise source search using the sound source search apparatus according to the present invention when the noise search target is an aircraft in flight or traveling. The sound source exploration device 30 includes a large number of microphones 3 on a flat surface 2 like the sound source exploration device 1 shown in FIG. The flat plate 4 that provides the surface 2 and the multiple microphones 3 define the sound pressure detection element array 5 in the sound source exploration device 30. In this system, it is assumed that the noise source search target is the aircraft 40 in flight or traveling, and the jet noise of the engine 41 is the true sound source 42. The sound source exploration device 30 includes a start signal generation device 31 that can optically capture the moment when the aircraft 40 passes through an explorable region as a reflection of a laser beam and output a capture signal.
[0046]
In the sound source exploration device 30, the sound pressure signals from all the microphones 3 are amplified by the amplifiers 32 under the control performed by the electronic computer 34 in response to the captured signal from the start signal generating device 31, and the amplified sound pressures are amplified. The signal is simultaneously recorded as a discretized signal by the simultaneous recording device 33. The recorded data is analyzed by the electronic computer 34, stored in the data storage device 35, and displayed on the result display device 36 such as a display. By analyzing the data input to the simultaneous recording device 33 in conjunction with the sound pressure signal from the microphone 3, the noise source 42 can be specified for each instantaneous position of the aircraft 40. The start signal generator 31 plays the role of a starter for starting simultaneous recording. However, the passage of the aircraft 40 is not limited to capturing using laser light, but optically using a CCD image and acoustically. In addition, the ultrasonic reflection can be detected by converting a pressure change accompanying the passage of the aircraft 40 into a voltage. The electronic computer 34 obtains the sound source sound pressure distribution of the virtual sound source estimated for the sound source 42 based on the sound pressure signal detected by each sound pressure detection element 3, and the sound pressure based on the evaluation function determined from the sound source sound pressure distribution. It functions as an arithmetic unit for the sound source exploration device 30 for obtaining the optimum arrangement on the surface 2 of the detection element 3.
[0047]
The sound source search according to the present invention has been described by taking the search for a noise source in an aircraft as an example, but the application of the present invention is not limited to this example, and the noise source search for a model or an actual machine in a wind tunnel experiment, an automobile, a truck, etc. Investigating noise sources of objects with moving sound sources such as vehicles, railway vehicles, noise sources of large equipment such as factories, plants, power generation facilities, buildings, etc., various rotary machines such as axial flow or centrifugal compressors, blowers, turbines, etc. It is clear that the present invention can be applied to noise source search for OA or noise abnormality diagnosis of OA equipment and home appliances.
[0048]
【The invention's effect】
As described above, the sound source search method and apparatus according to the present invention include a plurality of sound pressure detection elements in which the sound pressure of sound propagated from a sound source existing in space is distributed on the surface of a two-dimensional plane or a three-dimensional curved surface. The sound source sound pressure distribution of the virtual sound source estimated for the sound source is obtained based on the sound pressure signal detected by each sound pressure detecting element, and the sound pressure detecting element of the sound pressure detecting element is determined based on the evaluation function determined from the sound source sound pressure distribution. Since the optimum arrangement on the surface is determined, the sound source is detected by the sound pressure detecting element arranged in the optimum arrangement thus obtained, and the sound source sound pressure distribution of the virtual sound source is obtained. By distinguishing from the false sound source, it is possible to easily and efficiently determine the true sound source by the sound source search device. As described above, since it becomes possible to obtain highly accurate information as the three-dimensional information about the frequency band of the noise source, the position of the noise source, and the magnitude of the noise, the work load for noise cause determination is reduced, and the noise is reduced. More accurate and effective noise countermeasures such as reduction targets, degree of modification of noise source parts, availability of reduction devices, presence / absence of usage condition changes, scale of measures such as ambient soundproofing, whether or not to radically change mechanical devices It can be selected in a short time.
[Brief description of the drawings]
FIG. 1 is an explanatory diagram for explaining an outline of a sound source search method according to the present invention;
FIG. 2 is an explanatory diagram schematically illustrating an outline of selective use of a sound pressure detecting element in a sound source searching method and apparatus according to the present invention.
FIG. 3 is a schematic diagram showing an outline of an aspect of noise source search using a sound source search device according to the present invention when a noise search target is an aircraft in flight or traveling.
FIG. 4 is a schematic diagram showing aircraft noise sources.
FIG. 5 is a schematic diagram showing how a single microphone detects sound propagated from a plurality of noise sources of an object.
FIG. 6 is a schematic diagram showing a state in which the sound intensity of a noise source on an object is obtained from sounds detected by a plurality of microphones provided in the sound source search device.
FIG. 7 is a diagram showing a mathematical model for designing a microphone array by a genetic algorithm according to the present invention.
8 is a diagram showing combinations of microphone coordinates having the maximum fitness for different generations based on the mathematical model shown in FIG. 7;
FIG. 9 is a diagram corresponding to FIG. 8 showing the position of the sound source estimated by an individual having the maximum fitness in two-dimensional contour display (left side) and three-dimensional display (right side) for each different generation. It is.
10 is a diagram showing a result of sound pressure distribution when a sound source exploration was simulated using a microphone array obtained by an optimum design having the coordinates of the 51st generation shown in FIG. 8, and FIG. Dimensional contour display, (b) is the three-dimensional display of (a).
FIGS. 11A and 11B are diagrams showing the results of sound pressure distribution when sound source exploration is simulated using an existing cross-shaped microphone array, where FIG. 11A is a two-dimensional contour display, and FIG. 11B is a diagram of FIG. 3D display.
FIG. 12 is a diagram for explaining the outline of a conventional sound source search device;
[Explanation of symbols]
1,20,30

Exploration device

2,22 Surface
3,23 Sound pressure detection element (microphone) 4 Flat plate
5 Sound pressure detection element array
10, 40 Object (aircraft) 11 Noise surface
12, 25, 26, 42 Noise source (true sound source) 13 False sound source
14 Sound source sound pressure distribution (contour) of virtual sound source
31 Start signal generator 32 Amplifier
33 Simultaneous recording device 34 Electronic computer (arithmetic unit) 35 Data storage device 36 Result display device
ξ Sound pressure difference φ Sound source sound pressure gradient

Claims

The sound pressure of sound propagated from a sound source existing in the space is detected by a plurality of sound pressure detection elements distributed on the surface of a two-dimensional plane or a three-dimensional curved surface, and the sound pressure detected by each of the sound pressure detection elements A sound pressure distribution of the virtual sound source estimated for the sound source based on a signal is obtained, and a sound pressure difference and a sound pressure gradient in the sound source sound pressure distribution are determined on the surface of the sound pressure detecting element based on an evaluation function. A sound source exploration method comprising determining an optimal placement.

The evaluation function is a linear combination or product of a non-negative function for a sound pressure difference between a true sound source position and a false sound source position of the sound source and a non-negative function for a sound pressure gradient in the vicinity of the true sound source position obtained from the sound pressure distribution. The sound source search method according to claim 1, wherein the optimal arrangement of the sound pressure detection elements is determined as an arrangement coordinate of the sound pressure detection elements when maximizing the evaluation function.

In the calculation process of maximizing the evaluation function, a group is formed from a large number of individuals in which the arrangement coordinates of the sound pressure detection elements are arranged at a genetic locus, and selection, crossover, and mutation of the individuals in the group are performed. The sound source search method according to claim 2, wherein a genetic algorithm that repeats the transition through the process consisting of generations is applied, and the evaluation function is used as the fitness of the selection.

All the sound pressure signals detected by the sound pressure detection element are recorded, and the sound pressure signal used for obtaining the sound source sound pressure distribution among the recorded sound pressure signals is selected according to the frequency to be searched. The sound source search method according to any one of claims 1 to 3, wherein a method of determining the optimum arrangement by the evaluation function is applied when performing.

The sound pressure signal detected by the sound pressure detecting element is discretized based on the captured signal from the start signal generating device that receives a physical signal accompanying the passage of the object having the sound source and outputs the captured signal. 2. A sound source search method according to claim 1, further comprising starting simultaneous recording as a signal.

The sound source detection method according to claim 1, wherein the sound pressure detection element is a microphone or a pressure sensor.

A sound pressure detection element array having a surface of a two-dimensional plane or a three-dimensional curved surface and a plurality of sound pressure detection elements for detecting the sound pressure of sound propagated from a sound source existing in space on the surface; And obtaining the sound source sound pressure distribution of the virtual sound source estimated for the sound source based on the sound pressure signal detected by each sound pressure detecting element, and using the sound pressure difference and the sound pressure gradient in the sound source sound pressure distribution as variables A sound source exploration device comprising an arithmetic device for obtaining an optimum arrangement of the sound pressure detection elements on the surface based on a function.

The evaluation function is a linear combination or product of a non-negative function for a sound pressure difference between a true sound source position and a false sound source position of the sound source and a non-negative function for a sound pressure gradient in the vicinity of the true sound source position obtained from the sound pressure distribution. The sound source exploration device according to claim 7, wherein the optimum arrangement of the sound pressure detection elements is determined as an arrangement coordinate of the sound pressure detection elements when the evaluation function is maximized.

In the calculation process of maximizing the evaluation function, a group is formed from a large number of individuals in which the arrangement coordinates of the sound pressure detection elements are arranged at a genetic locus, and selection, crossover, and mutation of the individuals in the group are performed. The sound source search device according to claim 8, wherein a genetic algorithm that repeats the transition through the process for generations is applied, and the evaluation function is used as the fitness of the selection.

All the sound pressure signals detected by the sound pressure detection element are recorded, and the sound pressure signal used for obtaining the sound source sound pressure distribution among the recorded sound pressure signals is selected according to the frequency to be searched. The sound source exploration device according to claim 7, wherein a method for determining the optimum arrangement by the evaluation function is applied when performing the operation.

Since it further comprises a simultaneous recording device that simultaneously records the sound pressure signal detected by the sound pressure detection element as a discretized signal, and a start signal generation device that outputs a start signal for starting recording by the simultaneous recording device. The sound source exploration device according to claim 7.

The sound source detection device according to claim 7, wherein the sound pressure detection element is a microphone or a pressure sensor.