JP2760177B2 - 音源位置推定方法 - Google Patents
音源位置推定方法Info
- Publication number
- JP2760177B2 JP2760177B2 JP3230139A JP23013991A JP2760177B2 JP 2760177 B2 JP2760177 B2 JP 2760177B2 JP 3230139 A JP3230139 A JP 3230139A JP 23013991 A JP23013991 A JP 23013991A JP 2760177 B2 JP2760177 B2 JP 2760177B2
- Authority
- JP
- Japan
- Prior art keywords
- sound source
- source position
- spectrum
- virtual
- microphone
- 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.)
- Expired - Fee Related
Links
Landscapes
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Description
【発明の詳細な説明】
【０００１】
【産業上の利用分野】本発明は、複数のマイクロホンを
配置して音源位置を推定する音源位置推定方法に関する
ものである。
配置して音源位置を推定する音源位置推定方法に関する
ものである。
【０００２】
【従来の技術】マイクロホンアレイを用いた音源位置の
推定に際し、音源とマイクロホンの距離が比較的近距離
にある場合、音源から放射される音波を球面波とみなし
て音源位置を推定することが有効であることは文献「多
数センサによる音源位置の推定」日本音響学会誌４６巻
７号（１９９０）に詳しく述べられている。この従来の
音源位置推定方法では、音源が存在すると考えられる空
間内に複数の仮想音源位置を想定し、マイクロホンアレ
イを構成する各マイクロホンより得られる信号のスペク
トルについて先の仮想音源位置と各マイクロホンとの距
離に基づき位相と振幅を補正した後、すべてのマイクロ
ホンにおける補正後のスペクトルを加算平均する。加算
平均された平均音源スペクトルＺ（ω，ｉ）は（数１）
のように表わされる。
推定に際し、音源とマイクロホンの距離が比較的近距離
にある場合、音源から放射される音波を球面波とみなし
て音源位置を推定することが有効であることは文献「多
数センサによる音源位置の推定」日本音響学会誌４６巻
７号（１９９０）に詳しく述べられている。この従来の
音源位置推定方法では、音源が存在すると考えられる空
間内に複数の仮想音源位置を想定し、マイクロホンアレ
イを構成する各マイクロホンより得られる信号のスペク
トルについて先の仮想音源位置と各マイクロホンとの距
離に基づき位相と振幅を補正した後、すべてのマイクロ
ホンにおける補正後のスペクトルを加算平均する。加算
平均された平均音源スペクトルＺ（ω，ｉ）は（数１）
のように表わされる。
【０００３】
【数１】
【０００４】ここで、Ｍはマイクロホンアレイを構成す
るマイクロホンの総数、Ｘm（ω）はｍ番目のマイクロ
ホンから得られるスペクトルを角周波数ωの関数として
表わしたもの、ｉはすべての仮想音源位置に順番に付け
た番号、ｒimはｉ番目の仮想音源とｍ番目のマイクロホ
ンとの距離、ｃは音速である。ただし、各マイクロホン
より得られる信号のＳ／Ｎ比は、マイクロホンと実際の
音源の距離が遠くなるほど低下すると考えられるので、
（数１）における加算平均においては、各項にｒimの反
比例項を重み関数として乗じた後、平均を行なってい
る。
るマイクロホンの総数、Ｘm（ω）はｍ番目のマイクロ
ホンから得られるスペクトルを角周波数ωの関数として
表わしたもの、ｉはすべての仮想音源位置に順番に付け
た番号、ｒimはｉ番目の仮想音源とｍ番目のマイクロホ
ンとの距離、ｃは音速である。ただし、各マイクロホン
より得られる信号のＳ／Ｎ比は、マイクロホンと実際の
音源の距離が遠くなるほど低下すると考えられるので、
（数１）における加算平均においては、各項にｒimの反
比例項を重み関数として乗じた後、平均を行なってい
る。
【０００５】この平均音源スペクトルのパワー｜Ｚ
（ω，ｉ）｜は，仮想音源位置と実際の音源位置とが一
致したときに最大値をとり、この平均音源スペクトルの
パワーの最大値をもたらす仮想音源位置を実際の音源位
置と推定する。したがって、音源位置を推定するために
は、実際の音源が存在すると考えられる空間全体にわた
って、平均スペクトルのパワーを仮想音源位置の関数と
して表わす必要がある。
（ω，ｉ）｜は，仮想音源位置と実際の音源位置とが一
致したときに最大値をとり、この平均音源スペクトルの
パワーの最大値をもたらす仮想音源位置を実際の音源位
置と推定する。したがって、音源位置を推定するために
は、実際の音源が存在すると考えられる空間全体にわた
って、平均スペクトルのパワーを仮想音源位置の関数と
して表わす必要がある。
【０００６】
【発明が解決しようとする課題】しかしながら、上記従
来の音源位置推定方法では、音源が存在すると考えられ
る空間が広範囲に広がったり、音源位置の推定精度を高
くするために仮想音源位置の間隔を狭めたりする必要が
ある場合には、平均音源スペクトルを求めなければなら
ない仮想音源位置の数が幾何級数的に増加し、これに伴
い、音源位置を推定するための計算量が増加してしまう
という問題を有していた。
来の音源位置推定方法では、音源が存在すると考えられ
る空間が広範囲に広がったり、音源位置の推定精度を高
くするために仮想音源位置の間隔を狭めたりする必要が
ある場合には、平均音源スペクトルを求めなければなら
ない仮想音源位置の数が幾何級数的に増加し、これに伴
い、音源位置を推定するための計算量が増加してしまう
という問題を有していた。
【０００７】本発明は、上記従来の問題点を解決するも
のであり、仮想音源位置を想定する空間が拡大された
り、音源位置の推定精度を高める場合にも、計算量の極
端な増加を防止して推定効率を向上させることができる
ようにした音源位置推定方法を提供することを目的とす
る。
のであり、仮想音源位置を想定する空間が拡大された
り、音源位置の推定精度を高める場合にも、計算量の極
端な増加を防止して推定効率を向上させることができる
ようにした音源位置推定方法を提供することを目的とす
る。
【０００８】
【課題を解決するための手段】この目的を達成するため
に本発明による音源位置推定方法は、音源から発する信
号を空間内に配置された複数のマイクロホンで受信し、
この信号のスペクトル成分を各マイクロホンについて求
め、音源が存在すると想定される空間内に仮想した音源
位置とマイクロホンの位置との間の距離に基づき、各マ
イクロホンで受信した信号から求めたスペクトルの位相
と振幅を補正することによって音源のスペクトルを求
め、この音源スペクトルをマイクロホンについての加算
平均した平均音源スペクトルのパワーを仮想音源位置の
関数として計算することにより、真の音源位置を推定す
るに際し、仮想音源の存在想定空間を段階的に絞り込
み、かつ仮想音源位置の配置間隔を段階的に狭めながら
繰り返して音源位置を推定するようにしたものである。
に本発明による音源位置推定方法は、音源から発する信
号を空間内に配置された複数のマイクロホンで受信し、
この信号のスペクトル成分を各マイクロホンについて求
め、音源が存在すると想定される空間内に仮想した音源
位置とマイクロホンの位置との間の距離に基づき、各マ
イクロホンで受信した信号から求めたスペクトルの位相
と振幅を補正することによって音源のスペクトルを求
め、この音源スペクトルをマイクロホンについての加算
平均した平均音源スペクトルのパワーを仮想音源位置の
関数として計算することにより、真の音源位置を推定す
るに際し、仮想音源の存在想定空間を段階的に絞り込
み、かつ仮想音源位置の配置間隔を段階的に狭めながら
繰り返して音源位置を推定するようにしたものである。
【０００９】
【作用】したがって、本発明によれば、１回の音源位置
推定における仮想音源位置の数を少数にとどめることが
できるので、平均音源スペクトルを求める計算回数を減
少させることができ、音源位置を推定する空間が拡大さ
れたり、仮想音源位置の間隔が狭められたりする場合に
おいても、音源位置を推定するために必要な計算量の極
端な増加を抑えて推定の精度を維持することができる。
推定における仮想音源位置の数を少数にとどめることが
できるので、平均音源スペクトルを求める計算回数を減
少させることができ、音源位置を推定する空間が拡大さ
れたり、仮想音源位置の間隔が狭められたりする場合に
おいても、音源位置を推定するために必要な計算量の極
端な増加を抑えて推定の精度を維持することができる。
【００１０】
【実施例】以下、本発明の一実施例について、図面を参
照しながら説明する。
照しながら説明する。
【００１１】図２は本発明の一実施例における音源位置
推定方法を説明するため音源とマイクロホンアレイの配
置および仮想音源位置の存在範囲を示す図である。
推定方法を説明するため音源とマイクロホンアレイの配
置および仮想音源位置の存在範囲を示す図である。
【００１２】図２において、２１は５００〜５０００Ｈ
ｚの帯域で白色なスペクトルをもつ信号を発する音源、
２２ａ〜２２ｈは音源２１の設置された平面から５０ｃ
ｍ離れた平面上の半径６０ｃｍの円周上にほぼ等間隔に
配置され、円形マイクロホンアレイを構成する８個の無
指向性のマイクロホン、２３は音源が存在すると考えら
れ、仮想音源位置を配置する１ｍ四方の２次元平面の範
囲、２４はこの２次元空間に導入した直交座標軸ｘ、ｙ
を示す（以下、この座標系の単位はｃｍとする）。この
座標系において、仮想音源位置を配置する範囲は（数
２）により表わされる。
ｚの帯域で白色なスペクトルをもつ信号を発する音源、
２２ａ〜２２ｈは音源２１の設置された平面から５０ｃ
ｍ離れた平面上の半径６０ｃｍの円周上にほぼ等間隔に
配置され、円形マイクロホンアレイを構成する８個の無
指向性のマイクロホン、２３は音源が存在すると考えら
れ、仮想音源位置を配置する１ｍ四方の２次元平面の範
囲、２４はこの２次元空間に導入した直交座標軸ｘ、ｙ
を示す（以下、この座標系の単位はｃｍとする）。この
座標系において、仮想音源位置を配置する範囲は（数
２）により表わされる。
【００１３】
【数２】
【００１４】また、実際の音源位置は（０、０）と表わ
される。本実施例においては、各マイクロホン２２ａ〜
２２ｈから得られる信号は、サンプリング周波数１２ｋ
Ｈｚで同時サンプリングされた後、５１２ポイントのサ
ンプリングデータを用いてスペクトル成分Ｘm（ω）に
変換するものとする。ここで得られたＸm（ω）をもと
に、（数１）により平均音源スペクトルＺ（ω，ｉ）を
求め、後述のような周波数成分についてこのパワー｜Ｚ
（ω，ｉ）｜を加算平均し、この最大値を与える仮想音
源位置を推定音源位置とする。本実施例においては、
（数２）の範囲に、ｘ軸方向、Ｙ軸方向ともに２ｃｍ間
隔の格子点を想定し、これを仮想音源位置とする。
される。本実施例においては、各マイクロホン２２ａ〜
２２ｈから得られる信号は、サンプリング周波数１２ｋ
Ｈｚで同時サンプリングされた後、５１２ポイントのサ
ンプリングデータを用いてスペクトル成分Ｘm（ω）に
変換するものとする。ここで得られたＸm（ω）をもと
に、（数１）により平均音源スペクトルＺ（ω，ｉ）を
求め、後述のような周波数成分についてこのパワー｜Ｚ
（ω，ｉ）｜を加算平均し、この最大値を与える仮想音
源位置を推定音源位置とする。本実施例においては、
（数２）の範囲に、ｘ軸方向、Ｙ軸方向ともに２ｃｍ間
隔の格子点を想定し、これを仮想音源位置とする。
【００１５】以下、本発明実施例と上記従来例の差異を
理解しやすくするため、両者を対比しながら説明する。
理解しやすくするため、両者を対比しながら説明する。
【００１６】従来の音源位置推定方法では、図４の仮想
音源位置配置図に示すすべての格子点４１について平均
音源スペクトルＺ（ω，ｉ）を求める必要がある。本実
施例と同じ状況において、従来の音源位置推定方法によ
って音源位置の推定を行なおうとすると、２６０１個
（５１個×５１個）の格子点４１について平均音源スペ
クトルＺ（ω，ｉ）を求めることになる。図５は上記従
来例のように２ｃｍ間隔の２６０１個の格子点４１のす
べてについて、平均音源スペクトルＺ（ω，ｉ）を求
め、１０００〜１５００Ｈｚの周波数成分について加算
平均した結果を示す分布図である。これより音源位置が
（２、０）と推定され、実際の音源位置とよく一致して
いることがわかる。
音源位置配置図に示すすべての格子点４１について平均
音源スペクトルＺ（ω，ｉ）を求める必要がある。本実
施例と同じ状況において、従来の音源位置推定方法によ
って音源位置の推定を行なおうとすると、２６０１個
（５１個×５１個）の格子点４１について平均音源スペ
クトルＺ（ω，ｉ）を求めることになる。図５は上記従
来例のように２ｃｍ間隔の２６０１個の格子点４１のす
べてについて、平均音源スペクトルＺ（ω，ｉ）を求
め、１０００〜１５００Ｈｚの周波数成分について加算
平均した結果を示す分布図である。これより音源位置が
（２、０）と推定され、実際の音源位置とよく一致して
いることがわかる。
【００１７】一方、本発明実施例では、１回の位置推定
で平均音源スペクトルを求める仮想音源位置の数を少数
にとどめ、音源位置の推定を繰り返し、仮想音源位置を
想定する空間を段階的に狭めていく。
で平均音源スペクトルを求める仮想音源位置の数を少数
にとどめ、音源位置の推定を繰り返し、仮想音源位置を
想定する空間を段階的に狭めていく。
【００１８】図１（ａ）は本発明実施例における１回目
の音源位置推定の際の仮想音源位置の配置図、図１
（ｂ）は本発明実施例における２回目の音源位置推定の
際の仮想音源位置の配置図である。まず、（数２）の範
囲に図１（ａ）に示すように、１０ｃｍ間隔の１２１
（１１個×１１個）個の格子点１を想定し、これを初回
の音源位置推定における仮想音源位置とし、加算平均ス
ペクトルＺ（ω，ｉ）を求める。これを、５００〜１０
００Ｈｚの周波数成分について加算平均した結果が図３
（ａ）であり、これより最大の平均加算スペクトルのパ
ワーを与える点が推定音源位置として（０、０）と求ま
る。次に、図１（ｂ）に示すように、２回目の音源位置
推定のための仮想音源位置の範囲を、初回の音源位置推
定で求めた音源位置１１を中心に、１回目の音源位置推
定において隣接した８つの格子点１２ａ〜１２ｈで囲ま
れた空間（数３）となるように設定する。
の音源位置推定の際の仮想音源位置の配置図、図１
（ｂ）は本発明実施例における２回目の音源位置推定の
際の仮想音源位置の配置図である。まず、（数２）の範
囲に図１（ａ）に示すように、１０ｃｍ間隔の１２１
（１１個×１１個）個の格子点１を想定し、これを初回
の音源位置推定における仮想音源位置とし、加算平均ス
ペクトルＺ（ω，ｉ）を求める。これを、５００〜１０
００Ｈｚの周波数成分について加算平均した結果が図３
（ａ）であり、これより最大の平均加算スペクトルのパ
ワーを与える点が推定音源位置として（０、０）と求ま
る。次に、図１（ｂ）に示すように、２回目の音源位置
推定のための仮想音源位置の範囲を、初回の音源位置推
定で求めた音源位置１１を中心に、１回目の音源位置推
定において隣接した８つの格子点１２ａ〜１２ｈで囲ま
れた空間（数３）となるように設定する。
【００１９】
【数３】
【００２０】この範囲に２ｃｍ間隔の１２１個（１１個
×１１個）の格子点１３を想定し、これを２回目の音源
位置推定における仮想音源位置とし、平均音源スペクト
ルＺ（ω，ｉ）を求めた後、２０００〜２５００Ｈｚの
周波数成分について加算平均すると、図３（ｂ）に示す
ように、スペクトルのパワーの最大値を与える仮想音源
位置が（２、０）と求まり、実際の音源位置とよく一致
していることがわかる。つまり、２回の音源位置推定を
繰り返すことにより、合計２４２個の仮想音源位置につ
いてのみ平均音源スペクトルを求めるだけで、合計２６
０１個の仮想音源位置について平均音源スペクトルを求
める場合と同様の結果が得られた。
×１１個）の格子点１３を想定し、これを２回目の音源
位置推定における仮想音源位置とし、平均音源スペクト
ルＺ（ω，ｉ）を求めた後、２０００〜２５００Ｈｚの
周波数成分について加算平均すると、図３（ｂ）に示す
ように、スペクトルのパワーの最大値を与える仮想音源
位置が（２、０）と求まり、実際の音源位置とよく一致
していることがわかる。つまり、２回の音源位置推定を
繰り返すことにより、合計２４２個の仮想音源位置につ
いてのみ平均音源スペクトルを求めるだけで、合計２６
０１個の仮想音源位置について平均音源スペクトルを求
める場合と同様の結果が得られた。
【００２１】本実施例においては、２段階の音源位置推
定を行なっているが、この初回の音源位置推定におい
て、５００〜１０００Ｈｚという、２回目の推定に用い
た周波数帯域よりも低域の周波数帯域を用いているの
は、以下のような理由による。
定を行なっているが、この初回の音源位置推定におい
て、５００〜１０００Ｈｚという、２回目の推定に用い
た周波数帯域よりも低域の周波数帯域を用いているの
は、以下のような理由による。
【００２２】図３（ｂ）は本実施例より得られた２ｃｍ
間隔の仮想音源位置に対する平均音源スペクトルのパワ
ーを２０００〜２５００Ｈｚという高域の成分について
加算平均した結果である。図５と比較してわかるよう
に、高域のスペクトルを用いた場合には極大値の間隔が
互いに狭まっていることがわかる。したがって、１０ｃ
ｍ間隔の格子点１（図１（ａ）参照）を仮想音源位置と
する初回の音源位置推定において、高域のスペクトルを
用いてしまうと、平均音源スペクトルＺ（ω，ｉ）のパ
ワーは、仮想音源位置近傍における実音源の存在の確か
らしさを代表する値とはみなすことができなくなってし
まう。よって、本実施例のように仮想音源位置の間隔を
異ならせて段階的に音源位置推定を行なっていく場合に
は、仮想音源位置の配置間隔に対応して、平均音源スペ
クトルのうち、いずれの帯域を用いるかを決定する必要
がある。図３（ｂ）のように、真の音源位置以外に平均
音源スペクトルのパワーが極大値をもつことは、マイク
ロホン２２ａ〜２２ｈ（図２参照）の数や、用いるスペ
クトルの帯域との関係も含めて上記の文献「多数センサ
による音源位置の推定」において考察されている。
間隔の仮想音源位置に対する平均音源スペクトルのパワ
ーを２０００〜２５００Ｈｚという高域の成分について
加算平均した結果である。図５と比較してわかるよう
に、高域のスペクトルを用いた場合には極大値の間隔が
互いに狭まっていることがわかる。したがって、１０ｃ
ｍ間隔の格子点１（図１（ａ）参照）を仮想音源位置と
する初回の音源位置推定において、高域のスペクトルを
用いてしまうと、平均音源スペクトルＺ（ω，ｉ）のパ
ワーは、仮想音源位置近傍における実音源の存在の確か
らしさを代表する値とはみなすことができなくなってし
まう。よって、本実施例のように仮想音源位置の間隔を
異ならせて段階的に音源位置推定を行なっていく場合に
は、仮想音源位置の配置間隔に対応して、平均音源スペ
クトルのうち、いずれの帯域を用いるかを決定する必要
がある。図３（ｂ）のように、真の音源位置以外に平均
音源スペクトルのパワーが極大値をもつことは、マイク
ロホン２２ａ〜２２ｈ（図２参照）の数や、用いるスペ
クトルの帯域との関係も含めて上記の文献「多数センサ
による音源位置の推定」において考察されている。
【００２３】なお、本実施例においては、８つのマイク
ロホン２２ａ〜２２ｈを音源２１から５０ｃｍの高さの
半径６０ｃｍの円周上に配置しているが、マイクロホン
アレイを構成するマイクロホンの数およびその配置につ
いては任意に選定することができる。また、本実施例に
おいては、仮想音源の存在想定空間を２次元の平面とし
たが、これを１次元、あるいは３次元の空間としても適
用可能である。
ロホン２２ａ〜２２ｈを音源２１から５０ｃｍの高さの
半径６０ｃｍの円周上に配置しているが、マイクロホン
アレイを構成するマイクロホンの数およびその配置につ
いては任意に選定することができる。また、本実施例に
おいては、仮想音源の存在想定空間を２次元の平面とし
たが、これを１次元、あるいは３次元の空間としても適
用可能である。
【００２４】更に、本実施例では音源位置の推定を２回
繰り返しているが、これを２回以上繰り返すことによ
り、仮想音源の存在想定空間がより広範囲に広がった場
合や、より高精度の推定が望まれる場合に対応すること
が可能である。本発明は、このほか、その基本的技術思
想を逸脱しない範囲で種々変更することができる。
繰り返しているが、これを２回以上繰り返すことによ
り、仮想音源の存在想定空間がより広範囲に広がった場
合や、より高精度の推定が望まれる場合に対応すること
が可能である。本発明は、このほか、その基本的技術思
想を逸脱しない範囲で種々変更することができる。
【００２５】
【発明の効果】以上説明したように本発明によれば、音
源位置推定の計算を段階的に狭められた仮想音源の存在
空間に対して繰り返し行なうことにより、推定の精度を
劣化させることなく、音源位置推定に伴う計算量を削減
させることができ、推定効率を向上させることができ
る。
源位置推定の計算を段階的に狭められた仮想音源の存在
空間に対して繰り返し行なうことにより、推定の精度を
劣化させることなく、音源位置推定に伴う計算量を削減
させることができ、推定効率を向上させることができ
る。
【図１】（ａ）本発明の一実施例における音源位置推定
方法に用いる１回目の音源位置推定の際の仮想音源位置
の配置図 （ｂ）同実施例における音源位置推定方法に用いる２回
目の音源位置推定の際の仮想音源位置の配置図
方法に用いる１回目の音源位置推定の際の仮想音源位置
の配置図 （ｂ）同実施例における音源位置推定方法に用いる２回
目の音源位置推定の際の仮想音源位置の配置図
【図２】本発明の一実施例における音源位置推定方法を
説明するための音源とマイクロホンアレイの配置および
仮想音源位置の存在範囲を示す図
説明するための音源とマイクロホンアレイの配置および
仮想音源位置の存在範囲を示す図
【図３】（ａ）本発明の一実施例における音源位置推定
方法により１回目の音源位置を推定する際、５００〜１
０００Ｈｚの音源成分を用いて音源位置の推定を行なっ
た場合の加算平均スペクトルのパワー分布図 （ｂ）同実施例における音源位置推定方法により２回目
の音源位置を推定する際、２０００〜２５００Ｈｚの音
源成分を用いて音源位置の推定を行なった場合の加算平
均スペクトルのパワー分布図
方法により１回目の音源位置を推定する際、５００〜１
０００Ｈｚの音源成分を用いて音源位置の推定を行なっ
た場合の加算平均スペクトルのパワー分布図 （ｂ）同実施例における音源位置推定方法により２回目
の音源位置を推定する際、２０００〜２５００Ｈｚの音
源成分を用いて音源位置の推定を行なった場合の加算平
均スペクトルのパワー分布図
【図４】従来の音源位置推定方法に用いる仮想音源位置
の配置図
の配置図
【図５】従来の音源位置推定方法による１０００〜１５
００Ｈｚの音源成分の加算平均スペクトルのパワー分布
図
００Ｈｚの音源成分の加算平均スペクトルのパワー分布
図
１１ １回目の音源位置推定により音源位置と推定され
た点 １２ 格子点（１回目の音源位置推定により音源位置と
推定された点１１に隣接した８つの仮想音源位置） ２１ 実際の音源 ２２ マイクロホンアレイを構成する無指向性マイクロ
ホン ２３ 仮想音源位置の想定する２次元空間の範囲 ２４ 同２次元空間に導入したｘｙ直交座標軸
た点 １２ 格子点（１回目の音源位置推定により音源位置と
推定された点１１に隣接した８つの仮想音源位置） ２１ 実際の音源 ２２ マイクロホンアレイを構成する無指向性マイクロ
ホン ２３ 仮想音源位置の想定する２次元空間の範囲 ２４ 同２次元空間に導入したｘｙ直交座標軸
Claims (1)
- 【請求項１】音源から発する信号を空間内に配置された
複数のマイクロホンで受信し、この信号のスペクトル成
分を各マイクロホンについて求め、音源が存在すると想
定される空間内に仮想した音源位置とマイクロホンの位
置との間の距離に基づき、各マイクロホンで受信した信
号から求めたスペクトルの位相と振幅を補正することに
よって音源のスペクトルを求め、この音源スペクトルを
マイクロホンについての加算平均した平均音源スペクト
ルのパワーを仮想音源位置の関数として計算することに
より、真の音源位置を推定するに際し、仮想音源の存在
想定空間を段階的に絞り込み、かつ仮想音源位置の配置
間隔を段階的に狭めながら繰り返して音源位置を推定す
ることを特徴とする音源位置推定方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP3230139A JP2760177B2 (ja) | 1991-09-10 | 1991-09-10 | 音源位置推定方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP3230139A JP2760177B2 (ja) | 1991-09-10 | 1991-09-10 | 音源位置推定方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH0572025A JPH0572025A (ja) | 1993-03-23 |
| JP2760177B2 true JP2760177B2 (ja) | 1998-05-28 |
Family
ID=16903196
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP3230139A Expired - Fee Related JP2760177B2 (ja) | 1991-09-10 | 1991-09-10 | 音源位置推定方法 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP2760177B2 (ja) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102901950A (zh) * | 2012-09-20 | 2013-01-30 | 浙江工业大学 | 平面阵列识别声源三维坐标的方法 |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5867581A (en) * | 1994-10-14 | 1999-02-02 | Matsushita Electric Industrial Co., Ltd. | Hearing aid |
| JP5265327B2 (ja) * | 2008-12-15 | 2013-08-14 | 日東紡音響エンジニアリング株式会社 | 音響インピーダンスの計算方法及びシステム |
| JP5642027B2 (ja) * | 2011-07-06 | 2014-12-17 | 株式会社日立パワーソリューションズ | 異音診断装置および異音診断方法 |
| JP6100655B2 (ja) * | 2013-09-18 | 2017-03-22 | 株式会社小野測器 | 音源探査装置 |
| DE102018110759A1 (de) * | 2018-05-04 | 2019-11-07 | Sennheiser Electronic Gmbh & Co. Kg | Mikrofonarray |
| CN115150712B (zh) * | 2022-06-07 | 2024-06-18 | 中国第一汽车股份有限公司 | 一种车载麦克风系统及汽车 |
-
1991
- 1991-09-10 JP JP3230139A patent/JP2760177B2/ja not_active Expired - Fee Related
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102901950A (zh) * | 2012-09-20 | 2013-01-30 | 浙江工业大学 | 平面阵列识别声源三维坐标的方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| JPH0572025A (ja) | 1993-03-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US8577055B2 (en) | Sound source signal filtering apparatus based on calculated distance between microphone and sound source | |
| Huang et al. | Real-time passive source localization: A practical linear-correction least-squares approach | |
| US6760449B1 (en) | Microphone array system | |
| US8233353B2 (en) | Multi-sensor sound source localization | |
| KR101591220B1 (ko) | 공간적 전력 밀도에 기초하여 마이크 위치 결정을 위한 장치 및 방법 | |
| US20070274528A1 (en) | Acoustic Processing Device | |
| JP7235534B2 (ja) | マイクロホンアレイ位置推定装置、マイクロホンアレイ位置推定方法、およびプログラム | |
| JP2760177B2 (ja) | 音源位置推定方法 | |
| US12081949B2 (en) | Systems and methods for loudspeaker layout mapping | |
| JP5642339B2 (ja) | 信号分離装置及び信号分離方法 | |
| JP2020148909A (ja) | 信号処理装置、信号処理方法およびプログラム | |
| EP3414919A1 (en) | Microphone probe, method, system and computer program product for audio signals processing | |
| JP2013536477A (ja) | 到来方向推定値から曖昧性を解消する装置及び方法 | |
| EP1883046A1 (en) | Image interpolation method and apparatus therefor | |
| JP4293986B2 (ja) | 音場を表す方法及びシステム | |
| JP5060465B2 (ja) | 収音装置、収音方法、収音プログラム、記録媒体 | |
| JPWO2006030834A1 (ja) | 信号到来方向推定装置、信号到来方向推定方法、および信号到来方向推定用プログラム | |
| KR20190001170A (ko) | 신호의 도래각을 추정하는 방법 및 장치 | |
| JP2022132839A (ja) | 音響処理装置、音響処理システム、音響処理方法およびプログラム | |
| CN118918244A (zh) | 一种基于点云的建筑物三维重建方法及系统 | |
| JP6650245B2 (ja) | インパルス応答生成装置及びプログラム | |
| CN112116136B (zh) | 一种最短路径的生成方法及相关装置 | |
| JP2018170617A (ja) | 収音装置、プログラム及び方法 | |
| US20250220380A1 (en) | Multidimensional acoustic crosstalk cancellation filter interpolation | |
| US12192739B2 (en) | Interpolation of finite impulse response filters for generating sound fields |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| LAPS | Cancellation because of no payment of annual fees |