next up previous print clean
Next: SMOOTHING WITH BOX AND Up: LEAKY INTEGRATION Previous: Plots

Two poles

Integration twice is an operation with two poles. Specifically,
\begin{displaymath}
{1\over (1-Z)^2}
\ =\ 
(1+Z+Z^2+Z^3+ \cdots)
(1+Z+Z^2+Z^3+ \cdots)
\ =\ 
1 + 2Z + 3Z^2 + 4Z^3 + 5Z^4 + \cdots\end{displaymath} (9)
Notice that the signal is $(1,2,3,\cdots)$,which is a discrete representation of the function $f(t)=t\,{\rm step} (t)$.Figure 2 shows the result when the two integrations are leaky integrations.

 
leak2
leak2
Figure 2
A cascade of two leaky integrators.


view

We see the signal begin as t but then drop off under some weight that looks exponential. A second time-derivative filter $(-i\omega)^2$has an amplitude spectrum $\vert\omega^2\vert$.Likewise, a second integration has an amplitude spectrum $\vert 1/\omega^2\vert$,which is about what we see in Figure 2, except that at $\omega =0$leaky integration has rounded off the $\infty$.

Instead of allowing two poles to sit on top of each other (which would look like just one pole), I moved the pole slightly off $\Re \omega = 0$so that $\Re\omega\gt$.As in Figure [*], another pole is included (but not shown) at negative frequency. This extra pole is required to keep the signal real. Of course the two poles are very close to each other. The reason I chose to split them this way is to prepare you for filters where the poles are far apart.

EXERCISES:

  1. Show that multiplication by (1 - Z) in discretized time is analogous to time differentiation in continuous time. Show that dividing by (1 - Z) is analogous to integration. What are the limits on the integral?
  2. A simple feedback operation is $y_t = (1 - \epsilon)y_{t - 1} + x_t$.Give a closed-form expression for the output yt if xt is an impulse. Estimate the decay time $\tau$ of your solution (the time it takes for yt to drop to e-1y0)? For small $\epsilon$,say = 0.1, .001, or 0.0001, what is $\tau$?
  3. Find an analytic expression for the plot on the right side of Figure 1 as a function of $\omega$.Show that it is like $1/\vert\omega\vert$.
  4. In continuous time, the signal analogous to that in Figure 2 is te-t. What is the analogous frequency function?

next up previous print clean
Next: SMOOTHING WITH BOX AND Up: LEAKY INTEGRATION Previous: Plots
Stanford Exploration Project
10/21/1998