Free Books
   Free Books    Physical Audio Signal Processing  


Excitation Examples


Localized Displacement Excitations

Whenever two adjacent components of $ \underline{x}_K$ are initialized with equal amplitude, only a single $ W$-variable will be affected. For example, the initial conditions (for time $ n+1$)

\begin{eqnarray*} y_{n,m-1} &=& 1\\ y_{n-1,m} &=& 1 \end{eqnarray*}

will initialize only $ y^{-}_{n,m-1}$, a solitary left-going pulse of amplitude 1 at time $ n=0$, as can be seen from Eq.$ \,$(E.11) by adding the leftmost columns explicitly written for $ \mathbf{T}^{-1}$. Similarly, the initialization

\begin{eqnarray*} y_{n-1,m-2} &=& 1\\ y_{n,m-1} &=& 1 \end{eqnarray*}

gives rise to an isolated right-going pulse $ y^{+}_{n,m-1}$, corresponding to the leftmost column of $ \mathbf{T}^{-1}$ plus the first column on the left not explicitly written in Eq.$ \,$(E.11). The superposition of these two examples corresponds to a physical impulsive excitation at time 0 and position $ m-1$:

$\displaystyle y_{n-1,m-2}$ $\displaystyle =$ $\displaystyle 1$  
$\displaystyle y_{n,m-1}$ $\displaystyle =$ $\displaystyle 2$  
$\displaystyle y_{n-1,m}$ $\displaystyle =$ $\displaystyle 1$ (E.12)

Thus, the impulse starts out with amplitude 2 at time 0 and position $ m-1$, and afterwards, impulses of amplitude 1 propagate away to the left and right along the string.

In summary, we see that to excite a single sample of displacement traveling in a single-direction, we must excite equally a pair of adjacent colums in $ \mathbf{T}^{-1}$. This corresponds to equally weighted excitation of K-variable pairs the form $ (y_{n,m},y_{n-1,m\pm1})$.

Note that these examples involved only one of the two interleaved computational grids. Shifting over an odd number of spatial samples to the left or right would involve the other grid, as would shifting time forward or backward an odd number of samples.

Localized Velocity Excitations

Initial velocity excitations are straightforward in the DW paradigm, but can be less intuitive in the FDTD domain. It is well known that velocity in a displacement-wave DW simulation is determined by the difference of the right- and left-going waves [437]. Specifically, initial velocity waves $ v^{\pm}$ can be computed from from initial displacement waves $ y^\pm$ by spatially differentiating $ y^\pm$ to obtain traveling slope waves $ y'^\pm$, multiplying by minus the tension $ K$ to obtain force waves, and finally dividing by the wave impedance $ R=\sqrt{K\epsilon }$ to obtain velocity waves:

$\displaystyle v^{+}$ $\displaystyle =$ $\displaystyle -cy'^{+}= \frac{f^{{+}}}{R}$  
$\displaystyle v^{-}$ $\displaystyle =$ $\displaystyle \;cy'^{-}= -\frac{f^{{-}}}{R}, \protect$ (E.13)

where $ c=\sqrt{K/\epsilon }$ denotes sound speed. The initial string velocity at each point is then $ v(nT,mX)=v^{+}(n-m)+v^{-}(n+m)$. (A more direct derivation can be based on differentiating Eq.$ \,$(E.4) with respect to $ x$ and solving for velocity traveling-wave components, considering left- and right-going cases separately at first, and arguing the general case by superposition.)

We can see from Eq.$ \,$(E.11) that such asymmetry can be caused by unequal weighting of $ y_{n,m}$ and $ y_{n,m\pm1}$. For example, the initialization

\begin{eqnarray*} y_{n-1,m+1} &=& +1\\ y_{n-1,m} &=& -1 \end{eqnarray*}

corresponds to an impulse velocity excitation at position $ m+1/2$. In this case, both interleaved grids are excited.


More General Velocity Excitations

From Eq.$ \,$(E.11), it is clear that initializing any single K variable $ y_{n,m}$ corresponds to the initialization of an infinite number of W variables $ y^{+}_{n,m}$ and $ y^{-}_{n,m}$. That is, a single K variable $ y_{n,m}$ corresponds to only a single column of $ \mathbf{T}^{-1}$ for only one of the interleaved grids. For example, referring to Eq.$ \,$(E.11), initializing the K variable $ y_{n-1,m}$ to -1 at time $ n$ (with all other $ y_{n,m}$ intialized to 0) corresponds to the W-variable initialization

