Two quantities have been proposed as criteria for measuring the ``naturalness'' of synthetic reverberation: the time density and the frequency density . A good reverberator should provide high values of both densities, thus giving smooth, dense time and frequency responses.
The frequency density is defined as the average number of resonances per
Hertz. A general expression can be derived from the order of the
system (5), assuming that all the poles are distinct and no
In the prototype case where the delay lines all have the same length , we have
The time density is defined as the number of nonzero samples per second in the impulse response. In actual rooms, is an increasing function of time. In order to obtain dense reverberation after the early reflections (e.g., after 80 msec), it helps to use different delay lengths.
The actual positions of frequency peaks depend on the feedback matrix and the delay lengths. If the delay lengths are fixed, we can vary some time-frequency properties of the structure simply by varying the distribution of eigenvalues of the feedback matrix. The total length of the delay lines should be chosen in such a way that the frequency density, as determined by (32), is high enough. Then the matrix eigenvalues can be adjusted to avoid resonant peak clustering or other undesirable mode distributions.
It is interesting to discuss the effect of eigenvalues in the prototype case of equal delays. A uniform distribution of eigenvalues along the unit circle is optimum for the frequency response in the sense that it minimizes the maximum distance between peaks. However, it produces a highly repetitive time response. Conversely, clustering the eigenvalues around a point on the unit circle can be good for maximizing the length of time patterns, but the clustering of frequency peaks produces a poor reverberator amplitude response vs. frequency. We see from these considerations that there is a time-frequency tradeoff. This tradeoff can be addressed using circulant matrices.
A couple of examples of different eigenvalue distributions are given in
Fig. Fig. 3. The matrix used in Fig. Fig. 3b is
simply obtained by a right circular shift of the rows of the matrix
which is given by the junction of equal-impedance
waveguides and, as already stated, has eigenvalues only at and .
We can express as the product
is the right-shift matrix
As a side comment, we notice that is the scattering matrix of the circulator, a circuit device which can be used to obtain the multiplication of one-port scattering parameters .
The shape of the frequency response depends also on the zeros which were discussed in section IV. In particular, Theorem 2 provides a way of setting the zeros exactly over the poles in the prototype equal-delay case. We anticipated in section IV that the way to choose the vectors and indicated in Theorem 2 can be useful for getting a flat amplitude response at low frequencies when the delay lengths are slightly varied from the prototype case. Fig. Fig. 4 depicts the time and frequency responses for the CFDN using the same feedback matrix as in Fig. Fig. 3b, having , , and delay lengths . As we can see from Fig. Fig. 4, we are able to get a nearly flat amplitude response at low frequencies without losing the reverberating character of the time response. We believe that this is a good alternative to allpass filters which tend to have degenerate impulse responses when the poles approach the unit circle.