advanced26 min studyLesson 23 in Mathematics

Fourier analysis

Translate between position and frequency while keeping normalization, convergence, and boundaries explicit.

On this page

What you’ll understand

  • State a consistent Fourier transform convention and its inverse.
  • Explain why derivatives and convolution simplify in frequency space.
  • Separate finite approximations, L2 convergence, and pointwise reconstruction.

A complicated waveform can be easier to understand through how much of each oscillation it contains. Fourier analysis turns translation into phase, differentiation into multiplication, and many linear evolution equations into independent modes. These transformations underpin signal processing, quantum momentum space, and spectral numerical methods. Their success does not remove analytic obligations: normalization factors, boundary conditions, and the sense in which an infinite reconstruction converges are all part of the calculation.

Choose the convention before doing algebra

For a sufficiently well-behaved function on the real line, define its transform by integrating against a complex oscillation. We use angular frequency k, a minus sign in the forward exponential, and a factor 1/(2π) in the inverse. Other conventions distribute factors differently or use cycles per unit length. None is inherently more correct, but mixing them produces incorrect amplitudes. Start with Schwartz functions, whose smoothness and rapid decay make the manipulations below safe.

A periodic function uses a discrete Fourier series because allowable modes must fit the period. A function on the full line uses a continuous transform. A discrete Fourier transform is another object: a finite invertible linear map on sampled data. Interpreting that data as samples of a continuous signal requires sampling assumptions. Aliasing can identify distinct continuous frequencies on the same sample grid, so a high-resolution plot of finite data does not recover frequencies that were never distinguished.

f^(k)=∫Re−ikxf(x) dx,f(x)=12π∫Reikxf^(k) dk\widehat f(k)=\int_{\mathbb R}e^{-ikx}f(x)\,dx,\qquad f(x)=\frac1{2\pi}\int_{\mathbb R}e^{ikx}\widehat f(k)\,dk
These formulas hold directly for Schwartz functions. Broader L1/L2 classes require suitable inversion hypotheses or an L2 interpretation rather than automatic pointwise equality.↔ Scroll wide equations horizontally.

Operators become multipliers

Integration by parts gives the transform of a derivative: multiplication by ik, provided the boundary term vanishes. Two derivatives give −k². This makes constant-coefficient differential equations especially simple. The heat equation becomes a first-order ordinary differential equation for each k, and each mode decays at a rate proportional to k². Short wavelengths disappear quickly because they have large frequency, explaining the smoothing effect without following every point of the original profile.

Convolution combines shifted copies of one function weighted by another. Its transform is the product of the transforms under our convention. Consequently a translation-invariant linear filter can be described by a multiplier in frequency space. Division by a multiplier requires caution where it vanishes or becomes small. An inverse filter may amplify tiny measurement errors, and a Green function may need a distributional interpretation or boundary prescription. Algebraic simplification can reveal ill conditioning rather than cure it.

