JP6788272B2 - 音源の検出方法及びその検出装置 - Google Patents

Description

この発明は、音源の検出方法及びその検出装置に係り、詳しくは音響インテンシティーの測定結果を利用する音源の検出方法及びその検出装置に関するものである。

三次元空間における音の振る舞いを測定する手法としては、以下に例示する方法が知られている。
例えば、直交する３軸方向（ｘ，ｙ，ｚの軸方向）に対する音響インテンシティーを測定し、その測定結果を利用することにより、どの方向からどのタイミングでどのような音が到来しているかというようにコンサートホールの反射音などを解析する方法１（下記非特許文献１、図１５参照）がある。
一方、音響インテンシティーではなく音圧を測定することにより、コンサートホールの反射音などを解析する方法２（下記非特許文献２、図１６参照）と方法３（下記非特許文献３、図１７参照）もある。

上記方法１は、音響インテンシティーの測定結果を利用しているため測定精度が高くなるという長所がある。また、上記方法２，３は、有用な音源情報、すなわち反射を与えている音源いわゆる虚音源（図１６や図１７に示す〇印）が抽出されるため解りやすいという長所がある。
このことから、例えば、上記方法１の音響インテンシティーを測定した結果に対して、上記方法２，３のような反射音（虚音源）を抽出する手法が適用できれば、精度の良い有用な音源情報の解析が可能になりそうである。

ちなみに、上記方法２，３における虚音源の抽出に関する基本的な考え方は、マイクロホンで測定される音圧波形のピークの情報を用いて音源の存在確定とその座標推定を行うものである。
具体的には、図１８に例示するように、ｘ−ｙ平面空間（ｚ＝０）の異なる位置にある２つの音源Ｓ₁，Ｓ₂について、測定点（原点）Ｏの近傍の異なる位置に設置した３つのマイクロホンＰ₁〜Ｐ₃により測定した場合、その測定される音圧波形は、そのタイミングと振幅が互いに異なる。このため、マイクロホンＰ₁〜Ｐ₃で得られる各音圧波形におけるピークを検知することにより、音源の存在を確定し、また、その各音圧波形におけるピークのレベルと時間差の組み合わせから音源の座標を推定するようにしている。

そこで、本発明者は、上記方法２，３における音源の抽出手法ともいえる「ピーク検知」の手法を、上記方法１における音響インテンシティーの測定結果を利用した音源の解析方法に応用するという新たな方法４を案出し、その方法４の有用性を調べた。ここでは、その案出した音源情報（音源の位置）を抽出する方法について、「ピーク検知」を用いて音響インテンシティーから音源情報を検出する方法４と称することにする。

この「ピーク検知」を用いて音響インテンシティーから音源情報を検出する方法４は、以下に示す「音響インテンシティーの測定」、「音源の抽出」及び「音源座標の推定」という処理工程をこの順番で行うという構成を採用している。

＜音響インテンシティーの測定＞
（１） 受音点で直交する３軸方向（ｘ，ｙ，ｚ）の音響インテンシティーＩ_x，Ｉ_y，Ｉ_z（単位：Ｗ／ｍ²）を測定する。
（２） 音響インテンシティーＩ_x，Ｉ_y，Ｉ_zの合計値であるＩ_all＝√（Ｉ_x ²＋Ｉ_y ²＋Ｉ_z ²）を算出する（√（）は括弧内の値の平方根を示す。以下同様。）。
＜音源の抽出＞
（３） 上記合計値の音響インテンシティーの波形におけるピークが音源を表していると仮定する。
（４） 上記合計値のピークを時間窓幅Δｔ［秒］ごとに抽出する。
（５） （４）で抽出されたピークから音源の大きさ及び距離を算出する。
＜音源座標の推定＞
（６） （４）で抽出されたピークに対応する上記Ｉ_x，Ｉ_y，Ｉ_zにおけるピークを時間窓幅Δｔ［秒］ごとに抽出する。
（７） 上記Ｉ_x，Ｉ_y，Ｉ_zにおけるピークの符号付の大きさと（５）で算出した音源の距離とから音源の座標を推定する。

図１９に、この検出方法４の概念と各工程における結果の一例を示す。図１９中における符号Ｐは受音点、符号Ｓ₁は実際の音源を表している実音源、符号Ｓ₂〜Ｓ₇は反射音を表している虚音源をそれぞれ示す。
また、上記（１）〜（５）の詳細については、後述する実施の形態において説明する。
上記（７）における音源の座標は、Ｉ_x，Ｉ_y，Ｉ_zから抽出されたピークの符号付の大きさ±｜Ｉ_px｜，±｜Ｉ_py｜，±｜Ｉ_pz｜にて表される原点Ｏからの方位（±｜Ｉ_px｜，±｜Ｉ_py｜，±｜Ｉ_pz｜）に対して、上記（５）で算出した合計値Ｉ_allのピークが与える距離ほど離れた地点として推定される。

この検出方法４を試行したときの結果の一例を図２０に示す。
この結果の一例では、推定対象としている実音源Ｓ1を、図２０に示されるようにｙ軸上（ｘ＝０，ｚ＝０となる位置）に配置した設定としている。このため、このときのｘ方向とｚ方向における音響インテンシティーＩx，Ｉzは、本来Ｉx＝０、Ｉz＝０になるはずである。
ところが、測定の結果では、音響インテンシティーＩx，Ｉzに多少の残差波形が確認されている。
したがって、この「ピーク検知」を用いた検出方法４では、それらの残差波形のピークをｘ方向及びｚ方向の座標として誤って認識することになる。この結果、最終的には、図２０に示されるように、検出（推定）した音源Ｓ1の位置（黒丸）が、ｙ軸上にあるはずの実音源Ｓ1の正規の位置（白丸）からずれた場所に推定されてしまうことになる。このような結果は、音源位置の推定精度としては低い結果であると言わざるを得ない。

この検出方法４における音源位置の検出ずれは、上記方法２，３において音圧波形に対して開発されたピーク検知の手法を音響インテンシティーの波形に対してそのまま流用したことによる矛盾から生じた結果であると考えることができる。
つまり、音響インテンシティーは、一般にどの方向からどの程度の大きさの音が到来しているかを表す物理量である。このため、例えば、２つの音源Ｓ₁，Ｓ₂が、図２１に示されるようなｘ−ｙ平面空間（ｚ＝０）の各位置（特に音源Ｓ₂はｙ軸上の位置）にそれぞれ配置されている場合、音源Ｓ₂から原点Ｏに到来する音のｘ方向の音響インテンシティーＩ_xはＩ_x＝０となり、ピーク情報が存在しないことになる。このように音響インテンシティーの波形のピーク情報から音源の位置を検出しようとしても、音源の座標を得ることができない場合がある。
したがって、音響インテンシティーを対象とした解析においては、その波形のピーク情報から音源情報（位置）を得られない音源の位置が存在することがわかる。

R.W.Guy and A.Abdou, "A measurement system and method to investigate the directional characteristics of sound fields in enclosures," Noise Control Engineering Journal, 41, 8-18 (1994). Y. Yamasaki and T. Itou, "Measurement of spatial information in sound field by closely located four point microphone method,'' Journal of the Acoustical Society of Japan. (E), 10, 101-110 (1989). K.Sekiguchi, S.Kimura and T.Hanyu, "Analysis of sound field on spatial information using a four-channel microphone system based on regular tetrahedron peak method,'' Applied Acoustics, 37,305-323 (1992).

この発明は、音響インテンシティーの測定結果を利用して音源の位置を的確に検出することができる音源の検出方法及びその検出装置を提供するものである。

この発明（Ａ１）音源の検出方法は、
音場において受音点で直交する３軸方向（ｘ，ｙ，ｚ）の音響インテンシティーＩ_x，Ｉ_y，Ｉ_z［Ｗ／ｍ²］を測定し、それらの合計値であるＩ_all＝√（Ｉ_x ²＋Ｉ_y ²＋Ｉ_z ²）を算出する音響インテンシティーの測定工程と（√（）は括弧内の値の平方根を示す。以下同様。）、
前記測定工程で得た合計値の音響インテンシティー波形におけるピークが音源を表していると仮定して、そのピークを時間窓幅ごとに抽出し、その抽出されたピークから音源の大きさ及び距離を算出する音源の第１抽出工程と、
前記測定工程で得た合計値のピーク幅を推定して、そのピーク幅に対応する前記音響インテンシティーＩ_x，Ｉ_y，Ｉ_zにおけるピーク幅内の大きさの符号付の平均値を算出し、前記音響インテンシティーＩ_x，Ｉ_y，Ｉ_zにおけるピーク幅内の符号付の平均値と前記第１抽出工程で得た音源の距離とから音源の座標を推定する音源座標の第１推定工程と、
を備えるものである。

また、この発明（Ａ２）の音源の検出方法は、
音場において受音点で直交する３軸方向（ｘ，ｙ，ｚ）の音響インテンシティーＩ_x，Ｉ_y，Ｉ_z［Ｗ／ｍ²］を測定し、それらの合計値であるＩ_all＝√（Ｉ_x ²＋Ｉ_y ²＋Ｉ_z ²）を算出する音響インテンシティーの測定工程と、
前記測定工程で得た合計値の音響インテンシティーを音源インテンシティー成分と浮遊インテンシティー成分に区別する空間移動速度の閾値を設定し、前記合計値の音響インテンシティーから前記閾値以下の空間移動速度になる音源インテンシティー成分を切り出すか、又は、前記合計値の音響インテンシティーを音圧で除して得られる音響粒子速度の合計値から前記閾値以下の空間移動速度になる当該音響粒子速度の到来時間と対応する音源インテンシティー成分を前記合計値の音響インテンシティーから切り出し、その切り出された一片として含まれる音源インテンシティー成分の強さ及び到来時間の平均値を算出して音源の大きさ及び距離とする音源の第２抽出工程と、
前記第２抽出工程で切り出された音源インテンシティー成分の一片と同じ時間幅に含まれる音響インテンシティーＩ_x，Ｉ_y，Ｉ_zを切り出し、その切り出した音響インテンシティーＩ_x，Ｉ_y，Ｉ_zの符号付の大きさの平均値を算出し、前記音響インテンシティーＩ_x，Ｉ_y，Ｉ_zの符号付の平均値と前記第２抽出工程で得た音源の距離とから音源の座標を推定する音源座標の第２推定工程と、
を備えるものである。

