ConceptComplete

Fourier Analysis - Key Properties

The mathematical properties of Fourier transforms reveal deep connections between time and frequency domains, providing powerful tools for analyzing linear systems and wave phenomena.

Convolution Theorem and Linear Response

TheoremConvolution Theorem

The Fourier transform converts convolution to multiplication:

F[fāˆ—g](k)=F[f](k)ā‹…F[g](k)\mathcal{F}[f * g](k) = \mathcal{F}[f](k) \cdot \mathcal{F}[g](k)

where the convolution is defined as:

(fāˆ—g)(x)=āˆ«āˆ’āˆžāˆžf(x′)g(xāˆ’x′)dx′(f * g)(x) = \int_{-\infty}^{\infty} f(x')g(x - x')dx'

Similarly, multiplication in real space becomes convolution in Fourier space:

F[fā‹…g](k)=12Ļ€[F[f]āˆ—F[g]](k)\mathcal{F}[f \cdot g](k) = \frac{1}{2\pi}[\mathcal{F}[f] * \mathcal{F}[g]](k)

This theorem is fundamental in signal processing and physics. For a linear time-invariant system with impulse response h(t)h(t), the output for input f(t)f(t) is:

y(t)=(hāˆ—f)(t)ā€…ā€ŠāŸ¹ā€…ā€Šy~(ω)=h~(ω)f~(ω)y(t) = (h * f)(t) \implies \tilde{y}(\omega) = \tilde{h}(\omega)\tilde{f}(\omega)

The transfer function h~(ω)\tilde{h}(\omega) completely characterizes the system's frequency response.

ExampleRC Circuit as a Low-Pass Filter

For an RC circuit with input voltage Vin(t)V_{in}(t) and output Vout(t)V_{out}(t) across the capacitor:

τdVoutdt+Vout=Vin,τ=RC\tau\frac{dV_{out}}{dt} + V_{out} = V_{in}, \quad \tau = RC

Fourier transforming with F[dV/dt]=iωV~\mathcal{F}[dV/dt] = i\omega\tilde{V}:

V~out(ω)=11+iωτV~in(ω)=H(ω)V~in(ω)\tilde{V}_{out}(\omega) = \frac{1}{1 + i\omega\tau}\tilde{V}_{in}(\omega) = H(\omega)\tilde{V}_{in}(\omega)

The transfer function H(ω)=1/(1+iωτ)H(\omega) = 1/(1 + i\omega\tau) has magnitude:

∣H(ω)∣=11+ω2Ļ„2|H(\omega)| = \frac{1}{\sqrt{1 + \omega^2\tau^2}}

which falls off for ω>1/Ļ„\omega > 1/\tau, making this a low-pass filter with cutoff frequency ωc=1/Ļ„\omega_c = 1/\tau.

Uncertainty Relations

TheoremFourier Uncertainty Principle

For a function f(x)f(x) with Fourier transform f~(k)\tilde{f}(k), define the variances:

(Ī”x)2=∫x2∣f(x)∣2dx∫∣f(x)∣2dxāˆ’(∫x∣f(x)∣2dx∫∣f(x)∣2dx)2(\Delta x)^2 = \frac{\int x^2|f(x)|^2dx}{\int|f(x)|^2dx} - \left(\frac{\int x|f(x)|^2dx}{\int|f(x)|^2dx}\right)^2

(Ī”k)2=∫k2∣f~(k)∣2dk∫∣f~(k)∣2dkāˆ’(∫k∣f~(k)∣2dk∫∣f~(k)∣2dk)2(\Delta k)^2 = \frac{\int k^2|\tilde{f}(k)|^2dk}{\int|\tilde{f}(k)|^2dk} - \left(\frac{\int k|\tilde{f}(k)|^2dk}{\int|\tilde{f}(k)|^2dk}\right)^2

Then the uncertainty relation holds:

Ī”xā‹…Ī”k≄12\Delta x \cdot \Delta k \geq \frac{1}{2}

with equality for Gaussian wave packets.

In quantum mechanics, setting k=p/ā„k = p/\hbar gives Heisenberg's uncertainty principle:

Ī”xā‹…Ī”pā‰„ā„2\Delta x \cdot \Delta p \geq \frac{\hbar}{2}

This is not a statement about measurement but a fundamental property of wave functions: narrow localization in position requires broad spread in momentum, and vice versa.

