JP2021135202A

JP2021135202A - 音源探査システム、及びその音源探査方法並びに音源探査プログラム

Info

Publication number: JP2021135202A
Application number: JP2020032257A
Authority: JP
Inventors: 浩之細谷; Hiroyuki Hosoya; 健太郎林; Kentaro Hayashi; 貴裕中野; Takahiro Nakano
Original assignee: Mitsubishi Heavy Industries Ltd
Current assignee: Mitsubishi Heavy Industries Ltd
Priority date: 2020-02-27
Filing date: 2020-02-27
Publication date: 2021-09-13
Also published as: CN113311389A; EP3872515A1

Abstract

【課題】音源探査において汎用性を向上させることのできる音源探査システム、及びその音源探査方法並びに音源探査プログラムを提供することを目的とする。【解決手段】空間Ｓを囲むように配置され、音圧及び粒子速度を検出する複数のセンサ３と、センサ３で検出した音圧及び粒子速度に基づいて、空間Ｓ内における探査対象の周波数に応じた仮想位置の音圧を推定する推定部５と、仮想位置の音圧に基づいて、周波数に対応する音源探査を行う探査部６とを備える。【選択図】図１

Description

本開示は、音源探査システム、及びその音源探査方法並びに音源探査プログラムに関するものである。

音源探査は、アレイを構成するセンサを用いて行われる。アレイは、対象とする周波数に合わせてセンサが配置されている。

アレイではセンサが物理的に配置されるため、対応可能周波数が限定されてしまう。特許文献１では、アレイにおけるセンサの配置位置を物理的に可変とする機構が開示されている。

国際公開第２０１４／１５６２９２号

しかしながら、センサの配置位置を物理的に可変とする場合であっても、アレイの可変形態数によって対応可能周波数が限定されてしまう。このため、対象とする周波数に対して最適なアレイ配置として適切な計測を行うことができない可能性がある。

このため、汎用性のある音源探査システムが求められている。

本開示は、このような事情に鑑みてなされたものであって、音源探査において汎用性を向上させることのできる音源探査システム、及びその音源探査方法並びに音源探査プログラムを提供することを目的とする。

本開示の第１態様は、空間を囲むように配置され、音圧及び粒子速度を検出する複数のセンサと、前記センサで検出した音圧及び粒子速度に基づいて、前記空間内における探査対象の周波数に応じた仮想位置の音圧を推定する推定部と、前記仮想位置の音圧に基づいて、前記周波数に対応する音源探査を行う探査部と、を備える音源探査システムである。

本開示の第２態様は、空間を囲むように配置され、音圧及び粒子速度を検出する複数のセンサで検出した音圧及び粒子速度に基づいて、前記空間内における探査対象の周波数に応じた仮想位置の音圧を推定する工程と、前記仮想位置の音圧に基づいて、前記周波数に対応する音源探査を行う工程と、を有する音源探査方法である。

本開示の第３態様は、空間を囲むように配置され、音圧及び粒子速度を検出する複数のセンサで検出した音圧及び粒子速度に基づいて、前記空間内における探査対象の周波数に応じた仮想位置の音圧を推定する処理と、前記仮想位置の音圧に基づいて、前記周波数に対応する音源探査を行う処理と、をコンピュータに実行させるための音源探査プログラムである。

本開示によれば、音源探査において汎用性を向上させることができるという効果を奏する。

本開示の一実施形態に係る音源探査システムの概略構成を示す図である。本開示の一実施形態に係る情報処理部のハードウェア構成の一例を示した図である。本開示の一実施形態に係るセンサアレイの断面において音場を示す図である。本開示の一実施形態に係るセンサアレイの断面において仮想位置を設定した例を示す図である。本開示の一実施形態に係る音源可視化の例を示す図である。本開示の一実施形態に係る音源探査処理の手順の一例を示すフローチャートである。本開示の一実施形態に係る音源探査処理の変形例の手順の一例を示すフローチャートである。本開示の一実施形態に係る音源探査処理の変形例の手順の一例を示すフローチャートである。

以下に、本開示に係る音源探査システム、及びその音源探査方法並びに音源探査プログラムの一実施形態について、図面を参照して説明する。

図１は、本開示の一実施形態に係る音源探査システム１の概略構成を示す図である。図１に示すように、本実施形態に係る音源探査システム１は、センサアレイ２と、情報処理部４とを主な構成として備えている。

