Fourier Series Simulator

The core idea is that any periodic function with period $T$ can be constructed from a constant plus a sum of sines and cosines at harmonically related frequencies.

Here, $f(t)$ is the synthesized waveform, $N$ is the number of harmonics used (the main simulator control), $T$ is the period, and $a_n$ and $b_n$ are the Fourier coefficients that depend on the target waveform shape.

The coefficients are calculated from the original function $g(t)$ using integration. This is how the simulator "knows" the recipe for a square, triangle, or sawtooth wave.

For a square wave of amplitude $A$, this yields $b_n = \frac{4A}{n\pi}$ for odd $n$ (others zero). For a triangle wave, $a_n$ coefficients follow a $\frac{1}{n^2}$ pattern. This difference explains why the triangle wave converges much faster—its higher harmonic amplitudes decay rapidly.

Audio Synthesis & Music: Fourier synthesis is the foundation of subtractive synthesis in electronic music. A sawtooth wave, rich in harmonics, is generated and then filtered to create the sounds of synthesizers. The simulator directly shows the harmonic content that defines a sound's timbre.

Signal Processing & Compression: JPEG image and MP3 audio compression rely on a related technique (the Fourier transform). By identifying and keeping only the most significant frequency components (like the strong lower harmonics you see in the spectrum), file sizes can be drastically reduced with minimal perceived quality loss.

Vibration Analysis in Mechanical Engineering: In CAE, irregular vibrations from machinery are decomposed into their harmonic components via Fourier analysis. This allows engineers to identify the source of a vibration—for instance, pinpointing whether it's caused by an unbalanced fan at 30 Hz or bearing defect at 150 Hz—enabling targeted fixes.

Solving Partial Differential Equations: Fourier's original application was solving the heat equation. The temperature distribution in a metal rod can be expressed as a sum of sinusoidal modes. This "separation of variables" method is still a cornerstone for analytical solutions in heat transfer, acoustics, and quantum mechanics.

First, you might think that "increasing N always improves accuracy," but in reality, it's not that simple. While this holds true for smooth waveforms like the triangle wave, for waveforms with discontinuities like the square wave, the approximately 9% overshoot due to the Gibbs phenomenon does not disappear, even if N is increased to 100 or 1000. This is a mathematically proven convergence property. In practice, "how many terms to include" is decided by balancing acceptable error with computational cost. For instance, if a 9% distortion is critical in audio signal reproduction, approximating with a square wave itself might be inappropriate.

Next, note that the "spectrum" you see in the simulator is that of an ideal periodic waveform. The spectrum obtained by applying FFT to actual CAE data (e.g., time-history vibration data) contains noise and non-stationary components and will not show such clean discrete lines. For example, performing FFT on engine vibration data will show peaks corresponding to rotational speed, but there will always be broad noise in the background. It's crucial not to equate the clean spectrum from the tool directly with reality.

Finally, understand the "waveform properties" indicated by the decay rates of the coefficients, such as $1/n$ or $1/n^2$ . A slower decay (like $1/n$ for a square wave) means more high-frequency components and includes sharp changes. Conversely, a waveform with faster decay (like $1/n^2$ for a triangle wave) is smoother. This directly relates to practical decisions like mesh sizing or time-step settings in simulations. To simulate sharp phenomena like shock waves, where fast decay cannot be assumed, very fine meshes or high sampling frequencies become necessary.