Sign in

Not a member? | Forgot your Password?

Search Online Books

Search tips

Free Online Books

Free PDF Downloads

A Quadrature Signals Tutorial: Complex, But Not Complicated

Understanding the 'Phasing Method' of Single Sideband Demodulation

Complex Digital Signal Processing in Telecommunications

Introduction to Sound Processing

C++ Tutorial

Introduction of C Programming for DSP Applications

Fixed-Point Arithmetic: An Introduction

Cascaded Integrator-Comb (CIC) Filter Introduction


FFT Spectral Analysis Software

See Also

Embedded SystemsFPGA

Chapter Contents:

Search Physical Audio Signal Processing


Book Index | Global Index

Would you like to be notified by email when Julius Orion Smith III publishes a new entry into his blog?


Digital Waveguide Model

Figure C.3: Digital simulation of the ideal, lossless waveguide with observation points at $ x=0$ and $ x=3X=3cT$. (The symbol ``$ z^{-1}$'' denotes a one-sample delay.)

In this section, we interpret the sampled d'Alembert traveling-wave solution of the ideal wave equation as a digital filtering framework. This is an example of what are generally known as digital waveguide models [430,431,433,437,442].

The term $ y_r\left[(n-m)T\right]\isdef y^{+}(n-m)$ in Eq.$ \,$(C.16) can be thought of as the output of an $ m$-sample delay line whose input is $ y^{+}(n)$. In general, subtracting a positive number $ m$ from a time argument $ n$ corresponds to delaying the waveform by $ m$ samples. Since $ y^{+}$ is the right-going component, we draw its delay line with input $ y^{+}(n)$ on the left and its output $ y^{+}(n-m)$ on the right. This can be seen as the upper ``rail'' in Fig.C.3

Similarly, the term $ y_l\left[(n+m)T\right]\isdeftext y^{-}(n+m)$ can be thought of as the input to an $ m$-sample delay line whose output is $ y^{-}(n)$. Adding $ m$ to the time argument $ n$ produces an $ m$-sample waveform advance. Since $ y^{-}$ is the left-going component, it makes sense to draw the delay line with its input $ y^{-}(n+m)$ on the right and its output $ y^{-}(n)$ on the left. This can be seen as the lower ``rail'' in Fig.C.3.

Note that the position along the string, $ x_m = mX= m cT$ meters, is laid out from left to right in the diagram, giving a physical interpretation to the horizontal direction in the diagram. Finally, the left- and right-going traveling waves must be summed to produce a physical output according to the formula

$\displaystyle y(t_n,x_m) = y^{+}(n-m) + y^{-}(n+m). \protect$ (C.17)

We may compute the physical string displacement at any spatial sampling point $ x_m$ by simply adding the upper and lower rails together at position $ m$ along the delay-line pair. In Fig.C.3, ``transverse displacement outputs'' have been arbitrarily placed at $ x=0$ and $ x=3X$. The diagram is similar to that of well known ladder and lattice digital filter structures (§C.9.4,[297]), except for the delays along the upper rail, the absence of scattering junctions, and the direct physical interpretation. (A scattering junction implements partial reflection and partial transmission in the waveguide.) We could proceed to ladder and lattice filters by (1) introducing a perfectly reflecting (rigid or free) termination at the far right, and (2) commuting the delays rightward from the upper rail down to the lower rail [432,434]. The absence of scattering junctions is due to the fact that the string has a uniform wave impedance. In acoustic tube simulations, such as for voice [87,297] or wind instruments [195], lossless scattering junctions are used at changes in cross-sectional tube area and lossy scattering junctions are used to implement tone holes. In waveguide bowed-string synthesis (discussed in a later section), the bow itself creates an active, time-varying, and nonlinear scattering junction on the string at the bowing point.

Any ideal, one-dimensional waveguide can be simulated in this way. It is important to note that the simulation is exact at the sampling instants, to within the numerical precision of the samples themselves. To avoid aliasing associated with sampling, we require all waveshapes traveling along the string to be initially bandlimited to less than half the sampling frequency. In other words, the highest frequencies present in the signals $ y_r(t)$ and $ y_l(t)$ may not exceed half the temporal sampling frequency $ f_s
\isdeftext 1/T$; equivalently, the highest spatial frequencies in the shapes $ y_r(x/c)$ and $ y_l(x/c)$ may not exceed half the spatial sampling frequency $ \nu_s \isdeftext 1/X$.

A C program implementing a plucked/struck string model in the form of Fig.C.3 is available at

Previous: Sampled Traveling Waves
Next: Digital Waveguide Interpolation

Order a Hardcopy of Physical Audio Signal Processing

About the Author: Julius Orion Smith III
Julius Smith's background is in electrical engineering (BS Rice 1975, PhD Stanford 1983). He is presently Professor of Music and Associate Professor (by courtesy) of Electrical Engineering at Stanford's Center for Computer Research in Music and Acoustics (CCRMA), teaching courses and pursuing research related to signal processing applied to music and audio systems. See for details.


No comments yet for this page

Add a Comment
You need to login before you can post a comment (best way to prevent spam). ( Not a member? )