一方、この発明（Ｂ１）の音源の検出装置は、
音場において受音点で直交する３軸方向（ｘ，ｙ，ｚ）の音響インテンシティーＩ_x，Ｉ_y，Ｉ_z［Ｗ／ｍ²］を測定し、それらの合計値であるＩ_all＝√（Ｉ_x ²＋Ｉ_y ²＋Ｉ_z ²）を算出する音響インテンシティーの測定手段と、
前記測定手段で得た合計値の音響インテンシティー波形におけるピークが音源を表していると仮定して、そのピークを時間窓幅ごとに抽出し、その抽出されたピークから音源の大きさ及び距離を算出する音源の第１抽出手段と、
前記測定手段で得た合計値のピーク幅を推定して、そのピーク幅に対応する前記音響インテンシティーＩ_x，Ｉ_y，Ｉ_zにおけるピーク幅内の大きさの符号付の平均値を算出し、前記音響インテンシティーＩ_x，Ｉ_y，Ｉ_zにおけるピーク幅内の符号付の平均値と前記第１抽出手段で得た音源の距離とから音源の座標を推定する音源座標の第１推定手段と、
を備えるものである。

また、この発明（Ｂ２）の音源の検出装置は、
音場において受音点で直交する３軸方向（ｘ，ｙ，ｚ）の音響インテンシティーＩ_x，Ｉ_y，Ｉ_z［Ｗ／ｍ²］を測定し、それらの合計値であるＩ_all＝√（Ｉ_x ²＋Ｉ_y ²＋Ｉ_z ²）を算出する音響インテンシティーの測定手段と、
前記測定手段で得た合計値の音響インテンシティーを音源インテンシティー成分と浮遊インテンシティー成分に区別する空間移動速度の閾値を設定し、前記合計値の音響インテ
ンシティーから前記閾値以下の空間移動速度になる音源インテンシティー成分を切り出すか、又は、前記合計値の音響インテンシティーを音圧で除して得られる音響粒子速度の合計値から前記閾値以下の空間移動速度になる当該音響粒子速度の到来時間と対応する音源インテンシティー成分を前記合計値の音響インテンシティーから切り出し、その切り出された一片として含まれる音源インテンシティー成分の強さ及び到来時間の平均値を算出して音源の大きさ及び距離とする音源の第２抽出手段と、
前記第２抽出手段で切り出された音源インテンシティー成分の一片と同じ時間幅に含まれる音響インテンシティーＩ_x，Ｉ_y，Ｉ_zを切り出し、その切り出した音響インテンシティーＩ_x，Ｉ_y，Ｉ_zの符号付の大きさの平均値を算出し、前記音響インテンシティーＩ_x，Ｉ_y，Ｉ_zの符号付の平均値と前記第２抽出手段で得た音源の距離とから音源の座標を推定する音源座標の第２推定手段と、
を備えるものである。

上記発明Ａ１及びＢ１の音源の検出方法及び検出装置によれば、第１抽出工程（手段）と第１推定工程（手段）を備えない検出方法や検出装置に比べて、音響インテンシティーの測定結果を利用して音源の位置を的確に検出することができる。

また、上記発明Ａ２及びＢ２の音源の検出方法及び検出装置によれば、第２抽出工程（手段）と第２推定工程（手段）を備えない検出方法や検出装置に比べて、音響インテンシティーの測定結果を利用して音源の存在も含めてその位置をより的確に検出することができる。

実施の形態１に係る音源の検出装置の概要を示す説明図である。 実施の形態１に係る音源の検出方法の処理工程を示すフローチャートである。 図２の検出方法等の概念とその検出方法の一過程で得られる結果の一例を示す説明図である。 音響インテンシティーを測定するマイクロホンの設置関係を示す説明図である。 図２の検出方法等における音場の構成例とその検出方法による結果の一例を示す説明図である。 実施の形態２に係る音源の検出装置の概要を示す説明図である。 実施の形態２に係る音源の検出方法の処理工程を示すフローチャートである。 音源インテンシティー成分と浮遊インテンシティー成分の概念及び性質を示す説明図である。 図７の検出方法等の概念とその検出方法の一過程で得られる結果の一例を示す説明図である。 図７の検出方法等における音場の構成例とその検出方法による結果の一例を示す説明図である。 図７の検出方法等における音場の構成例とその検出方法（音響粒子速度の空間移動速度を採用した場合）による結果の一例を示す説明図である。 長所（ｅ１）に関する比較試験の条件とその結果を示す説明図である。 長所（ｅ２）に関する比較試験の条件とその結果を示す説明図である。 参考の比較試験の結果を示す説明図である。 従来の音解析方法１の概念とその方法１で得られる結果の一例を示す説明図である。 従来の音解析方法２で得られる結果の一例を示す説明図である。 従来の音解析方法３で得られる結果の一例を示す説明図である。 音解析方法２、３の音源抽出原理の一部を示す説明図である。 ピーク検知を用いて音響インテンシティーから音源情報を検出する方法４の概念とその一過程で得られる結果の一例を示す説明図である。 図１９の検出方法４における音場の構成例とその検出方法４による結果の一例を示す説明図である。 音響インテンシティーの波形のピーク情報から音源情報を得られない場合の一例を示す説明図である。 図２の検出方法において音源の抽出が正しくできないおそれのある場合の一例を示す説明図である。

以下、この発明を実施するための形態（以下、単に「実施の形態」という）について添付の図面を参照しながら説明する。

［実施の形態１］
図１及び図２は、実施の形態１に係る音源の検出方法及びその検出装置の概要を示すものである。図１はその音源の検出装置の概要を示し、図２はその音源の検出方法の概要（工程）を示している。

実施の形態１に係る音源の検出装置１Ａは、実施の形態１に係る音源の検出方法を実施することができる装置である。この検出方法は、前述した「ピーク検知」を用いて音響インテンシティーから音源情報を検出する方法４に対して、音響インテンシティーの波形におけるピークの大きさではなく所定のピーク幅内の平均値で処理して音源を抽出したうえで音源位置の推定を行うよう改良した点に特徴がある。この観点からすると、この検出方法は、言わば「改良型ピーク検知」を用いて音響インテンシティーから音源情報を検出する方法と称することができる。

そして、検出装置１Ａは、図１に示されるように、音響インテンシティーの測定手段２、音源の第１抽出手段３Ａ、及び音源位置の第１推定手段４Ａを主に備えている。また、検出装置１Ａは、必要に応じて、各手段２，３Ａ，４Ａにおける結果などの情報を表示するための液晶パネル等の表示手段（記憶手段も含む）５を備えることもある。
この検出装置１Ａにおける音響インテンシティーの測定手段２、音源の第１抽出手段３Ａ、及び音源位置の第１推定手段４Ａについては、例えば、主として検出方法の各工程を処理するためのプログラムや必要なデータからなるソフトウェアによって構成することができる。また、そのソフトウェアのプログラムやデータをパーソナルコンピュータ等の機器に格納・搭載して作動させるよう構築することにより送出装置１Ａが製作される。この点は、後述する実施の形態２に係る検出装置１Ｂについても同様である。

（音響インテンシティーの測定手段及び測定工程）
音響インテンシティーの測定手段２は、図１に示されるように、音場において受音点で直交する３軸方向（ｘ，ｙ，ｚ）の音響インテンシティーＩ_x，Ｉ_y，Ｉ_z［Ｗ／ｍ²］を測定する測定部２１と、その音響インテンシティーＩ_x，Ｉ_y，Ｉ_zの合計値であるＩ_all＝√（Ｉ_x ²＋Ｉ_y ²＋Ｉ_z ²）を算出する算出部２２を有している。
また、音源の検出方法では、図２に示されるように、この測定部２１の処理と算出部２２の処理とをステップ１０（以下、このような語句を単に「Ｓ１０」として表記することもある。）とステップ１１（Ｓ１１）としてこの順で行い、このステップ１０、１１で示す２つの処理工程を併せて音響インテンシティーの測定工程としている。

測定部２１は、音場における音源の情報を計測する測定機器から音響インテンシティーの測定に必要な情報を入力して音響インテンシティーＩ_x，Ｉ_y，Ｉ_zを算出する処理を行う部分である。測定機器としては、図１に示されるように、例えば、マイクロホン６、音響粒子速度センサー７、データ取得機器８等で主に構成される測定機器が適用される。
このうち、マイクロホン６のみを用いて音響インテンシティーを測定する場合は、図３の左側部分に示されるように、推定対象となる音源（実音源Ｓ₁、虚音源Ｓ₂などの複数の音源）が存在する音場（室内等の閉鎖空間など）においてマイクロホン６により音圧を測定する。次いで、マイクロホン６で測定された音圧の波形情報から、その測定時に設定される受音点Ｐで直交する３軸方向（ｘ，ｙ，ｚの方向）の音響インテンシティーＩ_x，Ｉ_y，Ｉ_zを入手するために必要な音情報（音圧と音響粒子速度）を得る。
また、マイクロホン６と音響粒子速度センサー７を用いて音響インテンシティーを測定する場合は、上記音情報の１つである音響粒子速度をマイクロホンの出力波形から算出するのではなく音響粒子速度センサー７により直接測定して得る。この場合、音響粒子速度センサー７は、上記３軸方向の軸上にそれぞれ設置する。
データ取得機器８は、マイクロホン６や音響粒子速度センサー７により取得した音源の情報を増幅して蓄積（記憶）する、マイクアンプ、記憶装置等で構成される機器である。

