ГОСТ Р 8.1036-2024. Национальный стандарт Российской Федерации. Государственная система обеспечения единства измерений. Стандарты частоты и (или) времени. Методы оценки основных метрологических характеристик с использованием сигналов, передаваемых глобальными навигационными спутниковыми системами
Актуально на 02.11.2026.
Как цитировать: "ГОСТ Р 8.1036-2024. Национальный стандарт Российской Федерации. Государственная система обеспечения единства измерений. Стандарты частоты и (или) времени. Методы оценки основных метрологических характеристик с использованием сигналов, передаваемых глобальными навигационными спутниковыми системами". ГОСТ Р 8.1036-2024. https://bukva.ai/docs/standart-gost-r-8-1036-2024
(утв. и введен в действие Приказом Росстандарта от 12.04.2024 № 447-ст)
Примечание
Документ введен в действие с 30.09.2024.
Утвержден и введен в действие Приказом Федерального агентства по техническому регулированию и метрологии от 12 апреля 2024 г. № 447-ст
НАЦИОНАЛЬНЫЙ СТАНДАРТ РОССИЙСКОЙ ФЕДЕРАЦИИ
ГОСУДАРСТВЕННАЯ СИСТЕМА ОБЕСПЕЧЕНИЯ ЕДИНСТВА ИЗМЕРЕНИЙ
СТАНДАРТЫ ЧАСТОТЫ И (ИЛИ) ВРЕМЕНИ
МЕТОДЫ ОЦЕНКИ ОСНОВНЫХ МЕТРОЛОГИЧЕСКИХ ХАРАКТЕРИСТИК С ИСПОЛЬЗОВАНИЕМ СИГНАЛОВ, ПЕРЕДАВАЕМЫХ ГЛОБАЛЬНЫМИ НАВИГАЦИОННЫМИ СПУТНИКОВЫМИ СИСТЕМАМИ
State system for ensuring the uniformity of measurements. Frequency and/or time standards. Methods for assessing basic metrological characteristics using signals transmitted by global navigation satellite
Дата введения 30 сентября 2024 года
Предисловие
1 РАЗРАБОТАН Федеральным государственным унитарным предприятием «Всероссийский научно-исследовательский институт физико-технических и радиотехнических измерений» (ФГУП «ВНИИФТРИ»)
2 ВНЕСЕН Техническим комитетом по стандартизации ТК 206 «Эталоны и поверочные схемы», подкомитетом ПК 12 «Эталоны и поверочные схемы в области измерений времени и частоты»
3 УТВЕРЖДЕН И ВВЕДЕН В ДЕЙСТВИЕ Приказом Федерального агентства по техническому регулированию и метрологии от 12 апреля 2024 г. № 447-ст
Правила применения настоящего стандарта установлены в статье 26 Федерального закона от 29 июня 2015 г. № 162-ФЗ «О стандартизации в Российской Федерации». Информация об изменениях к настоящему стандарту публикуется в ежегодном (по состоянию на 1 января текущего года) информационном указателе «Национальные стандарты», а официальный текст изменений и поправок — в ежемесячном информационном указателе «Национальные стандарты». В случае пересмотра (замены) или отмены настоящего стандарта соответствующее уведомление будет опубликовано в ближайшем выпуске ежемесячного информационного указателя «Национальные стандарты». Соответствующая информация, уведомление и тексты размещаются также в информационной системе общего пользования — на официальном сайте Федерального агентства по техническому регулированию и метрологии в сети Интернет (www.rst.gov.ru)
Введение
В современных науке и производстве, в таких высокотехнологических видах жизнедеятельности, как связь, транспорт, энергетика, медицина, сфера обеспечения обороноспособности и безопасности государства и другие, широко применяются стандарты частоты и (или) времени или прецизионные средства измерений на их основе.
В процессе эксплуатации таких средств измерений возникает необходимость оценки их характеристик на местах эксплуатации, что на практике, как правило, трудно реализуемо. В то же время развитие оборудования приема сигналов глобальных навигационных спутниковых систем и методов обработки полученных результатов измерений позволяет произвести оценку основных метрологических характеристик стандартов частоты и (или) времени без доставки к местам их размещения стандартов, имеющих более высокие показатели точности.
Данный стандарт направлен на развитие методов контроля характеристик стандартов частоты и (или) времени, снижение рисков выведения из строя высокоточного оборудования при транспортировке к месту проведения работ, а также на сокращение расходов на их содержание.
Объектом стандартизации являются методы оценки основных метрологических характеристик стандартов частоты и (или) времени с использованием сигналов, передаваемых глобальными навигационными спутниковыми системами.
1 Область применения
Настоящий стандарт распространяется на стандарты частоты и (или) времени, относительная погрешность меры частоты которых находится в пределах +/- 1,0·10-9, что соответствует уровню рабочих эталонов первого — третьего разрядов по государственной поверочной схеме для средств измерений времени и частоты <1>, и устанавливает методы оценки их основных метрологических характеристик с использованием сигналов, передаваемых глобальными навигационными спутниковыми системами.
<1> Приказ Федерального агентства по техническому регулированию и метрологии от 26 сентября 2022 г. № 2360 «Об утверждении государственной поверочной схемы для средств измерений времени и частоты».
2 Нормативные ссылки
В настоящем стандарте использованы нормативные ссылки на следующие стандарты:
ГОСТ 8.567 Государственная система обеспечения единства измерений. Измерения времени и частоты. Термины и определения
ГОСТ Р 8.736 Государственная система обеспечения единства измерений. Измерения прямые многократные. Методы обработки результатов измерений. Основные положения
Примечание - При пользовании настоящим стандартом целесообразно проверить действие ссылочных стандартов в информационной системе общего пользования — на официальном сайте Федерального агентства по техническому регулированию и метрологии в сети Интернет или по ежегодному информационному указателю «Национальные стандарты», который опубликован по состоянию на 1 января текущего года, и по выпускам ежемесячного информационного указателя «Национальные стандарты» за текущий год. Если заменен ссылочный стандарт, на который дана недатированная ссылка, то рекомендуется использовать действующую версию этого стандарта с учетом всех внесенных в данную версию изменений. Если заменен ссылочный стандарт, на который дана датированная ссылка, то рекомендуется использовать версию этого стандарта с указанным выше годом утверждения (принятия). Если после утверждения настоящего стандарта в ссылочный стандарт, на который дана датированная ссылка, внесено изменение, затрагивающее положение, на которое дана ссылка, то это положение рекомендуется применять без учета данного изменения. Если ссылочный стандарт отменен без замены, то положение, в котором дана ссылка на него, рекомендуется применять в части, не затрагивающей эту ссылку.
3 Термины и определения
В настоящем стандарте применены термины по ГОСТ 8.567.
4 Сокращения
В настоящем стандарте применены следующие сокращения:
ГНСС — глобальная навигационная спутниковая система (GPS, ГЛОНАСС, Galileo, Beidou и другие);
ГСВЧ — государственная служба времени, частоты и определения параметров вращения Земли;
ГЭТ 1 — государственный первичный эталон единиц времени, частоты и национальной шкалы времени;
МБМВ — международное бюро мер и весов;
НСП — неисключенная систематическая погрешность;
СКО — среднее квадратическое относительное отклонение;
СКДО — среднее квадратическое относительное двухвыборочное отклонение;
СПО — специальное программное обеспечение;
СЧВ — стандарты частоты и (или) времени;
ШВС — шкала времени системная ГНСС;
ЭВИ — эфемеридно-временная информация;
UTC(SU) — национальная шкала координированного времени (национальная шкала времени);
UTC(SU)ГЛО — реализация национальной шкалы времени UTC(SU), передаваемая системой ГЛОНАСС;
UTC(SU)СИ — локальная реализация национальной шкалы времени UTC(SU), формируемая эталоном, средством измерений или техническим устройством;
5 Методы оценки основных метрологических характеристик стандартов частоты и (или) времени с использованием сигналов, передаваемых глобальными навигационными спутниковыми системами
5.1 Методы оценки смещений (разности) ШВ СЧВ
В соответствии с требованиями федерального закона [1], а также постановления [2] ГСВЧ распространяет информацию о точном значении московского времени и календарной дате. Информация о точном значении московского времени и календарной дате формируется на основе UTC(SU), является официальной, общедоступной и обязательной для использования в Российской Федерации. Указанная информация, а также эталонные сигналы частоты и времени распространяются с использованием ГНСС ГЛОНАСС, спутниковых систем связи, радиосвязи и радиовещания.
СИ, применяемые в сфере государственного регулирования обеспечения единства измерений и обеспечивающие контроль погрешности синхронизации времени, должны прослеживаться к ГЭТ 1 с установлением обязательных метрологических требований к смещениям (разности) местных ШВ относительно UTC(SU).
При оценке смещений (разности) ШВ СЧВ следует применять поверенные приемники сигналов ГНСС, прошедшие процедуру калибровки внутренних задержек сигналов (задержек сигналов в приемном тракте) и имеющие соответствующие сертификаты с указанием значений погрешности (неопределенности) измерений при калибровке.
Основными методами сравнения ШВ по сигналам ГНСС являются метод «односторонней передачи» (абсолютный режим) и метод «одновременного наблюдения» (дифференциальный режим), описанные в [3], [4].
ГНСС являются сложными многокомпонентными системами, и погрешность координатно-временных определений с их применением зависит от большого количества факторов.
Основными источниками погрешности являются:
- погрешность ЭВИ;
- атмосферная задержка сигнала ГНСС (в ионосфере и тропосфере);
- условия приема сигнала ГНСС: блокирование искусственными и природными препятствиями, переотражения сигнала, радиочастотные помехи и др.;
- аппаратурные погрешности (калибровка ГНСС приемника, шум измерений и др.).
Погрешность, вносимую большинством из этих факторов, невозможно заранее спрогнозировать, поэтому она трактуется как случайная и оценивается статистическими методами. Для достоверной оценки случайной составляющей погрешности интервалы времени наблюдений должны составлять не менее одних суток, при этом интервалы времени наблюдений и интервалы времени измерения рекомендуется выбирать в соответствии с ГОСТ 8.567.
Границы (доверительные границы) случайной погрешности оценки смещений ШВ и суммарной погрешности измерений оцениваются в соответствии с ГОСТ Р 8.736. Составляющие НСП измерений более подробно рассмотрены в 5.1.1 - 5.1.3, 5.2.1 и 5.2.2.
5.1.1 Методы оценки смещений (разности) ШВ с применением приемников сигналов ГНСС, оснащенных выходом импульсных сигналов 1 Гц
Для реализации данного метода применяется один приемник сигналов ГНСС, размещаемый в месте дислокации исследуемого СЧВ. Обязательным условием в данном случае является наличие у приемника нормированной метрологической характеристики, характеризующей пределы допускаемых смещений формируемой ШВ (выходного импульсного сигнала 1 Гц) относительно UTC(SU).
Измерения выполняют по схеме, приведенной на рисунке А.1.
Для удобства обработки результатов измерений рекомендуется обеспечить их автоматическую фиксацию с применением СПО, считывающего показания измерителя временных интервалов (двухканального частотомера, работающего в режиме измерения интервалов времени) через имеющийся интерфейс (RS-232, USB-порт, LAN и т.д.).
Среднее значение смещений ШВ СЧВ относительно UTC(SU)СИ, приведенное к середине интервала времени наблюдений \Delta T_{[UTC(SU)^{\text{СИ}} - \text{СЧВ}]}(t_{i}) , вычисляют по формуле
\Delta T_{[UTC(SU)^{\text{СИ}} - \text{СЧВ}]}(t_{i}) = \Delta T\prime_{[UTC(SU)^{\text{СИ}} - \text{СЧВ}]}(t_{i}) + \tau_{K1} - \tau_{K2} , (1)
где \Delta T\prime_{[UTC(SU)^{\text{СИ}} - \text{СЧВ}]}(t_{i}) — показание измерителя интервалов времени в момент времени ti;
\tau_{K1} — значение задержки импульсного сигнала 1 Гц в кабеле № 1;
\tau_{K2} — значение задержки импульсного сигнала 1 Гц в кабеле № 2.
Доверительные границы НСП оценки смещений ШВ \Theta_{\Sigma}(P) в данном случае оценивают в соответствии с ГОСТ Р 8.736 и вычисляют по формуле
\Theta_{\Sigma}(P) = \pm k \cdot \sqrt{\sum_{i = 1}^{5} \Theta_{i}^{2}} , (2)
где k — коэффициент, определяемый принятой доверительной вероятностью;
\Theta_{1} — составляющая НСП, обусловленная пределами допускаемых смещений формируемой приемником сигналов ГНСС ШВ относительно UTC(SU);
\Theta_{2} — составляющая НСП, обусловленная пределами погрешности измерителя интервалов времени;
\Theta_{3} — составляющая НСП, обусловленная пределами погрешности определения координат антенны приемника сигналов ГНСС;
\Theta_{4} , \Theta_{5} — составляющие НСП, обусловленные пределами погрешности измерения задержки импульсного сигнала 1 Гц в кабелях № 1 и № 2.
5.1.2 Методы оценки смещений (разности) ШВ с использованием псевдодальномерных кодовых измерений
Для оценки разности ШВ по кодовым псевдодальномерным измерениям рекомендуется использование данных в формате CGGTTS или аналогичном. Файлы в формате CGGTTS содержат оценки расхождения ШВ, являющейся опорной ШВ для приемника и ШВС ГНСС по каждому видимому КА, полученные по различным типам сигналов. При расчете расхождений ШВ используют только данные непосредственных кодовых измерений псевдодальности и данные навигационного сообщения о разности между ШВС ГНСС и опорными ШВ [UTC(SU), UTC(USNO) и т.д.].
Также расхождение ШВ может быть рассчитано с использованием псевдодальномерных кодовых ГНСС измерений, содержащихся в файлах в формате RINEX [5] или подобном. Алгоритмы расчетов подробно описаны в [6] - [8], также существует несколько вариантов СПО, предназначенного для проведения расчетов по описанным алгоритмам, распространяемого как на коммерческой основе, так и свободно. Например, свободно распространяемое МБМВ ПО с открытым исходным кодом на языке Fortran под названием «R2CGGTTS».
Формат файлов CGGTTS стандартизован МБМВ и описан в [9], [10], специализированные ГНСС приемники для временных сличений, как правило, штатно формируют файлы в данном формате.
5.1.2.1 Абсолютный режим
При выполнении измерений применяется один приемник сигналов ГНСС, размещаемый в месте дислокации исследуемого СЧВ. В качестве опорных сигналов для приемника используются выходные сигналы исследуемого СЧВ: импульсный сигнал 1 Гц и высокочастотный синусоидальный сигнал (как правило, 5 или 10 МГц).
Измерения выполняют по схеме, приведенной на рисунке Б.1. Значения задержек сигналов в приемном тракте приемника и задержки импульсного сигнала 1 Гц в кабеле № 2, как правило, заносятся в файлы настроек приемников и учитываются при обработке результатов измерений автоматически.
Структура формата файлов CGGTTS и ему подобных допускает дискретизацию данных с интервалами времени измерения 5 и 16 мин, а также кратных этим значениям.
Значение смещения ШВ СЧВ относительно реализации национальной шкалы времени UTC(SU), передаваемой системой ГЛОНАСС, \Delta T_{[UTC(SU)^{\text{гло}} - \text{СЧВ}]^{s}}(t_{i}) вычисляют по формуле
\Delta T_{[UTC(SU)^{\text{гло}} - \text{СЧВ}]^{s}}(t_{i}) = \frac{\sum_{j = 1}^{M} \cdot T_{s}^{j}(t_{i})}{M} + \tau_{c}, (3)
где \Delta T_{s} — смещение ШВ СЧВ относительно UTC(SU), передаваемой системой ГЛОНАСС, по типу сигнала s в момент времени ti;
M — общее количество наблюдаемых КА в момент времени ti;
\Delta T_{s}^{j} — расхождение ШВ приемника со ШВ ГЛОНАСС по наблюдениям j-го КА, типу сигнала s (столбец REFSYS файла CGGTTS);
\tau_{c} — поправка к ШВ системы ГЛОНАСС относительно UTC(SU), передаваемая в кадре навигационного сообщения [11].
Доверительные границы НСП оценки смещений ШВ \Theta_{\Sigma}(P) в данном случае оценивают в соответствии с ГОСТ Р 8.736 и вычисляют по формуле
\Theta_{\Sigma}(P) = \pm k \cdot \sqrt{\sum_{i = 1}^{4} \Theta_{i}^{2}} , (4)
где k — коэффициент, определяемый принятой доверительной вероятностью;
\Theta_{1} — составляющая НСП, обусловленная погрешностью формирования информации о поправке к ШВ системы ГЛОНАСС относительно UTC(SU), передаваемой в кадре навигационного сообщения;
\Theta_{2} — составляющая НСП, обусловленная погрешностью измерений при калибровке приемника сигналов ГНСС;
\Theta_{3} — составляющая НСП, обусловленная пределами погрешности определения координат антенн приемника сигналов ГНСС;
\Theta_{4} — составляющая НСП, обусловленная пределами погрешности измерения задержки импульсного сигнала 1 Гц в кабеле № 2.
5.1.2.2 Дифференциальный режим
Дифференциальный режим позволяет существенно уменьшить составляющие погрешности за счет погрешности формирования ЭВИ, атмосферной задержки сигналов и нестабильности ШВС ГНСС.
Для реализации дифференциального режима при измерениях применяется два аналогичных приемника сигналов ГНСС, размещаемых в местах дислокации исследуемого СЧВ и СЧВ, имеющего более высокие показатели точности (далее — опорный СЧВ).
Для оценки смещений (разности) ШВ следует использовать только те типы сигналов ГНСС, для которых в процессе калибровки определены задержки сигналов в приемном тракте обоих приемников.
Измерения выполняют одновременно по схеме, приведенной на рисунке В.1. Для оценки характеристик используется измерительная информация двух приемников, для чего должен быть обеспечен информационный обмен.
Обязательным условием в данном случае является наличие у опорного СЧВ, нормированной метрологической характеристики, характеризующей пределы допускаемых смещений формируемой ШВ (выходного импульсного сигнала 1 Гц) относительно UTC(SU).
В зависимости от длины базовой линии (расстояния между исследуемым СЧВ и опорным СЧВ) могут применяться различные методы оценки разности ШВ.
Метод усреднения по одновременно наблюдаемым КА (метод Common-View)
Данный метод является наиболее предпочтительным при сравнении удаленных ШВ, поскольку позволяет существенно уменьшить составляющие погрешности, обусловленные погрешностью формирования ЭВИ и нестабильностью ШВС ГНСС. Недостатком данного метода является ограничение по расстоянию между исследуемым СЧВ и опорным СЧВ, поскольку с увеличением расстояния уменьшается количество одновременно наблюдаемых КА. Рекомендуемая длина базовых линий — до 2000 км.
При использовании метода Common-View значение смещения ШВ исследуемого СЧВ относительно ШВ опорного СЧВ \Delta T_{s}(t_{i}) вычисляют по формуле
\Delta T_{s}(t_{i}) = \frac{\sum_{j = 1}^{M} \Delta T_{\text{СЧВ},s}^{j}(t_{i}) - \Delta T_{\text{СЧВэ},s}^{j}(t_{i})}{M} , (5)
где \Delta T_{\text{СЧВ},s}^{j}(t_{i}) — смещение ШВ исследуемого СЧВ относительно ШВС ГНСС, передаваемой j-м КА, по типу сигнала s в момент времени ti;
\Delta T_{\text{СЧВэ},s}^{j}(t_{i}) — смещение ШВ опорного СЧВ относительно ШВС ГНСС, передаваемой j-м КА, по типу сигнала s в момент времени ti;
M — количество одновременно наблюдаемых КА.
Доверительные границы НСП оценки смещений ШВ исследуемого СЧВ относительно UTC(SU) \Theta_{\Sigma}(P) в данном случае оценивают в соответствии с ГОСТ Р 8.736 и вычисляют по формуле
\Theta_{\Sigma}(P) = \pm k \cdot \sqrt{\sum_{i = 1}^{7} \Theta_{i}^{2}} , (6)
где k — коэффициент, определяемый принятой доверительной вероятностью;
\Theta_{1} , \Theta_{2} — составляющие НСП, обусловленные погрешностью измерений при калибровке приемников сигналов ГНСС № 1 и № 2 соответственно;
\Theta_{3} , \Theta_{4} — составляющие НСП, обусловленные пределами погрешности определения координат антенн приемников сигналов ГНСС № 1 и № 2 соответственно;
\Theta_{5} , \Theta_{6} — составляющие НСП, обусловленные пределами погрешности измерения задержек импульсных сигналов 1 Гц в кабелях № 2 и № 3 соответственно;
\Theta_{7} — составляющая НСП, обусловленная пределами допускаемых смещений ШВ опорного СЧВ относительно UTC(SU).
1 В соответствии с [11] погрешность формирования информации о поправке к ШВ системы ГЛОНАСС относительно UTC(SU), передаваемой в кадре навигационного сообщения, находится в пределах +/- 1 мкс.
2 Если приемник способен принимать сигналы ГНСС в нескольких частотных диапазонах, можно практически полностью исключить ионосферную составляющую погрешности за счет формирования безыоносферной комбинации сигналов, однако в этом случае увеличивается случайная аппаратурная составляющая погрешности кодовых измерений.
3 Если оба приемника способны принимать сигналы ГНСС в нескольких частотных диапазонах, можно выбрать оптимальный тип сигнала в зависимости от длины базовой линии; например, для расстояний между приемниками менее 100 км рекомендуется использовать одночастотные измерения, при расстоянии более 300 км — безыоносферную комбинацию.
4 Для небольших базовых линий путь, который проходит сигнал ГНСС в атмосфере, достаточно схож, и воздействие атмосферы существенно уменьшается при формировании разности по одновременно наблюдаемым КА. При этом нескомпенсированная часть ионосферной задержки оказывается меньше, чем увеличение шума измерений из-за формирования безыоносферной комбинации.
Метод усреднения по всем наблюдаемым на каждом пункте КА (метод All-In-View)
Данный метод применяют при отсутствии или слишком малом количестве одновременно наблюдаемых КА. Рекомендуемая длина базовых линий — более 3000 км.
При использовании метода All-In-View значение смещения ШВ исследуемого СЧВ относительно ШВ опорного СЧВ \Delta T_{s}(t_{i}) вычисляют по формуле
\Delta T_{s}(t_{i}) = \frac{\sum_{j = 1}^{N} \Delta T_{\text{СЧВ},s}^{j}(t_{i})}{N} - \frac{\sum_{k = 1}^{M} \Delta T_{\text{СЧВэ},s}^{k}(t_{i})}{M} , (7)
где \Delta T_{\text{СЧВ},s}^{j} — смещение ШВ исследуемого СЧВ относительно ШВС ГНСС, передаваемой j-м КА, по типу сигнала s в момент времени ti;
\Delta T_{\text{СЧВэ},s}^{k} — смещение ШВ опорного СЧВ относительно ШВС ГНСС, передаваемой k-м КА, по типу сигнала s в момент времени ti;
N — количество КА, наблюдаемых в месте дислокации исследуемого СЧВ;
M — количество КА, наблюдаемых в месте дислокации опорного СЧВ.
Доверительные границы НСП оценки смещений ШВ исследуемого СЧВ относительно UTC(SU) \Theta_{\Sigma}(P) в данном случае оценивают в соответствии с ГОСТ Р 8.736 и вычисляют по формуле (6).
5.1.3 Методы оценки смещений (разности) ШВ с использованием псевдодальномерных кодовых и фазовых измерений
Если имеющиеся ГНСС приемники способны принимать кодовые и фазовые измерения как минимум в двух частотных диапазонах, для оценки смещений могут быть применены методы апостериорного высокоточного абсолютного позиционирования (англ. PPP — Precise Point Positioning). Существует большое количество алгоритмов высокоточного абсолютного позиционирования, отличающихся используемыми моделями воздействующих факторов, способом разрешения неоднозначности фазовых псевдодальномерных измерений, составом уточняемых параметров и другими нюансами (далее — PPP).
Основные принципы метода PPP и применяемые модели описаны в работах [12], [13], варианты алгоритмов с целочисленным разрешением фазовой неоднозначности — в работах [14], [15]. Существуют как коммерческие, так к свободно распространяемые реализации метода PPP с открытым исходным кодом (например, RTKLib, GPSTk, PPPLib и другие). Также есть ряд публично доступных онлайн-сервисов, которые позволяют любому желающему получить PPP решение, загрузив файл с псевдодальномерными измерениями в формате RINEX (например, NRCAN PPP, GAPS, APPS, magicPPP).
Общим для данных алгоритмов является использование в качестве исходной измерительной информации безыоносферных комбинаций кодовых и фазовых измерений, а также применение в решении высокоточных апостериорных орбит КА и поправок к ШВ КА. Это позволяет практически полностью исключить влияние ионосферы и космического сегмента ГНСС на точность получаемых оценок. В ходе решения одновременно с разностью ШВ также уточняются координаты фазового центра ГНСС антенны и другие параметры. В качестве высокоточных эфемерид могут применяться публично доступные данные центров анализа ГНСС информации или сводные продукты международной службы ГНСС (англ. IGS — International GNSS Service). В зависимости от точности оперативность продуктов IGS составляет от нескольких часов до двух недель после окончания соответствующих суток UTC.
Метод PPP аналогичен методу All-In-View и позволяет получить оценки расхождения ШВ двух пунктов вне зависимости от расстояния между ними. По аналогии с дифференциальным режимом сравнения ШВ по кодовым измерениям применяются два аналогичных приемника сигналов ГНСС, размещаемые в местах дислокации исследуемого СЧВ и опорного СЧВ. Измерения выполняют одновременно по схеме, приведенной на рисунке В.1. Для оценки характеристик используется измерительная информация двух приемников, для чего должен быть обеспечен информационный обмен. Независимо обрабатывается измерительная информация двух ГНСС приемников, после чего вычисляется разность их ШВ \Delta T(t_{i}) по следующей формуле:
\Delta T(t_{i}) = \Delta T_{\text{СЧВ}}(t_{i}) - \Delta T_{\text{СЧВэ}}(t_{i}) , (8)
где \Delta T_{\text{СЧВ}}(t_{i}) — смещение ШВ исследуемого СЧВ относительно опорной ШВ высокоточных эфемерид, оцененное с помощью PPP в момент времени ti;
\Delta T_{\text{СЧВэ}}(t_{i}) — смещение ШВ опорного СЧВ относительно опорной ШВ высокоточных эфемерид, оцененное с помощью PPP в момент времени ti.
Доверительные границы НСП оценки смещений ШВ исследуемого СЧВ относительно UTC(SU) \Theta_{\Sigma}(P) в данном случае оценивают в соответствии с ГОСТ Р 8.736 и вычисляют по формуле
\Theta_{\Sigma}(P) = \pm k \cdot \sqrt{\sum_{i = 1}^{5} \Theta_{i}^{2}} , (9)
где k — коэффициент, определяемый принятой доверительной вероятностью;
\Theta_{1} , \Theta_{2} — составляющие НСП, обусловленные погрешностью измерений при калибровке приемников сигналов ГНСС № 1 и № 2 соответственно;
\Theta_{3} , \Theta_{4} — составляющие НСП, обусловленные пределами погрешности измерения задержек импульсных сигналов 1 Гц в кабелях № 2 и № 3 соответственно;
\Theta_{5} — составляющая НСП, обусловленная пределами допускаемых смещений ШВ опорного СЧВ относительно UTC(SU).
5.2 Методы оценки относительной погрешности меры частоты (относительной разности частот)
При оценке относительной разности частот не предъявляются требования к калибровке приемника сигналов ГНСС, также нет необходимости измерять задержку импульсного сигнала 1 Гц до входа в приемник или в измеритель интервалов времени.
5.2.1 Методы оценки относительной погрешности меры частоты (относительной разности частот) с применением приемников сигналов ГНСС, оснащенных выходом высокочастотных синусоидальных сигналов
Для реализации данного метода используются приемники сигналов ГНСС, оснащенные выходом высокочастотных синусоидальных сигналов со стандартными действительными значениями частоты (1, 5, 10 МГц и т.д.). При выполнении измерений применяется один приемник сигналов ГНСС, размещаемый в месте дислокации исследуемого СЧВ.
Обязательным условием в данном случае является наличие у приемника нормированной метрологической характеристики, характеризующей пределы допускаемой относительной погрешности меры частоты (относительной разности частот) выходных высокочастотных синусоидальных сигналов.
Измерения выполняют по схеме, приведенной на рисунке Г.1. Как правило, частотные (фазовые) компараторы поставляются в комплекте с СПО, обеспечивающим оценку различных частотных характеристик, в том числе и относительной разности частот \Delta_{0}f_{j} , вычисляемой по формуле
\Delta_{0}f_{j} = (\frac{\Delta f}{f})\text{изм} , (10)
где (\frac{\Delta f}{f})\text{изм} — результат измерения относительной разности частот.
Доверительные границы НСП оценки относительной разности частот \Theta_{\Sigma}(P) в данном случае определяют характеристиками приемника сигналов ГНСС и частотного (фазового) компаратора, оценивают в соответствии с ГОСТ Р 8.736 и вычисляют по формуле
\Theta_{\Sigma}(P) = \pm k \cdot \sqrt{\sum_{i = 1}^{2} \Theta_{i}^{2}} , (11)
где k — коэффициент, определяемый принятой доверительной вероятностью;
\Theta_{1} — составляющая НСП, обусловленная пределами погрешности частотного (фазового) компаратора;
\Theta_{2} — составляющая НСП, обусловленная пределами допускаемой относительной погрешности меры частоты (относительной разности частот) выходных высокочастотных синусоидальных сигналов приемника сигналов ГНСС.
5.2.2 Методы оценки относительной погрешности меры частоты (относительной разности частот) на основе оценок разности ШВ
Исходными данными для оценки относительной погрешности меры частоты (относительной разности частот) являются значения смещений (разности) ШВ, вычисленные по формулам (3), (5), (7) или (8).
Значение относительной разности частот исследуемого СЧВ и опорного СЧВ \Delta_{0}f_{j} на интервале времени между конечным и начальным моментами времени tк и tн вычисляют по формуле
\Delta_{0}f_{j} = \frac{\Delta T(t_{\text{к}}) - \Delta T(t_{\text{н}})}{\tau_{\text{и}}} , (12)
где \Delta T(t_{\text{к}}) — смещение ШВ исследуемого СЧВ относительно опорной ШВ в момент времени tк, вычисленное по формулам (3), (5), (7) или (8);
\Delta T(t_{\text{н}}) — смещение ШВ исследуемого СЧВ относительно опорной ШВ в момент времени tн, вычисленное по формулам (3), (5), (7) или (8);
\tau_{\text{и}} — интервал времени измерения, вычисляемый по формуле \tau_{\text{и}} = \tau_{\text{к}} - \tau_{\text{н}} .
Определение НСП оценки относительной разности частот необходимо осуществлять на основе значений составляющих НСП, используемых для расчетов смещений ШВ. При этом следует исключить все составляющие, имеющие постоянное значение, а именно задержки сигналов ГНСС в приемном тракте и задержки импульсных сигналов 1 Гц в кабелях.
1 При оценке относительной разности частот с применением приемников сигналов ГНСС, оснащенных выходом высокочастотных синусоидальных сигналов, для повышения точности измерений рекомендуется проводить повторные измерения и из полученного ряда вычислять среднее арифметическое значение.
2 Рекомендации по выбору интервалов времени измерения и интервалов времени наблюдений в зависимости от требуемой точности измерений приведены в [16].
5.3 Методы оценки среднего относительного изменения (дрейфа) частоты
Исходными данными для оценки среднего относительного изменения (дрейфа) частоты являются значения относительной разности частот, исследуемого и опорного СЧВ. Среднее относительное изменение (дрейф) частоты \nu вычисляют по формуле
\nu = \frac{\Delta_{0}f_{\text{к}} - \Delta_{0}f_{\text{н}}}{\tau_{\text{н}}} , (13)
где \Delta_{0}f_{\text{н}} — относительная разность частот исследуемого и опорного СЧВ, определенная в начале интервала времени наблюдений;
\Delta_{0}f_{\text{к}} — относительная разность частот исследуемого и опорного СЧВ, определенная в конце интервала времени наблюдений;
\tau_{\text{н}} — интервал времени наблюдений, где \tau_{\text{н}} \geq 10\tau_{\text{и}} .
1 Как правило, среднее относительное изменение (дрейф) частоты нормируют при конкретно определенных интервалах времени измерения; например, среднее суточное относительное изменение частоты, в данном случае для оценки среднего относительного изменения (дрейфа) частоты рекомендуемый интервал времени наблюдений составляет 10 суток.
2 Для повышения точности измерений рекомендуется использование всех значений из полученного по формулам (10) и (12) ряда. В этом случае оценку среднего относительного изменения (дрейфа) частоты осуществляют в соответствии с приложением 5 [16].
5.4 Методы оценки относительной нестабильности частоты
Нестабильность меры частоты — это характеристика случайных и систематических изменений частоты во времени. Для количественного описания нестабильности меры частоты применяют ряд характеристик.
В ГОСТ 8.567 в качестве количественных характеристик нестабильности частоты приведены СКО и СКДО частоты.
Исходными данными для оценки нестабильности частоты являются значения относительной разности частот, исследуемого и опорного СЧВ.
СКО частоты \sigma вычисляют по формуле
\sigma = \sqrt{\frac{\sum_{j = 1}^{n} (\Delta_{0}f_{j} - \Delta_{0}\overline{f})^{2}}{n - 1}} , (14)
где \Delta_{0}f_{j} — относительная разность частот исследуемого и опорного СЧВ, определенная в момент времени tj при интервале времени измерения \tau_{\text{и}} ;
\Delta_{0}\overline{f} — средняя относительная разность частот исследуемого и опорного СЧВ на интервале времени наблюдений при \tau_{\text{н}} \geq 10\tau_{\text{и}} ;
n — число результатов измерения при \tau_{\text{н}} = n \cdot \tau_{\text{и}} .
СКДО частоты \sigma_{\text{у}}(\tau_{\text{и}}) вычисляют по формуле
\sigma_{\text{у}}(\tau_{\text{и}}) = \sqrt{\frac{\sum_{j = 1}^{n} (\Delta_{0}f_{j + 1} - \Delta_{0}f_{j})^{2}}{2 \cdot (n - 1)} ,} (15)
где \Delta_{0}f_{j} — относительная разность частот исследуемого и опорного СЧВ, определенная в момент времени tj при интервале времени измерения \tau_{\text{и}} ;
n — число результатов измерения при \tau_{\text{н}} = n \cdot \tau_{\text{и}} .
1 Результаты оценки нестабильности частоты зависят от интервалов времени измерения, выборки и наблюдения. Как правило, интервалы времени измерения и выборки принимают равными, а интервал времени наблюдения принимают таким, чтобы количество результатов измерений на интервале времени наблюдений составляло не менее 10.
2 Не рекомендуется применять термины «кратковременная» или «долговременная» нестабильность частоты без указания значений интервалов времени измерения (выборки) и времени наблюдений.
СХЕМА ИЗМЕРЕНИЙ ПРИ ОЦЕНКЕ СМЕЩЕНИЙ (РАЗНОСТИ) ШКАЛ ВРЕМЕНИ СТАНДАРТОВ ЧАСТОТЫ И (ИЛИ) ВРЕМЕНИ С ПРИМЕНЕНИЕМ ПРИЕМНИКОВ СИГНАЛОВ ГЛОБАЛЬНЫХ НАВИГАЦИОННЫХ СПУТНИКОВЫХ СИСТЕМ, ОСНАЩЕННЫХ ВЫХОДОМ 1 Гц В АБСОЛЮТНОМ РЕЖИМЕ
СХЕМА ИЗМЕРЕНИЙ ПРИ ОЦЕНКЕ СМЕЩЕНИЙ (РАЗНОСТИ) ШКАЛ ВРЕМЕНИ СТАНДАРТОВ ЧАСТОТЫ И (ИЛИ) ВРЕМЕНИ С ИСПОЛЬЗОВАНИЕМ ПСЕВДОДАЛЬНОМЕРНЫХ КОДОВЫХ ИЗМЕРЕНИЙ В АБСОЛЮТНОМ РЕЖИМЕ
СХЕМА ИЗМЕРЕНИЙ ПРИ ОЦЕНКЕ СМЕЩЕНИЙ (РАЗНОСТИ) ШКАЛ ВРЕМЕНИ СТАНДАРТОВ ЧАСТОТЫ И (ИЛИ) ВРЕМЕНИ С ИСПОЛЬЗОВАНИЕМ ПСЕВДОДАЛЬНОМЕРНЫХ КОДОВЫХ ИЗМЕРЕНИЙ В ДИФФЕРЕНЦИАЛЬНОМ РЕЖИМЕ
СХЕМА ИЗМЕРЕНИЙ ПРИ ОЦЕНКЕ ОТНОСИТЕЛЬНОЙ РАЗНОСТИ ЧАСТОТ С ПРИМЕНЕНИЕМ ПРИЕМНИКОВ СИГНАЛОВ ГЛОБАЛЬНЫХ НАВИГАЦИОННЫХ СПУТНИКОВЫХ СИСТЕМ, ОСНАЩЕННЫХ ВЫХОДОМ ВЫСОКОЧАСТОТНЫХ СИГНАЛОВ
БИБЛИОГРАФИЯ
| [1] | Федеральный закон от 3 июня 2011 г. № 107-ФЗ «Об исчислении времени» | |
| [2] | Постановление Правительства Российской Федерации от 23 марта 2001 г. № 225 «Об утверждении Положения о Государственной службе времени, частоты и определения параметров вращения Земли» | |
| [3] | Ф. Риле. Стандарты частоты. Принципы и приложения. Перевод с английского Н.Н. Колачевского, Москва, Физматлит, 2009 | |
| [4] | К. Одуан, Б. Гино. Измерение времени. Основы GPS, Москва, Техносфера, 2002 | |
| [5] | RINEX The Receiver Independent Exchange Format Version 3.03 International GNSS Service (IGS), RINEX Working Group and Radio Technical Commission for Maritime Services Special Committee 104 (RTCM-SC104), July 14th, 2015 Acknowledgment: RINEX Version 3.02 and 3.03 is based on RINEX Version 3.01 which was developed by: Werner Gurtner, Astronomical Institute of the University of Bern, Switzerland and Lou Estey, UNAVCO, Boulder Colorado, USA | |
| [6] | D.W. Allan, M.A. Weiss. Requires Adobe Acrobat ReaderAccurate Time and Frequency Transfer During Common-View of a GPS Satellite, 34th Annual Frequency Control Symposium. Pp. 334 — 346. May. 1980 | |
| [7] | V.S. Zhang, T.E. Parker, M.A. Weiss, F.M. Vannicola. Requires Adobe Acrobat Reader Multi-Channel GPS/GLONASS Common-View between NIST and USNO, IEEE International Frequency Control Symposium. Pp. 598 — 606. June 2000 | |
| [8] | M.A. Lombardi, L.M. Nelson, A.N. Novick, V.S. Zhang. Requires Adobe Acrobat ReaderTime and Frequency Measurements Using the Global Positioning System, Cal. Lab. Int. J. Metrology. Pp. 26 — 33 (July — September 2001) | |
| [9] | J. Azoubib, W. Lewandowski. CGGTTS GPS/GLONASS Data Format Version 02 | |
| [10] | P. Defraigne, G. Petit. CGGTTS-Version 2E: an extended standard for GNSS Time Transfer. Royal Observatory of Belgium, BIPM | |
| [11] | Глобальная навигационная спутниковая система ГЛОНАСС. Интерфейсный контрольный документ. Навигационный радиосигнал в диапазонах L1, L2. (Редакция 5.1), Москва, 2008 | |
| [12] | J.F. Zumberge, M.B. Heflin, D.C. Jefferson, M.M. Watkins, F.H. Webb. Precise point positioning for the efficient and robust analysis of GPS data from large networks, J. Geophys. Res., 1997, vol. 102, no. B3. Pp. 5005 — 5017 | |
| [13] | J. Kouba. A Guide to Using International GNSS Service (IGS) Products. 2009 | |
| [14] | P. Collins, S. Bisnath, F. Lahaye, and P. Héroux (2010). «Undifferenced GPS ambiguity resolution using the decoupled clock model and ambiguity datum fixing», NAVIGATION: Journal of The Institute of Navigation, Vol. 57, No. 2, Summer 2010. Pp. 123 — 135 | |
| [15] | S. Loyer, F. Perosanz, F. Mercier et al. Zero-difference GPS ambiguity resolution at CNES — CLS IGS Analysis Center. J Geod 86. Pp. 991 — 1003 (2012) | |
| [16] | Рекомендации по метрологии МИ 2188-92 | Государственная система обеспечения единства измерений. Меры частоты и времени. Методика поверки |