11.07.2014 Views

Flute acoustics: measurement, modelling and design - School of ...

Flute acoustics: measurement, modelling and design - School of ...

Flute acoustics: measurement, modelling and design - School of ...

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

CHAPTER 3. MEASURING ACOUSTIC IMPEDANCE 38<br />

subsequent row in (3.4) by the first row yields the first column <strong>of</strong> A,<br />

A j 1 = A 11 b j /b 1 , (3.7)<br />

in terms <strong>of</strong> A 11 . A 11 can be given any value without affecting impedance <strong>measurement</strong>s; it is<br />

usually set equal to cosh(ikx 1 ) (equivalent to assuming the first microphone has unity gain <strong>and</strong><br />

that the <strong>measurement</strong> duct is cylindrical).<br />

Taking pairs <strong>of</strong> rows from (3.5), we may eliminate p ′ <strong>and</strong> obtain linear equations in the unknowns<br />

A j 2 . For example, rows 1 <strong>and</strong> 2 combine to give b ′ 2 (A 11 + A 12 / ¯Z ′ ) = b ′ 1 (A 21 + A 22 / ¯Z ′ ).<br />

The elements A j 2 are determined by eliminating the pressure (p ′ or p ′′ ) from each pair <strong>of</strong> rows<br />

in (3.5) <strong>and</strong> (3.6) <strong>and</strong> solving the resulting system. Note that for all n > 2 the system is overdetermined<br />

(in the algebraic, noise-free sense) <strong>and</strong> for n = 2 is algebraically equivalent to the TMTC<br />

technique <strong>of</strong> Gibiat & Laloë.<br />

So, e.g. for a head with two microphones, calibrated with three known loads,<br />

[<br />

b<br />

′<br />

2<br />

/ ¯Z ′ −b ′ 1 / ¯Z ′<br />

][ ]<br />

A12<br />

b<br />

2 ′′/<br />

Z ¯′′<br />

−b<br />

1 ′′/<br />

Z ¯<br />

=<br />

′′ A 22<br />

[<br />

b<br />

′<br />

1<br />

A 21 − b ′ 2 A 11<br />

b ′′<br />

1 A 21 − b ′′<br />

2 A 11<br />

]<br />

. (3.8)<br />

The above-outlined calibration technique assumes very little about the geometry <strong>of</strong> the impedance<br />

head <strong>and</strong> the characteristics <strong>of</strong> the transducers. If the calibration is complete, then<br />

wall losses within the impedance head do not need to be taken into account explicitly.<br />

the multiple microphone technique with cylindrical waveguide <strong>and</strong> microphones attached at<br />

known distances from the reference plane, the calibration parameters may be recast in a more<br />

instructive form. If each microphone has an admittance <strong>of</strong> y j /Z 0 , then the pressure <strong>and</strong> upstream<br />

flow at microphone j are related to those at microphone j − 1 according to<br />

[<br />

p j<br />

Z 0 U + j<br />

]<br />

=<br />

[<br />

1 0<br />

−y j 1<br />

]<br />

T<br />

[<br />

p j −1<br />

Z 0 U + j −1<br />

]<br />

In<br />

, (3.9a)<br />

where T is the transfer matrix for a cylindrical pipe<br />

T =<br />

[<br />

cosh(ikd j ) sinh(ikd j )<br />

sinh(ikd j ) cosh(ikd j )<br />

]<br />

(3.9b)<br />

<strong>and</strong> d j = x j − x j −1 . For j = 1, p j −1 <strong>and</strong> U + are the pressure <strong>and</strong> flow at the reference plane<br />

j −1<br />

<strong>and</strong> d 1 = x 1 . The microphone signals are equal to κ j p j , where κ j is the gain <strong>of</strong> microphone<br />

j . Taking (3.9) <strong>and</strong> a calibrated matrix A, one can determine k, κ j for j = 1,...,n <strong>and</strong> y j for<br />

j = 1,...,n −1 (y n , the dimensionless admittance <strong>of</strong> the microphone closest to the source, cannot<br />

be determined). For a <strong>measurement</strong> set-up with a combination <strong>of</strong> pressure <strong>and</strong> flow transducers,<br />

or all flow transducers, the calibration proceeds in much the same way <strong>and</strong> an equation<br />

similar to (3.9) can be constructed.<br />

Sometimes a third calibration is unnecessary or impracticable. In these cases one may precalculate<br />

the complex wavenumber k, using a theory that accounts for viscothermal losses<br />

within the waveguide. For a given impedance head, a single set <strong>of</strong> three calibrations can determine<br />

the degree <strong>of</strong> confidence one may take in this assumption, <strong>and</strong> the errors involved in<br />

making it. The remaining elements <strong>of</strong> A are then found from (3.4) <strong>and</strong> (3.5) as described.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!