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

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

Иллюстрации

Показать все

Реферат

Область техники, к которой относится изобретение

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

Уровень техники

Спутниковые системы позиционирования, например, глобальная система позиционирования (GPS) использует спутники в космическом пространстве для определения положения объектов на земле. В GPS сигналы от спутников поступают на приемник GPS и используются для определения положения приемника GPS. В настоящее время для гражданских приемников GPS доступны два типа измерений GPS, соответствующие каждому каналу коррелятора, с захваченным сигналом спутника GPS. Два типа измерений GPS - это псевдодальность и интегрированная фаза несущей для двух несущих сигналов, L1 и L2, с частотами 1,5754 ГГц и 1,2276 ГГц или длинами волны 0,1903 м и 0,2442 м соответственно. Измерение псевдодальности (или измерение кода) является основной наблюдаемой GPS, которую могут делать приемники GPS всех типов. Оно использует коды C/A или P, модулирующие несущие сигналы. Измерение записывает истинное время, необходимое соответствующему коду для распространения от спутника к приемнику, т.е. время поступления сигнала на приемник по часам приемника минус время выхода сигнала со спутника по часам спутника.

Измерение фазы несущей получается путем интегрирования реконструированной несущей сигнала при поступлении его на приемник. Таким образом, измерение фазы несущей также является мерой разницы по времени распространения, определяемой временем выхода сигнала со спутника по часам спутника и временем его поступления на приемник по часам приемника. Однако поскольку начальное количество полных периодов при распространении между спутником и приемником, когда приемник начинает отслеживать фазу несущей сигнала, обычно не известно, разница по времени распространения может иметь ошибку, составляющую множественные периоды несущей, т.е. в измерении фазы несущей может присутствовать неопределенность, составляющая целое число периодов. Для разрешения этой неопределенности, составляющей целое число периодов, были разработаны различные подходы. Хороший обзор таких подходов изложен в Главе 7, разделе 7.8, "Ambiguity Fixing," в GPS Спутник Surveying, 3-е издание, Альфреда Лейка (Alfred Leick), издательство John Wiley & Sons, Нью-Йорк, 2004.

Таким образом, для навигации объекта приемник GPS, присоединенный к объекту, принимает измерения GPS от совокупности спутников, и эти измерения GPS обрабатываются компьютерной системой, подключенной к приемнику GPS, для выработки обновлений положения объекта. Для вычисления обновлений положения объекта дальность или расстояние между приемником GPS и каждым из совокупности спутников вычисляется путем умножения времени распространения сигнала на скорость света. Эти дальности обычно называются псевдодальностями (ложными дальностями), поскольку часы приемника обычно имеют значительную погрешность по времени, которая приводит к общему смещению в измеренной дальности. Это общее смещение из-за погрешности часов приемника можно разрешить совместно с пространственными координатами приемника GPS в порядке обычных навигационных расчетов. К ошибкам или шуму в расчетной дальности также могут приводить различные другие факторы, включая эфемеридную ошибку, ошибку хронирования часов спутника, атмосферные эффекты, шум приемника и ошибку многолучевого распространения. При автономной GPS-навигации, когда приемник GPS получает дальности на основании кода и/или фазы несущей в отношении совокупности спутников в поле зрения, без консультации с какой-либо опорной станцией, обычно бывает очень трудно снизить погрешности или шумы в отношении дальностей.

Для устранения или снижения этих погрешностей в приложениях GPS обычно используются дифференциальные операции. Дифференциальные операции GPS (DGPS) предусматривают базовую или опорную станцию, имеющую базовый или опорный приемник GPS, систему навигации, имеющую полевой или навигационный приемник GPS, присоединенный к объекту, в отношении которого осуществляется навигация, и канал связи между опорной станцией и системой навигации. Опорная станция обычно находится в известном положении, и измерения, полученные и/или корректировки, вычисленные на ней, передаются на и используются системой навигации для исключения большинства факторов ошибки. Дифференциальные операции с использованием измерений фазы несущей часто называют операциями позиционирования/навигации на основе кинематики в реальном времени (RTK).

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

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

