Просмотр сообщений

В этом разделе можно просмотреть все сообщения, сделанные этим пользователем.


Сообщения - andyp

Страницы: 1 [2] 3 4
16
Этот дистрибутив можно взять сайта AstDys, только он без танцев с бубном
работать не будет. Я уже в этом убедился.

17
Я собираюсь попробовать создать программу OrbSearch для расчета и уточнения элементов орбит астероидов.
Основная задача данного приложения - расчет как можно более точных элементов орбиты
на основе большого числа наблюдений, точность должна быть не хуже, чем у OrbFit.
Кроме того, данное приложение должно быть проще по использованию по сравнению
с OrbFit и свободно от ошибок, содержащихся в OrbFit.

Предполагается, что OrbSearch будет написан на C++ с использованием наработок,
которые были созданы в процессе работы над LinOccult.
Результаты работы OrbSearch должны быть тесно интегрированы с LinOccult.
OrbSearch будет создан на платформе UNIX (Linux) с возможностью портирования на другие платформы.
Лицензия GPL GNU Open Source.

Программа будет содержать два основных модуля расчета - определения достоверных
наблюдений из множества доступных и модуля уточнений элементов орбиты по этим наблюдениям.

Однако, соревноваться с OrbFit - задача на порядок сложнее, чем соревнование с
написанной на коленке WinOccult. Поэтому скорее всего, будет необходима разнообразная
помощь от максимально возможного числа людей. Если это реально, я начну разработку
более конкретного плана действий.

Андрей Плеханов

18
Я ставлю ту же дату, которая использовалась в FindOrb для
сравнения результатов.
Больше количество цифр, увы, OrbFit не выводит.
Вот уточнения для момента покрытия:

        Semimajor axis     =   3.11413634253370E+00 AU
        Eccentricity       =   0.097958494963160
        Inclination        =  24.7416299334460 deg
        Long. of node      =  32.5987625215677 deg
        Arg. of pericenter = 254.2884485838171 deg
        Mean anomaly       = 141.1157588755332 deg
        Epoch of elements  : MJD   53760.82903000 TDT (Jan 25, 2006, 19.896720 h)

0.082"    0.078"-152.2

19
Вот результаты вычисления элементов орбиты:

        Semimajor axis     =   3.11426988335481E+00 AU
        Eccentricity       =   0.097916548457838
        Inclination        =  24.7415480631589 deg
        Long. of node      =  32.5986098618435 deg
        Arg. of pericenter = 254.3016678368601 deg
        Mean anomaly       = 132.1672897593710 deg
        Epoch of elements  : MJD   53711.00000000 TDT (Dec  7, 2005,  0.000000 h)

0.099"    0.090"  52.6

20
Вот что получилось:

        Semimajor axis     =   3.49312285725626E+00 AU
        Eccentricity       =   0.047205229885712
        Inclination        =   5.4165710090750 deg
        Long. of node      =  79.4781701082467 deg
        Arg. of pericenter =   2.6478233820814 deg
        Mean anomaly       =  30.2149813640963 deg
        Epoch of elements  : MJD   53741.86300000 TDT (Jan  6, 2006, 20.712000 h)

Ошибка: 0.098"    0.079"  98.9

21
В процессе обработки полученных данных возникли мелкие вопросы.

1) Почему Максим (а вслед за ним Андрей) выбрали эпохой элементов 2006, январь 8.28? Не опечатка ли это? Ведь покрытие ожидается 2006, январь 8.690 - то есть около 16 ч 34.2 м UT.

При проведении расчётов я указывал эпоху элементов именно январь 8.28.

2) Почему у Андрея получилось два набора элементов с разным эллипсом ошибок (http://www.kadar.ru/forum/index.php/topic,206.msg1308.html#msg1308 и http://www.kadar.ru/forum/index.php/topic,206.msg1427.html#msg1427 )? Соответствует ли им один и тот же наблюдательный материал или разный?

3) Нельзя ли "попросить" программу Find_orb выводить бОльшее количество значащих цифр, чем в сообщении http://www.kadar.ru/forum/index.php/topic,206.msg1359.html#msg1359 - там их явно недостаточно?

Какие выводы следуют из рассмотрения представленного материала?

1) Неплохо было бы уточнить... координаты покрываемой звезды.

2) По сравнению с орбитой из файла astorb.dat элементы орбиты Максима "сместили" полосу на 260 км к северу; среднеквадратичная погрешность положения полосы при этом возросла со 141 км до 463 км. Вероятно, программа Find_orb непригодна для наших целей: её точности явно не хватает.

