Интерактивные Web-формы созданы для отладки расчетов, связанных с прогнозом
движения космического аппарата. Расчет положения планет, плотности атмосферы
и т.д. и т.п. смотрятся исключительно со стороны воздействия соответствующих сил
при движении в космическом пространстве.
В основу программ расчета положена библиотека ESTK (Earth Space ToolKit),
разработанная автором сайта с использованием накомпленного в
ИПМ им. М.В. Келдыша опыта решения
подобных задач. ЕSTK написана на Си и работает в операционных системах
Windows и Linux (+UNIX системы).
ESTK - не является FREE-software, т.к. используется в ряде
задач комерческого характера. Тем не менее, в разделе
ссылки можно найти ссылки на свободное программное
обеспечение. Отмечу, что его качество оставляет желать лучшего, да и зоопарк
языков программирования и операционных систем огромен.
По мере возможности и наличия времени автор будет добавлять интересные расчеты
и выставлять их в интерактивном виде. Прошу прощения за топорность форм, т.к.
написание скриптов, HTML и взаимодействие серверов для организации расчетов -
это не моя работа :)
Библиотека ESTK (Earth Space ToolKit) Ядро 41.2 [17 декабря 2020 г.]
На FTP
ftp://ssd.jpl.nasa.gov/pub/eph/planets/ содерджатся массивы ежегодников
планет. На этом же FTP можно взять исходные коды программ для расчета эфемерид
планет.
На
http://ssd.jpl.nasa.gov/horizons.cgi#top находится Web-форма для запроса эфемерид
ежегодника 405 на эфемеридное время.
В данной форме осуществляется расчет положения (км) и скорости (км/с) планет
на время UTC в системе координат (СК) J2000. Выходные данные этой формы
протестированы на совпадения с расчетами в JPL.
Перейти к web-форме планетарий
На
http://tai.bipm.org/iers/conv2003/conv2003_c5.html
содержатся коды программ перевода систем координат из J2000 в WGS-84 и наоборот.
Приведено 3 модели: IAU200A, IAU200B, IAU1980.
Программы перевода используют следующие обновляемые таблицы:
http://maia.usno.navy.mil/ser7/finals.all
http://maia.usno.navy.mil/ser7/finals2000A.data
http://maia.usno.navy.mil/ser7/tai-utc.dat
На
http://hpiers.obspm.fr/eop-pc/products/matrice/matrice.html находится
Web-форма для запроса матрицы перевода модели IAU2000A на заданное время UTC.
В данной форме осуществляется расчет матриц перевода из J2000 в WGS-84 и
наоборот на время UTC. Выходные данные этой формы
протестированы на совпадения с расчетами в форме, указанной выше.
Матрица записывается таким образом, что перевод осуществляется умножением
соответствующей матрицы на вектор-столбец.
Полезная ссылка на книгу
В.Е. Жаров Сферическая астрономия М: 2002.
Перейти к web-форме iers
Плотность атмосферы (кг/м^3) расчитывается по следующим моделям:
На
http://omniweb.gsfc.nasa.gov/vitmo/msis_vitmo.html находится Web-форма для
запроса параметров атмосферы MSISE-1990 на заданное время UTC, положение
объекта и по заданным индексам солнечной активности:
F10_7 | Индекс солнечной активности, равный
плотности потока радиоизлучения Солнца на волне длиной 10.7 см, 10^-22Вт/(м^2Гц) |
F81 | Средне-взвешенное за 81 сутки
значение индекса F10_7 , 10^-22Вт/(м^2Гц) |
Kp | Квазилогарифмический,
планетарный среднесуточный индекс геомагнитной активности |
В данной форме осуществляется расчет плотности атмосферы (кг/м^3) для
всех моделей на время UTC и по заданному вектору положения (км) в WGS-84
.
Перейти к web-форме атмосфера
Статистика доступа к solarinex.txt.
Для расчета атмосферы требуется наличие индексов
солнечной активности. Помимо истории, требуется и прогноз.
В Баллистическом центре ИПМ
эта проблема решена. Скачиваются данные и строится прогноз. На моем сайте
попытаюсь предоставить
продукт, включающий исторические данные с 2002 года
и прогноз с текущего дня на 27 суток.
На
ftp://ftp.swpc.noaa.gov/pub/27DO.txt дается 27-дневный прогноз:
F10_7 | Индекса солнечной активности, равного
плотности потока радиоизлучения Солнца на волне длиной 10.7 см, 10^-22Вт/(м^2Гц);
|
A | Planetary A Index; |
Kp | Квазилогарифмического,
планетарного среднесуточного индекса геомагнитной активности. |
На
ftp://ftp.swpc.noaa.gov/pub/indices/quar_DSD.txt дается ежеквартальный
индекс солнечной активности. (Реальные, не прогнозные данные)
Файлы с индексами ежедневно скачиваются м обрабатываются (делается собственный
прогноз),
результатом является файл со значениями
F10_7 | Индекс солнечной активности, равный
плотности потока радиоизлучения Солнца на волне длиной 10.7 см, 10^-22Вт/(м^2Гц>)
|
F81 | Средне-взвешенное за 81 сутки
значение индекса F10_7 , 10^-22Вт/(м^2Гц) |
Kp | Квазилогарифмический,
планетарный среднесуточный индекс геомагнитной активности |
В файле есть историческая часть и прогнозная. Файл обновляется раз в день.
Файл-результат можно скачать по ссылке с моего сайта:
ftp://ftp.kiam1.rssi.ru/pub/gps/solar/solarinex.txt.
На
ftp://ftp.ngdc.noaa.gov/STP/SOLAR_DATA/ находятся
индексы солнечной активности. (сам пока с ними не разбирался)
На
http://www.sidc.be/archive
Находятся измерения индексов от "SIDC - Solar Influences Data Center".
В БЦ ИПМ
эти файлы накапливаются. Сводный файл представлен на
ftp://ftp.kiam1.rssi.ru/pub/gps/solar/alldata.meu формат файла представлен ниже
1 столбец DATE - Дата UTC
2 столбец SUNSPOT - Количество пятен на Солнце
3 столбец WOLF - Число Вольфа по данным обсерватории Катания (IT9)
4 столбец F10_7 - Плотность потока излучения на волне длиной 10.7 см
5 столбец AK_H - A-индекс по данным обсерватории Chambon-la-Foret (F)
6 столбец AK_W - A-индекс по данным обсерватории Wingst (DL)
7 столбец AP - Планетарный А-индекс Ap
8 столбец ISN - Вчерашнее значение числа солнечных пятен (International Sunspot Number)
9 столбец ST - Количество станций наблюдения ISN
Символы ### означают отсутствие данных
Перейти к web-форме индексы