センサアレイ２は、複数のセンサ３を有している。センサ３は、空間Ｓを囲むように配置されている。本実施形態では、センサ３によって取り囲まれた空間Ｓは球状である場合を例として説明する。センサ３で囲んだ空間Ｓは閉空間となるため、センサ３は例えば４つ以上設けられることが好ましい。なお、空間Ｓについては、球状に限定されない。

配置された各センサ３は、それぞれ音圧及び粒子速度を検出する。音圧は、センサ３が配置された位置における音による圧力である。粒子速度は、センサ３が配置された位置における音響粒子の速度である。換言すると、粒子速度は、空気の振動である。また、後述するように粒子速度は、音圧勾配に相当する。そして、各センサ３で検出された音圧及び粒子速度は、情報処理部４へ出力される。

センサアレイ２として配置されるセンサ３は、図１に示すように所定の空間Ｓを取り囲むように配置されている。センサ３の配置位置や配置密度等は図１に限定されない。

センサアレイ２として配置されるセンサ３の間隔については、対象とする周波数で想定される波長の１／２程度に設定されることが好ましい。具体的には、対象とする最高周波数の波長の１／２がセンサ３の間隔の最大値として設定されることが好ましい。

情報処理部４は、センサ３からの情報に基づいて任意の位置の音圧を推定して音源探査を行う。音響探査を行う場合には、探査対象とする周波数に対応した位置（間隔）の音圧を計測する必要がある。しかしながら、物理的に配置したセンサ３の計測値をそのまま使用する場合には、音圧の計測位置が固定されてしまうため、対応可能な周波数が限定されてしまう。そこで、情報処理部４では、配置したセンサ３の計測値を用いて仮想的な位置（特にセンサ３を配置していない位置）の音圧を推定することによって、様々な周波数に対応して音響探査を行う。

図２は、本実施形態に係る情報処理部４のハードウェア構成の一例を示した図である。
図２に示すように、情報処理部４は、コンピュータシステム（計算機システム）であり、例えば、ＣＰＵ１１と、ＣＰＵ１１が実行するプログラム等を記憶するためのＲＯＭ（ＲｅａｄＯｎｌｙＭｅｍｏｒｙ）１２と、各プログラム実行時のワーク領域として機能するＲＡＭ（ＲａｎｄｏｍＡｃｃｅｓｓＭｅｍｏｒｙ）１３と、大容量記憶装置としてのハードディスクドライブ（ＨＤＤ）１４と、ネットワーク等に接続するための通信部１５とを備えている。なお、大容量記憶装置としては、ソリッドステートドライブ（ＳＳＤ）を用いることとしてもよい。これら各部は、バス１８を介して接続されている。

また、情報処理部４は、キーボードやマウス等からなる入力部や、データを表示する液晶表示装置等からなる表示部などを備えていてもよい。

なお、ＣＰＵ１１が実行するプログラム等を記憶するための記憶媒体は、ＲＯＭ１２に限られない。例えば、磁気ディスク、光磁気ディスク、半導体メモリ等の他の補助記憶装置であってもよい。

後述の各種機能を実現するための一連の処理の過程は、プログラムの形式でハードディスクドライブ１４等に記録されており、このプログラムをＣＰＵ１１がＲＡＭ１３等に読み出して、情報の加工・演算処理を実行することにより、後述の各種機能が実現される。なお、プログラムは、ＲＯＭ１２やその他の記憶媒体に予めインストールしておく形態や、コンピュータ読み取り可能な記憶媒体に記憶された状態で提供される形態、有線又は無線による通信手段を介して配信される形態等が適用されてもよい。コンピュータ読み取り可能な記憶媒体とは、磁気ディスク、光磁気ディスク、ＣＤ−ＲＯＭ、ＤＶＤ−ＲＯＭ、半導体メモリ等である。

図１に示されるように、情報処理部４は、推定部５と、探査部６とを備えている。

推定部５は、センサ３で検出した音圧及び粒子速度に基づいて、空間Ｓ内における探査対象の周波数に応じた仮想位置ＶＰの音圧を推定する。具体的には、推定部５は、センサ３によって検出した音圧及び粒子速度に基づいて、キルヒホッフ−ヘルムホルツ（Ｋｉｒｃｈｈｏｆｆ−Ｈｅｌｍｈｏｌｔｚ）の積分定理により、空間Ｓの内部の音圧を推定する。

キルヒホッフ−ヘルムホルツの積分定理とは、空間Ｓの境界（空間Ｓの外周面）から、空間Ｓの内部の音場を再生する。具体的には、キルヒホッフ−ヘルムホルツの積分定理は、以下の式（１）によって記述される。

