allgosts.ru07. МАТЕМАТИКА. ЕСТЕСТВЕННЫЕ НАУКИ07.040. Астрономия. Геодезия. География

ГОСТ Р 53608-2009 Глобальная навигационная спутниковая система. Методы и технологии выполнения геодезических и землеустроительных работ. Разрешение неоднозначности фазовых измерений псевдодальности. Основные положения

Обозначение:
ГОСТ Р 53608-2009
Наименование:
Глобальная навигационная спутниковая система. Методы и технологии выполнения геодезических и землеустроительных работ. Разрешение неоднозначности фазовых измерений псевдодальности. Основные положения
Статус:
Действует
Дата введения:
01/01/2011
Дата отмены:
-
Заменен на:
-
Код ОКС:
07.040

Текст ГОСТ Р 53608-2009 Глобальная навигационная спутниковая система. Методы и технологии выполнения геодезических и землеустроительных работ. Разрешение неоднозначности фазовых измерений псевдодальности. Основные положения

ФЕДЕРАЛЬНОЕ АГЕНТСТВО ПО ТЕХНИЧЕСКОМУ РЕГУЛИРОВАНИЮ И МЕТРОЛОГИИ

НАЦИОНАЛЬНЫЙ

СТАНДАРТ

РОССИЙСКОЙ

ФЕДЕРАЦИИ

ГОСТР

53608-

2009

Глобальная навигационная спутниковая система

МЕТОДЫ И ТЕХНОЛОГИИ ВЫПОЛНЕНИЯ ГЕОДЕЗИЧЕСКИХ И ЗЕМЛЕУСТРОИТЕЛЬНЫХ

РАБОТ

Разрешение неоднозначности фазовых измерений

псевдодальности

Основные положения

Издание официальное

3

м

tfi

Москва

Стакдартинформ

2010

ГОСТ Р 53608—2009

Предисловие

Цели и принципы стандартизации в Российской Федерации установлены Федеральным законом от 27 декабря 2002 г. N9164-ФЗ «О техническом регулировании», а правила применения национальных стандартов Российской Федерации — ГОСТ Р 1.0—2004 «Стандартизация в Российской Федерации. Основные положения»

Сведения о стандарте

1 РАЗРАБОТАН 29-м Научно-исследовательским институтом Министерства обороны Российской Федерации

2 ВНЕСЕН Техническим комитетом по стандартизации ТК 363 «Радионавигация»

3 УТВЕРЖДЕН И ВВЕДЕН В ДЕЙСТВИЕ Приказом Федерального агентства по техническому регулированию и метрологии от 15 декабря 2009 г. № 933-ст

4 8ВЕДЕН ВПЕРВЫЕ

Информация об изменениях к настоящему стандарту публикуется е ежегодно издаваемом информационном указателе «Национальные стандарты». а текст изменений и поправок в ежемесячно издаваемых информационных указателях «Национальные стандарты». В случаях пересмотра (замены) или отмены настоящего стандарта соответствующее уведомление будет опубликовано в ежемесячно издаваемом информационном указателе «Национальные стандарты». Соответствующая информация, уведомления и тексты размещаются также в информационной системе общего пользования — на официальном сайте Федерального агентства по техническому регулированию и метрологии в сети Интернет

© Стандартинформ, 2010

Настоящий стандарт не может быть полностью или частично воспроизведен, тиражирован и распространен в качестве официального издания без разрешения Федерального агентства по техническому регулированию и метрологии

ГОСТ Р 53608—2009

Содержание

1 Область применения............................................1

2 Обозначения и сокращения........................................1

3 Общие положения.............................................1

4 Алгоритм разрешения неоднозначности.................................2

4.1 Предварительная обработка фазовых измерений.........................2

4.2 Определение первых разностей фаз.................................2

4.3 Определение вторых разностей фаз.................................3

4.4 Учет корреляции вторых разностей фаз...............................4

4.5 Получение целочисленных параметров вторых разностей фаз..................5

in

ГОСТ Р 53608—2009

НАЦИОНАЛЬНЫЙ СТАНДАРТ РОССИЙСКОЙ ФЕДЕРАЦИИ

Глобальная навигационная спутниковая система

