JPH0237971B2 - ONBAHYOKAKEISOKUKI - Google Patents

ONBAHYOKAKEISOKUKI

Info

Publication number
JPH0237971B2
JPH0237971B2 JP9297383A JP9297383A JPH0237971B2 JP H0237971 B2 JPH0237971 B2 JP H0237971B2 JP 9297383 A JP9297383 A JP 9297383A JP 9297383 A JP9297383 A JP 9297383A JP H0237971 B2 JPH0237971 B2 JP H0237971B2
Authority
JP
Japan
Prior art keywords
sound
delay time
reverberation
time
value
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 - Lifetime
Application number
JP9297383A
Other languages
Japanese (ja)
Other versions
JPS59216022A (en
Inventor
Yoichi Ando
Kenji Shima
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.)
Mitsubishi Electric Corp
Original Assignee
Mitsubishi Electric 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 Mitsubishi Electric Corp filed Critical Mitsubishi Electric Corp
Priority to JP9297383A priority Critical patent/JPH0237971B2/en
Publication of JPS59216022A publication Critical patent/JPS59216022A/en
Publication of JPH0237971B2 publication Critical patent/JPH0237971B2/en
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01HMEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
    • G01H3/00Measuring characteristics of vibrations by using a detector in a fluid

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Description

【発明の詳細な説明】[Detailed description of the invention]

この発明は音楽などを聴く音場の音響特性を評
価しやすい物理量および心理量の形に変換して出
力できる計測器に関するものである。 従来から好ましい音や音場に関する研究が行な
われているが、今回行こなつた研究により、音響
装置の再生音の基準となるコンサートホールにお
ける良い音場を決めるために必要なパラメータと
しては、次に示すようなものがあることが、次第
に明らかになつて来た。 即ち、両耳の音の特徴を表わす重要な客観的な
パラメータとしては聴取音圧・第1反射音の遅れ
時間・後続残響音の残響時間および両耳間相互相
関係数の4つの要素て有ることが、一連の模擬音
場におけるプリフアレンス(聴感上の心地良さ)
の試験によつて明らかとなつた。 次に上記4つの要素に関して詳しく説明する。 まず、第1図に反射壁の存在する空間における
音源と人頭との関係を示す。図において11は人
頭、12は音源、13,14は反射壁である。こ
こで、音源信号をp(t)、音源から左耳および右
耳へのインパルス応答をそれぞれhl(t)・h∫(t)
としたとき、左耳および右耳の信号fl(t)とfr
(t)とはそれぞれ fl(t)=∫t pp(ν)hl(t−ν)dν=p(t)*h
l
(t) ……(式1a) fr(t)=∫t pp(ν)hr(t−ν)dν=p(t)*h
r
(t) ……(式1b) と表わされる。上記式の中の*印はコンポリユー
シヨンを示している。 なお図においてはhonのn=0は直接耳に入る
場合を、n=1は反射壁13で反射した音が耳に
入る場合を、n=2は反射壁14で反射した音が
耳に入る場合を、m=lは左の耳に音が入つてく
る場合を、m=rは右の耳に音が入つてくる場合
をそれぞれ示している。 ところで第1図においては、反射壁13,14
で反射した後耳に入る音は1つの反射壁に対して
2つしか示していないが、この反射が多数起こつ
ているとして、この反射壁13,14における反
射時のインパルス応答をWo(t)とすると、左の
耳及び右の耳へ達するインパルス応答hl(t)、hr
(t)はそれぞれ hl(t)=n=0 GoWo(t−△to)*hol(t) ……(式2a) hr(t)=n=0 GnWo(t−△to)*hor(t) ……(式2b) と表わすことができる。この(式2a)と(式2b)
を用いると上記の(式1a)と(式1b)はそれぞ
れ fl(t)=n=0 p(t)*GoWo(t−△to)*hol ……(式3a) fr(t)=n=0 p(t)*GoWo(t−△to)*hor ……(式3b) と表わさせる。 ここで音源p(t)が一様な輻射特性を持たな
い場合には、各々の方向別の輻射パターンを考慮
して、p(t)をpo(t)で置換できる。 ところで両耳に入つてくる音響信号の中の情報
として、独立で客観的な音響パラメータが含まれ
ており、その第1のパラメータとして音源信号p
(t)をあげることができる。その音源信号を用
いて、その長時間の自己相関関数ΦP(τ)を Φp(τ)= lim T→∞1/2T∫T -Tp′(t)p′(t+τ)dt ……(式4) と表わすことができる。ここでp′(t)=p(t)
*s(t)であり、s(t)は耳の感度に対応し、
理論上は中耳と外耳の特性で表わされるが、実用
上は聴感特性を近似したものとしてよく知られて
いるGフイルタのインパルス応答として表わすこ
とができる。当然ながら、p′(t)のパワーΦp
(o)で(式4)を割り算することによつて正規
化自己相関関数φp(τ)は φp(τ)=Φp(τ)/Φp(o)……(
式5) と表わすことができる。 第2図a,bは上記(式5)に対応する正規化
自己相関関数の測定値を図に示したものである。
第2図aはギーボンス(Gibbons)による音楽
“ローヤルパベーン(Royal Pavane)”に対応し
た測定値を表わしたものであり、この音楽を以後
音学Aと称する。第2図bはアーノルド
(Arnold)による音楽“シンホニエツタ、作品
48:ムーブメント、アレグロコンブリオ
(Synfornietta、Opus48:Movement、
Allegro con brio)”に対応した測定値を表わし
たものであり、この音学を以後音楽Bと称する。 次に、第2の客観的なパラメータとしては壁な
どの境界における反射によつて生ずるインパルス
応答があげられる。これは直接音と第一反射音と
の間の初期時間遅れに関係するとともに、初期反
射音や後続残響音さらに反射に基づくスペクトラ
ムの変化などにも関係するものである。 第3の客観的なパラメータとしては左右の耳へ
のインパルス応答hol(t),hor(t)があげられ
る。このインパルス応答は音の定位に重要な役割
を演じるものであり、それぞれ互いに独立な関係
にはない。そのことは、中央定位する信号の場合
には、hol(t)hor(t)となることから明らか
である。 次に、この2つのインパルス応答hol(t)とhor
(t)の間の相互の従属関係を導びき出す。 まず、両耳の信号fl(t),fr(t)の間の長時
間の両耳間相互関関数Φlr(τ)をfl(t),fr(t)
を用いて表わすと、 Φlr(τ)= lim T→∞1/2T∫T -Tfl(t)fr(t +τ)dt、|τ|≦|ms ……(式6) となる。 ところで、拡がり感あるいは方向のない感覚
は、両耳間相互相関が小さい値の場合に発生する
ものであり、定まつた方向からの信号だけの場合
には|τ|<1msにおいて、両耳間相互相関関
数は大きなピークを持つ。ここで|τ|<1ms
としているのは、両耳の信号fl(t)、fr(t)の
間の時間差は、両耳間距離と音速との関係より、
通常1ms以上にはなり得ないことによるもので
ある。 まず、直接音のみによる両耳間相互相関関数
Φ(o) lr(τ)はfl(t)=p(t)*hpl(t)、fr
t)=p
(t)*hpr(t)として(式6)に代入すれば得
られる。ところで、直接音のみの正規化自己相関
関数φ(o) lr(τ)は と表わすことができる。なおφ(o) lr(τ)は仮にhpl
(t)hpr(t)ならばほぼ1になる。ここでΦ(o) ll
(o)とΦ(o) rr(o)は左と右のそれぞれの耳におけ
る信号のτ=0の自己相関関数を示している。次
に直接波の自己相関関数が小さくなる時間以後に
別の反射音が直接音に加わるものとすれば、その
ときの正規化自己相関関数φ(N) lr(τ)は、Wo(t)
がデイラツクのデルタ関数δ(t)に等しい場合
には、 と表わされる。ここで、Φ(n) lr(τ)はn番目の反
射の両耳間相互関関数をΦ(n) ll(o)、Φ(n) rr(o)
はそ
れぞれの左、右の耳のn番めの反射音のτ=0に
おける自己相関関数を示している。なお、通常の
部屋で音源が正面にある場合には、両耳間相互相
関関数の最大値はτ=0に極めて近い所で得られ
る。 ところで、正規化両耳間相互相関の強さを
IACCとして IACC=|φlr(τ)|max、for|τ|1ms
……(式9) と定義する。なお、φlr(τ)=Φlr(τ)/√ll
(o)Φrr(o)であり、Φll(o)、Φrr(o)はそ

