ГОСТ Р 53608—2009
ет существенно уменьшить число определяемых параметров без ущерба для точности координатных
определений. Основными методами разрешения неоднозначности фазовых измерений псевдодаль
ности являются методы, базирующиеся на использовании вторых разностей фаз. Вторые разности
вычисляются с использованием первых разностей фаз.
3.4 Решения, получаемые по вторым разностям фаз. как правило, позволяют при нормальных
условиях наблюдений вычислитьзначения целочисленных параметров комбинаций измеренных фаз и
на этой основе определять относительные координаты определяемых пунктов с точностью порядка
первых единиц миллиметров.
3.5 Исходными данными, необходимыми для разрешения неоднозначности фазовых измерений,
являются:
- фазовые измерения надвух несущихчастотахдля всех наблюдаемых навигационныхспутников
ГЛОНАСС и GPS на интервале непрерывного слежения за фазой:
- абсолютные геоцентрические координаты исходного пункта;
- предварительные координаты определяемого пункта.
4 Алгоритм разрешения неоднозначности
4.1 Предварительная обработка фазовых измерений
В измеренные значения фаз несущих вводятся поправки за влияние тропосферы, ионосферы и
вращения Земли вокруг своей оси.
Поправки за влияние тропосферы вычисляются одним из существующих методов (например,
методы Саастомойнена иХопфилда). При длине измеряемогобазиса не более 5 км введение поправок
за тропосферу может быть пропущено.
Поправки за влиянио ионосферы могут вычисляться с использованием двухчастотных синхрон
ных измерений кодовых и фазовых псевдодальностей. Поправки за влияние ионосферы для фазовых
ПД численно равны поправкам, полученным для кодовых ПД. но знаки их противоположны (поправки к
кодовым ПД отрицательны, а поправки к фазовым ПД положительны). В случае измерений кодовых ПД
формулы для вычисления этих поправок соответственно на первой и второй частотах имеют вид:
f2- г
Ad™ =(d ,-d 2),
<
1
)
Г1 ’2
Ad™(d ,-d 2).
(
2
)
где f, , f2— значения первой и второй несущих частот;
dv d2 — значения измеренных кодовых ПД соответственно на первой и второй частотах.
Поправка на вращение Земли вокруг своей оси за время распространения сигнала от КА до пун
кта наблюдения (в метрах) вычисляется по формуле:
шз_
Дd*»
(3)
---------р sinA cosВ sinZ,
п
где <о— угловая скорость вращения Земли вокруг своей оси:
ае— большая полуось общеземного эллипсоида;
р — расстояние до КА (в метрах);
с — скорость света, равная 299792498 м/с;
В — широта пункта наблюдения;
A, Z — соответственно азимут и зенитное расстояние КА.
4.2 Определение первых разностей фаз
Исходным соотношением служит формула:
0 , ( И -
где;’— номер спутника:
номер момента времени измерения /,•;
1,2 — номера исходного и определяемого пунктов;
величина ПРФ;
т\ (i)’m2(i) —
измеренные значения фаз.
(4)
2