Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов

В соответствии с классической постановкой задачи опре­деления орбиты ИСЗ с использованием метода наименьших квад­ратов определим начальные условия движения ИСЗ на момент прохождения им экваториальной плоскости Земли. Схема опреде­ления орбиты ИСЗ приведена на рис. 3.15. Пусть орбита определя­ется по измерениям даль­ности с наземных измери­тельных пунктов — НИПов, число которых варьируется, .например, от одного до трех. Тогда на каждом витке по­лета ИСЗ может быть про­ведено от одного до трех се­ансов измерений. С целью выявления особенностей про­цедуры метода наименьших квадратов (МНК) примем достаточно простые матема-

Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru

Рис. 3.15. Схема организации экс­перимента по определению орбиты ИСЗ

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

Уравнения движения ИСЗ в центральном поле тяготения Зем­ли, в экваториальной системе координат, имеют вид

Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru

Расчетные (опорные) начальные условия для этой системы уравнений определяются расчетной орбитой выведения и точно­стью выведения ИСЗ носителем. В стохастической интерпретации применительно к байесовской постановке задачи это означает, что начальные условия для системы уравнений (3.156) представляют собой случайный вектор Х(0) с компонентами х(0) =х0,

Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru

с заданным математическим ожиданием Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru и корреляционной матрицей Рх (0).

Предположим, что наземный измерительный пункт (НИП) со­вершает плоское движение по окружности в плоскости, параллель­ной плоскости экватора, с угловой скоростью, равной угловой ско­рости Земли. В таком случае определены экваториальные коорди­наты НИП:

Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru

где Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru — координаты НИП в гринвичской системе координат; Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru — координаты НИП в абсолютной эква­ториальной системе координат;

Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru

Здесь Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru - момент, к которому «привязываются» начальные усло­вия xо; Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru — угловая скорость Земли; Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru —-прямое восхожде­ние гринвичского меридиана на момент Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru ;t и Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru — среднее солнеч­ное время. Модель измерения дальности — простейшая. Считаем, что с наземных измерительных пунктов измеряется дальность р, связанная с координатами ИСЗ и НИП соотношением

Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru

причем измерения осуществляются дискретно в моменты Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru и соп­ровождаются центрированными некоррелированными ошибками Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru со средним квадратичным отклонением Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru , так что можно записать

Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru

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

Обозначим Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru , где Х(0)—искомый «истин­ный» вектор начальных условий движения ИСЗ; Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru где Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru — расчетное значение дальности в момент Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru соответствующее движению ИСЗ из искомых «истинных» начальных условий Х(0). Тогда, линеаризуя зависимость дальности от начальных условий в окрестности точки Х(0), можно записать

Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru

где Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru — расчетное значение дальности в момент Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru ,соответству­ющее движению ИСЗ из начальных условий Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru . Приведенное выражение удобно переписать в виде

Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru


Производные Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru называются баллистическими.

Уравнения (3.158) называются условными, если измерения име­ют различную физическую природу или не являются равноточны­ми. В последнем случае вводят так называемые веса измерений

Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru

причем значение Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru — произвольное, и уравнения (3.158) принима­ют вид

Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru

Как уже указывалось, алгоритм метода наименьших квадратов определяется из условия минимизации суммы квадратов невязок Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru , взятых в случае неравноточных измерений с весами pi.

Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru

