Properties of the Fourier Transform

The Fourier transform earns its central place in physics not as an isolated integral but as a dictionary. Every elementary manipulation of a function in the spatial or time domain — sliding it along its axis, multiplying it by an oscillation, stretching its argument, differentiating it, or convolving two functions — corresponds to a clean, usually algebraic, operation on the transform. Sliding becomes a phase; stretching becomes reciprocal stretching; differentiation becomes multiplication by ikik; convolution becomes an ordinary product. Assemble these correspondences and the transform stops being a formula to evaluate and becomes a calculus: you translate a hard analytic problem into frequency space, solve it there by inspection, and translate the answer back. Throughout we adopt the symmetric convention f^(k)=12π∫f e−ikxdx\hat f(k)=\frac{1}{\sqrt{2\pi}}\int f\,e^{-ikx}\mathrm{d}x, which splits the 2π2\pi evenly between transform and inverse and makes the transform unitary. Every property below carries 2π2\pi factors that depend on this choice, so we flag the convention at every turn — a convolution here picks up a 2π\sqrt{2\pi} that vanishes in the e−2πiνxe^{-2\pi i\nu x} convention. The natural home for these rules is the Schwartz space S\mathcal S of rapidly decreasing smooth functions, where the transform is a bijection; by density and duality they extend to L2L^2 and to tempered distributions. Three dualities organize everything. Translation and modulation are mirror images: a shift in one domain is a phase in the other. Differentiation and multiplication-by-coordinate are mirror images: ∂x↔ik\partial_x\leftrightarrow ik and −ix↔∂k-ix\leftrightarrow\partial_k. And scaling is reciprocal — the algebraic seed of the uncertainty principle. None of these is a coincidence; each follows from the symmetric way xx and kk enter the kernel e−ikxe^{-ikx}. Mastering this table is the gateway to the convolution theorem, to Parseval–Plancherel energy conservation, and to the transform method for differential equations.

Key takeaways
  • The transform converts operations in physical space into algebra in frequency space: translation becomes a linear phase e−ikae^{-ika}, modulation a rigid frequency shift, differentiation multiplication by ikik.
  • Scaling is reciprocal — squeeze a function in space and its transform spreads in frequency, the analytic root of the uncertainty relation Δx Δk≥12\Delta x\,\Delta k\ge \tfrac12.
  • Convolution maps to a pointwise product (the convolution theorem); in the symmetric convention f∗g^=2π f^g^\widehat{f*g}=\sqrt{2\pi}\,\hat f\hat g.
  • Duality exchanges function and transform: applying the transform twice is the parity operator, so F4=id\mathcal F^4=\mathrm{id}, and every property comes paired.
  • Symmetry ties reality and evenness of ff to Hermitian symmetry and reality of f^\hat f; the unit Gaussian e−x2/2e^{-x^2/2} is its own transform and saturates Δx Δk=12\Delta x\,\Delta k=\tfrac12 at every width.
  • Key idea: the Fourier transform intertwines pairs of operations — translation ↔ phase, modulation ↔ shift, differentiation ↔ multiplication by ikik, scaling ↔ reciprocal scaling, convolution ↔ product — and each rule has a dual obtained by swapping the two domains.
The core insight

The Fourier transform intertwines pairs of operations: translation ↔ phase, modulation ↔ shift, differentiation ↔ multiplication by ikik, scaling ↔ reciprocal scaling, and convolution ↔ product. Each rule has a dual obtained by swapping the spatial and frequency domains, because xx and kk enter the kernel e−ikxe^{-ikx} symmetrically.

Where this lives today

The property table is the daily working language of imaging, communication, and instrumentation. Magnetic-resonance imaging acquires its data in the frequency domain: the scanner samples kk-space and the image is the inverse transform, so the translation theorem is the tool radiologists use to register images, and the scaling theorem sets the field-of-view/resolution trade-off. Every 4G, 5G, and Wi-Fi link runs on orthogonal frequency-division multiplexing — the modulation theorem applied thousands of times per symbol, with the FFT doing the transform in silicon. Gravitational-wave observatories dig signals out of noise by matched filtering — convolution against a template — and the convolution theorem turns that expensive integral into a cheap multiplication. The shift theorem's deepest consequence remains a live open problem: in X-ray crystallography and coherent diffractive imaging you measure only the magnitude ∣f^(k)∣|\hat f(k)|, while all the structure sits in the phase you cannot record. Recovering that lost phase — the phase-retrieval problem — is "position lives in the phase" turned into a decades-long computational challenge, and why every protein structure solved by crystallography is a triumph over the missing phase.

The physical picture