ぞれ左右の耳の信号音のτ=0における自己相関
関数を示している。 以上において、音源信号の長時間の自己相関関
数と、音が壁などで反射することによつて生ずる
複数のインパルス応答と、両耳の信号の相関関係
を示す両耳間相互相関関数とに関して述べたが、
つぎにこれらの物理量から、両耳の音の信号の持
徴を表わす重要な客観的なパラメータである4つ
の要素即ち、聴取音圧、第1反射音の遅れ時間、
後続残響音の残響時間、および両耳間相関関数の
値のそれぞれの最適値のもとめ方について述べ
る。 前述の音楽Aおよび音楽Bを用いて行つたプリ
フアレンステストの結果を第3図a,bに示す。
横軸はIACC、縦軸は聴取音圧〔単位、dBA〕で
ある。聴取音圧は自己相関関数のτ=0の値であ
ると考えることができ、グラフの横軸の値はプリ
フアレンスを示している。第3図a,bより明ら
かに最適な聴取音圧〔p〕pはIACCにさほど依存
しておらず、ややスローテンポの音楽Aの場合に
は77〜79dBA、テンポの速い音楽Bの場合には
79〜80dBAであることがわかる。いずれの場合
も聴取音圧の最適値は79dBA前後の値を取つて
いることがわかる。 次に、音楽とスピーチを用いてスピーカ再生時
の直接音と単一反射音からなる合成音場をプリフ
アレンスで評価した結果によれば、音源信号の正
規化された自己相関関数φp(τ)を求め、反射音
のレベルを直接音の±6dBに亘つて変化させたと
き、その反射音の最適遅れ時間は、|Φp(τ)|が
第1反射音のレベルG1の1/10に相当するその時
間に対応することが明らかになつた。そこで|
Φp(τ)|が第1反射音のレベルC1の1/10に相当
する時間τdを横軸に取り、プリフアレンスが最
大となる単一反射音の遅れ時間〔△t1pを縦軸に
して表わした図が第4図である。図中に示された
範囲は、プリフアレンスの最大値より0.1低い時
の遅れ時間を示したもので、図中の記号○は第1
反射音のレベルG1=6dB、●はG1=0dB、□は
G1=−6dBの時のそれぞれを示している。とくに
|Φp(τ)がΦp(o)の0.1倍になる時間をτpとす
るとG1=0dBの場合、τd=τpと表現できる。なお
前述の音楽AとBのτpは第2図aとbからも類推
できるようにそれぞれ127msと35msである。
ここで、第4図から明らかなように、図中の各点
に近似させて直線を引くことができ、その直線に
対応する両軸を見てみるとτdはプリフアレンス
が最大となる単一反射音の遅ね時間〔△t1pにほ
ぼ一致することがわかる。同時に|Φp(τ)|が
Φp(o)の0.1倍になる時間τpともほぼ一致する。
即ち、これを数式で表現すれば、最も好ましい第
一反射音の遅延時間〔△t1pは 〔△t1〕p=τp ……(式10) と表わすことができる。 なお、τ>τpにおいて、 |φp(τ)|≦KGc ……(式11) であり、この時K=const(=0.1)、C=const(=
1.0)、
The present invention relates to a measuring instrument that can convert and output the acoustic characteristics of a sound field in which music is listened to, into physical quantities and psychological quantities that are easy to evaluate. Research has been conducted on desirable sounds and sound fields, but the research conducted this time revealed that the following parameters are necessary to determine a good sound field in a concert hall, which is the standard for the sound reproduced by acoustic equipment. It gradually became clear that there was something to show for it. In other words, there are four important objective parameters that represent the characteristics of binaural sound: audible sound pressure, delay time of first reflected sound, reverberation time of subsequent reverberant sound, and interaural cross-correlation coefficient. That is, the preference (auditory comfort) in a series of simulated sound fields.
This was clarified by a test. Next, the above four elements will be explained in detail. First, Figure 1 shows the relationship between a sound source and a person's head in a space with reflective walls. In the figure, 11 is a human head, 12 is a sound source, and 13 and 14 are reflecting walls. Here, the sound source signal is p(t), and the impulse responses from the sound source to the left ear and right ear are h l (t) and h∫(t), respectively.
Then, the left ear and right ear signals f l (t) and f r
(t) is fl(t)=∫ t p p(ν)h l (t-ν)dν=p(t)*h
l
(t) ...(Formula 1a) fr(t)=∫ t p p(ν)h r (t-ν)dν=p(t)*h
r
(t) ... (Equation 1b) It is expressed as follows. The * mark in the above formula indicates composition. In the figure, n=0 of h on means that the sound directly enters the ear, n=1 means that the sound reflected by the reflective wall 13 enters the ear, and n=2 means that the sound reflected from the reflective wall 14 enters the ear. m=l indicates the case where the sound enters the left ear, and m=r indicates the case where the sound enters the right ear. By the way, in FIG. 1, the reflecting walls 13, 14
The number of sounds that enter the ear after being reflected by the walls is shown only two for one reflecting wall, but assuming that many such reflections occur, the impulse response when reflected at the reflecting walls 13 and 14 is written as W o (t ), the impulse responses reaching the left and right ears hl(t), h r
(t) is h l (t) = n=0 G o W o (t-△t o ) * h ol (t) ... (Equation 2a) h r (t) = n=0 It can be expressed as GnW o (t−△t o )*h or (t) (Formula 2b). This (Equation 2a) and (Equation 2b)
Using (Equation 1a) and (Equation 1b) above, fl (t) = n=0 p (t) * G o W o (t-△t o ) * h ol ... (Equation 3a) ) fr(t)= n=0 p(t)*G o W o (t-△t o )*h or ... (Formula 3b). Here, if the sound source p(t) does not have uniform radiation characteristics, p(t) can be replaced by p o (t), taking into consideration the radiation pattern for each direction. By the way, the information in the acoustic signal that enters both ears includes independent and objective acoustic parameters, and the first parameter is the sound source signal p.
(t) can be given. Using the sound source signal, the long-term autocorrelation function ΦP(τ) is calculated as Φp(τ)=lim T→∞1/2T∫ T -T p′(t)p′(t+τ)dt... 4) It can be expressed as Here p'(t)=p(t)
*s(t), where s(t) corresponds to the sensitivity of the ear,
Theoretically, it is expressed by the characteristics of the middle ear and the outer ear, but in practice it can be expressed as the impulse response of a G filter, which is well known as an approximation of auditory characteristics. Naturally, the power Φp of p′(t)
By dividing (Equation 4) by (o), the normalized autocorrelation function φp(τ) is obtained as φp(τ) = Φp(τ)/Φp(o)...(
Equation 5) can be expressed as: FIGS. 2a and 2b illustrate the measured values of the normalized autocorrelation function corresponding to the above (Equation 5).
FIG. 2a shows the measured values corresponding to the music "Royal Pavane" by Gibbons, which is hereinafter referred to as phonology A. Figure 2b is the music “Sinho Nietsta” by Arnold, Op.
48: Movement, Allegro con brio (Synfornietta, Opus48: Movement,
This acoustics is referred to as music B hereafter.Next, the second objective parameter is the impulse generated by reflection at boundaries such as walls. This is related to the initial time delay between the direct sound and the first reflected sound, as well as the initial reflected sound, subsequent reverberation, and changes in the spectrum due to reflection. Objective parameters of 3 include the impulse responses h ol (t) and h or (t) to the left and right ears.These impulse responses play an important role in sound localization and are independent of each other. This is clear from the fact that in the case of a centrally localized signal, h ol (t) h or (t). Next, the impulse responses of these two h ol (t) and h or
Derive the mutual dependency between (t). First, we define the long-term interaural correlation function Φ lr (τ) between the binaural signals fl(t) and fr(t) as fl(t) and fr(t).
When expressed using Φ lr (τ)=lim T→∞1/2T∫ T -T f l (t) fr (t + τ) dt, |τ|≦|ms (Equation 6). By the way, a sense of spaciousness or a sense of no direction occurs when the interaural cross-correlation is a small value, and in the case of only signals from a fixed direction, the interaural cross-correlation is The cross-correlation function has a large peak. Here |τ|<1ms
This is because the time difference between the binaural signals fl(t) and fr(t) is, from the relationship between the binaural distance and the speed of sound,
This is because the time normally cannot exceed 1 ms. First, the interaural cross-correlation function Φ (o) lr (τ) with only direct sound is fl (t) = p (t) * h pl (t), f r (
t)=p
It can be obtained by substituting (t)*h pr (t) into (Equation 6). By the way, the normalized autocorrelation function φ (o) lr (τ) for only the direct sound is It can be expressed as Note that φ (o) l r(τ) is temporarily h pl
If (t) h pr (t), it becomes approximately 1. Here Φ (o) ll
(o) and Φ (o) rr (o) show the τ=0 autocorrelation function of the signal at the left and right ears, respectively. Next, if another reflected sound is added to the direct sound after the time when the autocorrelation function of the direct wave becomes small, then the normalized autocorrelation function φ (N) lr (τ) is W o (t )
is equal to the Dirac delta function δ(t), then It is expressed as Here, Φ (n) lr (τ) is the interaural interaction function of the n-th reflex as Φ (n) ll (o), Φ (n) rr (o)
represents the autocorrelation function at τ=0 of the nth reflected sound of the left and right ears, respectively. Note that when the sound source is in front of the user in a normal room, the maximum value of the interaural cross-correlation function is obtained extremely close to τ=0. By the way, the strength of the normalized interaural cross-correlation is
IACC as IACC = |φ lr (τ) | max, for | τ | 1ms
...(Formula 9) is defined as follows. Note that φl r (τ)=Φ lr (τ)/√ ll
(o) Φ rr (o), and Φ ll (o) and Φ rr (o) respectively indicate the autocorrelation functions at τ=0 of the signal tones for the left and right ears. The above describes the long-term autocorrelation function of the sound source signal, the multiple impulse responses caused by sound reflecting off walls, etc., and the interaural cross-correlation function that indicates the correlation between binaural signals. However,
Next, from these physical quantities, we can determine four elements that are important objective parameters that represent the characteristics of binaural sound signals: audible sound pressure, delay time of the first reflected sound,
We will discuss how to find the optimal values for the reverberation time of the subsequent reverberant sound and the value of the interaural correlation function. The results of the preference test conducted using the aforementioned music A and music B are shown in FIGS. 3a and 3b.
The horizontal axis is IACC, and the vertical axis is audible sound pressure (unit, dBA). The listening sound pressure can be considered to be the value of τ=0 of the autocorrelation function, and the value on the horizontal axis of the graph indicates the preference. From Figures 3a and b, it is clear that the optimal listening sound pressure [p] p does not depend much on IACC; it is 77 to 79 dBA for music A with a rather slow tempo, and 79 dBA for music B with a fast tempo. teeth
It can be seen that it is 79-80dBA. It can be seen that the optimal value of the listening sound pressure is around 79 dBA in both cases. Next, according to the results of evaluating a synthesized sound field consisting of a direct sound and a single reflected sound when reproduced by a speaker using music and speech, it is found that the normalized autocorrelation function φ p (τ) of the sound source signal is When the reflected sound level is varied by ±6 dB of the direct sound, the optimal delay time of the reflected sound is |Φ p (τ)| is 1/10 of the first reflected sound level G 1 It has become clear that the period corresponds to that time. There|
The horizontal axis is the time τd at which Φ p (τ) | corresponds to 1/10 of the level C 1 of the first reflected sound, and the vertical axis is the delay time [△t 1 ] p of the single reflected sound at which the preference is maximum. FIG. 4 is a diagram showing the axis. The range shown in the figure indicates the delay time when it is 0.1 lower than the maximum value of the preference, and the symbol ○ in the figure indicates the delay time when the preference is 0.1 lower than the maximum value.
Reflected sound level G 1 = 6 dB, ● is G 1 = 0 dB, □ is
Each figure shows the case when G 1 = -6 dB. In particular, if τ p is the time during which |Φ p (τ) becomes 0.1 times Φ p (o), when G 1 =0 dB, it can be expressed as τ d = τ p . Note that the τ p of music A and B mentioned above are 127 ms and 35 ms, respectively, as can be inferred from Figure 2 a and b.
As is clear from Figure 4, a straight line can be drawn to approximate each point in the figure, and when looking at both axes corresponding to the straight line, τd is the single reflection with the maximum preference. It can be seen that the sound delay time [△t 1 ] almost matches p . At the same time, it almost coincides with the time τ p when |Φ p (τ)| becomes 0.1 times Φ p (o).
That is, if this is expressed mathematically, the most preferable delay time [Δt 1 ] p of the first reflected sound can be expressed as [Δt 1 ]p=τ p (Formula 10). In addition, when τ>τ p , |φ p (τ)|≦KG c ... (Formula 11), and in this case K=const (=0.1), C=const (=
1.0),

【式】である。正規化自己相関 関数Φp(τ)を用いたとき上記(式11)は τ>τpにおいて、|Φp(τ)|≦0.1
……(式12) と表現できる。 さらに、音源信号の自己相関関数は最適残響時
間とも密接な関係があり、この関係を第5図に示
した。縦軸は後続残響音の最適残響時間〔Tsub〕
pであり、横軸はτpである。ここで言う残響時間
とは直接音が60dB減衰する時間ではなく、残響
部の信号が60dB減衰する時間として表現してい
る。図中、音楽Aと音楽Bは先に述べたものと同
様であるが、音楽Eはモーツアルト(Mozart)
による“交響曲ハ長調K、551ジユピター第4楽
章、モルトアレグロ”であり、スピーチSは国木
田独歩の“利根川の瀬の音”の「空はどこまでも
青く、明るく高いのに驚きます。」(τp=12ms)
である。図から明らかなように、第5図に示され
た関係は〔Tsub〕p≒(23±10)τpなる関係で
ほぼ近似することができる。 次に、第6図は直接音と単一反射音からなる合
成音場において測定した結果を図に表わしたもの
であつて、横軸に反射音の到来方向ξを縦軸に正
規化されたプリフアレンスとIACCの値を示して
ある。この第6図からプリフアレンスはIACCの
値が減少するに従つて増加することがわかる。即
ち、プリフアレンスのスコアとIACCの値との間
の相関関数は負(−0.76:1%有意水準)であつ
て、このことはIACCがτ=0で最大値をとる場
合に成立するものである。最も効果的にIACCを
小さくするには正面から±(55゜±20゜)の範囲に
初期反射音が到来するようにすればよいことも図
より読み出すことができる。 次に、上述の4つのパラメータによるプリフア
レンスの尺度について述べる。ところでこのプリ
フアレンスの尺度は比較試験によつて求められた
ものであり、各々のパラメータはプリフアレンス
レンスの尺度に対して独立に影響をあたえる。こ
の結果、重畳の理が適用できるので各々の客観的
なパラメーターを最適値によつて正規化すること
により、コンサートホールなどにおいて得られた
音場のプリフアレンスのデータを一般化すること
が可能となる。以下に比較試験によつて得られた
プリフアレンスの尺度について説明する。 まず、聴取音圧の関数としてのプリフアレンス
の尺度S1を第7図に示す。この図において、最適
聴取音圧における尺度は零に設定されている。S1
の値は、最も好ましい聴取音圧(0dB)を中心と
して左右ほぼ対称であるが、弱い聴取音圧の方に
ずれた方が、強い音圧の方にずれたときよりもや
やプリフアレンスの尺度が良い傾向がある。この
ことを数式によつて表現すれば、 S1(LL)≦S1(−LL) ……(式13) となる。ここで、LL=20log(p/〔p〕p)であ
り、pは聴取音圧、〔p〕pは最適な聴取音圧であ
る。なお、図中○は音楽A,Xは音楽Bの値を示
している。 次に、直接音と第1反射音の間の遅れ時間(第
1反射音遅れ時間)の関数としてのプリフアレン
スの尺度S2を第8図に示す。この図の横軸は最も
好ましい時間遅れ〔△t1pによつて正規化されて
いる。ところで、第2反射音の遅れ時間は、その
最適値〔△t2pが〔△t2p1.8τpであることが知
られているので、測定を行うにあたつては第2反
射音の遅れ時間としてはこの最適値を用いた。も
ちろん、(式10)の条件があることは言うまでも
ない。図において、○,a,A,△のそれぞれは
音楽Aによる個々の測定結果であり、X,b,
B,△のそれぞれは音楽Bによる個々の測定結果
であり、Cは音楽C、Dは音楽D、□は音楽E、
●はスピーチSによる測定結果である。音楽A,
B,Eについてはすでに先に述べたが、音楽Cは
ハイドンの交響曲第102番変ロ長調第2楽章アダ
ージヨ(τp=65ms)であり、音楽Dはワーグ
ナーのジークフリートイデイル(Siegfried
Idyll)332小節(τp=40ms)である。なお第6
図と同様最適な第一反射音の遅れ時間における尺
度は零に設定されている。 さらに、後続残響音の残響時間を関数としたプ
リフアレンスの尺度S3を第9図に示す。図におい
て実線は反射音の全ての音圧
[Formula]. When using the normalized autocorrelation function Φ p (τ), the above (Equation 11) becomes |Φ p (τ)|≦0.1 when τ>τ p
...(Formula 12) Furthermore, the autocorrelation function of the sound source signal has a close relationship with the optimal reverberation time, and this relationship is shown in FIG. The vertical axis is the optimal reverberation time of the subsequent reverberant sound [Tsub]
p , and the horizontal axis is τ p . The reverberation time referred to here is not expressed as the time for the direct sound to attenuate by 60 dB, but as the time for the signal in the reverberation section to attenuate by 60 dB. In the figure, music A and music B are the same as those mentioned above, but music E is Mozart.
``Symphony in C major K, 551 Jupiter 4th movement, molto allegro'', and Speech S is ``Symphony in C major K, 551 Jupiter 4th movement, molto allegro'', and Speech S is Doppo Kunikida's ``Sound of the Tone River'', ``I'm surprised at how blue the sky is, how bright and high it is.'' (τp =12ms)
It is. As is clear from the figure, the relationship shown in FIG. 5 can be approximately approximated by the relationship [Tsub]p≈(23±10)τp. Next, Figure 6 shows the results measured in a synthetic sound field consisting of a direct sound and a single reflected sound, where the horizontal axis represents the arrival direction ξ of the reflected sound, and the vertical axis represents the normalized Preference and IACC values are shown. It can be seen from FIG. 6 that the preference increases as the value of IACC decreases. That is, the correlation function between the preference score and the IACC value is negative (-0.76: 1% significance level), and this holds true when IACC takes the maximum value at τ = 0. . It can also be read from the figure that the most effective way to reduce IACC is to make the early reflected sound arrive within a range of ± (55° ± 20°) from the front. Next, a preference measure using the above four parameters will be described. Incidentally, this scale of preference was determined through a comparative test, and each parameter independently influences the scale of preference. As a result, since the principle of superposition can be applied, it becomes possible to generalize sound field preference data obtained in concert halls etc. by normalizing each objective parameter by its optimal value. . The preference scale obtained in the comparative test will be explained below. First, a measure of preference S 1 as a function of listening sound pressure is shown in FIG. In this figure, the scale at optimal hearing sound pressure is set to zero. S 1
The values of are almost symmetrical around the most preferable listening sound pressure (0 dB), but the preference scale is slightly lower when the listening sound pressure is lower than when the listening pressure is lower than when the sound pressure is higher. There is a good trend. If this is expressed mathematically, S 1 (LL)≦S 1 (−LL) (Formula 13). Here, LL=20log(p/[p] p ), p is the audible sound pressure, and [p] p is the optimal audible sound pressure. Note that in the figure, ○ indicates the value of music A, and X indicates the value of music B. Next, the measure S 2 of preference as a function of the delay time between the direct sound and the first reflected sound (first reflected sound delay time) is shown in FIG. The horizontal axis of this figure is normalized by the most preferred time delay [Δt 1 ] p . By the way, the delay time of the second reflected sound is known to have an optimal value [△t 2 ] p of [△t 2 ] p 1.8τ p . This optimal value was used as the delay time of reflected sound. Of course, it goes without saying that the condition of (Equation 10) exists. In the figure, ○, a, A, and △ are individual measurement results for music A, and X, b,
Each of B and △ is an individual measurement result using music B, C is music C, D is music D, □ is music E,
● is the measurement result by Speech S. Music A,
I have already mentioned B and E earlier, but music C is Haydn's Symphony No. 102 in B flat major, second movement Adagio (τp = 65ms), and music D is Wagner's Siegfried.
Idyll) 332 bars (τ p = 40 ms). Furthermore, the 6th
As in the figure, the optimal delay time scale for the first reflected sound is set to zero. Furthermore, the preference scale S 3 as a function of the reverberation time of the subsequent reverberant sound is shown in FIG. In the figure, the solid line is the total sound pressure of reflected sound.