Минимизация (3.160) по Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru приводит, как известно (см. (3.144), к следующим уравнениям относительно Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru :

Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru

где Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru —блочная матрица 6×(N+l) вида Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru ;

Р- диагональная матрица (N+l) × (N+1) с элементами рi2;

Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru - вектор (N+1)×1.

Уравнение (3.161) называют нормальным уравнением. Обозна­чим А = НТРН и В = НТР. Элементы akn матрицы А и bk матрицы В можно записать в скалярной форме:

Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru

Где Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru частные производные от измеряемого параметра pi по определяемым начальным условиям. Их удобно записать в форме

Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru

В этом выражении Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru - частные производные от измеряемых . параметров по компонентам вектора состояния Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru , вычисленные в момент Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru - частные производные от компонент вектора состояния ИСЗ Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru по начальным условиям Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru , вычисленные в момент Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru . Существуют различные вычислительные алгоритмы ме­тода наименьших квадратов, отличающихся способом определения частных производных, входящих в (3.163). Используем для опреде­ления этих частных производных метод конечных разностей. Частные производные Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru определяются численным дифференцированием по следующей схеме. Задаются начальньши условиями

Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru

Здесь Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru —априорные расчетные (опорные) начальные условия,, используемые в качестве начального (нулевого) приближения в алгоритме наименьших квадратов. Вектор Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru при вычислении производной Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru формируется следующим образом: Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru , где Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru — априорное среднеквадратичное отклонение компоненты Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru вектора состояния ИСЗ, соответствующее априор­ной корреляционной матрице Px(0). Для остальных компонент Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru ; j — 2...6. При начальных условиях (3.164) численно ин­тегрируются уравнения движения ИСЗ (3.156) и рассчитываются значения Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru вектора состояния ИСЗ на все требуемые моменты i=l, 2, ..., N. Затем принимаются начальные условия

Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru

и при этих начальных условиях точно также рассчитываются зна­чения вектора состояния ИСЗ Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru , i=l, 2,..., N. Частные производ­ные Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru определяются по формуле

Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru

Аналогично рассчитываются частные производные компонент Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru текущего вектора состояния ИСЗ Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru по другим компонентам Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru вектора начальных условий Х(0). Частные произ­водные измеряемых параметров по компонентам текущего вектора состояния могут быть рассчитаны по их явным выражениям на ос­нове рассчитанных текущих значений компонент Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru вектора Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru . Заметим, что с целью уменьшения объема вычислений при опреде­лении производных Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru можно использовать более простые выражения

Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru

где Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru — значение k-й компоненты вектора состояния момент Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru при начальных условиях Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru .

С целью дальнейшего уменьшения затрат машинного времени при прогнозировании состояния ИСЗ на моменты Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru (i=1, 2,..., N) можно вместо численного интегрирования уравнений (3.156) воспользоваться приближенным аналитическим методом по следую­щей схеме. Вначале осуществляется переход от вектора состояния Х(0) с компонентами Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru в декартовых экваториаль­ных координатах к вектору элементов орбиты в следующей после­довательности. Определяются:

большая полуось орбиты

Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru

где

Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru

фокальный параметр

Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru

Эксцентриситет Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru

Наклонение Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru

долгота восходящего узла

Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru

где k — произвольный множитель (это равенство определяет cos Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru );

аргумент широты ИСЗ

Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru

истинная аномалия

Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru

аргумент перицентра

Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru

эксцентрическая аномалия

Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru

время прохождения перигея

Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru

период обращения ИСЗ

Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru

Определим эксцентрическую аномалию Еi, соответствующую» положению ИСЗ на требуемый расчетный момент Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru из уравнения Кеплера:

Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru

Это уравнение обычно решают методом итераций, представляя егов виде

Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru

Затем определяется истинная аномалия:

Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru

Далее определяются долгота, радиус-вектор, радиальная и трансверсальная составляющие скорости ИСЗ для расчетного мо­мента Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru .

Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru

Затем определяются компоненты вектора Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru в экваториальной системе координат:

Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru

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

1. Численно интегрируются уравнения движения (3.155) при начальных условиях Х(0) и определяется опорная орбита ИСЗ, а также компоненты Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru вектора состояния Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru , i=0, 1, 2, ..., N, для моментов Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru , в которые производятся измерения.

2. На каждый момент Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru , i—0, 1, 2,..., N, в который производят­ся измерения Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru , определяется расчетная дальность Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru .

3. Формируются разности Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru . При этом предполагает­ся, что выборка Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru известна. Одновременно иск­лючаются так называемые аномальные измерения, не удовлетворяющие условию

Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru

где Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru — допустимое отклонение измерения от расчетного, ко­торое задается заранее из технических соображений.

4. Для моментов измерений Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru , удовлетворяющих условию (3.165), вычисляются производные от координат Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru текущего вектора состояния Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru по начальным условиям Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru , k,j=1, 2, …, 6; i=0, 1, 2, ..., N.

5. Рассчитываются частные производные Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru от измеряемых параметров по компонентам xki вектора Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru и определяются значе­ния производных измеряемых параметров по начальным условиям Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru .

6. Производится формирование нормальных уравнений Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru .

7. Решается система нормальных уравнений и определяются поправки Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru к начальным условиям Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru .

Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru 8. Формируется новый вектор начальных условий Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru , и решение задачи повторяется вновь с п. 1 до тех пор, пока на r-й итерации не будет выполнено условие

Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru

причем значение Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru задается заранее.

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

Ниже приводятся результаты имитационного моделирования процесса определения орбиты ИСЗ методом наименьших квадра­тов. Схема моделирования в ос­новном соответствует рис. 3.2. Исходные данные при моделирова­нии были приняты следующими. Параметры расчетной (опорной) орбиты ИСЗ: период Т = 22260 с, эксцентриситет е=0,75, наклоне­ние i = 63°, долгота восходящего узла Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru = 84°. При этом экватори-

альные координаты ИСЗ на опорной орбите в момент t0 равны Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru = 892,2 км, Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru = 9877 км, Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru = -1,46 км, Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru = -2,82 км/с, Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru = = 4,27 км/с, Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru = 6,25 км/с. Требовалось уточнить значения компо­нент вектора Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru , Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru , Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru , Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru , Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru , Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru при условии, что они заданы с априорными разбросами Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru = 0,5 км, Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru = 20 м/с. В качестве измеряемого параметра принималась наклон­ная дальность до НИП, измеряемая со случайной ошибкой Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru ), име­ющей среднеквадратичное отклонение 25 м. Число НИП варьирова­лось от одного до трех. Геодезические координаты НИП приведены в табл. 3.1.

