JP2015173391A - 演算装置及び演算方法、並びにコンピュータプログラム及び記録媒体 - Google Patents

演算装置及び演算方法、並びにコンピュータプログラム及び記録媒体 Download PDF

Info

Publication number
JP2015173391A
JP2015173391A JP2014048817A JP2014048817A JP2015173391A JP 2015173391 A JP2015173391 A JP 2015173391A JP 2014048817 A JP2014048817 A JP 2014048817A JP 2014048817 A JP2014048817 A JP 2014048817A JP 2015173391 A JP2015173391 A JP 2015173391A
Authority
JP
Japan
Prior art keywords
sound
calculation
peripheral
detected
propagation direction
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
JP2014048817A
Other languages
English (en)
Inventor
慎一 菊池
Shinichi Kikuchi
慎一 菊池
快友 今西
Yasutomo Imanishi
快友 今西
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Pioneer Corp
Original Assignee
Pioneer Electronic Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Pioneer Electronic Corp filed Critical Pioneer Electronic Corp
Priority to JP2014048817A priority Critical patent/JP2015173391A/ja
Publication of JP2015173391A publication Critical patent/JP2015173391A/ja
Pending legal-status Critical Current

Links

Landscapes

  • Circuit For Audible Band Transducer (AREA)

Abstract

【課題】音響空間内に放射された音の流れを好適に検出する。
【解決手段】演算装置(1)は、音源(300)からの音を少なくとも2箇所(210,220)で検出する音検出手段(110)と、少なくとも2箇所の各々で音を検出した時間(T0,T1)に基づいて、音が音源から伝播する伝播方向に関連する演算を実行する演算手段(120)とを備える。このような演算装置によれば、音響空間内に放射された音の流れを、直感的且つ簡易的に知ることが可能となる。
【選択図】図6

Description

