The harmonic analysis of continuous waveforms

Written by Guy Fernando

Created Jan 2020 - Last modified Sep 2024

Jean-Baptiste Joseph Fourier (c. 1768 - c. 1830) was a French mathematician and physicist noted for his investigations into
heat flow, and scientific advisor to Napoleon. He is also recognized for discovering that any periodic motion can be described
as a infinite series of trigonometric functions, now commonly known as the Fourier series. He is lesser known for the discovery
of the greenhouse effect, having determined that the temperature on Earth would in fact be lower if not for the insulating
properties of the atmosphere. Believing that heat also held healing properties, he met his fate by tripping down the stairs in
his house possibly while being wrapped up in a blanket.

In 1807, Fourier submitted his paper to the Institut de France, claiming that any periodic signal could be represented by a series. The principle idea of a Fourier series expansion is to represent a function f(x) of period 2π as an infinite series of sinusoidal and cosinusoidal functions, given by the following equations.

\( \begin{aligned}
f(x) = \frac{a_{0}}{2} + a_{1}.cos \; x + a_{2}.cos \; 2x + a_{3}.cos \; 3x \; + ... +
\; b_{1}.sin \; x + b_{2}.sin \; 2x + b_{3}.sin \; 3x ...
\end{aligned} \)

\( \begin{aligned}
f(x) = \frac{a_{0}}{2} + \sum_{n=1}^ \infty (a_{n}.cos(nx) + b_{n}.sin(nx))
\;\;\;\;\;\;(1)
\end{aligned} \)

Where integrals are used to find the Fourier series coefficients, given by the following equations.

\( \begin{aligned}\;\;\;\;\;
a_{0} = \frac{1}{\pi } \int_{0}^{2\pi} f(x) \; dx
\;\;\;\;\;\;(2)
\end{aligned} \)

\( \begin{aligned}\;\;\;\;\;
a_{n} = \frac{1}{\pi } \int_{0}^{2\pi} f(x).cos(nx) \; dx
\;\;\;\;\;\;(3)
\end{aligned} \)

\( \begin{aligned}\;\;\;\;\;
b_{n} = \frac{1}{\pi } \int_{0}^{2\pi} f(x).sin(nx) \; dx
\;\;\;\;\;\;(4)
\end{aligned} \)

Note, all waveforms presented in this article are symmetric about the x-axis, and so will have no DC component.

\( \begin{aligned}
\therefore a_{0} = 0
\end{aligned} \)

Also note, all waveforms are asymmetric about the y-axis (odd functions), and so will have no cosine terms.

\( \begin{aligned}
\therefore a_{n} = 0
\end{aligned} \)

A sine wave or sinusoid is a mathematical curve that describes a smooth periodic oscillation.

By definition a sine wave is the first and the only term in a Fourier series. It having only the first harmonic that being
the fundamental tone.

\( \begin{aligned}
b_{1} = 1
\end{aligned} \)

\( \begin{aligned}
f(x) = A.sin(x)
\end{aligned} \)

A square wave is a non-sinusoidal periodic waveform in which the amplitude alternates at a steady frequency between
fixed minimum and maximum values, with the same duration at minimum and maximum.

A single period of the square wave between 0 and 2π can be broken down into two discrete intervals, and the integral
of each interval taken.

Using equation (4).

\( \begin{aligned}
b_{n} = \frac{1}{\pi} \int_0^\pi A.sin(nx) \; dx \; + \; \frac{1}{ \pi } \int_\pi^{2\pi} -A.sin(nx) \; dx
\end{aligned} \)

\( \begin{aligned}
\;\;\; = \frac{A}{\pi} \big\{ \left[\frac{-cos(nx)}{n}\right]_0^\pi - \left[\frac{-cos(nx)}{n}\right]_\pi^{2\pi}\big\}
\end{aligned} \)

\( \begin{aligned}
\;\;\; = \frac{A}{n\pi} \big\{ \left[-cos(n\pi)-(cos(0))\right] - \left[-cos(2n\pi)-(-cos(n\pi))\right] \big\}
\end{aligned} \)

\( \begin{aligned}
\;\;\; = \frac{A}{n\pi} \big\{ \left[1-cos(n\pi) \right] - \left[cos(n\pi)-1 \right] \big\}
\end{aligned} \)

\( \begin{aligned}
b_{n} = \frac{4A}{\pi} \left[\frac{1-cos(n\pi)}{2n}\right]
\end{aligned} \)

Substituting values of *b _{n}* back into equation (1).

\( \begin{aligned} f(x) = \frac{4A}{ \pi } \left[ sin(x) + \frac{1}{3}sin(3x) + \frac{1}{5}sin(5x) + \frac{1}{7}sin(7x) \; ... \right] \end{aligned} \)

The equation shows that a square wave contains only odd harmonics, with subsequent overtones decreasing in amplitude.

The sawtooth wave is a non-sinusoidal waveform. It is so named based on its resemblance to the teeth of a plain-toothed saw
with a zero rake angle. The convention is that a sawtooth wave ramps upward and then sharply drops.