Результаты статистического моделирования приведены в табл. 3.2. В ней даны поправки к начальным условиям в зависимости от номера приближения при определении орбиты по данным трех НИП.

Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru

3.3.4. Аналитический алгоритм априорной оценки точности с использованиемметода наименьших квадратов

На основе метода наименьших квадратов может быть построен высокоэффективный с точки зрения быстроты вычисле­ний алгоритм оценки точности определения орбит ИСЗ, близких к круговым. Обычные методы, применяемые для оценки точности, весьма трудоемки, так как требуют вычисления либо баллистиче­ских производных (3.162), если применяется метод наименьших квадратов, либо использования рекуррентных соотношений вида (3.43), (3.44), (3.45) для апостериорной корреляционной матрицы расширенного вектора состояния оцениваемой динамической си­стемы, если для определения орбиты используется байесовский ал­горитм, например, фильтр Калмана или фильтр второго порядка. Для определения орбит измерения производятся на отдельных участках полета ИСЗ, определяемых, в частности, условиями ра­диовидимости ИСЗ, если измерения проводятся с НИПов. При этом на участках радиовидимости НИПы работают в течение огра­ниченных отрезков времени, называемых сеансами траекторных измерений. Поэтому представляет интерес алгоритм априорной

Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru

Рис. 3.16. Орбитальная наклонная сис­тема координат

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

В таком случае вопрос о сходимости оценок к истинным значе­ниям компонент вектора состояния не возникает и можно ограни­читься упрощенными моделями движения ИСЗ и траекторных из­мерений. Можно ограничиться, в частности, использованием лине­аризованной в окрестности опорной орбиты моделью движения ИСЗ. Решению линеаризованной в окрестности опорной экватори­альной орбиты системе уравнений соответствует фундаментальная матрица (3.77). Однако, когда опорная круговая орбита является наклонной, линеаризованная система уравнений движения не до­пускает решения в квадратурах. С целью устранения этой трудно­сти перейдем от экваториальной к орбитальной системе координат, начало которой совпадает с центром масс Земли, оси Охо и Оуо лежат в плоскости опорной орбиты, ось Охо направлена в восходя­щий узел орбиты, ось Ozo нормальна к плоскости орбиты и образу­ет правую тройку с осями Охо и Оуо (рис. 3.16).

Введем в рассмотрение вектор состояния ИСЗ в сферических орбитальных координатах:

Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru

В таком случае, линеаризуя соотношения, связывающие вектор состояния Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru ИСЗ в декартовой экваториальной системе координат с вектором Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru , получаем

Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru

где

Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru

Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru

Здесь Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru — параметры опорной орбиты, Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru - орбитальная угловая скорость ИСЗ; и — аргумент широты в опорном дви­жении.