本発明は、例えば音響空間内に放射された音の流れを視覚的に把握するための演算を行う演算装置及び演算方法、並びにコンピュータプログラム及び記録媒体の技術分野に関する。
音の定位方向を知る手段として、音響空間内に放射された音の流れを視覚的に把握するという技術が知られている。例えば特許文献1及び2では、複数の検出点における音圧や粒子速度から放射エネルギー方向を演算するという技術が開示されている。
特開2005−236636号公報 特許第5181865号
しかしながら、上述した特許文献1及び2に記載された技術では、音を検出するマイクロフォン等の装置に高い感度が要求されてしまう。このため、音の流れを検出するためには、例えば専用のツールを用いることが前提となってしまい、例えばコストの増大等の技術的問題点が生ずる。
本発明が解決しようとする課題には上記のようなものが一例として挙げられる。本発明は、音響空間内に放射された音の流れを好適に検出することが可能な演算装置及び演算方法、並びにコンピュータプログラム及び記録媒体を提供することを課題とする。
上記課題を解決する演算装置は、音源からの音を少なくとも2箇所で検出する音検出手段と、前記少なくとも2箇所の各々で前記音を検出した時間に基づいて、前記音が前記音源から伝播する伝播方向に関連する演算を実行する演算手段とを備える。
上記課題を解決する演算方法は、音源からの音を少なくとも2箇所で検出する音検出工程と、前記少なくとも2箇所の各々で前記音を検出した時間に基づいて、前記音が前記音源から伝播する伝播方向に関連する演算を実行する演算工程とを備える。
上記課題を解決するコンピュータプログラムは、音源からの音を少なくとも2箇所で検出する音検出工程と、前記少なくとも2箇所の各々で前記音を検出した時間に基づいて、前記音が前記音源から伝播する伝播方向に関連する演算を実行する演算工程とをコンピュータに実行させる。
上記課題を解決する記録媒体は、上述したコンピュータプログラムが記録されている。
第1実施例に係る演算装置の全体構成を示す概略図である。 第1実施例に係る音源と基準マイク及び周辺マイクの位置関係を示す概念図(その1)である。 第1実施例に係る基準マイク及び周辺マイクで検出された信号を示すグラフ(その1)である。 第1実施例に係る音源と基準マイク及び周辺マイクの位置関係を示す概念図(その2)である。 第1実施例に係る基準マイク及び周辺マイクで検出された信号を示すグラフ(その2)である。 第2実施例に係る音源と基準マイク及び周辺マイクの位置関係を示す概念図である。 第2実施例に係る基準マイク及び周辺マイクで検出された信号を示すグラフである。 各周辺マイクに対応して演算される時間差ベクトルを示す概念図である。 各周辺マイクに対応して演算される時間差ベクトルの合成を示す概念図である。 表示部における表示例を示す平面図である。 第3実施例に係る畳み込み演算を示す概念図である。 第4実施例に係る平均値の演算方法を示す概念図である。
<1>
本実施形態に係る演算装置は、音源からの音を少なくとも2箇所で検出する音検出手段と、前記少なくとも2箇所の各々で前記音を検出した時間に基づいて、前記音が前記音源から伝播する伝播方向に関連する演算を実行する演算手段とを備える。
本実施形態の演算装置によれば、その動作時には、例えばスピーカ等の音源から放射される音が、音検出手段において検出される。本実施形態に係る音検出手段は、少なくとも2箇所において音を検出することが可能に構成されている。具体的には、例えば異なる位置に配置された2つ以上の音検出器(例えば、マイクロフォン等)から得られた音を示す信号を夫々取得可能に構成されている。
音検出手段で検出された音は、演算手段における演算に用いられる。本実施形態に係る演算手段は、少なくとも2箇所で音を検出した時間に基づいて、音が前記音源から伝播する伝播方向に関連する演算を実行する。なお、ここでの「伝播方向に関連する演算」とは、音の伝播方向を直接的に示すデータの演算を意味するだけでなく、音の伝播方向を含むデータや検出位置における音圧や位相等の演算を含む趣旨である。演算手段は、例えば一の箇所で音が検出された時間と、他の箇所で音が検出された時間との差分に基づいて、音源から放射された音がどのような方向に伝播しているのかを演算する。
上述した構成によれば、相異なる箇所で音が検出された時間を利用して伝播方向を演算できる(言い換えれば、音が各箇所で検出された時間さえ正確に知ることができれば伝播方向を演算できる)ため、例えば音圧や粒子速度を利用して伝播方向を検出する場合と比べて、音検出器に高い感度が求められない。よって、高価な音検出器を利用せずとも、簡単な処理で精度良く音の伝播方向を知ることができる。
以上説明したように、本実施形態に係る演算装置によれば、音響空間内に放射された音の流れを好適に検出することが可能である。
<2>
本実施形態に係る演算装置の一態様では、前記音検出手段は、基準位置と、該基準位置の周辺に位置する少なくとも1箇所の周辺位置とで夫々前記音を検出し、前記演算手段は、前記基準位置と前記周辺位置とで前記音を検出した時間差に基づいて、前記基準位置における前記伝播方向に関連する演算を実行する。
この態様によれば、音検出手段により、基準位置と少なくとも1箇所の周辺位置とで音が検出される。基準位置及び周辺位置で音が検出されると、演算手段では、基準位置で音が検出された時間と、周辺位置で音が検出された時間との時間差が演算される。そして、演算手段では更に、演算された時間差に基づいて、基準位置における音の伝播方向に関連する演算が実行される。
上述した構成によれば、基準位置と周辺位置との検出時間の差分に基づいて、基準位置における音の伝播方向を検出することができる。なお、周辺位置が2箇所以上設定されている場合には、基準位置で音が検出された時間との時間差は、周辺位置毎に演算されることになる。この場合には、例えば周辺位置毎の時間差ベクトルを合成して用いればよい。周辺位置の数を増加させることで、伝播方向をより正確に検出することが可能となる。
<3>
上述の如く基準位置と周辺位置とで音を検出する態様では前記音検出手段は、前記基準位置と、該基準位置と同一平面上の周辺位置である第1周辺位置及び第2周辺位置と、前記同一平面上でない周辺位置である第3周辺位置とで夫々音を検出し、前記演算手段は、前記基準位置と、前記第1周辺位置、前記第2周辺位置及び前記第3周辺位置とで前記音を検出した時間差に基づいて、前記基準位置における3次元の前記伝播方向に関連する演算を実行してもよい。
この場合、音検出手段により、基準位置と、第1、第2及び第3周辺位置とで音が検出される。第1周辺位置及び第2周辺位置は、基準位置と同一平面上に位置する周辺位置である。一方、第3周辺位置は、第1及び第2周辺位置が位置する平面と同一平面上に位置しない周辺位置である。なお、周辺位置は、第1、第2及び第3周辺位置以外に設定されていてもよく、他の周辺位置で音が検出された時間を後述の演算に利用することもできる。
基準位置、第1、第2及び第3周辺位置で音が検出されると、演算手段では、基準位置で音が検出された時間と、第1、第2及び第3周辺位置で音が検出された時間との時間差が演算される。そして、演算手段では更に、演算された時間差に基づいて、基準位置における音の伝播方向に関連する演算が実行される。ここで本態様では特に、上述したように、第1及び第2周辺位置は、基準位置と同一平面上に位置する周辺位置として、第3周辺位置は、第1及び第2周辺位置が位置する平面と同一平面上に位置しない周辺位置として設定されている。このため、演算手段は少なくとも相異なる3方向に関する時間差を演算でき、その結果として3次元での音の伝播方向に関する演算を実行することができる。
<4>
上述の如く基準位置と周辺位置とで音を検出する態様では、前記演算手段は、前記基準位置及び前記周辺位置の各々で検出した前記音を示す信号の第1ピーク又は最大値の時間差に基づいて、前記基準位置における前記伝播方向に関連する演算を実行してもよい。
この場合、第1ピーク(即ち、検出を開始してから最初に検出されるピーク)又は最大値によって基準位置と周辺位置とでの音の検出時間が比較される。従って、検出された音を示す信号に複数のピークが存在するような場合であっても、正確に検出時間の差分を演算できる。
<5>
上述の如く基準位置と周辺位置とで音を検出する態様では、前記演算手段は、前記基準位置及び前記周辺位置の各々で検出した前記音を示す信号に、所望の周波数に対応する正弦波を畳み込み演算し、前記基準位置における前記所望の周波数を有する音の前記伝播方向に関連する演算を実行してもよい。
この場合、予め所望の周波数に対応する正弦波を畳み込み演算しておくことにより、所望の周波数を有する音についての伝播方向に関連する演算が容易となる。具体的には、畳み込み演算により、音検出手段で検出された音のうち、所望の周波数を有する音に対応する信号のピークが顕著に現れる。よって、例えばピーク位置の差分を利用した伝播方向の算出が容易となる。
<6>
上述の如く所望の周波数に対応する正弦波の畳み込み演算を実行する態様では、前記演算手段は、前記基準位置及び前記周辺位置で前記音を検出した時間差の所定期間の平均値に基づいて、前記基準位置における前記伝播方向に関連する演算を実行してもよい。
この場合、音を検出した時間差は、所定期間に(即ち、複数のピークを利用して)複数回演算される。そして、複数の時間差の平均値に基づいて、伝播方向に関連する演算が実行される。このようにすれば、時間差について1つの第1ピーク又は最大値を用いる場合、空間に放射された音の過渡的な伝搬方向を演算することができるのに対し、空間に放射された音の定常的な音の伝搬方向を演算することができる。
<7>
本実施形態に係る演算装置の他の態様では、前記演算手段は、前記少なくとも2箇所の各々で前記音を検出した時間に加えて、前記少なくとも2箇所の各々で検出した前記音の音圧に基づいて、前記伝播方向に関連する演算を実行する。
この態様によれば、音が検出される時間差に加えて、検出された音の音圧(振幅)が演算に利用される。よって、音の伝播方向に加えて音の大きさを知ることが可能となり、伝播方向に関する情報を、伝播方向及び音の大きさを含む情報として算出することができる。このため、例えば伝播方向に関する情報を視覚化する場合等において、より適切な表示が可能となる。
<8>
本実施形態に係る演算装置の他の態様では、前記演算手段による演算結果に応じて、前記伝播方向を表示する表示手段を更に備える。
この態様によれば、演算手段によって演算された伝播方向を視覚的に把握することが可能となる。
<9>
上述の如く表示手段を更に備える態様では、前記表示手段は、前記伝播方向をベクトル表示してもよい。
この場合、伝播方向がベクトル表示されるため、より直感的に音の流れを把握することが可能となる。
<10>
本実施形態に係る演算方法は、音源からの音を少なくとも2箇所で検出する音検出工程と、前記少なくとも2箇所の各々で前記音を検出した時間に基づいて、前記音が前記音源から伝播する伝播方向に関連する演算を実行する演算工程とを備える。
本実施形態に係る演算方法によれば、上述した演算装置と同様に、音響空間内に放射された音の流れを好適に検出することが可能である。
なお、本実施形態に係る演算方法においても、上述した本実施形態に係る演算装置における各種態様と同様の各種態様を採ることが可能である。
<11>
本実施形態に係るコンピュータプログラムは、音源からの音を少なくとも2箇所で検出する音検出工程と、前記少なくとも2箇所の各々で前記音を検出した時間に基づいて、前記音が前記音源から伝播する伝播方向に関連する演算を実行する演算工程とをコンピュータに実行させる。
本実施形態に係るコンピュータプログラムによれば、上述した演算装置及び演算方法と同様に、音響空間内に放射された音の流れを好適に検出することが可能である。
なお、本実施形態に係るコンピュータプログラムにおいても、上述した本実施形態に係る演算装置における各種態様と同様の各種態様を採ることが可能である。
<12>
本実施形態に係る記録媒体は、上述した本実施形態に係るコンピュータプログラムが記録されている。
本実施形態に係る記録媒体によれば、記録されたコンピュータプログラムを実行することにより、音響空間内に放射された音の流れを好適に検出することが可能である。
以下では、図面を参照して本発明の実施例について詳細に説明する。
<第1実施例>
初めに、図1を参照しながら、第1実施例に係る演算装置100の構成について説明する。ここに図1は、第1実施例に係る演算装置の全体構成を示す概略図である。
図1において、第1実施例に係る演算装置100は、主な構成要素として、音検出部110及び演算部120を備えて構成されている。
音検出部110は、本発明の「音検出手段」の一例であり、接続された基準マイク210及び周辺マイク220において検出された音を示す信号を入力とする。なお、基準マイク210及び周辺マイク220は、例えば一般的なマイクロフォン等として構成されており、音源300から放射された音を検出する。音検出部110は、入力された信号(即ち、基準マイク210で検出された音を示す信号(以下、適宜「基準マイク信号」と称する)及び周辺マイク220で検出された音を示す信号(以下、適宜「周辺マイク信号」と称する))を直接、或いは各種処理を施して、演算部120に出力可能に構成されている。
演算部120は、本発明の「演算手段」の一例であり、時間差算出部121及び伝播方向算出部122を備えている。時間差算出部121は、音検出部110から入力された基準マイク信号及び周辺マイク信号を比較して、基準マイク210及び周辺マイク220において音が検出された時間差を算出し、伝播方向算出部122に出力する。伝播方向算出部122は、時間差算出部121で算出された時間差に基づいて、音源300から放射された音の伝播方向を算出する。演算部120は、算出した音の伝播方向を示すデータを、表示部400に出力可能に構成されている。
表示部400は、本発明の「表示手段」の一例であり、例えば液晶モニタ等のディスプレイとして構成されている。表示部400は、演算部120で算出された音の伝播方向を視覚的に把握可能な態様で表示可能とされている。
次に、第1実施例に係る演算装置100の動作について、図2から図5を参照して説明する。ここに図2及び図4は夫々、第1実施例に係る音源と基準マイク及び周辺マイクの位置関係を示す概念図である。また図3及び図5は夫々、第1実施例に係る基準マイク及び周辺マイクで検出された信号を示すグラフである。
図2に示すように、原点を基準にしてX方向で見た場合に、基準マイク210が周辺マイク220より原点に近い位置に配置されており、原点より後方に配置された音源300からの音がX方向の正方向に放射されている場合を考える。この場合、音源300から放射された音は、先ず基準マイク210において検出され、その後に周辺マイク220に検出されることになると考えられる。
ここで図3に示すように、音検出部110において、基準マイク信号及び周辺マイク信号が図に示すようなものとして検出されたとする。この場合、時間差算出部121は、原点に対して遠方に配置されたマイクから近傍に配置されたマイクの時間差を算出する。具体的には、原点に対して遠方に配置された周辺マイク信号の最大値が検出された時刻T1から、原点に対して近傍に配置された基準マイク信号の最大値が検出された時刻T0の時間差を算出する。なお、時間差算出部121は、最大値に代えて第1ピーク等の他のピークを比較することで時間差を算出してもよい。
図3の信号から算出された時間差は、T1−T0>0となる。この結果、基準マイク210が配置されている基準位置での音の伝播方向は、伝播方向算出部122により、図2に示す伝播方向ベクトルD1として算出される。即ち、音の伝播方向は、X方向で見た場合の正方向として算出される。
他方、図4に示すように、原点を基準にしてX方向で見た場合に、基準マイク210が周辺マイク220より原点に近い位置に配置されており、周辺マイク220より原点から離れて配置された音源300からの音がX方向の負方向に放射されている場合を考える。この場合、音源300から放射された音は、先ず周辺マイク220において検出され、その後に基準マイク210に検出されることになると考えられる。
ここで図5に示すように、音検出部110において、基準マイク信号及び周辺マイク信号が図に示すようなものとして検出されたとする。この場合、時間差算出部121は、原点に対して遠方に配置された周辺マイク信号の最大値が検出された時刻T1から、原点に対して近傍に配置された基準マイク信号の最大値が検出された時刻T0の時間差を算出する。すると算出された時間差は、T1−T0<0となる。この結果、基準マイク210が配置されている基準位置での音の伝播方向は、伝播方向算出部122により、図4に示す伝播方向ベクトルD2として算出される。即ち、音の伝播方向は、X方向で見た場合の負方向として算出される。
以上説明したように、第1実施例に係る演算装置100によれば、基準マイク210で音が検出される時間と、周辺マイク220で音が検出される時間との差分を利用して、音の伝播方向を算出することができる。第1実施例に係る演算装置では特に、音が検出されたタイミングさえ正確に検出できればよいため、例えば高感度のマイクロフォン等を使用せずとも好適に伝播方向を知ることができる。
<第2実施例>
次に、図6から図10を参照しながら、第2実施例に係る演算装置について説明する。ここに図6は、第2実施例に係る音源と基準マイク及び周辺マイクの位置関係を示す概念図である。また図7は、第2実施例に係る基準マイク及び周辺マイクで検出された信号を示すグラフであり、図8は各周辺マイクに対応して演算される時間差ベクトルを示す概念図である。更に図9は、各周辺マイクに対応して演算される時間差ベクトルの合成を示す概念図であり、図10は、表示部における表示例を示す平面図である。
なお、第2実施例は、上述した第1実施例と比べて一部の構成が異なるのみであり、その他の点については概ね同様である。このため、以下では既に説明した第1実施例と異なる部分について詳細に説明し、他の重複する部分については適宜説明を省略するものとする。
図6に示すように、第2実施例に係る演算装置100では、1つの基準マイク210に対し、複数の周辺マイク220が配置されている。具体的には、基準マイク210を囲うように、第1周辺マイク221、第2周辺マイク222、第3周辺マイク223、第4周辺マイク224が配置されている。この場合、音源300から放射された音は、先ず第1周辺マイク221において検出され、その後に基準マイク210に検出され、更にその後に第2周辺マイク222及び第3周辺マイク223に検出され、最後に第4周辺マイク224に検出されることになると考えられる。
ここで図7に示すように、音検出部110において、基準マイク信号及び各周辺マイク信号(具体的には、第1周辺マイク221で検出された音を示す第1周辺マイク信号、第2周辺マイク222で検出された音を示す第2周辺マイク信号、第3周辺マイク223で検出された音を示す第3周辺マイク信号、及び第4周辺マイク224で検出された音を示す第4周辺マイク信号)が図に示すようなものとして検出されたとする。なお、第2周辺マイク信号及び第3周辺マイク信号はほぼ同時刻に同様の信号として検出されるものであるため併せて図示している。
図8に示すように、時間差算出部121は、まず、各軸方向に対して基準マイク信号及び各周辺マイク信号の対応するピークが検出された時間差を算出する。具体的には、基準マイクに対しX方向に配置された第1周辺マイクと第4周辺マイクの時間差を算出する。基準マイク信号と第1周辺マイク信号との時間差は、X方向に対し遠方に配置された基準マイク信号から近傍に配置された第1周辺マイク信号の差分を算出する。すなわち、T0−T1>0となる。この結果、基準マイク210と第1周辺マイク221との関係での時間差ベクトルx1は、X方向での正方向のベクトルとして算出される。基準マイク信号と第4周辺マイク信号の場合、X方向に対して遠方に配置された第4周辺マイク信号から近傍に配置された基準マイク信号の差分を算出する。すなわち、T0−T4<0となる。この結果、基準マイク210と第4周辺マイク224との関係での時間差ベクトルx2は、X方向での正方向のベクトルとして算出される。次に、基準マイクに対しY方向に配置された第2周辺マイクと第3周辺マイクの時間差を算出する。基準マイク信号と第2周辺マイク信号との時間差は、Y方向に対して遠方に配置された基準マイク信号から近傍に配置された第2周辺マイク信号の差分を算出する。すなわち、T0−T2<0となる。この結果、基準マイク210と第2周辺マイク222との関係での時間差ベクトルy1は、Y方向での負方向のベクトルとして算出される。基準マイク信号と第3周辺マイク信号との時間差は、Y方向に対して遠方に配置された第3周辺マイク信号から近傍に配置された基準マイク信号の差分が算出される。すなわち、T3−T1>0となる。この結果、基準マイク210と第3周辺マイク223との関係での時間差ベクトルy2は、Y方向での正方向のベクトルとして算出される。
図9に示すように、伝播方向算出部122は、上記周辺マイク220毎に算出された時間差ベクトルx1、x2、y1及びy2を夫々合成して、基準位置での音の伝播方向ベクトルD3を算出する。具体的には、伝播方向ベクトルD3は、以下の数式(1)を用いて算出される。
D3=(x1+x2)+j(y1+y2) ・・・(1)
合成の結果、伝播方向ベクトルD3は、Y方向の成分が相殺され、X方向の正方向のベクトルとして算出される。
図10に示すように、このような伝播方向ベクトルの算出は、例えば複数箇所を基準位置として行われ、各基準位置において算出される。このようにして算出された複数の伝播方向ベクトルは、表示部400において、対応する基準位置に合わせて表示される。図に示す例では、ユーザ500の右前方にある音源300から放射される音の伝播方向が表示されている。また、各ベクトルの濃淡により音圧の強弱も表されている。具体的には、色が濃いほど音圧が強いことを示しており、音源300から遠くなるほど、音圧が弱くなっているのが分かる。このように、音の伝播方向に加えて音圧等の他のデータも併せて表示してもよい。
以上説明したように、第2実施例に係る演算装置によれば、複数の周辺マイク220を配置することにより、2次元平面上での伝播方向を知ることが可能である。ちなみに、第2実施例では、周辺マイクを同一平面上に配置したため、2次元での伝播方向が算出されている。これに対し、例えば、同一平面上ではない位置に他の周辺マイクを配置すれば、3次元での伝播方向も算出することが可能である。なお、3次元での伝播方向の算出についても、上述した方法と同様に算出できる(即ち、各周辺マイク信号との時間差ベクトルを算出し、それらを合成することで算出できる)ため、ここでの詳細な説明は省略する。
<第3実施例>
次に、図11を参照しながら、第3実施例に係る演算装置について説明する。ここに図11は、第3実施例に係る畳み込み演算を示す概念図である。
なお、第3実施例は、上述した第1実施例と比べて伝播方向の算出方法が異なるのみであり、その他の点については概ね同様である。このため、以下では既に説明した第1実施例と異なる部分について詳細に説明し、他の重複する部分については適宜説明を省略するものとする。第3実施例では、基準マイク210及び周辺マイク220の配置は、図2で示す場合と同様であるとする。
図11に示すように、第3実施例に係る演算装置100では、演算部120において、基準マイク信号及び周辺マイク信号の各々に対し、所望の周波数に対応する正弦波の畳み込み演算が行われる。このような畳み込み演算を行うことで、基準マイク信号及び周辺マイク信号の各々において、所望の周波数に対応する音のピークが顕著に現れる。よって、畳み込み演算後の基準マイク信号及び周辺マイク信号を用いて伝播方向を算出すれば、所望の周波数を有する音の伝播方向を知ることができる。具体的には、畳み込み演算後の基準マイク信号のピークが検出される時刻T0及び周辺マイク信号のピークが検出される時刻T1の時間差を利用して伝播方向ベクトルが算出される。
<第4実施例>
次に、図12を参照しながら、第4実施例に係る演算装置について説明する。ここに図12は、第4実施例に係る平均値の演算方法を示す概念図である。
なお、第4実施例は、上述した第4実施例と比べて伝播方向の算出方法が異なるのみであり、その他の点については概ね同様である。このため、以下では既に説明した第4実施例と異なる部分について詳細に説明し、他の重複する部分については適宜説明を省略するものとする。第4実施例でも、基準マイク210及び周辺マイク220の配置は、図2で示す場合と同様であるとする。
図12に示すように、第4実施例に係る演算装置100では、先ず第3実施例と同様に、基準マイク信号及び周辺マイク信号の各々に対し、所望の周波数に対応する正弦波の畳み込み演算が行われる。これにより、基準マイク信号及び周辺マイク信号の各々において、所望の周波数に対応する音のピークが顕著に現れることになる。
続いて、畳み込み演算後の基準マイク信号及び周辺マイク信号を用いて、所定期間における時間差の平均値が算出される。具体的には、先ず各ピークに対応する期間Δt毎に複数の時間差が算出される。即ち、第1ピーク同士の時間差Tは、基準マイク信号の第1ピークの検出時間T01及び周辺マイクの第1ピークの検出時間T11を用いて以下の数式(2)で算出される。
=T11−T01 ・・・(2)
同様に、第2ピーク同士の時間差Tは、基準マイク信号の第2ピークの検出時間T02及び周辺マイクの第1ピークの検出時間T12を用いて以下の数式(3)で算出される。
=T12−T02 ・・・(3)
そして、第nピーク同士の時間差Tは、基準マイク信号の第nピークの検出時間T0n及び周辺マイクの第1ピークの検出時間T1nを用いて以下の数式(4)で算出される。
=T1n−T0n ・・・(4)
これらn個の時間差T〜T算出されると、n個の時間差T〜Tの平均値Tが以下の数式(5)を用いて算出される。
=(T+T+・・・Tn−1+T)/n ・・・(5)
このようにして算出された時間差の平均値Tを用いれば、1つのピークの時間差のみを用いる場合と比べて、空間に放射された音の定常的な音の伝搬方向を演算することができる。
本発明は、上述した実施形態に限られるものではなく、特許請求の範囲及び明細書全体から読み取れる発明の要旨或いは思想に反しない範囲で適宜変更可能であり、そのような変更を伴う演算装置及び演算方法、並びにコンピュータプログラム及び記録媒体もまた本発明の技術的範囲に含まれるものである。
100 演算装置
110 音検出部
120 演算部
121 時間差算出部
122 伝播方向算出部
210 基準マイク
220 周辺マイク
300 音源
400 表示部
500 ユーザ