Think of the kernel e−ikxe^{-ikx} as a meter reporting how much of frequency kk a signal contains. Slide the whole signal by aa and you have not changed which frequencies are present, only their relative timing — so the magnitude spectrum ∣f^∣|\hat f| is untouched and only a phase e−ikae^{-ika} appears. All information about where a feature sits is banished into the phase. The first figure makes this literal: a pulse and its shifted copy have identical magnitude spectra, and the entire record of the displacement is a straight-line phase ramp of slope −a-a.

Figure
Translation ↔ phase. Shifting a pulse by a distance a leaves the magnitude spectrum exactly unchanged and adds only the linear phase e^(−ika), a ramp of slope −a. All information about position lives in the phase, never in the magnitude.
Shift a pulse: the magnitude spectrum is unchanged, the phase picks up a ramp −ka position domain magnitude spectrum shift a = 2 f(x) f(x − a) |f̂(k)| identical for both The phase carries the shift arg = 0 arg f̂(x) = 0 arg f̂(x − a) = −ka angular frequency k →

Multiplying the signal by an oscillation eik0xe^{ik_0x} does the dual thing: it relabels every frequency by +k0+k_0, so the whole spectrum slides to f^(k−k0)\hat f(k-k_0) — the mirror of translation, and the reason AM radio, heterodyne receivers, and frequency combs all reduce to one theorem in one domain or the other. The deepest property is scaling: squeeze a pulse in time and it must contain higher frequencies over a wider band, since you cannot localize a wave sharply in both domains at once. The reciprocal law f(ax)^=∣a∣−1f^(k/a)\widehat{f(ax)}=|a|^{-1}\hat f(k/a) is the algebraic engine of the uncertainty relation, and the Gaussian is the one shape that trades off optimally, saturating the bound at every width.

Figure
Scaling is reciprocal. Compressing a Gaussian to width σ = 0.5 gives a narrow pulse Δx = 0.35 but a broad spectrum Δk = 1.41, while σ = 1.5 gives a wide pulse Δx = 1.06 and a narrow spectrum Δk = 0.47. In every case the product Δx·Δk = 1/2 — the Gaussian saturates the uncertainty bound.
Narrow in x ⇒ wide in k, and the product Δx·Δk = 1/2 for every width position domain frequency domain σ = 0.5, Δx = 0.35 σ = 1.5, Δx = 1.06 wide spectrum Δk = 1.41 narrow Δk = 0.47 Δx · Δk = 1/2 always

Finally, differentiation turns into multiplication by ikik for a reason at the heart of linear physics: ∂xeikx=ik eikx\partial_x e^{ikx}=ik\,e^{ikx}, so the exponentials are eigenfunctions of the derivative. Because every translation-invariant operator is built from derivatives (equivalently, is a convolution), the transform diagonalizes all of them at once — which is why the Fourier basis is the natural coordinate system for constant-coefficient differential equations, filters, and the normal modes of a translationally symmetric system.

The core equations

Everything follows from the defining pair, stated in the symmetric convention where 1/2π1/\sqrt{2\pi} appears both ways and the transform is unitary on L2L^2.

f^(k)=12π∫−∞∞f(x) e−ikx dx,f(x)=12π∫−∞∞f^(k) eikx dk(1)\hat f(k)=\frac{1}{\sqrt{2\pi}}\int_{-\infty}^{\infty} f(x)\,e^{-ikx}\,\mathrm{d}x,\qquad f(x)=\frac{1}{\sqrt{2\pi}}\int_{-\infty}^{\infty}\hat f(k)\,e^{ikx}\,\mathrm{d}k \qquad (1)

where: f(x)f(x) — function of the spatial or time variable xx; f^(k)\hat f(k) — its transform, a function of angular frequency kk; kk — angular frequency (radians per unit length or time); 1/2π1/\sqrt{2\pi} — symmetric normalization. Read f^(k)\hat f(k) as the overlap of ff with the probe eikxe^{ikx}, and ff as a superposition of those probes. Sliding ff by a real distance aa multiplies the transform by a pure phase; the magnitude is untouched.

f(x−a)^(k)=e−ika f^(k)(2)\widehat{f(x-a)}(k)=e^{-ika}\,\hat f(k) \qquad (2)

where: aa — real shift in the spatial domain; e−ikae^{-ika} — linear phase factor of unit modulus. Reading: a translation is invisible to ∣f^∣|\hat f|, appearing only as a phase slope. Its dual: multiplying ff by a carrier eik0xe^{ik_0x} translates the whole spectrum.

eik0xf(x)^(k)=f^(k−k0)(3)\widehat{e^{ik_0 x}f(x)}(k)=\hat f(k-k_0) \qquad (3)

where: k0k_0 — carrier angular frequency; f^(k−k0)\hat f(k-k_0) — the spectrum shifted by k0k_0. Reading: modulation in xx is the exact mirror of translation in kk. Stretching the argument by aa stretches the spectrum reciprocally and rescales its height.