上記の式（１）において、Ωは空間Ｓを示しており、∂Ωは空間Ｓの境界面（センサ３が配置された外周面）を示しており、ｒｓは境界面（センサ３の位置）を示す位置ベクトルを示しており、ｒは空間Ｓ内の任意の位置を示す位置ベクトルを示しており、ｐは音圧を示しており、ｎは、∂Ω上における外向き法線ベクトルを示している。すなわち、ｐ（ｒｓ）は、センサ３の位置における音圧を示している。そして、ｐ（ｒ）は、空間Ｓ内における任意の位置の音圧を示している。また、∂ｐ／∂ｎは、音圧勾配を示している。音圧勾配は、粒子速度に相当する。具体的には、音圧勾配と粒子速度とｖの関係は、∂ｐ／∂ｎ＝ｊρωｖとして表される。なお、ｊは虚数単位であり、ρは媒質（空気）の密度であり、ωは角周波数である。Ｇは、以下の式（２）に示すようにグリーン関数を示している。

式（２）において、ｋは、波数である。式（１）及び式（２）によれば、センサ位置における音圧ｐ（ｒｓ）及び音圧勾配∂ｐ／∂ｎ（粒子速度）に基づいて、空間Ｓの内部の任意の位置ｒの音圧ｐ（ｒ）を導出することができる。このため、キルヒホッフ−ヘルムホルツの積分定理を用いることで、例えば図３に示すように、空間Ｓ内の音場（音圧状態）を再生することができる。

このように、空間Ｓ内部の音場を再生することができるため、推定部５は、空間Ｓ内において探査対象の周波数に応じた仮想位置ＶＰを設定する。仮想位置ＶＰとは、音圧を推定する位置である。探査対象の周波数の音の検知感度は、例えば検知位置や検知位置の間隔等に依存する。このため、仮想位置ＶＰについては、探索対象の周波数に対応して、適切に検知可能なように予め仮想位置ＶＰの設定数は配置位置、間隔等が設定されている。

なお、探査において指向性を考慮する場合には、仮想位置ＶＰは、探査対象の周波数及び探査方向に基づいて空間Ｓ内に複数設定されることが好ましい。このように仮想位置ＶＰを設定することで、音源の探査を対象とする方向である探査方向に対して指向性をもたせることができる。具体的には、探査方向における雑音に対する感度を向上させることができる。

例えば、仮想位置ＶＰは、図４のように空間Ｓ内に設定される。図４では、任意の周波数に応じて仮想位置ＶＰを設定した例である。このように仮想位置ＶＰが設定されると、式（１）を用いて各仮想位置ＶＰに対応する音圧が推定される。

このように、空間Ｓの周囲のセンサ３の計測結果を用いて空間Ｓの内部の任意の点の音圧を取得することができるため、仮想位置ＶＰを適宜変更することによってさまざまな周波数に対応して音圧を取得することが可能となる。

なお、仮想位置ＶＰについては、周波数や探査方向に限定されず、用途等の他の要因にも基づいて設定されることとしてもよい。

探査部６は、仮想位置ＶＰの音圧に基づいて、周波数に対応する音源探査を行う。具体的には、探査部６は、複数の仮想位置ＶＰの音圧に基づいてビームフォーミングを行い、音源探査を行う。各仮想位置ＶＰのそれぞれに届く信号に遅延などの処理を行うことで、等価的に非常に鋭い指向性を持たせ、特定の方向から到来する音波の強度を求める。そして、鋭い指向性の方向を変化させ、探査方向におけるエリアをスキャンして音源の位置及び強度を推定する。音源探査が行われると、探査方向におけるエリアがスキャンされるため、該エリア内における音源の位置及び強度が可視化可能となる。

図５は、探査対象をファン８として、センサアレイ２に対するファン８の方向を探査方向とし音源探査（音源可視化）を行った場合の例を示している。図５では、音源の強度を４段階として示している。図５の例では、ファン８におけるＰの位置が音源と特定することができる。このように音源探査が行われることによって、音源を特定することが可能となる。

次に、上述の情報処理部４による音源探査処理の一例について図６を参照して説明する。図６は、本実施形態に係る音源探査処理の手順の一例を示すフローチャートである。図６に示すフローは、例えば、使用者等によって音源探査の開始指示がされた場合に実行される。

まず、各センサ３から計測結果を取得する（Ｓ１０１）。

次に、空間Ｓ内に設定された仮想位置ＶＰにおける音圧を再生する（Ｓ１０２）。

