Эфемеридный расчет для заданного момента времени и конкретного спутника

Страницы работы

3 страницы (Word-файл)

Содержание работы

                        Задание 4 (последнее , 4 часа). Эфемеридный расчет

                                                      Общие замечания

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

            В рамках Задания 4 студенты знакомятся с содержанием jps-сообщения GE (эфемеридная информация) и самостоятельно выполняют эфемеридный расчет для заданного момента времени и конкретного спутника.

                                      Порядок выполнения Задания следующий.

1. На момент времени __часов __ минуты 00 секунд (Москва) вычислить координаты i-го спутника, где i определяется из таблицы:

№ ЭВМ

1

2

3

4

5

6

7

8

9

10

   I

2

3

15

16

18

21

22

26

29

31

2. Рассчитать дальность до спутника и результат сравнить с соответствующим измерением псевдодальности (в метрах).

3. Отличие результата вычислений и измерений объяснить.

4. Сформировать дифференциальную поправку к псевдодальности.

Замечание 1. Расчеты рекомендуется выполнить в среде MatLab (либо любой другой на языках Pascal, C++, Basic). На точности вычислений не экономить!

Замечание 2.  Определения эфемерид

                M0   Средняя аномалия на время привязки

                Dn   Отличие среднего движения от расчетного значения

                   e   Эксцентриситет

            (A)1/2  Квадратный корень большой полуоси

    (OMEGA)0   Долгота восходящего узла орбитальной плоскости на недельную эпоху ()

                  i0   Угол наклонения на время привязки

                  w   Аргумент перигея

OMEGADOT   Скорость изменения прямого восхождения ()

            IDOT   Скорость изменения угла наклонения

                Cuc     Амплитуда косинусной гармонической поправки к аргументу широты

                Cus     Амплитуда синусной гармонической поправки к аргументу широты

                Crc     Амплитуда косинусной гармонической поправки к радиусу орбиты

                Crs     Амплитуда синусной гармонической поправки к радиусу орбиты

                Cic     Амплитуда косинусной гармонической поправки к углу наклонения

                Cis     Амплитуда синусной гармонической поправки к углу наклонения

                 toe   Опорное время привязки эфемерид

            IODE   Идентификатор набора данных (эфемериды)

Замечание 3. Элементы временных вычислений.

Время излучения сигнала

                                                 t=tGPS-tRANGE,

где tGPS – GPS-время из сообщения GT (в секундах),

      tRANGE – псевдодальность (сек) из сообщения RC.

Cдвиг временной шкалы спутника

                                               ,

            af0, af1, af2 и toc, содержатся в сообщении GE.

Системное время GPS (без учета поправок на релятивистские эффекты и групповую задержку)

                                               tS=t-∆tSV

Замечание 4. Элементы координатных вычислений.

Большая полуось

Расчетное среднее движение - радиан/секунда

tk = tS - toe

Время от опорной эпохи эфемерид*

n = n0 +Dn*π

Скорректированное среднее движение

Средняя аномалия

Уравнение Кеплера для эксцентрической аномалии (может быть решено итерациями) итерациями. Для 0-й итерации Ek=Mk)

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

Fk = nk + ω*π

Аргумент широты

Возмущения второй гармоники

duk = Cus sin 2fk + Cuc cos 2fk

Коррекция аргумента широты

drk = Crc cos 2fk + Crs sin 2fk

Коррекция радиуса

dik = Cic cos 2fk + Cis sin 2fk

Коррекция наклонения

uk = Fk + duk

Скорректированный аргумент широты

Скорректированный радиус

Скорректированное наклонение

Координаты в орбитальной плоскости

Скорректированная долгота восходящего узла

Геоцентрические фиксированные (ECEF) координаты

* t системное время GPS на момент передачи, то есть системное время, скорректированное на время прохождения сигнала (дальность/скорость света). Кроме того, параметр tk должен выражать действительную полную разность между временем t и временем эпохи toe и учитывать начало/конец перехода недели, то есть, если tk  больше 302 400 секунд, то следует вычесть из tk 604 800 секунд, а если tk меньше минус 302 400 секунд, следует добавить к tk 604 800 секунд.

 X=xkcos(tRANGE)+ yksin(tRANGE)

 Y=-xksin(tRANGE)+ ykcos(tRANGE)                         ECEF-координаты с учетом вращения Земли

 Z=zk

Замечание 5. Расчет дальности и дифференциальной поправки.

Дальность

R=((X-XO)2+(Y-YO)3+(Z-ZO)2)1/2, где XO, YO, ZO – координаты приемника.

Дифференциальная поправка

                     ΔR=R- (tRANGE + ∆tSV – tNAVT –TGD + )*c,

где tRANGE – псевдодальность (сек) из сообщения RC,

      tNAVT  – оценка сдвига временной шкалы приемника из сообщения ТО,

      TGD    – поправка групповой задержки из сообщения GE,

     – релятивистская поправка (в секундах),

Замечание 6. Константы.

m - универсальный гравитационный параметр WGS-84, равный (3.986005 x 1014 м32);

c  - скорость света в вакууме, равная (2.99792458 x 108 м/с);

   

π= 3.1415926535898;

= 0.000072921151467 с-1.

Содержание отчета

1. Исходные данные.

2. Алгоритм эфемеридного расчета с выводом всех промежуточных вычислений.

3. Объяснение результатов.

Литература

1. Конспект лекций.

2. ICD-200 (основной документ, регламентирующий работу GPS), интернет.

Похожие материалы

Информация о работе