Mathematical Methods (Second Year) MT 2009: −π < x < π; −π < x < π; −π < x < 0 x 0 < x < π
Mathematical Methods (Second Year) MT 2009: −π < x < π; −π < x < π; −π < x < 0 x 0 < x < π
Mathematical Methods (Second Year) MT 2009: −π < x < π; −π < x < π; −π < x < 0 x 0 < x < π
A. Fourier Series
1. Sketch graphs of the following functions in the range −2π < x < 2π given that all the functions
are periodic with period 2π:
(a) f (x) = |x|, −π < x < π;
(b) f (x) = x, −π < x < π;
n
0 −π < x < 0
(c) f (x) =
x 0<x<π
2. State whether the following functions are even, odd or neither even nor odd:
(a) sin(x) (b) cos(x) (c) sin2 (x) (d) cos2 (x) (e) sin(2x) (f) cos(2x)
(g) x cos(x) (h) sin(x) + cos(x) (i) x + sin(x) (j) ex
0 a 2a
0 a 2a 4a
0 a 2a
Figure 2: Two different periodic continuations of f (x) with periods 2a and 4a respectively.
Both continuations coincide with f (x) on the interval [0, a] and both continuations can be ex-
pressed in terms of Fourier series. These Fourier series will be different, but both of them will
coincide with f (x) on the interval [0, a]! In the simplest case one has a function g(x) defined
on [0, L] with g(0) = g(L). Then there are two generic ways of continuing this function to
the interval: a symmetric continuation g1 (x) = g1 (−x) = g(x) and and antisymmetric one
g2 (x) = −g2 (−x) = g(x). The Fourier series for g1 (x) will only contain cosines and is called the
Fourier cosine series. The Fourier series for g2 (x) will only contain sines and is called the Fourier
sine series.
1
4. Find (a) the Fourier sine series and (b) the Fourier cosine series for the function
In order to determine the Fourier sine series, you need to extend the function to the interval
[−π, 0] by requiring that f (−x) = −f (x). Now determine the Fourier sine series for the resulting
odd function on the interval [−π, π]. This Fourier series involves only sines as the extended
function is odd. The Fourier sine series equals f (x) on the original interval [0, π].
5. Find (a) the Fourier sine series and (b) the Fourier cosine series for the function f (x) = x2 in
the interval 0 < x < 2.
7.∗ N.B. This problem is for students who have already taken the complex analysis short option
(and they are strongly advised to do it).
Evaluate the Fourier transforms of the following functions for a > 0 (use the reside theorem
where appropriate)
2a 1
(a) f (x) = a2 +x2
, (b) f (x) = exp(−a|x|), (c) f (x) = cosh(x) .
8.∗ A direct application of Fourier transform techniques is Fraunhofer diffraction. See the discussion
on p. 443 of Riley/Hobson/Bence. Consider a beam of monochromatic light of wavelength
λ incident on a two dimensional screen of width 2Y . The direction of the incident beam is
y
k’
Y
k θ
0
x
−Y
specified by the wave vector ~k. The magnitude of ~k is given by the wave number k = |~k| = 2π λ .
The essential quantity in the Fraunhofer diffraction pattern is the dependence of the observed
amplitude (and hence intensity) on the angle θ between the viewing direction (which we specify
through a vector ~k 0 ) and the direction of the incident beam ~k. We suppose that at the position
2
(0, y) the amplitude of the transmitted light is f (y) per unit length in the y-direction (note that
f (y) may be complex). The function f (y) is called an aperture function. Both the screen and
the beam are assumed infinite in the z-direction. The total light amplitude at a position
with x0 > 0 will be the superposition of all Huygens’ wavelets originating from the various parts
of the screen. For very large |~r0 these can be treated as plane waves to give
³ ´
Z Y exp i~k 0 · (~r0 − y~ey )
A(~r0 ) = dy f (y) . (3)
−Y |~r0 − y~ey |
³ ´
Here the factor exp i~k 0 · (~r0 − y~ey ) represents the phase change undergone bby the light in
travelling from the point y~ey on the screen to the point ~r0 and the denominator represents the
reduction in amplitude with distance. Writing ~k 0 = k cos(θ)~ex + k sin(θ)~ey and assuming that
|~r0 | À Y we can approximate |~r0 − y~ey | ≈ |~r0 | and hence approximate (3) as
Z ³ ´
exp i~k 0 · ~r0 Y
A(~r0 ) ≈ dy f (y) exp −iky sin(θ) (4)
|~r0 | −Y
Now using that f (y) = 0 for |y| > Y we can extende the integration limits to ±∞. Then the
intensity in direction θ is given by
2π ˜ 2π b
I(θ) = |A|2 = 2
|f (k sin(θ))|2 ≡ I(θ). (5)
|~r0 | |~r0 |2
b for
Use this to calculate I(θ)
(a) The “single finite slit” diffraction pattern, given by the aperture function
(b) The “two narrow slits” diffraction pattern, given by the aperture function
= 1 for (−a−b)
2 ≤ x ≤ −a+b
2
= 0 otherwise
where b < a. How is the result (c) related to those for (a) and (b)? Can you explain this in
terms of the convolution theorem?
d2 Φ
− K 2 Φ = f (x) , (6)
dx2
show that its solution can be written as
Z
1 ∞
f˜(k)
Φ(x) = − √ dk eikx . (7)
2π −∞ k2 + K 2
3
10. Show that
(a) The Fourier transform of f (ax) is a1 f˜(k/a);
(b) The Fourier transform of f (a + x) is eika f˜(k);
(c) The Fourier transform of eiqx f (x) is f˜(k − q);
df (x)
(d) The Fourier transform of dx is ikaf˜(k);
˜
(e) The Fourier transform of xf (x) is i dfdk(k) .
4V0 X 1 ³ nπy ´
V (x, y) = sin e−nπx/a .
π n a
n=1,3,...