Fourier Series
Fourier Series
Parseval’s Identity
L L
nπx nπx
Z Z
cos dx = 0 & sin dx = 0, (1)
−L L −L L
Z L
mπx nπx
sin cos dx = 0 (2)
−L L L
Z L (
mπx nπx 0, m 6= n
sin sin dx = (3)
−L L L L, m=n
Z L (
mπx nπx 0, m 6= n
cos cos dx = (4)
−L L L L, m=n
Thus, if we have
∞
a0 X n nπx nπx o
f (x) = + an cos + bn sin , x ∈ [0, 2L]
2 n=1
L L
2L
1 nπx
Z
an = f (x) cos dx, n = 0, 1, 2, . . . and (5)
L 0 L
2L
1 nπx
Z
bn = f (x) sin dx, n = 1, 2, . . . (6)
L 0 L
with
2L
1 nπx
Z
an = f (x) cos dx, n = 0, 1, 2, . . . and (8)
L 0 L
2L
1 nπx
Z
bn = f (x) sin dx, n = 1, 2, . . . (9)
L 0 L
Remark
The Fourier series of a function is defined whenever the integrals in (8)
and (9) have meaning.
MA 201 (2016): PDE
Applications
Parseval’s Identity
with
c+2L
1 nπx
Z
an = f (x) cos dx, n = 0, 1, 2, . . . and (11)
L c L
c+2L
1 nπx
Z
bn = f (x) sin dx, n = 1, 2, . . . (12)
L c L
Remark
The Fourier series of a function is defined whenever the integrals in (11)
and (12) have meaning.
MA 201 (2016): PDE
Applications
Parseval’s Identity
Gibbs’ Phenomenon
• This is about how the Fourier series of a piecewise continuously
differentiable periodic function behaves at a jump discontinuity. It is
named after the American mathematical physicist J. W Gibbs.
• Consider the following periodic function whose definition in one
period is:
1, 0 < x < 1,
f (x) = (13)
0, 1 < x < 2.
• This function can be represented as
∞
1 2 X sin(2n − 1)πx
f (x) = + , x ∈?. (14)
2 π n=1 2n − 1
• To see how well this infinite series represents the function, let us
truncate the series after N terms.
• Let the sum of these first N terms of the infinite series be denoted
by SN :
N
1 2 X sin(2n − 1)πx
SN = + . (15)
2 π n=1 2n − 1
MA 201 (2016): PDE
Applications
Parseval’s Identity
Gibbs’ Phenomenon
Gibbs’ Phenomenon
Gibbs’ Phenomenon
Gibbs’ Phenomenon
Gibbs’ Phenomenon
Gibbs’ Phenomenon
Definition
Let f (t) : [−L, L] → R be a piecewise C 1 -function. Define the adjusted function g (t)
as follows:
1
2 [f (t + ) + f (t − )], −L < t < L,
g (t) = (16)
1 + −
2 [f (−L ) + f (L )], t = ±L.
Note:
The above definition tells that g (t) coincides with f (t) at all points in (−L, L), where
f (t) is continuous, but g (t) is the average of the left-hand and right-hand limits of
f (t) at points of discontinuity in (−L, L).
Theorem
Let f : [−L, L] → R be a piecewise C 1 function and let g (t) be the
adjusted function as defined in (16). Then Fourier series of f (t) = g (t),
for all t ∈ [−L, L].
Note:
• The existence of Fourier series depends on the evaluation of Fourier
coefficients.
• On the other hand, convergence of the Fourier series is done via
adjusted function g .
• Thus, a highly oscillating function may be decomposed as sum of a
(possibly infinite) set of simple oscillating trigonometric functions.
• A discontinuous function is approximated by a continuous function.
A0 = h, An = 0, n = 1, 2, 3, . . . , (18)
2h
Bn = nπ , n odd, (19)
0, n even.
∞
A0 X h nπt nπt i
f (t) = + An cos + Bn sin ,
2 n=1
L L
∞
C0 X h nπt nπt i
g (t) = + Cn cos + Dn sin .
2 n=1
L L
∞ h
A0 X nπt nπt i
f (t)g (t) = g (t) + An g (t) cos + Bn g (t) sin .
2 n=1
L L
∞
" #
L L L L
A0 nπt nπt
Z Z X Z Z
f (t)g (t)dt = g (t)dt+ An g (t) cos dt + Bn g (t) sin dt
−L 2 −L n=1 −L L −L L
L ∞
1 A0 C0 X
Z
f (t)g (t)dt = + [An Cn + Bn Dn ].
L −L 2 n=1
L ∞
1 A20 X 2
Z
[f (t)]2 dt = + [An + Bn2 ].
L −L 2 n=1
• This result can obtained easily from the previous theorem by taking g (t) = f (t).
• The left-hand side represents the mean square value of f (t).
• It can, therefore, be thought of in terms of energy if f (t) represents a signal.
• What Parseval’s theorem states therefore is that the energy of a signal expressed
as a waveform is proportional to the sum of the squares of its Fourier coefficients.
P = EI = RI 2 .
∞
A0 X h nπt nπt i
I (t) = + An cos + Bn sin .
2 n=1
L L
L ∞
1 A20 1X 2
Z
Iav = I 2 (t)dt = + (A + Bn2 ).
2L −L 4 2 n=1 n
L ∞
1 A20 X 2
Z
[f (t)]2 dt = + [An + Bn2 ].
L −L 2 n=1
so that here L = π.
• We get
∞
1 π
2π 4 X 16
Z
t 4 dx = +
π −π 9 n=1
n4
giving us
∞
2π 4 2π 4 X 16
= +
5 9 n=1
n4
• Hence the result follows.
MA 201 (2016): PDE
Applications
Parseval’s Identity
Bernoulli’s Solution
• In [0, π], Bernoulli gave the solution of (21) as a series of the form
u = b1 sin x cos ct + b2 sin 2x cos 2ct + . . . . (24)
• When t = 0, we should have
φ(x) = b1 sin x + b2 sin 2x + . . . . (25)
This is possible and can be treated as the Fourier sine series of φ(x),
and which converges to φ(x).
Few Facts:
As we are expecting u(x, t) as a solution of (21), term-wise
differentiation of the series given by (24) should exist. In other sense, the
limit function u of the infinite series (24) should be differentiable.
Term-wise differentiation of an infinite series is not always possible.
In fact, infinite sum of continuous functions may have discontinuous
limit.
1 1 1 + x2
1+ 2
+ 2 2
+ ... = .
1+x (1 + x ) x2
MA 201 (2016): PDE
Applications
Parseval’s Identity
which fails at x = 0. In fact, the RHS series diverges for all x(?)
then ∞
X nπ n nπx nπx o
f 0 (x) = −An sin + Bn cos .
n=1
L L L
x x ∞ x
A0 nπτ nπτ o
Z Z X Z n
f (τ )dτ = dτ + An cos + Bn sin dτ.
−L −L 2 n=1 −L L L