次に、再生した音圧に基づいて、音源探査を行い、音源可視化を行う（Ｓ１０３）。

なお、上記例では、空間Ｓの外周面の音圧としてセンサ３の計測結果を用いてキルヒホッフ−ヘルムホルツの積分定理を適用する場合を説明したが、外周面において推定された音圧分布を用いてキルヒホッフ−ヘルムホルツの積分定理を適用することとしてもよい。例えば、空間Ｓを球体状として、センサ３が球面上に配置されている場合には、センサ３で検出した音圧及び粒子速度に基づいて、球面上の音圧分布を推定し、球面上の音圧分布に基づいて空間Ｓ内における仮想位置ＶＰの音圧を推定することとしてもよい。具体的には、球面調和関数（球面上でのフーリエ変換）を基底関数とする圧縮センシングを用いることによって、球面上に配置されたセンサ３の計測結果から球面上の音圧分布を表現（推定）することができる。すなわち、球面の音圧分布が、球面調和関数の重ね合わせで表現される。なお、球面上の音圧分布を推定することができれば、球面調和関数を基底関数とする圧縮センシングを用いる場合に限定されない。

このような場合には、音源探査処理のフローは図７のようになる。なお、図６のフローと同様の処理については同一符号を付している。すなわち、Ｓ１０１の後に各センサ３からの計測結果に基づいて、球面上の音圧分布を推定する（Ｓ２０２）。そして、音圧分布に基づいて、空間Ｓ内に設定された仮想位置ＶＰにおける音圧を再生する（Ｓ１０２）。このようにすることで、音圧の推定精度を向上させることができる。すなわち、仮想位置ＶＰの間隔を狭く設定することができ、より高い周波数にも対応することができる。

また、空間Ｓを直方体状として、センサ３が直方体上に配置されている場合には、センサ３で検出した音圧及び粒子速度に基づいて、直方体上の音圧分布を推定し、直方体上の音圧分布に基づいて空間Ｓ内における仮想位置ＶＰの音圧を推定することとしてもよい。具体的には、正弦波を基底関数とする圧縮センシングを用いることによって、直方体上に配置されたセンサ３の計測結果から直方体上の音圧分布を表現（推定）することができる。すなわち、直方体上の音圧分布が、波長が異なる正弦波の重ね合わせで表現される。なお、直方体上の音圧分布を推定することができれば、正弦波を基底関数とする圧縮センシングを用いる場合に限定されない。

このような場合には、音源探査処理のフローは図８のようになる。なお、図６のフローと同様の処理については同一符号を付している。すなわち、Ｓ１０１の後に各センサ３からの計測結果に基づいて、直方体上の音圧分布を推定する（Ｓ３０２）。そして、音圧分布に基づいて、空間Ｓ内に設定された仮想位置ＶＰにおける音圧を再生する（Ｓ１０２）。このようにすることで、音圧の推定精度を向上させることができる。すなわち、仮想位置ＶＰの間隔を狭く設定することができ、より高い周波数にも対応することができる。

以上説明したように、本実施形態に係る音源探査システム、及びその音源探査方法並びに音源探査プログラムによれば、空間Ｓを囲む複数のセンサ３によって検出した音圧及び粒子速度により、囲まれている空間Ｓの内部の音圧を推定することが可能となる。すなわち、探査対象である周波数に応じて空間Ｓ内に設定された位置である（センサ３が配置されていない位置）の音圧を推定することが可能となる。そして仮想位置ＶＰの音圧に基づくことで、周波数に対する音源探査を行うことができる。

つまり、センサ３の配置位置の変更を行うことなく、任意の周波数を探査対象として音源探査を行うことが可能となる。換言すると、周波数に応じたセンサアレイを仮想的に再現し、音源探査を行うことが可能となる。すなわち、音源探査において汎用性を向上させることができる。

本開示は、上述の実施形態のみに限定されるものではなく、発明の要旨を逸脱しない範囲において、種々変形実施が可能である。

以上説明した各実施形態に記載の音源探査システム、及びその音源探査方法並びに音源探査プログラムは例えば以下のように把握される。
本開示に係る音源探査システム（１）は、空間（Ｓ）を囲むように配置され、音圧及び粒子速度を検出する複数のセンサ（３）と、前記センサ（３）で検出した音圧及び粒子速度に基づいて、前記空間（Ｓ）内における探査対象の周波数に応じた仮想位置（ＶＰ）の音圧を推定する推定部（５）と、前記仮想位置（ＶＰ）の音圧に基づいて、前記周波数に対応する音源探査を行う探査部（６）と、を備える。

