特許 6100655 音源探査装置

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】6100655

(24)【登録日】2017年3月3日

(45)【発行日】2017年3月22日

(54)【発明の名称】音源探査装置

(51)【国際特許分類】

   G01H 3/00 20060101AFI20170313BHJP

【ＦＩ】

   G01H3/00 A

【請求項の数】4

【全頁数】10

(21)【出願番号】特願2013-192815(P2013-192815)

(22)【出願日】2013年9月18日

(65)【公開番号】特開2015-59794(P2015-59794A)

(43)【公開日】2015年3月30日

【審査請求日】2016年7月1日

(73)【特許権者】

【識別番号】000145806

【氏名又は名称】株式会社小野測器

(74)【代理人】

【識別番号】100094330

【弁理士】

【氏名又は名称】山田 正紀

(74)【代理人】

【識別番号】100079175

【弁理士】

【氏名又は名称】小杉 佳男

(72)【発明者】

【氏名】猿渡 克己

(72)【発明者】

【氏名】三津山 元基

【審査官】 田中 秀直

(56)【参考文献】

【文献】 特開平０５−０７２０２５（ＪＰ，Ａ）

【文献】 特開２０１３−０１５４６８（ＪＰ，Ａ）

【文献】 特開２０１１−２０３１６２（ＪＰ，Ａ）

【文献】 特開平０９−２５１２９９（ＪＰ，Ａ）

【文献】 特開２００６−２７９７４０（ＪＰ，Ａ）

【文献】 特開平０３−１７９２８５（ＪＰ，Ａ）

【文献】 特開平０４−１８４１８３（ＪＰ，Ａ）

【文献】 特開２００３−２０２２７６（ＪＰ，Ａ）

【文献】 特開平１１−１６８７９１（ＪＰ，Ａ）

【文献】 特開２００５−２８４４４３（ＪＰ，Ａ）

【文献】 特開平０８−１５２４６７（ＪＰ，Ａ）

(58)【調査した分野】（Int.Cl.，ＤＢ名）

Ｇ０１Ｈ １／００−１７／００

(57)【特許請求の範囲】

【請求項1】

配列された複数のマイクロホンを備えたプローブと、
前記プローブで受音して得た音信号に基づいて、空間上の関心領域内に離散的に設定した複数の初期の測定点の音圧を算出し隣接する測定点の音圧どうしの差分からなる差分音圧を算出して該差分音圧を閾値と比較し、さらに
閾値以上の差分音圧が得られた隣接する２つの測定点どうしの間に新たな測定点を設定して該新たな測定点の音圧を算出し、該新たな測定点を含む隣接する測定点の音圧どうしの差分からなる差分音圧を算出して閾値と比較する過程を、前記関心領域に関する全ての差分音圧が閾値以下となるまで繰り返す音源探査演算を実行する演算部とを有することを特徴とする音源探査装置。

【請求項2】

前記関心領域が、前記プローブ側から覗めたときに２次元的に広がる空間であって、前記演算部が、前記音源探査演算を、互いに異なる２方向それぞれについて独立に実行するものであることを特徴とする請求項１記載の音源探査装置。

【請求項3】

前記関心領域が、前記プローブ側から覗めたときに２次元的に広がる空間であって、前記演算部が、２次元差分演算を用いて前記新たな測定点を設定するものであることを特徴とする請求項１記載の音源探査装置。

【請求項4】

前記演算部は、前記複数の初期の測定点を、隣接する測定点どうしの間隔が関心周波数領域の最高周波数あるいは関心周波数の音波の波長の１／２以下の間隔となるように設定するものであることを特徴とする請求項１から３のうちいずれか１項記載の音源探査装置。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、空間上の関心領域内の音源位置を探査する音源探査装置に関する。

【背景技術】

【0002】

従来より、複数のマイクロホンで受音して得た音信号から音源位置を探査する技術が知られており、そのうちの典型的な技術の１つとしてビームフォーミング式の音源探査技術が存在する。

【0003】

このビームフォーミング式を採用して音源探査を行なうにあたっては、関心領域内の離散的な複数の測定点それぞれについて音圧を算出し、その音圧のピークの位置が音源位置として同定される（特許文献１参照）。

【0004】

