The Fourier series representation of a signal, as expressed for
a square wave by
st=
a
0
+∑
k
=1∞
a
k
cos2πktT+∑
k
=1∞
b
k
sin2πktT
s
t
a
0
k
1
a
k
2
k
t
T
k
1
b
k
2
k
t
T
(1)
says that the left and right sides are "equal". We need to
investigate equality through an example.
Let's find the spectrum of the square wave
sqt
sq
t
.
The expressions for the Fourier coefficients have the common
form
a
k
b
k
=2T∫0T2cos2πktTsin2πktTd
t
−2T∫T2Tcos2πktTsin2πktTd
t
a
k
b
k
2
T
t
0
T
2
2
k
t
T
2
k
t
T
2
T
t
T
2
T
2
k
t
T
2
k
t
T
(2)
The cosine coefficients
a
k
a
k
are all zero, and the sine coefficients are
b
k
={4πk if k is odd0 if k is even
b
k
4
k
k is odd
0
k is even
(3)
Thus, the Fourier series for the square wave is
sqt=∑k∈13…4πksin2πktT
sq
t
k
k
1
3
…
4
k
2
k
t
T
(4)
As we see in
Figure 1, the
Fourier series requires many more terms to provide the same
quality of approximation as we found with the
half-wave
rectified sinusoid.
We can verify that more terms are needed by considering the
power spectrum and the approximation error shown in
Figure 2.
This difference between the two Fourier series results because
the half-wave rectified sinusoid's Fourier coefficients are
proportional to
1k2
1
k
2
while those of the square wave are proportional to
1k
1
k
.
In short, the square wave's coefficients decay more slowly
with increasing frequency. Said another way, the square-wave's
spectrum contains more power at higher frequencies than does
the half-wave-rectified sinusoid.
Calculate the harmonic distortion for the square wave.
Total harmonic distortion in the square wave is
1−4π2=20%
1
4
2
20
%
.
When comparing the square wave to its Fourier series
representation it is not clear that the two are equal. The fact
that the square wave's Fourier series requires more terms for a
given representation accuracy is not important. However, close
inspection of Figure 3 does reveal a
potential issue: Does the Fourier series really equal the square
wave at all values of
t
t?
In particular, at each step-change in the square wave, the
Fourier series exhibits a peak followed by rapid
oscillations. As more terms are added to the series, the
oscillations seem to become more rapid and smaller, but the
peaks are not decreasing. Consider this mathematical question
intuitively: Can a discontinuous function, like the square wave,
be expressed as a sum, even an infinite one, of continuous ones?
One should at least be suspicious, and in fact, it can't be thus
expressed. This issue brought
Fourier
much criticism from the French Academy of Science (Laplace,
Legendre, and Lagrange comprised the review committee) for
several years after its presentation on 1807. It was not
resolved for also a century, and its resolution is interesting
and important to understand from a practical viewpoint.
The extraneous peaks in the square wave's Fourier series
never disappear; they are termed
Gibb's phenomenon after the American physicist
Josiah Willard Gibbs. They occur whenever the signal is
discontinuous, and will always be present whenever the signal
has jumps. Let's return to the question of equality; how can the
equal sign in the
definition of
the Fourier series
be justified? The partial answer is that pointwise—each and
every value of
t
t—equality
is not guaranteed. What mathematicians
later in the nineteenth century showed was that the rms error of
the Fourier series was always zero.
limit
K
→
∞
rms
ε
K
=0
K
rms
ε
K
0
(5)
What this means is that the difference between an actual signal
and its Fourier series representation may not be zero, but the
square of this quantity has
zero integral!
It is through the eyes of the rms value that we define equality:
Two signals
s
1
t
s
1
t
,
s
2
t
s
2
t
are said to be equal in the
mean square if
rms
s
1
−
s
2
=0
rms
s
1
s
2
0
.
These signals are said to be equal
pointwise if
s
1
t=
s
2
t
s
1
t
s
2
t
for all values of
tt. For Fourier
series, Gibb's phenomenon peaks have finite height and zero
width: The error differs from zero only at isolated
points--whenever the periodic signal contains
discontinuities--and equals about 9% of the size of the
discontinuity. The value of a function at a finite set of points
does not affect its integral. This effect underlies the reason
why defining the value of a discontinuous function, like we
refrained from doing in defining the
step function, at its
discontinuity is meaningless. Whatever you pick for a value has
no practical relevance for either the signal's spectrum or for
how a system responds to the signal. The Fourier series value
"at" the discontinuity is the average of the values on either
side of the jump.