【式】の場合の後続残響音の 残響時間を関数としたプリフアレンスの尺度S3
あり、破線はG=1.1の場合のそれである。図に
おいて、○,aは音楽Aについて、X,bは音楽
Bについて、□は音楽Eについて、●はスピーチ
Sについてのそれぞれ実験結果であり、いずれも
G=4.1の場合のものである。G=4.1はコンサー
トホールの後方にいる場合のような残響音が多い
場合に相当し、G=1.1は、ホールの前席におけ
るような直接音の多い場合に対応している。最も
好まれる残響時間におけるプリフアレンスの尺度
は零に設定されている。 次にIACCの関数としてのプリフアレンスの尺
度S4を第10図に示す。図において、○は音楽A
のとき、X,bは音楽Bのときの実験結果であ
る。IACCの性格上τ=0においてその最大値を
とれば、音像は正面方向にある。図から明らかな
ように、IACCが1に近づくとプリフアレンスの
尺度S4は急激に小さくなる。したがつて、できる
限りIACCは0.5より小さくなるようにする必要が
ある。 これまでにプリフアレンスの4つの尺度S1から
S4までについて述べたが、これらの4つの尺度に
ついてそれぞれ以下に示す近似式で近似すること
ができる。 まず、第7図に示した聴取音圧のプリフアレン
スの尺度S1は S1=−α1|X13/2 ……(式14) ただし、 X1=20log(p/〔p〕p
……(式15) と表わすことができる。なお、 α=0.07±0.03、X1>0 0.04±0.02、X1<0 ……(式16) 次に、第8図に示した第一反射音の遅れ時間の
プリフアレンスの尺度S2は S2=−α2|X23/2 ……(式17) ただし、 X2=log(△t1/〔△t1p) ……(式18) と表わすことができる。なお、 α2=1.42±0.6、X2>0 1.11±0.5、X2<0 ……(式19) である。 さらに、第9図に示した後続残響音の残響時間
のプリフアレンスの尺度S3は S3=−α3|X33/2 ……(式20) ただし、 X3=log(Tsub/〔Tsub〕)p)……(式21) と表わすことができる。なお、 α3=(0.74±0.25)G+ (0.45±0.15)、X3>0 −(0.42±0.14)G+ (2.36±0.79)、X3<0 ……(式22) であり もし、(式22)においてα3が負となるときには
α3=0……(式22′)とする。 残りの第10図に示したIACCのプリフアレン
スの尺度S4は S4=−α4X4 3/2 ……(式23) ただし、 X4=IACC ……(式24) と表わすことができる。なお、 α4=1.45±0.44 ……(式25) である。 重畳の理に基づいて、コンサートホールなどに
おけるプリフアレンスの全体の尺度Sとして S=4i=0 Si ……(式26) で表わすことができる。もちろん、Si(i=1、
2、3、4)は前述の4つのプリフアレンスの尺
度を示している。 このようにして得られたプリフアレンスの尺度
S及びSi(i=1、2、3、4)を用いることに
よつて音場の正確な評価を行うことができる。こ
の発明はコンサートホールや室内におけるプリフ
アレンスを測定するために音場における聴取音
圧、第1反射音の遅れ時間、後続残響音の残響時
間、および両耳間相互相関係数の4つのパラメー
タより、それぞれのプリフアレンスの尺度及び全
体のプリフアレンス尺度を求めることのできる測
定器を提供することが目的である。 (1) 本発明の音場評価計測器 以下、本発明の音場評価計測器の一実施例を
図を用いて説明する。 第11図及び第12図は本発明による音場評
価計測器の概略の構成を示したブロツク図であ
る。第12図は4つのパラメーターごとに示し
たブロツク図である。図において1は人頭また
はダミーヘツド、2l,2rは外耳道入口部に
挿入されたマイクロフオン、3,3l,3rは
前置増巾器、4は物理量解析器、5は比較器、
6は心理量変換器、7は総合評価器、8はプリ
フアレンスの出力端子、9はプリフアレンスを
記録するための記録器、10は音場評価計測
器、40は加算器、41は聴取音圧解析器、4
2は第1反射音遅れ時間解析器、43は後続残
響音残響時間解析器、44は両耳間相互相関関
数解析器、51は聴取音圧心理量変換器、52
は第1反射音遅れ時間心理量変換器、53は後
続残響音残響時間心理量変換器、54は両耳間
相互相関係数心理量変換器、y,y1,y2,y3
y4は比較データ入力端子である。その比較デー
タ入力端子y,y1,y2,y3,y4にはそれぞれの
最適値が入力しており、入力端子y1からは最適
聴取音圧(〔p〕p)に対応する信号が、入力端
子y2からは最適第1反射音遅れ時間(〔△t1p
に対応する信号が、入力端子y3からは最適後続
残響音残響時間(〔Tsub〕p)に対応する信号
が、入力端子y4からは最適両耳間相互相関係数
に対応する信号が入力する。 なお、第12図の聴取音圧解析器41・第1
反射音遅れ時間解析器42・後続残響音残響時
間解析器43・両耳間相互関関数解析器44は
第11図の物理量解析器4に対応する物理量解
析部に相当する。第12図の聴取音圧心理量変
換器51・第1反射音遅れ時間心理量変換器5
2・後続残響音残響時間心理量変換器53・両
耳間相互相関係数心理量変換器54は第11図
の比較器5と心理量変換器6とに対応する比較
及び心理量変換部に相当する。同様に総合評価
器7は総合評価部に相当する。 まず、最初に第11図を用いて発明の音場評
価計測器の一実施例の簡単な流れを説明する。
第1にコンサートホールやリスニングルームな
ど、音場評価を行い場所にダミーヘツドまたは
人頭1を設置する。その左右の外耳道入口に取
付けられたマイクロホン2l,2rによつて両
耳の音圧信号が検出され、前置増巾器3によつ
て増巾し加算されたのち、物理量解析器4によ
つて両耳の音圧信号から物理量即ち聴取音圧、
第1反射音遅れ時間、後続残響音残響時間およ
び両耳間相互相関係数IACCを求め、そのそれ
ぞれの値を比較器5において入力端子yより比
較データとして入力されるそれぞれの最適値と
比較し、次にすでに説明した(式14)(式17)
(式20)(式28)に基ずいた処理を行なう心理量
変換器6によつて心理量つまりプリフアレンス
に対応した量に変換し、その後総合評価器7に
おいて総合評価して音場におけるプリフアレン
スを出力端子8より出力するものであるこの計
測器によつてコンサートホールやリスニングル
ームの音場評価を行こなうことができる。 次に第11図に示した一実施例をさらに4つ
の個々のパラメータの流れまで示した第12図
により、詳細に説明する。 前述の通り、マイクロフオン2l,2rより
検出された音圧信号は、前置増巾器3l,3r
により増巾され、加算器40で加算された後、
物理量解析部に相当する聴取音圧解析器41、
第1反射音遅れ時間解析器42、後続残響音残
響時間解析器43のそれぞれに入力されそれぞ
れの解析器において聴取音圧p、第1反射音遅
れ時間△t1、後続残響音残響時間Tsubのそれ
ぞれが解析測定され出力される。また、この3
つの流れとは別に、加算器40で加算される前
の2つの前置増巾器3l,3rからの出力が、
物理量解析部に相当する両耳間相関関数解析器
44へ入力され、両耳間相互相関係数の最大値
IACCが解析され、出力される。次に、物理量
解析部に相当する4つの解析器41,42,4
3,44のそれぞれから出力される聴取音圧
p、第1反射音遅れ時間△t1、後続残響音残響
時間Tsub、最大両耳間相互相関係数IACCはそ
れぞれ、比較及び心理量変換部に相当する聴取
音圧心理量変換部51、第1反射音遅れ時間心
理量変換器52、後続残響音残響時間心理量変
換器53に入力され、他の比較データ入力端子
y1,y2,y3,y4のそれぞれから入力された最適
聴取音圧〔p〕p、最適第1反射音遅れ時間〔△
t1p、最適後続残響音残響時間〔Tsub〕p、最適
両耳間相互相関係数と比較される。最適両耳間
相互相関係数は0と考えてもよいが、実用上は
0.4〜0.5以下に設定しても問題ない場合が多い
ので、この実施例では0.4に設定している。 この比較及び心理量変換部における心理量変
換のための計算式は、聴取音圧心理量変換部5
1においては前述の(式14)〜(式16)、第1
反射音遅れ時間心理量変換器52においては前
述の(式17)〜(式19)、後続残響音残響時間
心理量変換器53においては前述の(式20)〜
(式22′)、両耳間相互相関係数心理量変換器5
4においては(式23)〜(式25)、の通りであ
つて、この比較及び心理変換部より得られる出
力であるところのプリフアレンスの尺度S1とS2
とS3とS4をマイコンプログラム等による計算処
理によつて求めることができる。しかし、あま
り精度を問題としないような場合は、下記に示
すような変換テーブルを用いて行つても良い。
It is the measure of preference S 3 as a function of the reverberation time of the subsequent reverberant sound in the case of [Equation], and the dashed line is that for the case of G=1.1. In the figure, ○ and a are experimental results for music A, X and b for music B, □ for music E, and ● for speech S, all of which are for the case of G=4.1. G=4.1 corresponds to a case where there is a lot of reverberant sound, such as when being at the rear of a concert hall, and G=1.1 corresponds to a case where there is a lot of direct sound, such as at the front seats of a hall. The preference scale at the most preferred reverberation time is set to zero. Next, the preference measure S 4 as a function of IACC is shown in FIG. In the diagram, ○ is music A
When , X, b are the experimental results for music B. Due to the nature of IACC, if the maximum value is taken at τ=0, the sound image is in the front direction. As is clear from the figure, as IACC approaches 1, the preference measure S 4 decreases rapidly. Therefore, it is necessary to make IACC smaller than 0.5 as much as possible. So far, four scales of preference S 1 to
Although up to S4 has been described, these four scales can be approximated by the approximate expressions shown below. First, the hearing sound pressure preference scale S 1 shown in Fig . 7 is S 1 = −α 1 | )
...(Formula 15) In addition, α=0.07±0.03, X 1 >0 0.04±0.02, X 1 <0 ... (Formula 16) Next, the measure of preference S 2 for the delay time of the first reflected sound shown in Fig. 8 is S 2 = −α 2 |X 2 | 3/2 (Formula 17 ) However, it can be expressed as : Note that α 2 =1.42±0.6, X 2 >0 1.11±0.5, X 2 <0 (Formula 19). Furthermore , the scale S 3 of the preference of the reverberation time of the subsequent reverberant sound shown in FIG. 9 is S 3 = −α 3 | Tsub]) p )... (Equation 21) In addition, α 3 = (0.74±0.25)G+ (0.45±0.15), X 3 >0 − (0.42±0.14 ) G+ (2.36±0.79), 22), when α 3 is negative, α 3 =0 (Formula 22'). The remaining IACC preference measure S 4 shown in FIG . 10 can be expressed as S 4 = −α 4 . Note that α 4 =1.45±0.44 (Formula 25). Based on the principle of superposition, the overall measure S of preference in a concert hall or the like can be expressed as S= 4i=0 Si (Formula 26). Of course, S i (i=1,
2, 3, and 4) indicate the four preference scales mentioned above. By using the preference scales S and S i (i=1, 2, 3, 4) obtained in this way, it is possible to accurately evaluate the sound field. This invention uses four parameters to measure preference in a concert hall or room: audible sound pressure in the sound field, delay time of the first reflected sound, reverberation time of the subsequent reverberant sound, and interaural cross-correlation coefficient. The object is to provide a measuring device with which individual preference measures and overall preference measures can be determined. (1) Sound field evaluation measuring instrument of the present invention Hereinafter, one embodiment of the sound field evaluating measuring instrument of the present invention will be described with reference to the drawings. FIGS. 11 and 12 are block diagrams showing the general structure of a sound field evaluation measuring instrument according to the present invention. FIG. 12 is a block diagram showing each of the four parameters. In the figure, 1 is a human head or dummy head, 2l, 2r are microphones inserted into the entrance of the external auditory canal, 3, 3l, 3r are preamplifiers, 4 is a physical quantity analyzer, 5 is a comparator,
6 is a psychological quantity converter, 7 is a comprehensive evaluator, 8 is a preference output terminal, 9 is a recorder for recording preferences, 10 is a sound field evaluation measuring device, 40 is an adder, 41 is audible sound pressure analysis vessel, 4
2 is a first reflected sound delay time analyzer, 43 is a subsequent reverberation sound reverberation time analyzer, 44 is an interaural cross-correlation function analyzer, 51 is a listening sound pressure psychological quantity converter, 52
is a first reflected sound delay time psychological quantity converter, 53 is a subsequent reverberation sound reverberation time psychological quantity converter, 54 is an interaural cross-correlation coefficient psychological quantity converter, y, y 1 , y 2 , y 3 ,
y4 is a comparison data input terminal. The respective optimum values are input to the comparison data input terminals y, y 1 , y 2 , y 3 , y 4 , and a signal corresponding to the optimum hearing sound pressure ([p] p ) is input from the input terminal y 1 . However, from input terminal y 2 , the optimal first reflected sound delay time ([△t 1 ] p )
A signal corresponding to the optimal subsequent reverberant sound reverberation time ([Tsub] p ) is input from the input terminal y 3 , and a signal corresponding to the optimal interaural cross-correlation coefficient is input from the input terminal y 4 . do. It should be noted that the hearing sound pressure analyzer 41 and the first
The reflected sound delay time analyzer 42, the subsequent reverberant sound reverberation time analyzer 43, and the interaural correlation function analyzer 44 correspond to a physical quantity analyzer corresponding to the physical quantity analyzer 4 in FIG. Hearing sound pressure psychological quantity converter 51 and first reflected sound delay time psychological quantity converter 5 in FIG.
2. The subsequent reverberation sound reverberation time psychological quantity converter 53 and the interaural cross-correlation coefficient psychological quantity converter 54 are the comparison and psychological quantity converting sections corresponding to the comparator 5 and the psychological quantity converter 6 in FIG. Equivalent to. Similarly, the comprehensive evaluator 7 corresponds to a comprehensive evaluation section. First, a simple flow of an embodiment of the sound field evaluation measuring device of the invention will be explained using FIG. 11.
First, a sound field is evaluated and a dummy head or head 1 is installed in a location such as a concert hall or listening room. The sound pressure signals of both ears are detected by the microphones 2l and 2r installed at the entrances of the left and right external auditory canals, amplified and added by the preamplifier 3, and then processed by the physical quantity analyzer 4. From the sound pressure signals of both ears, the physical quantity, that is, the audible sound pressure,
The first reflected sound delay time, the subsequent reverberant sound reverberation time, and the interaural cross-correlation coefficient IACC are determined, and the respective values are compared with the respective optimal values input as comparison data from the input terminal y in the comparator 5. , then the already explained (Equation 14) (Equation 17)
The psychological quantity converter 6, which performs processing based on (Equation 20) and (Equation 28), converts it into a psychological quantity, that is, a quantity corresponding to the preference, and then comprehensively evaluates it in the comprehensive evaluator 7 to determine the preference in the sound field. This measuring instrument, which is output from the output terminal 8, can be used to evaluate the sound field of a concert hall or listening room. Next, the embodiment shown in FIG. 11 will be explained in detail with reference to FIG. 12, which shows the flow of four individual parameters. As mentioned above, the sound pressure signals detected from the microphones 2l and 2r are transmitted to the preamplifiers 3l and 3r.
After being amplified by and added by an adder 40,
an audible sound pressure analyzer 41 corresponding to a physical quantity analysis section;
It is input to each of the first reflected sound delay time analyzer 42 and the subsequent reverberant sound reverberation time analyzer 43, and in each analyzer, the audible sound pressure p, the first reflected sound delay time Δt 1 , and the subsequent reverberant sound reverberation time Tsub are input. Each is analyzed, measured, and output. Also, these 3
Apart from the two flows, the outputs from the two preamplifiers 3l and 3r before being added by the adder 40 are
The maximum value of the binaural cross-correlation coefficient is input to the binaural correlation function analyzer 44, which corresponds to the physical quantity analysis unit.
IACC is parsed and output. Next, four analyzers 41, 42, 4 corresponding to the physical quantity analysis section are
The audible sound pressure p, first reflected sound delay time Δt 1 , subsequent reverberation sound reverberation time Tsub, and maximum interaural cross-correlation coefficient IACC output from each of It is input to the corresponding listening sound pressure psychological quantity converter 51, first reflected sound delay time psychological quantity converter 52, subsequent reverberation sound reverberation time psychological quantity converter 53, and other comparison data input terminals.
Optimal hearing sound pressure [p] p input from each of y 1 , y 2 , y 3 , and y 4 , optimal first reflected sound delay time [△
t 1 ] p , the optimal subsequent reverberation sound reverberation time [Tsub] p , and the optimal interaural cross-correlation coefficient are compared. The optimal binaural cross-correlation coefficient can be considered to be 0, but in practice
In many cases, there is no problem even if the value is set to 0.4 to 0.5 or less, so in this embodiment, it is set to 0.4. The calculation formula for this comparison and the psychological quantity conversion in the psychological quantity converting section is the hearing sound pressure psychological quantity converting section 5.
1, the above-mentioned (Equations 14) to (Equations 16), the first
The reflected sound delay time psychological quantity converter 52 uses the above-mentioned (Formula 17) to (Formula 19), and the subsequent reverberation sound reverberation time psychological quantity converter 53 uses the above-mentioned (Formula 20) to
(Formula 22'), interaural cross-correlation coefficient psychological quantity converter 5
In 4, (Formula 23) to (Formula 25) are as follows, and the preference scales S 1 and S 2 are the outputs obtained from this comparison and psychological conversion section.
, S 3 and S 4 can be obtained through calculation processing using a microcomputer program or the like. However, if accuracy is not a big issue, a conversion table as shown below may be used.

【表】【table】

【表】【table】

【表】 上記変換テーブル,はそれぞれ聴取音圧
心理量変換部51と第1反射音遅れ時間心理量
変換器52の変換テーブルの一例である。この
変換テーブル,には聴取音圧pと最適聴取
音圧〔p〕p、第1反射遅れ時間△t1と最適第1
反射遅れ時間との値に対応したそれぞれのプリ
フアレンスの尺度S1,S2が記憶されている。こ
こでは聴取音圧心理量変換部51と第1反射音
遅れ時間心理量変換器52の変換テーブルしか
示していないが、後続残響音残響時間心理量変
換器53と両耳間相互相関関係数心理量変換器
54に関しても同様な変換テーブルを作成して
行なうことができる。 このようにして比較及び心理量変換器の変換
器51〜54から出力されたプリフアレンス尺
度S1〜S4は、総合評価器7に入力させられ、総
合評価を行なつた結果として得られた全体のプ
リフアレンスの尺度Sが出力端子8に出力され
る。この出力は記録器9に記録される。ただ
し、この一実施例においては記録器9を含んで
いるが、かならずしも含んでいる必要はなくま
た、なくても良いし、表示器に変ても良いこと
はいうまでもない。 つぎに、第12図の音場評価計測器10に用
いられている物理量解析部の聴取音圧解析器4
1、第1反射音遅れ時間解析器42、後続残響
音残響時間解析器43および両耳間相互相関関
数解析器44のそれぞれの構成の一例について
第13図a〜eを用いて説明する。 第13図aは相関計を変形したM.R.シユレ
ーダー(M.R.Schroeder)による2乗積分形残
響計を用いた第1反射音遅れ時間解析器42あ
るいは後続残響音残響時間解析器43の基本構
成を示した図であつて、音圧に対応した出力も
含まれているという聴取音圧解析器41の機能
を有するものである。第13図bは第1反射音
遅れ時間解析器42あるいは後続残響音残響時
間解析器43の第1反射音または残響時間に対
応したコードを出力するために用いられるプラ
イオリテイエンコーダの一例を示した図であ
る。第13図cとdは聴取音圧解析器の他の構
成例を示した図である。第13図eは両耳間相
互相関係数を求める回路の構成例を示す図であ
る。 第13図aにおいて、400は2乗積分形残
響計を用いた第1反射音遅れ時間解析器あるい
は後続残響音残響時間解析器として用いられる
回路の入力端子、401,402,403は遅
延回路、404〜411は乗算器、412〜4
15は積分器、416〜418は加算器、41
9〜421は比較器、422は自己相関係数
φ0の出力端子、423〜425はそれぞれ第
1反射音の遅れ時間または後続残響音残響時間
の情報を出力するO1〜Oo-1の出力端子、42
6は減衰器、第13図bにおいて、430は出
力端子423〜425から出力されるO1
Oo-1が入力するプライオリテイーエンコーダ、
x1〜xo-1は第1反射音遅れ時間または後続残響
音残響時間に対応したコード出力、第13図
c,dにおいて、435は整流回路、436は
低減通過フイルタ、437は乗算器、第13図
eにおいて、441は両耳間相互相関係数を求
める回路の左側入力信号のLinの入る入力端
子、442は同じく右側入力信号Rinの入る入
力端子、443,444は加算器、445〜4
48は自乗平均回路RMS、449は演算回路、
450は相互相関係数に対応した出力をとり出
す出力端子である。 なお第13図aの回路においては、遅延回
路、加算器及び比較器をそれぞれ(n−1)
個、遅延回路からの出力信号が入力する乗算
器、積分器、積分器からの出力信号が入力する
と乗算器及び出力端子をそれぞれn個備えてい
るが、図面上では(n−1)個のものは代表し
て3個、n個のものは代表して4個だけ猫いて
あり、残りのものは省略してあり、図中に猫い
たものに連続符号を取つてある。したがつて2
番目の遅延回路402の次に(n−1)番目の
遅延回路403として3番目から(n−2)番
目の遅延回路には便宜上、符号を取つていな
い。 まず、第13図aについて説明する。入力端
子400に信号が与えられるとその入力信号
は、乗算器404〜407に与えられ、また遅
延回路401にも与えられる。遅延回路401
の出力は乗算器405へ与えられるとともに次
の遅延回路402にも与えられ順次遅れた信号
を(n−1)番目の遅延回路403の出力まで
を作り出す。それぞれの遅延回路401〜40
2の出力は乗算器405〜407によつて掛算
される。乗算器404は入力信号そのものの自
乗演算を行つている。次に、その乗算器404
〜407の出力信号はそれぞれ積分器412〜
415において積分されその出力信号はさらに
乗算器408〜411によつて自乗される。な
お、積分器412の出力信号は乗算器408へ
の入力とは別にそのまま端子422に出力され
る。この信号は自己相関係数φ(o)つまり入
力信号のパワーの自乗であり、音圧信号に対応
した情報を含んでいる。パワーの自乗の処理は
単に絶対値化としても等価な機能を持たせるこ
とができる。ところで乗算器408〜411の
出力信号は、つまり、自己相関係数の自乗φ2 p
2 1,……,φ2 o-1の信号は加算器416〜418
に加えられ、加算器416〜417の出力はそ
れぞれ次の加算器に加えられ、(n−1)番目
の加算器418まで行なわれる。ところで、乗
算器408からの出力信号と、加算器416〜
417からの出力信号とは、それぞれ、比較回
路419〜421において、最後尾の加算器4
18の出力信号を減衰器426で減衰させた後
の信号と比較され、出力端子423〜425に
それぞれの出力信号O1〜423〜Oo-1425
を出力する。なお減衰器426の減衰比は第1
反射音遅れ時間を求める場合には(0.1)2に、
後続残響音残響時間を求める場合には
(0.001)2に設定されたものが使用される。これ
はすでに最適遅れ時間検出器として知られたも
のと同様のものであつて、自己相関関数の値が
1/10あるいは1/1000になる遅れ時間を求めるも
のである。また遅延回路401から遅延回路4
03までの遅延時間の合計は第1反射音遅れ時
間及び後続残響音残響時間の長さ程度、つまり
それぞれ100〜200ms及び3〜5sec程度必要で
ある。当然ながら、この時間の経過に伴なつ
て、積分回路412〜415はリセツトされる
必要がある。ただし、リセツト回路は図示して
いない。 次に第13図bに示すプライオリテイーエン
コーダ430は、第13図aに示した比較器4
19〜421の出力信号O1423〜Oo-142
5をエンコードするもので、第13図aの減衰
器426の減衰比(0.1)2または(0.001)2に対
応した遅れ時間τ(=τp)(最初に減衰比以下に
小なくなつた遅れ時間)に対応したコードx1
xo-1を表示、又は出力する。 次の第13図cとdとは聴取音圧解析器とし
て用いる他の回路例を示す図であるが、これの
場合も第13図aの出力端子422から出力信
号φ0得るのと同様に、聴取音圧を知ることが
できる。この構成は、第13図cに示した通り
低域通過フイルタ436と整流回路435又は
自乗乗算器437を組み合わせた回路である。 次の第13図eの両耳間相互相関係数を求め
る回路519は正規化相関係数を求めるための
回路方式の1例を用いている。まず入力端子4
41から入力信号Linが、もう一方の入力端子
442から入力信号Rinがそれぞれ入力される
と、入力信号Lin,Rinはそれぞれ直接自乗平
均回路445,448に入りそれぞれの自乗平
均値γ,δとなり出力される。また入力信号
LinとRinとの間の和信号を加算器443より、
差信号を一方の信号を負にして他方の信号に加
える加算器444より、それぞれ取り出して自
乗平均回路446,447に入力して、同様に
和信号、差信号をそれぞれの自乗平均値α,β
を出力する。このようにして得られた4つの自
乗平均値の信号α,β,γ,δは演算回路44
9において演算処理つまり(α2−β2)/(4γδ)
を行ない、出力端子450に相互相関係数に対
応した出力φlr(o)つまり両耳間相互相関係数
(IACC)が得られる。ただし、第13図eの
回路を用いる場合には、入力信号LinとRinが
入力端子441,442に入力する前に0〜1
msの遅延回路を設け、その遅れ時間を変化さ
せ、その最大値を求めるように構成すると、前
述の(式9)に対応する回路構成となり、その
回路より(式9)に忠実な値が求められること
になる。そのための遅延回路は第13図e中に
は図示していない。なお、より精度の高い両耳
間相互相関係数(IACC)の演算を行うには前
述の(式7)〜(式9)までの式を用いて数値
演算を実行するとよい。 ここまでに説明してきた一実施例においては
聴取音圧解析器41、第1反射音遅れ時間解析
器42及び後続残響音残響時間解析器43への
入力は、ダミーヘツドまたは人頭1の左右の耳
に入つてくる2つの音響信号の加算信号を用い
てきたが、通常、その2つの信号間の遅れ時間
は1ms以下であり、これらの3つの解折器4
1〜43の解析結果に大きい差異を与えること
はないので、左右いずれか一方の音響信号を使
用しても十分な解析結果を得ることはできる。
これまでに説明してきた第1発明の一実施例で
あるところ音場評価計測器10の最適聴取音圧
〔p〕pは音源の種類によつても異なるが、それ
らはほぼ(79±5)dB程度であるので、ここ
では代表して79dBを用いた。また、最適第1
反射音遅れ時間〔△t1pと最適後続残響音残響
時間〔Tsub〕pとは実際の音源そのものから第
1反射音遅れ時間解析器42により最初に減衰
比以下に小さくなつた遅れ時間τpを求めてこれ
を最適第1反射音遅れ時間〔△t1pとすると良
く、その(23±10)倍をして最適後続残響音残
響時間〔Tsub〕pとして用いればよい。この一
実施例では〔Tsub〕p=23τpを用いた。 以上のようにこの本発明の音場評価計測器に
よれば1個または2個のマイクロホンと、この
マイクロホンから得られた音響信号を増巾する
ための増巾器と、この増巾信号から聴取音圧p
を計測するための聴取音圧計測手段と、同様に
増巾信号から第1反射音遅れ時間△t1を計測す
るための第1反射音遅れ時間計測手段と、また
同様に増巾信号から後続残響音残響時間Tsub
を計測するための後続残響音残響時間計測手段
と、さらに同じく増巾信号から両耳間相互相関
係数IACCを計測するための両耳間相互相関係
数計測手段と、上記聴取音圧pと最適聴取音圧
〔p〕pとの比の10を底とする対数の20倍の値の
絶対値の3/2乗の値に、上記聴取音圧が最適聴
取音圧より大きいときには−(0.07±0.03)を、
小さいときには−(0.04±0.02)を乗じてなる
値をプリフアレンス(聴感上の心理良さ)の尺
度S1として出力する、つまり、前述の(式14)
〜(式16)に対応する出力を出す第1の比較変
換手段と、上記第1反射音遅れ時間△t1と最適
第1反射音遅れ時間〔△t1pとの比の10を底と
する対数の絶対値の3/2乗の値に、上記第1反
射音遅れ時間が最適第1反射音遅れ時間より大
きいときには−(1.42±0.6)を、小さいときに
は−(1.11±0.5)を乗じてなる値をプリフアレ
ンスの尺度S2として出力する、つまり、前述の
(式17)〜(式19)に対応する出力を出す第2
の比較変換手段と、上記後続残響音残響時間
Tsubと最適後続残響音残響時間〔Tsub〕pとの
比の10を底とする対数の絶対値の3/2乗の値に、
上記後続残響音残響時間が最適後続残響音残響
時間より大きいときには、反射音全体の音圧の
−(0.74±0.25)倍の値と−(0.45±0.15)との
加算値を、小さいときには反射音全体の音圧の
−(0.42±0.14)倍の値と−(2.36±0.79)との
加算値を乗じたのちの値が負なるときにはその
ままの値を、正なるときには零をプリフアレン
スの尺度S3として出力する、つまり、前述の
(式20)〜(式22′)に対応する出力を出す第3
の比較変換手段と、上記両耳間相互相関係数の
3/2乗の値に−(1.45±0.44)を乗じた値をプリ
フアレンスの尺度S4として出力する、つまり前
述の(式23)〜(式25)に対応する出力を出す
第4の比較変換手段と、第1から第4までの比
較変換手段より出力されたプリフアレンスの尺
度S1,S2,S3,S4を総合評価してプリフアレン
スの全体の尺度Sを出力する総合評価手段とを
備えているので正確で忠実な音場評価を行なう
ことができるという効果を奏するものである。 (2) 本発明の応用例 次に音場における聴取音圧、第1反射音の遅
れ時間、後続残響音の残響時間、および両耳間
相互相関係数の4つの物理的パラメータに基づ
いて、音場を補正し、より好ましい音場を作る
ことのできる本発明を応用した第1、第2の音
響装置の各一実施例を図を用いて説明する。 (2‐1) 第1の音響装置 第14図aは音場を補正し、より好ましい
音場を作ることのできる第1の音響装置の一
実施例の構成ブロツク図であつて、第14図
bは第14図aの音響装置に含まれる音場拡
大装置の1例、第14図cは第14図aの音
響装置に含まれる残響装置の他の例を示した
ものである。 まず第14図aにおいて、501は入力端
子INR、502は入力端子INL、503a〜
503dは加算器、504a,504bは残
響装置、505a,505b,505r,5
05lは減衰器、506a,506bは遅延
回路、507は1種の分周回路であるスケー
ラ、508a,508bはアツプダウンカウ
ンタ、509はデイジタル−アナログ変換
器、510a,510cは第1反射音遅れ時
間解析器、510bは後続残響音残響時間解
析器、511a〜551cはエンコーダ、5
12a,512b,513a,513bは比
較器、514a,514b,514cは平滑
回路、515a,515bは目標値を設定す
るための可変抵抗器、516は音場拡大装
置、517r,517lはスピーカを駆動す
るためのパワーアンプ、518r,518l
はスピーカ、519は両耳間相互相関係数
(IACC)計算器、520は乗算器、531,
532は残響装置504a,504bから音
場拡大装置516への入力端子inr,inl、5
21はCLOCK信号入力端子、533,53
4は音場拡大装置516の出力端子outr,
outl、535は比較器513bから音場拡大
装置516への入力端子attmである。第1
4図bにおいて、530a〜530dは加算
器、536は減衰器、537は低域通過フイ
ルター、538は高域通過フイルター、53
9は移相回路又は遅延回路である。第14図
cにおいて、504cはM.R.シユレーダー
(Schroeder)とB.S.アタル(Atal)による
残響器、540a〜540iは加算器、54
1a,451bは位相反転回路、542a,
542bはボルテージフオロー、543a〜
543fは遅延回路、544a〜544fは
減衰器、545は減衰器504cの入力端
子、546は減衰器504cの出力端子であ
る。なお、図中、他の図面中の符号と同一符
号は同一または相当部分を示している。 次に第14図a,b,cを用いて、第1の
音響装置の一実施例の動作を説明する。ま
ず、入力端子INR501とINL502にレコ
ード、テープ、マイク、ラジオ等からの音響
信号が与えられる。この音響信号は加算器5
03cによつて加算され、第1反射音遅れ時
間解析器510cに入力される。この第1反
射音遅れ時間解析器510cは先に述べた第
12図の42、つまり第13図aに示された
回路と同じ回路によつて構成されており、こ
の場合、第1反射音遅れ時間を測定できるよ
うに減衰比は1/10に設定されている。この
第1反射音遅れ時間解析器511の出力はエ
ンコーダ511に入力され、エンコーダ51
1からは、音源にとつて最適な第1反射音の
遅れ時間〔△t1pをτpとして求めた出力
〔(△t1p=τp)が出る。このエンコーダ51
1cからの出力は比較器512aに入力され
るとともに、乗算器520にも入力され、23
倍されて最適後続残響音残響時間〔Tsub〕p
(=23゜τp)として比較器512bに与えられ
る。 一方、リスニングルームなどの音場に置か
れたダミーヘツドまたは人頭1の左右の耳の
所に設けられたマイクロフオン2r,2lに
入つてきた信号を前置増巾器3l,3rで増
巾し、その出力を加算器503dによつて加
算したのち、先と同様に第12図の42と同
じ第1反射音遅れ時間解析器510aに入力
し、かつ、後続残響音残響時間解析器510
bにも入力される。この後続残響音残響時間
解析器510bの構成は第12図の43にお
いて述べたものと同じ回路構成のものであ
る。次に第1反射音遅れ時間解析器510a
の出力はエンコーダ511aに入力され、音
場の影響を付与された後の第1反射音の遅れ
時間△t1に対応した信号が得られる。この
後、比較器512aに入力される。また、後
続残響音残響時間解析器510bの出力はエ
ンコーダ511bに入力され、音場の影響を
付与された後の後続残響音残響時間Tsubに
対応した信号が得られ、この信号は比較器5
12bに入力される。ところで、比較器51
2aにおいて、エンコーダ511cからの最
適第1反射音遅れ時間〔△t1pと、エンコー
ダー511aからの音場の影響を付与された
後の第1反射音の遅れ時間△t1とが比較さ
れ、もし最適第1反射音遅れ時間〔△t1p
比して第1反射音遅れ時間△t1の方が大きい
場合には、第1反射音遅れ時間△t1を小さく
するために、比較器512aからカウンタ5
08aにカウントダウンの信号CDを与え、
カウンタ508aの内容を小さくし、そのた
めに、カウンタ508aの内容に基づいて分
周を行うスケーラ507においては分周する
値が小さくなり、入力端子521からの
CLOCK信号を分周して得られるDELAY
CLOCK信号の周波数が高くなり、その
DELAY CLOCK信号により残響回路504
a,504bの中のBBD(バケツト・ブリゲ
ート・デイバイス)などによつて構成される
遅延回路506a,506bの遅延時間が小
さくなるというように動作するものである。
しかし、もし最適第1反射音遅れ時間〔△
t1pに比して第1反射音遅れ時間△t1の方が
小さい場合には比較器512aからカウンタ
508aにカウントアツプの信号CUを与え、
カウンタ508aの内容が大きくなり、スケ
ーラ507において分割する値が大きくな
り、そのためDELAY CLOCK信号の周波数
が低くなり、遅延回路506a,506bの
遅延時間が大きくなるというように上述の場
合と反対の動作をする。 また、比較器512bにおいては、乗算器
520の最適後続残響音残響時間〔Tsub〕p
と、エンコーダ511bからの音場の影響を
付与された後の後続残響音残響時間Tsubと
が比較され、もし、最適後続残響音残響時間
〔Tsub〕pに対して、後続残響音残響時間
Tsubが大きい場合には、カウンタ508b
にカウントアツプの信号CUが送出され、カ
ウンタ508bの内容がデイジタル−アナロ
グ変換器509によつてアナログ値に変換さ
れ、平滑回路514cを経て、残響回路50
4a,504bのそれぞれの減衰器505
a,505bに与えられ、その信号によつ
て、減衰器505a,505bにおける減衰
値が大きくなり、そのために後続残響音残響
時間Tsubを小さくするというように動作す
る。しかし、もし、最適後続残響音残響時間
〔Tsub〕pに対して後続残響音残響時間Tsub
が小さい場合には、カウンタ508bにカウ
ントダウンの信号CDが送出され、カウンタ
508bの内容がデイジタル−アナログ変換
器509によつてアナログ値に変換され、平
滑回路514cを経て減衰器505a,50
5bに与えられ、その信号によつて減衰器5
05a,505bにおける減衰値が小さくな
り、そのために残響時間Tsubを大きくする
というように上述の場合と反対の動作を行な
う。 ところで前置増巾器3l,3rからの加算
される前のそれぞれ増巾信号は、両耳間相互
相関係数計算器519に入力され、両耳間相
互相関係数IACCが求められてアナログ化さ
れ、次の平滑回路514bに入力される。次
に比較器513bにおいて、平滑回路514
bからの信号と、両耳間相互相関係数IACC
の目標値を設定する可変抵抗器515bから
の目標電圧信とを比較する。もし、平滑回路
514bからの信号、つまり両耳相互相関係
数IACCに対応する信号の方が、可変抵抗器
515bからの目標信号より大きい場合は、
比較器513bからの出力信号は小さくなつ
て、音場拡大装置516の入力端子att in5
35に入力され、音場拡大装置516はより
両耳相互相関係数IACCを小さくするという
ように動作する。しかし、もし、両耳相互相
関係数IACCに対応する信号の方が、可変抵
抗器515bからの目標信号より小さい場合
は、比較器513bからの出力信号は大きく
なつて、音場拡大装置の入力端子att in53
5に入力され、音場拡大装置516はより両
耳相互相関係数IACCを大きくするというよ
うに上述の場合と反対の動作をする。 ところで加算器503dで加算された加算
信号を入力された第1反射音遅れ時間解析器
510aは、聴取音圧に対応した音圧信号
φ0を出力しており、この音圧信号φ0を平滑
回路514aを通した後、比較器513aに
入力させる。この比較器513aにおいて、
音量の目標値を設定する可変抵抗器515a
からの目標となる電圧信号と、平滑回路51
4aからの信号とを比較して、もし、聴取音
圧に対応した平滑回路514aからの出力信
号の方が、可変抵抗515aからの音量の目
標値よりも大きい場合は、比較器513aか
らの出力が大きくなり、減衰器505r,5
05lの減衰が大きくなつて聴取音圧を下げ
るというように動作する。しかも、もし、聴
取音圧に対応する信号の方が、音量の目標値
よりも小さい場合は、比較器513aからの
出力が小さくなり、減衰器505r,505
lの減衰が小さくなつて聴取音圧を上げると
いうように上述の場合と反対の動作をする。 なお、第14図a内に出てきた平滑回路5
14a,514b,514cはいずれも、急
激な信号の変化による雑音の発生を防止する
ためと、残響回路504a,504b内の遅
延回路506a,506bの遅延時間も急激
な変化を起こさないようにするために用いら
れている。 次に、音場拡大装置516の働きを、第1
4図bに示す音場拡大装置を用いた場合で説
明する。 まず残響装置504a,504bから入力
端子inr531及びinl532に入つた2つの
信号の差成分が、加算器530aによつて得
られ、減衰器536を経て低域通過フイルタ
537を通つた後、クロストークする位相が
逆相になるような位相で、加算器530b,
530cにより入力端子inr531から出力
端子outr533までと、入力端子inl532
から出力端子outl534までのそれぞれの主
経路に加算される。それとは別に、低域通過
フイルタ537の出力はさらに、高域通過フ
イルタ538を経て、移相回路又は遅延回路
539において遅延または移送回転が与えら
れた後、同じく、クロストークする位相が逆
相となるような位相で、上述のそれぞれの主
経路に加算器530d,530eによつて加
算され、その加算後の信号は出力端子outr5
33とoutl534にそれぞれ出力される。こ
の音場拡大装置516の減衰器536を変化
することによつて両耳間相互相関係数IACC
を変化させることができる。またこの減衰器
536は、入力端子att in535からの信号
の値つまり、比較器513bからの出力信号
が大きくなると減衰が大きくなるように設定
されている。 さらに第14図cに示してあるのはM.R.
シユレーダー(Schroeder)とB.S.アタル
(Atal)による残響器504cで、第14図
aに示してある残響装置504a,504b
の代りにこの残響器504cを用いれば、よ
り自然な残響感が得られる。 なお、上記一実施例において、マイクロホ
ン2r,2lをダミーヘツドまたは人頭1の
外耳道入口の音圧を得るように述べて来た
が、両耳間相互相関係数IACCの値を得る必
要がないときは1本のマイクロホンでもよい
し、両耳間相互関係数IACCを求める場合で
もダミーヘツドなどに取付けられたものでな
い2本のマイクロホンによつて行つてよい場
合もある。音響信号については2チヤンネル
のステレオ信号のように示したが、必ずし
も、2チヤンネル信号の必要のない場合もあ
ることは言うまでもない。 以上のように、第1の音響装置において
は、入力端子INR,INLとこの入力端子INR
INLから入つてきた音響信号より最適第1反
射音遅れ時間〔△t1pを計測する最適第1反
射音遅れ時間計測手段と、その最適第1反射
音遅れ時間〔△t1pの(23±10)倍に対応す
る最適後続残響音残響時間〔Tsub〕pを出力
する最適後続残響音残響時間出力手段と、音
場におかれたマイクロフオンからの音場信号
より音場における聴取音圧φ0を計測する聴
取音圧計測手段と、その音場信号より第1反
射音遅れ時間△t1を計測する第1反射音遅れ
時間計測手段と、その音場信号より後続残響
音残響時間Tsubを計測する後続残響音残響
時間計測手段と、その音場信号より両耳間相
互相関係数IACCを計測する両耳間相互相関
係数計測手段と、上記最適第1反射音遅れ時
間〔△t1pと上記第1反射音遅れ時間△t1
を比較し、その差に応じた信号を出力する第
1の比較手段と、上記最適後続残響音残響時
間〔Tsub〕pと上記後続残響音残響時間Tsub
を比較し、その差に応じた信号を出力する第
2の比較手段と、上記第1の比較手段の出力
信号によつて遅延時間が変化し、上記第2の
比較手段の出力信号によつて残響時間が変化
し、かつ上記入力端子から入力する音響信号
に残響音を付加する残響手段と、あらかじめ
目標とする両耳間相互相関係数の値を設定で
きる両耳間相互相関係数設定手段と、上記の
両耳間相互相関係数IACCと上記の両耳間相
互相関係数設定手段の設定値とを比較し、そ
の差に応じた信号を出力する第3の比較手段
と、上記残響手段の出力信号を受け、上記の
第3の比較手段の出力信号に応じて、音場の
両耳間相互相関係数IACCを変化することの
できる出力を有する音場拡大手段と、あらか
じめ目標とする聴取音圧の値を設定できる聴
取音圧設定手段と、上記の音場の聴取音圧
φ0と上記聴取音圧設定手段の設定値を比較
し、その差に応じた信号を出力する第4の比
較手段と、上記音場拡大手段の出力を入力と
して受け、上記第4の比較手段の出力信号に
応じて、減衰率を変化できる減衰手段と、上
記減衰手段からの信号を増巾し、空間に音響
信号を放射する電気音響変換手段とを備えて
いるので、音場を補正し、より好ましい音場
を作ることができるというだけでなく、この
第1の音響装置内で音源について最適な第1
反射音遅れ時間〔△t1pと最適残響時間
〔Tsub〕pが得られ、使用者の好みに応じて聴
取音圧や両耳間相互相関係数IACCの値を変
えることができるという効果を奏するもので
ある。 (2‐2) 第2の音響装置 次に、比較的反射の少ない和室や反射が多
くても、残響時間が極度に短い車室内のよう
な音場における音場装置においては、音場に
おける聴取音圧、第1反射音遅れ時間、後続
残響音残響時間の3つの物理的パラメータの
みに基づいて音場を補正しだけで、十分なよ
り好ましい音場を作ることができる。 なお、前述の第1の音響装置をここで対象
としている音場での音響装置として用いるこ
とができるのは言うまでもないが、使用音場
がこのように小さな部屋、比較的反射の少な
い和室や車室内に限定された場合、両耳間相
互相関係数IACCのパラメータを使用しなく
ても第2の音響装置は前述の第1の音響装置
と同一の効果を得られるばかりでなく、音響
装置がコンパクトになり、しかも安価で、そ
の上、回路構成が前述の第1の音響装置より
簡単になるため、作業性が良く、しかも耐久
性が良くなるなどの点で第1の音響装置より
以上の効果を奏することができる。 そこで、音場における聴取音圧、第1反射
音遅れ時間及び後続残響音残響時間の3つの
物理的パラメータを用いて音場を補正し、よ
り好ましい音場を作ることのできるこの第2
の音響装置の一実施例を図を用いて説明す
る。 第15図が第2の音響装置の一実施例のブ
ロツク図である。 図において451〜453は絶対値化回
路、511dはエンコーダ、551は入力端
子INR、552は入力端子INL、553は積
分リセツト信号入力端子である。554は第
1反射音遅れ時間解析器で前述の第13図a
の回路図及び第14図aの第1反射音遅れ時
間解析器510a,510c、後続残響音残
響時間解析器510bのものと内部構成は多
少変えているが、機能的にはほとんど同じ
で、通常の自己相関器の構成に減衰器426
と比較器419を付加したことを特徴として
いる。このような構成の場合も、前述の第1
反射音遅れ時間解析器510a,510cあ
るいは後続残響音残響時間解析器510bと
同様に動作する。なお、図中第13図及び第
14図a中の符号と同一符号は、同一または
相当部分を示している。 まず入力端子INR551及びINL552か
ら音響信号が入力されと、その2つの信号は
加算器503cによつて加算せられ、第1反
射音遅れ時間解析器544に入力される。こ
の第1反射音遅れ時間解析器554におい
て、第14図aの第1反射音遅れ時間解析器
510a,510cと同様な動作を行ない、
そこからの出力は前述の第14図aの場合と
同様にエンコーダ511dに入力され、音響
信号にとつて最適の第1反射音遅れ時間〔△
t1pが出力される。この最適第1反射音遅れ
時間〔△t1pに対応した遅延時間が得られる
ようにスケーラ507にエンコーダ511d
からの信号が与えられる。スケーラ507は
入力端子521からのCLOCK信号を分周
し、DELAY CLOCK信号を作り、残響回路
504a,504bの遅延回路506a,5
06bにDELAY CLOCK信号が送られる。
音響信号は入力端子INR551とINL552
からそれぞれ残響回路504aと504bに
入力され残響音が付与されて減衰器505
r,505lに送られる。一方、第1反射音
遅れ時間解析器554の音圧に対応した信号
φ0が、平滑回路514aを経て、比較器5
13aに与えられ、音圧の目標値を設定する
可変抵抗器515aの電圧と比較されて、も
し、音圧が、その目標値より小さい場合に
は、比較器513aの出力は小さくなり、減
衰器505r,505lはその比較器513
aより送られてきた信号によつて減衰量を小
さくするというように動作する。しかし、も
し音圧が、その目標値より大きい場合には、
比較器513aの出力は大きくなり、減衰器
505r,505lはその比較器513aよ
り送られてきた信号によつて減衰量を大きく
するというよに上述の場合と反対の動作をす
る。そのようにして得られた減衰器505
r,505lの出力はパワーアンプ517
r,517lを経て、スピーカ518r,5
18lによつて再生される。 このとき、残響時間の最適値である最適後
続残響音残響時間〔Tsub〕pは、先に述べた
ように最適第1反射音遅れ時間〔△t1p
(23±10)倍であることが望ましい。残響回
路504a,504bのような形で残響器が
構成されているときには、後続残響音残響時
間Tsub=−ε・△t1/log(g)となること
が知られている。なおgは減衰率である。こ
こでTsub=(23±10)△t1なる関係を用いれ
ば、減衰率g=0.588〜0.811の範囲にしなけ
ればならないことがわかる。なお、後続残響
音残響時間Tsub=23△t1のときは減衰率g
=0.741である。即ち、このように減衰率g
を設定すれば、最適第1反射音残響時間〔△
t1pと最適後続残響音残響時間〔Tsub〕pが得
られる。もちろん、このような処理を処して
も、音場における影響が付加されるが、和室
など比較的反射の少ない空間や、自動車車室
など反射が多くても、残響時間が極度に短い
空間の場合には大変効果がある。 ところでこの第2の音響装置においては、
小さな部屋や比較的反射の少ない和室や車室
内のような音場において、第1の音響装置と
同一又はそれ以上の効果を奏し、より簡単な
音響装置を提供しようとするものであるが、
この第2の音響装置は、小さな部屋や車室内
以外の音場で使用した場合、第1の音響装置
よりは多少性能はおちるが、それでも、聴取
音圧、第1反射音遅れ時間及び後続残響音残
響時間の3つの物理的パラメータにより音場
を補正しているので、従来の音響装置よりは
るかに好ましい音場を作ることがで、その
上、第2の音響装置の方が第1の音響装置よ
り、コンパクトで安価で回路構成が簡単で、
作業性が良く耐久性も良いという利点があ
る。 以上のように、この第2の音響装置におい
ては、入力端子INR、INLと、この入力端子
INR,INLから入つてきた音響信号から最適
第1反射音遅れ時間を計測する最適第1反射
音遅れ時間計測手段と、上記音響信号から聴
取音圧に対応する信号を計測する聴取音圧対
応信号計測手段と、上記最適第1反射音遅れ
時間に対応する信号を制御信号として受け、
上記音響信号を入力し、上記音響信号に上記
最適第1反射音遅れ時間に対応した遅延時間
を付加するとともに、その後続残響音残響時
間が上記最適第1反射音遅れ時間の(23±
10)倍に設定された残響手段と、あらかじめ
目標とする聴取音圧の値を設定できる聴取音
圧設定手段と、上記聴取音圧対応信号計測手
段の出力である聴取音圧と上記聴取音圧設定
手段の設定値とを比較し、その差に対応する
信号を出力する比較手段と、上記残響手段の
出力信号が入力され、上記比較手段の出力信
号に応じて、上記残響手段から入力された音
響信号に与える減衰率を変化させることので
きる減衰手段と、上記減衰手段からの信号を
増巾し、空間に音響信号を放射する電気音響
変換手段とを備えているので、音場を補正
し、より好ましい音場を作るだけでなく、こ
の第2の音響装置内で、最適な第1反射音遅
れ時間と最適後続残響音残響時間が得られ、
使用者の好みに応じて聴取音圧を変えること
ができるという効果を奏するものである。 以上のように、この発明によれば、音楽などを
聞く音場の音響特性を測定して物理量および心理
量の音場における聴取音圧、第1反射音遅れ時
間、後続残響音残響時間、両耳間相互相関係数の
4つパラメータで評価するための計測器を提供す
ることができとという効果を奏する。
[Table] The above conversion tables are examples of conversion tables for the audible sound pressure psychological quantity converter 51 and the first reflected sound delay time psychological quantity converter 52, respectively. This conversion table contains the hearing sound pressure p, the optimal hearing sound pressure [p] p , the first reflection delay time △t 1 , and the optimal first
Preference scales S 1 and S 2 corresponding to the reflection delay time are stored. Although only the conversion tables for the audible sound pressure psychological quantity converter 51 and the first reflected sound delay time psychological quantity converter 52 are shown here, the following reverberation sound reverberation time psychological quantity converter 53 and the interaural cross-correlation coefficient psychological quantity converter 52 are shown. A similar conversion table can be created for the quantity converter 54 as well. The preference scales S 1 to S 4 outputted from the converters 51 to 54 of the comparison and psychological quantity converter in this way are input to the comprehensive evaluator 7, and the overall result obtained as a result of the comprehensive evaluation is A preference measure S is output to the output terminal 8. This output is recorded on the recorder 9. However, although the recording device 9 is included in this embodiment, it is not necessarily necessary to include it, and it goes without saying that it may be omitted or may be replaced with a display device. Next, the audible sound pressure analyzer 4 of the physical quantity analysis section used in the sound field evaluation measuring instrument 10 in FIG.
1. An example of the configuration of each of the first reflected sound delay time analyzer 42, subsequent reverberant sound reverberation time analyzer 43, and binaural cross-correlation function analyzer 44 will be explained using FIGS. 13a to 13e. FIG. 13a shows the basic configuration of the first reflected sound delay time analyzer 42 or the subsequent reverberant sound reverberation time analyzer 43 using a square integral type reverberation meter by MR Schroeder, which is a modified correlator. It has the function of the listening sound pressure analyzer 41 that also includes an output corresponding to sound pressure. FIG. 13b shows an example of a priority encoder used to output a code corresponding to the first reflected sound or reverberation time of the first reflected sound delay time analyzer 42 or the subsequent reverberant sound reverberation time analyzer 43. It is a diagram. FIGS. 13c and 13d are diagrams showing other configuration examples of the audible sound pressure analyzer. FIG. 13e is a diagram showing an example of the configuration of a circuit for calculating the interaural cross-correlation coefficient. In FIG. 13a, 400 is an input terminal of a circuit used as a first reflected sound delay time analyzer or a subsequent reverberant sound reverberation time analyzer using a square integral reverberation meter, 401, 402, and 403 are delay circuits; 404-411 are multipliers, 412-4
15 is an integrator, 416 to 418 are adders, 41
9 to 421 are comparators, 422 is an output terminal for the autocorrelation coefficient φ 0 , and 423 to 425 are O 1 to O o-1 output terminals that output information on the delay time of the first reflected sound or the reverberation time of the subsequent reverberant sound, respectively. Output terminal, 42
6 is an attenuator, and in FIG. 13b, 430 is O 1 to output from output terminals 423 to 425.
O o-1 input priority encoder,
x 1 to x o-1 are code outputs corresponding to the first reflected sound delay time or the subsequent reverberant sound reverberation time. In FIGS. 13c and d, 435 is a rectifier circuit, 436 is a reduction pass filter, 437 is a multiplier, In FIG. 13e, 441 is an input terminal into which the left input signal Lin of the circuit for calculating the interaural cross-correlation coefficient is input, 442 is an input terminal into which the right input signal Rin is also input, 443 and 444 are adders, and 445 to 4
48 is a root mean square circuit RMS, 449 is an arithmetic circuit,
450 is an output terminal that takes out an output corresponding to the cross-correlation coefficient. In the circuit of FIG. 13a, each of the delay circuit, adder, and comparator is (n-1).
The output signal from the delay circuit is input to the multiplier, the integrator is input, and the output signal from the integrator is input to the multiplier and output terminal. Three representative objects are shown, and only four representative objects are shown.The rest are omitted, and consecutive symbols are given to the objects shown in the figure. Therefore 2
The (n-1)th delay circuit 403 next to the (n-1)th delay circuit 402, and the third to (n-2)th delay circuits are not labeled for convenience. First, FIG. 13a will be explained. When a signal is applied to input terminal 400, the input signal is applied to multipliers 404 to 407 and also to delay circuit 401. Delay circuit 401
The output is given to the multiplier 405 and also given to the next delay circuit 402 to produce successively delayed signals up to the output of the (n-1)th delay circuit 403. Each delay circuit 401-40
The outputs of 2 are multiplied by multipliers 405-407. Multiplier 404 performs a square operation on the input signal itself. Next, the multiplier 404
~407 output signals are respectively integrators 412~
The output signal is integrated at 415 and further squared by multipliers 408-411. Note that the output signal of the integrator 412 is output as is to the terminal 422, apart from being input to the multiplier 408. This signal is an autocorrelation coefficient φ(o), that is, the square of the power of the input signal, and contains information corresponding to the sound pressure signal. The processing of squaring the power can also have an equivalent function simply by converting it into an absolute value. By the way, the output signals of the multipliers 408 to 411 are the square of the autocorrelation coefficient φ 2 p ,
2 1 , ..., φ 2 o-1 signals are sent to adders 416 to 418
The outputs of adders 416 to 417 are added to the next adder, respectively, up to the (n-1)th adder 418. By the way, the output signal from the multiplier 408 and the adders 416 to
The output signal from 417 is the last adder 4 in comparator circuits 419 to 421, respectively.
The output signals of 18 are compared with the signals after being attenuated by an attenuator 426, and the respective output signals O 1 - 423 - O o-1 425 are output to the output terminals 423 - 425.
Output. Note that the attenuation ratio of the attenuator 426 is the first
When calculating the reflected sound delay time, (0.1) 2 is used.
When determining the subsequent reverberation sound reverberation time, the value set to (0.001) 2 is used. This is similar to what is already known as an optimal delay time detector, and is used to find the delay time at which the value of the autocorrelation function becomes 1/10 or 1/1000. Also, from the delay circuit 401 to the delay circuit 4
The total delay time up to 03 is required to be about the length of the first reflected sound delay time and the subsequent reverberant sound reverberation time, that is, about 100 to 200 ms and 3 to 5 seconds, respectively. Naturally, as this time elapses, the integrating circuits 412-415 need to be reset. However, the reset circuit is not shown. Next, the priority encoder 430 shown in FIG. 13b is connected to the comparator 4 shown in FIG. 13a.
19~421 output signal O 1 423~ O o-1 42
5, and the delay time τ (=τ p ) corresponding to the attenuation ratio (0.1) 2 or (0.001) 2 of the attenuator 426 in FIG. 1 code corresponding to time)...
x Display or output o-1 . The following Figures 13c and d are diagrams showing other examples of circuits used as audible sound pressure analyzers, and in this case as well, the output signal φ 0 is obtained from the output terminal 422 in Figure 13a , you can know the audible sound pressure. This configuration is a combination of a low-pass filter 436 and a rectifier circuit 435 or a square multiplier 437, as shown in FIG. 13c. A circuit 519 for determining the interaural cross-correlation coefficient shown in FIG. 13e uses an example of a circuit system for determining the normalized correlation coefficient. First, input terminal 4
When the input signal Lin is input from 41 and the input signal Rin is input from the other input terminal 442, the input signals Lin and Rin directly enter the root mean square circuits 445 and 448, respectively, and output their respective root mean square values γ and δ. be done. Also the input signal
The sum signal between Lin and Rin is output from the adder 443,
The difference signal is taken out from an adder 444 that makes one signal negative and adds it to the other signal, and inputs it to root mean square circuits 446 and 447, and similarly, the sum signal and the difference signal are converted to respective root mean square values α and β.
Output. The four root mean square value signals α, β, γ, and δ thus obtained are sent to the calculation circuit 44.
9, the calculation process is (α 2 − β 2 )/(4γδ)
As a result, an output φ lr (o) corresponding to the cross-correlation coefficient, that is, an interaural cross-correlation coefficient (IACC) is obtained at the output terminal 450. However, when using the circuit shown in FIG. 13e, the input signals Lin and Rin are 0 to 1
If a delay circuit of ms is provided, the delay time is varied, and the maximum value is determined, the circuit configuration corresponds to the above-mentioned (Equation 9), and a value faithful to (Equation 9) can be obtained from that circuit. It will be done. A delay circuit for this purpose is not shown in FIG. 13e. Note that in order to calculate the interaural cross-correlation coefficient (IACC) with higher accuracy, numerical calculations may be performed using the above-mentioned equations (7) to (9). In the embodiment described so far, the inputs to the audible sound pressure analyzer 41, the first reflected sound delay time analyzer 42, and the subsequent reverberation sound reverberation time analyzer 43 are input to the dummy head or the left and right ears of the human head 1. Although the sum signal of two incoming acoustic signals has been used, the delay time between the two signals is usually less than 1 ms, and these three decomposers 4
Since there is no large difference in the analysis results of Nos. 1 to 43, sufficient analysis results can be obtained even if one of the left and right acoustic signals is used.
In the embodiment of the first invention described above, the optimum audible sound pressure [p] p of the sound field evaluation measuring instrument 10 varies depending on the type of sound source, but it is approximately (79±5). Since it is about dB, 79 dB is used here as a representative. Also, the optimal first
Reflected sound delay time [△t 1 ] p and optimal subsequent reverberant sound reverberation time [Tsub] p are the delay time τ that first becomes smaller than the attenuation ratio by the first reflected sound delay time analyzer 42 from the actual sound source itself. It is best to find p and use it as the optimal first reflected sound delay time [△t 1 ] p , or multiply it by (23±10) and use it as the optimal subsequent reverberation sound reverberation time [Tsub] p . In this example, [Tsub] p =23τ p was used. As described above, the sound field evaluation measuring instrument of the present invention includes one or two microphones, an amplification device for amplifying the acoustic signal obtained from the microphone, and an amplification signal that can be heard from the amplified signal. sound pressure p
audible sound pressure measuring means for measuring the audible sound pressure, and a first reflected sound delay time measuring means for measuring the first reflected sound delay time Δt 1 from the amplified signal; Reverberation sound reverberation time Tsub
a subsequent reverberation sound reverberation time measuring means for measuring the reverberation time; furthermore, a binaural cross-correlation coefficient measuring means for measuring the interaural cross-correlation coefficient IACC from the amplified signal; Optimum hearing sound pressure [ p ] If the above hearing sound pressure is greater than the optimal hearing sound pressure, - (0.07 ±0.03),
When it is small, the value obtained by multiplying by -(0.04±0.02) is output as the preference (psychological quality in terms of hearing) S 1 , that is, the above (Equation 14)
A first comparison conversion means that outputs an output corresponding to ~(Equation 16), and a ratio of the first reflected sound delay time △t 1 and the optimum first reflected sound delay time [△t 1 ] p to the base of 10. If the first reflected sound delay time is larger than the optimal first reflected sound delay time, set -(1.42±0.6) to the value of the absolute value of the logarithm to the 3/2 power, and if it is smaller, set -(1.11±0.5) to the value. The second function outputs the multiplied value as the preference measure S2 , that is, outputs corresponding to the above (Equation 17) to (Equation 19).
Comparison conversion means and the above subsequent reverberant sound reverberation time
The value is the 3/2 power of the absolute value of the logarithm to the base 10 of the ratio between Tsub and the optimal subsequent reverberation sound reverberation time [Tsub] p .
When the reverberation time of the subsequent reverberant sound is larger than the optimal reverberation time of the subsequent reverberation sound, the sum of -(0.74±0.25) times the sound pressure of the entire reflected sound and -(0.45±0.15) is used; Preference scale S3 In other words, the third function outputs the output corresponding to (Equation 20) to (Equation 22') above.
and outputs the value obtained by multiplying the value of the 3/2 power of the above-mentioned interaural cross-correlation coefficient by -(1.45±0.44) as the preference measure S4 , that is, the above-mentioned (Equation 23) ~ Comprehensively evaluates the preference scales S 1 , S 2 , S 3 , and S 4 output from the fourth comparison conversion means that outputs the output corresponding to (Equation 25) and the first to fourth comparison conversion means. Since the present invention is provided with comprehensive evaluation means for outputting the overall scale S of the preference, it is possible to perform an accurate and faithful sound field evaluation. (2) Application example of the present invention Next, based on the four physical parameters of the audible sound pressure in the sound field, the delay time of the first reflected sound, the reverberation time of the subsequent reverberant sound, and the interaural cross-correlation coefficient, DESCRIPTION OF THE PREFERRED EMBODIMENTS One embodiment of each of a first and second acoustic device to which the present invention is applied, which can correct a sound field and create a more preferable sound field, will be described with reference to the drawings. (2-1) First acoustic device FIG. 14a is a configuration block diagram of an embodiment of the first acoustic device that can correct the sound field and create a more preferable sound field. 14b shows one example of the sound field expansion device included in the audio device of FIG. 14a, and FIG. 14c shows another example of the reverberation device included in the audio device of FIG. 14a. First, in FIG. 14a, 501 is the input terminal INR , 502 is the input terminal INL , 503a to
503d is an adder, 504a, 504b are reverberation devices, 505a, 505b, 505r, 5
05l is an attenuator, 506a and 506b are delay circuits, 507 is a scaler which is a type of frequency dividing circuit, 508a and 508b are up-down counters, 509 is a digital-to-analog converter, and 510a and 510c are first reflected sound delay times. Analyzer, 510b is a subsequent reverberation sound reverberation time analyzer, 511a to 551c are encoders, 5
12a, 512b, 513a, 513b are comparators; 514a, 514b, 514c are smoothing circuits; 515a, 515b are variable resistors for setting target values; 516 is a sound field expansion device; 517r, 517l drive speakers. Power amplifier for 518r, 518l
is a speaker, 519 is an interaural cross-correlation coefficient (IACC) calculator, 520 is a multiplier, 531,
532 are input terminals inr, inl, 5 from the reverberation devices 504a, 504b to the sound field expansion device 516.
21 is the CLOCK signal input terminal, 533, 53
4 is the output terminal outr of the sound field expansion device 516;
outl, 535 is an input terminal attm from the comparator 513b to the sound field expansion device 516. 1st
4b, 530a to 530d are adders, 536 is an attenuator, 537 is a low pass filter, 538 is a high pass filter, and 53
9 is a phase shift circuit or a delay circuit. In FIG. 14c, 504c is a reverberator made by MR Schroeder and BS Atal, 540a to 540i are adders, and 54
1a, 451b are phase inversion circuits, 542a,
542b is voltage follow, 543a~
543f is a delay circuit, 544a to 544f are attenuators, 545 is an input terminal of attenuator 504c, and 546 is an output terminal of attenuator 504c. Note that the same reference numerals in the drawings as those in other drawings indicate the same or corresponding parts. Next, the operation of one embodiment of the first acoustic device will be explained using FIGS. 14a, 14b, and 14c. First, an acoustic signal from a record, tape, microphone, radio, etc. is applied to the input terminals INR 501 and INL 502. This acoustic signal is sent to the adder 5
03c and input to the first reflected sound delay time analyzer 510c. This first reflected sound delay time analyzer 510c is constituted by the circuit 42 in FIG. 12 described above, that is, the same circuit as that shown in FIG. 13a. The damping ratio is set to 1/10 so that time can be measured. The output of this first reflected sound delay time analyzer 511 is input to the encoder 511, and the encoder 51
1 outputs an output [(Δt 1 ] p = τ p ) obtained by setting the delay time [Δt 1 ] p of the first reflected sound that is optimal for the sound source as τ p . This encoder 51
The output from 1c is input to the comparator 512a, and is also input to the multiplier 520, 23
Multiplied by optimal subsequent reverberation sound reverberation time [Tsub] p
(=23°τ p ) is applied to the comparator 512b. On the other hand, the signals coming into the dummy head placed in a sound field such as a listening room or the microphones 2r and 2l installed at the left and right ears of the human head 1 are amplified by preamplifiers 3l and 3r. , the outputs are added by an adder 503d, and then input to the first reflected sound delay time analyzer 510a, which is the same as 42 in FIG.
It is also input to b. The configuration of this subsequent reverberation sound reverberation time analyzer 510b is the same circuit configuration as that described at 43 in FIG. 12. Next, the first reflected sound delay time analyzer 510a
The output is input to the encoder 511a, and a signal corresponding to the delay time Δt 1 of the first reflected sound after being affected by the sound field is obtained. Thereafter, it is input to the comparator 512a. Further, the output of the subsequent reverberation sound reverberation time analyzer 510b is input to the encoder 511b, and a signal corresponding to the subsequent reverberation sound reverberation time Tsub after being affected by the sound field is obtained, and this signal is sent to the comparator 5
12b. By the way, comparator 51
2a, the optimal first reflected sound delay time [Δt 1 ] p from the encoder 511c and the delay time Δt 1 of the first reflected sound after being affected by the sound field from the encoder 511a are compared. , if the first reflected sound delay time △t 1 is larger than the optimal first reflected sound delay time [△t 1 ] p , then in order to make the first reflected sound delay time △t 1 smaller, , comparator 512a to counter 5
Give countdown signal CD to 08a,
The content of the counter 508a is made smaller, and therefore, in the scaler 507 that performs frequency division based on the content of the counter 508a, the value to be divided becomes smaller, and the value from the input terminal 521 becomes smaller.
DELAY obtained by dividing the CLOCK signal
The frequency of the CLOCK signal increases and
The reverberation circuit 504 is activated by the DELAY CLOCK signal.
The delay time of the delay circuits 506a and 506b formed by BBDs (bucketed bridge devices) in the circuits a and 504b is reduced.
However, if the optimal first reflection delay time [△
t 1 ] If the first reflected sound delay time Δt 1 is smaller than p , a count-up signal CU is given from the comparator 512a to the counter 508a,
The contents of the counter 508a become larger, the value to be divided by the scaler 507 becomes larger, the frequency of the DELAY CLOCK signal becomes lower, and the delay time of the delay circuits 506a and 506b becomes larger. do. In addition, in the comparator 512b, the optimal subsequent reverberation sound reverberation time [Tsub] p of the multiplier 520
and the subsequent reverberant sound reverberation time Tsub after being given the influence of the sound field from the encoder 511b, and if the subsequent reverberant sound reverberation time Tsub is the optimal subsequent reverberant sound reverberation time
If Tsub is large, counter 508b
A count-up signal CU is sent out, and the contents of the counter 508b are converted to an analog value by the digital-to-analog converter 509, and then passed through the smoothing circuit 514c to the reverberation circuit 50
Attenuators 505 for each of 4a and 504b
The signals are applied to the attenuators 505a and 505b to increase the attenuation value in the attenuators 505a and 505b, thereby reducing the subsequent reverberation time Tsub. However, if the optimal subsequent reverberant sound reverberation time [Tsub] p , the subsequent reverberant sound reverberation time Tsub
is small, a countdown signal CD is sent to the counter 508b, and the contents of the counter 508b are converted to an analog value by the digital-to-analog converter 509, and then passed through the smoothing circuit 514c to the attenuators 505a and 50.
5b, and the signal causes the attenuator 5 to
The attenuation values at 05a and 505b become smaller, and therefore the reverberation time Tsub is increased, which is the opposite operation to the above case. By the way, the amplified signals from the preamplifiers 3l and 3r before being added are input to the binaural cross-correlation coefficient calculator 519, where the binaural cross-correlation coefficient IACC is calculated and converted into analog signals. and is input to the next smoothing circuit 514b. Next, in the comparator 513b, the smoothing circuit 514
Signal from b and interaural cross-correlation coefficient IACC
is compared with a target voltage signal from variable resistor 515b which sets a target value of . If the signal from the smoothing circuit 514b, that is, the signal corresponding to the binaural cross-correlation coefficient IACC, is larger than the target signal from the variable resistor 515b,
The output signal from the comparator 513b becomes smaller and is input to the input terminal att in5 of the sound field expansion device 516.
35, the sound field expansion device 516 operates to further reduce the binaural cross-correlation coefficient IACC. However, if the signal corresponding to the binaural cross-correlation coefficient IACC is smaller than the target signal from the variable resistor 515b, the output signal from the comparator 513b becomes larger and is input to the sound field expansion device. terminal att in53
5, the sound field expansion device 516 performs an operation opposite to the above case, such as increasing the binaural cross-correlation coefficient IACC. By the way, the first reflected sound delay time analyzer 510a, which receives the addition signal added by the adder 503d, outputs a sound pressure signal φ 0 corresponding to the audible sound pressure, and smoothes this sound pressure signal φ 0 . After passing through the circuit 514a, the signal is input to the comparator 513a. In this comparator 513a,
Variable resistor 515a for setting the target volume value
The target voltage signal from the smoothing circuit 51
If the output signal from the smoothing circuit 514a corresponding to the listening sound pressure is larger than the target volume value from the variable resistor 515a, the output from the comparator 513a is compared with the signal from the comparator 513a. becomes larger, and the attenuator 505r, 5
The attenuation of 05l increases and the listening sound pressure is lowered. Moreover, if the signal corresponding to the listening sound pressure is smaller than the target volume value, the output from the comparator 513a becomes smaller, and the attenuators 505r and 505
The operation is opposite to the above case, as the attenuation of l becomes smaller and the listening sound pressure is increased. Furthermore, the smoothing circuit 5 shown in Fig. 14a
14a, 514b, and 514c are all used to prevent noise from occurring due to sudden signal changes, and to prevent sudden changes in the delay times of delay circuits 506a and 506b in the reverberation circuits 504a and 504b. It is used in Next, the function of the sound field expansion device 516 is
A case will be explained in which the sound field enlarging device shown in FIG. 4b is used. First, the difference component between the two signals input from the reverberation devices 504a and 504b to the input terminals inr531 and inl532 is obtained by the adder 530a, and after passing through the attenuator 536 and the low-pass filter 537, crosstalk occurs. Adders 530b,
530c from input terminal inr531 to output terminal outr533 and input terminal inl532
to the output terminal outl534. Apart from that, the output of the low-pass filter 537 is further passed through a high-pass filter 538 and given a delay or transfer rotation in a phase shift or delay circuit 539, after which the crosstalking phase is also reversed. Adders 530d and 530e add the signal to each of the above-mentioned main paths with a phase such that
33 and outl 534, respectively. By changing the attenuator 536 of this sound field expansion device 516, the interaural cross-correlation coefficient IACC
can be changed. Further, this attenuator 536 is set so that the attenuation increases as the value of the signal from the input terminal att in 535, that is, the output signal from the comparator 513b increases. Furthermore, what is shown in Figure 14c is MR.
Reverberation device 504a, 504b shown in FIG. 14a with reverberation device 504c by Schroeder and BS Atal.
If this reverberator 504c is used instead of , a more natural feeling of reverberation can be obtained. In the above embodiment, it has been described that the microphones 2r and 2l are used to obtain the sound pressure at the dummy head or at the entrance of the external auditory canal of the human head 1, but when it is not necessary to obtain the value of the interaural cross-correlation coefficient IACC, may be performed using one microphone, or even when determining the interaural correlation coefficient IACC, it may be possible to use two microphones that are not attached to a dummy head or the like. Although the acoustic signal is shown as a two-channel stereo signal, it goes without saying that there may be cases where a two-channel signal is not necessary. As described above, in the first acoustic device, the input terminals INR , INL and the input terminals INR, INL ,
Optimal first reflected sound delay time measuring means for measuring the optimal first reflected sound delay time [△t 1 ] p from the acoustic signal coming in from IN L , and the optimal first reflected sound delay time [△t 1 ] p Optimum subsequent reverberation sound reverberation time output means that outputs an optimum subsequent reverberation sound reverberation time [Tsub] corresponding to (23±10) times p , and a sound field signal from a microphone placed in the sound field. Audible sound pressure measuring means for measuring audible sound pressure φ 0 , first reflected sound delay time measuring means for measuring a first reflected sound delay time Δt 1 from the sound field signal, and subsequent reverberation sound from the sound field signal. a subsequent reverberation sound reverberation time measuring means for measuring the reverberation time Tsub; an interaural cross-correlation coefficient measuring means for measuring the interaural cross-correlation coefficient IACC from the sound field signal; and the above-mentioned optimal first reflected sound delay time. [△t 1 ] p and the first reflected sound delay time △t 1 , and outputs a signal according to the difference; and the optimum subsequent reverberation sound reverberation time [Tsub] p and Above subsequent reverberation sound reverberation time Tsub
and a second comparison means that outputs a signal according to the difference, and the delay time changes depending on the output signal of the first comparison means, and the delay time changes depending on the output signal of the second comparison means. a reverberation means for changing the reverberation time and adding reverberation sound to the acoustic signal inputted from the input terminal; and an interaural cross-correlation coefficient setting means for setting a target interaural cross-correlation coefficient value in advance. and third comparison means for comparing the binaural cross-correlation coefficient IACC with the setting value of the binaural cross-correlation coefficient setting means and outputting a signal according to the difference, and the reverberation a sound field expanding means having an output capable of receiving the output signal of the means and changing the interaural cross-correlation coefficient IACC of the sound field according to the output signal of the third comparison means; a listening sound pressure setting means that can set a value of the listening sound pressure to be heard; 4 comparison means, attenuation means that receives the output of the sound field expansion means as input and can change the attenuation rate according to the output signal of the fourth comparison means, and amplifies the signal from the attenuation means. , an electroacoustic transducer that radiates an acoustic signal into the space, so it is not only possible to correct the sound field and create a more preferable sound field, but also to optimize the sound source within this first acoustic device. The first
The effect is that the reflected sound delay time [△t 1 ] p and the optimal reverberation time [Tsub] p can be obtained, and the values of the listening sound pressure and interaural cross-correlation coefficient IACC can be changed according to the user's preference. It is something that plays. (2-2) Second acoustic device Next, sound field devices for sound fields such as Japanese-style rooms with relatively few reflections and car interiors where even if there are many reflections, the reverberation time is extremely short, Correcting the sound field based on only three physical parameters: sound pressure, first reflected sound delay time, and subsequent reverberant sound reverberation time is sufficient to create a more favorable sound field. It goes without saying that the above-mentioned first acoustic device can be used as an acoustic device in the sound field targeted here, but the sound field used is a small room like this, a Japanese-style room with relatively little reflection, or a car. When confined indoors, the second acoustic device can not only obtain the same effect as the first acoustic device without using the interaural cross-correlation coefficient IACC parameter, but also It is more compact and inexpensive, and has a simpler circuit configuration than the first acoustic device, so it is easier to work with and has better durability than the first acoustic device. It can be effective. Therefore, the sound field is corrected using three physical parameters: the audible sound pressure in the sound field, the first reflected sound delay time, and the subsequent reverberation sound reverberation time, and this second method can create a more preferable sound field.
An embodiment of the acoustic device will be described with reference to the drawings. FIG. 15 is a block diagram of an embodiment of the second acoustic device. In the figure, 451 to 453 are absolute value conversion circuits, 511d is an encoder, 551 is an input terminal INR , 552 is an input terminal INL , and 553 is an integral reset signal input terminal. 554 is the first reflected sound delay time analyzer, which is shown in Fig. 13a mentioned above.
Although the circuit diagram and the internal configuration of the first reflected sound delay time analyzers 510a, 510c and the subsequent reverberant sound reverberation time analyzer 510b in FIG. Attenuator 426 in the autocorrelator configuration of
It is characterized by the addition of a comparator 419. In the case of such a configuration, the above-mentioned first
It operates similarly to the reflected sound delay time analyzers 510a, 510c or the subsequent reverberation sound reverberation time analyzer 510b. Note that the same reference numerals as those in FIG. 13 and FIG. 14a indicate the same or equivalent parts. First, acoustic signals are input from the input terminals INR 551 and INL 552, and the two signals are added by the adder 503c and input to the first reflected sound delay time analyzer 544. This first reflected sound delay time analyzer 554 performs the same operation as the first reflected sound delay time analyzers 510a and 510c in FIG. 14a,
The output therefrom is input to the encoder 511d as in the case of FIG. 14a described above, and the optimum first reflected sound delay time [△
t 1 ] p is output. An encoder 511d is applied to the scaler 507 so as to obtain a delay time corresponding to this optimal first reflected sound delay time [△t 1 ] p .
A signal is given from The scaler 507 divides the frequency of the CLOCK signal from the input terminal 521 to create a DELAY CLOCK signal, and outputs the signal to the delay circuits 506a and 5 of the reverberation circuits 504a and 504b.
A DELAY CLOCK signal is sent to 06b.
The audio signals are input terminals IN R 551 and IN L 552.
are input to reverberant circuits 504a and 504b respectively, and reverberant sound is added to the attenuator 505.
r, sent to 505l. On the other hand, the signal φ 0 corresponding to the sound pressure of the first reflected sound delay time analyzer 554 is passed through the smoothing circuit 514a to the comparator 5.
If the sound pressure is smaller than the target value, the output of the comparator 513a becomes smaller, and the output of the attenuator 513a becomes smaller. 505r and 505l are the comparators 513
It operates by reducing the amount of attenuation by the signal sent from a. However, if the sound pressure is greater than the target value,
The output of the comparator 513a increases, and the attenuators 505r and 505l increase the amount of attenuation by the signal sent from the comparator 513a, which is the opposite operation to the above case. Attenuator 505 thus obtained
The output of r, 505l is the power amplifier 517
r, 517l, speaker 518r, 5
Regenerated by 18l. At this time, the optimal subsequent reverberation sound reverberation time [Tsub] p , which is the optimal value of the reverberation time, is (23 ± 10) times the optimal first reflected sound delay time [△t 1 ] p , as stated earlier. This is desirable. It is known that when a reverberator is configured as the reverberation circuits 504a and 504b, the subsequent reverberation sound reverberation time Tsub=-ε·Δt 1 /log(g). Note that g is the attenuation rate. Here, by using the relationship Tsub=(23±10) Δt1 , it can be seen that the attenuation rate g must be in the range of 0.588 to 0.811. In addition, when the subsequent reverberation sound reverberation time Tsub = 23△t 1 , the attenuation rate g
=0.741. That is, the attenuation rate g
By setting the optimal first reflected sound reverberation time [△
t 1 ] p and the optimal subsequent reverberation sound reverberation time [Tsub] p are obtained. Of course, even with such processing, effects in the sound field will be added, but in spaces with relatively little reflection, such as Japanese-style rooms, or spaces with a lot of reflection, such as car cabins, the reverberation time is extremely short. is very effective. By the way, in this second audio device,
The present invention aims to provide a simpler acoustic device that has the same or better effect than the first acoustic device in a small room or a sound field with relatively little reflection, such as a Japanese-style room or a car interior.
When this second acoustic device is used in a sound field other than a small room or a vehicle interior, the performance is somewhat lower than that of the first acoustic device, but it still has a high audible sound pressure, first reflected sound delay time, and subsequent reverberation. Since the sound field is corrected using three physical parameters, namely the sound reverberation time, it is possible to create a much more favorable sound field than with conventional sound devices, and the second sound device is more effective than the first sound field. It is more compact, cheaper, and has a simpler circuit configuration than other devices.
It has the advantages of good workability and good durability. As described above, in this second audio device, the input terminals IN R and IN L and this input terminal
Optimal first reflected sound delay time measuring means for measuring the optimal first reflected sound delay time from the acoustic signals coming in from IN R and IN L ; and an audible sound pressure for measuring a signal corresponding to the audible sound pressure from the acoustic signals. a corresponding signal measuring means, and receiving a signal corresponding to the optimum first reflected sound delay time as a control signal;
The above acoustic signal is input, a delay time corresponding to the above optimum first reflected sound delay time is added to the above acoustic signal, and the subsequent reverberation sound reverberation time is (23±) of the above optimum first reflected sound delay time.
10) Reverberation means set to double, audible sound pressure setting means that can set a target audible sound pressure value in advance, and audible sound pressure that is the output of the audible sound pressure corresponding signal measuring means and the audible sound pressure. a comparison means that compares the set value of the setting means and outputs a signal corresponding to the difference; and an output signal of the reverberation means is inputted, and the output signal of the reverberation means is inputted from the reverberation means in accordance with the output signal of the comparison means. Since it is equipped with attenuation means that can change the attenuation rate given to the acoustic signal and electroacoustic conversion means that amplifies the signal from the attenuation means and radiates the acoustic signal into space, the sound field can be corrected. In addition to creating a more favorable sound field, an optimal first reflected sound delay time and an optimal subsequent reverberant sound reverberation time are obtained within this second acoustic device,
This has the effect that the listening sound pressure can be changed according to the user's preference. As described above, according to the present invention, the acoustic characteristics of the sound field in which music is listened to are measured, and the physical and psychological quantities of the listening sound pressure in the sound field, the first reflected sound delay time, the subsequent reverberant sound reverberation time, and the like are measured. This has the advantage that it is possible to provide a measuring instrument for evaluating the interaural cross-correlation coefficient using four parameters.