ここで音源信号を正確に探査しようとすると、関心領域内を細かく分割して多数の測定点を設定し、それら多数の測定点それぞれについて音圧を算出する必要がある。ただし、この場合、演算時間が長くかかり、リアルタイム性が失われるおそれがある。一方、リアルタイム性を向上させようとして測定点どうしの間隔を空けて粗い分割とすると、音源位置の正確性が失われる結果となる。

【先行技術文献】

【特許文献】

【0005】

【特許文献1】特開平５−７２０２５号公報

【発明の概要】

【発明が解決しようとする課題】

【0006】

本発明は、音源探査のリアルタイム性の確保と音源位置同定の正確性を高度に両立させた音源探査装置を提供することを目的とする。

【課題を解決するための手段】

【0007】

上記目的を達成する本発明の音源探査装置は、
配列された複数のマイクロホンを備えたプローブと、
上記プローブで受音して得た音信号に基づいて、空間上の関心領域内に離散的に設定した複数の初期の測定点の音圧を算出し隣接する測定点の音圧どうしの差分からなる差分音圧を算出してその差分音圧を閾値と比較し、さらに
閾値以上の差分音圧が得られた隣接する２つの測定点どうしの間に新たな測定点を設定してその新たな測定点の音圧を算出し、その新たな測定点を含む隣接する測定点の音圧どうしの差分からなる差分音圧を算出して閾値と比較する過程を、関心領域に関する全ての差分音圧が閾値以下となるまで繰り返す音源探査演算を実行する演算部とを有することを特徴とする。

【0008】

本発明の音源探査装置は、関心領域内の、閾値を越える差分音圧が得られた領域のみ新たな測定点が設定されて詳細な音圧分布が測定される。このため、少ない演算量で音源位置が正確に同定される。

【0009】

ここで、本発明の音源探査装置において、上記関心領域が、プローブ側から覗めたときに２次元的に広がる空間であって、上記演算部が、音源探査演算を、互いに異なる２方向それぞれについて独立に実行するものであってもよい。

【0010】

この場合、シンプルな演算で済む。

【0011】

あるいは、本発明の音源探査装置において、上記関心領域が、プローブ側から覗めたときに２次元的に広がる空間であって、上記演算部が、２次元差分演算を用いて新たな測定点を設定するものであってもよい。

【0012】

この場合、２次元的な差分音圧が得られる。

【0013】

また、本発明の音源探査装置において、上記演算部は、複数の初期の測定点を、隣接する測定点どうしの間隔が関心周波数領域の最高周波数あるいは関心周波数の音波の波長の１／２以下の間隔となるように設定するものであることが好ましい。

【0014】

この場合、小さな差分音圧が得られた２つの測定点どうしの間に存在する音圧のピークの見逃しが防止される。

【発明の効果】

【0015】

以上の本発明によれば、音源探査のリアルタイム性の確保と音源位置同定の正確性を高度に両立させることができる。

【図面の簡単な説明】

【0016】

【図1】本発明の一実施形態としての音源探査装置のブロック図である。

【図2】プローブを構成する複数のマイクロホンの配列の一例を示した図である。

【図3】ビームフォーミングの原理説明図である。

【図4】演算部で実行される音源探査演算の演算内容を示すフローチャートである。

【図5】音圧分布および測定点の１次元断面を示した模式図である。

【図6】差分音圧の配列を示した図である。

【図7】２次元的な関心領域内の音圧分布を可視化して示した図である。

【図8】関心領域内の差分音圧計算の他の例を示した図である。

【発明を実施するための形態】

【0017】

以下、本発明の実施の形態を説明する。

【0018】

図１は、本発明の一実施形態としての音源探査装置のブロック図である。

【0019】

この音源探査装置１は、プローブ１０と演算部２０とを有する。

【0020】

プローブ１０は、配列された複数のマイクロホン１２（図２参照）を備えている。詳細は後述する。

【0021】

演算部２０は、プローブ１０での受音により得られた音信号に基づいて、音源が含まれていると考えられる関心領域Ｄ内に離散的に設定された測定点Ｐの音圧を算出し、さらに音源探査のための演算を行なう。関心領域の測定点Ｐの設定は、初期的には、測定者や上位の装置等から設定される関心周波数ｆ_０あるいは関心周波数領域ｆ_１〜ｆ_２（ｆ_１＜ｆ_２）の最高周波数ｆ_２の音波の波長の、例えば１／２の間隔となるように設定される。この演算部２０による演算内容の詳細についても後述する。

