Какое программное обеспечение подходит для интегрирования радиальной скорости по времени [закрыто]

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

Проблема в том, что я не знаю, какую программу/программное обеспечение мне следует использовать для построения наиболее подходящей линии, проходящей через точки необработанных данных радиальной скорости по фазе. Более того, мне нужно интегрировать эту линию наилучшего соответствия, а также построить интегрированную кривую.

Может ли кто-нибудь порекомендовать хорошее программное обеспечение, которое я могу использовать?

(Я пробовал интеграцию с программой под названием Logger Pro , но она не увенчалась успехом.)

Похоже, вам нужен какой-то язык программирования, который может справиться с этим, например Mathematica, R, Python, Matlab и т. д. Существует множество вариантов. Это просто зависит от того, к какому из них у вас есть доступ и вам удобно учиться.
Я согласен, вам, вероятно, придется программировать его самостоятельно. Я не знаком с Matlab или R, но для python есть пакет astropy , в котором могут быть некоторые строительные блоки для начала работы, поэтому вам не придется делать все с нуля. На первый взгляд, в нем есть инструменты для подгонки модели и построения графиков.
Используйте правило Симпсона или правило трапеции. Можно было бы тривиально сделать в электронной таблице.
Я голосую за то, чтобы закрыть этот вопрос как не по теме, потому что он касается не астрономии, а численных методов.
Я голосую за то, чтобы закрыть этот вопрос как не относящийся к теме, потому что этот вопрос относится к Stack Exchange рекомендаций по программному обеспечению.

Ответы (1)

Я не уверен, что вы на самом деле пытаетесь сделать. Радиальная скорость обычно не измеряется, она рассчитывается по прогнозируемой скорости. Но если я поверю вам на слово, что ваши необработанные данные — это радиальная скорость, то это все упростит.

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

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

Таким образом, задача интегрирования становится задачей простой геометрии (у вас есть высокий тонкий прямоугольник с треугольником наверху для каждого интервала). Давайте подумаем о звезде с 30-дневным периодом. Это ~ 40 000 минут, поэтому я бы, вероятно, использовал 5 или 8-минутные интервалы ... начиная с минуты 0, какова площадь под кривой от t = 0 до t = 8 минут? затем с 8 до 16, вплоть до 43 200-й минуты. Сумма этих площадей и есть общая площадь. Электронная таблица Excel может легко сделать это.

На самом деле, если вы также используете 5-минутные интервалы и сравните результаты, вы увидите, насколько лучше (надеюсь, НЕ намного) меньший интервал ... аналитическое интегрирование доводится до предела интервала, приближающегося к нулю, но в реальной жизни , численное интегрирование необходимо почти для всех «интересных» задач. Вот 2 ссылки, кстати:

Слишком сложно.