A single period of the sawtooth wave between 0 and 2π can be broken down into two discrete intervals, and the integral
of each interval taken.

Using equation (4).

\( \begin{aligned}
b_{n} = \frac{1}{\pi} \int_0^\pi \big(\frac{Ax}{\pi}\big) .sin(nx) \; dx \; + \;
\frac{1}{ \pi } \int_\pi^{2\pi} \big(\frac{Ax}{\pi}-2A\big) .sin(nx) \; dx
\end{aligned} \)

Using integration by parts.

\( \begin{aligned}
\;\;\; = \frac{A}{\pi} \big\{ \left[\frac{sin(nx)-nx.cos(nx)}{ {\pi}n^{2} }\right]_0^\pi +
\left[\frac{sin(nx)-n(x-2\pi).cos(nx)}{{\pi}n^{2}}\right]_\pi^{2\pi}\big\}
\end{aligned} \)

After simplifying, for values of n > 0.

\( \begin{aligned}
b_{n} = \frac{2A}{\pi} \left[\frac{-(-1)^{n} }{n}\right]
\end{aligned} \)

Substituting values of *b _{n}* back into equation (1).

\( \begin{aligned} f(x) = \frac{2A}{ \pi } \left[ sin(x) - \frac{1}{2}sin(2x) + \frac{1}{3}sin(3x) - \frac{1}{4}sin(4x) \; ... \right] \end{aligned} \)

The equation shows that a sawtooth wave contains all harmonics, with subsequent overtones decreasing in amplitude.

A triangle wave is a non-sinusoidal waveform named for its triangular shape. It is a periodic, piecewise linear, continuous
real function.

A single period of the triangle wave between 0 and 2π can be broken down into three discrete intervals, and the integral
of each interval taken.

Using equation (4).

\( \begin{aligned}
b_{n} = \frac{1}{\pi}\int_0^{\frac{\pi}{2}}\big(\frac{2Ax}{\pi}\big).sin(nx)\;dx \; + \;
\frac{1}{\pi}\int_{\frac{\pi}{2}}^{\frac{3\pi}{2}}\big(\frac{-2Ax}{\pi}+2A\big).sin(nx)\;dx \; + \;
\frac{1}{\pi}\int_{\frac{3\pi}{2}}^{2\pi}\big(\frac{2Ax}{\pi}-4A\big).sin(nx)\;dx
\end{aligned} \)

Using integration by parts.

\( \begin{aligned}
\;\;\; = \frac{A}{\pi}\big\{\left[\frac{2(sin(nx)-nx.cos(nx))}{ {\pi}n^{2} }\right]_0^{\frac{\pi}{2}} +
\left[\frac{-2(sin(nx)+({\pi}x-nx).cos(nx))}{{\pi}n^{2}}\right]_{ \frac{\pi}{2} }^{\frac{3\pi}{2}} +
\left[\frac{2(sin(nx)+(2{\pi}n-nx).cos(nx))}{{\pi}n^{2}}\right]_{\frac{3\pi}{2}}^{2\pi}\big\}
\end{aligned} \)

After simplifying, for values of n > 0.

\( \begin{aligned}
b_{n} = \frac{8A}{\pi^2}\left[\frac{(-(-1)^{n}(1-cos(n\pi)))}{2n^2}\right]
\end{aligned} \)

Substituting values of *b _{n}* back into equation (1).

\( \begin{aligned} f(x) = \frac{8A}{\pi^{2}}\left[ sin(x) - \frac{1}{9}sin(3x) + \frac{1}{25}sin(5x) - \frac{1}{49}sin(7x) \; ... \right] \end{aligned} \)

The equation shows that a triangle wave contains only odd harmonics, with subsequent overtones decreasing in amplitude.

Now that we have calculated the coefficients of each Fourier series term, we can plot the series using javascript to see how close it resembles the ideal waveform it represents. Recall that to generate square, sawtooth or triangle waveforms, an infinite number of terms is required, however below, only the first four terms are plotted.

Notice that the plotted waveform resembles the ideal waveform, but is far from perfect. You can add additional terms by typing
them into the "series" textbox above. This will illustrate that when more terms are included, the closer the waveform
approximation approaches the ideal. Change the "waveform" dropdown list above to show different waveform approximations.

The sine, square, sawtooth and triangle waveforms produced by common audio signal generators is described here, however the
Fourier series for any periodic waveform may be calculated using this method. See this example of an online
Audio Signal Generator here on i4cy.com.

A signal synthesised by a signal generator at its upper frequency limit is often imperfect, and appears similar to a Fourier
series missing its higher harmonic terms. This in reality is mainly due to the limited frequency response of analogue electronic
circuitry, and in the case of digital electronic synthesis, the digital to analogue converter's Nyquist limit.

Widely available computing tools make it far easier to learn about the Fourier series compared to when I was taught the subject
as a student. Even though relatively involved calculus is required, for example when solving expressions that require "integration
by parts", free tools such WolframAlpha or
Symbolab can be used to evaluate such steps without knowledge of advanced
mathematics.

This website is powered using ultra low power green locally based servers.

Copyright © i4cy 2000-2024. All rights reserved.