f(ax)^(k)=1∣a∣ f^ ⁣(ka),a≠0(4)\widehat{f(ax)}(k)=\frac{1}{|a|}\,\hat f\!\left(\frac{k}{a}\right),\qquad a\neq 0 \qquad (4)

where: aa — real nonzero dilation of the argument; 1/∣a∣1/|a| — amplitude factor. Reading: narrowing in xx (a>1a>1) widens in kk — the uncertainty trade-off. Differentiation becomes multiplication by ikik; the nn-th derivative a degree-nn polynomial.

f′(x)^(k)=ik f^(k),f(n)^(k)=(ik)n f^(k)(5)\widehat{f'(x)}(k)=ik\,\hat f(k),\qquad \widehat{f^{(n)}}(k)=(ik)^n\,\hat f(k) \qquad (5)

where: f(n)f^{(n)} — the nn-th derivative, assumed to decay so boundary terms vanish; (ik)n(ik)^n — the polynomial multiplier replacing the operator. Reading: calculus in xx becomes arithmetic in kk. The dual rule trades multiplication by xx for differentiation in kk.

−ix f(x)^(k)=ddkf^(k)(6)\widehat{-ix\,f(x)}(k)=\frac{\mathrm{d}}{\mathrm{d}k}\hat f(k) \qquad (6)

where: −ix-ix — multiplication operator in xx; d/dk\mathrm{d}/\mathrm{d}k — differentiation in kk. Reading: the pair (x,∂x)(x,\partial_x) maps to (i∂k,ik)(i\partial_k,ik) — the Fourier fingerprint of the canonical commutator. Convolution and product are dual, with a convention-dependent 2π\sqrt{2\pi}.

f∗g^(k)=2π f^(k) g^(k),fg^(k)=12π (f^∗g^)(k)(7)\widehat{f*g}(k)=\sqrt{2\pi}\,\hat f(k)\,\hat g(k),\qquad \widehat{fg}(k)=\frac{1}{\sqrt{2\pi}}\,(\hat f*\hat g)(k) \qquad (7)

where: (f∗g)(x)=∫f(y)g(x−y) dy(f*g)(x)=\int f(y)g(x-y)\,\mathrm{d}y — the convolution; 2π\sqrt{2\pi} — artifact of the symmetric convention. Reading: smearing in one domain is a pointwise product in the other. Reality and evenness of ff constrain the transform.

f^(k)‾=f^(−k)  (f real);f^(−k)=f^(k)  (f even)(8)\overline{\hat f(k)}=\hat f(-k)\ \ (f\ \text{real});\qquad \hat f(-k)=\hat f(k)\ \ (f\ \text{even}) \qquad (8)

where:   ⋅  ‾\overline{\;\cdot\;} — complex conjugation. Reading: real ff has a Hermitian spectrum; real-and-even ff has a real-and-even spectrum, i.e. a cosine transform. The master duality: transforming twice reflects the argument.

F[f^](x)=f(−x) ⟹ F2=P,F4=id(9)\mathcal F[\hat f](x)=f(-x)\ \Longrightarrow\ \mathcal F^2=\mathcal P,\quad \mathcal F^4=\mathrm{id} \qquad (9)

where: F\mathcal F — the transform operator; P\mathcal P — the parity operator f(x)↦f(−x)f(x)\mapsto f(-x). Reading: the transform is a fourth root of the identity, so its eigenvalues are {1,i,−1,−i}\{1,i,-1,-i\} and its eigenbasis is the Hermite functions.

Derivation ADifferentiation becomes multiplication and why exponentials diagonalize every translation-invariant operator

The goal of this derivation is to prove the most consequential entry in the table, f′^=ik f^\widehat{f'}=ik\,\hat f, by one integration by parts, and then draw out the payoff most treatments leave implicit: because the exponentials are eigenfunctions of ∂x\partial_x, the Fourier transform is the change of basis that simultaneously diagonalizes every constant-coefficient (translation-invariant) linear operator, turning differential equations into algebra and filters into pointwise multiplication.

Step 1. Start from the definition applied to f′f', and integrate by parts with u=e−ikxu=e^{-ikx}, dv=f′ dx\mathrm{d}v=f'\,\mathrm{d}x: f′^(k)=12π∫f′(x)e−ikxdx=12π[fe−ikx]−∞∞+ik2π∫f(x)e−ikxdx\widehat{f'}(k)=\frac{1}{\sqrt{2\pi}}\int f'(x)e^{-ikx}\mathrm{d}x = \frac{1}{\sqrt{2\pi}}\big[f e^{-ikx}\big]_{-\infty}^{\infty}+\frac{ik}{\sqrt{2\pi}}\int f(x)e^{-ikx}\mathrm{d}x. Why: integration by parts is the one move that transfers the derivative off ff and onto the known kernel e−ikxe^{-ikx}, where differentiating just brings down a factor.

Step 2. For f∈Sf\in\mathcal S (indeed for any f→0f\to 0 at infinity) the boundary term vanishes, leaving f′^(k)=ik f^(k)\widehat{f'}(k)=ik\,\hat f(k). Why: the whole rule hinges on the surface term dying — for functions that do not decay, such as constants or plane waves, the naive rule fails and one must pass to tempered distributions.

Step 3. Iterate nn times: f(n)^(k)=(ik)nf^(k)\widehat{f^{(n)}}(k)=(ik)^n\hat f(k), so any operator L(∂x)=∑mcm∂xmL(\partial_x)=\sum_m c_m\partial_x^m transforms to multiplication by the polynomial L(ik)L(ik). Why: this converts a constant-coefficient linear ODE into a polynomial equation you can divide through — the subject of the third derivation below.

Step 4. The reason behind all of this is ∂xeikx=ik eikx\partial_x e^{ikx}=ik\,e^{ikx}: the exponentials are eigenfunctions of ∂x\partial_x. Since a translation-invariant operator commutes with all shifts, it is diagonalized in the same eigenbasis, and the transform is that diagonalizing map. Why: this is why the Fourier basis — not any other — is the privileged coordinate system for translationally symmetric physics, from crystal band structure to the modes of a uniform string. The rule f(n)^=(ik)nf^\widehat{f^{(n)}}=(ik)^n\hat f is verified against the defining integral, and its dual (6)(6) follows by differentiating under the integral sign.

Go deeper · Derivation B · 1 minScaling implies the uncertainty trade-off
Derivation BScaling implies the uncertainty trade-off

The goal here is to turn the reciprocal scaling law into a hard inequality — the Heisenberg–Weyl bound Δx Δk≥12\Delta x\,\Delta k\ge\tfrac12 — by combining the differentiation rule, Plancherel's identity, the Cauchy–Schwarz inequality, and one integration by parts, and to identify the Gaussian as the unique shape that saturates it. This is the analytic core from which the quantum position–momentum relation follows by the substitution k=p/ℏk=p/\hbar. Take ff real, normalized ∫∣f∣2dx=1\int|f|^2\mathrm{d}x=1, centered at the origin, with spreads Δx2=∫x2∣f∣2dx\Delta x^2=\int x^2|f|^2\mathrm{d}x and Δk2=∫k2∣f^∣2dk\Delta k^2=\int k^2|\hat f|^2\mathrm{d}k.

Step 1. By the differentiation rule f′^=ikf^\widehat{f'}=ik\hat f and Plancherel, ∫∣f′∣2dx=∫k2∣f^∣2dk=Δk2\int|f'|^2\mathrm{d}x=\int k^2|\hat f|^2\mathrm{d}k=\Delta k^2. Why: this rewrites the frequency spread as the energy of the derivative in xx-space, so both factors of the bound live in the same domain.

Step 2. Apply Cauchy–Schwarz to the pair xfxf and f′f': ∣∫xf f′‾ dx∣2≤(∫x2∣f∣2)(∫∣f′∣2)=Δx2 Δk2\big|\int x f\,\overline{f'}\,\mathrm{d}x\big|^2\le\big(\int x^2|f|^2\big)\big(\int|f'|^2\big)=\Delta x^2\,\Delta k^2. Why: Cauchy–Schwarz converts a product of two spreads into a bound on their overlap, with equality only when f′∝xff'\propto xf — whose solution is the Gaussian.

Step 3. Evaluate the real part by parts, using ddx∣f∣2=2 Re(f‾f′)\tfrac{\mathrm d}{\mathrm dx}|f|^2=2\,\mathrm{Re}(\overline f f'): Re∫xf f′‾ dx=12∫x ddx∣f∣2 dx=−12\mathrm{Re}\int x f\,\overline{f'}\,\mathrm{d}x=\tfrac12\int x\,\tfrac{\mathrm d}{\mathrm dx}|f|^2\,\mathrm{d}x=-\tfrac12. Why: the boundary term 12[x∣f∣2]\tfrac12[x|f|^2] vanishes for decaying ff, and normalization pins the remaining integral to 11.

Step 4. Since ∣⋅∣2≥(Re)2|\cdot|^2\ge(\mathrm{Re})^2, combine Steps 2 and 3: 14≤Δx2Δk2\tfrac14\le\Delta x^2\Delta k^2, hence Δx Δk≥12\Delta x\,\Delta k\ge\tfrac12. Why: equality requires both Cauchy–Schwarz saturation (f′∝xff'\propto xf) and reality, forcing the Gaussian f∝e−x2/2σ2f\propto e^{-x^2/2\sigma^2}, which hits the bound for every σ\sigma as the second figure shows. The constant 12\tfrac12, and the Gaussian saturation for σ=0.5,1.0,1.5\sigma=0.5,1.0,1.5 giving Δx Δk=0.5\Delta x\,\Delta k=0.5 in each case, were verified independently; with k=p/ℏk=p/\hbar this is exactly Δx Δp≥ℏ/2\Delta x\,\Delta p\ge\hbar/2.

Go deeper · Derivation C · 1 minSolving a screened Poisson equation by transform
Derivation CSolving a screened Poisson equation by transform

The goal of this derivation is to show the property table doing real work end to end: take the screened Poisson (Yukawa) equation −y′′+λ2y=g-y''+\lambda^2 y=g, use the differentiation rule to turn it into algebra in kk-space, invert with the convolution theorem while tracking every 2π\sqrt{2\pi} the symmetric convention demands, and recover the exponential Green's function — a result that must be convention-independent, and is. Assume y→0y\to0 at infinity and λ>0\lambda>0.

Step 1. Transform both sides and use linearity term by term: −y′′^+λ2y^=g^-\widehat{y''}+\lambda^2\hat y=\hat g. Why: linearity of the integral lets the transform pass through each term separately, which is the whole reason the method works on linear equations.

Step 2. Apply the differentiation rule twice, y′′^=(ik)2y^=−k2y^\widehat{y''}=(ik)^2\hat y=-k^2\hat y, giving (k2+λ2)y^=g^(k^2+\lambda^2)\hat y=\hat g and hence y^(k)=g^(k)/(k2+λ2)\hat y(k)=\hat g(k)/(k^2+\lambda^2). Why: the second derivative has become multiplication by −k2-k^2, so the differential equation is now a one-line division — no integration in kk-space at all.

Step 3. Invert with the convolution theorem. Writing h^(k)=1/(k2+λ2)\hat h(k)=1/(k^2+\lambda^2), the product h^g^\hat h\hat g equals 12πh∗g^\tfrac{1}{\sqrt{2\pi}}\widehat{h*g} from (7)(7), so y=12π(h∗g)y=\tfrac{1}{\sqrt{2\pi}}(h*g) with h=F−1[h^]h=\mathcal F^{-1}[\hat h]. Why: the symmetric 2π\sqrt{2\pi} must be carried here — dropping it, a common slip, corrupts the Green's-function prefactor.

Step 4. Evaluate hh by contour integration: ∫eikx/(k2+λ2) dk=(π/λ)e−λ∣x∣\int e^{ikx}/(k^2+\lambda^2)\,\mathrm{d}k=(\pi/\lambda)e^{-\lambda|x|}, so h(x)=12π⋅πλe−λ∣x∣=π2 1λe−λ∣x∣h(x)=\tfrac{1}{\sqrt{2\pi}}\cdot\tfrac{\pi}{\lambda}e^{-\lambda|x|}=\sqrt{\tfrac{\pi}{2}}\,\tfrac{1}{\lambda}e^{-\lambda|x|}. Why: this symmetric-convention inverse is not 12λe−λ∣x∣\tfrac{1}{2\lambda}e^{-\lambda|x|} — that form belongs to the physicist convention, and confusing the two is the classic source of a spurious 2π2\pi.

Step 5. Combine: y=12π⋅π2 1λ∫e−λ∣x−x′∣g(x′) dx′=12λ∫−∞∞e−λ∣x−x′∣g(x′) dx′y=\tfrac{1}{\sqrt{2\pi}}\cdot\sqrt{\tfrac{\pi}{2}}\,\tfrac{1}{\lambda}\int e^{-\lambda|x-x'|}g(x')\,\mathrm{d}x'=\dfrac{1}{2\lambda}\displaystyle\int_{-\infty}^{\infty} e^{-\lambda|x-x'|}g(x')\,\mathrm{d}x'. Why: the two 2π\sqrt{2\pi} factors cancel to leave 12λ\tfrac{1}{2\lambda} — the convolution of the source with the Yukawa Green's function G=12λe−λ∣x−x′∣G=\tfrac{1}{2\lambda}e^{-\lambda|x-x'|}, which must be convention-independent, and is. The pair 1/(k2+λ2)↔π/2 λ−1e−λ∣x∣1/(k^2+\lambda^2)\leftrightarrow \sqrt{\pi/2}\,\lambda^{-1}e^{-\lambda|x|} (symmetric convention) and the final 12λ\tfrac{1}{2\lambda} prefactor were both checked numerically.

Limiting cases

How the table specializes in the key regimes:

CaseResultMeaning
Pure shift, no scalingf(x−a)^(k)=e−ikaf^(k)\widehat{f(x-a)}(k)=e^{-ika}\hat f(k)The magnitude ∣f^∣|\hat f| is invariant under translation; every bit of position information lives in the phase.
Self-dual Gaussianf=e−x2/2⇒f^=e−k2/2f=e^{-x^2/2}\Rightarrow\hat f=e^{-k^2/2}In the symmetric convention the unit Gaussian is its own transform; scaling forces it to saturate Δx Δk=12\Delta x\,\Delta k=\tfrac12.
Real and evenf real, even⇒f^ real, evenf\ \text{real, even}\Rightarrow\hat f\ \text{real, even}Reality plus evenness collapse the symmetry rules and the transform reduces to a cosine transform; real-odd gives imaginary-odd.
Distributional limit1^(k)=2π δ(k)\widehat{1}(k)=\sqrt{2\pi}\,\delta(k)The rules extend to tempered distributions: an infinitely wide constant has an infinitely sharp spectrum, the extreme of the scaling law.
Narrow-pulse limitδ(x−a)^(k)=12πe−ika\widehat{\delta(x-a)}(k)=\tfrac{1}{\sqrt{2\pi}}e^{-ika}A point in space is a flat, unit-magnitude spectrum with a pure phase ramp — the shift theorem carried to a delta.

Worked examples

1 · A shifted, modulated Gaussian (confidence builder)

From the self-dual pair f(x)=e−x2/2↔f^(k)=e−k2/2f(x)=e^{-x^2/2}\leftrightarrow\hat f(k)=e^{-k^2/2}, find the transform of g(x)=ei5xf(x−3)g(x)=e^{i5x}f(x-3) and locate where its magnitude peaks.

Step 1. Handle the shift first with the translation theorem (2)(2): the transform of f(x−3)f(x-3) is e−3ikf^(k)=e−3ike−k2/2e^{-3ik}\hat f(k)=e^{-3ik}e^{-k^2/2} — magnitude untouched, only a phase ramp of slope −3-3.

Step 2. Now the modulation theorem (3)(3): multiplying by the carrier ei5xe^{i5x} shifts the whole spectrum by k0=5k_0=5, replacing kk by k−5k-5 everywhere. Thus g^(k)=e−3i(k−5) e−(k−5)2/2\hat g(k)=e^{-3i(k-5)}\,e^{-(k-5)^2/2}. Result: g^(k)=e−3i(k−5)e−(k−5)2/2\hat g(k)=e^{-3i(k-5)}e^{-(k-5)^2/2}.

Step 3. The magnitude ∣g^(k)∣=e−(k−5)2/2|\hat g(k)|=e^{-(k-5)^2/2} is a Gaussian peaking at k=5k=5 with value 11: modulation slid the spectral peak from 00 to 55, while the shift contributed only the phase e−3i(k−5)e^{-3i(k-5)}, which vanishes at the peak. The spectrum peaks at k=5k=5 with ∣g^(5)∣=1|\hat g(5)|=1; position lives in the phase, carrier frequency in the peak location.

2 · Rectangular pulse, the sinc, and scaling made numerical

A rectangular pulse f(x)=1f(x)=1 for ∣x∣≤T|x|\le T and 00 otherwise has transform f^(k)=2/π sin⁡(kT)/k\hat f(k)=\sqrt{2/\pi}\,\sin(kT)/k. For T=1T=1, find the zero-frequency value and the first spectral zero, then predict what happens when the pulse is halved to T=12T=\tfrac12.

Step 1. Zero-frequency value: as k→0k\to0, sin⁡(kT)/k→T\sin(kT)/k\to T, so f^(0)=2/π T\hat f(0)=\sqrt{2/\pi}\,T — the pulse area divided by 2π\sqrt{2\pi}. For T=1T=1 this is 2/π=0.798\sqrt{2/\pi}=\mathbf{0.798}. f^(0)=0.798\hat f(0)=0.798.

Step 2. First zero: sin⁡(kT)=0\sin(kT)=0 first at kT=πkT=\pi, so k=π/Tk=\pi/T, giving a main lobe of half-width π/T\pi/T. For T=1T=1 the spectrum first crosses zero at k=π=3.14k=\pi=\mathbf{3.14}. First zero at k=π≈3.14k=\pi\approx 3.14.

Step 3. Halve the pulse, T→12T\to\tfrac12: this is the scaling law (4)(4) with a=2a=2. The first zero moves to k=π/T=2π≈6.28k=\pi/T=2\pi\approx\mathbf{6.28}, twice as far out, and the amplitude f^(0)\hat f(0) halves to 0.3990.399. A pulse squeezed by two in space has a spectrum stretched by two and lowered by two — the reciprocal scaling of the uncertainty trade-off, made arithmetic. The Gaussian version of this trade-off is exact rather than approximate: for fσ=e−x2/2σ2f_\sigma=e^{-x^2/2\sigma^2} one finds Δx=σ/2\Delta x=\sigma/\sqrt2 and Δk=1/(σ2)\Delta k=1/(\sigma\sqrt2), so Δx Δk=12\Delta x\,\Delta k=\tfrac12 at every width (as the second figure confirms for σ=0.5,1.0,1.5\sigma=0.5,1.0,1.5). Halving the pulse doubles the spectral width to k=2πk=2\pi and halves the peak; the Gaussian ties the uncertainty bound at every width.

Physics of the portrait

The property table is a portrait of one man's idea outrunning its own rigor. Fourier asserted, to an unconvinced Paris Academy, that any function could be built from sines and cosines — a claim so bold his examiners refused it publication for years. The transform and its algebra are the mature form of that assertion, and it took another century of Dirichlet, Riemann, Plancherel, and finally Schwartz to put every property on a rigorous foundation.

Origins · 1 minJoseph Fourier and how this was discovered
JF
Joseph Fourier
1768–1830 · heat, trigonometric series, and the algebra of the transform

Jean-Baptiste Joseph Fourier was an orphan, a revolutionary who narrowly escaped the guillotine, and the scientific governor of Napoleon's occupied Egypt before he was the mathematician we name the transform for. His decisive work grew out of a concrete problem — how heat diffuses through a solid — and in his 1807 memoir to the Institut de France he made the audacious claim that an arbitrary function could be represented as a sum of sines and cosines. Lagrange, Laplace, and Legendre balked: that a single formula could reproduce a curve with corners and jumps seemed absurd, and the memoir went unpublished. Fourier pressed on, and his 1822 Théorie analytique de la chaleur laid out both the series for periodic functions and the integral transform on the whole line, complete with the operational rules — how the transform responds to differentiation and translation — that turned it into a method for the heat equation. He was right about far more than he could prove; the convergence questions his claim raised drove a century of analysis and, through Dirichlet, Riemann, and Lebesgue, essentially created the modern theory of functions. The uncertainty principle latent in his scaling rule, the energy identity Plancherel proved unitary in 1910, and the extension to the Dirac delta that Schwartz's distributions legitimized in the 1940s all descend from the one insight Fourier saw first: that the language of frequency turns the calculus of heat, waves, and diffusion into algebra.

History
1807
— Fourier submits his memoir on the propagation of heat to the Institut de France, claiming any function expands in sines and cosines
1822
— Théorie analytique de la chaleur lays out the series, the integral transform, and the operational rules
1910
— Plancherel proves the L² energy identity, making the transform a unitary map
1933
— Wiener's generalized harmonic analysis extends the transform to signals of infinite energy
1945
— Schwartz's theory of distributions extends every property to the Dirac delta and its derivatives

Going further

Texts. Stone & Goldbart, Mathematics for Physics, §2.2 · Arfken, Weber & Harris, Mathematical Methods for Physicists, 7th ed., Ch. 20 · Hassani, Mathematical Physics, 2nd ed., Ch. 8 · Byron & Fuller, Mathematics of Classical and Quantum Physics, Ch. 5 · Riley, Hobson & Bence, Mathematical Methods for Physics and Engineering, 3rd ed., §13.1 · Folland, Fourier Analysis and Its Applications, Ch. 7 (rigorous, on Schwartz space).

Provenance. Conventions and property formulas follow Stone & Goldbart §2.2 (symmetric 1/2π1/\sqrt{2\pi}). A verifier should confirm: the self-dual pair e−x2/2↔e−k2/2e^{-x^2/2}\leftrightarrow e^{-k^2/2} and the symmetric-convention pair 1/(k2+λ2)↔π/2 λ−1e−λ∣x∣1/(k^2+\lambda^2)\leftrightarrow\sqrt{\pi/2}\,\lambda^{-1}e^{-\lambda|x|} (not 12λe−λ∣x∣\tfrac{1}{2\lambda}e^{-\lambda|x|}, which is the physicist-convention form); the uncertainty constant 12\tfrac12 with Gaussian saturation Δx Δk=0.500\Delta x\,\Delta k=0.500 for σ=0.5,1.0,1.5\sigma=0.5,1.0,1.5; the rectangular-pulse values f^(0)=2/π=0.798\hat f(0)=\sqrt{2/\pi}=0.798 and first zero at k=πk=\pi; and F4=id\mathcal F^4=\mathrm{id}. All 2π2\pi factors are convention-dependent and stated explicitly.

Prerequisites. The Fourier Transform · Fourier Series and Periodic Functions · The Dirac Delta Function. Leads to. The Convolution Theorem · Parseval's and Plancherel's Theorems · Solving Differential Equations by Transform Methods · Tempered Distributions and Fourier Transforms.

Ask the tutor: Ask me to run any signal through the property table step by step — shift it, scale it, modulate it, or feed it to a differential equation — and I will show you the transform emerge one rule at a time.

What does "duality" actually mean, and why is F⁴ the identity?

Duality says the forward and inverse transforms differ only by the sign in the exponent, so transforming twice reproduces the inversion integral at −x-x: F2f=f(−x)\mathcal F^2 f = f(-x), the parity operator. Parity twice is the identity, so F4=id\mathcal F^4=\mathrm{id}. The consequence is beautiful: since F4=1\mathcal F^4=1, its eigenvalues can only be the fourth roots of unity {1,i,−1,−i}\{1,i,-1,-i\}, and the eigenfunctions are the Hermite functions Hn(x)e−x2/2H_n(x)e^{-x^2/2} with eigenvalue (−i)n(-i)^n. The Gaussian is the n=0n=0 case, eigenvalue 11, which is why it is self-dual.

Why does differentiation turn into multiplication by ik?

Because the complex exponentials eikxe^{ikx} are eigenfunctions of the derivative operator: ∂xeikx=ik eikx\partial_x e^{ikx}=ik\,e^{ikx}. The transform decomposes ff into these exponentials, and differentiating ff just multiplies each coefficient by ikik. Concretely, one integration by parts moves the derivative off ff and onto the kernel e−ikxe^{-ikx}, with the boundary term vanishing for any decaying function. The payoff is enormous: every constant-coefficient differential operator becomes a polynomial in kk, so a differential equation becomes algebra you solve by division — the mechanism behind the transform method for the heat, wave, and Helmholtz equations.

Is this uncertainty relation the same as Heisenberg's in quantum mechanics?

It is the same inequality, and the quantum version is a special case. The mathematical statement Δx Δk≥12\Delta x\,\Delta k\ge\tfrac12 holds for any function and its transform — a theorem about waves, with no physics beyond the reciprocal scaling of the transform. Quantum mechanics adds one ingredient: de Broglie's relation p=ℏkp=\hbar k. Substituting k=p/ℏk=p/\hbar gives Δx Δp≥ℏ/2\Delta x\,\Delta p\ge\hbar/2, the Heisenberg uncertainty principle. So the physical content — that you cannot simultaneously localize position and momentum — is inherited entirely from the mathematics of Fourier pairs, with ℏ\hbar supplying only the unit conversion.

Why is the magnitude spectrum unchanged when I shift a signal?

Because shifting a signal changes only the relative timing of its frequency components, not how much of each is present. The transform of a translate is the original times a pure phase e−ikae^{-ika}, and ∣e−ika∣=1|e^{-ika}|=1, so ∣f^(k)∣|\hat f(k)| is identical for a signal and all its shifts. The record of where the signal sits is stored entirely in the phase — the linear ramp of slope −a-a in the first figure. This is a deep fact with practical teeth: in crystallography and coherent imaging you measure the magnitude but not the phase, and because position lives in the phase, recovering an image from magnitude alone (the phase-retrieval problem) is genuinely hard.

Why does a product become a convolution and vice versa?

Because the transform diagonalizes convolution. A convolution f∗gf*g is a translation-invariant smearing, and translation-invariant operations are exactly the ones the exponentials diagonalize — so in the frequency domain, where those exponentials are the coordinates, convolution collapses to pointwise multiplication (times a 2π\sqrt{2\pi} in the symmetric convention). Read the other way, multiplying two signals mixes their frequencies, and that mixing is a convolution of their spectra. This makes the transform indispensable computationally: an O(N2)O(N^2) convolution becomes an O(N)O(N) multiply, and the FFT gets you there in O(Nlog⁡N)O(N\log N) — why matched filters, DSP, and large-integer multiplication all run through the Fourier domain.

When do these rules fail, and what fixes them?

The rules assume enough decay for the boundary terms in integration by parts to vanish and the integrals to converge — the Schwartz space S\mathcal S is the clean setting where every property holds exactly. They fail for functions that do not decay: a constant, a plane wave, or the delta have no ordinary Fourier integral. The fix is the theory of tempered distributions, which extends the transform to these objects by duality against Schwartz test functions. There the differentiation and translation rules survive verbatim, and you get the crucial pairs 1^=2π δ\widehat{1}=\sqrt{2\pi}\,\delta and δ^=1/2π\widehat{\delta}=1/\sqrt{2\pi} — the content of the follow-on article on tempered distributions.

Before this

Where this leads