\begin{eqnarray*} y^{+}_{n,m-(2\mu+1)}&=&+1, \quad \mu =0,1,2,\cdots\\ y^{-}_{n,m-(2\mu+1)}&=&-1, \quad \mu =0,1,2,\cdots \end{eqnarray*}

with all other W variables being initialized to zero. In view of earlier remarks, this corresponds to an impulsive velocity excitation on only one of the two subgrids. A schematic depiction from $ \mu = m-5$ to $ m+5$ of the W variables at time $ n$ is as follows:

\begin{displaymath}\begin{array}{crrrrr\vert rrrrrrc} \cdots & 1 & 0 & 1 & 0 & 1... ... \cdots\\ & & & & & & m & & & & & \mu & \rightarrow \end{array}\end{displaymath} (E.14)

Below the solid line is the sum of the left- and right-going traveling-wave components, i.e., the corresponding K variables at time $ n$. The vertical lines divide positions $ \mu=m-1$ and $ \mu=m$. The initial displacement is zero everywhere at time $ n$, consistent with an initial velocity excitation. At times $ \nu=n+1,n+2,n+3,n+4$, the configuration evolves as follows:

\begin{displaymath}\begin{array}{crrrrr\vert rrrrrrc} \cdots & 0 & 1 & 0 & 1 & 0... ... 0 & 0 & 0 & 0 & 0 & 1 & 0 & 0 & 0 & 0 & 0 & \cdots \end{array}\end{displaymath} (E.15)

\begin{displaymath}\begin{array}{crrrrr\vert rrrrrrc} \cdots & 1 & 0 & 1 & 0 & 1... ... 0 & 0 & 0 & 0 & 1 & 0 & 1 & 0 & 0 & 0 & 0 & \cdots \end{array}\end{displaymath} (E.16)

\begin{displaymath}\begin{array}{crrrrr\vert rrrrrrc} \cdots & 0 & 1 & 0 & 1 & 0... ... 0 & 0 & 0 & 1 & 0 & 1 & 0 & 1 & 0 & 0 & 0 & \cdots \end{array}\end{displaymath} (E.17)

\begin{displaymath}\begin{array}{crrrrr\vert rrrrrrc} \cdots & 1 & 0 & 1 & 0 & 1... ... 0 & 0 & 1 & 0 & 1 & 0 & 1 & 0 & 1 & 0 & 0 & \cdots \end{array}\end{displaymath} (E.18)

The sequence $ [\dots,1,0,1,0,1,\dots]$ consists of a dc (zero-frequency) component with amplitude $ 1/2$, plus a sampled sinusoid of amplitude $ 1/2$ oscillating at half the sampling rate $ f_s=1/T$. The dc component is physically correct for an initial velocity point-excitation (a spreading square pulse on the string is expected). However, the component at $ f_s/2$ is usually regarded as an artifact of the finite difference scheme. From the DW interpretation of the FDTD scheme, which is an exact, bandlimited physical interpretation, we see that physically the component at $ f_s/2$ comes about from initializing velocity on only one of the two interleaved subgrids of the FDTD scheme. Any asymmetry in the excitation of the two interleaved grids will result in excitation of this frequency component.

Due to the independent interleaved subgrids in the FDTD algorithm, it is nearly always non-physical to excite only one of them, as the above example makes clear. It is analogous to illuminating only every other pixel in a digital image. However, joint excitation of both grids may be accomplished either by exciting adjacent spatial samples at the same time, or the same spatial sample at successive times instants.

In addition to the W components being non-local, they can demand a larger dynamic range than the K variables. For example, if the entire semi-infinite string for $ m<0$ is initialized with velocity $ 2X/T$, the initial displacement traveling-wave components look as follows:

\begin{displaymath}\begin{array}{crrrrrr\vert rrrrrc} \cdots & 6 & 5 & 4 & 3 & 2... ... 0 & 0 & 0 & 0 & 0 & 0 & 0 & 0 & 0 & 0 & 0 & \cdots \end{array}\end{displaymath} (E.19)

