Sign in

username:

password:



Not a member?

Search Online Books



Search tips

Free Online Books

Ads

Chapters

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?

  

Digitizing Elementary Reflectances via the Bilinear Transform

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

$\displaystyle s = c \frac{1-z^{-1}}{1+z^{-1}}$ (Q.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$ (Q.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.$ \,$(Q.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