Когда ошибки более или менее скорректированы с использованием автономной системы, системы RTK и/или WADGPS, компьютерная система, подключенная к навигационному приемнику GPS, обычно вычисляет обновление положения приемника один раз в каждом из ряда периодов. Было разработано несколько разных методов обновления положения приемника GPS с использованием измерения фазы несущей GPS. Примерами этих методов являются процессы наименьших квадратов и фильтра Кальмана. Метод наименьших квадратов подробно рассмотрен в Главе 4, "Least Squares Adjustments," в GPS Спутник Surveying, 3-е издание, Альфреда Лейка (Alfred Leick), издательство John Wiley & Sons, Нью-Йорк, 2004. Фильтр Кальмана кратко описан в конце главы. В традиционной системе навигации GPS положение приемника GPS обычно обновляется с частотой 1 герц, т.е. одно обновление в секунду. Для объектов со средней и высокой динамикой часто желательно, чтобы обновления положения вычислялись с более высокой частотой. Однако, ввиду сложности полного процесса наименьших квадратов или фильтра Кальмана, вычислительных возможностей компьютерной системы может не хватить, если обновления положения будут вычисляться с частотой, значительно превышающей 1 герц. Поэтому необходим метод спутниковой навигации, позволяющий значительно увеличить частоту, с которой производятся обновления точного положения, и одновременно минимизировать сложность вычислений.

Сущность изобретения

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

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

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

Согласно другому варианту осуществления настоящего изобретения способ навигации объекта осуществляется в так называемом режиме работы «максимальной точности», в котором корректировки положения, получаемые на низкой частоте, не используют для коррекции обновлений положения, получаемых на высокой частоте. Вместо этого высокоскоростной процесс выполняют в кадре относительного положения, привязанном к опорной точке. Для минимизации роста ошибки со временем в высокоскоростном процессе, изменения в измерениях фазы несущей, используемых для вычисления обновлений положения, регулируют с использованием всей имеющейся информации, включая, но без ограничения: (1) измерения фазы несущей на двух частотах для устранения эффектов рефракции в ионосфере; (2) тщательно смоделированные тропосферные эффекты; (3) корректировки часов и орбиты спутника, обеспечиваемые широкомасштабной или глобальной спутниковой системой позиционирования; и т.д. Этот режим работы очень полезен в случаях, когда нужны обновления относительного положения вместо обновлений абсолютного положения. Преимущество использования этого режима работы проявляется в системе RTK. Он позволяет системе RTK, которая обеспечивает навигацию относительно положения базового приемника RTK, продолжать работать в режиме относительной навигации, даже когда канал связи между мобильным приемником и базовой станцией, который обеспечивает корректировки измерений фазы несущей, прерывается.

Краткое описание чертежей

Фиг.1 - блок-схема спутниковой системы навигации, имеющей компьютерную систему, подключенную к приемнику GPS.

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

Фиг.3 - логическая блок-схема, более подробно иллюстрирующая аспект способа навигации согласно одному варианту осуществления настоящего изобретения.

Фиг.4 - логическая блок-схема, более подробно иллюстрирующая другой аспект способа навигации согласно одному варианту осуществления настоящего изобретения.

Фиг.5 - логическая блок-схема, более подробно иллюстрирующая аспект способа навигации согласно альтернативному варианту осуществления настоящего изобретения.

Фиг.6 - логическая блок-схема, более подробно иллюстрирующая другой аспект способа навигации согласно одному варианту осуществления настоящего изобретения.

Фиг.7 - логическая блок-схема, иллюстрирующая способ редактирования измерений GPS согласно одному варианту осуществления настоящего изобретения.

Описание

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

Для облегчения дифференциальных операций система 100 также может быть подключена к опорной станции 120 посредством канала связи 124. Опорная станция 120 обеспечивает наблюдаемые GPS, измеренные на ней и/или корректировки ошибки GPS, вычисленные на ней. В широкомасштабных или глобальных приложениях система 100 может быть подключена к одному или нескольким центральным концентраторам 130, поддерживающим связь с группой опорных станций (не показана) по радио- и/или спутниковым каналам связи 134. Концентратор(ы) 130 принимает(ют) наблюдаемые GPS от группы опорных станций и вычисляет(ют) корректировки, которые передаются в систему 100.

Согласно одному варианту осуществления настоящего изобретения система 100 включает в себя центральный процессор (CPU) 140, память 148, совокупность входных портов 153, 154 и 155, один или несколько выходных портов 156 и необязательный пользовательский интерфейс 158, связанные между собой одной или несколькими шинами связи 152. Входные порты 154 служат для приема данных от приемника 110 GPS, опорной станции 120 и/или концентратора 130 соответственно и выходной(ые) порт(ы) 156 можно использовать для вывода результатов вычислений. Альтернативно, результаты вычислений можно отображать на устройстве дисплея пользовательского интерфейса 158.

