The GNSS integer ambiguities: estimation and validation
The GNSS integer ambiguities: estimation and validation
The GNSS integer ambiguities: estimation and validation
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
(appendix A.1):<br />
Qy = C sd<br />
pφ ⊗ D T D (2.60)<br />
with the transformation matrix D T as defined in (2.48).<br />
2.4.1 Variance component <strong>estimation</strong><br />
<strong>The</strong> stochastic model as specified by the vc-matrix of the observations can be given as:<br />
Qy = σ 2 Gy<br />
(2.61)<br />
with σ 2 the variance factor of unit weight, <strong>and</strong> Gy the cofactor matrix. If the variance<br />
factor is assumed to be unknown, it can be estimated a posteriori. For that purpose<br />
equation (2.61) can be generalized into:<br />
Qy =<br />
p<br />
α<br />
σ 2 αGα<br />
(2.62)<br />
Hence, the stochastic model is written as a linear combination of unknown factors <strong>and</strong><br />
known cofactor matrices, (Tiberius <strong>and</strong> Kenselaar 2000). <strong>The</strong> factors σ 2 α can represent<br />
variances as well as covariances, for example variance factors for the code <strong>and</strong> phase<br />
data, a variance for each satellite-receiver pair, <strong>and</strong> the covariance between code <strong>and</strong><br />
phase data. <strong>The</strong> problem of estimating the unknown (co-)variance factors is referred to<br />
as variance component <strong>estimation</strong>.<br />
2.4.2 Elevation dependency<br />
Until now the st<strong>and</strong>ard deviations of a certain observation type were chosen equal for<br />
all satellites, although it might be more realistic to assign weights to the observations of<br />
different satellites depending on their elevation. That is because increased signal attenuation<br />
due to the atmosphere <strong>and</strong> multipath effects will result in less precise observations<br />
from satellites at low elevations. <strong>The</strong>refore, Euler <strong>and</strong> Goad (1991) have suggested to<br />
use elevation-dependent weighting:<br />
σ 2 p s r = (qs r) 2 σ 2 p (2.63)<br />
with σp the st<strong>and</strong>ard deviation of code observations to satellites in the zenith, <strong>and</strong><br />
q s r = 1 + a · exp{ −εsr }, (2.64)<br />
ε0<br />
where a is an amplification factor, ε s r the elevation angle, <strong>and</strong> ε0 a reference elevation<br />
angle. <strong>The</strong> values of a <strong>and</strong> ε0 depend on the receiver that is used, <strong>and</strong> on the observation<br />
type – the function q in equation (2.64) may be different for code <strong>and</strong> phase observations<br />
as well as for observations on different frequencies.<br />
Experimental results have shown that an improved stochastic model can be obtained if<br />
elevation dependent weighting is applied, see e.g. (Liu 2002).<br />
<strong>GNSS</strong> stochastic model 21