f′^(k)=ikf^(k),f∗g^(k)=f^(k)g^(k),u^(k,t)=e−κk2tu^(k,0)\widehat{f'}(k)=ik\widehat f(k),\qquad\widehat{f*g}(k)=\widehat f(k)\widehat g(k),\qquad\widehat u(k,t)=e^{-\kappa k^2t}\widehat u(k,0)
Take Schwartz data for the first identities. The last formula solves the full-line heat equation u_t=κu_xx for κ>0 and t≥0 with suitable initial data.↔ Scroll wide equations horizontally.

Energy can converge when a pointwise picture is subtle

Plancherel's theorem extends the transform from a convenient dense class to L2 while preserving its norm up to the chosen normalization. This means square-integrable signals can be handled even when the forward integral is not an ordinary absolutely convergent integral. Their transforms, like the original L2 elements, are defined up to changes on null sets. Equality of integrated energy does not automatically prescribe a value at every point or guarantee pointwise convergence of arbitrary partial reconstructions.

For periodic piecewise smooth data, Fourier partial sums approach the midpoint of a jump and exhibit Gibbs overshoot near it. Increasing the number of modes narrows the troublesome region without eliminating the limiting relative overshoot. Smoother summation methods, such as Fejér averaging, change convergence behavior. This is useful when interpreting field truncations and spectral simulations: the approximation rule and the function class both matter. A finite spectral cutoff is a controlled representation choice, not an exact description of every scale.

∫R∣f(x)∣2 dx=12π∫R∣f^(k)∣2 dk\int_{\mathbb R}|f(x)|^2\,dx=\frac1{2\pi}\int_{\mathbb R}|\widehat f(k)|^2\,dk
Plancherel holds for L2 functions using the transform's L2 extension. The factor 1/(2π) follows from the stated angular-frequency convention.↔ Scroll wide equations horizontally.

A heat mode is an eigenfunction

On a periodic interval of length 2π, solve ut=κuxxu_t=\kappa u_{xx} with u(x,0)=3cos⁡(2x)u(x,0)=3\cos(2x) and κ>0.

  1. The periodic basis contains cos(nx) for integer n. The initial profile has only the n=2 cosine mode, so try u(x,t)=a(t)cos⁡(2x)u(x,t)=a(t)\cos(2x).

  2. Two x derivatives produce uxx=−4a(t)cos⁡(2x)u_{xx}=-4a(t)\cos(2x). Comparing the coefficients gives a′(t)=−4κa(t)a'(t)=-4\kappa a(t) with a(0)=3.

  3. The scalar equation gives a(t)=3e−4κta(t)=3e^{-4\kappa t}. Its amplitude decays without changing the wavelength because the linear equation does not couple distinct modes.

  4. The squared L2 norm starts at 9π and decays as 9πe−8κt9\pi e^{-8\kappa t}. This is consistent with diffusion dissipating the nonconstant part of the profile rather than conserving its energy.

The result

u(x,t)=3e−4κtcos⁡(2x)u(x,t)=3e^{-4\kappa t}\cos(2x). The mode with twice the frequency decays four times as fast as the n=1 mode.

Explore the model

Interactive lab

Interactive lab · spectral

Build a profile from finitely many modes

Change mode count N and coefficient decay p in the finite sum Σ sin(nx−t)/nᵖ. More modes add shorter wavelengths, while larger p suppresses them more strongly. This explicit family demonstrates truncation and scale; it does not transform arbitrary input or prove reconstruction at jumps.

Build a profile from finitely many modesA finite sum of sin(nx − t)/nᵖ. Playback changes t; increasing the cutoff includes finer spatial scales. This is not a proof about a continuum limit or renormalization. Current readings: Sum at x = π/2: 0.867; Smallest wavelength: 2π / 5; Highest-mode coefficient: 0.2. All parameters are adjustable with the labeled controls below.01.573.144.716.28−2.56−1.2801.282.56x · radiansFinite sine sum
Finite sumHighest included mode
Sum at x = π/2
0.867
Smallest wavelength
2π / 5
Highest-mode coefficient
0.2
t = 0 s

A finite sum of sin(nx − t)/nᵖ. Playback changes t; increasing the cutoff includes finer spatial scales. This is not a proof about a continuum limit or renormalization. Parameter changes pause playback and recompute the model at the current elapsed time. Reset restores both parameters and returns time to zero. Playback starts paused. Adjust either parameter, then choose Play or Step.

Look forA mode expansion explains scale, but its truncation and convergence rule must remain visible.

Check your understanding

With the stated convention, what is the Fourier transform of the second derivative of a Schwartz function?

Try it yourself

A periodic heat profile begins as 2+sin⁡x+sin⁡(3x)2+\sin x+\sin(3x). Write its later profile and identify its long-time limit.

Take a moment to make a prediction or work it out on paper. Open a hint when you need one.

Give me a hint

Treat the constant as frequency zero and evolve each nonzero mode independently.

Compare with the solution

The zero mode has no second derivative and stays 2. The other amplitudes decay at rates κ and 9κ, giving u(x,t)=2+e−κtsin⁡x+e−9κtsin⁡(3x)u(x,t)=2+e^{-\kappa t}\sin x+e^{-9\kappa t}\sin(3x). As t→∞, both oscillations vanish and the profile approaches its spatial average, 2. Periodicity supplies the boundary conditions behind this mode basis.

Take it with you

  • Keep one transform convention throughout a calculation.
  • Frequency multipliers expose smoothing, inverse problems, and singular modes.
  • L2 convergence and pointwise reconstruction are different mathematical promises.

Read further

Original sources and deeper explanations to follow at your own pace.

Optional progress, saved in this browser when available. Revisit any lesson at any time.

Lesson not marked complete.