マイクロホン６としては、音響インテンシティーの測定法に応じた種類や数のものが使用される。例えば、全指向性マイクロホンを使用するＰ−Ｐ法の場合にはそのマイクロホンを４つ又は６つ使用し、単一指向性マイクロホンを使用するＣ−Ｃ法の場合にはそのマイクロホンを６つ使用する。
測定部２１では、音響インテンシティーの測定に際して、上記いずれかの種類のマイクロホン６を単独で使用するか、又は、マイクロホン６と音響粒子速度センサー７とを組み合わせて使用することができる。例えば、その測定に際して全指向性マイクロホンを４つ用いるＰ−Ｐ法を適用する場合は、４つの小型の全指向性マイク６ａ，６ｂ，６ｃ，６ｄを、図４に例示するように、軸方向の原点（Ｏ）とする位置とその原点から所定の距離（間隔）ｋずつ離れた直交３軸（ｘ軸、ｙ軸、ｚ軸）上の３つの位置とに配分して配置されるよう、図示しない治具により支持した構成にしたうえで使用する。図４に示す例では、マイク６ａを原点の位置、マイク６ｂをｘ軸上、マイク６ｃをｙ軸上、マイク６ｄをｚ軸上に配置している。また、このときの間隔ｋは、例えば、１０ｍｍ〜３０ｍｍの範囲内で設定される。

測定は、例えば全指向性マイクロホンを４つ用いるＰ−Ｐ法を適用する場合、以下のように行われる。
はじめに、上記構成の全指向性マイク６ａ，６ｂ，６ｃ，６ｄを測定対象の音場に設置する。この際、マイク６ａが測定時に音場に設定される直交３軸の原点Ｏの位置に、マイク６ｃが実音源Ｓ₁と向き合う正面の方向になるｙ軸上の距離ｋの位置（ｘ＝０，ｙ＝ｋ，ｚ＝０）に、マイク６ｂが音源Ｓ₁に対して水平右側の方向になるｘ軸上の距離ｋの位置（ｘ＝ｋ，ｙ＝０，ｚ＝０）に、マイク６ｄが鉛直上方の方向になるｚ軸上の距離ｋの位置（ｘ＝０，ｙ＝０，ｚ＝ｋ）に、それぞれ配置された状態で設置される。距離ｋについては、例えば１５ｍｍに設定した。
次いで、音場における実音源Ｓ₁から測定信号を再生して、上記状態で設置された全指向性マイク６ａ，６ｂ，６ｃ，６ｄにより測定信号を録音して得られる音圧波形又はその音圧波形から算出されたインパルス応答をデータ取得機器８に保存する。
実施の形態１では、音場として、図５の左側の部分に示されるように、反射音が非常に小さく実音源Ｓ₁のみが音源として検知されるように壁面、天井面及び床面に吸音材を敷設した部屋（閉鎖空間）を採用している。この部屋は、ｘ，ｙ，ｚ方向の各寸法が４ｍ，３．３ｍ，４ｍである大きさのものである。
実音源Ｓ₁としては、スピーカーを使用した。また、そのスピーカーは、直交３軸におけるｙ軸上の２ｍの位置（ｘ＝０、ｙ＝２ｍ、ｚ＝０の座標位置）に、スピーカーの正面部分が原点Ｏに向く状態で設置した。直交３軸の原点（Ｏ）は、閉鎖空間の床面のほぼ中心の位置から１．１ｍの高さに設定した。

そして、測定部２１では、次の手順により音響インテンシティーＩ_x，Ｉ_y，Ｉ_zが得られる。
まず、音情報の音圧については、上記全指向性マイクロホンを４つ用いるＰ−Ｐ法を適用した場合、その４つのマイク６ａ，６ｂ，６ｃ，６ｄにより測定した時刻ｔの音圧をｐ_a(t)，ｐ_b(t)，ｐ_c(t)，ｐ_d(t) としたとき、以下に示すｘ，ｙ，ｚ方向の時刻ｔにおける音圧の平均ｐ_x(t) ，ｐ_y(t) ，ｐ_z(t)として算出される。
ｐ_x(t)=（ｐ_a(t)＋ｐ_b(t)）／２
ｐ_y(t)=（ｐ_a(t)＋ｐ_c(t)）／２
ｐ_z(t)=（ｐ_a(t)＋ｐ_d(t)）／２
次に、音情報の音響粒子速度については、以下の数式で表されるｘ，ｙ，ｚ方向の時刻ｔにおける音響音響粒子速度Ｕ_x(t)，Ｕ_y(t)，Ｕ_z(t)として算出される。
最後に、音響インテンシティーは「音響インテンシティー＝音圧×音響粒子速度」として表されることから、上記音圧及び音響粒子速度の測定結果を用いて、以下の数式で表されるｘ，ｙ，ｚ方向の時刻ｔにおける音響インテンシティーＩ_x(t)，Ｉ_y(t)，Ｉ_z(t)として求められる。
Ｉ_x(t) ＝ ｐ_x(t)×ｕ_x(t)，Ｉ_y(t) ＝ ｐ_y(t)×ｕ_y(t)，Ｉ_z(t) ＝ ｐ_z(t)×ｕ_z(t)

算出部２２は、測定部２１で得た音響インテンシティーＩ_x，Ｉ_y，Ｉ_zのデータを用いてそれらの合計値であるＩ_allを下記の数式により算出する処理を行う部分である。この合計値Ｉ_allの単位も「Ｗ／ｍ²」である。
Ｉ_all＝√（Ｉ_x ²＋Ｉ_y ²＋Ｉ_z ²）

この測定手段２における測定部２１及び算出部２２で得られる音響インテンシティーＩ_x，Ｉ_y，Ｉ_zと合計値Ｉ_allの結果については、例えば、表示手段５において、図３に右側の部分に示されるような波形として表示することができる。図３等に示される横軸ｔは測定時の測定信号を再生してからの経過時間を示し、その横軸ｔと直交する上下方向（縦軸）は波形の振幅を示す。また、図３の右側部分に示す波形は、コンサートホールで測定したときに得られた音響インテンシティーの波形の一例である。

（音源の第１抽出手段及び第１抽出工程）
第１抽出手段３Ａは、図１に示されるように、測定手段２（工程）で得た合計値Ｉ_allの音響インテンシティー波形におけるピークを時間窓幅ごとに抽出する抽出部３１と、その抽出されたピークから音源の大きさ及び距離を算出する算出部３２を有している。
また、音源の検出方法では、図２に示されるように、この抽出部３１の処理と算出部３２の処理とをステップ１２、１３（Ｓ１２、Ｓ１３）としてこの順で行い、このステップ１２、１３で示す２つの処理工程を併せて音源の第１抽出工程としている。

抽出部３１は、合計値Ｉ_allの音響インテンシティー波形におけるピークが音源を表していると仮定したうえで、図３に示すように、そのピークを時間窓幅Δｔごとに抽出する処理を行う部分である。
時間窓幅Δｔは、連続した長い時間軸の波形から或る一定のデータ列を切り取り、そのデータ列の中からピークを見つける際の切り取るときの時間の長さに相当するものある。これにより、ピークは、時間窓幅Δｔの間隔ごとに１つずつ抽出されることになる。
また、時間窓幅Δｔは、解析帯域幅Δｆ（Ｈｚ）に反比例する値である（Δｔ∝Δｆ^-1）。例えば、解析帯域幅Δｆが４０００[Ｈｚ]（１００[Ｈｚ]〜４１００[Ｈｚ]）のときは時間窓幅Δｔが１．３[ｍｓ]となるのに対し、解析帯域幅Δｆが２０００[Ｈｚ]（１００[Ｈｚ]〜２１００[Ｈｚ]）と狭くなるときはその時間窓幅Δｔが３．５[ｍｓ]となり、長い値になる。
図３では、時間窓幅Δｔの間隔ごとに抽出した各ピーク（図中に「・」で示している部分）にＳ₁，Ｓ₂，Ｓ₃，Ｓ₄，・・・の記号をそれぞれ付しているが、この記号は推定対象の音場における複数の音源Ｓ₁，Ｓ₂，Ｓ₃，Ｓ₄，・・・に対するピークであることを示している。

算出部３２は、その抽出されたピークから音源の大きさ及び距離を算出する処理を行う部分である。
ここで、音源の大きさとは、実音源Ｓ₁や反射音（虚音源Ｓ₂，Ｓ₃，…）の強さ、すなわち実音源Ｓ₁や虚音源Ｓ₂，Ｓ₃，…から到来する音の音量を表すものである。また、音源の距離は、実音源Ｓ₁や虚音源Ｓ₂，Ｓ₃，…から受音点Ｐまでの距離を表すものである。この音源の距離は、音速で除すれば実音源や反射音の到来時間を表す値となる。
そして、音源の大きさ及び距離は次のようにして算出される。
音源の大きさについては、各ピーク位置における合計値Ｉ_allの大きさ（図３の縦軸の値）（［Ｗ/ｍ²］）として算出する。また、音源の距離については、各ピーク位置の時間［ｓ］（図３の横軸の値）に音速［ｍ/ｓ］を乗じた値（［ｍ］）として算出する。