Память 148 может включать в себя высокоскоростную оперативную память и может включать в себя энергонезависимую память большой емкости, например, одно или несколько запоминающих устройств на основе магнитных дисков. Память 148 также может включать в себя память большой емкости, размещенную дистанционно от центрального процессора 140. В памяти 148 предпочтительно хранится операционная система 162, база данных 170 и программы или процедуры 164 приложения GPS, в том числе процедуры 166, реализующие способ навигации с использованием последовательных изменений в измерениях фазы несущей согласно одному варианту осуществления настоящего изобретения. Операционная система 162 и прикладные программы и процедуры 164, хранящиеся в памяти 148, должны выполняться CPU 140 компьютерной системы 100. В памяти 148 предпочтительно также хранятся структуры данных, используемые при выполнении процедур 166 приложения GPS, например, измерения и корректировки GPS, а также другие структуры данных, рассматриваемые в этом документе.

Операционная система 162 может являться, но без ограничения, встроенной операционной системой, UNIX, Solaris или Windows 95, 98, NT 4.0, 2000 или XP. В более общем случае операционная система 162 имеет процедуры и инструкции для передачи, обработки, доступа, сохранения и поиска данных.

По причинам, указанным ниже, в памяти 148 также может храниться исполнимый модуль реального времени (RTX) 167, который является компьютерной программой для многозадачных операций в реальном времени. Согласно одному варианту осуществления настоящего изобретения RTX 167 позволяет встраивать операционную систему 162 в процедуры 166 для обеспечения множественных потоков, что позволяет разным задачам в процедурах 166 выполняться "квази-одновременно", то есть может казаться, что разные задачи выполняются в одно и то же время и что система 100 одновременно производит разные операции. Это позволяет процедурам 166 включать в себя два или более параллельных задания или процесса, которые выполняются в разных потоках. RTX 167 управляет началом и окончанием каждого потока и позволяет потокам взаимодействовать друг с другом. RTX 167 также позволяет ставить данные в очередь на поток, передавать данные между потоками и приводить процесс к последовательному виду, упорядочивая события.

Кроме того, RTX 167 поддерживает стандартное управление множественными потоками, что позволяет выполнению потока ожидать событий, которые инициируются событиями в другом потоке. Событие - это состояние, которое может быть задано или очищено на потоке. Когда потоку предписано ожидать одно или несколько событий, поток приостанавливается, пока все события не будут заданы. Это значительно упрощает синхронизацию и связь между потоками. Благодаря RTX 167 выполнение потока осуществляется на основании приоритета. Поток с более высоким приоритетом выполняется прежде потока с более низким приоритетом. Среди потоков с одинаковым приоритетом потоки выполняются по циклическому принципу. Каждому потоку выделяется квант времени для выполнения. В качестве RTX 167 можно использовать готовый коммерческий RTX. Примеры таких коммерчески доступных исполнимых модулей реального времени включают в себя CMX-RTX от CMX Systems, Inc., Concurrent Real time Executive (CORTEX) от Australian Real Time Embedded Systems (ARTESYS) и Nucleus RTX, от Accelerated Technology Inc.

Как указано пунктирной линией 105 на фиг.1, в некоторых вариантах осуществления приемник 110 GPS и частично или полностью компьютерная система 100 интегрированы в единое устройство, в едином корпусе, например, переносное, ручное или даже карманное устройство отслеживания положения, или автомобильную или иную мобильную систему позиционирования и/или навигации. В других вариантах осуществления приемник 110 GPS и компьютерная система 100 не интегрированы в единое устройство.

