next up previous print clean
Next: FILTERS IN PARALLEL Up: Spectrum and phase Previous: Phase of a rational

ROBINSON'S ENERGY-DELAY THEOREM

Here we will see that a minimum-phase filter

has less energy delay than any other one-sided filter with the same spectrum.  More precisely, the energy summed from zero to any time t for the minimum-phase wavelet is greater than or equal to that of any other wavelet with the same spectrum.

Here is how I prove Robinson's energy-delay theorem: compare two wavelets, $F_{\rm in}$ and $F_{\rm out}$,that are identical except for one zero, which is outside the unit circle for $F_{\rm out}$and inside for $F_{\rm in}$.We can write this as
\begin{eqnarray}
F_{\rm out}(Z) &= & (b + sZ)\, F(Z) \ F_{\rm in}(Z) &= & (s + bZ)\, F(Z) \end{eqnarray} (38)
(39)
where b is bigger than s, and F is arbitrary but of degree n. Proving the theorem for complex-valued b and s is left as an exercise. Notice that the spectrum of b+sZ is the same as that of s+bZ. Next, tabulate the terms in question.

t $F_{\rm out}$ $F_{\rm in}$ $F^2_{\rm out} - F^2_{\rm in}$ $\sum^t_{k = 0} (F^2_{\rm out} - F^2_{\rm in})$
  bf0 sf0 $(b^2 - s^2)\, {f_0}^2 $ $(b^2 - s^2)\, {f_0}^2 $
1 bf1 + sf0 sf1 + bf0 $(b^2 - s^2)\, ({f_1}^2 - {f_0}^2)$ $(b^2 - s^2)\, {f_1}^2$
$\vdots$ $\vdots$      
k bfk + sfk -1 sfk + bfk -1 $(b^2 - s^2)\, ({f_k}^2 - f_{k - 1}^2)$ $(b^2 - s^2)\, {f_k}^2$
$\vdots$ $\vdots$      
n + 1 sfn bfn (b2 - s2) (-fn2)  

The difference, which is given in the right-hand column, is always positive. An example of the result is shown in Figure 20.

 
robinson
Figure 20
Total energy
[2] percentage versus time. `x' for minimum-phase wavelet. `o' for nonminimum phase.

robinson
view burn build edit restore

Notice that (s+bZ)/(b+sZ) is an all-pass filter. Multiplying by an all-pass filter does not change the amplitude spectrum but instead introduces a zero and a pole. The pole could cancel a preexisting zero, however. To sum up, multiplying by a causal/anticausal all-pass filter can move zeros inside/outside the unit circle. Each time we eliminate a zero inside the unit circle, we cause the energy of the filter to come out earlier. Eventually we run out of zeros inside the unit circle, and the energy comes out as early as possible.

EXERCISES:

  1. Repeat the proof of Robinson's minimum-energy-delay theorem for complex-valued b, s, and fk. (HINT: Does $F_{\rm in} = (s + bZ)\, F$ or $F_{\rm in} = (\overline{s} + \overline{b} Z)F$?)

next up previous print clean
Next: FILTERS IN PARALLEL Up: Spectrum and phase Previous: Phase of a rational
Stanford Exploration Project
10/21/1998