（音源位置の第１推定手段及び第１推定工程）
音源位置の第１推定手段４Ａは、図１に示されるように、測定手段２及びその測定工程で得た合計値Ｉ_allのピーク幅を推定する推定部４１と、そのピーク幅に対応する音響インテンシティーＩ_x，Ｉ_y，Ｉ_zにおけるピーク幅内の符号付の大きさの平均値を算出する算出部４２と、音響インテンシティーＩ_x，Ｉ_y，Ｉ_zにおけるピーク幅内の符号付の大きさの平均値と第１抽出手段３Ａ及び工程で得た音源の距離とから音源の座標を推定
する推定部４３とを有している。
また、音源の検出方法では、図２に示されるように、この推定部４１の処理と算出部４２の処理と推定部４３の処理とをステップ１４〜１６（Ｓ１４〜Ｓ１６）としてこの順で行い、このステップ１４〜１６で示す３つの処理工程を併せて音源位置の第１推定工程としている。

推定部４１は、音響インテンシティーの合計値Ｉ_allのピーク幅を推定する処理を行う部分である。
このピーク幅の推定は、例えば、合計値の音響インテンシティー波形において、時間窓幅Δｔの間隔ごとに抽出される各ピークの前後で盛り上がる波形部分の始点と終点の間をピーク幅として推定する。具体的には、上記音響インテンシティー波形において、時間窓幅Δｔの間隔ごとの各ピークから時間軸上を前後に移動させたときに波形の変化量（傾きなど）が再度緩やかになる前後２箇所の時刻を解析し、その前後２箇所の時刻で挟まれる波形部分をピーク幅として得ることができる。
この推定部４１の推定処理により、図３の右側部分の上部や図５の中央部に示されるように、一点鎖線で挟まれる部分が、時間窓幅Δｔの間隔ごとの各ピークに対するピーク幅としてそれぞれ求められる。

算出部４２は、推定部４１で推定した各ピーク幅に対応する音響インテンシティーＩ_x，Ｉ_y，Ｉ_zにおける各ピーク幅内の大きさの符号付の平均値を算出する処理を行う部分である。
各ピーク幅に対応する音響インテンシティーＩ_x，Ｉ_y，Ｉ_zにおける各ピーク幅は、図３の右側部分の下部側に示されるように、音響インテンシティーＩ_x，Ｉ_y，Ｉ_zの各波形のうち合計値Ｉ_allの音響インテンシティーの波形で推定された２本の破線で挟まれる部分と同じ範囲（時刻間）で挟まれる波形部分になる。
また、各ピーク幅内の音響インテンシティーＩ_x，Ｉ_y，Ｉ_zの符号付の大きさの平均値（±｜Ｉ_{x_ave}｜，±｜Ｉ_{y_ave}｜，±｜Ｉ_{z_ave}｜）は、例えば、ピーク幅の開始時間をｔ₁、その終了時間をｔ₂とし、時刻ｔのときのＩ_x，Ｉ_y，Ｉ_zをそれぞれＩ_x（ｔ），Ｉ_y（ｔ），Ｉ_z（ｔ）と表記すると、以下の数式に基づいてそれぞれ算出される。
この算出部４２の算出処理により、図３の右側部分の下部側に示されるように、Ｉ_x，Ｉ_y，Ｉ_zの各音響インテンシティー波形の各ピーク幅における符号付の平均値（図中に「○」で示される部分）がそれぞれ求められる。

推定部４３は、算出部４２で得た音響インテンシティーＩ_x，Ｉ_y，Ｉ_zにおけるピーク幅内の符号付の平均値と第１抽出手段３Ａ及びその工程で得た音源の距離とから音源の座標を推定する処理を行う部分である。
音源の座標は、例えば、あるピークに対応するピーク幅内のｘ，y，ｚ方向の音響イ
ンテンシティー波形の符号付の平均値（±｜Ｉ_{x_ave}｜，±｜Ｉ_{y_ave}｜，±｜Ｉ_{z_ave}｜）にて表される原点からの方位（±｜Ｉ_{x_ave}｜，±｜Ｉ_{y_ave}｜，±｜Ｉ_{z_ave}｜）に対して、第１推定工程で算出した合計値Ｉ_allのピークが与える距離ほど離れた地点として推定される。このときの合計値Ｉ_allのピークが与える距離ほど離れた地点とは、例えば、音源の第１抽出工程等で算出した合計値Ｉ_allのピークが与える距離がＬ［ｍ］であり、Ｉ_x，Ｉ_y，Ｉ_zから算出された方位が（＋２，−３，＋４）であるとした場合、直交３軸の原点Ｏから座標（＋２，−３，＋４）の方向に、原点Ｏから距離（長さ）Ｌ［ｍ］に相当する分だけ離れた位置を意味する。
この算出部４２の算出処理により、図５の右側部分に一部の結果が示されるように、推定対象の実音源Ｓ₁について推定した座標（図中の黒丸）は、実音源Ｓ₁の正規の位置（白丸）であるｙ軸上の２ｍの位置（ｘ＝０、ｚ＝０）とほぼ一致していることがわかる。

（音源位置の検出結果）
図５の右側部分に示されている一部の結果から、実施の形態１に係る音源の検出方法（その検出装置１Ａ）、すなわち別称の「改良型ピーク検知」を用いて音響インテンシティーから音源情報を検出する方法によれば、次のことがわかる。
すなわち、この検出方法では、ｘ方向及びｚ方向の音響インテンシティーＩ_x，Ｉ_zに測定誤差による残差波形が存在するものの（図５の中央部分）、前述した「ピーク検知」を用いて音響インテンシティーから音源情報を検出する方法４とは異なり、ピーク幅の大きさの平均値を適用しているため、ｘ方向及びｚ方向の音響インテンシティーＩ_x，Ｉ_zがいずれも、その残差波形における特定の値に支配されることがなく、ほぼゼロの値（座標）として認識できていることがわかる。

この結果、このときに検出（推定）した音源Ｓ₁の位置（黒丸）は、図５の右側部分に一部の結果として示されているとおり、ｙ軸上にあるはずの音源Ｓ₁の正規の位置（白丸）とほぼ一致した場所に推定されるようになる。つまり、この検出方法によれば、上記検出する方法４の場合（図２０参照）に比べて、音源Ｓ1の位置が精度よくｙ軸上に推定されており、音源の座標の推定精度が改善されていることがわかる。
したがって、この検出方法及びその検出装置１Ａでは、音源の位置を的確に検出することができる。

［実施の形態２］
図６及び図７は、実施の形態２に係る音源の検出方法及びその検出装置の概要を示すものである。図６はその音源の検出装置の概要を示し、図７はその音源の検出方法の概要（工程）を示している。

実施の形態２に係る音源の検出装置１Ｂは、実施の形態２に係る音源の検出方法を実施することができる装置である。この検出方法は、前述した「ピーク検知」を用いて音響インテンシティーから音源情報を検出する方法４や実施の形態１に係る検出方法とは異なり、後述する理由からピーク検知をまったく用いることなく、音響インテンシティー又は音響粒子速度の空間移動速度に着目して音源を抽出したうえで音源位置の推定を行うよう工夫した点に特徴がある。この観点からすると、この検出方法は、言わば「速度検知」を用いて音響インテンシティーから音源情報を検出する方法と称することができる。

そして、検出装置１Ｂは、図６に示されるように、音響インテンシティーの測定手段２、音源の第２抽出手段３Ｂ、及び音源位置の第２推定手段４Ｂを主に備えている。また、検出装置１Ｂは、実施の形態１に係る検出装置１Ａの場合と同様に、必要に応じて、各手段２，３Ｂ，４Ｂにおける結果などの情報を表示するための液晶パネル等の表示手段（記憶手段も含む）５を備えることもある。

ここで、実施の形態２に係る音源の検出方法及びその検出装置１Ｂが完成するに至るまでの背景について、以下に説明する。

まず着想のきっかけは、前述した実施の形態１に係る検出方法及びその検出装置１Ａでは、例えば、音響インテンシティーの波形が激しく変化する測定結果や逆にほとんど変化しない測定結果に対して音源自体の存在を正しく認識できない可能性があることを発見したことにある。
この認識できない可能性がある理由は、実施の形態１に係る検出方法及びその検出装置１Ａでは音源の存在を判別するための音源の第１抽出工程（第１抽出手段３Ａ）において音響インテンシティーの波形のピークを検知する手法（改良型ピーク検知）を採用している点に影響しているものと考えられる。

具体的には、例えば、図２２に例示するように、直交３軸方向における原点（Ｏ）上の音源Ｓ₀やその原点近傍の音源Ｓ₁などが存在すると、原理上、その音響インテンシティー波形の振幅及びその変化がゼロもしくは極小になる（同図の左側下部及び右側下部に例示する結果を参照。）。また、その音源Ｓ₀、音源Ｓ₁に対する実際の測定を行った場合には、その音響インテンシティーの波形が多くの誤差を含む暴れた波形になることが多い。つまり、原点やその近傍の位置に存在する音源Ｓ0、Ｓ₁については、その音響インテンシティー波形のピークそのものを正確に抽出することが難しくなる。
したがって、このような原点やその近傍の位置にある音源（Ｓ₀，Ｓ₁）に対しては、前述した「改良型ピーク検知」の検出方法を用いても、その音源の抽出が正しくできない可能性がある。

そこで、本発明者は、このような実状に鑑みて鋭意研究した結果、ピークを検知する手法を用いない音源の抽出方法として、音響インテンシティー又は音響粒子速度の空間移動速度に着目して音源を抽出する手法が有効であることを見い出し、最終的に「速度検知」を用いて音響インテンシティーから音源の位置を検出する方法及びその検出装置を完成させるに至った。

そして、この検出方法及びその検出装置１Ｂは、以下の仮定を前提にすることで成立している。
すなわち、音響インテンシティーには、次の２種類の成分が存在すると仮定する。
１．音源を示しているもの。以下、これを「音源インテンシティー成分」と呼ぶこととする。
２．音源間をさまよっているもの。以下、これを「浮遊インテンシティー成分」と呼ぶこととする。
また、上記２種類のインテンシティー成分は、次の性質を有すると仮定する。
３．音源インテンシティー成分は、空間移動速度が遅い。
４．浮遊インテンシティー成分は、空間移動速度が速い。
図８は、この仮定した内容を参考までに図案化して示している。