На фиг.2 показана логическая блок-схема, иллюстрирующая способ 200 навигации объекта с использованием последовательных разностей в измерениях фазы несущей согласно одному варианту осуществления настоящего изобретения. Способ 200 реализован в процедурах 167 и может выполняться системой 100. Согласно фиг.2 способ 200 включает в себя процесс 201 инициализации и два одновременных процесса, высокоскоростной процесс 210 и низкоскоростной процесс 220. Процесс 201 инициализации используется для вычисления начального положения объекта и других начальных параметров, необходимых высокоскоростному процессу 210. Высокоскоростной процесс 210 включает в себя ряд процессов 215 распространения положения, в каждом из которых вычисляется обновление положения объекта в одном из ряда малых периодов. Низкоскоростной процесс 220 включает в себя ряд процессов 225 коррекции положения, в каждом из которых вычисляется коррекция положения объекта в одном из ряда больших периодов. Согласно фиг.2B в каждой большом периоде T_m может быть несколько, например 10, малых периодов t_mn (m = 0, 1, 2, 3,... и n=0, 1, 2, 3,...). В зависимости от режима работы, малые периоды могут также совпадать с большими периодами, что рассмотрено ниже.

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

где x - вектор коррекции состояния, выражающий коррекцию состояния процесса, управляемого в дискретном времени, которое в этом случае может включать в себя положение приемника и время по часам приемника; z представляет значение обновления измерения, которое задается разностью между измерением относительно спутника и ожидаемым значением для измерения, которое вычисляется из первоначально оцененного состояния; n представляет шум в измерениях; h - вектор чувствительности измерения, который характеризует чувствительность измерения к изменению состояния. Вектор h сформирован путем разложения уравнения, связывающего псевдодальность измерения с положением приемника GPS в ряд Тейлора. Элементы h включают в себя первые производные обновления измерения относительно вектора коррекции. Вектор коррекции состояния x включает в себя, по меньшей мере, одну коррекцию положения приемника GPS. Он также может включать в себя коррекцию часов приемника GPS. В случае измерений фазы несущей вектор коррекции состояния может включать в себя неизвестный коэффициент неопределенности, который можно разрешить с использованием традиционного метода разрешения неопределенности фазы несущей. Для облегчения нижеследующего рассмотрения предположим, что вектор состояния является четырехэлементным вектором, т.е. состояние включает в себя только положение приемника и время по часам приемника.

Уравнение (1) можно разложить, чтобы связать вектор коррекции состояния с множеством измерений в отношении множественных спутников в общем периоде:

В этом уравнении z - вектор, состоящий из обновлений в отношении множественных спутников, H - матрица, состоящая из чувствительностей измерения в отношении множественных спутников, x - по-прежнему вектор состояния и n - вектор шума измерения, включающий в себя множество значений шума измерения, связанных с обновлениями в z. Матрица чувствительности измерения H зависит от геометрии спутников, которая относится к совокупности геометрических соотношений между приемником и спутниками. Обновления измерения в векторе z часто именуются префиксными остатками. Решая уравнение (2) в наименьших квадратах, получим:

где верхний индекс T представляет операцию транспонирования и верхний индекс "-1" представляет операцию обращения матрицы.

Альтернативой решению уравнения (2) с использованием уравнения (3) является решение методом взвешенных наименьших квадратов, которое имеет вид:

где W - ковариационная матрица измерений, диагональные элементы которой выражают среднеквадратичные отклонения шума измерения в векторе шума n и недиагональные элементы представляют смешанный момент второго порядка между измерениями. Поскольку смешанный момент второго порядка между измерениями обычно предполагается равным нулю, недиагональные элементы W обычно равны нулю.

Для простоты описания в следующих рассмотрениях используется простое уравнение в средних квадратах, уравнение (3). Уравнение (3) можно дополнительно упростить:

или

где A = (H^TH)^-1 и B = AH^T.

Иногда бывает полезно строить остаточную матрицу чувствительности, S, которая отображает обновления, z, или префиксные остатки, в постфиксные остатки, соответствующие измерениям GPS в обновлениях z. Постфиксные остатки выражаются как элементы в остаточном векторе Δ:

где

где I - квадратная единичная матрица, ранг которой равен количеству измерений или количеству элементов в z.

В процессе инициализации 201 коррекция положения в векторе коррекции, x, прибавляется к первоначально оцененному положению для обеспечения скорректированной оценки положения приемника. Коррекция часов приемника часто рассматривается как мешающий параметр и не обновляется. Это возможно, поскольку зависимость от часов приемника является линейной и большие ошибки в этом значении не влияют на решение. Поскольку уравнения дальности, уравнения (1)-(7), являются нелинейными, может потребоваться повторное выполнение всех вычислений, если оценка начального положения имеет большую ошибку. Таким образом, процесс 201 инициализации можно повторять с использованием того же множества измерений, которое получено в первом периоде, или, альтернативно, его можно повторять по нескольким последовательным периодам измерения. Таким образом, несколько больших периодов в процессе 201 инициализации можно взять до получения точного начального положения и ассоциированных матриц A и H (или B и S).

