Page  00000001 Digital Waveguide Modeling of W~oodwind Toneholes Gary P. Scavone garyeccrma.stanford.edu Julius 0. Smith, HII j oseccrma.stanford.edu Center for Computer Research in M~usic and Acoustics (OCR MA ) Department of M1~usic, Stanford University Stanford, Califonmia 94305 USA Ab str act This paper demonstrates a digital waveguide implementation of a six-hole woodwind tonehole lattice, based on theory and measurements published by Keefe (1981). Woodwind tonehole transmission-matrix parameters are converted to traveling-wave scattering parameters suitable for digital waveguide implementation. Second-order digital filters are designed to approximate the reflection and transmission transfer functions implied by the Keefe data. In this way, the tonehole is implemented by a two-port scattering junction which accounts for both series and shunt complex impedances. Alternatively, the tonehole can be implemented as a one-multiply, one-filter, threeport scattering junction. The results of a digital waveguide six-hole flute bore implementation using both models are compared to Keefe (1990), with excellent agreement. In this way, the best available acoustic theory regarding toneholes is efficiently and accurately simulated in discrete-time. 1 Two-Port Tonehole M~odel The fundamental acoustic properties of toneholes have been extensively studied and reported by Keefe (1981, 1990). The model described by Keefe (1990) is an accurate representation for a tonehole unit, assuming adjacent tonehole interactions are negligible. In this description, acoustic variables at the tonehole junction are related by a transmission matrix of series and shunt impedance parameters. K~eefe's original derivation of the tonehole parameters was based on a symmetric T section, as shown in Fig. 1 (Keefe, 1981). The series impedance terms, Za, result from an analysis of symmetric pressure distribution, or a pressure anti-node, at the tonehole, for which pressure is symmetric and equal across the junction. The transmission matrix which results under this analysis is given by [ J~] = 1+ 2Z, z-1 Za (i + tZ) [P2 (1) obtained by cascading the three matrices which correspond to the three impedance terms. Based on the approximation that IZa/ZJI < 1, Eq. (1) can be reduced to the form [ Lp:]=[Zj;' ~][~]: (2) which is the basic tonehole unit cell given by K~eefe for transmission-matrix calculations. The values of Za and Z5 vary according to whether the tonehole is open (o) or closed (c) as Fig. 1: T section transmission-line representation of the tonehole. an analysis of anti-symmetric pressure distribution, or a pressure node, at the tonehole junction. In this Z~'= Zo(a/b)2 (jkt6 + &) Zb)= -jZo(a/b)2cot(ktj, 4jO) = joab2t) z&") = -~~l)~~ (3a) (3b) (3c) (3d) case, volume flow is symmetric and equal across the Definitions and descriptions of the various parameters junction. The shunt impedance term, Z5, results from in Eqs. (3a) (3d) can be found in (Keefe, 1990).

