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 [434,435,437,441,445].
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
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 http://ccrma.stanford.edu/~jos/pmudw/.