The helical coordinate

INVERSE FILTERS AND OTHER FACTORIZATIONS

Mathematics sometimes seems a mundane subject, like when it does the ``accounting'' for an engineer. Other times, as with the study of causality and spectral factorization, it brings unexpected amazing new concepts into our lives. There are many little-known, fundamental ideas here; a few touched on next.

Start with an example. Consider a mechanical object. We can strain it and watch it stress or we can stress it and watch it strain. We feel knowledge of the present and past stress history is all we need to determine the present value of strain. Likewise, the converse, history of strain should tell us the stress. We could say there is a filter that takes us from stress to strain; likewise, another filter takes us from strain to stress. What we have here is a pair of filters that are mutually inverse under convolution. In the Fourier domain, one is literally the inverse of the other. What is remarkable is that in the time domain, both are causal. They both vanish before zero lag $\tau=0$ .

Not all causal filters have a causal inverse. The best known name for one that does is ``minimum-phase filter.'' Unfortunately, this name is not suggestive of the fundamental property of interest, ``causal with a causal (convolutional) inverse.'' I could call it CCI. An example of a causal filter without a causal inverse is the unit delay operator--with $Z$ -transforms, the operator $Z$ . If you delay something, you cannot get it back without seeing into the future, which you are not allowed to do. Mathematically, $1/Z$ cannot be expressed as a polynomial (actually, a convergent infinite series) in positive powers of $Z$ .

Physics books do not tell us where to expect to find transfer functions that are CCI. I think I know why they do not. Any causal filter has a ``sharp edge'' at zero time lag where it switches from nonresponsiveness to responsiveness. The sharp edge might cause the spectrum to be large at infinite frequency. If so, the inverse filter is small at infinite frequency. Either way, one of the two filters is unmanageable with Fourier transform theory, which (you might have noticed in the mathematical fine print) requires signals (and spectra) to have finite energy. Finite energy means the function must get really small in that immense space on the $t$ -axis and the $\omega$ axis. It is impossible for a function to be small and its inverse be small. These imponderables become manageable in the world of Time Series Analysis (discretized time axis).

The helical coordinate