СПОСОБ ОПРЕДЕЛЕНИЯ ВЕКТОРА СОСТОЯНИЯ ПОДВИЖНОГО ОБЪЕКТА ПРИ ПОМОЩИ НАВИГАЦИОННЫХ КОСМИЧЕСКИХ АППАРАТОВ КОСМИЧЕСКОЙ НАВИГАЦИОННОЙ СИСТЕМЫ Российский патент 1996 года по МПК G01C21/00 

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

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

При этом возможно использование сформированной поправки как для коррекции координат, так и для коррекции навигационного параметра.

Недостатком известного способа является наличие погрешности в определении координат на уровне 4 6 метров.

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

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

Наряду с перечисленными выше недостатками снижение точности при использовании НКА КНС Navstar вызвано наличием режима ограниченного доступа - S/A Selective Availability введенного министерством обороны США искусственного загрубления навигационных сигналов открытого доступа.

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

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

Способ осуществляется следующим образом.

При использовании навигационной аппаратуры потребителя фирмы Trimble Navig Ltd возможно получение в специальном режиме работы результатов измерений псевдодальностей и доплеровского сдвига частоты несущей сигнала на частоте 1575,42 МГц по всем НКА, сигналы которых сопровождались приемным устройством на текущий момент времени.

В стандартном режиме работы шестиканальная НАП указанных выше моделей одновременно имеет возможность сопровождать сигналы 6 НКА и может периодически или в режиме мультиплексирования следить за 8 НКА, находящимися в зоне радиовидимости ее антенного устройства. Это позволяет существенно упростить реализацию условия идентичности НКА в рабочих созвездиях для эталонного пункта и подвижного объекта, чей вектор состояния определяется. Так, при записи на эталонном пункте навигационных параметров практически по всем видимым НКА можно гарантировать, что четыре или как минимум три из них войдут в состав рабочего созвездия второго комплекта НАП, установленного на подвижном объекте при расстоянии между пунктами по крайней мере до 500 километров.

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

Координаты эталонного пункта при этом известны с высокой точностью. Они определяются заранее, например, при помощи геодезической привязки.

Вычисление положения подвижного объекта при этом производится с использованием следующих выражений:


где оценка скорректированного вектора положения подвижного объекта на момент времени t_i;
вектор положения подвижного объекта, полученный на момент времени t_i в результате решения навигационной задачи в НАП, установленной на подвижном объекте;
вектор дифференциальной коррекции координат подвижного объекта;
H_п(t_i) градиентная матрица, рассчитываемая с учетом определенных координат подвижного объекта;
вектор дифференциальной коррекции разности псевдодальностей, измеренных на подвижном объекте.

Разности псевдодальностей и на эталонном пункте, и на подвижном объекте формируются следующим образом: один из спутников рабочего созвездия выбирается ведущим обычно это самый "высокий" изо всех видимых НКА и измерения по этому спутнику вычитаются из измерений по остальным спутникам созвездия.

Вектор дифференциальной коррекции определяется в соответствии со следующими выражениями:

ΔX_к(t_i) = X_к(t_i)-X_p
ΔY_к(t_i) = Y_к(t_i)-Y_p
ΔZ_к(t_i) = Z_к(t_i)-Z_p
где результаты измерений псевдодальности по сигналам k-го и ведущего (0) НКА соответственно, полученные на эталонном пункте и подвижном объекте,
результаты по дальности до k-го и ведущего (0) НКА соответственно, полученные на подвижном объекте при помощи вычислителя, входящего в состав НАП, установленной на подвижном объекте,
X_k, Y_k, Z_k координаты k-го НКА в геоцентрической системе координат,
X_p, Y_p, Z_p координаты подвижного объекта или эталонного пункта в геоцентрической системе координат.

Выражения для градиентной матрицы при этом носят следующий вид:
H_к.1 = ΔX_к(R_к-ΔX_o)R_o
H_к.2 = ΔY_к(R_к-ΔY_o)R_o
H_к.3 = ΔZ_к(R_к-ΔZ_o)R_o
Использование в качестве навигационных параметров, по которым решается навигационная задача по определению вектора состояния подвижного объекта, разностей псевдодальностей позволяет упростить вычисление обратной градиентной матрицы, а также скомпенсировать влияние различных величин уходов временной школы НАП на эталонном пункте и подвижном объекте.

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

R_x(t_i), R₀(t_i) расчетные дальности до k-го и ведущего НКА, полученные на эталонном пункте.

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

