In this article we take a look at group velocity and at the extraction of the envelope of a wave packet using the ideas of the Hilbert transform.

**Interference of two waves**

A single sinusoidal wave is infinite in extent and periodic in space and time. When waves interact, the dynamics are more interesting. The simplest case is the superposition of two waves. Assume the two components have equal amplitudes and approximately equal wavenumbers and frequencies:

The two components move with phase speeds and . We write the mean values and and the differences and . Then by elementary trigonometry we have

The second factor here (the carrier wave) represents a wave with wavenumber moving with phase speed , close to the phase speeds of the two components. The first factor is slowly varying in space. It is called the *envelope*, and has wavenumber and frequency . It moves with a speed , called the *group velocity*,

This velocity may be radically different from the phase velocity . More generally, we can write the group velocity as

**Group velocity of Rossby waves**

Let us look at planetary waves in the atmosphere. For nondivergent quasigeostrophic flow on a beta plane the *Rossby wave* phase speed is

where is the mean zonal flow, is the wavenumber and is constant. Because of the minus sign, the phase speed of the waves is *westward* relative to the mean flow. We compute the group velocity:

We have the surprising result that the group velocity is directed towards the east (relative to the mean flow), opposite to the phase velocity. The figure below is a Hovmöller diagram, showing the wave amplitude as a function of longitude (horizontal) and time (downward axis). The blue arrow shows the slow progress of an individual wave maximum. The red arrow shows the much faster propagation of the wave group.

Group velocity is of immense importance in weather forecasting. The large wave-like disturbances in the atmosphere at middle latitudes travel at the phase speed . But the energy travels at the group velocity , which is often *much greater*. Wave minima or troughs are commonly linked to stormy weather. Through the action of group velocity, a new storm can appear “spontaneously” downstream of an existing chain of storms. The propagation of energy is more rapid than the movement of the individual storms.

**Extraction of the envelope **

The envelope of a wave packet may be extracted using ideas based on the Hilbert transform. For full details, see Bracewell (2000, pp. 359-367). Let be a function on a periodic domain . We perform the following operations in sequence:

- Compute the Fourier coefficients: .
- Set the coefficients to zero for negative index: where is the Heavyside sequence.
- Compute the inverse transform: .
- Double and take the absolute value: .

In words, we calculate the Fourier series, throw away the negative frequencies, invert, double and take the absolute value.

A simple example illustrates the technique. Suppose . There are just two non-vanishing terms in the Fourier series: . Elimination of the negative frequency part leaves and twice the absolute value of this is , as expected.

We consider a signal with two wave packets, each with a different carrier frequency and envelope. These are shown in the figure below. The envelope calculated as above is also shown.

The envelope extraction may be combined with low-pass or band-pass filtering by replacing the Heavyside function by a suitable masking function, for example

which eliminates all components except in the frequency band . We show below the effects of band-passing the signal with wavenubers in the range (top) and (bottom), isolating the two packets.

**The Hilbert transform**

The theoretical explanation of the envelope extraction method is given in Bracewell (pp. 359-367). The Hilbert transform of a function is defined as

where the Cauchy principal value of the integral is intended. For a real signal, there is an associated complex signal called *the analytical signal*, defined by

and the amplitude of this signal contains the information about the envelope.

References

Bracewell, R, 2000: *The Fourier Transform and its Applications.* Third Edn., McGraw-Hill, New York. 636pp.

You must be logged in to post a comment.