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

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


Сообщения - andyp

Страницы: [1] 2 3 4
1
Вопрос к автору программы. Проект заброшен? Планируется дальнейшее развитие? <tlen kb портирование последней версии под Windows?

Проект не заброшен, но пока не развивается,
я использую для расчета покрытий.
В принципе я реализовал то, что хотел.
Если есть желающие добавлять новые фичи -
проблем нет, проект открытый.

2
Выложил первую версию программы - OrbSearch 1.0.0 Beta.
Исходники можно найти на моей страничке:

http://andyplekhanov.narod.ru

и на

http://sourceforge.net/projects/orbsearch

Начальная орбита вычисляется с помощью метода Гаусса по трем наблюдениям,
а затем уточняется с помощью нахождения минимума функции отклонения Simplex методом.
В отличии от известных программ упор делается на статистическую обработку и вычисление
неопределености обриты путем вычисления нескольких вариантов.

В программе используются библиотеки из LinOccult и GNU scientific library.

Вот несколько сравнений результатов с OrbFit:

Астероид 49:

OrbSearch (best fitted):

  Epoch (y m d)             2008/05/25.25 (54611.25)
  Semi-major axis (a)         3.09061843 AU
  Eccentricity (e)            0.23131746
  Inclination (i)             3.18124886 deg
  Long. of ascending node   286.13984314 deg
  Argument of perihelion    110.03371310 deg
  Mean anomaly (M)        240.29977619 deg

OrbFit:

        Semimajor axis     =   3.09061933737767E+00 AU
        Eccentricity       =   0.231307001433158
        Inclination        =   3.1814838945372 deg
        Long. of node      = 286.1301224096144 deg
        Arg. of pericenter = 110.0438989194008 deg
        Mean anomaly       = 240.3045776164618 deg
        Epoch of elements  : MJD   54611.28400000 TDT (May 25, 2008,  6.816000
h)

Астероид 165:

OrbSearch (best fitted):

  Epoch (y m d)             2008/05/20.00 (54606.00)
  Semi-major axis (a)         3.12572643 AU
  Eccentricity (e)            0.08338424
  Inclination (i)            11.23946034 deg
  Long. of ascending node   302.65132928 deg
  Argument of perihelion    350.83658219 deg
  Mean anomaly (M)        308.39276761 deg

OrbFit:

        Semimajor axis     =   3.12571555818036E+00 AU
        Eccentricity       =   0.083388649509875
        Inclination        =  11.2393034099170 deg
        Long. of node      = 302.6507359851024 deg
        Arg. of pericenter = 350.8433943595194 deg
        Mean anomaly       = 308.4035998068666 deg
        Epoch of elements  : MJD   54606.09000000 TDT (May 20, 2008,  2.160000
h)

3
Основной каркас программы создан.
Реализован модуль вычисления предварительной
орбиты по трем точкам.
Результат совпадает с OrbFit.

4
Вышла новая версия программы LinOccult 2.0.0 Alpha.
Скачать ее можно здесь:

http://sourceforge.net/projects/linoccult/

или здесь:

http://andyplekhanov.narod.ru/soft/soft.htm

Основные изменения:
- Все библиотеки и само приложение переделаны под "чистый" C++ с использованием namespace, потоков и STL
- Убран устаревший код
- В математическую библиотеку вынесена работа с векторами и матрицами, интегрирование и полиномы Чебышева
- Изменена структура библиотеки ввода-вывода из/в различные астрономические форматы.
Добавлены новые форматы - для использования в OrbSearch
- Выделены две новые библиотеки: для работы с различными астрономическими базами данных и
астрономическими алгоритмами
- Исправлены некоторые мелкие ошибки
- Доделаны или улучшены некоторые процедуры

5
А исходник у этой программы есть ?
Хотелось бы использовать его в  OrbSearch.

6
Если надо очень точно посчитать, можно посмотреть исходники LinOccult:
http://astro-archive.prao.ru/soft/soft.php


7
Хороший сборник.
Жаль не могу забрать экземплярчик.

8
A mne Preston prislal UCAC2 po pochte.
Po obichnoi pochte ne po electronnoi.

9
Если делать все по-взрослому, то надо печатать сборник докладов.

10
Мой друг Женя Малахов портировал LinOccult
под Windows. Эту версию можно найти на моем сайте:

http://andyplekhanov.narod.ru/soft/linoccult/linoccult_1.1.3.zip

Перед использованием почитайте ReadMeRus.pdf.
Если Window будет показывать кривые шрифты
в документации, попробуйте ее распечатать.

11
Я думаю, что в данном случае ошибся Стив.
Вот что он мне написал:

You found a bug in my code :)  Thanks...

Obscode 809 is La Silla. This is one of the sources that I treat as "high
accuracy".  But my code should only reduce the rms for specific instances
and it was reducing the rms for ALL La Silla observations... Opps!

I have fixed this now so it shouldn't happen in future orbit fits.

Поэтому никто это покрытие и не наблюдал.

Андрей

12
А координаты наблюдательной точки будут какие ?

13
Надо сделать как правильно.
Проблема в том, что до 1972 года время наблюдений
определялось в UT1, а после в UTC.
Надо эти времена переводить в ET(TDT).
Как это делается для UTC я знаю, а для UT1 нет.
Я даже не знаю в чем отличие UT1 и UTC.

14
Я думаю, что наблюдения должны фильтроваться по следующему принципу

1) Улучшаем орбиту по всем доступным наблюдениям.
2) Определяем максимальную невязку.
3) Если эта невязка больше некоторого порога, то исключаем это наблюдение (не все наблюдения, которые больше порога, а только одно наблюдение).
4) Снова, улучшаем орбиту и переходим к п.2.

В общем, это достаточно стандартная процедура фильтрации наблюдений.


Да, это стандартный метод.
Но он не всегда подходит.

Тут проблема в том, что прежде чем улучшать орбиту, надо ее вычислить.
А для вычисления надо выбрать три наблюдения (какие ?) из множества доступных.

У меня немного другая идея - попробовать определить несколько орбит
по различным комбинациям трех наблюдений и из всех доступных нам наблюдений выбрать те,
которые попали на большее число этих орбит.

15
Вот есть задача разобрать как работают некоторые процедуры
перевода времени
UT1 -- TDT -- TAI -- UTC
в OrbFit.

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