Hilbert transform and Fourier series

A few days ago I wrote about the Hilbert transform and gave as an example that the Hilbert transform of sine is cosine. We'll bootstrap that example to find the Hilbert transform of any periodic function from its Fourier series.

The Hilbert transform of a function f(t) is a function f_H(x) defined by

where the integral is interpreted in the sense of the Cauchy principal value, the limit as the singularity is approach symmetrically from both sides.

The Hilbert transform shifts and scales conveniently. Shifting a function by any amount h shifts its transform by the same amount. And scaling a function by any amount k > 0 scales its transform the same way. That is, we have the following transform pairs.

Now since the Hilbert transform of sin(t) is cos(x), the Hilbert transform of sin(t + /2) must be cos(x + /2). But sin(t + /2) is cos(t), and cos(x + /2) is sin(x), so the Hilbert transform of cos(t) is -sin(x). In this case, the Hilbert transform has the same pattern as differentiation.

Now if > 0 the scaling rule tells us the Hilbert transform of sin(t) is cos(x) and the Hilbert transform of cos(x) is -sin(x). Here the analogy with differentiation breaks down because differentiation would bring out a factor of from the chain rule [1].

Putting these facts together, if we have a function f written in terms of a Fourier series

then its Hilbert transform is

In other words, we replace the bs with as and the as with -bs. [2]

Notice that there's no b₀ term above. In signal processing terminology, there's no DC offset. In general a Fourier series has a constant term, and the Hilbert transform of a constant is 0. So again like differentiation, constants go away.

If there is no DC offset, then applying the Hilbert transform to f twice gives -f. If there is a DC offset, applying the Hilbert transform to f twice gives -f with the DC offset removed.

Opposite sign convention

Unfortunately there are two definitions of the Hilbert transform in common use: the one at the top of this post and its negative. What changes if we use the other convention?

We noted above that applying the Hilbert transform to f twice gives -f. This means that the inverse transform is the negative transform [3]. In symbols, if H is the Hilbert transform operator, then -H^2f = -Hf and so H^-1 = -H. So the disagreement over whether to include a negative sign in the definition of the Hilbert transform amounts to a disagreement over which to call the forward transform and which to call the inverse transform.

The shifting and scaling rules apply to both definitions of the transform. But with the opposite sign convention, the Hilbert transform of sine is negative cosine and the Hilbert transform of cosine is sine. So our bottom line becomes replace the as with bs and the bs with -as" in the Fourier series.

Footnotes

[1] Incidentally, the Hilbert transform commutes with differentiation. That is, the transform of the derivative is the derivative of the transform.

[2] This is an example of parallel replacement. We replace a_n, with -b_n and b_n, with a_n at the same time.

[3] For signals with no DC offset. Otherwise the Hilbert transform is not invertible.

The post Hilbert transform and Fourier series first appeared on John D. Cook.