【0022】

図２は、プローブを構成する複数のマイクロホンの配列の一例を示した図である。

【0023】

ここでは、中心点に、このプローブ１０の向きを視認するためのカメラ１１が配置されており、その中心点を取り巻くように３６個のマイクロホン１２が配列されている。本実施形態のプローブ１０は、４個のマイクロホン１２ａ，１２ｂ，１２ｃ，１２ｄを、それらのマイクロホンどうしの間隔が不均等となるように配置し、それら４個のマイクロホン１２ａ，１２ｂ，１２ｃ，１２ｄからなる組を、円周方向に９等分した各領域それぞれに配置している。

【0024】

図３は、ビームフォーミングの原理説明図である。ビームフォーミング自体は広く知られた技術であり、ここではその概要を説明するにとどめる。

【0025】

ここでは、マイクロホン１２として無指向性のマイクロホンが使われている。

【0026】

図３（Ａ）は、プローブ１０の正面に音源が存在する場合の模式図である。

【0027】

音源がプローブ１０の正面に存在すると、音源から各マイクロホン１２までの距離は互いに等しく、それら複数のマイクロホン１２での受音で得られた音信号を互いに加算すると位相が等しいために大きな振幅の信号が得られることになる。

【0028】

図３（Ｂ）は、プローブ１０が向いた方向に対する斜めの方向に音源が存在する場合の模式図である。

【0029】

音源がプローブ１０に対し斜めの位置にあると、音源と各マイクロホン１２との間の音伝播路の長さがそれぞれ異なる。すなわち、ある１つの時刻に音源から発せられた音が各マイクロホンに到達する時刻が異なる。言い換えると、各マイクロホン１２で同時に受音された音信号は位相が異なっている。したがってそれら複数のマイクロホンでの受音で得られた音信号を互いに加算すると、互いに打ち消し合って小さな信号となってしまう。すなわち、複数のマイクロホン１２で得られた音信号を互いに加算すると図３（Ｃ）に示すように、全体として指向性を待つことになる。

【0030】

ところで、今度は複数のマイクロホンでの受音で得られた音信号を互いの位相を調整して加算すると、複数の測定点Ｐ（図１参照）のうちの、その位相の調整の仕方に応じた測定点Ｐに指向の方向を合わせることができる。

【0031】

ここで、音源から発せられる音が高周波になると経路差と周期が合致し音源が存在しない位置に音源が存在するかのような信号（いわゆるゴースト）が得られることが知られている。図２に示すマイクロホン１２どうしの間隔を不均等にするのは、このゴーストの発生を抑えるためである。例えば、マイクロホン１２を等間隔に配置した場合の測定可能周波数が５００Ｈｚ〜３０００Ｈｚである場合に、図２に示すように不均等に配置することにより、測定可能周波数を５００Ｈｚ〜１０ｋＨｚとすることができる。

【0032】

図４は、演算部で実行される音源探査演算の演算内容を示すフローチャートである。

【0033】

図５は、音圧分布および測定点の１次元断面を示した模式図である。図５の横軸は関心領域内の座標位置［ｍ］、縦軸は音圧［Ｐａ］を表わしている。

【0034】

さらに、図６は、差分音圧の配列を示した図である。

【0035】

この図６の横軸は、図５の横軸と同じ座標位置［ｍ］、縦軸は差分の音圧［Ｐａ］を表わしている。

【0036】

音源探査にあたっては、先ず、初期測定点が設定される（ステップＳ０１）。ここでは図１に示すように、演算部２０に、関心のある周波数ｆ_０あるいは関心のある周波数帯域（ｆ_１〜ｆ_２）のうちの最高周波数ｆ_２が設定され、演算部２０は隣接する測定点どうしの間隔がその関心周波数ｆ_０あるいは、最高周波数ｆ_２の音の波長の１／２の間隔となるように、初期の測定点が設定される。

【0037】

図５には、この初期の測定点Ｐ１，Ｐ２，・・・，Ｐ１５の音圧が四角形の記号で示されている。

【0038】

次にプローブ１０を構成する複数のマイクロホン１２での受音で得られた音信号を取得し（ステップＳ０２）、その音信号に基づいて、ステップＳ０１で設定された初期測定点それぞれについての音圧を算出し（ステップＳ０３）、さらに隣接する測定点の音圧どうしの差分からなる差分音圧を算出する（ステップＳ０４）。