МЕТОДЫ И ТЕХНОЛОГИИ ВЫПОЛНЕНИЯ ГЕОДЕЗИЧЕСКИХ И ЗЕМЛЕУСТРОИТЕЛЬНЫХ РАБОТ

Разрешение неоднозначности фазовых измерений псевдодальности

Основные положения

Global navigation satellite system. Methods and technologies of geodetic and cadastral works execution. Ambiguity resolution for phase pseudo-ranges. 8asic principles

Дата введения — 2011—01—01

1 Область применения

Настоящий стандарт распространяется на методы и технологии выполнения геодезических и землеустроительных работ с использованием аппаратуры потребителей глобальной навигационной спутниковой системы.

Настоящий стандарт устанавливает методы определения целочисленных параметров измеренных фазовых псевдодальностей до космических аппаратов спутниковых систем типа GPS. выполняемую при высокоточных определениях относительных координат пунктов с использованием наблюдений космических аппаратов этих систем, при выполнении геодезических и землеустроительных работ.

2 Обозначения и сокращения

8 настоящем стандарте применяют следующие сокращения:

ГЛОНАСС — глобальная навигационная спутниковая система Российской Федерации:

КА — космический аппарат;

МНК — метод наименьших квадратов;

ПД — псевдодальность:

СКП — среднее квадратическое отклонение:

GPS — глобальная навигационная спутниковая система Соединенных Штатов Америки.

3 Общие положения

3.1 Величина фазовой псевдодальности выражается в виде целого числа длин волн и дробной части длины волны. Дробная часть фазы может быть получена с высокой точностью, характеризуемой СКП на уровне первых единиц миллиметров. Прямое измерение целочисленной части фазы не представляется возможным, вследствие чего возникает неоднозначность фазовых измерений, требующая разрешения.

3.2 Разрешение неоднозначности фазовых измерений псевдодальности представляет собой процесс нахождения целого числа циклов (волн), укладывающихся в расстоянии между навигационным спутником и приемником при обработке фазовых измерений.

3.3 Для решения этой задачи с точностью, требуемой для выполнения геодезических и землеустроительных работ, применяют специальные методы математической обработки фазовых измерений, использующие при обработке не сами измерения, а их разности: первые, вторые и третьи. Это позволя-

Иэдание официальное

1

ГОСТ Р 53608—2009

ет существенно уменьшить число определяемых параметров без ущерба для точности координатных определений. Основными методами разрешения неоднозначности фазовых измерений псеедодаль* ности являются методы, базирующиеся на использовании вторых разностей фаз. Вторые разности вычисляются с использованием первых разностей фаз.

3.4 Решения, получаемые по вторым разностям фаз, как правило, позволяют при нормальных условиях наблюдений вычислить значения целочисленных параметров комбинаций измеренных фаз и на этой основе определять относительные координаты определяемых пунктов с точностью порядка первых единиц миллиметров.

3.5 Исходными данными, необходимыми для разрешения неоднозначности фазовых измерений, являются:

• фазовые измерения на двух несущих частотах для всех наблюдаемых навигационных спутников ГЛОНАСС и GPS на интервале непрерывного слежения за фазой:

• абсолютные геоцентрические координаты исходного пункта;

• предварительные координаты определяемого пункта.

4 Алгоритм разрешения неоднозначности

4.1 Предварительная обработка фазовых измерений

В измеренные значения фаз несущих вводятся поправки за влияние тропосферы, ионосферы и вращения Земли вокруг своей оси.

Поправки за влияние тропосферы вычисляются одним из существующих методов (например, методы Саастомойнена и Хопфилда). При длине измеряемого базиса не более 5 км введение поправок за тропосферу может быть пропущено.

Поправки за влияние ионосферы могут вычисляться с использованием двухчастотных синхрон* ных измерений кодовых и фазовых псевдодальностей. Поправки за влияние ионосферы для фазовых ПД численно равны поправкам, полученным для кодовых ПД. но знаки их противоположны (поправки к кодовым ПД отрицательны, а поправки к фазовым ПД положительны). 8 случае измерений кодовых ПД формулы для вычисления этих поправок соответственно на первой и второй частотах имеют вид:

Ad

Ad

(d,-d2).

(Ф - d2).

d>

<2>