Sampling Theorem

TheoremNyquist-Shannon Sampling Theorem

If a function f(t)f(t) has Fourier transform f~(ω)\tilde{f}(\omega) that vanishes for āˆ£Ļ‰āˆ£>ωmax⁔|\omega| > \omega_{\max} (band-limited), then f(t)f(t) is completely determined by its samples at intervals Ī”t=Ļ€/ωmax⁔\Delta t = \pi/\omega_{\max}:

f(t)=āˆ‘n=āˆ’āˆžāˆžf(nĪ”t) sinc(tāˆ’nĪ”tĪ”t)f(t) = \sum_{n=-\infty}^{\infty} f(n\Delta t)\,\text{sinc}\left(\frac{t - n\Delta t}{\Delta t}\right)

where sinc(x)=sin⁔(Ļ€x)/(Ļ€x)\text{sinc}(x) = \sin(\pi x)/(\pi x). The Nyquist frequency is ωN=Ļ€/Ī”t=2ωmax⁔\omega_N = \pi/\Delta t = 2\omega_{\max}.

This theorem is crucial for digital signal processing: signals must be sampled at least twice per period of the highest frequency component to avoid aliasing.

ExampleAliasing in Digital Audio

Human hearing ranges to about 20 kHz. CD audio samples at 44.1 kHz, giving a Nyquist frequency of 22.05 kHz. Frequencies above this fold back into the audible range (aliasing), so an anti-aliasing filter removes frequencies above 20 kHz before sampling.

Multi-dimensional Fourier Transforms

DefinitionFourier Transform in $\mathbb{R}^n$

For a function f(x)f(\mathbf{x}) on Rn\mathbb{R}^n:

f~(k)=∫Rnf(x)eāˆ’ikā‹…xdnx\tilde{f}(\mathbf{k}) = \int_{\mathbb{R}^n} f(\mathbf{x})e^{-i\mathbf{k} \cdot \mathbf{x}}d^nx

f(x)=1(2Ļ€)n∫Rnf~(k)eikā‹…xdnkf(\mathbf{x}) = \frac{1}{(2\pi)^n}\int_{\mathbb{R}^n} \tilde{f}(\mathbf{k})e^{i\mathbf{k} \cdot \mathbf{x}}d^nk

The Laplacian becomes multiplication: F[āˆ‡2f](k)=āˆ’k2f~(k)\mathcal{F}[\nabla^2 f](\mathbf{k}) = -k^2\tilde{f}(\mathbf{k}) where k=∣k∣k = |\mathbf{k}|.

ExampleGreen's Function for Poisson Equation

To solve āˆ‡2Ļ•(x)=āˆ’Ļ(x)/ϵ0\nabla^2\phi(\mathbf{x}) = -\rho(\mathbf{x})/\epsilon_0, Fourier transform both sides:

āˆ’k2Ļ•~(k)=āˆ’Ļ~(k)/ϵ0-k^2\tilde{\phi}(\mathbf{k}) = -\tilde{\rho}(\mathbf{k})/\epsilon_0

Thus Ļ•~(k)=ρ~(k)/(k2ϵ0)\tilde{\phi}(\mathbf{k}) = \tilde{\rho}(\mathbf{k})/(k^2\epsilon_0). For a point charge ρ(x)=qĪ“3(x)\rho(\mathbf{x}) = q\delta^3(\mathbf{x}), ρ~(k)=q\tilde{\rho}(\mathbf{k}) = q:

Ļ•(x)=qϵ0(2Ļ€)3∫eikā‹…xk2d3k=q4πϵ0∣x∣\phi(\mathbf{x}) = \frac{q}{\epsilon_0(2\pi)^3}\int \frac{e^{i\mathbf{k} \cdot \mathbf{x}}}{k^2}d^3k = \frac{q}{4\pi\epsilon_0|\mathbf{x}|}

recovering Coulomb's law.

RemarkFast Fourier Transform (FFT)

The discrete Fourier transform of NN points naively requires O(N2)O(N^2) operations. The FFT algorithm reduces this to O(Nlog⁔N)O(N\log N) by exploiting symmetries, revolutionizing digital signal processing, crystallography, and numerical PDE solvers.

These properties make Fourier analysis an indispensable tool for understanding oscillatory and wave phenomena across all areas of physics.