したがって、この検出方法及びその検出装置１Ｂにおいては、以上の仮定のもとに、音響インテンシティーの波形から速度検知により音源を抽出するように構成している。
以下、この検出方法及びその検出装置１Ｂの説明を続ける。

（音響インテンシティーの測定手段及び測定工程）
音響インテンシティーの測定手段２は、図６に示されるように、音場において受音点で直交する３軸方向（ｘ，ｙ，ｚ）の音響インテンシティーＩ_x，Ｉ_y，Ｉ_z［Ｗ／ｍ²］を測定する測定部２１と、その音響インテンシティーＩ_x，Ｉ_y，Ｉ_zの合計値であるＩ_all＝√（Ｉ_x ²＋Ｉ_y ²＋Ｉ_z ²）を算出する算出部２２を有している。
この測定部２１と算出部２２は、実施の形態１に係る検出装置１Ａにおける測定部２１及び算出部２２と同じ構成からなるものである。

また、音源位置の検出方法では、図７に示されるように、この測定部２１の処理と算出部２２の処理とをステップ２０とステップ２１としてこの順で行い、このステップ２０、２１で示す２つの処理工程を併せて音響インテンシティーの測定工程としている。
この音響インテンシティーの測定工程についても、実施の形態１に係る検出方法における音響インテンシティーの測定工程と同じ構成になっている。

この測定部２１による測定の概念を図９の左側の部分に例示する。
また、この測定手段２における測定部２１及び算出部２２で得られる音響インテンシティーＩ_x，Ｉ_y，Ｉ_zと合計値Ｉ_allの結果については、実施の形態１の場合（図３）と同様に、例えば、表示手段５において、図９の右側の部分に示されるような波形として表示することができる。また、図９の右側部分に示す波形は、コンサートホールで測定したときに得られた音響インテンシティーの波形の一例である。
さらに、この測定手段２及び測定工程では、実施の形態１の測定手段及び測定工程に関して説明したとおり、音響インテンシティーを得るために必要な音情報の１つである音響粒子速度Ｕ_x，Ｕ_y，Ｕ_z［ｍ／ｓ］が得られるため、その合計値であるＵ_all＝√（Ｕ_x ²＋Ｕ_y ²＋Ｕ_z ²）［ｍ／ｓ］も同時に算出することができる。

（音源の第２抽出手段及び第２抽出工程）
第２抽出手段３Ｂは、図６に示されるように、測定手段２及びその測定工程で得た合計値Ｉ_allの音響インテンシティーを前記仮定の音源インテンシティー成分と浮遊インテンシティー成分に区別する空間移動速度の閾値Ｖｓを設定する設定部３５と、合計値Ｉ_allの音響インテンシティーから閾値Ｖｓ以下の空間移動速度になる音源インテンシティー成分を切り出すか、又は、合計値Ｕ_allの音響粒子速度から閾値Ｖｓ以下の空間移動速度になる当該音響粒子速度の到来時間と対応する音源インテンシティー成分を合計値Ｉ_allの音響インテンシティーから切り出す切り出し部３６と、その切り出された一片として含まれる音源インテンシティー成分の強さ及び到来時間の平均値を算出して音源の大きさ及び距離とする算出部３７とを有している。
また、音源の検出方法では、図７に示されるように、この設定部３５の処理と切り出し部３６の処理と算出部３７の処理とをステップ２２〜２４（Ｓ２２〜Ｓ２４）としてこの順で行い、このステップ２２〜２４で示す３つの処理工程を併せて音源の第２抽出工程としている。

設定部３５は、合計値Ｉ_allの音響インテンシティーから特に音源インテンシティー成分を切り出すために適切な閾値Ｖｓを設定する処理を行う部分である。
このような閾値Ｖｓは、厳密には測定環境（部屋等の音場の大きさ、音源の距離など）や解析条件によって変化するが、一般的に音の到来時間ｔ［秒］に比例した値となる。
また、この閾値Ｖｓについては、例えば到来時間ｔ［秒］の関数で与えられる以下の数式で求めることができる。
Ｖｓ＝（ｔ×１０００＋１０）^1.5×４０ ［ｍ／秒］
上記式中の「１０００」は 到来時間ｔの単位［秒］から［ミリ秒］への変換を示し、
「１．５」は到来時間ｔに比例する次元を示し、「４０」は到来時間ｔを空間移動速度の閾値Ｖｓに変換する際の比率を示す。
実施の形態２では、受音点Ｐ（原点Ｏ）から実音源Ｓ₁までの距離を２［ｍ］にしているため、実音源Ｓ₁からの音の到来時間ｔは、ｔ＝２／３４４＝０．００５８１４［秒］となる。「３４４」は、気温２０℃のときの音速［ｍ／ｓ］を示す。
したがって、このときの閾値Ｖｓは、上記数式からＶｓ＝（０．００５８１４×１０００＋１０）^1.5×４０＝２，５１６［ｍ／秒］となる。

切り出し部３６は、合計値Ｉ_allの音響インテンシティーから閾値Ｖｓ以下になる音源インテンシティー成分を切り出す処理を行う部分である。
この音源インテンシティー成分の切り出しは、以下のように行われる。
はじめに、合計値Ｉ_allの音響インテンシティー又は音響粒子速度の合計値Ｕ_allから閾値Ｖｓ以下の空間移動速度になる音響インテンシティー又は音響粒子速度の到来時間、すなわち音源インテンシティー成分の到来時間を抽出する。
この音源インテンシティー成分の到来時間の抽出に際しては、まず、音響インテンシティーの合計値Ｉ_allにおける或る時刻の座標から次の時刻の座標への移動距離、又は音響
粒子速度の合計値Ｕ_allにおける或る時刻の座標から次の時刻の座標への移動距離を算出する。例えば、サンプリング周波数が４８，０００[Ｈｚ]の条件で測定を実施した場合、ある時刻から次の時刻への最小の時間間隔は１／４８，０００[秒]となる。
続いて、上記移動距離を移動時間で除して音響インテンシティー又は音響粒子速度の空間移動速度を算出する。この算出で得られる空間移動速度が閾値Ｖｓ以下の小さい値になった場合、そのときの時刻を「ｔ_s」とする。また、上記音響インテンシティー又は音響粒子速度の空間移動速度の算出を同様に続ける。この続けた算出で得られる空間移動速度が閾値Ｖｓよりも大きい値になった場合、そのときの時刻を「ｔ_e」とする。このときの時刻ｔ_sから時刻ｔ_eまでの時間が上記抽出すべき音源インテンシティー成分の到来時間になる。
最後に、この抽出された上記時刻ｔ_sから上記時刻ｔ_eまでの到来時間に含まれる音響インテンシティーの合計値Ｉ_allを切り出す。
ここで、この切り出された合計値Ｉ_all群を「合計値Ｉ_allの音響インテンシティーの一片」と呼ぶこととする。また、このときの合計値Ｉ_allの音響インテンシティーの一片は、全て音源インテンシティー成分で構成されており、或る一つの音源を表していることになる。

この切り出し部３６で切り出した音源インテンシティー成分の結果について、図９の右側の上部に例示する。
図９においては、合計値Ｉ_allの音響インテンシティーの波形のうち２本の破線で挟まれる波形部分（実線の両方向矢印で示す範囲：時間幅）が、切り出された音源インテンシティー成分に該当する。またこの場合、図９においては、合計値Ｉ_allの音響インテンシティーの波形のうち音源インテンシティー成分どうし間における波形部分（破線の両方向矢印で示す範囲：時間幅）が、閾値Ｖｓよりも大きい空間移動速度になる浮遊インテンシティー成分に該当することになる。
そして、この切り出したＩ_allの音響インテンシティーの一片（図９における１つの実線両方向矢印で示される波形部分、換言すれば音源インテンシティー成分の部分）が、１つの音源を表しているものと考えられる。
図９においては、切り出された音響インテンシティーの一片（音源インテンシティー成分）のそれぞれに対して大きさ及び時間の平均値を算出した値を「〇」の白丸で示している。そして、これらの白丸が、抽出した音源Ｓ₁，Ｓ₂，Ｓ₃，Ｓ₄，・・・となる。

算出部３７は、切り出し部３６で切り出した合計値Ｉ_allの音響インテンシティーの一片として含まれる音源インテンシティー成分の強さ及び到来時間の平均値を算出する処理を行う部分である。
ここで、上記強さとは、図９に示されるように音源インテンシティー成分の波形部分の振幅（図中の縦軸の値）をいう。また、上記到来時間とは、図９に示されるように音源インテンシティー成分の情報がある音源から受音点に到来するまでの所要時間（図中の横軸の値）をいう。
このときの算出は、以下に示すように行われる。
すなわち、音響インテンシティーの合計値Ｉ_allの切り出し開始時間をｔ₁，その切り出し終了時間をｔ₂、時刻ｔのときの合計値Ｉ_allをＩ_all（ｔ）と表記すると、強さの平均値Ｉ_{all_ave}と到来時間の平均値ｔ_aveは、次の数式に基づいてそれぞれ算出される。

算出部３７では、このようにして算出された音源インテンシティー成分の強さ及び到来時間の平均値を、最終的に各音源（Ｓ1，Ｓ2，Ｓ3，Ｓ4，・・・）の大きさ及び距離とす
る。ここでいう音源の大きさ及び距離とは、実施の形態１における音源の第１抽出手段３Ａ及びその抽出工程で算出されるときの音源の大きさ及び距離と同じ概念のものである。
図９には、各音源（Ｓ1，Ｓ2，Ｓ3，Ｓ4，・・・）の大きさ（強度の平均値）を「○」の白丸で示している。

