PHY431 Slides FourierTransforms OptionalReview
PHY431 Slides FourierTransforms OptionalReview
1
= [exp(−iω / 2) − exp(iω/2)]
−iω
1 exp(iω / 2) − exp(−iω/2)
=
(ω/2) 2i
sin(ω/2) F(w)
=
(ω/2)
F (ω ) = sinc(ω/2) Imaginary
Component = 0
w
Example: the Fourier Transform of a
Gaussian, exp(-at2), is itself!
∞
F {exp(−at 2 )} = ∫ − ) exp(−iωt ) dt
2
exp( at
−∞
∝ exp(−ω / 4a)
2
The details are a HW problem!
0
t
0
w
Fourier Series & The Fourier Transform
The spectrum
∫ ∫
1
f (t ) = F (ω ) exp(iω t ) dω F (ω ) = f (t ) exp(−iω t ) dt
2π
−∞ −∞
Source: Prof. Rick Trebino, Georgia Tech
What do we hope to achieve with the
Fourier Transform?
We desire a measure of the frequencies present in a wave. This will
lead to a definition of the term, the spectrum.
Time
It will be nice if our measure also tells us when each frequency occurs.
Lord Kelvin on Fourier’s theorem
Lord Kelvin
Joseph Fourier
Fourier was
obsessed with the
physics of heat and
developed the
Fourier series and
transform to model
heat-flow problems.
sin(3wt)
Here, we write a
square wave as
a sum of sine
waves.
sin(5wt)
Any function can be written as the
sum of an even and an odd function.
E(-x) = E(x)
E ( x) ≡ [ f ( x) + f (− x)] / 2
O( x) ≡ [ f ( x) − f (− x)] / 2
O(-x) = -O(x)
f ( x) = E ( x) + O( x)
Fourier Cosine Series
Because cos(mt) is an even function (for all m), we can write an even
function, f(t), as:
∑
1
f(t) = F cos(mt)
π m
m =0
And where we’ll only worry about the function f(t) over the interval
(–π,π).
The Kronecker delta function
⎧1 if m = n
δ m,n ≡⎨
⎩0 if m ≠ n
Finding the coefficients, Fm, in a Fourier Cosine Series
∞
∑
1
Fourier Cosine Series: f (t ) = F cos(mt )
π m =0
m
To find Fm, multiply each side by cos(m’t), where m’ is another integer, and integrate:
π ∞ π
∑
1
∫
−π
f (t ) cos(m ' t ) dt =
π ∫ F
m=0 −π
m cos(mt ) cos(m ' t ) dt
π
⎧π if m = m '
But:
∫
−π
cos(mt ) cos(m ' t ) dt = ⎨
⎩ 0 if m ≠ m '
≡ π δ m,m '
π ∞
∑F π δ
1
So:
∫
−π
f (t ) cos(m ' t ) dt =
π m=0
m m,m '
ß only the m’ = m term contributes
Fm =
∫
−π
f (t ) cos(mt ) dt ß yields the
coefficients for
any f(t)!
Fourier Sine Series
∑
∞
1
f (t) = Fm′ sin(mt)
π
m= 0
where we’ll only worry about the function f(t) over the interval (–
π,π).
Finding the coefficients, F’m, in a Fourier Sine Series
∞
∑
1
Fourier Sine Series: f (t ) = F ′ sin(mt )
π m =0
m
To find Fm, multiply each side by sin(m’t), where m’ is another integer, and integrate:
π ∞ π
∑
1
But:
∫
−π
f (t ) sin(m ' t ) dt =
π ∫ F ′ sin(mt ) sin(m ' t ) dt
m =0 −π
m
π
⎧π if m = m '
∫
−π
sin(mt ) sin(m ' t ) dt = ⎨
⎩ 0 if m ≠ m '
≡ π δ m,m '
π ∞
So:
∑
1
∫ f (t ) sin(m ' t ) dt
−π
=
π
F′ π δ
m =0
m m,m ' ß only the m’ = m term contributes
∞ ∞
1 1
f (t ) =
π
∑ F
m =0
m cos(mt ) +
π
∑ F ′ sin(mt )
m =0
m
where
Fm =
∫ f (t) cos(mt) dt and Fm′ =
∫ f (t) sin(mt) dt
We can plot the coefficients of a Fourier Series
1
Fm vs. m
.5
0
5
10
25
15
20
30
m
We really need two such plots, one for the cosine series and another
for the sine series.
Discrete Fourier Series vs.
Continuous Fourier Transform
Fm vs. m
Let the integer
m become a F(m)
real number
and let the
coefficients,
Fm, become a
function F(m).
Again, we really need two such plots, one for the cosine series and
another for the sine series.
The Fourier Transform
Consider the Fourier coefficients. Let’s define a function F(m) that
incorporates both cosine and sine series coefficients, with the sine
series distinguished by making it the imaginary component:
F(m) º Fm – i F’m =
∫ f (t ) cos(mt ) dt − i
∫ f (t ) sin(mt ) dt
∫
The Fourier
F (ω ) = f (t ) exp ( −iω t ) dt Transform
−∞
F(w) is called the Fourier Transform of f(t). It contains equivalent
information to that in f(t). We say that f(t) lives in the time domain,
and F(w) lives in the frequency domain. F(w) is just another way of
looking at a function or wave.
The Inverse Fourier Transform
The Fourier Transform takes us from f(t) to F(w).
How about going back?
∑ ∑
1 1
f (t ) = Fm cos(mt ) + Fm' sin(mt )
π m =0
π m =0
∞
1 Inverse
f (t ) =
2π ∫
−∞
F (ω ) exp(iω t ) dω Fourier
Transform
The Fourier Transform and its Inverse
F (ω ) =
∫ f (t ) exp(−iωt ) dt
−∞
FourierTransform
∞
1
f (t ) =
2π ∫ F (ω ) exp(iωt ) dω
−∞
Inverse Fourier Transform
E (t ) → F {E (t )} or: E (t ) → E%(ω )
2
S (ω ) ≡ F {E (t )}
1
= [exp(−iω / 2) − exp(iω/2)]
−iω
1 exp(iω / 2) − exp(−iω/2)
=
(ω/2) 2i
sin(ω/2) F(w)
=
(ω/2)
F (ω ) = sinc(ω/2) Imaginary
Component = 0
w
Example: the Fourier Transform of a
Gaussian, exp(-at2), is itself!
∞
F {exp(−at 2 )} = ∫ − ) exp(−iωt ) dt
2
exp( at
−∞
∝ exp(−ω / 4a)
2
The details are a HW problem!
0
t
0
w
The Dirac delta function
d(t)
⎧∞ if t = 0
δ (t ) ≡ ⎨
⎩ 0 if t ≠ 0
t
⎧∞ if t = 0
The Dirac delta function δ (t ) ≡ ⎨
⎩ 0 if t ≠ 0
fm(t) = m exp[-(mt)2]/√p
d(t)
f3(t)
f2(t)
f1(t)
t
Dirac d-function Properties
d(t)
∫ δ (t ) dt = 1
−∞
t
∞ ∞
∫ δ (t − a) f (t ) dt = ∫ δ (t − a) f (a) dt = f (a)
−∞ −∞
∫ exp(±iωt ) dt = 2π δ (ω )
−∞
∞
∫ exp[±i(ω − ω ′)t ]
−∞
dt = 2π δ (ω − ω ′)
The Fourier Transform of d(t) is 1.
∞
d(t) 1
0
t w
∞
And the Fourier Transform of 1 is 2pd(w):
∫ 1 exp(−iωt ) dt = 2π δ (ω )
−∞
1 2pd(w)
t 0
w
The Fourier transform of exp(iw0 t)
∞
F {exp(iω0 t )}
∞
=
∫
−∞
exp(iω0 t ) exp(−i ω t ) dt
=
∫
−∞
exp(−i [ω − ω0 ] t ) dt = 2π δ (ω − ω0 )
exp(iw0t)
F {exp(iw0t)}
Im t
0
Re w
t 0 w0
0
F {cos(ω0t )} =
∫ cos(ω t ) exp(−i ω t ) dt
−∞
0
∞
1
=
2 ∫
−∞
[exp(i ω0 t ) + exp(−i ω0 t )] exp(−i ω t ) dt
∞ ∞
1 1
=
2 ∫ exp(−i [ω − ω ]t ) dt
−∞
0 +
2 ∫ exp(−i [ω + ω ] t ) dt
−∞
0
= π δ (ω − ω0 ) + π δ (ω + ω0 )
cos(w0t) F {cos(ω0t )}
t
0 w
-w0
0
+w0
Fourier Transform Symmetry Properties
Expanding the Fourier transform of a function, f(t):
∞
F (ω ) =
∫ [Re{ f (t)} + i Im{ f (t)}] [cos(ωt) − i sin(ωt)] dt
−∞ ∞
Even functions of w
Odd functions of w
The Modulation Theorem:
The Fourier Transform of E(t) cos(w0 t)
∞
F {E (t ) cos(ω0t )}
∞
=
∫
−∞
E (t ) cos(ω0t ) exp(−i ω t ) dt
1
∞
=
∫ E (t ) ⎡⎣exp(i ω t ) + exp(−i ω t )⎤⎦ exp(−i ω t ) dt
2
−∞
0
∞
0
1 1
2∫ 2∫
= E (t ) exp(−i [ω − ω ] t ) dt +0 E (t ) exp(−i [ω + ω ] t ) dt 0
−∞ −∞
1 % 1 %
F {E (t ) cos(ω0t )} = E (ω − ω0 ) + E (ω + ω0 )
2 2
Example:
E (t ) cos(ω0t ) F {E(t ) cos(ω0t )}
E(t) = exp(-t2)
t
-w0 w
0
w0
Scale Theorem
The Fourier transform
of a scaled function, f(at): F { f (at )} = F (ω /a) / a
∞
Proof: F { f (at )} =
∫
−∞
f (at ) exp( −iω t ) dt
=
∫ f (u) exp(−i [ω /a] u) du / a
−∞
= F (ω /a) / a
If a < 0, the limits flip when we change variables, introducing a
minus sign, hence the absolute value.
f(t) F(w)
The Scale
Theorem Short
pulse
in action
t w
The shorter
Medium-
the pulse, length
the broader pulse
the spectrum! t w
g(t) G(w)
F {a f (t ) + b g (t )} = w
t
aF { f (t )} + bF {g (t )}
F(w) +
f(t)+g(t) G(w)
Change variables : u = t − a
∞
∫
−∞
f (u ) exp(−iω[u + a ])du
∞
= exp(−iω a ) ∫ f (u ) exp(−iωu )du
−∞
= exp(−iω a ) F (ω )
Fourier Transform with respect to space
∞
F (k ) = ∫ f ( x) exp(−ikx) dx x
−∞
F {f(x)} = F(k)
∫∫
= f(x,y) exp[-i(kxx+kyy)] dx dy y
x
F (2){f(x,y)}
If f(x,y) = fx(x) fy(y),
f(0)
Effective width ≡ Area / height:
∞ Dteff
1 (Abs value is
Δteff ≡ ∫
f (0) −∞
f (t ) dt unnecessary
for intensity.)
0 t
∞ ∞
1 1 2π f (0)
Δω ≥ ∫
F (0) −∞
F (ω ) d ω = ∫
F (0) −∞
F (ω ) exp(iω [0]) d ω =
F (0)