本開示に係る音源探査システム（１）によれば、空間（Ｓ）を囲む複数のセンサ（３）によって検出した音圧及び粒子速度により、囲まれている空間（Ｓ）の内部の音圧を推定することが可能となる。すなわち、探査対象である周波数に応じて空間（Ｓ）内に設定された位置である仮想位置（センサ（３）が配置されていない位置）の音圧を推定することが可能となる。そして仮想位置の音圧に基づくことで、周波数に対する音源探査を行うことができる。つまり、センサ（３）の配置位置の変更を行うことなく、任意の周波数を探査対象として音源探査を行うことが可能となる。換言すると、周波数に応じたセンサアレイを仮想的に再現し、音源探査を行うことが可能となる。すなわち、音源探査において汎用性を向上させることができる。

本開示に係る音源探査システム（１）は、前記推定部は、前記センサによって検出した前記空間の境界における音圧及び粒子速度に基づいて、前記空間の内部の音圧を推定することとしてもよい。

本開示に係る音源探査システム（１）によれば、センサによって検出した空間の境界（外周面）における音圧及び粒子速度に基づくことによって、空間の内部の音圧を推定することが可能となる。

本開示に係る音源探査システム（１）は、前記推定部（５）は、前記空間（Ｓ）の周囲に配置された前記センサ（３）によって検出した音圧及び粒子速度に基づいて、キルヒホッフ−ヘルムホルツの積分定理により、前記空間（Ｓ）の内部の音圧を推定することとしてもよい。

本開示に係る音源探査システム（１）によれば、キルヒホッフ−ヘルムホルツの積分定理を用いることで、空間（Ｓ）の周囲に配置されたセンサ（３）によって検出した音圧及び粒子速度に基づいて、空間（Ｓ）の内部の音圧を推定することができる。

本開示に係る音源探査システム（１）は、前記仮想位置（ＶＰ）は、探査対象の周波数及び探査方向に基づいて前記空間（Ｓ）内に複数設定されることとしてもよい。

本開示に係る音源探査システム（１）によれば、仮想位置（ＶＰ）は、探査対象の周波数及び探査方向に基づいて空間（Ｓ）内に複数設定されることで、センサ（３）の配置位置を変更することなく、探査対象の周波数及び探査方向に対応して音源探査を行うことができる。

本開示に係る音源探査システム（１）は、前記探査部（６）は、複数の前記仮想位置（ＶＰ）の音圧に基づいてビームフォーミングを行い、音源探査を行うこととしてもよい。

本開示に係る音源探査システム（１）によれば、複数の仮想位置（ＶＰ）の音圧に基づいてビームフォーミングを行うことで、音源探査を行うことができる。

本開示に係る音源探査システム（１）は、前記センサ（３）は、球面上に配置されており、前記推定部（５）は、前記センサ（３）で検出した音圧及び粒子速度に基づいて、前記球面上の音圧分布を推定し、前記球面上の音圧分布に基づいて前記空間（Ｓ）内における前記仮想位置（ＶＰ）の音圧を推定することとしてもよい。

本開示に係る音源探査システム（１）によれば、センサ（３）が球面上に配置されている場合には、センサ（３）で検出した音圧及び粒子速度に基づいて、球面上の音圧分布を推定することができる。そして、センサ（３）で検出した音圧及び粒子速度に基づいて、球面上の音圧分布を推定することで、音圧の推定精度を向上させることができる。すなわち、仮想位置（ＶＰ）の間隔を狭く設定することができ、より高い周波数にも対応することができる。

本開示に係る音源探査システム（１）は、前記推定部は、球面調和関数を基底関数として、前記球面上の音圧分布を推定することとしてもよい。

本開示に係る音源探査システム（１）によれば、球面調和関数を基底関数とすることで球面上の音圧分布を推定することができる。

本開示に係る音源探査システム（１）は、前記センサ（３）は、直方体上に配置されており、前記推定部（５）は、前記センサ（３）で検出した音圧及び粒子速度に基づいて、前記直方体上の音圧分布を推定し、前記直方体上の音圧分布に基づいて前記空間（Ｓ）内における前記仮想位置（ＶＰ）の音圧を推定することとしてもよい。