（音源位置の第２推定手段及び第２推定工程）
音源位置の第２推定手段４Ｂは、図６に示されるように、第２抽出手段３Ｂ及びその抽出工程で切り出された音源インテンシティー成分の一片と同じ時間幅に含まれる音響インテンシティーＩ_x，Ｉ_y，Ｉ_zを切り出す切り出し部４５と、その切り出した音響インテンシティーＩ_x，Ｉ_y，Ｉ_zの大きさの平均値を算出する算出部４６と、音響インテンシティーＩ_x，Ｉ_y，Ｉ_zの大きさの平均値と第２抽出手段３Ｂ及びその抽出工程で得た音源の距離とから音源の座標を推定する推定部４７とを有している。
また、音源の検出方法では、図７に示されるように、この切り出し部４５の処理と算出部４６の処理と推定部４７の処理とをステップ２５〜２７（Ｓ２５〜Ｓ２７）としてこの順で行い、このステップ２５〜２７で示す３つの処理工程を併せて音源位置の第２推定工程としている。

切り出し部４５は、切り出された音源インテンシティー成分の一片と同じ時間幅に含まれる音響インテンシティーＩ_x，Ｉ_y，Ｉ_zを切り出す処理を行う部分である。このときの音源インテンシティー成分の一片と同じ時間幅は、図９の右側部分に示される各２本の破線で挟まれる部分（実線の両方向矢印で示す範囲）の時刻幅となる。また、このときの音響インテンシティーＩ_x，Ｉ_y，Ｉ_zの切り出しは、Ｉ_x，Ｉ_y，Ｉ_zの各音響インテンシティーの波形のうち上記各時間幅に挟まれる波形部分の情報を抽出することになる。

算出部４６は、上記切り出し部４５で切り出した音響インテンシティーＩ_x，Ｉ_y，Ｉ_zの符号付の大きさの平均値を算出する処理を行う部分である。
このときの大きさの平均値は、以下に示すように算出される。
すなわち、音響インテンシティーの合計値Ｉ_allの切り出し開始時間をｔ₁、切り出し終了時間をｔ₂、時刻ｔのときのＩ_x，Ｉ_y，Ｉ_zの値をＩ_x(t)，Ｉ_y(t)，Ｉ_z(t)と表記すると、時刻ｔのときのＩ_x，Ｉ_y，Ｉ_zの符号付の大きさの平均値（±｜Ｉ_x__ave｜，±｜Ｉ_{y_ave}｜，±｜Ｉ_{z_ave}｜）は、次の数式に基づいてそれぞれ算出される。
この算出部４６の算出処理により、図９の右側部分の下部に示されるように、上記時間幅におけるＩ_x，Ｉ_y，Ｉ_zの各音響インテンシティー（各音源インテンシティー）の大きさの平均値（図中に「○」の白丸で示される部分）がそれぞれ求められる。図９においてＩ_x，Ｉ_y，Ｉ_zの各音響インテンシティーの波形に付記された＋ｘ₁、＋ｙ₁、＋ｚ₁は、音源Ｓ₁に対応するＩ_x，Ｉ_y，Ｉ_zの符号付の平均値として、それぞれ＋符号の値が得られたことを示す。

推定部４７は、算出部４６で得た音響インテンシティーＩ_x，Ｉ_y，Ｉ_zの大きさの平均値と第２抽出手段３Ｂ及びその抽出工程で得た音源の距離とから音源の座標を推定する処理を行う部分である。
音源の座標の推定は、例えば、音源Ｓ₁に対応するx，y，z方向のＩ_x，Ｉ_y，Ｉ_zの符号付の平均値として、それぞれ＋ｘ₁、＋ｙ₁、＋ｚ₁が得られている場合、原点からの方位（＋ｘ₁、＋ｙ₁、＋ｚ₁）に対して、算出部３７で算出した音源Ｓ₁の到来時間の平均t_aveに音速を乗した値、すなわち音源距離ほど離れた地点として推定される。
この推定部４７の推定処理により、図１０の右側部分に一部の結果が示されるように、
推定対象の実音源Ｓ₁の推定した座標（図中の黒丸）は、実音源Ｓ₁の正規の位置（白丸）であるｙ軸上の２ｍの位置（ｘ＝０、ｚ＝０）とほぼ一致していることがわかる。

（音源位置の検出結果１）
図１０の右側部分に示されている一部の結果から、実施の形態２に係る音源の検出方法（その検出装置１Ｂ）、すなわち別称の「速度検知」を用いて音響インテンシティーから音源情報を検出する方法によれば、次のことがわかる。
すなわち、この検出方法では、実施の形態１に係る音源の検出方法等における検出結果（図５）と比較すると、ｘ方向における位置の推定に若干の誤差が発生しているものの、ｘ方向及びｚ方向の音響インテンシティーＩ_x，Ｉ_zの音源インテンシティーに該当する波形部分の大きさの値（平均値）がいずれもほぼゼロの値（座標）として認識できていることがわかる。

この結果、このときに検出（推定）した音源Ｓ₁の位置（黒丸）は、図１０の右側部分に一部の結果として示されるとおり、ｙ軸上にあるはずの音源Ｓ₁の正規の位置（白丸）とほぼ一致した場所に推定されるようになる。つまり、この検出方法によっても、実施の形態１に係る音源の検出方法等の場合と同様に、音源Ｓ₁の位置が精度よくｙ軸上に推定されていることがわかる。
したがって、この検出方法及びその検出装置１Ｂによっても、音源の位置を的確に検出することができる。

（音源位置の検出結果２）
また参考までに、音源の第２抽出手段３Ｂ又はその第２抽出工程において音響粒子速度の合計値Ｕ_allの空間移動速度と閾値Ｖｓとの比較により合計値Ｉ_allの音響インテンシティーから音源インテンシティー成分を切り出した後、音源座標の第２推定手段４Ｂ又はその第２推定工程において音源の座標を推定した場合の結果の一例を、図１１に示す。このときの測定は、上記音響インテンシティーの測定時と同じ条件で行っている。

このように音響粒子速度の合計値Ｕ_allを参照して合計値Ｉ_allの音響インテンシティーから音源インテンシティー成分を切り出した後に音源位置の検出を行った場合でも、図１１の右側部分に示される一部の結果から明らかなように、ｙ軸上にあるはずの音源Ｓ₁の正規の位置（白丸）とほぼ一致した場所に推定されており、これにより音源Ｓ₁の位置が精度よくｙ軸上に推定されていることがわかる。
したがって、実施の形態２に係る音源の検出方法及びその検出装置１Ｂは、音響インテンシティーの空間移動速度だけでなく音響粒子速度の空間移動速度を利用した場合でも、音源の位置を的確に検出することが可能になる。

そして、実施の形態２に係る音源の検出方法（その検出装置１Ｂ）は、図５や図１０及び図１１に示す結果からも明らかなように、実施の形態１に係る音源の検出方法（その検出装置１Ａ）とともに、音源の位置を的確に検出するための有用な手法（又は手段）であるといえる。
その一方で、実施の形態２に係る検出方法（その検出装置１Ｂ）は、音響インテンシティー波形の振幅情報を参照しないことから、以下の長所（ｅ１）、（ｅ２）を更に有しており、その振幅情報を参照する実施の形態１に係る音源の検出方法（その検出装置１Ａ）よりも、音響インテンシティー波形に対して原理的に相性の良い手法（又は手段）であると考えられる。
この点を考慮すると、実施の形態２に係る検出方法（その検出装置１Ｂ）は、実施の形態１に係る検出方法（その検出装置１Ａ）に比べて、音響インテンシティーの測定結果を利用して音源の位置をより的確に検出できることになる。

（ｅ１） 実施の形態２に係る音源の検出方法（その検出装置１Ｂ）は、音響インテンシティーの波形の振幅情報ではなく速度情報を参照するため、例えば、受音点Ｐ（原点）の位置やその近傍の位置に音源が存在し、その音源の音響インテンシティーの波形振幅がほぼゼロの値や、もしくは多くの測定誤差を含む乱れた振幅特性となるような音場に対しても、その音源の位置を推定することが可能である。
（ｅ２） 実施の形態２に係る音源の検出方法（その検出装置１Ｂ）は、音響インテンシティーの波形の振幅情報ではなく速度情報を参照するため、雑音の多い環境の音場において測定された乱れた振幅特性を有する音響インテンシティーからでも、その音源の位置を推定することが可能である。

＜長所（ｅ１）に関する比較試験＞
図１２は、受音点（原点）の近傍位置（ｘ, ｙ, ｚ＝−１５［mm］，０，０）に設置された実音源Ｓ₁の位置の推定を実施したときの試験内容及びその結果を示す。
試験を行った音場は、ｘ，ｙ，ｚ方向の各寸法が２．７ｍ，４．５ｍ，２．５ｍである大きさからなる居間であって、一般的な部屋としての音の響きや反射音が生じる閉鎖空間である。
そして、この閉鎖空間のほぼ中央の位置に、実施の形態１における音響インテンシティーの測定に関する説明のなかで示した構成のマイクロホン６（全指向性マイク６ａ〜６ｄ）を設置し、図１２の左側部分に示すように、原点に対応させる１つのマイク（６ａ：図４）の位置を受音点（Ｐ）として直交３軸方向（ｘ，ｙ，ｚ）を設定した。
また、推定対象の音源Ｓ₁として、スピーカーを使用した。この音源Ｓ₁のスピーカーは、そのスピーカー正面の中心部（○で示す白丸）が上記近傍位置に存在するよう設置した。試験は、音源Ｓ₁のスピーカーから周波数掃引音（スイープ正弦波）を再生し、それをマイクロホン６（全指向性マイク６ａ〜６ｄ）にて録音して得られた音圧波形よりインパルス応答を算出することで行うこととした。

