A method was previously reported [8] for converting the continuous-time woodwind tonehole model of Keefe [2] to a discrete-time two-port scattering junction model for implementation in the digital waveguide (DW) domain. The results using this technique are shown in Figure 1 and compared with reproduced results using the technique of Keefe [3] for a simple flute air column with six toneholes. Discrepancies between the two methods are mainly evident in early closed hole reflections. Keefe's results were calculated for a frequency range of 10 kHz and subsequently smoothed in the time-domain with a hamming window [5]. By trial and error, a lowpass filter with a cutoff frequency around 4 kHz was found to best reproduce Keefe's results. The DW results were obtained at a sampling rate of 44.1 kHz and then lowpass filtered to a 10 kHz bandwidth to correspond with the calculations of [3]. Further lowpass filtering is inherent from the Lagrangian, delay-line length interpolation technique used in this model [6]. Because such filtering is applied at different locations along the air column and is dependent on the particular fractional delay length modeled, a cumulative effect is difficult to accurately determine. As diagrammed in Figure 2, that portion of the signal reflected at the first tonehole is affected by only two interpolation filters, that at the second tonehole reflection is affected by four filtering operations, etc. Thus, *early* reflections in the DW model results are less lowpass filtered than the results of [3]. It should be noted that each fractional delay interpolation filter in this implementation can be combined with a lossy propagation filter, which models lumped thermoviscous losses along its corresponding segment of the air column and which is also given by a lowpass frequency response. In this way, the inaccuracies inherent in low-order delay length interpolation filters can often be minimized. Alternately, higher-order interpolation filters can be used which introduce minimal frequency magnitude distortion.

For the purpose of real-time modeling, the two-port implementation has a particular disadvantage: the two lumped characterizations of the tonehole as either closed or open cannot be efficiently unified into a single tonehole model. While it is possible to develop a cross-fading/interpolation scheme to simulate ``half-holing'', this would require that two simultaneous models be run to simulate just one tonehole. It is preferable to have one model with adjustable parameters to simulate the various states of the tonehole, from closed to open and all states in between.

To this end, it is best to consider a distributed model of the tonehole, such that ``fixed'' portions of the tonehole structure are separated from the ``variable'' component. The junction of the tonehole branch with the main air column of the instrument can be modeled in the DW domain using a three-port scattering junction, as described in [8]. This method inherently models only the shunt impedance term of the Keefe tonehole characterization, however, the negative length correction terms implied by the series impedances can be approximated by adjusting the delay line lengths on either side of the three-port scattering junction. The other ``fixed'' portion of the tonehole is the short branch segment itself, which is modeled in the DW domain by appropriately sized delay lines. This leaves only the characterization of the open/closed tonehole end. A simple inertance model of the open hole end offers the most computationally efficient solution. The impedance of the open end is then given by

(1) |

where is the characteristic impedance of the tonehole branch waveguide and is the speed of sound. An appropriate discrete-time filter implementation for can be obtained using the conformal bilinear transform from the -plane to the -plane [7, pp. 415-430], with the result

(3) |

(4) |

Download tonehole.pdf

Copyright ©

Center for Computer Research in Music and Acoustics (CCRMA), Stanford University

[Automatic-links disclaimer]