где оцениваемая величина на k-том шаге сглаживания;
ее оценки на k-том и предыдущем шагах;
К весовой коэффициент.

Предварительная фильтрация указанным алгоритмом корректирующих поправок разностно-дальномерных измерений вида (δDR(t_i)) не дает положительных результатов, возможно из-за временной декорреляции, вызванной влиянием Selec- tive Availability, особенно при интервалах между измерениями более 10 секунд.

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

Для вычисления вектора V_п движущейся НАП использовалось следующее выражение:

где оценка скорректированного вектора скорости НАП на привязываемом маршруте на момент времени;
вектор дифференциальной коррекции разностей радиальных квазискоростей, измеренных НАП на привязываемом маршруте.

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

K 1, 2, 3.

где результаты измерений радиальных квазискоростей по доплеровским сигналам k-го и ведущего (0) НКА, соответственно полученные на опорном и привязываемом пунктах;
V_к(t_i), V₀(t_i) результаты расчета радиальной скорости k-го и ведущего (0) НКА, соответственно полученные на опорном пункте;
результаты расчета скорости k-го и ведущего (0) НКА соответственно, полученные на привязываемом пункте;
производные координат k-го НКА в геоцентрической системе координат;
производные координат привязываемого или опорного пунктов в геоцентральной системе координат.

В этом случае вектор дифференциальной коррекции положения движущейся НАП определяется с учетом вычисленного вектора согласно следующему выражению:

при условии K_c + K 1.

где X_п(t_i) скорректированный по разностно-дальномерным измерениям дифференциальным методом вектора положения НАП, установленной на движущемся объекте.

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

Кроме того, так как алгоритм интегрирования измерений доплеровского сдвига несущей последователен, то в значениях псевдодальности, полученных с помощью интеграла от доплеровского сдвига частоты на начальном участке времени содержится остаточный шум из-за недостаточного числа измерений, по которым формируется среднее значение псевдодальности на момент времени t₀.

Поэтому до тех пор, пока позволяет одинаковость условий (например, до тех пор, пока не изменилось рабочее созвездие или не изменились эфемериды НКА) производится вычисление значений интегрированного доплера на текущий момент времени и осредняются значения псевдодальности, приведенные на начальный момент времени. При смене рабочего созвездия на подвижном пункте или на эталонном пункте вновь выбирается рабочее созвездие не менее, чем из трех НКА и начинается формирование псевдодальности на момент tj+1o