После вычисления начального положения и ассоциированных матриц A и H (или B и S) начинаются высокоскоростной процесс 210 и низкоскоростной процесс 220. В низкоскоростном процессе 220 вычисляется коррекция положения приемника и матрицы A и H (или B и S) в каждом большом периоде с использованием процесса 225 коррекции положения, который включает в себя вычисление наименьших квадратов, подобно описанному выше согласно уравнениям (1)-(7). Высокоскоростной процесс 210 распространяет положение приемника 110 GPS вперед по времени путем вычисления обновленного положения приемника 110 GPS в каждом из ряда малых периодов после процесса инициализации с использованием процесса 215 распространения положения. Согласно одному варианту осуществления настоящего изобретения осуществлением процессов 210 и 220 управляет RTX 167, который начинает два отдельных потока после процесса инициализации 201 для выполнения процессов 210 и 220 соответственно. Поскольку процессы 225 осуществляются с более низкой частотой, чем процессы 215 в высокоскоростном процессе 210, и поскольку осуществление каждого процесса 225 может зависеть от получения тех или иных результатов вычислений из высокоскоростного процесса, что более подробно объясняется ниже, более высокий приоритет может назначаться потоку, в котором выполняется высокоскоростной процесс 210, и более низкий приоритет может назначаться потоку, в котором выполняется низкоскоростной процесс 220. RTX также управляет постановкой в очередь данных на каждом потоке и передачей данных между потоками, что более подробно объясняется ниже.

Поскольку измерения фазы несущей обычно выполняются с точностью менее одного сантиметра, их можно использовать для распространения положения приемника вперед по времени с очень малым ростом ошибки. Согласно фиг.3 процесс распространения положения 215 в текущем малом периоде в высокоскоростном процессе может включать в себя этап 310, на котором наиболее недавно вычисленная коррекция положения приемника и матрицы A и H (или B и S) получается из очереди, связанной с низкоскоростным процессом 220. Матрицы A и H (или B и S) могут нуждаться в обновлении для учета вкладов от одного или нескольких спутников, не включенных в вычисление матриц посредством низкоскоростного процесса 220, или для удаления вкладов от одного или нескольких спутников, которые привели к ненадежным измерениям. Поскольку обновление матриц не является обязательным, эта часть этапа 310 рассматривается ближе к концу описания.

Процесс 215 дополнительно включает в себя этап 320, на котором скорость приемника 110 GPS используется для вычисления приблизительного распространения позиции. Согласно одному варианту осуществления настоящего изобретения приблизительное обновление положения приемника в каждом малом периоде в высокоскоростном процессе 210 сначала вычисляется путем прибавления к положению приемника, вычисленному в предыдущем малом периоде, начальной коррекции, которая является произведением скорости и времени, истекшего после предыдущего периода. Получаем:

где представляет скорость и Δt - время, прошедшее между двумя малыми периодами или интервал малого периода, и представляет расчетное положение приемника 110 GPS. Этап 320 включен в процесс 215, чтобы гарантировать, что наилучшая оценка положения используется при высокоскоростном вычислении коррекции благодаря поддержанию малости обновлений (префиксных остатков). На самом деле, этот этап не требуется, поскольку распространение положения с использованием изменения фазы несущей, согласно описанному ниже, должно автоматически учитывать отсутствие начальной коррекции, вычисленной на этапе 320. Кроме того, наибольший вклад в префиксные остатки обычно дает ошибка часов приемника, которую можно включать, а можно и не включать в вычисления на этапе 320.

