Next  |  Prev  |  Up  |  Top  |  Index  |  JOS Index  |  JOS Pubs  |  JOS Home  |  Search

Digitizing Elementary Reflectances by Bilinear Transform

Going to discrete time via the bilinear transform means making the substitution

$\displaystyle s = c \frac{1-z^{-1}}{1+z^{-1}}$ (F.11)

where $ c>0$ is an arbitrary real constant, usually taken to be $ c=2/T$ .

Solving for $ z^{-1}$ gives us the inverse bilinear transform:

$\displaystyle z^{-1}= \frac{1-s/c}{1+s/c} \protect$ (F.12)

In this case, we see that setting $ c=1$ further simplifies our universal reflectances in the digital domain:

Note that this choice of $ c$ is also the only one that eliminates delay-free paths in the fundamental elements. This allows them to be used as building blocks for explicit finite difference schemes.

We may still obtain the above results using the more typical value $ c=2/T$ (instead of $ c=1$ ) in the bilinear transform. From Eq.$ \,$ (F.12), it is clear that changing $ c$ amounts to a linear frequency scaling of $ s=j\omega$ . Such a scaling may be compensated by choosing the waveguide (port) impedances to be $ R_L = Lc = 2L/T$ (instead of $ R_L=L$ ) for the inductor, and $ R_C = T/(2C)$ (instead of $ 1/C$ ) for the capacitor.


Next  |  Prev  |  Up  |  Top  |  Index  |  JOS Index  |  JOS Pubs  |  JOS Home  |  Search

[How to cite this work]  [Order a printed hardcopy]  [Comment on this page via email]

``Physical Audio Signal Processing'', by Julius O. Smith III, W3K Publishing, 2010, ISBN 978-0-9745607-2-4.
Copyright © 2014-03-23 by Julius O. Smith III
Center for Computer Research in Music and Acoustics (CCRMA),   Stanford University
CCRMA