Алгоритм коррекции данных дальномерных измерений по фазе С/A кода КНС GPS Navstar с помощью интегрирования измерений доплеровского сдвига частоты несущей сигналов НКА имеет следующий вид:
ID_к(t_i) ID_к(t_i-1) + (PV_к(t_i) PV₀(t_i)) • (t_i-t_i-1), (12)
dPRO_к(t_i) (PR_к(t_i PR_o(t_i)) ID_к(t_i), (13)
dPRO_к(t_ik) (dPRO_к(t_i-1) • (i-1) + dPRO_к(t_i))/i, (14)
где ID_к(t_i), ID_к(t_i-1) интеграл от доплеровского сдвига частоты, соответствующий счетному значению разности псевдодальностей пары НКА на t_i и t_i-1 моменты времени для k-го НКА;
dPRO_к(t_i) разница псевдодальностей k-го и ведущего НКА на момент времени t₀, полученная в результате измерений в момент t_i;
dPRO_к(t_i), dPRO_к(t_i-1) осредненное значение dPRO_к(t_i).

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

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

По полученным разностям псевдодальностей на момент времени путем решения классической навигационной задачи (см. например, 1) определяется для каждого момента времени вектор состояния подвижного объекта.

На фиг. 1 представлены на плоскости х, y в системе координат Гаусса-Крюгера результаты одной из серий навигационных определений привязываемого пункта до и после проведения дифференциальной коррекции. Для наглядности точки фиксации положения на плоскости соединены отрезками последовательно по времени их получения.

Результаты фильтрации скорректированных координат с учетом значений информации о векторе V_п согласно формуле (II) проиллюстрированы на фиг. 2 в сравнении с результатами координат только по разностям псевдодальностей без фильтрации.

На фиг. 3 приведены результаты коррекции координат, определяемых по измерениям псевдодальности по фазе C/A кода НКС GPS Navstar с помощью интегрирования измерений доплеровского сдвига частоты несущей сигналов НКА.

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

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

Для проведения предварительного тестирования результатов применения дифференциальных алгоритмов при обработке информации снимаемой с НАП было применено программное обеспечение PFINDER фирмы Trimble Navigation (Pfinder User Guide-Trimble Navigation Ltd 1993), позволяющее обрабатывать аналогичные данные и обеспечивающее точность (2dMRS) 2-5 м. Результат применения этой программы приведен фиг. 7.

Качественный сравнительный анализ результатов, представленных на фиг. 3, 4 и фиг. 7, а также данные статистической обработки, полученные при моделировании, позволяют сделать вывод о достаточной степени эффективности используемых простейших алгоритмов дифференциальной обработки навигационных данных, получаемых следующими моделями НАП: TRANSPAK-11, ACUTIME, Sveesix.

Можно сделать следующие выводы:
информация, поступающая с НАП Trimble Navigation, является достаточной для реализации двух основных методов дифференциальной коррекции навигационных измерений;
использование метода поддержки дальномерных измерений по фазе С/A-кода с помощью интегрирования информации о доплеровском сдвиге частоты несущей позволяет снизить погрешность результатов дифференциальной обработки до 2 м даже без использования специальных процедур фильтрации.

Изобретение относится к области навигации, в частности, космической. Техническим эффектом использования изобретения является повышение быстродействия навигационных измерений при сохранении точности. Технический результат достигается тем, что при измерении вектора состояния определяемого объекта с использованием рабочего созвездия НКА навигационная задача решается на основе разности псевдодальностей и доплеровского сдвига частоты на временном интервале, значение разностей псевдодальностей на начало интервала формируется путем осреднения значений разностей псевдодальностей, полученных на всем интервале, привязка ко времени осуществляется путем решения навигационно-временной задачи, а текущие значения разностей псевдодальностей формируются с учетом осредненного значения на начало интервала и интеграла разности доплеровского сдвига частоты на этом интервале. 3 з.п. ф-лы, 7 ил.

1. Способ определения вектора состояния объекта при помощи навигационных космических аппаратов (НКА) космической навигационной системы (КНС), заключающийся в измерении на определенном объекте и эталонном пункте с привязкой ко времени псевдодальностей и доплеровского сдвига частот по навигационным сигналам НКА, образующим одну группу НКА-рабочее созвездие в составе КНС, и определении на эталонном пункте поправки к навигационному параметру НКА, передаче полученной поправки на определяемый объект и определении вектора состояния определяемого объекта с учетом измеренных навигационных параметров и поправок к ним, отличающийся тем, что для измерения выбирают созвездие НКА с числом не менее трех, на эталонном пункте и определяемом объекте выбирают один НКА-ведущий, формируют разности псевдодальностей до ведущего НКА и каждого НКА рабочего созвездия, формируют разности доплеровского сдвига частоты по ведущему НКА и каждому НКА рабочего созвездия, формируют разность псевдодальностей на начальный момент t₀ по текущим значениям разности псевдодальностей и текущим значениям разности доплеровского сдвига частоты как разности текущих значений разностей псевдодальностей и интеграла от разностей доплеровского сдвига частоты, при этом t₀ определяется по смене созвездий, а текущее значение вектора состояния определяемого объекта формируют с учетом осредненных значений разности псевдодальностей, сформированных как сумма осредненных значений разностей псевдодальностей на момент t₀ и интеграла от разности доплеровского сдвига частоты на интервале от начального момента времени t₀ до момента, на который определяется вектор состояния, при изменении рабочего созвездия НКА операции формирования текущего значения вектора состояния повторяют. 2. Способ по п.1, отличающийся тем, что при наличии в рабочем созвездии трех НКА запоминают последнее значение высоты, полученное на определяемом объекте путем решения навигационной задачи при наличии созвездия, состоящего из четырех и более НКА, и используют запомненное значение высоты при решении навигационной задачи совместно с значениями параметров, полученных на основе текущих навигационных измерений для определения вектора состояния. 3. Способ по п.1, отличающийся тем, что при наличии в рабочем созвездии трех НКА, используется значение высоты, получаемое от внешнего источника информации при решении навигационной задачи совместно с значениями параметров, полученных на основе текущих навигационных измерений для определения вектора состояния. 4. Способ по п. 2, отличающийся тем, что при появлении в рабочем созвездии четвертого НКА полученный на предыдущем измерении вектор состояния определяемого объекта уточняется с учетом решения навигационной задачи по информации от не менее чем четырех НКА.