Decay Time is Q Periods Next  |  Prev  |  Up  |  Top  |  Index  |  JOS Index  |  JOS Pubs  |  JOS Home  |  Search


Decay Time is Q Periods

Another well known rule of thumb is that the $ Q$ of a resonator is the number of ``periods'' under the exponential decay of its impulse response. More precisely, we will show that, for $ Q\gg 1/2$, the impulse response decays by the factor $ e^{-\pi}$ in $ Q$ cycles, which is about 96 percent decay, or -27 dB.

The impulse response corresponding to Eq. (C.8) is found by inverting the Laplace transform of the transfer function $ H(s)$. Since it is only second order, the solution can be found in many tables of Laplace transforms. Alternatively, we can break it up into a sum of first-order terms which are invertible by inspection (possibly after rederiving the Laplace transform of an exponential decay, which is very simple). Thus we perform the partial fraction expansion of Eq. (C.8) to obtain

$\displaystyle H(s) = \frac{g_1}{s-p_1} + \frac{g_2}{s-p_2}
$

where $ p_i$ are given by Eq. (C.9), and some algebra gives
$\displaystyle g_1$ $\displaystyle =$ $\displaystyle -g\frac{p_1}{p_2-p_1}$ (C.12)
$\displaystyle g_2$ $\displaystyle =$ $\displaystyle g\frac{p_2}{p_2-p_1}$ (C.13)

as the respective residues of the poles $ p_i$.

The impulse response is thus

$\displaystyle h(t) = g_1 e^{p_1t} + g_2 e^{p_2t}.
$

Assuming a resonator, $ Q>1/2$, we have $ p_2 = \overline{p}_1$, where $ p_1=\sigma_p +j\omega_p = -\alpha + j\omega_d$ (using notation of the preceding section), and the impulse response reduces to

$\displaystyle h(t) = g\,e^{-\alpha t} \left[\cos(\omega_p t) - \frac{\alpha}{\omega_p}
\sin(\omega_p t)\right]
= g\,e^{-\alpha t} \cos(\omega_p t + \phi),
$

where $ \phi\isdef \sin^{-1}(\alpha/\omega_p)$.

We have shown so far that the impulse response $ h(t)$ decays as $ e^{-\alpha t}$ with a sinusoidal radian frequency $ \omega_p=\omega_d$ under the exponential envelope. After Q periods at frequency $ \omega_p$, time has advanced to

$\displaystyle t_Q = Q\frac{2\pi}{\omega_p}
\approx \frac{2\pi Q}{\omega_0},
= \frac{\pi}{\alpha},
$

where we have used the definition Eq. (C.7) $ Q\isdef \omega_0/(2\alpha)$. Thus, after $ Q$ periods, the amplitude envelope has decayed to

$\displaystyle e^{-\alpha t_Q} = e^{-\pi} \approx 0.043\dots
$

which is about 96 percent decay. The only approximation in this derivation was

$\displaystyle \omega_p = \sqrt{\omega_0^2 - \alpha^2} \approx \omega_0
$

which holds whenever $ \alpha\ll\omega_0$, or $ Q\gg 1/2$.


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

[How to cite this work] [Order a printed hardcopy]

``Introduction to Digital Filters with Audio Applications'', by Julius O. Smith III, (August 2006 Edition).
Copyright © 2007-02-02 by Julius O. Smith III
Center for Computer Research in Music and Acoustics (CCRMA),   Stanford University
CCRMA  [Automatic-links disclaimer]