Conditions for Losslessness

Digital Waveguide Networks

When in the proof is diagonal and positive, a physical waveguide interpretation always exists with . A generalized waveguide interpretation exists for all via vector transformers [28, p. 55 sec. 4] in which acts as an ideal transformer (in the classical network theory sense) on the vector of all waveguide variables. If denotes the vector of physical pressures at the junction and denotes the physical volume velocities, then we have that the junction power, defined as is invariant with respect to insertion of a vector transformer (similarity transformation applied to the scattering matrix).

It can be quickly verified^{3} that all scattering matrices arising from the parallel
intersection of physical waveguides possess one eigenvalue equal to
and eigenvalues equal to [28]. In the case of physical
waveguides of equal impedances, the eigenvector associated with the
eigenvalue corresponds to equal incoming waves, while an eigenvector
associated with the eigenvalue corresponds to equal incoming waves on
branches, and a large opposite wave on the remaining branch which
pulls the junction pressure to zero. Adding a vector transformer to the
parallel -branch scattering junction gives scattering matrices of the
form
with
. To obtain more general
eigenvalue signatures
, a combination of series and
parallel junctions must be used. Finally, to reach the most general
complex case, we must admit complex eigenvalues.

We now consider the relation between unitary FDN feedback matrices and waveguide scattering junctions. As can be seen from comparing (27) to (true for any unitary matrix), we see that unitary corresponds to a scattering junction in which the total complex power is given by the ordinary norm of the incoming or outgoing traveling waves. Since the physical power associated with an incoming wave vector is , where in the absence of a vector transformer , we see that unitary corresponds to a scattering junction joining waveguides of equal wave impedance, i.e., . Since Householder reflections comprise only a subset unitary matrices, we see that a unitary FDN matrix corresponds to a transformer-coupled parallel/series waveguide junction in which all branch admittances are the same. In the more general (unnormalized) case in which the branch impedances are different, i.e., , we obtain (using a vector transformer) the larger class of scattering matrices which preserve an elliptic norm as induced by a positive-definite (or Hermitian) generalized junction impedance.

Since, as discussed above, only a subset of all -by- unitary matrices is given by a physical junction of waveguides, the unitary FDN point of view yields lossless systems outside the scope of those suggested by multiport scattering theory. On the other hand, since only normalized waveguide junctions exhibit unitary scattering matrices, the DWN approach gives rise to new classes of FDNs. Moreover, by considering more than one scattering junction, the DWN approach suggests new classes of network topologies following physical analogies. Similarly, FDN matrices can be partitioned to embed several FDN subsystems into larger FDN systems.

Formally, every DWN can be expressed as an FDN by collecting all of its delay lines into a diagonal delay matrix as in (4), and finding the matrix which computes the delay-line inputs from the delay-line outputs. Therefore, every waveguide network yields a feedback matrix for consideration in the FDN framework. Conversely, every real FDN can be expressed as a single-junction waveguide network using an ideal vector transformer at the junction.

Theorem 1 characterizes lossless FDN feedback matrices as those
having eigenvalues on the unit circle, where the definition of losslessness
was given by (27). It remains to be shown that satisfying
(27) implies that the poles of the corresponding
FDN are all on the unit circle. To this end, recall the form of the
state-transition matrix (8), and define the extended
generalized admittance

(28) |

This equation shows that is lossless if and only if is lossless, and its eigenvalues are on the unit circle by Theorem 1. But from (17), we have that the eigenvalues of are the poles of the corresponding FDN. Therefore, the FDN is lossless and its impulse response consists only of non-increasing and non-decaying modes. Desired decay characteristics versus frequency for obtaining a specific artificial reverberator can now be controlled separately by means of attenuation coefficients and filters, as indicated in section II.

Finite-Wordlength Effects

Conditions for Losslessness

Digital Waveguide Networks

``Circulant and Elliptic Feedback Delay Networks for Artificial Reverberation'', by Davide Rocchesso and Julius O. Smith III, preprint of version in

Download PDF version (cfdn.pdf)

Download compressed PostScript version (cfdn.ps.gz)

(Browser settings for best viewing results)

Copyright ©

Center for Computer Research in Music and Acoustics (CCRMA), Stanford University

(automatic links disclaimer)