試験は、この条件のもとで、前述した「ピーク検知」を用いた検出方法４と「改良型ピーク検知」を用いた実施の形態１に係る検出方法と「速度検知」を用いた実施の形態２に係る検知方法を、比較するためにそれぞれ実施した。

この比較試験の結果を図１２の中央部及び右側部分に示す。
図１２中の［Ａ］が「ピーク検知」を用いた検出方法４の結果、［Ｂ］が「改良型ピーク検知」を用いた実施の形態１に係る検出方法の結果、［Ｃ］が「速度検知」を用いた実施の形態２に係る検知方法の結果を示している。この［Ａ］、［Ｂ］及び［Ｃ］については、これ以後の比較試験の結果においても同様である。

はじめに、ｙ＝0，ｚ＝0の位置からの音源情報を表す音響インテンシティーＩ_y, Ｉ_zに関して確認する。この点については、理論上はゼロの値になるはずであるところ、実際の測定では図１２の中央部に示されるように［Ａ］、［Ｂ］及び［Ｃ］のいずれの場合であっても僅かなピークを有する波形となる結果が得られている。
これにより、いずれの場合でも音源の存在を確認するために利用する合計値Ｉ_allの音響インテンシティーの波形に対してピーク値が得られていることになるため、「速度検知」を用いた実施の形態２に係る検知方法だけではなく、「ピーク検知」を用いた検出方法４や「改良型ピーク検知」を用いた実施の形態１に係る検出方法であっても、その音源の座標を認識できていることがわかる。

これは、音源Ｓ₁のスピーカーがある程度の大きさを有するものであるため、ｙ＝０、ｚ＝０のピンポイントからだけではなく、そのピンポイントの周囲からも音が供給されていることになり、この結果、音源の位置をｙ＝０、ｚ＝０とは少しずれた場所として認識しているものと考えられる。
しかし、この音源Ｓ₁のスピーカーを極小の音源にした場合は、音響インテンシティーＩ_y, Ｉ_z の測定結果がいずれもゼロの値に近づくことになるため、波形のピーク情報を得ることができず、上記「ピーク検知」や「改良型ピーク検知」を用いた検出方法では、その極小の音源の位置を推定することが難しくなると考えられる。

次に、ｘ＝−１５[mm]というｘ＝０の近傍位置からの音源情報を表す音響インテンシティーＩ_xに関して確認する。この点については、理論上は若干の負のピークを有する音響インテンシティーＩ_xの波形になるはずであるところ、実際の測定では図１２の中央部に示されるように［Ａ］、［Ｂ］及び［Ｃ］のいずれの場合であっても、波形が負の側にとどまらず正の側にも振れる測定誤差を多く含む乱れた波形となる結果が得られている。

以上のことから、最終的な音源Ｓ₁の位置の推定結果は、図１２の右側に示されるような結果として得られた。

まず、「ピーク検知」を用いた検出方法４の場合［Ａ］には、音源Ｓ₁のｘ方向における正規の座標位置が正しくは負側の値になることに対して、ｘ方向の座標位置について誤差成分である正の値として誤認識した推定結果になっている。これは、音響インテンシティーＩ_xの測定波形において最も大きいピークが正の側に現れていることから、その正の側にあるピークを音源の座標情報として検知している結果であると考えられる。
また、「改良型ピーク検知」を用いた実施の形態１に係る検出方法の場合［Ｂ］には、
ピーク幅内の大きさの平均処理より誤差成分が緩和されているが、ｘ方向の座標位置について負の値として認識するまでには至っていない推定結果になっている。これは、音響インテンシティーＩ_xの波形のピーク幅内にある負側の正しい成分が平均処理により正側の誤差成分により相殺されてしまい、ｘ方向の座標位置について負の方向まで引き戻すことができていない結果であると考えられる。

これに対して、「速度検知」を用いた実施の形態２に係る検出方法の場合［Ｃ］には、ｘ方向の座標位置について負側の値として認識されている推定結果になっている。これは、合計値Ｉ_allの音響インテンシティーの波形のうちからゆっくりと移動している最初の波形部分（音源インテンシティーの部分）のみを音源情報として切り出しているため、ｘ方向の音響インテンシティーＩ_x の波形における正しい成分である負の成分だけを切り出すことができている結果であると考えられる。

以上総合すると、「速度検知」を用いた実施の形態２に係る検出方法のみが、受音点の近傍位置に音源Ｓ₁があるような音場であっても、その音源Ｓ₁の存在を正しく検出したうえでその音源Ｓ₁の位置について正しく推定することが可能であることが明らかである。

＜長所（ｅ２）に関する比較試験＞
図１３は、雑音の多い環境の音場で音源の位置の推定を実施したときの試験内容及びその結果を示す。
試験を行った音場は、上記長所（ｅ１）に関する比較試験のときの音場と同じ条件の居間（閉鎖空間）である。そして、この閉鎖空間のほぼ中央の位置に、実施の形態１における音響インテンシティーの測定に関する説明のなかで示した構成のマイクロホン６（全指向性マイク６ａ〜６ｄ）を設置し、図１３の左端部分に示すように、原点に対応させる１つのマイクロホン（６ａ：図４）の位置を受音点（Ｐ）として直交３軸方向（ｘ，ｙ，ｚ）を設定した。
また、試験では、推定対象の音源Ｓ₀として実音源となるスピーカーを使用した。この音源Ｓ₀のスピーカーは、上記音場においてスピーカー正面の中心部が、受音点からの方位角が４５°（正面＋ｙ方向を０°として時計回りの方向になす角度）、仰角が１７．６°となる１．５８ｍ離れた位置に設置した。試験では、このスピーカーから上記長所（ｅ１）に関する比較試験のときと同じ周波数掃引音（スイープ正弦波）を再生した。
さらに、試験では、雑音源Ｓ_nとして音源Ｓ₀とは異なるスピーカーを使用した。この雑音源Ｓ_nのスピーカーは、そのスピーカーの正面中心部が、受音点近傍から方位角が０°（正面＋ｙ方向）、仰角が０°となる１．５ｍ離れた位置に設置した。試験では、このスピーカーから、広域雑音としてのピンクノイズ（パワーが周波数に反比例する雑音）を再生した。

そして、試験は、上記音場において実音源Ｓ₀のスピーカーと雑音源Ｓ_nのスピーカーとからそれぞれ再生された音をマイクロホン６（全指向性マイク６ａ〜６ｄ）で録音し、その録音で得られる音圧波形からインパルス応答を算出することで行った。そして、このときの音源位置の検出における解析は、周波数が１００Ｈｚ〜４ｋＨｚの範囲で行った。
ちなみに、このときの音場は、実音源Ｓ₀の再生音と雑音源Ｓ_nの再生音の測定時におけるＳＮ比がほぼ０［ｄＢ］となり、実音源Ｓ₀の再生音が聴覚（聴感）で確認できないほど大きな雑音が存在している環境であった（図１３の左端部分の中程に示す周波数特性を参照）。
また、試験は、比較するためにこの条件のもとで、前述した「ピーク検知」を用いた検出方法４と「改良型ピーク検知」を用いた実施の形態１に係る検出方法と「速度検知」を用いた実施の形態２に係る検知方法とによりそれぞれ実施した。

この比較試験の結果を図１３の左端部分よりも右側の部分に示す。
はじめに、上記いずれの検出方法の場合（［Ａ］、［Ｂ］及び［Ｃ］）においても、３軸方向の音響インテンシティーＩ_x，Ｉ_y，Ｉ_zがピークの判別が困難なほど振幅特性が大きく乱れた波形になる結果になっている。
しかし、合計値Ｉ_allの音響インテンシティーの波形で雑音による振幅の乱れが多少軽減されていることから、上記いずれの検出方法の場合（［Ａ］、［Ｂ］及び［Ｃ］）においても、音源の存在自体は認識できている。ただし、推定された音源の座標に関しては、図１３の左端を除く右側の下部に示されるように、それぞれの検出方法で互いに異なる結果になっている。
つまり、「ピーク検知」を用いた検出方法４の場合［Ａ］と「改良型ピーク検知」を用いた実施の形態１に係る検出方法の場合［Ｂ］には、推定した音源Ｓ₁の座標位置が推定対象である実音源Ｓ₀の正規の位置からずれた結果になっている。これに対して、「速度検知」を用いた実施の形態２に係る検知方法の場合［Ｃ］には、推定した音源Ｓ₁の座標位置が実音源Ｓ₀の正規の位置とほぼ一致した結果になっている。
したがって、「速度検知」を用いた実施の形態２に係る検知方法が、最も高い精度で実音源Ｓ₀の位置を推定できることがわかる。

以上のことから、「速度検知」を用いた実施の形態２に係る検出方法は、雑音の多い環境の音場であっても、その音源Ｓ₀の存在を正しく検出したうえでその音源Ｓ₀の位置について正しく推定することが可能であることが明らかである。

ちなみに、図１３に示す比較試験１の結果においては、「ピーク検知」を用いた検出方法４の場合［Ａ］の方が、「改良型ピーク検知」を用いた実施の形態１に係る検出方法の場合［Ｂ］よりも精度の高い推定がされている。
これは、幅特性が大きく乱れた音響インテンシティーＩ_x，Ｉ_y，Ｉ_zの測定波形に対しては、ピーク幅内における大きさの平均化の手法を用いる「改良型ピーク検知」を用いた実施の形態１に係る検出方法では、真の音源情報が雑音情報で相殺されてしまうためであると考えられる。