【図面の簡単な説明】[Brief explanation of drawings]

第1図は反射壁と音源と人頭との関係を示す
図、第2図a,bは実際の音楽を用いて測定した
正規化自己相関関数の例を示す図、第3図a,b
は音楽Aおよび音楽Bを用いて行つたプリフアレ
ンステストの結果を示す図、第4図はτdと〔△
t1pとの関係を示す図、第5図はτpと〔Tsub〕p
の関係を示す図、第6図は正規化されたプリフア
レンスとIACCの値の関係を示す図、第7図は相
対聴取音圧とプリフアレンスの尺度S1との関係を
示す図、第8図は直接音と第1反射音の間の遅れ
時間とプリフアレンスの尺度S2との関係を示す
図、第9図は後続残響音の残響時間とプリフアレ
ンスの尺度S3との関係を示す図、第10図は
IACCとプリフアレンスの尺度S4との関係を示す
図、第11図第12図はこの発明による音場評価
計測器の一実施例を示す概略の構成を示すブロツ
ク図、第13図a,bは2乗積分形残響計を用い
た第1反射音遅れ時間解析器あるいは後続残響音
残響時間解析器の基本構成を示した図、第13図
c,dは聴取音圧解析器の他の構成例を示す図、
第13図eは両耳間相互相関係数を求める回路の
構成例を示す図、第14図aはこの発明を応用し
た音響装置の一実施例の構成ブロツク図、第14
図bとcは音場拡大装置と残響器の構成図、第1
5図はこの発明を応用した音響装置の他の実施例
の構成ブロツク図である。 図において、1は人頭またはダミーヘツド、2
r,2lはマイクロフオン、3,3r,3lは前
置増巾器、4は物理量解析器、5は比較器、6は
心理量変換器、7は総合評価器、8は出力端子、
9は記録器、10は音場評価計測器、41は聴取
音圧解析器、42,510a,510c,554
は第1反射音遅れ時間解析器、43,510bは
後続残響音残響時間解析器、44は両耳間相互相
関関数解析器、51は聴取音圧心理量変換器、5
2は第1反射音遅れ時間心理量変換器、53は後
続残響音残響時間心理量変換器、54は両耳間相
互相関係数心理量変換器、451〜453は絶対
値化回路、501,502,521,531,5
32,535,551,552は入力端子、50
3a〜503dは加算器、504a,504bは
残響装置、505a,505b,505r,50
5lは減衰器、506a,506bは遅延回路、
507はスケーラ、508a,508bはアツプ
ダウンカウンタ、509はデイジタル−アナログ
変換器、511a〜511dはエンコーダ、51
2a,512b,513a,513bは比較器、
514a〜514cは平滑回路、515a,51
5bは可変抵抗器、516は音場拡大装置、51
7r,517lはパワーアンプ、518r,51
8lはスピーカ、519は両耳間相互相関係数計
算器、520は乗算器、533,534は出力端
子である。なお、図中、同一符号は同一又は相当
部分を示す。
Figure 1 is a diagram showing the relationship between a reflective wall, a sound source, and a human head. Figure 2 a, b is a diagram showing an example of a normalized autocorrelation function measured using actual music. Figure 3 a, b
Figure 4 shows the results of the preference test conducted using music A and music B.
t 1 ] p , Figure 5 is a diagram showing the relationship between τ p and [Tsub] p , Figure 6 is a diagram showing the relationship between normalized preference and IACC value, and Figure 7 is a diagram showing the relationship between τ p and [Tsub] p. Figure 8 shows the relationship between the relative audible sound pressure and the preference measure S1 , Figure 8 shows the relationship between the delay time between the direct sound and the first reflected sound and the preference measure S2, and Figure 9 shows the relationship between the delay time between the direct sound and the first reflected sound and the preference measure S2 . The figure shows the relationship between the reverberation time of the subsequent reverberant sound and the preference scale S 3 .
Figures 11 and 12 are diagrams showing the relationship between IACC and preference scale S4 , and Figures 13a and 13b are block diagrams showing a schematic configuration of an embodiment of the sound field evaluation instrument according to the present invention. A diagram showing the basic configuration of a first reflected sound delay time analyzer or a subsequent reverberant sound reverberation time analyzer using a square-integral type reverberation meter, and Figures 13c and d are other configuration examples of the audible sound pressure analyzer. A diagram showing
FIG. 13e is a diagram showing an example of the configuration of a circuit for calculating the interaural cross-correlation coefficient, FIG. 14a is a configuration block diagram of an embodiment of an acoustic device to which the present invention is applied,
Figures b and c are block diagrams of the sound field expansion device and reverberator, the first
FIG. 5 is a block diagram of another embodiment of an audio device to which the present invention is applied. In the figure, 1 is a human head or dummy head, 2
r, 2l are microphones, 3, 3r, 3l are preamplifiers, 4 is a physical quantity analyzer, 5 is a comparator, 6 is a psychological quantity converter, 7 is a comprehensive evaluator, 8 is an output terminal,
9 is a recorder, 10 is a sound field evaluation measuring device, 41 is a listening sound pressure analyzer, 42, 510a, 510c, 554
5 is a first reflected sound delay time analyzer, 43, 510b is a subsequent reverberation sound reverberation time analyzer, 44 is an interaural cross-correlation function analyzer, 51 is an audible sound pressure psychological quantity converter, 5
2 is a first reflected sound delay time psychological quantity converter; 53 is a subsequent reverberation sound reverberation time psychological quantity converter; 54 is a binaural cross-correlation coefficient psychological quantity converter; 451 to 453 are absolute value conversion circuits; 501; 502,521,531,5
32, 535, 551, 552 are input terminals, 50
3a to 503d are adders, 504a, 504b are reverberation devices, 505a, 505b, 505r, 50
5l is an attenuator, 506a and 506b are delay circuits,
507 is a scaler, 508a and 508b are up-down counters, 509 is a digital-to-analog converter, 511a to 511d are encoders, 51
2a, 512b, 513a, 513b are comparators,
514a to 514c are smoothing circuits, 515a, 51
5b is a variable resistor, 516 is a sound field expansion device, 51
7r, 517l are power amplifiers, 518r, 51
8l is a speaker, 519 is an interaural cross-correlation coefficient calculator, 520 is a multiplier, and 533 and 534 are output terminals. In addition, in the figures, the same reference numerals indicate the same or corresponding parts.

