08.06.2013 Views

Bernese GPS Software Version 5.0 - Bernese GNSS Software

Bernese GPS Software Version 5.0 - Bernese GNSS Software

Bernese GPS Software Version 5.0 - Bernese GNSS Software

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

6.5 Preprocessing Phase Observations<br />

The change of the receiver clock (δk(t2) − δk(t1)) is computed using the linear combination<br />

selected in option “Frequency for the solution” (usually the ionosphere-free combination L3)<br />

for each epoch-difference.<br />

Now, we check whether the no-cycle slip hypothesis b1 = b2 = 0 holds. The residual in L3<br />

(ionosphere-free) linear combination is computed as<br />

r3 = β1 r1 + β2 r2 , where β1 = f2 1<br />

f 2 1 − f2 2<br />

The following condition should be met:<br />

<br />

|r3| ≤ 3g<br />

and β2 = − f2 2<br />

f 2 1 − f2 2<br />

. (6.5)<br />

(β1σ1) 2 + (β2σ2) 2 (6.6)<br />

When screening baseline files the factor g = √ 8 = √ 2 3 is due to triple-differencing (two<br />

receivers, two satellites, two epochs). When screening zero-difference files the factor is g =<br />

√ 4 = √ 2 2 instead (one receiver, two satellites, two epochs).<br />

Eqns. (6.3) resp. (6.4) allow us to compute the change of ionospheric refraction between the<br />

epochs t1 and t2<br />

I ij<br />

kℓ (t2) − I ij<br />

kℓ (t1)<br />

independently for both carriers (we assume b1 = b2 = 0 at present). The mean value m is<br />

computed as<br />

m = 1<br />

<br />

r1 +<br />

2<br />

f2 2<br />

f2 <br />

r2<br />

1<br />

(6.7)<br />

We check whether the condition<br />

m ≤ Mion<br />

(6.8)<br />

is met. The value of Mion (option “Maximum ionospheric change from epoch to epoch”) and the<br />

a priori RMS errors of the zero difference observables σ1 and σ2 (options “Sigma for L1/L2<br />

observations”) are input parameters in panel “MAUPRP 8: Cycle Slip Detection/Correction”. If<br />

conditions (6.6) and (6.8) hold, the no-cycle-slip hypothesis is accepted as true.<br />

In the opposite case a search over the values b1 and b5 is performed. All combinations<br />

<br />

r1<br />

b1p = NINT<br />

λ1<br />

<br />

r1<br />

b5q = NINT<br />

λ1<br />

<br />

+ p , p = −J1,...,−1,0,1,... ,J1<br />

− r2<br />

<br />

+ q , q = −J5,...,−1,0,1,... ,J5<br />

λ2<br />

(NINT = nearest integer) are formed and the “corrected” residuals<br />

(6.9)<br />

r1p = r1 − b1pλ1 , r2pq = r2 − (b1p − b5q)λ2 (6.10)<br />

are tested in the same way as the original residuals r1 and r2. The user has to specify<br />

the search ranges J1 and J5 (see “Search width to find L1/L5 cycle slip correction” in panel<br />

“MAUPRP 8: Cycle Slip Detection/Correction”). If one combination of r1p, r2pq meets the nocycle-slip<br />

hypothesis, this pair is assumed to be the cycle slip correction.<br />

If no combination for r1p, r2pq is found meeting the no-cycle-slip hypothesis a new ambiguity<br />

parameter should be introduced. But introducing too many ambiguity parameters would<br />

<strong>Bernese</strong> <strong>GPS</strong> <strong>Software</strong> <strong>Version</strong> <strong>5.0</strong> Page 119

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!