Здравствуйте, недавно столкнулся с такой проблемой. Нужно проверить работоспособность алгоритма наведения работающего на данных GPS и, желательно, на метеорологических спутниках, но о них нашел лишь данные в формате TLE. Так вот, подскажите, пожалуйста, незнающему , как нам эти данные преобразовать в нужный формат, чтобы проверить работоспособность алгоритма. Нужны такие вот данные: синусная поправка радиуса орбиты отличие главного движения от расчетного средняя аномалия косинусная поправка аргумента долготы эксцентриситет синусная поправка аргумента широты квадратный корень главной полуоси опорное время привязки эфемерид косинусная поправка к углу наклонения долгота восходящего узла орбитальной плоскости синусная поправка к углу наклонения угол наклонения косинусная поправка радиуса орбиты аргумент перигея скорость изменения прямого восхождения скорость изменения угла наклонения оценка групповой дифференциальной задержки поправка времени, мс поправка времени, мс/мс2 поправка времени, мс/мс
либо если кто знает, как из этих TLE можно узнать положение спутника в определенный момент времени, либо на промежутке.
Начинать можно сразу с Frequently Asked Questions for Revisiting Spacetrack Report #3 AIAA 2006-6753 ссылка на эту страничку идёт с Revisiting Spacetrack Report #3 AIAA 2006-6753 Собственно, саму книгу Fundamentals of Astrodynamics and Applications by David Vallado можно найти в интернете, правда, только первую редакцию 1997 года. См. также The OrbitTools Libraries NORAD SGP4/SDP4 Implementations in C++ and C# by Michael F. Henry З.Ы. Вы всё-таки посмотрите, как правильно по-русски называются термины, о которых вы пишите в первом посте.
Помню, что на сайте Вадима Чазова было что-то про преобразование TLE, но не мог вспомнить где. Сейчас нашёл. На страничке Нечто из астрометрии и небесной механики есть ссылка на Прогноз орбитального движения космического аппарата и в описании сказано: "Программная реализация алгоритма. Модуль unfornor.pas. procedure FromNoradToDesCart. procedure FromDesCartToNorad"
Здравствуйте ещё раз, но вдруг кто-то сможет подсказать именно с этой задачей, без вмешательства SGP4 ? Используя алгоритм GPS http://www.gps.gov/technical/icwg/IS-GPS-200H.pdf стр. 159 для метеоспутника NOAA 19 никак не получается верный ответ положения спутника и углов к нему. Вот пример tle файла: NOAA 19 1 33591U 09005A 16031.48171577 .00000103 00000-0 81022-4 0 9997 2 33591 99.0223 347.2628 0013044 245.7495 114.2314 14.12033971359709 В формулу отсюда лишь беру(с переводом в полукруг): io 9-13 0.55012388888888888888888888888889 угол наклонения OMEGAo 18-25 1.9292377777777777777777777777778 долгота восходящего e 27-33 0.013044 эксцентриситет w 35-42 1.365275 аргумент перигея Mo 44-51 0.63461888888888888888888888888889 средняя аномалия и отдельно квадратный корень главной полуоси=2690,9106265351883663254695944981 м1/2 остальные же значения беру за нули, т.к. такие явные поправки нам не нужны,да и неоткуда их брать может кто-то подсказать, что я не так делаю?
Вас только аватарка от моей ругани спасает Я понимаю, что фраза для вас мало что значит. Но не это сейчас главное. Хорошо. Начальные условия задачи вы привели. Но что с чем сравнивали? Что приняли за эталон и что получили вы? На какой интервал времени строили прогноз? Как получили значение большой полуоси орбиты? Кеплеровы элементы орбиты. Это как сказать. Это счислимые величины.
Сравнивал со значениями из Orbitron и сайта http://www.satview.org/forec.php?sat_id=33591U&dmy=01022016&hms=1433 Ну вот их значения я как раз и принимал,как те,к которым нужно стремиться. Рассчитывал на день вперед. То есть данные получены 31.01.16 в 11:33 Считаю на то же время, но уже 1ого февраля. У Орбитрона: аз 335 и ел -35,4 У меня: -18.2879 66.7167 насчет большой оси , есть значения его основное на сайте, что выше вписал, но и просто складывал радиус Земли + расстояние до орбиты спутника (870км) насчет того, что данные нельзя использовать в других формулах, то есть теперь можно считать лишь по SGP4/SDP4 ?( но ведь в GPS нам эти поправки нужны лишь для более точного определения (до метра),а тут это не важно.
Если у вас есть среднее движение, то значит есть и большая полуось. Третий закон Кеплера. В нашем случае — n2a3 = μ Пятая строчка в Table 30-II. Elements of Coordinate System (part 1 of 2) стр. 163. Не обязательно. Можно перейти от осреднённых элементов к оскулирующим (osculor в переводе с латыни – поцелуй. Ох уж эти французы!) и считать, как вашей душе угодно. Как это сделать, написано у Вадима Чазова. Ссылку приводил выше. Так какая нужна точность? Может быть вам можно будет ограничиться только вековыми возмущениями от второй зональной гармоники J2? --- Сообщения объединены, 1 фев 2016, Оригинальное время сообщения: 1 фев 2016 --- аз – азимут, а ел – высота? Так может быть дело не в орбите, а в переходе от инерциальной СК (или гринвичской СК) к топоцентрической?
