A1 Time-Frequency Analysis: David Murray
A1 Time-Frequency Analysis: David Murray
A1 Time-Frequency Analysis: David Murray
A1 Time-Frequency Analysis
David Murray
david.murray@eng.ox.ac.uk
www.robots.ox.ac.uk/∼dwm/Courses/2TF
Hilary 2011
A1 2011 2/1
Course Contents
Topic 5. Introduction
It follows that
Z +T /2 ∞ ∞
1 X X
∗ −imωt
Ave power = Cn einωt Cm e dt
T −T /2 n=−∞ m=−∞
A1 2011 5/1
Parseval revisited
1
FT [f (t)g(t)] = F (ω) ∗ G(ω) .
2π
Writing out fully, this is
Z ∞ Z ∞
1
f (t)g(t)e−iωt dt = F (p)G(ω − p) dp .
−∞ 2π −∞
♣ Example
Q: Determine the energy in the signal f (t) = u(t)e−t (i) in the time
domain, and (ii) by determining the energy spectral density and
integrating over frequency.
A: Part (i): To find the total energy in the time domain
Z ∞ Z ∞ ∞
2 exp(−2t)
|f (t)| dt = exp(−2t)dt = dt
−∞ 0 −2 0
1 1
= 0− =
−2 2
Part (ii): In the frequency domain
Z ∞ Z ∞
F (ω) = u(t) exp(−t) exp(−iωt)dt = exp(−t(1 + iω))dt
−∞ 0
∞
exp(−t(1 + iω)) 1
= − =
(1 + iω)
0 (1 + iω)
1
Hence the energy spectral density (ESD) is |F (ω)|2 =
1 + ω2
A1 2011 13 / 1
♣ Example /ctd
Integration of ESD over all frequency f (not ω remember!!) gives
the total energy of
Z ∞ Z ∞
1 1
|F (ω)|2 df = dω
−∞ 2π −∞ 1 + ω2
Substitute tan θ = ω
Z ∞ Z π/2
1 1
|F (ω)|2 df = sec2 θdθ
−∞ 2π −π/2 1 + tan2 θ
Z π/2
1
= dθ
2π −π/2
1 1
= π =
2π 2
Which agrees with our earlier result.
A1 2011 14 / 1
Correlation, intro
Low
High High
g(t)
t
1 2 3 4 5 6
Auto-Correlation
Auto-correlation Rff , describes how self-similar a signal is.
Sum the products of the signal f (t) and a copy of the signal at a
shifted time f (t + τ ).
then sum
f(t +τ )
f(t)
f(t)
A1 2011 17 / 1
Now take the pairs of numbers at random from the function f (t).
Our result shows that there is no rearrangement, random or ordered,
of
R the function Rvalues into φ(t) that would make
f (t)φ(t)dt > f (t)2 dt.
This statement and proof are only valid for finite energy signals
Corollary of Wiener-Khinchin
If f (t) and g(t) are in the input and output of a system with
transfer function H(ω), then
Cross-correlation
Symmetries
The cross-correlation does not in general have a definite
reflection symmetry. However, we can swap the functions,
Rfg (τ ) = Rgf (τ ).
Independent signals
If Rfg (0) = 0, the signal f (t) and g(t) have no dependence on
one another.
A1 2011 26 / 1
♣ Example
[Q] Determine Rfg (τ ) f(t) g(t)
1 1
[A] Start by sketching g(t + τ ) as
function of t (four sketches t t
2a 2a 3a 4a
needed).
Red f (t): f = 0 then f = t/2a then f = 0
Blue g(t + τ ): g = 0 then g = t/a − (2 + τ /a) then g = 1 then g = 0
1 1 1 1
t t t t
2a+ τ 4a+ τ 2a+ τ 4a+ τ 2a+ τ 4a+ τ 2a+ τ 4a+ τ
0 2a 0 2a 0 2a 0 2a
Left-most non-zero config valid for 0 ≤ 4a + τ ≤ a, so that
For −4a ≤ τ ≤ −3a:
∞ 4a+τ
(4a + τ )2
Z Z
t
Rfg (τ ) = f (t)g(t + τ )dt = · 1 dt =
−∞ 0 2a 4a
A1 2011 27 / 1
♣ Example /ctd
1 1 1 1
t t t t
2a+ τ 4a+ τ 2a+ τ 4a+ τ 2a+ τ 4a+ τ 2a+ τ 4a+ τ
0 2a 0 2a 0 2a 0 2a
For −a ≤ τ ≤ 0:
Z ∞ Z 2a t
t
τ
Rfg (τ ) = f (t)g(t + τ )dt = · − 2+ dt
−∞ 2a+τ 2a a a
The integrals and finding the maximum is a DIY exercise ...
A1 2011 28 / 1
Application
First
R ∞ why is it2finite power? Because one cannot evaluate
−∞
| sin ω0 t| dt. It is not “square-integrable”.
has the same period as its parent function, and large when τ = 0
τ
increasing
A1 2011 32 / 1
Summary
An energy signal f (t) has finite total energy and zero average power
Z ∞
ETot = |f (t)|2 dt 0 < ETot < ∞ .
−∞
Summary
A power signal f (t) is one where the total energy is infinite, and we
consider average power
Z T
1
PAve = lim |f (t)|2 dt 0 < PAve < ∞ .
T →∞ 2T −T