12.540 Principles of The Global Positioning System: Prof. Thomas Herring

12.
540 Principles of the Global

Positioning System
Lecture 11
Prof. Thomas Herring
03/17/03
12.540 Lec 11
Statistical approach to estimation

Summary
Look at estimation from statistical point of
view
Propagation of covariance matrices
Sequential estimation
03/17/03
12.540 Lec 11
Statistical approach to estimation

Examine the multivariate Gaussian distribution:
1
- (x- m )T V -1 (x- m )
1
2
Multivariant f (x) =
e
(2p ) n V
Minimize
(x - m)T V-1 (x - m) gives largest probability density
By minimizing the argument of the exponential in the

probability density function, we maximize the
likelihood of the estimates (MLE).
This is just weighted least squares where the weight
matrix is chosen to be the inverse of the covariance
matrix of data noise
03/17/03
12.540 Lec 11
Data covariance matrix

If we use the inverse of the covariance matrix of the
noise in the data, we obtain a MLE if data noise is
Gaussian distribution.
How do you obtain data covariance matrix?
Difficult question to answer completely
Issues to be considered:
Thermal noise in receiver gives on component
Multipath could be treated as a noise-like quantity
Signal-to-noise ratio of measurements allows an estimate of
the noise (discussed later in course).
In-complete mathematical model of observables can
sometimes be treated as noise-like.
Gain of GPS antenna will generate lower SNR at low elevation
angles
03/17/03
12.540 Lec 11

In practice in GPS (as well as many other fields), the
data covariance matrix is somewhat arbitrarily chosen.
Largest problem is temporal correlations in the
measurements. Typical GPS data set size for 24hours of data at 30 second sampling is 8x2880=23000
phase measurements. Since the inverse of the
covariance matrix is required, fully accounting for
correlations requires the inverse of 23000x23000
matrix.
To store the matrix would require, 4Gbytes of memory
Even if original covariance matrix is banded (ie.,
correlations over a time short compared to 24-hours),
the inverse of banded matrix is usually a full matrix
03/17/03
12.540 Lec 11

