The Fourier Transform on the Real Line

Part 13, Chapter 13: Fourier Series and Transforms

Learning objectives

  • Write the Fourier transform f^(ξ)\hat{f}(\xi) and the inverse Fourier transform
  • Apply the convolution theorem and the time-shift property
  • State and use Plancherel's theorem (the continuous analog of Parseval)
  • Recognize the Gaussian as the canonical fixed point f^=f\hat{f} = f

What replaces a Fourier series when the function is not periodic but lives on all of mathbbR\mathbb{R}? The discrete index nn becomes a continuous variable xi\xi, the sum becomes an integral, and the formulas are nearly identical, with one immensely useful new property. Convolution in space (which describes blurring, smoothing, the action of any linear time-invariant filter) becomes pointwise multiplication in frequency. That single fact, the convolution theorem, is why the Fourier transform is the central tool of signal processing, image analysis, quantum mechanics, and partial differential equations.

Definition and inversion

For a function f:mathbbRtomathbbCf: \mathbb{R} \to \mathbb{C} in L1L^1 (so the integrals exist), define the Fourier transform:

hatf(xi)=displaystyleintinftyinftyf(x)e2piixix,dx\hat{f}(\xi) = \displaystyle\int_{-\infty}^{\infty} f(x) e^{-2\pi i \xi x} \, dxinftyinftyf(x)e2piixix,dx

The inverse is symmetric:

f(x)=displaystyleintinftyinftyhatf(xi)e2piixix,dxif(x) = \displaystyle\int_{-\infty}^{\infty} \hat{f}(\xi) e^{2\pi i \xi x} \, d\xiinftyinftyhatf(xi)e2piixix,dxi

(Caveat: different books use different normalisations. The convention here uses 2pi2\pi in the exponent and no leading factors; some texts put 1/sqrt2pi1/\sqrt{2\pi} in front of both formulas, others put 1/(2pi)1/(2\pi) in front of one. Any consistent choice works; just stick to one.)

The core properties

Linearity: widehatalphaf+betag=alphahatf+betahatg\widehat{\alpha f + \beta g} = \alpha \hat{f} + \beta \hat{g}. The Fourier transform is a linear operator on L1L^1 (and extends to other spaces).

Time shift: if g(x)=f(xa)g(x) = f(x - a), then hatg(xi)=e2piiaxihatf(xi)\hat{g}(\xi) = e^{-2\pi i a \xi} \hat{f}(\xi). Shifting in time multiplies the transform by a phase factor, the magnitude is unchanged.

Modulation: if g(x)=e2piiaxf(x)g(x) = e^{2\pi i a x} f(x), then hatg(xi)=hatf(xia)\hat{g}(\xi) = \hat{f}(\xi - a). Multiplying by a complex exponential in time shifts the spectrum, this is how AM radio works.

Scaling: if g(x)=f(alphax)g(x) = f(\alpha x) with alpha>0\alpha > 0, then hatg(xi)=(1/alpha)hatf(xi/alpha)\hat{g}(\xi) = (1/\alpha) \hat{f}(\xi / \alpha). Compressing in time stretches in frequency, and vice versa, the uncertainty principle in disguise.

The convolution theorem

Define the convolution of ff and gg:

(f*g)(x)=displaystyleintinftyinftyf(t)g(xt),dt(f * g)(x) = \displaystyle\int_{-\infty}^{\infty} f(t) g(x - t) \, dtinftyinftyf(t)g(xt),dt

The convolution theorem states:

widehatf*g=hatfcdothatg\widehat{f * g} = \hat{f} \cdot \hat{g}

Convolution (the messy, expensive integral) becomes pointwise multiplication (cheap and parallel) in the frequency domain. That is the whole reason the Fast Fourier Transform exists in production code: you transform, multiply, and inverse-transform, an O(NlogN)O(N \log N) procedure that would otherwise cost O(N2)O(N^2).

Plancherel's theorem

The continuous analog of Parseval:

displaystyleintinftyinftyf(x)2,dx=displaystyleintinftyinftyhatf(xi)2,dxi\displaystyle\int_{-\infty}^{\infty} |f(x)|^2 \, dx = \displaystyle\int_{-\infty}^{\infty} |\hat{f}(\xi)|^2 \, d\xiinftyinftyhatf(xi)2,dxi

The Fourier transform preserves the L2L^2 norm, so it extends from L1capL2L^1 \cap L^2 to a unitary operator on all of L2(mathbbR)L^2(\mathbb{R}).