3) По сравнению с орбитой из файла astorb.dat оба набора элементов орбиты (от 9 декабря и от 18 декабря 2005 г.) Андрея "сдвинули" полосу на одно и то же расстояние: 57 км к югу. Погрешность положения полосы в первом случае получилась 129 км, во втором - 136 км.

4) Вероятность покрытия в середине полосы в результате уточнения возросла несущественно: с 8.7% до 9.6%. Следовательно, надо проводить больше астрометрических наблюдений в процессе подготовки к покрытию, причём желательно - на разных инструментах (чтобы уменьшить влияние систематических ошибок инструмента). Обрабатывать полученные снимки следует с максимальной тщательностью (так, как это делает Игорь Измайлов), чтобы "извлечь" из них максимально точную информацию. Желательно проведение исследований искажений, вносимых применяемыми инструментами.

Максим! Завершён ли цикл работ по наблюдению Zhangguoxi?

Если завершён, то не следует ли попытаться отправить материалы также и Стиву Престону - чтобы посмотреть, как ему удастся его обработать (может быть, он сумеет сильнее сократить погрешность)? Нет ли у Вас готового файла, в котором собраны все наблюдения Zhangguoxi для обработки программой Find_orb (в формате MPC) - в окончательном варианте?

Всего хорошего,
Алексей Лосюк

Эпоху для расчетов я взял такую-же, что была в FindOrb для сравнения полученных элементов.
А разные элементы и ошибки получились потому, что элементы рассчитаны для разных эпох.
Наблюдательный материал один и тот же.

22
Я рассчитал элементы орбиты с помощью OrbFit.
Вот что получилось:

     Propagated orbital elements for object 3028 (computed with least squares fit):
        Semimajor axis     =   3.01926921249583E+00 AU
        Eccentricity       =   0.024687685580059
        Inclination        =   9.5071891791671 deg
        Long. of node      = 189.8183951669222 deg
        Arg. of pericenter =  12.0328241475572 deg
        Mean anomaly       = 268.1498521874252 deg
        Epoch of elements  : MJD   53743.28000000 TDT (Jan  8, 2006,  6.720000 h)

0.094"    0.072"  95.9

23
А в OziExplorer это выглядит так. Кстати можно смотреть на Землю с любой точки,
не только с астероида.

26
Хочу уточнть по поводу карт в LinOccult'е.
Я этой проблемой занимался и вот что на настоящий момент есть:
1. Карты полосы покрытия на земном шаре LinOccult рисует, но это глобальные
карты, для наблюдателя не подходят. Причем это 3D карты выполненные
с помощью программы PovRay.
2. Для более детальных карт LinOccult выводит ценральную полосу покрытия
в виде GPS точек, пригодных для использования программой OziExplorer.
Здесь проблема в том, что нужны привязанные к местности карты для
точки наблюдения. Для Москвы и области есть прекрасные двухкилометровки,
для других регионов не знаю.
3. Я нашел несколько геоинформационных систем, которые отображают
различную информацию на картах, но освоить их очень непросто.

27
Да в общем ничего сложного
если обычная программа улучшения орбиты - работы на пару месяцев

Нет, обычную делать неинтересно да и смысла нет. Хотелось бы посоревноваться с
OrbFit и FindOrb за счет более тщательного анализа уже существующих данных.
И сделать программу более удобной в использовании. И согласовать с LinOccult.
Это не на два месяца.

28
Для LinOccult эпоха может быть любой, а для WinOccult надо только
определенную эпоху. Спросите Алексея Лосюка по этому поводу.

Я рассчитал полосу покрытия по элементам из OrbFit:

29
Рассчитал элементы орбиты в OrbFit.
Вот что получилось:

Propagation of orbital elements:
     Propagated orbital elements for object 3028 (computed with least squares fit):
        Semimajor axis     =   3.01933940705500E+00 AU
        Eccentricity       =   0.024676952444124
        Inclination        =   9.5071866357262 deg
        Long. of node      = 189.8185919235370 deg
        Arg. of pericenter =  11.9804586638075 deg
        Mean anomaly       = 262.1370294976466 deg
        Epoch of elements  : MJD   53711.00100000 TDT (Dec  7, 2005,  0.024000 h)

Ошибка:
0.087"    0.068"  95.4

Максим чем вы переводите из .rwo в MPC ?

30
Интересная информация. Я последнее время задумываюсь над тем, не начать ли проект по
написанию программы расчета элементов орбит астероидов аналогичную OrbFit и FindOrb
но более удобную, с использованием улучшенных алгоритмов и более подходящую для предсказаний
покрытий. Что по этому поводу думает астрономическая общественность ?

Страницы: 1 [2] 3 4