このため、音響インテンシティー波形の振幅情報を利用する「ピーク検知」を用いた検出方法４や「改良型ピーク検知」を用いた実施の形態１に係る検出方法では、音場の雑音や測定環境や音源の位置条件の違いによって、音源位置の推定結果の優劣が入れ替わることもある。
これに対して、「速度検知」を用いた実施の形態２に係る検出方法では、音場の雑音や測定環境や音源の位置条件の違いに左右されず、安定した精度で音源位置の推定が可能である。

＜参考の比較試験＞
図１４は、コンサートホールにおける反射音、すなわち虚音源の分布について、参考までに「ピーク検知」を用いた検出方法４と「改良型ピーク検知」を用いた実施の形態１に係る検出方法と「速度検知」を用いた実施の形態２に係る検知方法とで音源の位置を推定したときの比較試験の結果を示す。

この試験では、コンサートホールのステージに設置した全指向性スピーカー（実音源）から周波数掃引音（スイープ正弦波）を再生し、実音源から観客席に設置した受音点へのインパルス応答を測定し、上記Ｐ−Ｐ法により３軸方向の音響インテンシティーを算出した。そして、その音響インテンシティーの測定結果（図１４の左端部分に示す結果）を同様に利用して上記各検出方法により音源（虚音源）の位置をそれぞれ推定した。この推定された複数の音源の位置を、図１４の左端を除く右側部分のｘ−ｙ平面における円形の内外に併せて示している。

この図１４に示す結果から、上記３つの音源の検出方法では、虚音源の分布として様相が大幅に異なる分布結果が得られることがわかる。このことは、同じコンサートホールであるにもかかわらず、音源の検出方法が変わると、そのコンサートホールを異なる響きを有する空間として評価してしまうことになる。このうち「速度検知」を用いた実施の形態２に係る検知方法の場合［Ｃ］には、最も多くの虚音源の位置を推定することができている。
以上のことから、空間音響の分析を行う際には、精度の高い音源の検出手法を選択して適用することが重要であることがわかる。また、このような観点からすると、上記３つの検出方法のなかでは、「速度検知」を用いた実施の形態２に係る検知方法の場合が最も精度の高い音源位置の検出ができるといえる。

［他の実施の形態］
実施の形態１等では、音響インテンシティーの測定方法としてＰ−Ｐ法の一例を示したが、その測定方法としては、他にも例えば、以下に示す方法を適用してもよい。

例えば、全指向性マイクロホンを使用して測定する別のＰ−Ｐ法を適用してもよい。すなわち、４つの全指向性マイクロホンを正四面体の頂点に配置して測定する方法及びその機器（小野測器社：MI-331）や、２つの対向するマイクロホンからなるマイクロホンペアを直交３軸に配置して計６つの全指向性マイクホンを使用して測定する方法及びその機器（ＲＩＯＮ社：SI-33I、Ｇ.Ｒ.Ａ.Ｓ社 50VI-1など）である。
また、単一指向性マイクを使用して測定するＣ−Ｃ法を適応してもよい（特開２００８−２４９７０２号公報など）。このＣ−Ｃ法を適用する場合は、２つの対向するマイクロホンの出力の和にて音圧を算出し、それらの出力の差分にて音響粒子速度を算出した後、その算出した音圧と音響粒子速度との乗算により音響インテンシティーを得る。

さらに、粒子速度を直接測定することで音響インテンシティーを得る手法もある。この手法を適用する場合は、マイクホンと音響粒子速度センサーを組み合わせたＰＵセンサー（Microflown Technologies社：USP）から得られる音圧と音響粒子速度を乗算することにより音響インテンシティーを得ることができる。なお、Ｐ−Ｐ法やＣ−Ｃ法では、マイクロホンで得られた音圧から近似的に算出される音響粒子速度を用いて音響インテンシティーを得ることができる。

１Ａ，１Ｂ…音源の検出装置
２ …音響インテンシティーの測定手段
３Ａ…音源の第１抽出手段
３Ｂ…音源の第２抽出手段
４Ａ…音源座標の第１推定手段
４Ｂ…音源座標の第２推定手段
Ｉ_x，Ｉ_y，Ｉ_z…直交する３軸方向の音響インテンシティー
Ｉ_all…音響インテンシティーの合計値
Ｕ_x，Ｕ_y，Ｕ_z…直交する３軸方向の音響粒子速度
Ｕ_all…音響粒子速度の合計値

Claims (4)

1. 音場において受音点で直交する３軸方向（ｘ，ｙ，ｚ）の音響インテンシティーＩ_x，Ｉ_y，Ｉ_z［Ｗ／ｍ²］を測定し、それらの合計値であるＩ_all＝√（Ｉ_x ²＋Ｉ_y ²＋Ｉ_z ²）を算出する音響インテンシティーの測定工程と（√（）は括弧内の値の平方根を示す。以下同様。）、
   前記測定工程で得た合計値の音響インテンシティー波形におけるピークが音源を表していると仮定して、そのピークを時間窓幅ごとに抽出し、その抽出されたピークから音源の大きさ及び距離を算出する音源の第１抽出工程と、
   前記測定工程で得た合計値のピーク幅を推定して、そのピーク幅に対応する前記音響インテンシティーＩ_x，Ｉ_y，Ｉ_zにおけるピーク幅内の大きさの符号付の平均値を算出し、前記音響インテンシティーＩ_x，Ｉ_y，Ｉ_zにおけるピーク幅内の符号付の平均値と前記第１抽出工程で得た音源の距離とから音源の座標を推定する音源座標の第１推定工程と、
   を備える音源の検出方法。
2. 音場において受音点で直交する３軸方向（ｘ，ｙ，ｚ）の音響インテンシティーＩ_x，Ｉ_y，Ｉ_z［Ｗ／ｍ²］を測定し、それらの合計値であるＩ_all＝√（Ｉ_x ²＋Ｉ_y ²＋Ｉ_z ²）を算出する音響インテンシティーの測定工程と、
   前記測定工程で得た合計値の音響インテンシティーを音源インテンシティー成分と浮遊インテンシティー成分に区別する空間移動速度の閾値を設定し、前記合計値の音響インテンシティーから前記閾値以下の空間移動速度になる音源インテンシティー成分を切り出すか、又は、前記合計値の音響インテンシティーを音圧で除して得られる音響粒子速度の合計値から前記閾値以下の空間移動速度になる当該音響粒子速度の到来時間と対応する音源インテンシティー成分を前記合計値の音響インテンシティーから切り出し、その切り出された一片として含まれる音源インテンシティー成分の強さ及び到来時間の平均値を算出して音源の大きさ及び距離とする音源の第２抽出工程と、
   前記第２抽出工程で切り出された音源インテンシティー成分の一片と同じ時間幅に含まれる音響インテンシティーＩ_x，Ｉ_y，Ｉ_zを切り出し、その切り出した音響インテンシティーＩ_x，Ｉ_y，Ｉ_zの符号付の大きさの平均値を算出し、前記音響インテンシティーＩ_x，Ｉ_y，Ｉ_zの符号付の平均値と前記第２抽出工程で得た音源の距離とから音源の座標を推定する音源座標の第２推定工程と、
   を備える音源の検出方法。
3. 音場において受音点で直交する３軸方向（ｘ，ｙ，ｚ）の音響インテンシティーＩ_x，Ｉ_y，Ｉ_z［Ｗ／ｍ²］を測定し、それらの合計値であるＩ_all＝√（Ｉ_x ²＋Ｉ_y ²＋Ｉ_z ²）を算出する音響インテンシティーの測定手段と、
   前記測定手段で得た合計値の音響インテンシティー波形におけるピークが音源を表していると仮定して、そのピークを時間窓幅ごとに抽出し、その抽出されたピークから音源の大きさ及び距離を算出する音源の第１抽出手段と、
   前記測定手段で得た合計値のピーク幅を推定して、そのピーク幅に対応する前記音響インテンシティーＩ_x，Ｉ_y，Ｉ_zにおけるピーク幅内の大きさの符号付の平均値を算出し、前記音響インテンシティーＩ_x，Ｉ_y，Ｉ_zにおけるピーク幅内の符号付の平均値と前記第１抽出手段で得た音源の距離とから音源の座標を推定する音源座標の第１推定手段と、
   を備える音源の検出装置。
4. 音場において受音点で直交する３軸方向（ｘ，ｙ，ｚ）の音響インテンシティーＩ_x，Ｉ_y，Ｉ_z［Ｗ／ｍ²］を測定し、それらの合計値であるＩ_all＝√（Ｉ_x ²＋Ｉ_y ²＋Ｉ_z ²）を算出する音響インテンシティーの測定手段と、
   前記測定手段で得た合計値の音響インテンシティーを音源インテンシティー成分と浮遊インテンシティー成分に区別する空間移動速度の閾値を設定し、前記合計値の音響インテ
   ンシティーから前記閾値以下の空間移動速度になる音源インテンシティー成分を切り出すか、又は、前記合計値の音響インテンシティーを音圧で除して得られる音響粒子速度の合計値から前記閾値以下の空間移動速度になる当該音響粒子速度の到来時間と対応する音源インテンシティー成分を前記合計値の音響インテンシティーから切り出し、その切り出された一片として含まれる音源インテンシティー成分の強さ及び到来時間の平均値を算出して音源の大きさ及び距離とする音源の第２抽出手段と、
   前記第２抽出手段で切り出された音源インテンシティー成分の一片と同じ時間幅に含まれる音響インテンシティーＩ_x，Ｉ_y，Ｉ_zを切り出し、その切り出した音響インテンシティーＩ_x，Ｉ_y，Ｉ_zの符号付の大きさの平均値を算出し、前記音響インテンシティーＩ_x，Ｉ_y，Ｉ_zの符号付の平均値と前記第２抽出手段で得た音源の距離とから音源の座標を推定する音源座標の第２推定手段と、
   を備える音源の検出装置。