Methods on handling temporal correlations:
If measurements correlated over say 5-minute period, then
use samples every 5-minutes (JPL method)
Use full rate data, but artificially inflate the noise on each
measurement so that equivalent to say 5-minute sampling (ie.,
sqrt(10) higher noise on the 30-second sampled values
(GAMIT method)
When looking a GPS results, always check the data noise
assumptions (discussed more near end of course).
Assuming a valid data noise model can be developed,

what can we say about noise in parameter estimates?
03/17/03
12.540 Lec 11
Propagation of covariances
Given a data noise covariance matrix, the
characteristics of expected values can be
used to determine the covariance matrix of
any linear combination of the measurements.
Given linear operation : y = Ax with Vxx as
covariance matrix of x
Vyy =< yy T >=< Axx T A T >= A < xx T > A T
Vyy = AVxx A T
03/17/03
12.540 Lec 11
Propagation of covariance
Propagation of covariance can be used for
any linear operator applied to random
variables whose covariance matrix is already
known.
Specific examples:
Covariance matrix of parameter estimates from
least squares
Covariance matrix for post-fit residuals from least
squares
Covariance matrix of derived quantities such as
latitude, longitude and height from XYZ coordinate
estimates.
03/17/03
12.540 Lec 11
Covariance matrix of parameter estimates

Propagation of covariance can be applied to the
weighted least squares problem:
-1
-1
x = (A T Vyy A) -1 A T Vyy y
-1
-1
-1
-1
< x x T >= (A T Vyy A) -1 A T Vyy < yy T > Vyy A(A T Vyy A) -1

-1
Vx x = (A T Vyy A) -1
Notice that the covariance matrix of parameter

estimates is a natural output of the estimator if ATV-1A
is inverted (does not need to be)
03/17/03
12.540 Lec 11
Covariance matrix of estimated

parameters
Notice that for the rigorous estimation, the inverse of the data
covariance is needed (time consuming if non-diagonal)
To compute to parameter estimate covariance, only the
covariance matrix of the data is needed (not the inverse)
In some cases, a non-rigorous inverse can be done with say a
diagonal covariance matrix, but the parameter covariance matrix
is rigorously computed using the full covariance matrix. This is a
non-MLE but the covariance matrix of the parameters should be
correct (just not the best estimates that can found).
This techniques could be used if storage of the full covariance
matrix is possible, but inversion of the matrix is not because it
would take too long or inverse can not be performed in place.
03/17/03
12.540 Lec 11
10
Covariance matrix of post-fit residuals

Post-fit residuals are the differences between the
observations and the values computed from the
estimated parameters
Because some of the noise in the data are absorbed
into the parameter estimates, in general, the post-fit
residuals are not the same as the errors in the data.
In some cases, they can be considerably smaller.
The covariance matrix of the post-fit residuals can be
computed using propagation of covariances.
03/17/03
12.540 Lec 11
11
Covariance matrix of post-fit residuals

This can be computed using propagation on
covariances: e is the vector of true errors, and
v is vector of residuals
y = Ax + e
x = (A T Vyy-1 A) -1 A T Vyy-1 y
T -1
-1 T -1
v = y - Ax = I - A(A Vyy A) A Vyy e Eqn 1

144424443
Amount error reduced
Vvv =< vvT >= Vyy - A(A T Vyy-1 A) -1 A T
03/17/03
12.540 Lec 11
12
Post-fit residuals
Notice that we can compute the compute the
covariance matrix of the post-fit residuals (a large
matrix in generate)
Eqn 1 on previous slide gives an equation of the form
v=Be; why can we not compute the actual errors with
e=B-1v?
B is a singular matrix which has no unique inverse
(there is in fact one inverse which would generate the
true errors)
Note: In this case, singularity does not mean that
there is no inverse, it means there are an infinite
number of inverses.
03/17/03
12.540 Lec 11
13
Example
Consider the case shown below: When a rate
of change is estimated, the slope estimate will
absorb error in the last data point particularly
as Dt increases. (Try this case yourself)
6
Example of fitting slope to non-uniform data distribution
Data
Dt
3
2
1
0
03/17/03
Postfit error bar very small;

slope will always pass close
to this data point
Postfit error bar

somewhat reduced
0.0
10.0
20.0
Time
30.0
12.540 Lec 11
40.0
50.0
14
Covariance of derived quantities

Propagation of covariances can be used to determine
the covariance of derived quantities. Example latitude,
longitude and radius. q is co-latitude, l is longitude, R
is radius. DN, DE and DU are north, east and radial
changes (all in distance units).
Geocentric Case :
DN -cos(q )cos( l ) -cos(q ) sin( l) sin(q ) DX

cos( l )
0 DY
DE = - sin( l)
DU
X /R
Y /R
Z / R DZ
14444444244444443
A matrix for use in propagation from Vxx
03/17/03
12.540 Lec 11
15
Estimation in parts/Sequential estimation

A very powerful method for handling large
data sets, takes advantage of the structure of
the data covariance matrix if parts of it are
uncorrelated (or assumed to be uncorrelated).
V1
0
0
03/17/03
0
V2
0
-1
V1 -1
0
0
0
-1
0 = 0
0
V2
-1
0
V3
0
V3
12.540 Lec 11
16
Since the blocks of the data covariance matrix can be
separately inverted, the blocks of the estimation
(ATV-1A) can be formed separately can combined
later.
Also since the parameters to be estimated can be
often divided into those that effect all data (such as
station coordinates) and those that effect data a one
time or over a limited period of time (clocks and
atmospheric delays) it is possible to separate these
estimations (shown next page).
03/17/03
12.540 Lec 11
17
Sequential estimation with division of global and local parameters.
V is covariance matrix of new data (uncorrelated with priori
parameter estimates), Vxg is covariance matrix of prior parameter
estimates with estimates xg and xl are local parameter estimates,
xg+ are new global parameter estimates.
y A g
=
x g I
A l x g

0 x l
-1
T
-1
-1
x +g (A g T V-1 A g + Vxg
A g V-1 A l A g T V-1 y + Vxg
xg
)
=

T
T
T
A l V-1 A g
A l V-1 A l
A l V-1 y
x l
03/17/03
12.540 Lec 11
18
As each block of data is processed, the local
parameters, xl, can be dropped and the
covariance matrix of the global parameters xg
passed to the next estimation stage.
Total size of adjustment is at maximum the
number of global parameters plus local
parameters needed for the data being
processed at the moment, rather than all of
the local parameters.
03/17/03
12.540 Lec 11
19
Summary
We examined the way covariance matrices
and be manipulated
Estimation from a statistical point of view
Sequential estimation.
Next class continue with sequential estimation
in terms of Kalman Filtering.
Reminder: Paper topic and outline due
Wednesday.
03/17/03
12.540 Lec 11
20
10

12.540 Principles of The Global Positioning System: Prof. Thomas Herring

Uploaded by

Copyright:

Available Formats

12.540 Principles of The Global Positioning System: Prof. Thomas Herring

Uploaded by

Document Information

Original Description:

Original Title

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

12.540 Principles of The Global Positioning System: Prof. Thomas Herring

Uploaded by

Copyright:

Available Formats

12.

540 Principles of the Global

Statistical approach to estimation

Statistical approach to estimation

(x - m)T V-1 (x - m) gives largest probability density

By minimizing the argument of the exponential in the

Data covariance matrix

Data covariance matrix

Data covariance matrix

Assuming a valid data noise model can be developed,

Covariance matrix of parameter estimates

< x x T >= (A T Vyy A) -1 A T Vyy < yy T > Vyy A(A T Vyy A) -1

Notice that the covariance matrix of parameter

Covariance matrix of estimated

Covariance matrix of post-fit residuals

Covariance matrix of post-fit residuals

v = y - Ax = I - A(A Vyy A) A Vyy e Eqn 1

Example of fitting slope to non-uniform data distribution

Postfit error bar very small;

Postfit error bar

Covariance of derived quantities

Estimation in parts/Sequential estimation

You might also like