TheoremComplete

Fourier Analysis - Applications

Fourier methods provide powerful techniques for solving partial differential equations, analyzing signals, and understanding quantum phenomena through spectral decomposition.

TheoremSolution of PDEs via Separation of Variables

For linear PDEs with constant coefficients, separation of variables combined with Fourier series yields complete solutions. Consider the heat equation on [0,L][0, L] with boundary conditions u(0,t)=u(L,t)=0u(0,t) = u(L,t) = 0 and initial condition u(x,0)=f(x)u(x,0) = f(x):

βˆ‚uβˆ‚t=Ξ±βˆ‚2uβˆ‚x2\frac{\partial u}{\partial t} = \alpha\frac{\partial^2 u}{\partial x^2}

The solution is:

u(x,t)=βˆ‘n=1∞bneβˆ’Ξ±(nΟ€/L)2tsin⁑(nΟ€xL)u(x,t) = \sum_{n=1}^{\infty} b_n e^{-\alpha(n\pi/L)^2 t}\sin\left(\frac{n\pi x}{L}\right)

where bn=2L∫0Lf(x)sin⁑(nΟ€x/L)dxb_n = \frac{2}{L}\int_0^L f(x)\sin(n\pi x/L)dx are the Fourier sine coefficients of f(x)f(x).

Each mode decays exponentially with rate Ξ±(nΟ€/L)2\alpha(n\pi/L)^2, so higher frequency components (larger nn) decay faster, leading to smoothing of initial data.

ExampleQuantum Particle in a Box

The time-independent Schrâdinger equation for a particle in an infinite square well V(x)=0V(x) = 0 for 0<x<L0 < x < L, V(x)=∞V(x) = \infty otherwise:

βˆ’β„22md2ψdx2=Eψ,ψ(0)=ψ(L)=0-\frac{\hbar^2}{2m}\frac{d^2\psi}{dx^2} = E\psi, \quad \psi(0) = \psi(L) = 0

has solutions:

ψn(x)=2Lsin⁑(nΟ€xL),En=n2Ο€2ℏ22mL2\psi_n(x) = \sqrt{\frac{2}{L}}\sin\left(\frac{n\pi x}{L}\right), \quad E_n = \frac{n^2\pi^2\hbar^2}{2mL^2}

The general time-dependent wave function is:

Ξ¨(x,t)=βˆ‘n=1∞cnψn(x)eβˆ’iEnt/ℏ\Psi(x,t) = \sum_{n=1}^{\infty} c_n \psi_n(x)e^{-iE_n t/\hbar}

where cn=∫0Lψn(x)Ψ(x,0)dxc_n = \int_0^L \psi_n(x)\Psi(x,0)dx from the initial condition. This is a Fourier sine series in both spatial and energy representations.

TheoremGreen's Function via Fourier Transform

For the differential operator L=βˆ’βˆ‡2+m2\mathcal{L} = -\nabla^2 + m^2 in R3\mathbb{R}^3, the Green's function G(x)G(\mathbf{x}) satisfying:

(βˆ’βˆ‡2+m2)G(x)=Ξ΄3(x)(-\nabla^2 + m^2)G(\mathbf{x}) = \delta^3(\mathbf{x})

can be found by Fourier transform. With G~(k)\tilde{G}(\mathbf{k}):

(k2+m2)G~(k)=1β€…β€ŠβŸΉβ€…β€ŠG~(k)=1k2+m2(k^2 + m^2)\tilde{G}(\mathbf{k}) = 1 \implies \tilde{G}(\mathbf{k}) = \frac{1}{k^2 + m^2}

Inverting (using spherical coordinates with kβ‹…x=krcos⁑θ\mathbf{k} \cdot \mathbf{x} = kr\cos\theta):

G(x)=1(2Ο€)3∫eikβ‹…xk2+m2d3k=14Ο€βˆ£x∣eβˆ’m∣x∣G(\mathbf{x}) = \frac{1}{(2\pi)^3}\int \frac{e^{i\mathbf{k} \cdot \mathbf{x}}}{k^2 + m^2}d^3k = \frac{1}{4\pi|\mathbf{x}|}e^{-m|\mathbf{x}|}