В терминах компонент вектора Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru можно определить фундамен­тальную матрицу линеаризованной системы уравнений движения ИСЗ. Формально эта матрица, которую обозначим Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru , совпа­дает с фундаментальной матрицей (3.77).

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

Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru

Считаем, что измеряемым параметром является наклонная даль­ность р.

В таком случае фундаментальная матрица, соответствующая вектору Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru , имеет вид

Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru

Нелинейное соотношение, связывающее измерения с вектором состояния оцениваемой системы, после линеаризации в окрестно­сти опорной орбиты принимает вид

Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru

где Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru — отклонение результата измерения от расчетного значе­ния измеряемого параметра, вычисленного на опорной орбите; Hi — матрица 1×10 вида

Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru

Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru - случайная ошибка измерений с характеристиками

Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru

Предположим, что известна диагональная априорная корреля­ционная матрица расширенного вектора состояния Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru . Эта корреляционная матрица может быть учтена при обработке изме­рений по методу наименьших квадратов, исходя из следующих со­ображений. Поскольку нас интересуют не оценки, а их апостери­орные дисперсии, можно условно считать, что наличие априорной корреляционной матрицы Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru 10×10 эквивалентно проведению 10 дополнительных измерений дальности с точностью, характери­зуемой матрицей Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru . В таком случае матрица Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru , характеризу­ющая всю совокупность измерений в обрабатываемом сеансе, име­ет вид

Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru

Где блок Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru представляет собой диагональную матрицу (N+ 1) × (N+ 1) с элементами Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru .

Апостериорная корреляционная матрица вектора Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru в резуль­тате обработки сеанса из (N+1) измерений согласно (3.147) при­мет вид

Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru

Где

Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru

Соотношение (3.169) удобно записать в виде Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru , причем, учитывая структуру (3.167) матрицы Н и (3.166) матрицы Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru , это выражение можно привести к виду

Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru

ИЛИ

Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru

где R — матрица 10×10 вида

Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru

Таким образом, слагаемое R в правой части (3.171) отражает вклад измерений в уточнение оценок компонент вектора Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru . В ре­зультате приходим к следующей процедуре оценки точности, опи­сывающейся рекуррентными лишь по отношению к сеансам соотно­шениями:

Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru

где i, i+1— номера сеансов; Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru — фундаментальная матрица вида (3.166), связывающая момент окончания предшест­вующего сеанса измерений с моментом начала очередного сеанса. Алгоритм расчета по соотношениям (3.172) существенно упро­щается, если предположить, что в пределах одного сеанса частные производные, входящие в матрицу Н, связывающую измеряемые и оцениваемые величины, не зависят от номера измерения, т. е. по­стоянны. В таком случае для элементов матрицы R записываются выражения достаточно простого вида:

Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru

В этих выражениях (N+1) —число измерений в сеансе. Все элементы R11,..., R1010 в конечном счете определяются девятью ко­эффициентами, которые могут быть рассчитаны заранее:

Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru

Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru

Рис. 3.17. Сравнение аналитического алгоритма и имитационного моделирования

Здесь Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru — интервал между измерениями; Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru — время от начала сеанса до j-го измерения.

Аналогичные выражения могут быть получены и для случая, когда измеряются два параметра: дальность р и скорость изменения дальности Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru .

Эффективность рассмотренного алгоритма может быть показа­на на следующем примере. Пусть рассматривается задача оценки точности определения орбиты ИСЗ, находящегося на наклонной круговой 12-часовой орбите, по данным измерений дальности с трех НИПов.

На рис. 3.17 приведены кривые, показывающие результаты рас­четов оценки точности с помощью аналитического алгоритма (кривая 1) и рекуррентных соотношений фильтра Калмана (кри­вая 2). По оси абсцисс отложено время в сутках, а по оси ординат - нормированное cреднеквадратичное отклонение Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru - оценки компоненты Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru вектора состояния ИСЗ. Здесь Определение орбиты искусственного спутника Земли с использованием метода наименьших квадратов - student2.ru - априорное среднеквадратичное отклонение компоненты. Рисунок де­монстрирует близость результатов, полученных различными мето­дами, при этом время проведения расчетов с помощью аналитического алгоритма примерно в 60 раз меньше.

Наши рекомендации