Claims (12)

  1. 音源からの音を少なくとも2箇所で検出する音検出手段と、
    前記少なくとも2箇所の各々で前記音を検出した時間に基づいて、前記音が前記音源から伝播する伝播方向に関連する演算を実行する演算手段と
    を備える演算装置。
  2. 前記音検出手段は、基準位置と、該基準位置の周辺に位置する少なくとも1箇所の周辺位置とで夫々前記音を検出し、
    前記演算手段は、前記基準位置と前記周辺位置とで前記音を検出した時間差に基づいて、前記基準位置における前記伝播方向に関連する演算を実行する
    ことを特徴とする請求項1に記載の演算装置。
  3. 前記音検出手段は、前記基準位置と、該基準位置と同一平面上の周辺位置である第1周辺位置及び第2周辺位置と、前記同一平面上でない周辺位置である第3周辺位置とで夫々音を検出し、
    前記演算手段は、前記基準位置と、前記第1周辺位置、前記第2周辺位置及び前記第3周辺位置とで前記音を検出した時間差に基づいて、前記基準位置における3次元の前記伝播方向に関連する演算を実行する
    ことを特徴とする請求項2に記載の演算装置。
  4. 前記演算手段は、前記基準位置及び前記周辺位置の各々で検出した前記音を示す信号の第1ピーク又は最大値の時間差に基づいて、前記基準位置における前記伝播方向に関連する演算を実行することを特徴とする請求項2又は3に記載の演算装置。
  5. 前記演算手段は、前記基準位置及び前記周辺位置の各々で検出した前記音を示す信号に、所望の周波数に対応する正弦波を畳み込み演算し、前記基準位置における前記所望の周波数を有する音の前記伝播方向に関連する演算を実行することを特徴とする請求項2から4のいずれか一項に記載の演算装置。
  6. 前記演算手段は、前記基準位置及び前記周辺位置で前記音を検出した時間差の所定期間の平均値に基づいて、前記基準位置における前記伝播方向に関連する演算を実行する
    ことを特徴とする請求項5に記載の演算装置。
  7. 前記演算手段は、前記少なくとも2箇所の各々で前記音を検出した時間に加えて、前記少なくとも2箇所の各々で検出した前記音の音圧に基づいて、前記伝播方向に関連する演算を実行することを特徴とする請求項1から6のいずれか一項に記載の演算装置。
  8. 前記演算手段による演算結果に応じて、前記伝播方向を表示する表示手段を更に備えることを特徴とする請求項1から7のいずれか一項に記載の演算装置。
  9. 前記表示手段は、前記伝播方向をベクトル表示することを特徴とする請求項8に記載の演算装置。
  10. 音源からの音を少なくとも2箇所で検出する音検出工程と、
    前記少なくとも2箇所の各々で前記音を検出した時間に基づいて、前記音が前記音源から伝播する伝播方向に関連する演算を実行する演算工程と
    を備える演算方法。
  11. 音源からの音を少なくとも2箇所で検出する音検出工程と、
    前記少なくとも2箇所の各々で前記音を検出した時間に基づいて、前記音が前記音源から伝播する伝播方向に関連する演算を実行する演算工程と
    をコンピュータに実行させることを特徴とするコンピュータプログラム。
  12. 請求項11に記載のコンピュータプログラムが記録されていることを特徴とする記録媒体。