This is the Yukawa potential, describing screened Coulomb interactions with screening length 1/m1/m.

ExampleElectromagnetic Wave Propagation in Plasma

In a cold plasma, the dispersion relation for electromagnetic waves is:

Ο‰2=Ο‰p2+c2k2\omega^2 = \omega_p^2 + c^2k^2

where Ο‰p=ne2/(Ο΅0m)\omega_p = \sqrt{ne^2/(\epsilon_0 m)} is the plasma frequency. A wave packet E(x,0)=E0eβˆ’x2/(2Οƒ2)E(x,0) = E_0 e^{-x^2/(2\sigma^2)} decomposes as:

E(x,t)=12Ο€βˆ«E~0(k)ei(kxβˆ’Ο‰(k)t)dkE(x,t) = \frac{1}{2\pi}\int \tilde{E}_0(k)e^{i(kx - \omega(k)t)}dk

The group velocity vg=dω/dk=c2k/ωv_g = d\omega/dk = c^2k/\omega differs from the phase velocity vp=ω/kv_p = \omega/k, causing dispersion: the packet spreads as different frequencies travel at different speeds.

TheoremSpectral Theorem for Self-Adjoint Operators

In quantum mechanics, observables are represented by self-adjoint operators. The spectral theorem states that any self-adjoint operator A^\hat{A} can be diagonalized:

A^=∫aβ€‰βˆ£a⟩⟨aβˆ£β€‰dΞΌ(a)\hat{A} = \int a\,|a\rangle\langle a|\,d\mu(a)

where ∣a⟩|a\rangle are eigenstates with eigenvalue aa, and dμ(a)d\mu(a) is a spectral measure. For discrete spectra:

A^=βˆ‘nan∣an⟩⟨an∣\hat{A} = \sum_n a_n|a_n\rangle\langle a_n|

The wave function expands as ∣ψ⟩=βˆ‘ncn∣an⟩|\psi\rangle = \sum_n c_n|a_n\rangle (Fourier series in the eigenbasis), with measurement probabilities ∣cn∣2|c_n|^2 satisfying βˆ‘n∣cn∣2=1\sum_n |c_n|^2 = 1 (Parseval).

ExampleTime-Frequency Analysis and Wavelets

The Gabor transform (short-time Fourier transform) uses windowed Fourier basis:

Gf(t,Ο‰)=∫f(Ο„)w(Ο„βˆ’t)eβˆ’iωτdΟ„G_f(t,\omega) = \int f(\tau)w(\tau - t)e^{-i\omega\tau}d\tau

where w(t)w(t) is a window function (often Gaussian). This provides time-frequency localization but has fixed resolution Ξ”t⋅Δωβ‰₯1/2\Delta t \cdot \Delta\omega \geq 1/2.

Wavelet transforms use dilated and translated basis functions:

Wf(a,b)=∫f(t)Οˆβˆ—(tβˆ’ba)dtW_f(a,b) = \int f(t)\psi^*\left(\frac{t-b}{a}\right)dt

where aa is the scale and bb is the position. This provides variable time-frequency resolution, crucial for analyzing signals with features at multiple scales (e.g., gravitational wave detection, image compression).

RemarkFourier Analysis in Quantum Field Theory

In quantum field theory, fields are expanded in Fourier modes (creation and annihilation operators):

Ο•(x,t)=∫d3k(2Ο€)312Ο‰k[akei(kxβˆ’Ο‰kt)+ak†eβˆ’i(kxβˆ’Ο‰kt)]\phi(x,t) = \int \frac{d^3k}{(2\pi)^3}\frac{1}{\sqrt{2\omega_k}}\left[a_k e^{i(kx - \omega_k t)} + a_k^\dagger e^{-i(kx - \omega_k t)}\right]

where aka_k annihilates a particle with momentum kk. The vacuum ∣0⟩|0\rangle satisfies ak∣0⟩=0a_k|0\rangle = 0, and multi-particle states are built by applying ak†a_k^\dagger. This is a Fock space construction based on Fourier decomposition.

These applications demonstrate that Fourier analysis is not merely a mathematical technique but the natural language for describing linear systems, waves, and quantum phenomena across physics.