где /,, /j — значения первой и второй несущих частот; of,. d2 — значения измеренных кодовых ПД соответственно на первой и второй частотах. Поправка на вращение Земли вокруг своей оси за время распространения сигнала от КА до пун* кта наблюдения (в метрах) вычисляется по формуле:

<оз_

д<*»р =--р sinA cosB sinZ. (3)

где to*- угловая скорость вращения Земли вокруг своей оси; ас — большая полуось общеэемного эллипсоида; р — расстояние до КА (в метрах); с — скорость света, равная 299792498 м/с; в— широта пункта наблюдения;

A. Z — соответственно азимут и зенитное расстояние КА.

4.2 Определение первых разностей фаз Исходным соотношением служит формула:

(4}

где; — 1.2 —

номер спутника;

номер момента времени измерения /•; номера исходного и определяемого пунктов; величина ПРФ; измеренные значения фаз.

2

ГОСТ Р 53608—2009

С использованием ряда допущений, практически не влияющих на точность решения, величину S(j. i) представляют в виде

0i(/. i) = /{Й/'( + 8f' - t>2. + tjv1 + N{2 + е/.

N' - N' - N>. "12 "2 ,vi *

(«)

(6)

где/— несущая частота;

# j и т*._ поправка к шкале времени и интервал времени прохождения сигналов от/-го спутника до

т о*го пункта в момент f, (п = 1.2);

WJ — целочисленный параметр ПРФ при / = 1;

е/_случайная погрешность измерения первой разности фаз.

Формула (5) нелинейна относительно поправки ЙГ, в том смысле, что параметр т' (0 должен быть получен на истинный момент приема сигнала (f, + S,}.

Для существенного сокращения числа искомых параметров и повышения точности относитель-ных определений поправки к временным шкалам пунктов представляют в виде:

V = Vi (?)

где 6,f, — поправка к временной шкале М-ro пункта в момент первого измерения; са| — коэффициент линейного ухода временной шкалы М-ro пункта;

М — номер пункта (М = 1.2);

д/, — промежуток времени между текущим и начальным моментом измерений.

4.3 Определение вторых разностей фаз Исходным соотношением служит формула

где £>(/, к, i) — вторая разность фаз на момент (.

С учетом полученных первых разностей фаз вторые разности фаз определяют по формуле

02</,М = 'К -й/‘.)-/с уъ + (ЙГ*2 -p*c1)J + fc 1[_p*v +p*v(8f*b -р*.с 1)J- W

-6(£) + fc ^ p'. + p^(8f' -p^.c ’)] + /c I p^ + p'(8*' -p'.c !>] + ^2 +e*;

N14 */V* -N' ,v12 fV12 '*12'

<t0)

где N^2 — целочисленный параметр ВРФ. p'. и p*. — мгновенные дальности /-го и к-го спутников от л-го пункта в момент t,\

Pi и р* — скорость изменения дальности /• го и к-то спутников от л-го пункта в момент f,;

/— частота несущей; с— скорость света;

е*/ х случайная погрешность вторых разностей фаз в момент I. ■ •

Линеаризация позволяет получить уравнение погрешностей ВРФ. в котором неизвестными являются следующие параметры:

• поправки к предварительным координатам двух пунктов (шесть параметров);

• временные параметры двух пунктов (четыре неизвестных);

• целочисленные параметры, число которых на единицу меньше числа наблюденных спутников.

Систему уравнений погрешностей ВРФ на моменты f,.....tN описывают матрицей коэффициен

тов и вектором свободных членов.

Матрица коэффициентов представляют в виде

ч

1

<11)

А.

где N — число моментов измерений ВРФ;

Ау.....Аы — блоки идентичной структуры, зависящие от момента измерений и числа наблюдаемых

спутников.

3

ГОСТ Р 53608—2009

Например, для четырех спутников каждый из блоков Ал.....Аы определяется по формуле

®2г

а% бз> Сэ, of3i «з* fy

'V

е

4*

9ъ ^2. Ръ 9а

0з» h з, рэ, <73,

9*. ^4i Pai 9*

1 О О О 1 О 0 0 1

/8 1..

. N,

(12)

где ак1 и ckl — коэффициенты при поправках к начальным моментам измерений первого и второго пунктов;