Опять же согласно фиг.3 процесс 215 дополнительно включает в себя этап 330, на котором изменения в измерениях фазы несущей в течение интервала времени между текущим и предыдущим малыми периодами используются для вычисления изменения положения приемника в течение того же интервала времени. Обновление скорости на этапе 320 является хорошей оценкой положения приемника в текущем малом периоде. Но оно не учитывает никакого немоделированного шума ускорения и скорости, который можно скорректировать с использованием изменений в измерениях фазы несущей. Для вычисления изменения положения с использованием изменений в измерениях фазы несущей используется уравнение (5a) с матрицами A и H из процесса инициализации (если текущий малый период является одним из нескольких первых малых периодов после процесса инициализации) или из высокоскоростного процесса. Более эффективно использовать матрицу B и уравнение (5b). Если же используется матрица B, то матрица S и уравнение (6) также можно использовать для обеспечения альтернативного способа компенсации ухода на период или потери сигнала в ходе высокоскоростной обработки, что рассмотрено ниже в связи со способом компенсации потери сигнала спутника согласно одному варианту осуществления настоящего изобретения. Матрица S обычно весьма нечувствительна к положению пользователя. Поэтому необходимо производить повторные вычисления, только если было пройдено значительное расстояние после последнего процесса 225 коррекции положения.

Используются ли матрицы A и H или B и S, поскольку они вычисляются в процессах 225 коррекции положения один раз за каждый большой период, их можно повторно использовать для вычисления высокочастотных обновлений положения в последовательности малых периодов с достаточной точностью и, таким образом, не нужно повторно вычислять в каждом малом периоде. Это существенно снижает вычислительную нагрузку в высокоскоростном процессе 210, в котором единственными конкретными значениями, которые нужно повторно вычислять в каждом малом периоде для реализации уравнения (5a) или уравнений (5b) и (6), являются элементы вектора обновлений, z. Для каждого спутника, участвующего в вычислении матриц:

где ϕ используется для представления измерения фазы несущей, приведенной к соответствующей длине волны в метрах, и ρ используется для представления теоретической дальности между спутником и приемником. Нижний индекс, i, используется для представления текущего периода, и нижний индекс, i-1, используется для представления предыдущего периода. Теоретическую дальность ρ_i можно вычислять с использованием положения приемника, вычисленного в малом периоде i-1, которое можно обновлять или не обновлять с использованием уравнения (8). Теоретическую дальность ρ_i-1 можно просто вычислять с использованием положения приемника, вычисленного в малом периоде i-1. Оно не обновляется с использованием уравнения (8).

Процесс 200 может выполняться в разных режимах работы. Согласно одному варианту осуществления настоящего изобретения процесс 200 выполняется в режиме, именуемом режимом "максимальной доступности" или "высокочастотного распространения". В этом режиме работы не нужно использовать никаких поправок на эффекты рефракции в ионосфере или тропосфере, и измерения в отношении спутников с возвышением в несколько градусов могут быть включены при вычислении обновлений в уравнении (9). Эффекты рефракции можно игнорировать, поскольку изменения в эффектах рефракции в течение большого интервала (т.е. промежутка времени, соответствующего каждому большому периоду), который обычно составляет около 1 секунды, обычно ниже шума измерения фазы несущей. Таким образом, измерения фазы несущей L1 можно использовать без применения ионосферной поправки из измерений фазы несущей L2. Поскольку сигнал L1 сильнее сигнала L2 и измерения L1 легче получить в пограничных условиях, результат, полученный с использованием только измерений L1, может быть более надежным. Альтернативно, когда доступны измерения L2, можно также использовать среднее для измерения фазы несущей L1 и соответствующего L2 (переведенных в метры до усреднения). Это дает немного меньший шум в расчетном положении, поскольку шумы в измерениях L1 и L2 до некоторой степени устраняются в ходе процесса усреднения.

Таким образом, обновления, вычисленные с использованием уравнения (9), используются в уравнении (5a) или (5b) для вычисления изменения x состояния, которое включает в себя изменение положения приемника . Коррекцию скорости также можно генерировать из изменения положения , просто разделив изменение положения на интервал малого периода. Таким образом, коррекция скорости задается как:

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

Процесс 215 дополнительно включает в себя этап 340, на котором изменение положения, генерируемое из изменений в измерениях фазы несущей, используется для обновления положения приемника с использованием:

где представляет коррекцию положения приемника . В некоторых вариантах осуществления положение приемника определяется в декартовых координатах xyz и относительно геоцентрической, привязанной к земле, системы координат. Если определен в системе координат, оси которой направлены на север, восток, вверх, изменение положения нужно сначала умножить на соответствующую матрицу поворота, R, которая вычисляется в низкоскоростном процессе и передается высокоскоростному процессу совместно с матрицами A и H (или B и S). В этом случае

где представляет изменение положения в системе координат, оси которой направлены на север, восток, вверх, и представляет положение приемника в декартовых координатах.

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