Where this shows up
  • MRI k-space reconstruction: An MR scanner measures the Fourier transform of the spin-density map of the slice being imaged. The acquired raw data live in k-space, and the image is reconstructed by the inverse Fourier transform, this is literally the inversion formula at industrial scale.
  • Heisenberg uncertainty principle: If ff is well-concentrated in xx, then hatf\hat{f} must be spread out in xi\xi. The product of the two spreads is bounded below by a constant, the scaling property is the elementary version, and the rigorous statement is the variance inequality sigmaxsigmaxigeq1/(4pi)text(withoure2piixixtextconvention;equalshbar/2textinphysicsnotation)\sigma_x \sigma_\xi \geq 1/(4\pi)\text{ (with our }e^{-2\pi i\xi x}\text{ convention; equals }\hbar/2\text{ in physics notation)}xigeq1/(4pi)text(withoure2piixixtextconvention;equalshbar/2textinphysicsnotation).
  • Spectral analysis of seismic data: Earthquake records are non-periodic. To estimate the dominant frequencies (and hence subsurface structure), seismologists compute hatf(xi)2|\hat{f}(\xi)|^2, the power spectrum, via FFT.
  • Heat equation on the line: The heat equation ut=uxxu_t = u_{xx}xx on mathbbR\mathbb{R} is solved by Fourier transforming in xx. The PDE becomes the ODE hatut=4pi2xi2hatu\hat{u}_t = -4\pi^2 \xi^2 \hat{u}t=4pi2xi2hatu, whose solution is hatu(xi,t)=e4pi2xi2thatu_0(xi)\hat{u}(\xi, t) = e^{-4\pi^2 \xi^2 t} \hat{u}_0(\xi), a Gaussian multiplier in the frequency domain.

    Pause and think: The Gaussian f(x)=epix2f(x) = e^{-\pi x^2} satisfies hatf=f\hat{f} = f (the standard Gaussian is its own Fourier transform). What is hatf(0)\hat{f}(0), and what classical integral does it equal?

    Try it

    • Compute hatf(xi)\hat{f}(\xi) for the rectangular pulse f(x)=1f(x) = 1 on x<1/2|x| < 1/2, zero elsewhere. (Answer: textsinc(xi)=sin(pixi)/(pixi)\text{sinc}(\xi) = \sin(\pi \xi)/(\pi \xi).)
    • If hatf(xi)=exi\hat{f}(\xi) = e^{-|\xi|}, find hatf(0)\hat{f}(0). What does it represent? (Hint: it equals intf(x),dx\int f(x) \, dx.)
    • Apply the time-shift property to compute the Fourier transform of f(x)=textrect(x3)f(x) = \text{rect}(x - 3), where textrect\text{rect} is the standard rectangular pulse from the first exercise.
    • If intf2,dx=7\int |f|^2 \, dx = 7, what is inthatf2,dxi\int |\hat{f}|^2 \, d\xi? (Plancherel.)
    • Use the convolution theorem: if hatf(xi)=exi\hat{f}(\xi) = e^{-|\xi|} and hatg(xi)=1/(1+xi2)\hat{g}(\xi) = 1/(1 + \xi^2), write down widehatf*g\widehat{f * g}.

    A trap to watch for

    The most common Fourier-transform error is using the wrong convention. The formulas above use e2piixixe^{-2\pi i \xi x} in the forward transform, but many physics books use eiomegaxe^{-i \omega x} instead. That switch changes the inversion constant (you pick up a 1/(2pi)1/(2\pi)) and changes the form of Plancherel's constant. Always check the convention printed on page 1 of whatever reference you are using, and never copy formulas between books without verifying their conventions match.

    What you now know

    You can compute Fourier transforms of basic signals, apply the convolution theorem to convert a hard convolution into easy multiplication, and use Plancherel to track energy across the transform. The next chapter applies all of this Fourier machinery (plus separation-of-variables techniques) to solve the three canonical PDEs: Laplace, heat, and wave.

    Mark section complete →

    References

    • Garrity, T. (2002). All the Mathematics You Missed: But Need to Know for Graduate School. Cambridge University Press, ch. 13.
    • Stein, E. M., Shakarchi, R. (2003). Fourier Analysis: An Introduction. Princeton University Press, ch. 5.
    • Folland, G. B. (1992). Fourier Analysis and Its Applications. Wadsworth & Brooks/Cole, ch. 7-8.
    • Bracewell, R. N. (1999). The Fourier Transform and Its Applications (3rd ed.). McGraw-Hill, ch. 6-9.
    • Strichartz, R. S. (2003). A Guide to Distribution Theory and Fourier Transforms. World Scientific, ch. 2.

    This page is prerendered for SEO and accessibility. The interactive widgets above hydrate on JavaScript load.