bkj и dkl — коэффициенты при скоростях ухода шкал времени первого и второго пунктов; e*j’ Кг &к>и hk>' Р*у 9kl — коэффициенты при поправках к предварительным координатам первого и

второго пунктов.

Коэффициенты akj и ск-вычисляются по формулам

а*. =

Коэффициенты Ьк1 и dkl вычисляются по формулам:

-f;x

- ‘W „

dwS^riaftr1ft,,xpi -р*х

Коэффициенты eki, fki, дк1 вычисляются по формулам:

<ЧО)

ДХ,*

ек, ~

гЧХ|)

Pi*. ’

Г. -

*£»

ДУ*

•к)

Р1

П. 2

30)

- А2/

3Z,*

ski

р/

р?

Коэффициенты hki, pkj. qkj вычисляются no формулам:

- <ХО)

_ ДХ2/

. лх!

Pjf

Pzt

.

-AV

ДУ/

’ ТО)

Ра

Р2

- *0)

- А#

AZj

W r,' P$

(13)

(14)

(15)

(16)

Коэффициенты при целочисленных параметрах ВРФ образуют единичную матрицу.

Вектор свободных членов системы уравнений погрешностей ВРФ определяется как вектор, элементы которого представляют собой разности ВРФ. вычисляемых по предварительным значениям

определяемых параметров, и измеренных значений ВРФ на моменты.....tN.

Полученную систему уравнений поправок решают по МНК с учетом корреляции ВРФ.

4.4 Учет корреляции вторых разностей фаз

Корреляцию ВРФ учитывают путем использования в решении ковариационной матрицы ВРФ О. Данная матрица имеет блочно-диагональную структуру. Диагональные блоки являются квадратными матрицами. Их число равно числу N моментов измерений. Размер блока выражается числом, на едини

ГОСТ Р 53608—2009

цу меньшим числа М наблюдаемых спутников. Например, при М = 4 каждый из указанных блоков имеет вид

г4 2 21 2 4 2 ,2 2 4

J- 1.....N.

4.5 Получение целочисленных параметров вторых разностей фаз

Определяемые параметры ВРФ получают из решения системы нормальных уравнений.

Матрицу системы нормальных уравнений формируют в соответствии с формулой

В - АТРА, (18)

где Ат — матрица, транспонированная по отношению к матрице А;

Р — весовая матрица, обратная по отношению к ковариационной матрице О.

Вектор свободных членов системы нормальных уравнений формируют по формуле:

C = ATPL. (18)

где L — вектор свободных членов уравнений погрешностей ВРФ.

Из решения системы нормальных уравнений получают:

• временные параметры:

• поправки к координатам исходного и определяемого пунктов;

• целочисленные параметры ВРФ.

Вследствие погрешностей измерений, плохих условий наблюдений навигационных спутников (многолучевость, искусственные и естественные препятствия наблюдений и т. д.) фактические значения целочисленных параметров, как правило, не являются целыми числами. Поэтому полученные значения этих параметров должны быть округлены до целого в ближайшую сторону. С этими новыми значениями необходимо перевычислить вектор свободных членов уравнений погрешностей и получить новое решение с учетом поправок к шкалам времени пунктов и скоростям их изменения, а также поправок к координатам пунктов.

5

ГОСТ Р 53608—2009

УДК 629.783:1528.2+528.344+523.34.13]:006.354 ОКС 07.040 Э50 ОКСТУ 6801

Ключевые слова: глобальные навигационные спутниковые системы, геодезические и землеустроительные работы, неоднозначность, фазовая псеадодальность

6

Редактор £.С. Котлярова Технический редактор В.Н. Прусакова Корректор М.В. Буйная Компьютерная верстка Л.А. Круговой

Сдано е набор 12.11.2010. Подписаноа печать 14.12.2010. Формат 60 ■ 64у<£ Бумага офсетная. Гарнитура Ариял. Печать офсетная. Уеп. печ. я. 1.40. Уч.-иэд. я. 0.50. Тираж 89 экэ. Зак. 1030.

. 12399S Москва. Гранатный пер.. 4. info^gostinfo.iu

Набрано во на ПЭВМ.

Отпечатано о филиале — тип. «Московский печатник». 105062 Москва. Лялин пер., в.