Use of Kalman Filters in Time and Frequency Analysis
Use of Kalman Filters in Time and Frequency Analysis
frequency analysis
John Davis
1st May 2011
Overview of Talk
Simple example of Kalman filter
Discuss components and operation of filter
Models of clock and time transfer noise and
deterministic properties
Examples of use of Kalman filter in time and
frequency
Simple Example
1.5E-08
1.0E-08
Offset
5.0E-09
0.0E+00
0
10
12
-5.0E-09
-1.0E-08
Date
Truth
Measurement
14
16
18
20
1.0E-08
5.0E-09
0.0E+00
0
10
12
14
-5.0E-09
-1.0E-08
Truth
Measurement
Filter Estimate
16
18
20
Error
2.0E-09
0.0E+00
0
10
12
14
-2.0E-09
-4.0E-09
-6.0E-09
Date
Error
Uncertainty (lower)
Uncertainty (upper)
16
18
20
Offset
8.0E-09
6.0E-09
4.0E-09
2.0E-09
0.0E+00
-2.0E-09
10
12
Date
Measurement
Filter Estimate
14
16
18
20
x (tn ) ( ) x (tn1 )
( )
spacing
x (tn1 )
Estimate
of the state vector at iteration n-1
Offset
8.0E-09
6.0E-09
4.0E-09
2.0E-09
0.0E+00
-2.0E-09
10
12
14
Date
Measurement
Filter Estimate
16
18
20
Uncertainty Extrapolation
6.0E-09
4.0E-09
Error / Uncertainty
2.0E-09
0.0E+00
0
10
12
14
16
18
-2.0E-09
-4.0E-09
-6.0E-09
Error
Date
Uncertainty (upper)
Uncertainty Extrapolation (upper)
Uncertainty (lower)
Uncertainty Extrapolation (lower)
20
x (tn ) x (tn )
Design Matrix
Offset
8.0E-09
6.0E-09
4.0E-09
2.0E-09
0.0E+00
-2.0E-09
Measurement
Filter Estimate
10
12
Date
State Vector Extrapolation
14
16
18
20
P (t ) I K (t n ) H P (t )
K (tn )
Kalman Gain
Design Matrix
H
P (tn ) P(tn ) Parameter Covariance Matrix, at iteration
n, including and not including measurement n
respectively
I Identity Matrix
P (t ) I K (tn ) H P (t ) I K (tn ) H
K (tn )
Kalman Gain
Design Matrix
H
P (tn ) P(tn ) Parameter Covariance Matrix, at
iteration n, including and not including measurement
n respectively
I Identity Matrix
Error / Uncertainty
4.0E-09
2.0E-09
0.0E+00
0
10
12
14
16
18
-2.0E-09
-4.0E-09
-6.0E-09
Date
Error
Uncertainty Extrapolation (lower)
Uncertainty update (upper)
Uncertainty (lower)
Uncertainty Extrapolation (upper)
Uncertainty (upper)
Uncertainty Update (lower)
20
Q LQ (tn ) LQ (tn )
spacing
R LR (tn ) L (tn )
T
R
magnitude
-12
ADEV simulation
HDEV simulation
MDEV simulation
ADEV theory
HDEV theory
MDEV theory
z
-13
S ig m a
10
10
10
-14
-15
10
10
10
tau
10
10
10
Deterministic Properties
Linear frequency drift
Use three state vector components, time offset,
normalised frequency offset and linear frequency drift
2 2
0 1
0
0
1
x e
White Noise
n
Random Walk Noise
n
n 1
Markov Noise
0 k 1
n
n 1
Flicker noise, combination of Markov processes with
continuous range of relaxation parameters k
Markov noise processes and flicker noise have
memory
Kalman filters tend to work very well with Markov
noise processes
May only construct an approximation to flicker noise
in a Kalman filter
x x e
x kx e
x 10
-6
T im e O f s e t
4
2
0
-2
-4
-6
6
6.5
7
MJD
7.5
WPM Noise
-10
T im e O f s e t
0.5
-0.5
10
MJD
15
20
x 10
-11
T im e O f s e t
6
4
2
0
-2
-4
10
10.2
10.4
MJD
10.6
10.8
11
S ig m a
10
10
10
10
10
-12
ADEV Measurements
HDEV Measurements
MDEV Measurements
ADEV Theory
HDEV Theory
MDEV Theory
-13
-14
-15
-16
-17
10
10
10
tau
10
10
x 10
-10
1.5
1
Time Ofset
0.5
0
-0.5
-1
-1.5
-2
-2.5
10
MJD
15
20
S ig m a
10
10
10
10
-12
ADEV Measurements
HDEV Measurements
MDEV Measurements
ADEV Theory
HDEV Theory
MDEV Theory
-13
-14
-15
-16
10
10
10
tau
10
10
x 10
-10
T im e O f s e t
-1
-2
10
MJD
15
20
10
S ig m a
10
10
10
10
ADEV Measurements
HDEV Measurements
MDEV Measurements
ADEV Theory
HDEV Theory
MDEV Theory
-13
-14
-15
-16
10
10
10
tau
10
10
Clock Noise
x 10
WFM Noise
-8
T im e Ofset
-5
-10
-15
10
MJD
15
20
S ig m a
10
10
10
-11
ADEV Measurements
HDEV Measurements
MDEV Measurements
ADEV Theory
HDEV Theory
MDEV Theory
-12
-13
10
10
10
tau
10
10
x 10
RWFM Noise
-7
-1
T im e O fse t
-2
-3
-4
-5
-6
-7
-8
-9
0
10
MJD
15
20
-13
S ig m a
10
-12
10
10
ADEV Measurements
HDEV Measurements
MDEV Measurements
ADEV Theory
HDEV Theory
MDEV Theory
-14
-15
10
10
10
tau
10
10
FFM Noise
2
x 10
FFM Noise
-9
Time Ofset
0
-1
-2
-3
-4
-5
-6
-7
10
MJD
15
20
-14
S ig m a
10
ADEV Measurements
HDEV Measurements
MDEV Measurements
ADEV Theory
HDEV Theory
MDEV Theory
10
-15
10
10
10
tau
10
10
Periodic Instabilities
-5
1.622
x 10
Input Data
T im e O f s e t
1.6215
1.621
1.6205
1.62
1.6195
1.619
0
Day
D
GPS CV
T im e O f s e t
0.5
-0.5
-1
-1.5
MJD
10
11
12
13
10
10
10
10
-11
ADEV
HDEV
MDEV
-12
-13
-14
10
10
10
Tau
10
10
ADEV Measurements
HDEV Measurements
MDEV Measurements
ADEV Theory
HDEV Theory
MDEV Theory
-12
ADEV
10
-11
10
10
-13
-14
10
10
10
tau
10
10
( y Hx )
Residuals
Residual Variance (Theory)
V HP H
Time Ofset
0.5
0
-0.5
-1
-1.5
5.5
6.5
MJD
7.5
-9
Kalman
Filter Residuals (real data)
x 10
Time Ofset
1
0.5
0
-0.5
-1
-1.5
2.5
3.5
MJD
4.5
Computational Issues
x (t ) ( ) x (t
n 1
P (t ) ( ) P (t ) ( ) Q(tn )
n 1
ADEV = PED/Tau
ADEV Measurements
HDEV Measurements
MDEV Measurements
ADEV Theory
HDEV Theory
MDEV Theory
-12
ADEV
10
-11
10
10
-13
-14
10
10
10
tau
10
10
-11.4
-11.6
lo g P E D / ta u
-11.8
-12
-12.2
-12.4
-12.6
-12.8
-13
-13.2
-13.4
2.5
3.5
4
log Tau
4.5
5.5
10
10
10
10
-11
ADEV
HDEV
MDEV
-12
-13
-14
10
10
10
Tau
10
10
-11.4
-11.6
Add plot
lo g P E D / ta u
-11.8
-12
-12.2
-12.4
-12.6
-12.8
-13
-13.2
-13.4
2.5
3.5
4
log Tau
4.5
5.5
x 10
-8
T i m e -O f s e t (s )
5.5
5
4.5
4
3.5
3
2.5
5.4405
5.441
5.4415
MJD
5.442
5.4425
x 10
x 10
-8
T i m e -O f s e t (s )
4.5
3.5
5.4405
5.441
5.4415
MJD
5.442
5.4425
x 10
x 10
-9
T im e -O f s e t ( s )
4
2
0
-2
-4
-6
-8
-10
-12
5.4405
5.441
5.4415
MJD
5.442
5.4425
x 10
Summary
Title of Presentation
Name of Speaker
DateThe National Measurement System delivers world-class
measurement science & technology through these organisations