本開示に係る音源探査システム（１）によれば、センサ（３）が直方体上に配置されている場合には、センサ（３）で検出した音圧及び粒子速度に基づいて、直方体上の音圧分布を推定することができる。そして、センサ（３）で検出した音圧及び粒子速度に基づいて、直方体上の音圧分布を推定することで、音圧の推定精度を向上させることができる。すなわち、仮想位置（ＶＰ）の間隔を狭く設定することができ、より高い周波数にも対応することができる。

本開示に係る音源探査システム（１）は、前記推定部は、正弦波を基底関数として、前記直方体上の音圧分布を推定することとしてもよい。

本開示に係る音源探査システム（１）によれば、正弦波を基底関数とすることで直方体上の音圧分布を推定することができる。

本開示に係る音源探査方法は、空間（Ｓ）を囲むように配置され、音圧及び粒子速度を検出する複数のセンサ（３）で検出した音圧及び粒子速度に基づいて、前記空間（Ｓ）内における探査対象の周波数に応じた仮想位置（ＶＰ）の音圧を推定する工程と、前記仮想位置（ＶＰ）の音圧に基づいて、前記周波数に対応する音源探査を行う工程と、を有する。

本開示に係る音源探査プログラムは、空間（Ｓ）を囲むように配置され、音圧及び粒子速度を検出する複数のセンサ（３）で検出した音圧及び粒子速度に基づいて、前記空間（Ｓ）内における探査対象の周波数に応じた仮想位置（ＶＰ）の音圧を推定する処理と、前記仮想位置（ＶＰ）の音圧に基づいて、前記周波数に対応する音源探査を行う処理と、をコンピュータに実行させる。

１：音源探査システム
２：センサアレイ
３：センサ
４：情報処理部
５：推定部
６：探査部
８：ファン
１１：ＣＰＵ
１２：ＲＯＭ
１３：ＲＡＭ
１４：ハードディスクドライブ
１５：通信部
１８：バス
Ｓ：空間
ＶＰ：仮想位置

Claims

空間を囲むように配置され、音圧及び粒子速度を検出する複数のセンサと、
前記センサで検出した音圧及び粒子速度に基づいて、前記空間内における探査対象の周波数に応じた仮想位置の音圧を推定する推定部と、
前記仮想位置の音圧に基づいて、前記周波数に対応する音源探査を行う探査部と、
を備える音源探査システム。
前記推定部は、前記センサによって検出した前記空間の境界における音圧及び粒子速度に基づいて、前記空間の内部の音圧を推定する請求項１に記載の音源探査システム。
前記推定部は、キルヒホッフ−ヘルムホルツの積分定理により、前記空間の内部の音圧を推定する請求項２に記載の音源探査システム。
前記仮想位置は、探査対象の周波数及び探査方向に基づいて前記空間内に複数設定される請求項１から３のいずれか１項に記載の音源探査システム。
前記探査部は、複数の前記仮想位置の音圧に基づいてビームフォーミングを行い、音源探査を行う請求項１から４のいずれか１項に記載の音源探査システム。
前記センサは、球面上に配置されており、
前記推定部は、前記センサで検出した音圧及び粒子速度に基づいて、前記球面上の音圧分布を推定し、前記球面上の音圧分布に基づいて前記空間内における前記仮想位置の音圧を推定する請求項１から５のいずれか１項に記載の音源探査システム。
前記推定部は、球面調和関数を基底関数として、前記球面上の音圧分布を推定する請求項６に記載の音源探査システム。
前記センサは、直方体上に配置されており、
前記推定部は、前記センサで検出した音圧及び粒子速度に基づいて、前記直方体上の音圧分布を推定し、前記直方体上の音圧分布に基づいて前記空間内における前記仮想位置の音圧を推定する請求項１から５のいずれか１項に記載の音源探査システム。
前記推定部は、正弦波を基底関数として、前記直方体上の音圧分布を推定する請求項８に記載の音源探査システム。
空間を囲むように配置され、音圧及び粒子速度を検出する複数のセンサで検出した音圧及び粒子速度に基づいて、前記空間内における探査対象の周波数に応じた仮想位置の音圧を推定する工程と、
前記仮想位置の音圧に基づいて、前記周波数に対応する音源探査を行う工程と、
を有する音源探査方法。
空間を囲むように配置され、音圧及び粒子速度を検出する複数のセンサで検出した音圧及び粒子速度に基づいて、前記空間内における探査対象の周波数に応じた仮想位置の音圧を推定する処理と、
前記仮想位置の音圧に基づいて、前記周波数に対応する音源探査を行う処理と、
をコンピュータに実行させるための音源探査プログラム。