and the variables evolve forward in time as follows:

\begin{displaymath}\begin{array}{crrrrrr\vert rrrrrc} \cdots & 7 & 6 & 5 & 4 & 3... ... 2 & 2 & 2 & 2 & 2 & 2 & 1 & 0 & 0 & 0 & 0 & \cdots \end{array}\end{displaymath} (E.20)

\begin{displaymath}\begin{array}{crrrrrr\vert rrrrrc} \cdots & 8 & 7 & 6 & 5 & 4... ... 4 & 4 & 4 & 4 & 4 & 3 & 2 & 1 & 0 & 0 & 0 & \cdots \end{array}\end{displaymath} (E.21)

\begin{displaymath}\begin{array}{crrrrrr\vert rrrrrc} \cdots & 9 & 8 & 7 & 6 & 5... ... 6 & 6 & 6 & 6 & 5 & 4 & 3 & 2 & 1 & 0 & 0 & \cdots \end{array}\end{displaymath} (E.22)

Thus, the left semi-infinite string moves upward at a constant velocity of 2, while a ramp spreads out to the left and right of position $ \mu=m$ at speed $ c$, as expected physically. By Eq.$ \,$(E.9), the corresponding initial FDTD state for this case is

\begin{eqnarray*} y_{n,\mu} &=&0, \quad \mu\in{\bf Z}\\ y_{n-1,m-1} &=& -1,\\ y_{n-1,\mu} &=& -2, \quad \mu<m-1, \end{eqnarray*}

where $ {\bf Z}$ denotes the set of all integers. While the FDTD excitation is also not local, of course, it is bounded for all $ \mu $.

Since the traveling-wave components of initial velocity excitations are generally non-local in a displacement-based simulation, as illustrated in the preceding examples, it is often preferable to use velocity waves (or force waves) in the first place [447].

Another reason to prefer force or velocity waves is that displacement inputs are inherently impulsive. To see why this is so, consider that any physically correct driving input must effectively exert some finite force on the string, and this force is free to change arbitrarily over time. The ``equivalent circuit'' of the infinitely long string at the driving point is a ``dashpot'' having real, positive resistance $ 2R=2\sqrt{K\epsilon }$. The applied force $ f(t)$ can be divided by $ 2R$ to obtain the velocity $ v(t)$ of the string driving point, and this velocity is free to vary arbitrarily over time, proportional to the applied force. However, this velocity must be time-integrated to obtain a displacement $ y(t)$. Therefore, there can be no instantaneous displacement response to a finite driving force. In other words, any instantaneous effect of an input driving signal on an output displacement sample is non-physical except in the case of a massless system. Infinite force is required to move the string instantaneously. In sampled displacement simulations, we must interpret displacement changes as resulting from time-integration over a sampling period. As the sampling rate increases, any physically meaningful displacement driving signal must converge to zero.

Additive Inputs

Instead of initial conditions, ongoing input signals can be defined analogously. For example, feeding an input signal $ u_n$ into the FDTD via

$\displaystyle y_{n,m-1}$ $\displaystyle =$ $\displaystyle y_{n,m-1} + u_{n-1}$  
$\displaystyle y_{n,m}$ $\displaystyle =$ $\displaystyle y_{n,m} + 2u_n$  
$\displaystyle y_{n,m+1}$ $\displaystyle =$ $\displaystyle y_{n,m+1} + u_{n-1} \protect$ (E.23)

corresponds to physically driving a single sample of string displacement at position $ m$. This is the spatially distributed alternative to the temporally distributed solution of feeding an input to a single displacement sample via the filter $ H(z)=1-z^{-2}$ as discussed in [223].

Interpretation of the Time-Domain KW Converter

As shown above, driving a single displacement sample $ y_{n,m}$ in the FDTD corresponds to driving a velocity input at position $ m$ on two alternating subgrids over time. Therefore, the filter $ H(z)=1-z^{-2}$ acts as the filter $ H(z)=1-z^{-1}$ on either subgrid alone--a first-order difference. Since displacement is being simulated, velocity inputs must be numerically integrated. The first-order difference can be seen as canceling this integration, thereby converting a velocity input to a displacement input, as in Eq.$ \,$(E.23).

Next Section:
State Space Formulation
Previous Section:
State Transformations