7.2 The Fourier transform

The Fourier series of 7.1 handles functions that repeat with period TT. Most signals in real applications don’t repeat — a footstep, a single struck note, a Gaussian pulse, the pressure history of a passing aircraft. For these we need the Fourier transform, which is what the Fourier series becomes when the period goes to infinity.

This lesson develops the transform, surveys its operational properties (linearity, time-shift, scaling, differentiation), tabulates the most-useful pairs, and ends with the uncertainty principle that links the widths of ff and f~\tilde f.

Limit of a Fourier series

Take a function f(t)f(t) that is zero outside [T/2,T/2][-T/2, T/2] — a transient of finite extent. Make a periodic extension of ff with period TT. As long as TT is much larger than the support of ff, the periodic extension looks just like ff inside one period.

The Fourier series of the periodic extension is

f(t)  =  n=cnei2πnt/T,cn  =  1TT/2T/2f(t)ei2πnt/Tdt.f(t) \;=\; \sum_{n=-\infty}^{\infty} c_n\, e^{i\, 2\pi n t / T}, \qquad c_n \;=\; \frac{1}{T} \int_{-T/2}^{T/2} f(t)\, e^{-i\, 2\pi n t / T}\, dt.

Define the continuous frequency ωn2πn/T\omega_n \equiv 2\pi n / T. Adjacent values of ωn\omega_n differ by Δω=2π/T\Delta \omega = 2\pi / T. As TT \to \infty, Δω0\Delta \omega \to 0 and the discrete set {ωn}\{\omega_n\} fills up the real line.

From series to transform: the T → ∞ limit Derivation

Let f~(ω)Tcn\tilde f(\omega) \equiv T\, c_n when ω=ωn\omega = \omega_n. (Take the limit-preserving rescaling: as TT grows, the coefficients shrink like 1/T1/T, so multiplying by TT keeps the quantity finite.)

The coefficient formula becomes

f~(ωn)  =  T/2T/2f(t)eiωntdt.\tilde f(\omega_n) \;=\; \int_{-T/2}^{T/2} f(t)\, e^{-i \omega_n t}\, dt.

In the limit TT \to \infty this is the Fourier transform:

  f~(ω)    f(t)eiωtdt.  \boxed{\;\tilde f(\omega) \;\equiv\; \int_{-\infty}^{\infty} f(t)\, e^{-i \omega t}\, dt.\;}

Conversely, the series formula

f(t)  =  ncneiωnt  =  nf~(ωn)eiωntΔω2πf(t) \;=\; \sum_n c_n\, e^{i \omega_n t} \;=\; \sum_n \tilde f(\omega_n)\, e^{i \omega_n t}\, \frac{\Delta \omega}{2\pi}

(using 1/T=Δω/(2π)1/T = \Delta \omega / (2\pi)) becomes, in the limit, a Riemann sum that converges to the inverse Fourier transform:

  f(t)  =  12πf~(ω)eiωtdω.  \boxed{\;f(t) \;=\; \frac{1}{2\pi} \int_{-\infty}^{\infty} \tilde f(\omega)\, e^{i \omega t}\, d\omega.\;}

The discrete sum over harmonics has become a continuous integral over frequency. The period TT has dropped out; ω\omega is now a continuous variable, not a discrete set of harmonics.

The factor of 1/(2π)1/(2\pi) on the inverse transform is a convention. Different textbooks split it differently:

All three are mathematically equivalent; they shuffle factors of 2π2\pi around. Pick one and stick with it. When reading a textbook or paper, always check which convention it uses before applying any formula.

A function and its transform

f(t) — time|f̃(ω)| — frequency
f(t):

Time width Δt ≈ 1.00 · Frequency width Δω ≈ 1.00 · product Δt·Δω ≈ 1.00 (uncertainty principle: Δt·Δω ≥ ½).

Pick a function family — Gaussian, rectangle, windowed cosine, one-sided exponential — slide the width parameter, and watch f(t)f(t) and f~(ω)|\tilde f(\omega)| negotiate. The reciprocal-width relationship is universal: making ff narrower in time makes f~\tilde f broader in frequency, and vice versa. The Gaussian is the unique (up to scaling) function that’s its own Fourier transform — the “fixed point” of the transform operator, and the function that saturates the uncertainty bound.

Operational properties

The transform’s algebraic properties are what make it indispensable for solving linear differential equations. A short list (proofs are routine substitutions):

| Property | Time-domain | Frequency-domain | |---|---|---| | Linearity | af+bga f + b g | af~+bg~a \tilde f + b \tilde g | | Time shift | f(tt0)f(t - t_0) | eiωt0f~(ω)e^{-i \omega t_0}\, \tilde f(\omega) | | Frequency shift / modulation | eiω0tf(t)e^{i \omega_0 t}\, f(t) | f~(ωω0)\tilde f(\omega - \omega_0) | | Time scaling | f(at)f(a t) | 1af~(ω/a)\frac{1}{\lvert a \rvert}\, \tilde f(\omega / a) | | Differentiation | f˙(t)\dot f(t) | iωf~(ω)i \omega\, \tilde f(\omega) | | Higher derivative | f(n)(t)f^{(n)}(t) | (iω)nf~(ω)(i \omega)^n\, \tilde f(\omega) | | Multiplication by tt | tf(t)t f(t) | idf~/dωi\, d\tilde f / d\omega | | Complex conjugation | f(t)f^*(t) | f~(ω)\tilde f^*(-\omega) |

The differentiation property is the cash value of the entire transform for differential equations. A constant-coefficient linear ODE in time becomes a polynomial equation in ω\omega. A constant-coefficient linear PDE in (x,t)(x, t) becomes an algebraic equation in (k,ω)(k, \omega). The transform converts calculus into algebra.

