ConceptComplete

Fourier Methods - Examples and Constructions

Fourier methods provide explicit solutions to fundamental PDEs and enable sophisticated analysis of dispersion, decay, and regularity properties.

ExampleSolving the Wave Equation

For the wave equation utt=c2βˆ‡2uu_{tt} = c^2\nabla^2 u with initial data u(x,0)=f(x)u(x,0) = f(x), ut(x,0)=g(x)u_t(x,0) = g(x):

Taking Fourier transform in space: βˆ‚2u^βˆ‚t2=βˆ’4Ο€2c2∣ξ∣2u^\frac{\partial^2\hat{u}}{\partial t^2} = -4\pi^2c^2|\xi|^2\hat{u}

Solution: u^(ΞΎ,t)=f^(ΞΎ)cos⁑(2Ο€c∣ξ∣t)+g^(ΞΎ)2Ο€c∣ξ∣sin⁑(2Ο€c∣ξ∣t)\hat{u}(\xi,t) = \hat{f}(\xi)\cos(2\pi c|\xi|t) + \frac{\hat{g}(\xi)}{2\pi c|\xi|}\sin(2\pi c|\xi|t)

Inverse transforming gives the solution. The oscillatory nature in frequency space reflects wave propagation in physical space.

ExampleSchrΓΆdinger Equation

The linear SchrΓΆdinger equation iβ„βˆ‚Οˆβˆ‚t=βˆ’β„22mβˆ‡2ψi\hbar\frac{\partial\psi}{\partial t} = -\frac{\hbar^2}{2m}\nabla^2\psi becomes in Fourier space: iβ„βˆ‚Οˆ^βˆ‚t=2Ο€2ℏ2m∣ξ∣2ψ^i\hbar\frac{\partial\hat{\psi}}{\partial t} = \frac{2\pi^2\hbar^2}{m}|\xi|^2\hat{\psi}

giving ψ^(ΞΎ,t)=ψ^0(ΞΎ)eβˆ’2Ο€2iβ„βˆ£ΞΎβˆ£2t/m\hat{\psi}(\xi,t) = \hat{\psi}_0(\xi)e^{-2\pi^2i\hbar|\xi|^2t/m}.

The dispersion relation Ο‰=β„βˆ£ΞΎβˆ£22m\omega = \frac{\hbar|\xi|^2}{2m} shows different frequencies travel at different speeds, causing wave packet spreading.

ExampleFast Fourier Transform (FFT)

Computing discrete Fourier transforms naively requires O(N2)O(N^2) operations. The FFT algorithm reduces this to O(Nlog⁑N)O(N\log N) by recursively decomposing the transform: f^k=βˆ‘j=0Nβˆ’1fjeβˆ’2Ο€ijk/N\hat{f}_k = \sum_{j=0}^{N-1} f_j e^{-2\pi ijk/N}

This computational breakthrough enables real-time signal processing and efficient numerical PDE solvers.

Remark

Gibbs Phenomenon: Fourier series of discontinuous functions exhibit oscillations near jumps that don't disappear as more terms are added. The overshoot is approximately 9% of the jump height, reflecting the difficulty of approximating discontinuities with smooth basis functions.

ExamplePlancherel Theorem Applications

For the Dirichlet problem on the half-space {xn>0}\{x_n > 0\}: βˆ‡2u=0,u(xβ€²,0)=g(xβ€²)\nabla^2 u = 0, \quad u(x', 0) = g(x')

Taking Fourier transform in xβ€²=(x1,…,xnβˆ’1)x' = (x_1, \ldots, x_{n-1}): βˆ‚2u^βˆ‚xn2=4Ο€2βˆ£ΞΎβ€²βˆ£2u^\frac{\partial^2\hat{u}}{\partial x_n^2} = 4\pi^2|\xi'|^2\hat{u}

Solution: u^(ΞΎβ€²,xn)=g^(ΞΎβ€²)eβˆ’2Ο€βˆ£ΞΎβ€²βˆ£xn\hat{u}(\xi', x_n) = \hat{g}(\xi')e^{-2\pi|\xi'|x_n}

Inverse transforming gives the Poisson integral formula for the half-space.

These examples demonstrate how Fourier methods reduce PDEs to ODEs or algebraic equations, enabling explicit solutions and detailed analysis of solution behavior.