The helical coordinate

The spectral factorization concept

Interesting questions arise when we are given a spectrum and find ourselves asking how to find a filter that has that spectrum. Is the answer unique? We'll see not. Is there always an answer that is causal? Almost always, yes. Is there always an answer that is causal with a causal inverse (CwCI)? Almost always, yes.

Let us have an example. Consider a filter like the familiar time derivative $(1,-1)$ except let us downweight the $-1$ a tiny bit, say $(1,-\rho)$ where $0«\rho<1$ . Now the filter $(1,-\rho)$ has a spectrum $(1-\rho Z)(1-\rho/Z)$ with autocorrelation coefficients $(-\rho, 1+\rho^2,-\rho)$ that look a lot like a second derivative, but it is a tiny bit bigger in the middle. Two different waveforms, $(1,-\rho)$ and its time reverse both have the same autocorrelation. Spectral factorization could give us both $(1,-\rho)$ and $(\rho,-1)$ but we always want the one that is CwCI. The bad one is weaker on its first pulse. Its inverse is not causal. Below are two expressions for the filter inverse to $(\rho,-1)$ , the first divergent (filter coefficients at infinite lag are infinitely strong), the second convergent but noncausal.

$$\frac{1}{ \rho -Z} = \frac{ 1}{\rho}\ ( 1 +Z/\rho +Z^2/\rho^2+ \cdots)$$ (14)

$$\frac{1}{ \rho -Z} = \frac{-1}{ Z}\ ( 1 + \rho/Z + \rho^2/Z^2 + \cdots)$$ (15)

(Please multiply each equation by $\rho -Z$ and see it reduce to $1=1$ ).

So we start with a power spectrum and we should find a CwCI filter with that energy spectrum. If you input to the filter an infinite sequence of random numbers (white noise) you should output something with the original power spectrum.

We easily inverse Fourier transform the square root of the power spectrum getting a symmetrical time function, but we need a function that vanishes before $\tau=0$ . On the other hand, if we already had a causal filter with the correct spectrum we could manufacture many others. To do so all we need is a family of delay operators to convolve with. A pure delay filter does not change the spectrum of anything. Same for frequency-dependent delay operators. Here is an example of a frequency-dependent delay operator: First convolve with (1,2) and then deconvolve with (2,1). Both these have the same amplitude spectrum so their ratio has a unit amplitude (and nontrivial phase). If you multiply $(1+2Z)/(2+Z)$ by its Fourier conjugate (replace $Z$ by $1/Z$ ) the resulting spectrum is 1 for all $\omega$ .

Anything whose nature is delay is death to CwCI. The CwCI has its energy as close as possible to $\tau=0$ . More formally, my first book, FGDP, proves that the CwCI filter has for all time $\tau$ more energy between $t=0$ and $t=\tau$ than any other filter with the same spectrum.

Spectra can be factorized by an amazingly wide variety of techniques, each of which gives you a different insight into this strange beast. They can be factorized by factoring polynomials, by inserting power series into other power series, by solving least squares problems, by taking logarithms and exponentials in the Fourier domain. I've coded most of them and still find them all somewhat mysterious.

Theorems in Fourier analysis can be interpreted physically in two different ways, one as given, the other with time and frequency reversed. For example, convolution in one domain amounts to multiplication in the other. If we were to express the CwCI concept with reversed domains, instead of saying the ``energy comes as quick as possible after $\tau=0$ '' we would say ``the frequency function is as close to $\omega=0$ as possible.'' In other words, it is minimally wiggly with time. Most applications of spectral factorization begin with a spectrum, a real, positive function of frequency. I once achieved minor fame by starting with a real, positive function of space, a total magnetic field $\sqrt{H_x^2 +H_z^2}$ measured along the $x$ -axis and I reconstructed the magnetic field components $H_x$ and $H_z$ that were minimally wiggly in space (FGDP p.61).

The helical coordinate