JP2760177B2 - Sound source location estimation method - Google Patents
Sound source location estimation methodInfo
- 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
Description
【発明の詳細な説明】DETAILED DESCRIPTION OF THE INVENTION
【０００１】[0001]
【産業上の利用分野】本発明は、複数のマイクロホンを
配置して音源位置を推定する音源位置推定方法に関する
ものである。BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to a sound source position estimating method for estimating a sound source position by arranging a plurality of microphones.
【０００２】[0002]
【従来の技術】マイクロホンアレイを用いた音源位置の
推定に際し、音源とマイクロホンの距離が比較的近距離
にある場合、音源から放射される音波を球面波とみなし
て音源位置を推定することが有効であることは文献「多
数センサによる音源位置の推定」日本音響学会誌４６巻
７号（１９９０）に詳しく述べられている。この従来の
音源位置推定方法では、音源が存在すると考えられる空
間内に複数の仮想音源位置を想定し、マイクロホンアレ
イを構成する各マイクロホンより得られる信号のスペク
トルについて先の仮想音源位置と各マイクロホンとの距
離に基づき位相と振幅を補正した後、すべてのマイクロ
ホンにおける補正後のスペクトルを加算平均する。加算
平均された平均音源スペクトルＺ（ω，ｉ）は（数１）
のように表わされる。2. Description of the Related Art When estimating the position of a sound source using a microphone array, if the distance between the sound source and the microphone is relatively short, it is effective to estimate the position of the sound source by regarding sound waves emitted from the sound source as spherical waves. This is described in detail in the document "Estimation of sound source position by multiple sensors", Journal of the Acoustical Society of Japan, Vol. 46, No. 7 (1990). In this conventional sound source position estimating method, a plurality of virtual sound source positions are assumed in a space where a sound source is considered to exist, and the spectrum of a signal obtained from each microphone constituting a microphone array is compared with the previous virtual sound source position and each microphone. After correcting the phase and the amplitude on the basis of the distances, the corrected spectra in all microphones are added and averaged. The averaged average sound source spectrum Z (ω, i) obtained by addition and averaging is (Equation 1)
It is represented as
【０００３】[0003]
【数１】 (Equation 1)
【０００４】ここで、Ｍはマイクロホンアレイを構成す
るマイクロホンの総数、Ｘm（ω）はｍ番目のマイクロ
ホンから得られるスペクトルを角周波数ωの関数として
表わしたもの、ｉはすべての仮想音源位置に順番に付け
た番号、ｒimはｉ番目の仮想音源とｍ番目のマイクロホ
ンとの距離、ｃは音速である。ただし、各マイクロホン
より得られる信号のＳ／Ｎ比は、マイクロホンと実際の
音源の距離が遠くなるほど低下すると考えられるので、
（数１）における加算平均においては、各項にｒimの反
比例項を重み関数として乗じた後、平均を行なってい
る。Here, M is the total number of microphones constituting the microphone array, X m (ω) is a spectrum obtained from the m-th microphone as a function of angular frequency ω, and i is the position of all virtual sound sources. gave the order number, r im is the distance between the i-th virtual sound source and the m-th microphone, c is the speed of sound. However, since the S / N ratio of the signal obtained from each microphone is considered to decrease as the distance between the microphone and the actual sound source increases,
In averaging in equation (1), after multiplied by the inverse term of r im as the weighting function sections are performed average.
【０００５】この平均音源スペクトルのパワー｜Ｚ
（ω，ｉ）｜は，仮想音源位置と実際の音源位置とが一
致したときに最大値をとり、この平均音源スペクトルの
パワーの最大値をもたらす仮想音源位置を実際の音源位
置と推定する。したがって、音源位置を推定するために
は、実際の音源が存在すると考えられる空間全体にわた
って、平均スペクトルのパワーを仮想音源位置の関数と
して表わす必要がある。The power of the average sound source spectrum | Z
(Ω, i) | takes the maximum value when the virtual sound source position and the actual sound source position match, and estimates the virtual sound source position that gives the maximum value of the power of the average sound source spectrum as the actual sound source position. Therefore, in order to estimate the sound source position, it is necessary to represent the power of the average spectrum as a function of the virtual sound source position over the entire space where the real sound source is considered to exist.
【０００６】[0006]
【発明が解決しようとする課題】しかしながら、上記従
来の音源位置推定方法では、音源が存在すると考えられ
る空間が広範囲に広がったり、音源位置の推定精度を高
くするために仮想音源位置の間隔を狭めたりする必要が
ある場合には、平均音源スペクトルを求めなければなら
ない仮想音源位置の数が幾何級数的に増加し、これに伴
い、音源位置を推定するための計算量が増加してしまう
という問題を有していた。However, in the above-mentioned conventional sound source position estimating method, the space where the sound source is considered to be present is widened widely, and the interval between the virtual sound source positions is narrowed in order to increase the accuracy of the sound source position estimation. The number of virtual sound source positions for which the average sound source spectrum must be obtained increases exponentially, and the amount of calculation for estimating the sound source position increases accordingly. Had.
【０００７】本発明は、上記従来の問題点を解決するも
のであり、仮想音源位置を想定する空間が拡大された
り、音源位置の推定精度を高める場合にも、計算量の極
端な増加を防止して推定効率を向上させることができる
ようにした音源位置推定方法を提供することを目的とす
る。The present invention solves the above-mentioned conventional problems, and prevents an excessive increase in the amount of calculation even when the space for assuming the virtual sound source position is expanded or the accuracy of estimating the sound source position is increased. It is an object of the present invention to provide a sound source position estimating method capable of improving the estimation efficiency.
【０００８】[0008]
【課題を解決するための手段】この目的を達成するため
に本発明による音源位置推定方法は、音源から発する信
号を空間内に配置された複数のマイクロホンで受信し、
この信号のスペクトル成分を各マイクロホンについて求
め、音源が存在すると想定される空間内に仮想した音源
位置とマイクロホンの位置との間の距離に基づき、各マ
イクロホンで受信した信号から求めたスペクトルの位相
と振幅を補正することによって音源のスペクトルを求
め、この音源スペクトルをマイクロホンについての加算
平均した平均音源スペクトルのパワーを仮想音源位置の
関数として計算することにより、真の音源位置を推定す
るに際し、仮想音源の存在想定空間を段階的に絞り込
み、かつ仮想音源位置の配置間隔を段階的に狭めながら
繰り返して音源位置を推定するようにしたものである。In order to achieve this object, a sound source position estimating method according to the present invention comprises receiving a signal emitted from a sound source by a plurality of microphones arranged in a space,
The spectrum component of this signal is obtained for each microphone, and based on the distance between the virtual sound source position and the microphone position in the space where the sound source is assumed to be present, the phase of the spectrum obtained from the signal received by each microphone is calculated. The spectrum of the sound source is obtained by correcting the amplitude, and the power of the average sound source spectrum obtained by averaging this sound source spectrum with respect to the microphone is calculated as a function of the virtual sound source position. Is narrowed down step by step, and the sound source position is repeatedly estimated while the arrangement interval of the virtual sound source position is gradually narrowed.
【０００９】[0009]
【作用】したがって、本発明によれば、１回の音源位置
推定における仮想音源位置の数を少数にとどめることが
できるので、平均音源スペクトルを求める計算回数を減
少させることができ、音源位置を推定する空間が拡大さ
れたり、仮想音源位置の間隔が狭められたりする場合に
おいても、音源位置を推定するために必要な計算量の極
端な増加を抑えて推定の精度を維持することができる。Therefore, according to the present invention, the number of virtual sound source positions in one sound source position estimation can be reduced to a small number, so that the number of calculations for obtaining the average sound source spectrum can be reduced, and the sound source position can be estimated. Even when the space to be used is enlarged or the interval between the virtual sound source positions is narrowed, it is possible to suppress the extreme increase in the amount of calculation required for estimating the sound source position and maintain the estimation accuracy.
【００１０】[0010]
【実施例】以下、本発明の一実施例について、図面を参
照しながら説明する。An embodiment of the present invention will be described below with reference to the drawings.
【００１１】図２は本発明の一実施例における音源位置
推定方法を説明するため音源とマイクロホンアレイの配
置および仮想音源位置の存在範囲を示す図である。FIG. 2 is a diagram showing the arrangement of a sound source and a microphone array and the existence range of a virtual sound source position for explaining a sound source position estimating method in one embodiment of the present invention.
【００１２】図２において、２１は５００〜５０００Ｈ
ｚの帯域で白色なスペクトルをもつ信号を発する音源、
２２ａ〜２２ｈは音源２１の設置された平面から５０ｃ
ｍ離れた平面上の半径６０ｃｍの円周上にほぼ等間隔に
配置され、円形マイクロホンアレイを構成する８個の無
指向性のマイクロホン、２３は音源が存在すると考えら
れ、仮想音源位置を配置する１ｍ四方の２次元平面の範
囲、２４はこの２次元空間に導入した直交座標軸ｘ、ｙ
を示す（以下、この座標系の単位はｃｍとする）。この
座標系において、仮想音源位置を配置する範囲は（数
２）により表わされる。In FIG. 2, 21 is 500-5000H.
a sound source that emits a signal having a white spectrum in the z band,
22a to 22h are 50c from the plane on which the sound source 21 is installed.
Eight omnidirectional microphones 23 that are arranged at substantially equal intervals on a circumference with a radius of 60 cm on a plane separated by m and constitute a circular microphone array, 23 is considered to have a sound source, and arranges a virtual sound source position 24 is a rectangular coordinate axis x, y introduced into this two-dimensional space.
(Hereinafter, the unit of this coordinate system is cm). In this coordinate system, the range in which the virtual sound source position is arranged is represented by (Equation 2).
【００１３】[0013]
【数２】 (Equation 2)
【００１４】また、実際の音源位置は（０、０）と表わ
される。本実施例においては、各マイクロホン２２ａ〜
２２ｈから得られる信号は、サンプリング周波数１２ｋ
Ｈｚで同時サンプリングされた後、５１２ポイントのサ
ンプリングデータを用いてスペクトル成分Ｘm（ω）に
変換するものとする。ここで得られたＸm（ω）をもと
に、（数１）により平均音源スペクトルＺ（ω，ｉ）を
求め、後述のような周波数成分についてこのパワー｜Ｚ
（ω，ｉ）｜を加算平均し、この最大値を与える仮想音
源位置を推定音源位置とする。本実施例においては、
（数２）の範囲に、ｘ軸方向、Ｙ軸方向ともに２ｃｍ間
隔の格子点を想定し、これを仮想音源位置とする。The actual sound source position is expressed as (0, 0). In the present embodiment, each of the microphones 22a to 22a to
The signal obtained from 22h has a sampling frequency of 12k
After simultaneous sampling at Hz, it is assumed that the data is converted into a spectral component X m (ω) using 512 points of sampling data. On the basis of X m (ω) obtained here, an average sound source spectrum Z (ω, i) is obtained by (Equation 1), and this power | Z
(Ω, i) | is added and averaged, and the virtual sound source position giving the maximum value is set as the estimated sound source position. In this embodiment,
In the range of (Equation 2), grid points at 2 cm intervals in both the x-axis direction and the Y-axis direction are assumed, and these are defined as virtual sound source positions.
【００１５】以下、本発明実施例と上記従来例の差異を
理解しやすくするため、両者を対比しながら説明する。Hereinafter, in order to make it easier to understand the difference between the embodiment of the present invention and the above-mentioned conventional example, the description will be made while comparing both examples.
【００１６】従来の音源位置推定方法では、図４の仮想
音源位置配置図に示すすべての格子点４１について平均
音源スペクトルＺ（ω，ｉ）を求める必要がある。本実
施例と同じ状況において、従来の音源位置推定方法によ
って音源位置の推定を行なおうとすると、２６０１個
（５１個×５１個）の格子点４１について平均音源スペ
クトルＺ（ω，ｉ）を求めることになる。図５は上記従
来例のように２ｃｍ間隔の２６０１個の格子点４１のす
べてについて、平均音源スペクトルＺ（ω，ｉ）を求
め、１０００〜１５００Ｈｚの周波数成分について加算
平均した結果を示す分布図である。これより音源位置が
（２、０）と推定され、実際の音源位置とよく一致して
いることがわかる。In the conventional sound source position estimating method, it is necessary to obtain an average sound source spectrum Z (ω, i) for all grid points 41 shown in the virtual sound source position arrangement diagram of FIG. In the same situation as in the present embodiment, when trying to estimate the sound source position by the conventional sound source position estimation method, the average sound source spectrum Z (ω, i) is obtained for 2601 (51 × 51) grid points 41. Will be. FIG. 5 is a distribution diagram showing a result of obtaining an average sound source spectrum Z (ω, i) for all 2601 lattice points 41 at 2 cm intervals as in the above-described conventional example, and averaging the frequency components of 1000 to 1500 Hz. is there. From this, the sound source position is estimated to be (2, 0), and it can be seen that the sound source position matches the actual sound source position well.
【００１７】一方、本発明実施例では、１回の位置推定
で平均音源スペクトルを求める仮想音源位置の数を少数
にとどめ、音源位置の推定を繰り返し、仮想音源位置を
想定する空間を段階的に狭めていく。On the other hand, in the embodiment of the present invention, the number of virtual sound source positions for which the average sound source spectrum is obtained by one position estimation is kept to a small number, the estimation of the sound source positions is repeated, and the space in which the virtual sound source positions are assumed is stepwise. Narrow down.
【００１８】図１（ａ）は本発明実施例における１回目
の音源位置推定の際の仮想音源位置の配置図、図１
（ｂ）は本発明実施例における２回目の音源位置推定の
際の仮想音源位置の配置図である。まず、（数２）の範
囲に図１（ａ）に示すように、１０ｃｍ間隔の１２１
（１１個×１１個）個の格子点１を想定し、これを初回
の音源位置推定における仮想音源位置とし、加算平均ス
ペクトルＺ（ω，ｉ）を求める。これを、５００〜１０
００Ｈｚの周波数成分について加算平均した結果が図３
（ａ）であり、これより最大の平均加算スペクトルのパ
ワーを与える点が推定音源位置として（０、０）と求ま
る。次に、図１（ｂ）に示すように、２回目の音源位置
推定のための仮想音源位置の範囲を、初回の音源位置推
定で求めた音源位置１１を中心に、１回目の音源位置推
定において隣接した８つの格子点１２ａ〜１２ｈで囲ま
れた空間（数３）となるように設定する。FIG. 1A is a layout diagram of virtual sound source positions at the time of the first sound source position estimation in the embodiment of the present invention.
(B) is an arrangement diagram of virtual sound source positions at the time of the second sound source position estimation in the embodiment of the present invention. First, as shown in FIG.
Assuming (11 × 11) lattice points 1 as virtual sound source positions in the initial sound source position estimation, an averaged spectrum Z (ω, i) is obtained. This is
FIG. 3 shows the result of averaging the frequency components of 00 Hz.
(A), and the point at which the maximum average addition spectrum power is given is determined as (0, 0) as the estimated sound source position. Next, as shown in FIG. 1B, the range of the virtual sound source position for the second sound source position estimation is centered on the sound source position 11 obtained by the first sound source position estimation. Is set to be a space (Equation 3) surrounded by eight adjacent grid points 12a to 12h.
【００１９】[0019]
【数３】 (Equation 3)
【００２０】この範囲に２ｃｍ間隔の１２１個（１１個
×１１個）の格子点１３を想定し、これを２回目の音源
位置推定における仮想音源位置とし、平均音源スペクト
ルＺ（ω，ｉ）を求めた後、２０００〜２５００Ｈｚの
周波数成分について加算平均すると、図３（ｂ）に示す
ように、スペクトルのパワーの最大値を与える仮想音源
位置が（２、０）と求まり、実際の音源位置とよく一致
していることがわかる。つまり、２回の音源位置推定を
繰り返すことにより、合計２４２個の仮想音源位置につ
いてのみ平均音源スペクトルを求めるだけで、合計２６
０１個の仮想音源位置について平均音源スペクトルを求
める場合と同様の結果が得られた。In this range, 121 (11 × 11) grid points 13 at 2 cm intervals are assumed, and these are set as virtual sound source positions in the second sound source position estimation, and the average sound source spectrum Z (ω, i) is calculated. After the calculation, the averaging of the frequency components of 2000 to 2500 Hz is performed, and as shown in FIG. 3B, the virtual sound source position giving the maximum value of the spectrum power is obtained as (2, 0), It can be seen that they match well. In other words, by repeating the sound source position estimation twice, only the average sound source spectrum is obtained for only a total of 242 virtual sound source positions, and a total of 26
The same result as in the case of obtaining the average sound source spectrum for the 01 virtual sound source positions was obtained.
【００２１】本実施例においては、２段階の音源位置推
定を行なっているが、この初回の音源位置推定におい
て、５００〜１０００Ｈｚという、２回目の推定に用い
た周波数帯域よりも低域の周波数帯域を用いているの
は、以下のような理由による。In this embodiment, two-stage sound source position estimation is performed. In this first sound source position estimation, a frequency band of 500 to 1000 Hz, which is lower than the frequency band used for the second estimation, is used. Is used for the following reasons.
【００２２】図３（ｂ）は本実施例より得られた２ｃｍ
間隔の仮想音源位置に対する平均音源スペクトルのパワ
ーを２０００〜２５００Ｈｚという高域の成分について
加算平均した結果である。図５と比較してわかるよう
に、高域のスペクトルを用いた場合には極大値の間隔が
互いに狭まっていることがわかる。したがって、１０ｃ
ｍ間隔の格子点１（図１（ａ）参照）を仮想音源位置と
する初回の音源位置推定において、高域のスペクトルを
用いてしまうと、平均音源スペクトルＺ（ω，ｉ）のパ
ワーは、仮想音源位置近傍における実音源の存在の確か
らしさを代表する値とはみなすことができなくなってし
まう。よって、本実施例のように仮想音源位置の間隔を
異ならせて段階的に音源位置推定を行なっていく場合に
は、仮想音源位置の配置間隔に対応して、平均音源スペ
クトルのうち、いずれの帯域を用いるかを決定する必要
がある。図３（ｂ）のように、真の音源位置以外に平均
音源スペクトルのパワーが極大値をもつことは、マイク
ロホン２２ａ〜２２ｈ（図２参照）の数や、用いるスペ
クトルの帯域との関係も含めて上記の文献「多数センサ
による音源位置の推定」において考察されている。FIG. 3 (b) shows a 2 cm sample obtained from this embodiment.
This is the result of averaging the power of the average sound source spectrum with respect to the virtual sound source positions at intervals for components in the high range of 2000 to 2500 Hz. As can be seen from comparison with FIG. 5, when the spectrum in the high band is used, the intervals between the maximum values are narrowed to each other. Therefore, 10c
If the spectrum of the high frequency band is used in the initial sound source position estimation using the lattice points 1 at m intervals (see FIG. 1A) as the virtual sound source position, the power of the average sound source spectrum Z (ω, i) becomes This cannot be regarded as a value representing the likelihood of the existence of the real sound source near the virtual sound source position. Therefore, in the case where the sound source position estimation is performed stepwise with different intervals between the virtual sound source positions as in the present embodiment, any one of the average sound source spectra corresponding to the arrangement interval of the virtual sound source positions is used. It is necessary to decide whether to use a band. As shown in FIG. 3B, the fact that the power of the average sound source spectrum has a maximum value other than the true sound source position includes the number of microphones 22a to 22h (see FIG. 2) and the relationship with the spectrum band used. Thus, the above-mentioned document "Estimation of sound source position by multiple sensors" is considered.
【００２３】なお、本実施例においては、８つのマイク
ロホン２２ａ〜２２ｈを音源２１から５０ｃｍの高さの
半径６０ｃｍの円周上に配置しているが、マイクロホン
アレイを構成するマイクロホンの数およびその配置につ
いては任意に選定することができる。また、本実施例に
おいては、仮想音源の存在想定空間を２次元の平面とし
たが、これを１次元、あるいは３次元の空間としても適
用可能である。In the present embodiment, the eight microphones 22a to 22h are arranged on a circumference having a height of 50 cm from the sound source 21 and a radius of 60 cm. However, the number and arrangement of the microphones constituting the microphone array are described. Can be arbitrarily selected. Further, in the present embodiment, the assumed space of the virtual sound source is a two-dimensional plane, but this space can be applied as a one-dimensional or three-dimensional space.
【００２４】更に、本実施例では音源位置の推定を２回
繰り返しているが、これを２回以上繰り返すことによ
り、仮想音源の存在想定空間がより広範囲に広がった場
合や、より高精度の推定が望まれる場合に対応すること
が可能である。本発明は、このほか、その基本的技術思
想を逸脱しない範囲で種々変更することができる。Furthermore, in the present embodiment, the estimation of the sound source position is repeated twice. However, by repeating the estimation twice or more, the space where the virtual sound source is assumed to be spread becomes wider, or the estimation with higher accuracy is performed. Can be dealt with when is desired. The present invention can be variously modified without departing from the basic technical idea.
【００２５】[0025]
【発明の効果】以上説明したように本発明によれば、音
源位置推定の計算を段階的に狭められた仮想音源の存在
空間に対して繰り返し行なうことにより、推定の精度を
劣化させることなく、音源位置推定に伴う計算量を削減
させることができ、推定効率を向上させることができ
る。As described above, according to the present invention, the calculation of the sound source position estimation is repeatedly performed on the existence space of the virtual sound source which is narrowed stepwise, so that the estimation accuracy is not degraded. The amount of calculation involved in sound source position estimation can be reduced, and the estimation efficiency can be improved.
【図１】（ａ）本発明の一実施例における音源位置推定
方法に用いる１回目の音源位置推定の際の仮想音源位置
の配置図 （ｂ）同実施例における音源位置推定方法に用いる２回
目の音源位置推定の際の仮想音源位置の配置図FIG. 1A is a layout diagram of virtual sound source positions in a first sound source position estimation used in a sound source position estimation method according to an embodiment of the present invention; FIG. Diagram of virtual sound source position for sound source position estimation
【図２】本発明の一実施例における音源位置推定方法を
説明するための音源とマイクロホンアレイの配置および
仮想音源位置の存在範囲を示す図FIG. 2 is a diagram illustrating an arrangement of a sound source and a microphone array and an existing range of a virtual sound source position for describing a sound source position estimating method according to an embodiment of the present invention.
【図３】（ａ）本発明の一実施例における音源位置推定
方法により１回目の音源位置を推定する際、５００〜１
０００Ｈｚの音源成分を用いて音源位置の推定を行なっ
た場合の加算平均スペクトルのパワー分布図 （ｂ）同実施例における音源位置推定方法により２回目
の音源位置を推定する際、２０００〜２５００Ｈｚの音
源成分を用いて音源位置の推定を行なった場合の加算平
均スペクトルのパワー分布図FIG. 3A is a diagram illustrating a case where a first sound source position is estimated by a sound source position estimating method according to an embodiment of the present invention;
Power distribution diagram of averaging spectrum when sound source position is estimated using a sound source component of 000 Hz (b) When estimating the sound source position for the second time by the sound source position estimating method in the embodiment, a sound source of 2000 to 2500 Hz Power distribution diagram of averaging spectrum when sound source position is estimated using components
【図４】従来の音源位置推定方法に用いる仮想音源位置
の配置図FIG. 4 is a layout diagram of virtual sound source positions used in a conventional sound source position estimating method.
【図５】従来の音源位置推定方法による１０００〜１５
００Ｈｚの音源成分の加算平均スペクトルのパワー分布
図FIG. 5 shows 1000 to 15 obtained by a conventional sound source position estimation method.
Power distribution diagram of averaging spectrum of 00Hz sound source component
１１ １回目の音源位置推定により音源位置と推定され
た点 １２ 格子点（１回目の音源位置推定により音源位置と
推定された点１１に隣接した８つの仮想音源位置） ２１ 実際の音源 ２２ マイクロホンアレイを構成する無指向性マイクロ
ホン ２３ 仮想音源位置の想定する２次元空間の範囲 ２４ 同２次元空間に導入したｘｙ直交座標軸11 Point estimated as the sound source position by the first sound source position estimation 12 Grid point (eight virtual sound source positions adjacent to point 11 estimated as the sound source position by the first sound source position estimation) 21 Actual sound source 22 Microphone array 23 A range of a two-dimensional space assumed for a virtual sound source position 24 An xy orthogonal coordinate axis introduced in the two-dimensional space
Claims (1)
複数のマイクロホンで受信し、この信号のスペクトル成
分を各マイクロホンについて求め、音源が存在すると想
定される空間内に仮想した音源位置とマイクロホンの位
置との間の距離に基づき、各マイクロホンで受信した信
号から求めたスペクトルの位相と振幅を補正することに
よって音源のスペクトルを求め、この音源スペクトルを
マイクロホンについての加算平均した平均音源スペクト
ルのパワーを仮想音源位置の関数として計算することに
より、真の音源位置を推定するに際し、仮想音源の存在
想定空間を段階的に絞り込み、かつ仮想音源位置の配置
間隔を段階的に狭めながら繰り返して音源位置を推定す
ることを特徴とする音源位置推定方法。A signal emitted from a sound source is received by a plurality of microphones arranged in a space, and a spectral component of the signal is obtained for each microphone. A virtual sound source position and a microphone in a space where a sound source is assumed to exist are provided. The spectrum of the sound source is obtained by correcting the phase and the amplitude of the spectrum obtained from the signal received by each microphone based on the distance from the position of the microphone, and the power of the average sound source spectrum obtained by averaging the sound source spectrum with respect to the microphone Is calculated as a function of the virtual sound source position, when estimating the true sound source position, the space where the virtual sound source is supposed to exist is narrowed down stepwise, and the sound source position is repeatedly reduced while the arrangement interval of the virtual sound source positions is gradually narrowed. A sound source position estimating method characterized by estimating a sound source position.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP3230139A JP2760177B2 (en) | 1991-09-10 | 1991-09-10 | Sound source location estimation method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP3230139A JP2760177B2 (en) | 1991-09-10 | 1991-09-10 | Sound source location estimation method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH0572025A JPH0572025A (en) | 1993-03-23 |
| JP2760177B2 true JP2760177B2 (en) | 1998-05-28 |
Family
ID=16903196
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP3230139A Expired - Fee Related JP2760177B2 (en) | 1991-09-10 | 1991-09-10 | Sound source location estimation method |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP2760177B2 (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102901950A (en) * | 2012-09-20 | 2013-01-30 | 浙江工业大学 | Method for recognizing three-dimensional coordinates of sound sources via planar arrays |
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 (en) * | 2008-12-15 | 2013-08-14 | 日東紡音響エンジニアリング株式会社 | Method and system for calculating acoustic impedance |
| JP5642027B2 (en) * | 2011-07-06 | 2014-12-17 | 株式会社日立パワーソリューションズ | Abnormal sound diagnosis apparatus and abnormal sound diagnosis method |
| JP6100655B2 (en) * | 2013-09-18 | 2017-03-22 | 株式会社小野測器 | Sound source exploration device |
| DE102018110759A1 (en) * | 2018-05-04 | 2019-11-07 | Sennheiser Electronic Gmbh & Co. Kg | microphone array |
| CN115150712A (en) * | 2022-06-07 | 2022-10-04 | 中国第一汽车股份有限公司 | Vehicle-mounted microphone system and automobile |
-
1991
- 1991-09-10 JP JP3230139A patent/JP2760177B2/en not_active Expired - Fee Related
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102901950A (en) * | 2012-09-20 | 2013-01-30 | 浙江工业大学 | Method for recognizing three-dimensional coordinates of sound sources via planar arrays |
Also Published As
| Publication number | Publication date |
|---|---|
| JPH0572025A (en) | 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 | |
| JP3863306B2 (en) | Microphone array device | |
| KR101591220B1 (en) | Apparatus and method for microphone positioning based on a spatial power density | |
| EP3096539A1 (en) | Sound processing device and method, and program | |
| KR101456866B1 (en) | Method and apparatus for extracting the target sound signal from the mixed sound | |
| US8165317B2 (en) | Method and system for position detection of a sound source | |
| JP2015042989A (en) | Multi-sensor sound source localization | |
| JP2760177B2 (en) | Sound source location estimation method | |
| US8705319B2 (en) | Apparatus and method for resolving an ambiguity from a direction of arrival estimate | |
| US11632626B2 (en) | Audio encoding device and method | |
| US20200400489A1 (en) | Directional acoustic sensor, and methods of adjusting directional characteristics and attenuating acoustic signal in specific direction using the same | |
| JP5642339B2 (en) | Signal separation device and signal separation method | |
| JP4293986B2 (en) | Method and system for representing a sound field | |
| JP2020148909A (en) | Signal processor, signal processing method and program | |
| JP5060465B2 (en) | Sound collection device, sound collection method, sound collection program, recording medium | |
| KR101958337B1 (en) | The method and apparatus for estimating the direction of arrival of a signal | |
| JP6835425B2 (en) | Building height calculation device, building height calculation method, and program | |
| US11076250B2 (en) | Microphone array position estimation device, microphone array position estimation method, and program | |
| JP2007274250A (en) | Spatial angular distribution estimating method and receiver | |
| JP2017201747A (en) | Signal processing device, signal processing method, and signal processing program | |
| JP2022132839A (en) | Audio processing apparatus, audio processing system, audio processing method, and program | |
| JP4920270B2 (en) | Signal arrival direction estimation apparatus and method, signal separation apparatus and method, and computer program | |
| JP6728961B2 (en) | Received power estimation device, received power estimation method, and received power estimation program | |
| JP5713933B2 (en) | Sound source distance measuring device, acoustic direct ratio estimating device, noise removing device, method and program thereof | |
| EP4354904A1 (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 |