Numerical Examples of the Diagonalisation of the Equations (GPS)

As discussed in Sect. 8.3.7, a normal equation can be diagonalised and the related observation equation can be formed.

For the linearised observation equation (cf. Eq. 8.38)

tmp1E-1034

the least squares normal equation can be written as (cf. Eqs. 8.39 and 8.40)

tmp1E-1035

The normal Eq. a2.2 can be diagonalised as (cf. Eq. 8.41)

tmp1E-1036


The above diagonalisation process can be repeated r – 1 times to the second normal equation of Eq. a2.4, so that the second equation of Eq. a2.4 can be fully diagonalised and Eq. a2.4 can be represented as:

tmp1E-1037

wheretmp1E-1038is a diagonal matrix, r is the dimension oftmp1E-1039is a vector.

Normal Eq. a2.4 related observation equation is (cf. Eq. 8.43)

tmp1E-1042

where E is an identity matrix, andtmp1E-1043are residual vectors, which have the same property as V in Eq. a2.1.

By similarly repeating the above process r – 1 times to the observation equation of X2 (i.e., the second equation of Eq. a2.8), then Eq. a2.8 turns out to have a form of

tmp1E-1045

wheretmp1E-1046is in a form of a diagonal matrix where all elements are vectors of dimension r, P is a diagonal matrix of P,L’is a vector of L, and 1*2 is a residual vector that has the same property as V in Eq. a2.1. Equation a2.11 is the observation equation of normal Eq. a2.7.

Numerical examples to illustrate the diagonalisation process of the normal equation and observation equation are given below.

1. The Case of Two Variables

For the observation equation (where o is set to 1, which does not affect all results)

tmp1E-1048

the least squares normal equation is

tmp1E-1049

Because

tmp1E-1050

Eq. a2.13 is diagonalised as

tmp1E-1051

The solutiontmp1E-1052of Eq. a2.14 is the same as that of Eq. a2.13. Furthermore, to form the equivalent observation equation, there are

tmp1E-1054

thus, the observation equation related to Eq. a2.14 is

tmp1E-1055

The normal equation of the observation Eq. a2.15 is exactly the same as Eq. a2.14. This numerical example shows that the normal equation and the related observation equation can be diagonalised.

2. The Case of Three Variables

For the observation equation (where o is set to 1, which does not affect all results)

tmp1E-1056

the least squares normal equation is

tmp1E-1057

Because

tmp1E-1058

Eq. a2.17 is diagonalised as

tmp1E-1059

Thetmp1E-1060related normal equation can be further diagonalised. Because of

tmp1E-1063

Eq. a2.19 is further diagonalised as

tmp1E-1064

The solutiontmp1E-1065of Eq. a2.20 is the same as that of Eqs. a2.17 and a2.19. Furthermore, to form the equivalent observation equation of Eq. a2.19, there are

tmp1E-1067

 

 

 

tmp1E-1068

thus, the observation equation related to Eq. a2.19 is

tmp1E-1069

The X2 and X3 related observation equation can be further diagonalised as follows. Because

tmp1E-1070

the observation equation related to Eq. a2.20 is

tmp1E-1071

The normal Eq. a2.17 and its related observation Eq. a2.16 are fully diagonalised as Eqs. a2.20 and a2.22, respectively. These numerical examples show that the normal equation and the related observation equation can be diagonalised as described in Sect. 8.3.7.

Next post:

Previous post: