Transfer Function Method of Measuring In-Duct Acoustic Properties. I. Theory
Transfer Function Method of Measuring In-Duct Acoustic Properties. I. Theory
Theory
J. Y. Chung and D. A. Blaser
Citation: The Journal of the Acoustical Society of America 68, 907 (1980); doi: 10.1121/1.384778
View online: http://dx.doi.org/10.1121/1.384778
View Table of Contents: http://asa.scitation.org/toc/jas/68/3
Published by the Acoustical Society of America
A transfer-matrix approach for estimating the characteristic impedance and wave numbers of limp and rigid
porous materials
The Journal of the Acoustical Society of America 107, (2000); 10.1121/1.428404
Influence of errors on the two‐microphone method for measuring acoustic properties in ducts
The Journal of the Acoustical Society of America 79, (1998); 10.1121/1.393542
Transfer function method for measuring characteristic impedance and propagation constant of porous
materials
The Journal of the Acoustical Society of America 86, (1998); 10.1121/1.398241
On the multiple microphone method for measuring in-duct acoustic properties in the presence of mean flow
The Journal of the Acoustical Society of America 103, (1998); 10.1121/1.421289
Transfer function method of measuring in-duct acoustic
properties. I. Theory
J. Y. Chungand D. A. Blaser
EngineeringMechanicsDepartment,GeneralMotorsResearchLaboratories,Warren,Michigan48090
(Received10 December1979;acceptedfor publication5 June1980)
ACOUSTIC MICROPHONE
#1 MICROPHONE#2 and
DRIVER
907 J. Acoust. Soc. Am. 68(3), Sept. 1980 0001-4966/80/090907-07500.80 (D 1980 AcousticalSocietyof America 907
the incident and the reflected waves, respectively,
evaluated between the two microphone locations. .= R•i,i(r/)
e'•"'½"
dr/ [h,,.(r)
-hi(r)]e' •"'½*dr
h,2: the impulsive response correspondingto the
combined incident and reflected waves evaluated be- : oo
a,• ,•(r/)e'•"'•drl 0 [hr(r)-h,,.(r)]e'•'-•'dr,
tween the two microphone locations.
(14)
It should be noted that for convenience, the system or
input and output implicit in Eq. (6)was made opposite
to the actual physical condition (see Fig. I for the S•,,•(f)f •[h•,.(r)
-h,(r)]e
-'"•'r
dr
order of the microphone locations relative to the di- s,,•,'(f):f;[h•(r)-h•,.(r)]e
-•"•dr ' (15)
rection of propagation of the reflected wave). It can
be shown, however, that the impulsive response so whereSt• t•(f), the Fourier transform of R,• t•(r/) is the
auto-spectral density of the incident pressure compo-
defined (i.e., hr) is related to that obtainableby re-
versing the system input and output through a simple nent at the first microphone. Similarly, St• tr(f) is the
cross spectral density between the incident and reflect-
convolution integral i.e., the convolution of the two
quantities is a delta function. ed pressure components at the same microphone loca-
tion. The quantity at the left side of Eq. (15) is de-
We find from Eqs. (1), (2), and (7) that fined as the complex reflection coefficient* at the first
microphone location, or
p•.•(t)+p•,(t) = h•.(r)[p•i(t - •) +p•r(t- •)]d•. (8)
(16)
Using the relations in Eqs. (5) and (6), Eq. (8) be- Using the relation in Eq. (3), one has
comes,
E{Pt,(t)'Pt•(t+r)}
= rt(rl)E{Pti(t)Pt,(t+r-rl)}drl.
(17)
foøøp•i(t-
r)[h•,.(r)-hi(r)]dr It can be shown from Eq. (17) that
=fo•p,r(t-
r)[h•(r)-h,,.(r)]dr. (9) S•i•(f)= r,(rl)e.•,.,•,r
S• •(f)
I' (18)
Multiplying both sides of Eq. {0) by p,•(t- a), then tak- Thusfrom Eqs. (16) and (18), it is seenthatR,(f) is
ing expectedvalues6 we have the following relation: the Fourier transform of r,(r/), or
o•E{P,i(t
- r)p,•(t
- a)}[h,,.(r)
- )= •o
Hr(f hr(r)
e'•2'/*
o
St•,(f)
dr=str
tr(f), (21)
and
=fo
•E{P,r(t
- r)p,i(t
- a)}[h•{,)-
h,,.(r)]dr. (lO)
H,:(f): h,,.(r)
e'•"'•'dr:
Sz•.(f) (22)
But, S,,(f) '
where Hi, Hr, and Hxaare the acoustical transfer func-
E{p,•(t - r)p•i(t- a)}=R,i•(a- r) (11)
tions corresponding to the impulsive responses hi, hr,
and
and hx,., respectively.
E{p•(t- r)P•i(t- a)}=R•i•(a- r), (12) From Eqs. (15) to (22) the expression of the complex
where R,i,i and R,•r are the auto- and cross-correla- reflection coefficient becomes,
tion functions, respectively. SubstitutingEqs. {11) R•.(f)= [H,:(f) - Hi(f) ]/[Hr(f) -H•.•.(/)]. (23)
and (19.)into Eq. (10)andtakingFourier transforms,
Equation (23) is the basis of the new transfer function
we obtain,
technique for determining the acoustic impedance and
other acoustic parameters. From the definition in Eq.
f.•ofo•OR,i
,,(a-
r)[h,,.(r)-hi(r)]
e'•"•drda (16), the magnitudeof Rt representsthe amplituderatio
of the reflected and incident pressure components at the
first microphone location.
= It• ,•(a- r)[h•(r)- h•,.(r)]e'•"*• drda.
The relation between R• and the reflection coefficient
at the second microphone, R,. can be found from the
following relations:
Letting a- r =r/and da=dr/, then Eq. (13) can be re-
written as, From Eqs. (4)and (6), it is shown,
908 J.Acoust.
Soc.Am.,Vol.68,No.3,September
1980 J.Y.Chung
andD.A.Blaser:
In-duct
acoustic
properties.
I 908
r•.(r)p•.i(t- r)dr= hr(r)P•r(t- r)dr. (24)
o•ø{p•i(t
- r)[hr(r)
-h•2(r)]
+P•,(t
- r)
Substituting Eqs. (3)and (5) into Eq. (24), one has x[h,(r)-h•2(r)]
-p•(t- r)[h,O')-hi(r)]}dr=O, (32)
fo
• fo•r2(r)hi(•)P•i(t-
r-•)d•dr or
R•i •i(•) e'•2'•rd• Multiplying both sides of Eq. (34) by p•(t- a) and taking
expected values, we have, in terms of auto- and cross-
correlation functions,
x hi(•) e'j2'•td• = Rxi•(•7)e'•2'sc
o•R,,(e
- r)[h•(r)-h•2(r)]dr
x rx(•)e'•2'S•d• h•(r)e'•2'•*dr. (27) = R•i(a- r)[hr(r)-h•(r)]dr. (35)
Thus,
Taking the Fourier transform of Eq. (35), one shows,
Sxixi(f)'R2(f)'Hi(f)=Sx• x,(f)'Rx(f)'Hr(f) Hr -H•2 =(S• •i/S•)(Hr-Hi). (36)
A similar derivation leads to
or
u,(f) =e (30)
s =m(X/2), rn = 1, 2, 3, .... (39)
where k is the wavenumber and s is the microphone
Equation (39) indicates that the reflection coefficient
spacing. Eqs. (28) to (30) show that the magnitude of
cannot be determined from Eq. (23) at discrete fre-
the reflection coefficient is independent of the location
quency points for which the microphone spacing is an
at which it is measured. Physically, this should be
integer multiple of the half-wavelength of sound. In
obvious since the magnitude of the reflection coefficient
order to avoid these points up to a frequency fro, the
represents the ratio of the reflected and the incident
microphone spacing s must be chosen such that
sound power components in the tube. Provided the
losses at the tube wall are negligible, these components
[accordingto Eq. (39)],
are not functions of position along the tube. s •<c/2f.,, (40)
The expression of the reflection coefficient in Eq. (23) where c is the speed of sound. In order to determine
becomes indeterminate when H,-H•2 = 0. This singular the reflection coefficient R on the surface of a test
condition can be further investigated as follows: material which is not at the microphone location, Eq.
fo•ø{p•(t-
r)[h•2(r)-
h•
(r)]-p•,(t-
r)
,
where
(44)
1+R.I_201Ogxol
[I+R•.
=20loglo Ht[
+101Oglo(•)
=20loglo
H,-//•2- 201Oglol
Ht[+10loglo ,
II. THE COMPLEX ACOUSTIC IMPEDANCE (51)
Using Eqs. (23) and (44) and the standard relation where H u 12and H a12are the transfer functions measured
between the complex acoustic impedance and the re- at the upstream and downstream locations, respectively,
flection coefficient, one can obtain the complex acoustic also
impedance on the cross-sectional area of the tube at a
distance l from the first microphone, or
z/pc = (1 +R)/(1-R) is the transfer function of the acoustic element, which
canbe determinedindependently
from H •2andH •.,
Hi•sin(k/) sin[k(/-s)] (45) within the range of linear acoustic propagation. It
=Jcos[k(/-s•
- H12
cos(kli" shouldbe noted that the final form for TL in Eq. (51) is
with its real and imaginary parts being, valid only if the nondimensional wavenumberks [see
Eqs. (29) and (30)] is the same for both upstream and
x/pc=(Re(HI•.)
sin[k(2/- s)] - 1/2{sin[2k(l- s)] downstream' measurements.
+I sin(2kl)})/Ha, (46)
and
w, =s,, .A,/(pcl 1 +Ra12). (50) Thus usingRt in Eq. (23), the soundabsorptioncoef-
ficient shown in Eq. (56) can be rewritten as,
The transmission loss can then be defined using Eqs.
(49) and (50) as, a = 1- - . (sv)
912 J. Acoust.
Soc.Am.,Vol.68, No.3, September
1980 J.Y. Chung
andD.A. Blaser:In-duct
acoustic
properties.
I 912
SHerethe auto-spectrumis definedas the auto spectral density ence Function Method," J. Acoust. Soc. Am. 62, 388-395
multiplied by the frequency resolution bandwidth. (1977).
sJ.Y. Chung,"Cross-SpectralMethodof MeasuringAcoustic l•Thecross spectraldensityS•2in Eq. (66)is notcontaminated.
Intensity Without Error Caused by Instrument Phase Mis- Also the contaminations of the auto-spectral densities are
match," J. Acoust. Soc. Am. 64, 1613-1616 (1978). assumed to be mutually uncorrelated.
10j. y. Chung, "The Rejectionof Flow Noise Usinga Coher-
913 J. Acoust.
Soc.Am.,Vol.68, No.3, September
1980 J.Y. ChungandD. A. Blaser:In-ductacoustic
properties.
I 913