Useful Fourier-transform pairs

| f(t)f(t) | f~(ω)\tilde f(\omega) | |---|---| | Rectangular pulse, width TT centred at 0 | Tsinc(ωT/2)T\, \mathrm{sinc}(\omega T / 2) | | Gaussian, et2/2σ2e^{-t^2 / 2 \sigma^2} | σ2πeσ2ω2/2\sigma \sqrt{2\pi}\, e^{-\sigma^2 \omega^2 / 2} | | One-sided exponential decay, eαtθ(t)e^{-\alpha t}\, \theta(t) | 1/(α+iω)1 / (\alpha + i \omega) | | Two-sided exponential, eαte^{-\alpha \lvert t \rvert} | 2α/(α2+ω2)2\alpha / (\alpha^2 + \omega^2) — a Lorentzian | | Pure sinusoid, cos(ω0t)\cos(\omega_0 t) | π[δ(ωω0)+δ(ω+ω0)]\pi\, [\delta(\omega - \omega_0) + \delta(\omega + \omega_0)] | | Pure sinusoid, sin(ω0t)\sin(\omega_0 t) | iπ[δ(ωω0)δ(ω+ω0)]-i\pi\, [\delta(\omega - \omega_0) - \delta(\omega + \omega_0)] | | Delta function, δ(t)\delta(t) | 11 — flat spectrum | | Constant, 11 | 2πδ(ω)2\pi\, \delta(\omega) — pure DC | | Sign function, sgn(t)\mathrm{sgn}(t) | 2/(iω)2 / (i \omega) | | Heaviside step, θ(t)\theta(t) | πδ(ω)+1/(iω)\pi \delta(\omega) + 1/(i\omega) |

A few observations worth pulling out:

Fourier transform of a decaying exponential (damped oscillation) Worked Example

Let f(t)=eαtθ(t)f(t) = e^{-\alpha t}\,\theta(t) with α=3s1\alpha = 3\,\text{s}^{-1} (a sound that decays with time constant 1/31/3 s). Compute f~(ω)\tilde f(\omega):

f~(ω)=0eαteiωtdt=0e(α+iω)tdt=1α+iω.\tilde f(\omega) = \int_0^\infty e^{-\alpha t}\, e^{-i\omega t}\, dt = \int_0^\infty e^{-(\alpha + i\omega)t}\, dt = \frac{1}{\alpha + i\omega}.

The magnitude spectrum is f~(ω)=1/α2+ω2=1/9+ω2|\tilde f(\omega)| = 1/\sqrt{\alpha^2 + \omega^2} = 1/\sqrt{9 + \omega^2}. At ω=0\omega = 0: f~=1/3|\tilde f| = 1/3. At ω=3\omega = 3: f~=1/(32)0.236|\tilde f| = 1/(3\sqrt{2}) \approx 0.236, half-power point. The bandwidth (half-power) is Δω=α=3rad/s\Delta\omega = \alpha = 3\,\text{rad/s}, confirming that faster decay \Rightarrow broader spectrum.

Fourier transform of a Gaussian pulse Worked Example

Let f(t)=et2/(2σ2)f(t) = e^{-t^2/(2\sigma^2)} with σ=2ms\sigma = 2\,\text{ms} (a brief acoustic click). Using the Gaussian FT pair:

f~(ω)=σ2πeσ2ω2/2=22π×103e2×106ω2.\tilde f(\omega) = \sigma\sqrt{2\pi}\, e^{-\sigma^2 \omega^2/2} = 2\sqrt{2\pi}\times 10^{-3}\, e^{-2\times10^{-6}\,\omega^2}.

The spectral width (where the amplitude falls to 1/e1/e) is Δω=1/σ=500rad/s\Delta\omega = 1/\sigma = 500\,\text{rad/s}, i.e. Δf80Hz\Delta f \approx 80\,\text{Hz}. Doubling σ\sigma to 4 ms halves the bandwidth to 40Hz\sim 40\,\text{Hz} — the uncertainty relation at work.

The uncertainty principle

A function and its Fourier transform cannot both be sharply localised. With suitable definitions of widths Δt\Delta t and Δω\Delta \omega (typically the standard deviations of f2|f|^2 and f~2|\tilde f|^2),

  ΔtΔω    12.  \boxed{\;\Delta t \cdot \Delta \omega \;\geq\; \frac{1}{2}.\;}

Equality holds only for Gaussians; every other function strictly exceeds it. The proof is a Cauchy–Schwarz argument on the inner product of tf(t)t f(t) with f˙(t)\dot f(t).

This is the same inequality as Heisenberg’s ΔxΔp/2\Delta x \cdot \Delta p \geq \hbar / 2 in quantum mechanics, via p=kp = \hbar k and momentum-eigenstates-are-plane-waves. The uncertainty principle is not a fact about quantum physics; it is a fact about the Fourier transform applied to a function-space inner product. Quantum mechanics inherits it because its states live in such a space.

The practical consequence for acoustics: a short pulse (small Δt\Delta t) has a wide spectrum (large Δω\Delta \omega). A sharply defined pitch (small Δω\Delta \omega) requires listening over a long window (large Δt\Delta t). This is the central design tension of spectrogram analysis (Sound 8.2) and is what limits how precisely a click and a pitch can simultaneously be discerned.

What we use this for

The transform is the workhorse of frequency-domain physics:

The next lesson, 7.3 — Convolution and Parseval, develops the two algebraic identities that turn the Fourier transform from a calculation tool into the foundation of linear-systems theory.