Claims (1)

【特許請求の範囲】 1 マイクロホンと、上記マイクロホンの音響信
号を増巾するための増巾器と、上記増巾器からの
増巾信号から聴取音圧を計測するための聴取音圧
計測手段と、上記増巾信号から第1反射音遅れ時
間を計測するための第1反射音遅れ時間計測手段
と、上記増巾信号から後続残響音残響時間を計測
するための後続残響音残響時間計測手段と、上記
増巾信号から両耳間相互相関係数を計測するため
の両耳間相互相関係数計測手段と、上記聴取音圧
と最適聴取音圧との比の10を底とする対数の20倍
の値の絶縁値の3/2乗の値に、上記聴取音圧が最
適聴取音圧より大きいときには−(0.07±0.03)
を、小さいときには−(0.04±0.02)を乗じてな
る値をプリフアレンス(聴感上の心理良さ)の尺
度S1として出力する第1の比較変換手段と、上
記第1反射音遅れ時間と最適第1反射音遅れ時間
との比の10を底とする対数の絶対値の3/2乗の値
に、上記第1反射音遅れ時間が最適第1反射音遅
れ時間より大きいときには−(1.42±0.6)を、小
さいときには−(1.11±0.5)を乗じてなる値をプ
リフアレンスの尺度S2として出力する第2の比較
変換手段と、上記後続残響音残響時間と最適後続
残響音残響時間との比の10を底とする対数の絶対
値の3/2乗の値に、上記後続残響音残響時間が最
適後続残響音残響時間より大きいときには、反射
音全体の音圧の−(0.74±0.25)倍の値と−(0.45
±0.15)との加算値を、小さいときには反射音全
体の音圧の−(0.42±0.14)倍の値と−(2.36±
0.79)との加算値を乗じたのちの値が負なるとき
にはそのままの値を、正なるときは零をプリフア
レンスの尺度S3として出力する第3の比較変換手
段と、上記両耳間相互相関係数の3/2乗の値に−
(1.45±0.44)を乗じた値をプリフアレンスの尺
度S4として出力する第4の比較変換手段と、上記
のプリフアレンスの尺度S1,S2,S3,S4を総合評
価してプリフアレンスの全体の尺度Sを出力する
総合評価手段とを具備したことを特徴とする音場
評価計測器。 2 マイクロホンはダミーヘツドまたは人頭の外
耳道入口に取付けられたことを特徴とする特許請
求の範囲第1項記載の音場評価計測器。 3 最適第1反射音遅れ時間と最適後続残響音残
響時間とは第1反射音遅れ時間計測手段におい
て、それぞれ増巾信号の正規化自己相関関数が
0.1となる時間と、その(23±10)倍の時間とに
なるように設定したことを特徴とする特許請求の
範囲第1項または第2項記載の音場評価計測器。 4 総合評価手段において、プリフアレンスの尺
度S1とS2とS3とS4とを加算した値をプリフアレン
スの全体の尺度Sとして出力することを特徴とす
る特許請求の範囲第1項ないし第3項のいずれか
に記載の音場評価計測器。
[Claims] 1. A microphone, an amplifier for amplifying the acoustic signal of the microphone, and an audible sound pressure measuring means for measuring the audible sound pressure from the amplified signal from the amplifier. , a first reflected sound delay time measuring means for measuring a first reflected sound delay time from the amplified signal, and a subsequent reverberant sound reverberation time measuring means for measuring a subsequent reverberant sound reverberation time from the amplified signal. , an interaural cross-correlation coefficient measuring means for measuring an interaural cross-correlation coefficient from the amplified signal, and a logarithm of 20 to the base 10 of the ratio between the audible sound pressure and the optimum audible sound pressure. When the above hearing sound pressure is greater than the optimum hearing sound pressure, the value is the 3/2 power of the insulation value, which is twice the value - (0.07±0.03)
is multiplied by -(0.04±0.02) when it is small, and outputs the value obtained by multiplying the value by -(0.04±0.02) as a measure of preference (psychological quality in terms of hearing) S1, and the first reflected sound delay time and the optimal first reflection. If the first reflected sound delay time is greater than the optimal first reflected sound delay time, set -(1.42±0.6) to the 3/2 power of the absolute value of the logarithm to the base 10 of the ratio to the sound delay time. , when it is small, a value obtained by multiplying by -(1.11±0.5) is output as the preference scale S2 , and 10 of the ratio of the above-mentioned subsequent reverberant sound reverberation time and the optimal subsequent reverberant sound reverberation time. When the reverberation time of the subsequent reverberant sound is larger than the optimal reverberation time of the subsequent reverberant sound, the value of the absolute value of the logarithm to the power of 3/2 is -(0.74±0.25) times the sound pressure of the entire reflected sound. −(0.45
±0.15) and -(0.42±0.14) times the total sound pressure of the reflected sound when it is small, and -(2.36±
0.79), and outputs the same value as it is when the value is negative, and outputs zero when it is positive as the preference scale S3 , and the above-mentioned interaural mutual correlation. − to the value of the 3/2 power of the number
(1.45±0.44) and outputs the value multiplied by (1.45±0.44) as the preference scale S 4 , and comprehensively evaluates the above preference scales S 1 , S 2 , S 3 , and S 4 to obtain the overall A sound field evaluation measuring instrument characterized by comprising: a comprehensive evaluation means for outputting a scale S of . 2. The sound field evaluation measuring instrument according to claim 1, wherein the microphone is attached to a dummy head or to the entrance of the external auditory canal of a human head. 3. The optimal first reflected sound delay time and the optimal subsequent reverberant sound reverberation time are determined by the normalized autocorrelation function of the amplified signal in the first reflected sound delay time measuring means.
The sound field evaluation measuring instrument according to claim 1 or 2, characterized in that the time is set to be 0.1, and the time is set to be (23±10) times that time. 4. Claims 1 to 3, characterized in that the comprehensive evaluation means outputs the sum of the preference scales S1 , S2 , S3, and S4 as the overall preference scale S. A sound field evaluation measuring instrument as described in any of the paragraphs.
JP9297383A 1983-05-24 1983-05-24 ONBAHYOKAKEISOKUKI Expired - Lifetime JPH0237971B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP9297383A JPH0237971B2 (en) 1983-05-24 1983-05-24 ONBAHYOKAKEISOKUKI

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP9297383A JPH0237971B2 (en) 1983-05-24 1983-05-24 ONBAHYOKAKEISOKUKI