【0039】

図６（Ａ）は、この１回目の計算における差分音圧を示している。

【0040】

次にこの差分音圧と閾値とを比較する（ステップＳ０５）。そして閾値以上の差分音圧が存在していたら（ステップＳ０６）、その閾値を越える差分音圧が得られた、隣接する２つの測定点どうしの間に、新たな測定点を設定する（ステップＳ０７）。

【0041】

本実施形態では、隣接する２つの測定点の中点に新たな測定点が設定される。図５では、この新たな測定点は、○印で示される「２回目の測定点」Ｐ２．５，Ｐ４．５，Ｐ５．５，Ｐ６．５，Ｐ７．５，Ｐ９．５，Ｐ１０．５，Ｐ１１．５である。

【0042】

次にそれらの新たな測定点についての音圧が算出され（ステップＳ０８）。さらに、それらの新たな測定点を含む隣接する測定点の音圧どうしの差分からなる差分音圧が算出され（ステップＳ０９）、閾値と比較される（ステップＳ１０）。

【0043】

図６（Ｂ）には、１回目の計算における差分音圧（図６（Ａ））に重ねて、この２回目の計算における差分音圧が示されている。

【0044】

次にステップＳ０６に戻り、ここでは、図６（Ｂ）に示すように、未だ閾値を越える差分音圧が存在するため、さらに新たな測定点（図５に示す例では、Ｐ５．２５，Ｐ６．２５）が設定され（ステップＳ０７）、それらの新たな測定点（Ｐ５．２５，Ｐ６．２５）の音圧が算出され（ステップＳ０８）、さらに、それらの新たな測定点（Ｐ５．２５，Ｐ６．２５）を含む隣接する測定点の音圧どうしの差分からなる差分音圧が算出されて（ステップＳ０９）、閾値と比較される（ステップＳ１０）。

【0045】

図６（Ｃ）には、１回目および２回目の計算における差分音圧に重ねて、この３回目の計算による差分音圧が示されている。

【0046】

閾値と比較した結果（ステップＳ１０）、今度は閾値を越える差分音圧は存在しなくなり（ステップＳ０６）、この音源探査を終了する。

【0047】

ここでは、簡単に示すため１次元的な演算について説明したが、例えば２次元的に広がる関心領域内の互いに直交するｘ方向とｙ方向とのそれぞれについて上記の演算が行なわれて、その関心領域内の音圧分布が算出される。

【0048】

図７は、２次元的な関心領域内の音圧分布を可視化して示した図である。

【0049】

この図７におけるハッチングの違いは、例えば色の違いとして表現される。上記のようにして算出した音圧は、例えばこの図７に示すように表示することにより、ユーザにその情報が提供される。この図７には、音源を表わす音圧のピークＰＫが３箇所にあらわれている。図２に示すプローブ１０にはカメラ１１が備えられており、関心領域の画像を重ね合わせて表示することで、音源位置の情報が一層分かり易く提供される。

【0050】

ここでは、上記の通り、差分音圧が大きい領域についてのみ細かな間隔で音圧を算出しているため、演算量の低減、すなわちリアルタイム性の向上と、音源位置固定の正確性が高い次元で両立する。

【0051】

図８は、関心領域内の差分音圧計算の他の例を示した図である。

【0052】

ここに示す例では、前述のようにして２次元的に離散している各測定点Ｐ（図１参照）についての音圧を算出した後、以下に説明するようにして差分音圧が算出される。

【0053】

図８（Ａ）は、ｘ方向の差分演算を行なうフィルタを示す図である。

【0054】

関心領域内の任意の測定点（ｘ，ｙ）において、その測定点と、その測定点に対しｘ方向、ｙ方向の前後に隣接する各測定点の音圧に図８（Ａ）の各係数を掛けてその合計値を算出する。ここでは、この合計値をＬｘとする。

【0055】

図８（Ｂ）は、ｙ方向の差分演算を行なうフィルタを示す図である。

【0056】

関心領域内の任意の測定点（ｘ，ｙ）において、その測定点と、その測定点に対しｘ方向、ｙ方向の前後に隣接する各測定点の音圧に図８（Ｂ）の係数を掛けて、その合計値を算出する。ここでは、この合計値をＬｙとする。

【0057】

次いで、差分音圧Δとして、

【0058】

【数1】

【0059】

