booklet format - inaf iasf bologna
booklet format - inaf iasf bologna
booklet format - inaf iasf bologna
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
A.A. 2011/2012<br />
Temporal Data Analysis<br />
2.6 Convolution — Page 29<br />
Gaussian frequency distribution<br />
Let us assume we have f (t) = cosω 0 t, and the frequency ω 0 is not precisely defined, but is<br />
Gaussian distributed:<br />
What we are really measuring is<br />
P(ω) = 1<br />
σ 1 ω 2<br />
2π e− 2 σ 2<br />
∫ +∞<br />
f ¯<br />
1<br />
(t) =<br />
−∞ σ 1 ω 2<br />
2π e− 2 σ 2 cos(ω − ω 0 )t dω (A.1)<br />
i.e. a convolution integral in ω 0 . Instead of calculating this integral directly, we use the inverse of<br />
the Convolution Theorem (2.49), thus saving work and gaining higher enlightenment. But watch<br />
it! We have to handle the variables carefully. The time t in (A.1) has nothing to do with the<br />
Fourier trans<strong>format</strong>ion we need in (2.49). And the same is true for the integration variable ω.<br />
Therefore, we rather use t 0 and ω 0 for the variable pairs in (2.49). We identify:<br />
F (ω 0 ) =<br />
1<br />
σ 1 ω 2 0<br />
2π e− 2 σ 2<br />
1<br />
2π G(ω 0) = cosω 0 t<br />
The inverse Fourier transform of F (ω 0 ) and G(ω 0 ) are<br />
Finally we get:<br />
f (t 0 ) = 1 1 2π e− 2 σ2 t0<br />
2<br />
[ δ(t0 − t)<br />
g (t 0 ) = 2π<br />
2<br />
h(t 0 ) = e − 1 2 σ2 t 2 0<br />
[ δ(t0 − t)<br />
Now the only thing left is to Fourier transform h(t 0 ):<br />
¯ f (t) ≡ H(ω 0 ) =<br />
∫ +∞<br />
e − 1 2 σ2 t0<br />
2<br />
−∞<br />
= e − 1 2 σ2 t 2 cosω 0 t<br />
2<br />
[ δ(t0 − t)<br />
2<br />
+ δ(t ]<br />
0 + t)<br />
2<br />
+ δ(t ]<br />
0 + t)<br />
2<br />
+ δ(t ]<br />
0 + t)<br />
e −iω 0t 0<br />
dt 0<br />
2<br />
Now, this was more work than we had originally thought it would be. But look at what we have<br />
gained in insight!<br />
M.Orlandini 123