JP2014048817A 2014-03-12 2014-03-12 演算装置及び演算方法、並びにコンピュータプログラム及び記録媒体 Pending JP2015173391A (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2014048817A JP2015173391A (ja) 2014-03-12 2014-03-12 演算装置及び演算方法、並びにコンピュータプログラム及び記録媒体

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2014048817A JP2015173391A (ja) 2014-03-12 2014-03-12 演算装置及び演算方法、並びにコンピュータプログラム及び記録媒体

Related Child Applications (1)

Application Number Title Priority Date Filing Date
JP2018212849A Division JP2019023662A (ja) 2018-11-13 2018-11-13 演算装置

Publications (1)

Publication Number Publication Date
JP2015173391A true JP2015173391A (ja) 2015-10-01

Family

ID=54260464

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2014048817A Pending JP2015173391A (ja) 2014-03-12 2014-03-12 演算装置及び演算方法、並びにコンピュータプログラム及び記録媒体

Country Status (1)

Country Link
JP (1) JP2015173391A (ja)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107450052A (zh) * 2016-05-31 2017-12-08 华为终端(东莞)有限公司 一种室内导航定位的方法、系统以及相关设备

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS59105575A (ja) * 1982-12-08 1984-06-18 Matsushita Electric Ind Co Ltd 音方向検出方式
JPH06178399A (ja) * 1992-12-11 1994-06-24 Matsushita Electric Ind Co Ltd 音像定位係数算出装置
US20120026837A1 (en) * 2010-07-28 2012-02-02 Empire Technology Development Llc Sound direction detection
JP2012129873A (ja) * 2010-12-16 2012-07-05 Chubu Electric Power Co Inc 指定領域からの伝播音の再生方法とその装置
JP2012150059A (ja) * 2011-01-20 2012-08-09 Chubu Electric Power Co Inc 音源推定方法及び音源推定装置
WO2014013924A1 (ja) * 2012-07-19 2014-01-23 日東紡音響エンジニアリング株式会社 騒音識別装置及び騒音識別方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS59105575A (ja) * 1982-12-08 1984-06-18 Matsushita Electric Ind Co Ltd 音方向検出方式
JPH06178399A (ja) * 1992-12-11 1994-06-24 Matsushita Electric Ind Co Ltd 音像定位係数算出装置
US20120026837A1 (en) * 2010-07-28 2012-02-02 Empire Technology Development Llc Sound direction detection
JP2012129873A (ja) * 2010-12-16 2012-07-05 Chubu Electric Power Co Inc 指定領域からの伝播音の再生方法とその装置
JP2012150059A (ja) * 2011-01-20 2012-08-09 Chubu Electric Power Co Inc 音源推定方法及び音源推定装置
WO2014013924A1 (ja) * 2012-07-19 2014-01-23 日東紡音響エンジニアリング株式会社 騒音識別装置及び騒音識別方法

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107450052A (zh) * 2016-05-31 2017-12-08 华为终端(东莞)有限公司 一种室内导航定位的方法、系统以及相关设备
CN107450052B (zh) * 2016-05-31 2020-01-21 华为终端有限公司 一种室内导航定位的方法、系统以及相关设备

Similar Documents

Publication Publication Date Title
JP5250089B2 (ja) 音方向検出
JP4584951B2 (ja) 音源分離装置および音源分離方法
US9552644B2 (en) Motion analysis method and apparatus
US10321234B2 (en) Signal processing device and signal processing method
JP6806816B2 (ja) 距離差−周波数分析を用いた配管の漏洩感知装置および方法
JPWO2015136823A1 (ja) 目標抽出システム、目標抽出方法、情報処理装置およびその制御方法と制御プログラム
JP2015219138A (ja) 音源探査装置、音源探査方法、および音源探査プログラム
WO2015194137A1 (ja) 位置決定装置、漏洩検知システム、位置決定方法及びコンピュータ読み取り可能記録媒体
Wu Reconstructing transient acoustic radiation from an arbitrary object with a uniform surface velocity distribution
JP2015173391A (ja) 演算装置及び演算方法、並びにコンピュータプログラム及び記録媒体
Padois et al. On the use of geometric and harmonic means with the generalized cross-correlation in the time domain to improve noise source maps
JP2019023662A (ja) 演算装置
Mimani et al. Multiple line arrays for the characterization of aeroacoustic sources using a time-reversal method
JP2022033205A (ja) 演算装置
JP2020183963A (ja) 演算装置
JP2023126485A (ja) 演算装置
JPWO2021079496A5 (ja)
US10107355B2 (en) Active damping device and design method
JP2012149916A (ja) 過渡回復電圧測定装置、過渡回復電圧測定方法及び過渡回復電圧測定プログラム
US10156493B2 (en) Position determination device, position determination system, position determination method, and computer-readable recording medium
JP6415975B2 (ja) 音源方向推定装置
JP2017173141A (ja) ソーナー装置、聴音出力方法及びプログラム
JP6933303B2 (ja) 波源方向推定装置、波源方向推定方法、およびプログラム
US20240310206A1 (en) Quicksand amount observation system, quicksand amount observation apparatus, quicksand amount observation method, and computer-readable medium
JP6689568B2 (ja) 超音波測定装置及び超音波測定方法

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20170207

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20171117

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20171219

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20180117

A02 Decision of refusal

Free format text: JAPANESE INTERMEDIATE CODE: A02

Effective date: 20180814