を計算する（図８（Ｃ））。この差分音圧Δと閾値との大小を比較して、閾値を越える差分音圧Δが得られたときは、さらにｘ方向、ｙ方向双方に２分割することにより新たな測定点を設定して、その新たな測定点について音圧を計算する。この場合、１回につき新たな測定点は４点設定されることになる。

【0060】

ここでは、測定点を矩形の升目で表現しており、１回目の測定点をＰ_{ｉ，ｊ（ｉ，ｊ＝０，１，２）}とし、新たに設定された４つの測定点をＰ_ｘ，ｙ、Ｐ_{ｘ＋１，ｙ、}Ｐ_{ｘ，ｙ＋１}、Ｐ_{ｘ＋１，ｙ＋１}とする。これら４つの新たな測定点Ｐ_ｘ，ｙ、Ｐ_{ｘ＋１，ｙ}、Ｐ_{ｘ，ｙ＋１}、Ｐ_{ｘ＋１，ｙ＋１に}ついて音圧が計算される。

【0061】

上記の図８（Ａ），（Ｂ）の計算を行なうにあたっては、ｘ方向、ｙ方向に隣接する８つの測定点の音圧と、中央にある自分自身の音圧との合計９つの測定点の音圧が必要となる。

【0062】

ここでは、Ｐ_ｘ，ｙの測定点について、図８（Ａ），（Ｂ）の計算を行なうものとして説明する。このとき、測定点Ｐ_ｘ，ｙ、その周囲の８つの測定点とを含む、図８（Ｄ）に斜線を付した９つの測定点Ｐ_{ｘ−１，ｙ−１}、Ｐ_{ｘ，ｙ−１}、Ｐ_{ｘ＋１，ｙ−１}、Ｐ_{ｘ−１，ｙ}、Ｐ_ｘ，ｙ、Ｐ_{ｘ＋１，ｙ}、Ｐ_{ｘ−１，ｙ＋１}、Ｐ_{ｘ，ｙ＋１}、Ｐ_{ｘ＋１，ｙ＋１}の音圧が計算に使われる。

【0063】

ここでは、４つの新たな測定点Ｐ_ｘ，ｙ、Ｐ_{ｘ＋１，ｙ}、Ｐ_{ｘ，ｙ＋１}Ｐ_{ｘ＋１，ｙ＋１}の音圧は計算済であるが、残りの５つの測定点Ｐ_{ｘ−１，ｙ−１}、Ｐ_{ｘ，ｙ−１}、Ｐ_{ｘ＋１，ｙ−１}、Ｐ_{ｘ−１，ｙ}、Ｐ_{ｘ−１，ｙ＋１}の音圧のデータは存在しない。そこでここでは、それら５つの測定点Ｐ_{ｘ−１，ｙ−１}、Ｐ_{ｘ，ｙ−１}、Ｐ_{ｘ＋１，ｙ−１}、Ｐ_{ｘ−１，ｙ}、Ｐ_{ｘ−１，ｙ＋１}の音圧として、それぞれ、１回目の測定点Ｐ_０，０、Ｐ_1，０、Ｐ_１，０、Ｐ_０，１、Ｐ_０，１の音圧が計算に使われる。このような音圧データの代用を行なっても、Ｐ_{ｘ−１，ｙ−１}、Ｐ_{ｘ，ｙ−１}、Ｐ_{ｘ＋１，ｙ−１}、Ｐ_{ｘ−１，ｙ}、Ｐ_{ｘ−１，ｙ＋１}の音圧を計算してその計算した音圧を使った場合と最終的にほとんど同じ精度の計算結果が得られている。

【0064】

本発明では、このように２次元的な差分演算を採用して差分音圧を算出することによって新たな測定点を設定してもよい。

【0065】

以上の通り、本実施形態によれば、音源探査のリアルタイム性の確保と音源位置の正確性を高度に両立させることができる。

【0066】

尚、上述の説明では、関心周波数ｆ_０あるいは、最高周波数ｆ_２の音の波長の１／２の間隔となるように、初期の測定点を設定する旨、説明したが、必ずしも１／２の間隔に設定する必要はなく、それよりも狭い間隔に設定してもよい。

【符号の説明】

【0067】

１ 音源探査装置
１０ プローブ
１１ カメラ
１２，１２ａ，１２ｂ，１２ｃ，１２ｄ マイクロホン
２０ 演算部

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】