
P1: OSO/OVY P2: OSO/OVY QC: OSO/OVY T1: OSO
MHDQ256-Ch09 MHDQ256-Smith-v1.cls December 17, 2010 19:47
LT (Late Transcendental)
CONFIRMING PAGES
620 CHAPTER 9
..
Infinite Series 9-90
Knowing the signal is of the form A sin Bt, you would use the
data to try to solve for A and B. In this case, you don’t have
enough information to guarantee getting the right values for
A and B. Prove this by finding several values of A and B with
B = 8π that match the data. A famous result of H. Nyquist
from 1928 states that to reconstruct asignal of frequency f you
need at least 2 f samples.
45. The energy of a signal f (x) on the interval [−π, π]is
defined by E =
1
π
!
π
−π
[ f (x)]
2
dx.If f (x) has a Fourier
series f (x) =
a
0
2
+
∞
k=0
[a
k
coskx + b
k
sinkx], show that
E = A
2
0
+ A
2
1
+ A
2
2
+···, where A
k
=
a
2
k
+ b
2
k
. The se-
quence {A
k
} is called the energy spectrum of f (x).
46. Carefully examine the graphs in Figure 9.46. There is a Gibbs
phenomenon at x = 0. Does it appear that the size of the
Gibbsovershootchanges as thenumberof terms increases? We
examine that here. For the partial sum F
n
(x) as defined in
example9.1,itcanbe shownthattheabsolutemaximumoccurs
at
π
2n
. Evaluate F
n
π
2n
for n = 4, n = 6 and n = 8. Show
thatforlargen, the size of thebumpis
F
n
π
2n
− f
π
2n
≈
0.09. Gibbs showed that, in general, the size of the bump at a
jump discontinuity is about 0.09 times the size of the jump.
47. Some fixes have been devised to reduce the Gibbs phe-
nomenon. Define the σ-factors by σ
k
=
sin
kπ
n
kπ
n
for
k = 1, 2,...,n and consider the modified Fourier sum
a
0
2
+
n
k=0
[a
k
σ
k
coskx + b
k
σ
k
sinkx]. For example 9.1, plot the
modified sums for n = 4 and n = 8 and compare to Figure
9.46: f (x) =
1, −π<x < 0
−1, 0 < x <π
, F
2n−1
has critical point
at π/2n and lim
n→∞
F
2n−1
π
2n
=
2
π
π
0
sin x
x
dx ≈ 1.18.
EXPLORATORY EXERCISES
1. Suppose that you wanted to approximatea waveform with sine
functions (no cosines), as in the music synthesizer problem.
Such a Fourier sine series will be derivedin this exercise.You
essentially use Fourier series with a trick to guarantee sine
terms only. Start with your waveform as a function defined on
theinterval[0, l],forsomelengthl.Thendefineafunctiong(x)
that equals f (x)on[0, l] and thatis an odd function. Showthat
g(x) =
f (x), if 0 ≤ x ≤ l
− f (−x), if −l < x < 0
works. Explain why the
Fourier series expansion of g(x)on[−l, l] would contain sine
terms only. This series is the sine series expansion of f (x).
Showthefollowinghelpfulshortcut:the sine seriescoefficients
are
b
k
=
1
l
l
−l
g(x)sin
kπ
l
dx =
2
l
l
0
f (x) sin
kπ
l
dx.
Then compute the sine series expansion of f (x) = x
2
on
[0, 1] and graph the limit function on [−3, 3]. Analogous to
the above, develop a Fourier cosine series and find the cosine
series expansion of f (x) = x on [0, 1].
2. Fourier series area part of thefield of Fourier analysis, which
is central to many engineering applications. Fourier analysis
includes the Fourier transforms (and the FFT or Fast Fourier
Transform) and inverse Fourier transforms, to which you will
get a brief introduction in this exercise. Given measurements
of a signal (waveform), the goal is to construct the Fourier
series of a function. To start with a simple version of the prob-
lem, suppose the signal has the form f (x) =
a
0
2
+ a
1
cosπ x +
a
2
cos2π x +b
1
sinπ x +b
2
sin2π x and you have the mea-
surements f (−1) = 0, f (−
1
2
) = 1, f (0) = 2, f (
1
2
) = 1 and
f (1) = 0. Substitutingintothegeneralequationfor f (x),show
that f (−1) =
a
0
2
− a
1
+ a
2
= 0. Similarly,
a
0
2
− a
2
− b
1
= 1,
a
0
2
+ a
1
+ a
2
= 2,
a
0
2
− a
2
+ b
1
= 1, and
a
0
2
− a
1
+ a
2
= 0.
Note that the first and last equations are identical and that
b
2
never appears in an equation. Thus, you have four equa-
tionsand four unknowns.Solvethe equations. Youshould con-
clude that f (x) = 1 + cos π x + b
2
sinπ x, with no informa-
tionaboutb
2
.Fortunately,thereisan easier wayof determining
the Fourier coefficients. Recall that a
k
=
!
1
−1
f (x) cos kπ xdx
and b
k
=
!
1
−1
f (x) sin kπ xdx. You can estimate the in-
tegral using function values at x =−1/2, x = 0, x = 1/2
and x = 1. Find a version of a Riemann sum ap-
proximation that gives a
0
= 2, a
1
= 1, a
2
= 0 and b
1
= 0.
What value is given for b
2
? Use this Riemann sum
rule to find the appropriate coefficients for the data
f (−
3
4
) =
3
4
, f (−
1
2
) =
1
2
, f (−
1
4
) =
1
4
, f (0) = 0, f (
1
4
)=−
1
4
,
f (
1
2
) =−
1
2
, f (
3
4
) =−
3
4
and f (1) =−1. Compare to the
Fourier series of exercise 9.
3. Fourier series have been used extensively in processing
digital information, including digital photographs as well as
music synthesis. A digital photograph stored in “bitmap” for-
mat can be thought of as three functions f
R
(x, y), f
G
(x, y)
and f
B
(x, y). For example, f
R
(x, y) could be the amount of
red content in the pixel that contains the point (x, y). Briefly
explain what f
G
(x, y) and f
B
(x, y) would represent and how
the three functions could be combined to create a color pic-
ture. A sine series for a function f (x) on the interval [0, L]is
∞
k=1
b
k
sin
kπ x
L
where b
k
=
2
L
!
L
0
f (x) sin
kπ x
L
dx. Describe
what a sine series for a function f (x, y) with 0 ≤ x ≤ L and
0 ≤ y ≤ M would look like. If possible, take your favorite
photograph in bitmap format and write a program to find
Fourier approximations. The accompanying images were cre-
ated in this way. The first three images show Fourier approxi-
mationswith2,10and50terms, respectively.Noticethatwhile
the 50-term approximation is fairly sharp, there are some rip-
ples(or“ghosts”) outliningthetwopeople;the ripplesaremore
obviousinthe10-termimage.Brieflyexplainhowtheseripples
relate to the Gibbs phenomenon.
2 terms 10 terms 50 terms