Related Child Applications (2)

Application Number Title Priority Date Filing Date
JP2016968A Division JPH039226A (en) 1990-01-26 1990-01-26 Acoustic device
JP2016969A Division JPH039227A (en) 1990-01-26 1990-01-26 Acoustic device

Publications (2)

Publication Number Publication Date
JPS59216022A JPS59216022A (en) 1984-12-06
JPH0237971B2 true JPH0237971B2 (en) 1990-08-28

Family

ID=14069345

Family Applications (1)

Application Number Title Priority Date Filing Date
JP9297383A Expired - Lifetime JPH0237971B2 (en) 1983-05-24 1983-05-24 ONBAHYOKAKEISOKUKI

Country Status (1)

Country Link
JP (1) JPH0237971B2 (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5658483B2 (en) * 2010-05-26 2015-01-28 日本放送協会 Reality estimation device and program thereof
JP5150004B1 (en) * 2012-08-09 2013-02-20 リオン株式会社 Noise observation apparatus and noise observation method
CN105675125B (en) * 2016-01-31 2019-02-01 柳州市展虹科技有限公司 A kind of complicated stable sound field acoustic pressure test device intelligence control system

Also Published As

Publication number Publication date
JPS59216022A (en) 1984-12-06

Similar Documents

Publication Publication Date Title
JP4726875B2 (en) Audio signal processing method and apparatus
KR101370365B1 (en) A method of and a device for generating 3D sound
JP5524237B2 (en) Method and apparatus for applying echo to multi-channel audio signals using spatial cue parameters
EP0762804B1 (en) Three-dimensional acoustic processor which uses linear predictive coefficients
US20050216211A1 (en) Impulse response collecting method, sound effect adding apparatus, and recording medium
US7286946B2 (en) Transmission characteristic measuring device transmission characteristic measuring method, and amplifier
JPH02280199A (en) Reverberation device
JP2014529101A (en) Parametric encoder for encoding multi-channel audio signals
JP4234103B2 (en) Apparatus and method for determining impulse response and apparatus and method for providing speech
JP5627440B2 (en) Acoustic apparatus, control method therefor, and program
JPH0237971B2 (en) ONBAHYOKAKEISOKUKI
JPH0454439B2 (en)
JPH07334181A (en) Sound reverberation generating device
JP5801405B2 (en) Estimation of synthesized speech prototypes
JPS63234699A (en) Sound field correcting device
JPH0453480B2 (en)
Soeta et al. Spatial audio quality factors for headphone evaluation
JPH09327100A (en) Headphone reproducing device
JPH04295727A (en) Impulse-response measuring method
JPH02193500A (en) Correcting method for acoustic frequency characteristic
JPH07203597A (en) Headphone reproducing device
JPH05268693A (en) Sound field reproduction method
JPS63157600A (en) Sound field correction system
JPS6110182Y2 (en)
JPH0894731A (en) Method and device for detecting direction of sound source