JPH07334177A - Sound field analyzing device - Google Patents
Sound field analyzing deviceInfo
- Publication number
- JPH07334177A JPH07334177A JP6126471A JP12647194A JPH07334177A JP H07334177 A JPH07334177 A JP H07334177A JP 6126471 A JP6126471 A JP 6126471A JP 12647194 A JP12647194 A JP 12647194A JP H07334177 A JPH07334177 A JP H07334177A
- Authority
- JP
- Japan
- Prior art keywords
- time
- sound
- calculating
- sound intensity
- short
- 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
Links
Landscapes
- Stereophonic System (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
Description
【０００１】[0001]
【産業上の利用分野】本発明は、音楽ホール，劇場等に
おいて、音場の拡散状態を把握するための音場解析装置
に関する。BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to a sound field analysis device for grasping a diffusion state of a sound field in a music hall, a theater or the like.
【０００２】[0002]
【従来の技術】以下、従来のこの種の音場解析装置、特
に、音場の指向拡散度の算出手段について説明する。図
５は従来の指向拡散度算出手段の概略を示しており、図
中、51は方向別の音の強さの測定手段で、通常パラボラ
型のマイクロホンを用い、方向別の音の強さを測定し、
これを方向別の音の強さ52として次の指向拡散度の算出
手段53へ送る。2. Description of the Related Art A conventional sound field analysis apparatus of this type, in particular, a means for calculating the directional diffusion degree of a sound field will be described below. FIG. 5 shows an outline of a conventional directional diffusion degree calculating means. In the figure, 51 is a sound intensity measuring means for each direction, which normally uses a parabolic microphone to measure the sound intensity for each direction. Measure
This is sent to the next directional diffusion degree calculating means 53 as the sound intensity 52 for each direction.
【０００３】この指向拡散度の算出手段53における算出
式は(数１)で表され、これに基づき指向拡散度54((数
１)中のｄ)が算出されることになる。The calculation formula in the directional diffusion degree calculating means 53 is expressed by (Equation 1), and the directional diffusion degree 54 (d in (Equation 1)) is calculated based on this.
【０００４】[0004]
【数１】 [Equation 1]
【０００５】[0005]
【発明が解決しようとする課題】しかしながら、この従
来の音場解析装置により求められる指向拡散度は、方向
別の音の強さのみを入力としており、時間を考慮してい
ないために、音場の拡散状態の時間変化を把握すること
ができないという問題を有していた。However, since the directional diffusion degree obtained by this conventional sound field analysis apparatus only inputs the sound intensity of each direction and does not consider time, the sound field is not considered. There was a problem that it was not possible to grasp the time change of the diffusion state of.
【０００６】本発明は、上記従来の音場解析装置の問題
点を解決するものであり、音場の拡散状態を時間を追っ
て定量的に把握することができる音場解析装置を提供す
るものである。The present invention solves the above-mentioned problems of the conventional sound field analysis apparatus, and provides a sound field analysis apparatus capable of quantitatively grasping the diffusion state of the sound field over time. is there.
【０００７】[0007]
【課題を解決するための手段】本発明は、上記目的を達
成するために、２次元的に等角度おきに、または３次元
的に等立体角度おきに測定した各方向別インパルス応答
信号を入力とし、各方向別インパルス応答信号を時間分
割する短時間窓による時間分割手段と、方向毎の各短時
間窓内の音の強さを算出する単位時間あたりの音の強さ
の算出手段と、これを入力として各時間窓毎の短時間指
向拡散度を算出する指向拡散度の算出手段と、前記算出
手段による算出結果のグラフ化手段とを備えた装置と、
及び、方向毎の各短時間窓内の音の強さを入力として、
音の強さの方向間平均を算出する手段と、その平均値を
用いた音の強さの基準化を行う手段と、グラフ化手段と
を更に備えた装置とよりなるものである。In order to achieve the above object, the present invention inputs an impulse response signal for each direction, which is measured two-dimensionally at equal angles or three-dimensionally at equal cubic angles. The time division means for time division of each direction impulse response signal by a short time window, and means for calculating the sound intensity per unit time for calculating the sound intensity in each short time window for each direction, An apparatus including a directional diffusion degree calculating means for calculating a short-time directional diffusion degree for each time window with this as an input, and a graphing means of a calculation result by the calculating means,
And, inputting the sound intensity in each short time window for each direction,
The apparatus further comprises means for calculating an average of sound intensities between directions, means for normalizing the sound intensity using the average value, and graphing means.
【０００８】[0008]
【作用】したがって、本発明によれば、短時間指向拡散
度の時間変化を解析できるほか、方向別反射音レベルの
時間変化を解析でき、音場の拡散状態を時間を追って定
量的に把握することができる。Therefore, according to the present invention, it is possible to analyze the temporal change of the directional diffusion degree for a short time and the temporal change of the reflected sound level for each direction, and to quantitatively grasp the diffused state of the sound field over time. be able to.
【０００９】[0009]
【実施例】図１は本発明の音場解析装置の一実施例を示
す概略ブロック構成図である。図１において、１は音場
解析装置で全体の具体的手段を示すもので、方向別イン
パルス応答信号の測定手段２によって測定した２次元的
に等角度おき、または３次元的に等立体角度おきの各方
向別インパルス応答信号３を入力信号とする。この入力
信号３は、前記各方向別インパルス応答信号３を時間分
割する短時間窓による時間分割手段４と、方向毎の各短
時間窓内の音の強さを算出する単位時間あたりの音の強
さの算出手段５により、各短時間窓内の音の強さとして
出力される。ここで本発明の目的である指向拡散度の算
出を行うために、前記の単位時間あたりの音の強さの算
出手段５の出力を入力として各時間窓毎の短時間指向拡
散度を算出する指向拡散度の算出手段６と、グラフ化手
段７が設けられており、これにより短時間指向拡散度の
時間変化を解析することができるのである。さらに方向
別反射音レベルの時間変化を解析可能にするため、前記
単位時間あたりの音の強さの算出手段５の出力を入力と
する音の強さの方向間平均の算出手段８と、その平均値
を用いた音の強さの基準化手段９と、グラフ化手段10と
を設け、これにより方向別反射音レベルの時間変化を解
析できるようになり、音場の拡散状態を時間を追って定
量的に把握することができるようになったのである。1 is a schematic block diagram showing an embodiment of a sound field analyzing apparatus according to the present invention. In FIG. 1, reference numeral 1 denotes a sound field analysis apparatus as a whole, which is two-dimensionally equiangularly or three-dimensionally equidistantly three-dimensionally measured by a direction-specific impulse response signal measuring means 2. The impulse response signal 3 for each direction is used as an input signal. The input signal 3 includes a time dividing means 4 by a short time window for time-dividing the impulse response signal 3 for each direction, and a sound per unit time for calculating the sound intensity in each short time window for each direction. The intensity calculating means 5 outputs the intensity of the sound in each short time window. In order to calculate the directional diffusivity which is the object of the present invention, the short-time directional diffusivity for each time window is calculated by using the output of the means 5 for calculating the sound intensity per unit time as an input. The directional diffusion degree calculating means 6 and the graphing means 7 are provided so that the temporal change of the directional diffusion degree can be analyzed. Further, in order to make it possible to analyze the temporal change of the reflected sound level for each direction, a means 8 for calculating an average of sound intensities with the output of the sound intensity calculator 5 per unit time as an input, By providing means 9 for standardizing sound intensity using an average value and means 10 for graphing, it becomes possible to analyze the time change of the reflected sound level for each direction, and the diffusion state of the sound field can be tracked over time. It became possible to grasp it quantitatively.
【００１０】次に、上記構成において、その動作を図
１，図２，図３，図４を使用して説明する。Next, the operation of the above structure will be described with reference to FIGS. 1, 2, 3 and 4.
【００１１】まず、図１において、方向別インパルス応
答信号の測定手段２によって２次元的に等角度おきに、
または３次元的に等立体角度おきに測定した各方向別イ
ンパルス応答信号３が音場解析装置１に入力される。図
４にこの方向別インパルス応答信号の例を示してある
が、その具体的な方向別インパルス応答信号の測定手段
２は、例えば特開平４−94296号公報(特願平２−212533
号)に記載されているようなグローバルマイクロホン
や、一般的な指向性マイクロホン、またはマイクロホン
アレイによる指向性受音方法を用いることが考えられ
る。これらにより得られた各方向別インパルス応答信号
３に対し、短時間窓による時間分割手段４により各方向
別インパルス応答信号３を時間分割する。時間窓の大き
さは音場の平均自由路を音波が通過する時間以下となる
ように設定すればよく、図４では20msの時間窓となって
いる。次に単位時間あたりの音の強さの算出手段５は
(数２)のようになっており、これにより方向毎の各短時
間窓内の音の強さを算出する。First, in FIG. 1, by means of the impulse response signal measuring means 2 for each direction, two-dimensionally equiangularly,
Alternatively, the impulse response signal 3 for each direction measured three-dimensionally at equal stereoscopic angles is input to the sound field analysis device 1. FIG. 4 shows an example of the impulse response signal for each direction. A concrete measuring means 2 for the impulse response signal for each direction is disclosed in, for example, Japanese Patent Application Laid-Open No. 4-94296 (Japanese Patent Application No. 2-212533).
It is conceivable to use a directional sound receiving method using a global microphone as described in No.), a general directional microphone, or a microphone array. The impulse response signal 3 for each direction thus obtained is time-divided by the time division means 4 using a short time window. The size of the time window may be set so as to be equal to or less than the time taken for the sound wave to pass through the mean free path of the sound field. In FIG. 4, the time window is 20 ms. Next, the means 5 for calculating the sound intensity per unit time is
It is as in (Equation 2), and the sound intensity in each short time window for each direction is calculated by this.
【００１２】[0012]
【数２】 [Equation 2]
【００１３】このようにして計算された各方向の短時間
窓内の音の強さＥiを入力として、（数３）に基づく指
向拡散度の算出手段６では、各時間窓毎の短時間指向拡
散度ｄを、With the sound intensities E i in the short time windows in each direction calculated in this way as input, the directional diffusion degree calculating means 6 based on (Equation 3) calculates the short time for each time window. The directional diffusion d
【００１４】[0014]
【数３】 [Equation 3]
【００１５】のように算出し、グラフ化手段７により短
時間指向拡散度の時間変化グラフ11として出力する。図
２はこの短時間指向拡散度の時間変化のグラフ11の出力
例である。It is calculated as described above, and is output as a time change graph 11 of the short-term directional diffusion degree by the graphing means 7. FIG. 2 is an output example of the graph 11 of the time variation of the short-term directional diffusion coefficient.
【００１６】短時間指向拡散度は０から1.0の値をと
り、音場が完全拡散状態に移行するに従って1.0に近づ
く。The short-time directional diffusion coefficient takes a value of 0 to 1.0, and approaches 1.0 as the sound field shifts to the completely diffused state.
【００１７】次に、単位時間あたりの音の強さの算出手
段５の出力である方向毎の短時間窓内の音の強さＥiを
入力とし、(数４)に基づく音の強さの方向間平均の算出
手段８により、音の強さＥiの方向間平均Ｅmを算出す
る。Next, the sound intensity E i in the short-time window for each direction, which is the output of the sound intensity per unit time calculating means 5, is input, and the sound intensity based on (Equation 4) is used. An inter-direction average E m of the sound intensities E i is calculated by the inter-direction average calculation means 8.
【００１８】[0018]
【数４】 [Equation 4]
【００１９】算出された方向間の音の強さの平均値Ｅm
を用い、音の強さの基準化手段９において、(数５)に基
づき全方向の音の強さを基準化する。The average value E m of the sound intensities between the calculated directions
Using, the sound intensity standardizing means 9 standardizes the sound intensity in all directions based on (Equation 5).
【００２０】[0020]
【数５】 [Equation 5]
【００２１】次にグラフ化手段10により方向別反射音レ
ベルの時間変化のグラフ12として出力する。図３は図４
に示すように２次元的に等角度で18方向のインパルス応
答信号を測定したときの方向別反射音レベルの時間変化
のグラフの出力例である。Next, the graphing means 10 outputs the graph 12 of the time-dependent change of the reflected sound level for each direction. 3 is shown in FIG.
FIG. 7 is an output example of a graph of a time-dependent change in reflected sound level for each direction when impulse response signals in 18 directions are two-dimensionally equiangularly measured as shown in FIG.
【００２２】このように、方向別反射音エネルギーを時
間ごとに算出することにより、反射音エネルギーの方向
間のばらつきを時間を追って観察することができる。な
お、この拡散が良いほどグラフが円に近くなり、完全拡
散の場合、グラフは1.0の値の真円になる。As described above, by calculating the reflected sound energy for each direction for each time, it is possible to observe the variation in the reflected sound energy between the directions over time. Note that the better this diffusion is, the closer the graph becomes to a circle, and in the case of perfect diffusion, the graph becomes a perfect circle with a value of 1.0.
【００２３】これら図２及び図３に示したグラフによっ
て音場の拡散状態を時間を追って定量的に把握すること
ができるのである。The diffusion state of the sound field can be quantitatively grasped over time by the graphs shown in FIGS.
【００２４】[0024]
【発明の効果】本発明は上記実施例から明らかなよう
に、各方向別インパルス応答信号の時間分割手段と、単
位時間あたりの音の強さの算出手段と、この算出手段に
基づき短時間指向拡散度を算出する指向拡散度の算出手
段と、グラフ化手段とにより、短時間指向拡散度の時間
変化を解析できるほか、前記単位時間あたりの音の強さ
の算出手段に基づく前記音の強さの方向間平均の算出手
段，基準化手段，グラフ化手段等を付加することによ
り、方向別反射音レベルの時間変化を解析でき、音場の
拡散状態を時間を追って定量的に把握することができる
という効果を有する。As is apparent from the above-described embodiment, the present invention provides time division means for impulse response signals for each direction, means for calculating sound intensity per unit time, and directing for a short time based on this means. The directional diffusion degree calculating means for calculating the diffusion degree and the graphing means can analyze the time variation of the directional diffusion degree for a short time, and the sound intensity based on the calculating means of the sound intensity per unit time. It is possible to analyze the time variation of the reflected sound level for each direction by adding a means for calculating the inter-direction average, a standardization means, a graphing means, etc., and to quantitatively grasp the diffusion state of the sound field over time. It has the effect that
【図１】本発明音場解析装置の一実施例を示す概略ブロ
ック構成図である。FIG. 1 is a schematic block configuration diagram showing an embodiment of a sound field analysis device of the present invention.
【図２】短時間指向拡散度の時間変化の出力例を示すグ
ラフである。FIG. 2 is a graph showing an output example of a temporal change in short-term directional diffusion rate.
【図３】方向別反射音レベルの時間変化の出力例を示す
グラフである。FIG. 3 is a graph showing an output example of a temporal change in reflected sound level for each direction.
【図４】入力する方向別インパルス応答信号の例を示す
図である。FIG. 4 is a diagram showing an example of a direction-specific impulse response signal to be input.
【図５】従来例における指向拡散度算出の手段を示す概
略図である。FIG. 5 is a schematic diagram showing a means for calculating a directional diffusion degree in a conventional example.
１…音場解析装置、 ２…方向別インパルス応答信号の
測定手段、 ３…方向別インパルス応答信号、 ４…短
時間窓による時間分割手段、 ５…単位時間あたりの音
の強さの算出手段、 ６…指向拡散度の算出手段、
７，10…グラフ化手段、 ８…音の強さの方向間平均の
算出手段、 ９…音の強さの基準化手段、11…短時間指
向拡散度の時間変化のグラフ、 12…方向別反射音レベ
ルの時間変化のグラフ。DESCRIPTION OF SYMBOLS 1 ... Sound field analysis device, 2 ... Directional impulse response signal measuring means, 3 ... Directional impulse response signal, 4 ... Short time window time division means, 5 ... Sound intensity per unit time calculation means, 6 ... Means for calculating directional diffusion degree,
7, 10 ... Graphing means, 8 ... Sound intensity direction average calculation means, 9 ... Sound intensity standardizing means, 11 ... Short-time directional diffusion degree time change graph, 12 ... Direction-dependent A graph of the change over time of the reflected sound level.
───────────────────────────────────────────────────── フロントページの続き (51)Int.Cl.6 識別記号 庁内整理番号 ＦＩ 技術表示箇所 Ｈ０４Ｓ 7/00 Ｚ ─────────────────────────────────────────────────── ─── Continuation of the front page (51) Int.Cl. 6 Identification code Internal reference number FI technical display H04S 7/00 Z
Claims (2)
的に等立体角度おきに測定した各方向別インパルス応答
信号を入力とし、各方向別インパルス応答信号を時間分
割する短時間窓による時間分割手段と、方向毎の各短時
間窓内の音の強さを算出する単位時間あたりの音の強さ
の算出手段と、前記算出手段により算出された単位時間
あたりの音の強さを入力として各時間窓毎の短時間指向
拡散度を算出する指向拡散度の算出手段と、前記指向拡
散度の算出手段による算出結果のグラフ化手段とよりな
ることを特徴とする音場解析装置。1. A short time window for inputting an impulse response signal for each direction measured at two-dimensionally equiangular intervals or three-dimensionally for every three-dimensional equiangular angle and dividing the impulse response signals for each direction by time. Time division means, means for calculating the sound intensity per unit time for calculating the sound intensity in each short time window for each direction, and the sound intensity per unit time calculated by the calculating means A sound field analysis device comprising: a directional diffusion degree calculating means for calculating a short-time directional diffusion degree for each time window as an input; and a graphing means of a calculation result by the directional diffusion degree calculating means.
的に等立体角度おきに測定した各方向別インパルス応答
信号を入力とし、各方向別インパルス応答信号を時間分
割する短時間窓による時間分割手段と、方向毎の各短時
間窓内の音の強さを算出する単位時間あたりの音の強さ
の算出手段と、前記算出手段により算出された単位時間
あたりの音の強さを入力として音の強さの方向間平均を
算出する算出手段と、その平均値を用いた音の強さの基
準化手段と、前記基準化された音の強さのグラフ化手段
とよりなることを特徴とする音場解析装置。2. A short-time window for inputting impulse response signals for each direction measured at two-dimensionally equiangular intervals or for three-dimensionally equidistant angular intervals and time-dividing the impulse response signals for each direction. Time division means, means for calculating the sound intensity per unit time for calculating the sound intensity in each short time window for each direction, and the sound intensity per unit time calculated by the calculating means Comprising a calculating means for calculating the direction average of the sound intensity as an input, a sound intensity standardizing means using the average value, and a graphing means of the standardized sound intensity. Sound field analysis device characterized by.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP6126471A JPH07334177A (en) | 1994-06-08 | 1994-06-08 | Sound field analyzing device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP6126471A JPH07334177A (en) | 1994-06-08 | 1994-06-08 | Sound field analyzing device |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JPH07334177A true JPH07334177A (en) | 1995-12-22 |
Family
ID=14936047
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP6126471A Pending JPH07334177A (en) | 1994-06-08 | 1994-06-08 | Sound field analyzing device |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPH07334177A (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2013167607A (en) * | 2012-02-17 | 2013-08-29 | Nippon Telegr & Teleph Corp <Ntt> | Acoustic field diffusible indicator calculation device, method therefor, and program |
| CN105371945A (en) * | 2015-10-10 | 2016-03-02 | 同济大学 | Method and device for removing noise in hall pulse response measurement sound signals |
-
1994
- 1994-06-08 JP JP6126471A patent/JPH07334177A/en active Pending
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2013167607A (en) * | 2012-02-17 | 2013-08-29 | Nippon Telegr & Teleph Corp <Ntt> | Acoustic field diffusible indicator calculation device, method therefor, and program |
| CN105371945A (en) * | 2015-10-10 | 2016-03-02 | 同济大学 | Method and device for removing noise in hall pulse response measurement sound signals |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP1273927B1 (en) | Environmental noise monitoring | |
| Peterson | Simulating the response of multiple microphones to a single acoustic source in a reverberant room | |
| Argentieri et al. | A survey on sound source localization in robotics: From binaural to array processing methods | |
| Kotus et al. | Detection and localization of selected acoustic events in acoustic field for smart surveillance applications | |
| CN110186546A (en) | Hydrophone sensitivity free field wide band calibration method based on pink noise | |
| Nolan et al. | Isotropy in decaying reverberant sound fields | |
| Abdou et al. | Spatial information of sound fields for room‐acoustics evaluation and diagnosis | |
| Lovedee-Turner et al. | Three-dimensional reflector localisation and room geometry estimation using a spherical microphone array | |
| Morgenstern et al. | Design framework for spherical microphone and loudspeaker arrays in a multiple-input multiple-output system | |
| JPH07334177A (en) | Sound field analyzing device | |
| Bernstein et al. | Lateralization produced by envelope-based interaural temporal disparities of high-frequency, raised-sine stimuli: Empirical data and modeling | |
| Zhao et al. | Design and evaluation of a prototype system for real-time monitoring of vehicle honking | |
| EP2717263B1 (en) | Method, apparatus, and computer program product for categorical spatial analysis-synthesis on the spectrum of a multichannel audio signal | |
| Tervo et al. | Acoustic reflection path tracing using a highly directional loudspeaker | |
| US20120090394A1 (en) | System and method for acoustic characterization of solid materials | |
| Szwoch et al. | Detection of the incoming sound direction employing MEMS microphones and the DSP | |
| Farina et al. | 3D impulse response measurements on S. Maria del Fiore Church, Florence, Italy | |
| Oldham et al. | Computer applications in building and environmental acoustics | |
| Müller et al. | Measurement, Simulation and Auralization of Indoor Road Traffic Noise | |
| Götz et al. | Common-slope modeling of late reverberation in coupled rooms | |
| Hess et al. | Acoustical evaluation of virtual rooms by means of binaural activity patterns | |
| Massarani et al. | Measurement of Diffuse-Field Sound Absorption Coefficient of Materials Using the Two Microphones Method | |
| CN109270492A (en) | A kind of Selection of Regularization Parameter for big holographic distance | |
| O’Dwyer et al. | Machine learning for sound source elevation detection | |
| Günel et al. | Wavelet packet based analysis of sound fields in rooms using coincident microphone arrays |