Fourier Series & Epicycles Wave Synthesizer — 3D Lab Simulation
Decomposition of periodic functions into harmonic sinusoidal oscillations
Theoretical Foundation & Historical Context
In 1822, Joseph Fourier proved that any arbitrary periodic function can be decomposed into an infinite summation of sinusoidal and cosinusoidal harmonics. By chaining rotating circles (epicycles) of calculated radii and frequencies, the tip synthesizes the target wave in real-time.
Law: Fourier Expansion Theorem & Orthogonality of Trigonometric Systems | Scientist: Jean-Baptiste Joseph Fourier | Year: 1822
Key Learning Objectives & Competencies
- ✓ Decompose arbitrary periodic waveforms into an infinite convergent sum of orthogonal sine and cosine harmonics.
- ✓ Calculate Fourier coefficients a_n and b_n using orthogonality integration over one period.
- ✓ Demonstrate the Gibbs phenomenon: 8.95% persistent overshoot at step discontinuities as N approaches infinity.
Mathematical Formulations & Differential Equations
f(t) = a₀/2 + ∑ [ aₙ·cos(nωt) + bₙ·sin(nωt) ]
Infinite sum of harmonics with fundamental frequency ω = 2π / T
f_sq(t) = (4/π) · ∑ [ sin((2k-1)ωt) / (2k-1) ]
Contains exclusively odd harmonic components (n = 1, 3, 5, ...)
lim |f_N(t₀⁺) - f(t₀⁺)| ≈ 0.089490 · Δf
Persistent ~8.95% overshoot at discontinuities regardless of N
Physical Parameters & Units of Measurement
| Symbol | Physical Quantity | SI Unit | Default |
|---|---|---|---|
| N | Number of harmonic terms synthesized | — | — |
| ω | Fundamental angular frequency (rad/s) | — | — |
| bₙ | Amplitude of n-th harmonic corresponding to epicycle radius | — | — |
Real-World & Industrial Engineering Applications
Digital audio compression (MP3, AAC) and discrete cosine transform image compression (JPEG).
Digital Signal Processing (DSP) for bandpass filtering, acoustic noise cancellation, and equalization.
Quantum mechanics wave packet decomposition and Fourier-transform infrared spectroscopy (FTIR).
3-Step Controlled Experimental Protocol
- Step 1: Baseline Setup — Initialize the simulation in standard conditions and observe undisturbed system dynamics.
- Step 2: Parameter Variation — Gradually modulate physical parameters using the control panel sliders.
- Step 3: Verification & Phase Space Analysis — Compare phase space trajectories and conservation metrics against analytical ODE solutions.
Frequently Asked Questions (FAQ)
What types of functions can be represented by a Fourier series?
Any periodic function satisfying Dirichlet's conditions: it must be integrable over one period, have a finite number of local extrema, and have a finite number of isolated discontinuities.
Why does the Gibbs phenomenon occur at square wave edges?
Because continuous sine and cosine functions cannot converge uniformly at a jump discontinuity. The Fourier approximation always overshoots the jump by approximately 8.95% on both sides regardless of how many harmonics are added.
What is the difference between a Fourier Series and a Fourier Transform?
A Fourier Series applies strictly to periodic signals, decomposing them into discrete harmonic frequencies. The Fourier Transform generalizes this to non-periodic continuous signals, producing a continuous frequency spectrum.
Related Simulations & Experiments
Nonlinear dynamical system and deterministic chaos theory
Hypotenuse square area conservation & Generalized Law of Cosines
Calculus of variations: Fastest descent trajectory is a cycloid, not a straight line
Interactive Simulation: Interactive real-time 60 FPS WebGL simulation is active in your browser. Rotate, zoom, and inspect physical apparatus components with full camera freedom.