Page  00000002 To render these relationships in the digital waveguide domain, it is necessary to transform the planewave physical variables of pressure and volume velocity to traveling-wave variables as + + P_ P, - Pý 1 - (4) Ui -1 Z P+- P~ ( where Zo is the characteristic impedance of the cylindrical bore, which is equal on both sides of the tonehole. Waveguide pressure variables on both sides of the tonehole are then related by P[721 2(5) + where+ where method (Smith, 1983, p. 46), in conjunction with the equation-error method. That is, a one-pole model H (z) was first fit to the continuous-time response, H(eJi). Subsequently, the inverse error spectrum, H](eIJ)/H(eIJ) was modeled with a two-pole digital filter, H2 (z). The discrete-time approximation to H(eJi) was then given by H^i (z) /H2(z). The first step of this design process captures the peaks of the spectral envelope, while the second step models the "dips" in the spectrum. These particular calculations were performed for a tonehole of radius b = 4.765 mm, minimum tonehole height t, = 3.4 mm, tonehole radius of curvature re = 0.5 mm, and air column radius a = 9.45 mm. The results of Keefe (1981) were experimentally calibrated for frequencies less than about 5 kHz, so that the continuous-time responses evident in the figures are purely theoretical above this limit. Therefore, the discrete-time filter design process was weighted to produce better matching at low frequencies. n- = V+ ZaZs - Z, ZaZs + 2ZoZs + Z' 2ZoZs ZaZs + 2ZoZs + Z2 ' calculated using Eqs. (1) and (4) and then making appropriate simplifications for |Za/Zs| <K 1. Figure 2 depicts the waveguide tonehole two-port scattering junction in terms of these reflectances and transmittances. This structure is analogous to the fourmultiply Kelly-Lochbaum scattering junction (Kelly and Lochbaum, 1962). Frequency (kHz) 51 1 i----.----L--------4------- p+(n nn Open Hole - S Closed Hole - - Continuous-Time Responses - Discrete-Time Responses Pi Fig. 2: Digital waveguide tonehole two-port scattering junction. For the implementation of the reflectances and transmittances given by Eqs. (6a) - (6b) in the digital waveguide structure of Fig. 2, it is necessary to convert the continuous-time filter responses to appropriate discrete-time representations. In this study, use is made of an equation-error minimization technique (Smith, 1983) which matches both frequency response magnitude and phase. This technique is implemented in MATLAB~ by the function invfreqz. Figure 3 plots the responses of second-order discretetime filters designed to approximate the continuoustime magnitude and phase characteristics of the reflectances for closed and open toneholes. The openhole discrete-time filter was designed using Kopec's -51 0 2 4 6 8 10 12 14 16 18 20 Frequency (kHz) Fig. 3: Two-port tonehole junction closed-hole and open-hole reflectances, derived from Keefe (1981) shunt and series impedance parameters. (top) Reflectance magnitude; (bottom) Reflectance phase. Figure 4 plots the reflection function calculated for a six-hole flute bore, as described in (Keefe, 1990). The upper plot was calculated using Keefe's frequency-domain transmission matrices, such that the reflection function was determined as the inverse Fourier transform of the corresponding reflection coefficient. This response is equivalent to that provided by Keefe (1990), though scale factor discrepancies exist due to differences in open-end reflection models and lowpass filter responses. The lower plot was calculated from a digital waveguide model using twoport tonehole scattering junctions. Differences between the continuous- and discrete-time results are

Page  00000003 0.15 0.1 -- 0.05 - 0 0.05 - -0.1 0 1 2 3 4 5 6 7 8 9 1( Time (ms) be modeled by another waveguide and the reflection/transmission characteristics at its end by an appropriate digital filter. This tonehole model is then attached to the appropriate branch of the three-port junction. The bore characteristic admittance Yo is equal on either side of the junction, while the real tonehole characteristic admittance is Yoth. The three-port scattering junction equations for pressure traveling-wave components can be determined as Pa(t) = rop(t) + [ + ro] P(t) -2roph(t) (7a) 0 1 2 3 4 5 6 7 8 9 10 Time (ms) Fig. 4: Reflection functions for note G (three finger holes closed, three finger holes open) on a simple flute [see (Keefe, 1990)]. (top) Transmission-line calculation; (bottom) Digital waveguide two-port tonehole implementation. most apparent in early, high-frequency, closed-hole reflections. The continuous-time reflection function was low-pass filtered to remove time-domain aliasing effects incurred by the inverse Fourier transform operation and to better correspond with the plots of (Keefe, 1990). By trial and error, a lowpass filter with a cutoff frequency around 4 kHz was found to produce the best match to Keefe's results. The digital waveguide result was obtained at a sampling rate of 44.1 kHz and then lowpass filtered to a 10 kHz bandwidth, corresponding to that of (Keefe, 1990). Further lowpass filtering is inherent from the first-order Lagrangian, delay-line length interpolation technique used in this model (Valimaki, 1995). Because such filtering is applied at different locations along the "bore," a cumulative effect is difficult to accurately determine. The first tonehole reflection is affected by only two interpolation filters, while the second tonehole reflection is affected by four of these filtering operations. This effect is most responsible for the minor discrepancies apparent in the plots. 2 Three-Port Tonehole Model A tonehole junction may also be represented in the digital waveguide context by a lossless three-port junction. The three-port junction models sound wave interaction at the intersection of the air column and tonehole, as determined by conservation of volume flow and continuity of pressure. Wave propagation within the tonehole itself can subsequently p+(t) -= [l+ro]p+(t)+ropb(t) -2ropth(t) (7b) pt(t) = [1+ ro]p+(t) + [1+ ro]P (t) - [1+ 2ro]pth(t), (7c) where -Yoth -Zo rYoth + 2Yo Zo + 2Zoth (8) A one-multiply form of the three-port scattering equations is given by Pa(t) = Pb(t) + (9a) Pb(t) = pt(t)+w (9b) p (t) = p (t) + p-(t) -pth(t) + w, (9c) where S= ro [P+(t) + p(t) - 2p-t(t)] (10) An implementation of these equations in shown in Fig. 5. Pth-(n) Pth+(n) p,-(n) b-(n) Fig. 5: Tonehole three-port scattering junction implementation in one-multiply form. To complete the digital waveguide three-port tonehole implementation, it is necessary to determine an appropriate model for the tonehole section itself, and then attach this model to the junction. It is

Page  00000004 possible to implement the tonehole structure as a short, fractional delay, digital waveguide and apply an appropriate reflectance at its end. Depending on the tonehole geometry, the reflectance at the end of an open tonehole may be determined from either a flanged or unflanged (Levine and Schwinger, 1948) pipe approximation. The far end of a closed tonehole is appropriately modeled by an infinite impedance (or a pressure reflection without inversion). Given typical tonehole heights, however, a lumped reflectance model of the tonehole, which accounts for both the propagation delay and end reflection is more appropriate and easily implemented with a single low-order digital filter. In this sense, incoming tonehole pressure p4h(t) is calculated from the outgoing tonehole pressure pth s) andu the lumped tonehle diving poI 11 int reflectance, while the corresponding pressure radiated from the open tonehole is given by convolumittance. Figure 6 plots the reflection function obtained for the six-hole flute bore implemented using digital waveguide three-port tonehole junctions. The lumped open-hole reflectance incorporates an unflanged characteristic, while the closed-hole reflectance which best matches the Keefe (1990) data includes no propagation delay within the side branch. Alternatively, the lumped tonehole reflectance filters can be designed from the shunt impedance parameters of Eqs. (3a) and (3b), thus taking advantage of the data of Keefe (1981). The digital waveguide three-port tonehole junction implementation presented here corresponds to the two-port model when series impedance terms are neglected. In general, the series impedance terms are much less critical to the model performance than the shunt impedance, which is demonstrated by the similarity of the results for both implementations. Further, the series terms have more influence on closed-hole results than those for open holes (K~eefe, 1981). Tim (ns Fi. 6 elcinfnto frnt trefne hoe lsd trefne olsoe)o asml lt guide threfe-port junction to neole implhementaion. 3 Conclusions Current theoretical models of woodwind finger holes can be accurately implemented in the digital waveguide domain. The two-port tonehole waveguide implementation requires four second-order filtering operations per tonehole (details regarding a one-filter form are to be published in the proceedings of the 1997 IEEE Workshop on Applications of Signal Processing to Audio and Acoustics). The three-port implementation requires one multiply and one filtering operation. The results for both implementations are very similar, despite the fact that the three-port model neglects the series impedance terms. A more complete and detailed analysis of this topic, unconstrained by page number limitations, can be found in (Scavone, 1997). References Keefe, D. H. (1981). Woodwind Tome-hole Acoustics and the Spectrum Transformation Function. Ph.D. thesis, Case Western Reserve University. Keefe, D. H. (1990). Woodwind air column models. J. Acoust. Soc. Am., 88(1):35-51. K~elly, Jr., J. L. and Lochbaum, C. C. (1962). Speech synthesis. In Proc. Fourth Jut. Congress on Acoustics, pp. 1-4, Copenhagen, Denmark. Paper G42. Levine, H. and Schwinger, J. (1948). On the radiation of sound from an unflanged circular pipe. Phys. Rev., 73(4):383-406. Scavone, G. P. (1997). An. Acoustic Anaylsis of Single-Reed Woodwlinrd Instrumenrts with an Emphasis on Design and Performance Issues and Digital Waveguide Modeling Technriques. Ph.D. thesis, Music Dept., Stanford University. Available as CCRMA Technical Report No. STAN-M-100 or from ftp://ccrmaftp.stanford.edu/pub/Publications/Theses/ay ScavoneThesis/. Smith, J. 0. (1983). Techniques for Digital Filter Design and System Identification with Application to the Vriolin. Ph.D. thesis, Elec. Eng. Dept., Stanford University. Vhlimhki, V. (1995). Discrete-Time M~odeling of Acoustic Tubes Using Fractional Delay Filters. Ph.D. thesis, Helsinki University of Technology, Faculty of Electrical Engineering, Laboratory of Acoustic and Audio Signal Processing, Espoo, Finland, Report no. 37.