Qué software es apropiado para integrar la velocidad radial con respecto al tiempo [cerrado]

Soy un estudiante de secundaria y actualmente estoy tratando de usar el método Baade-Wesselink para derivar de forma independiente la relación período-luminosidad de las variables cefeidas. Noté que para obtener la variación del radio lineal de una cefeida, necesito integrar su velocidad radial durante su período.

El problema es que no sé qué programa/software debo usar para trazar la línea de mejor ajuste que pasa por los puntos de los datos sin procesar de la velocidad radial sobre la fase. Aún más, la cuestión es que necesito integrar esa línea de mejor ajuste y trazar también la curva integrada.

¿Alguien puede recomendar un buen software que pueda usar?

(He intentado la integración con un programa llamado Logger Pro pero no tuvo éxito).

Parece que desea algún tipo de lenguaje de programación que pueda manejar esto, como Mathematica, R, Python, Matlab, etc. Hay numerosas opciones. Solo depende de a cuál tenga acceso y se sienta cómodo aprendiendo.
Estoy de acuerdo, probablemente tendrás que programarlo tú mismo. No sé Matlab o R, pero para python está el paquete astropy que podría tener algunos componentes básicos para que comiences, así que no tendrás que hacer todo desde cero. De un vistazo rápido, tiene herramientas de trazado y ajuste de modelos.
Utilice la regla de Simpson o la regla del trapecio. Podría hacerse trivialmente en una hoja de cálculo.
Voto para cerrar esta pregunta como fuera de tema porque no se trata de astronomía, se trata de métodos numéricos.
Estoy votando para cerrar esta pregunta como fuera de tema porque esta pregunta pertenece al Intercambio de pila de recomendaciones de software.

Respuestas (1)

No estoy seguro de lo que realmente estás tratando de hacer. La velocidad radial generalmente no se mide, se calcula a partir de la velocidad proyectada. Pero si le tomo la palabra, que sus datos sin procesar son la velocidad radial, entonces eso simplifica las cosas.

Si conoces el período, PAG , y sabes el v Y si conoces la relación teórica entre ellos v = F ( pag ) entonces puede usar mínimos cuadrados lineales para calcular los parámetros ajustables de F ( pag ) al trazar v (medido) vs v (calculado).

Cualquier ecuación de "mejor ajuste" es capaz de interpolar valores de v así que todo lo que necesita es dividir el período en pequeños intervalos y aproximar el área bajo la curva usando sus dos puntos finales (rara vez es necesario usar ajustes parabólicos que también usan los puntos medios de esos pequeños intervalos).

Por lo tanto, el problema de la integración se convierte en uno de geometría simple (tiene un rectángulo alto y delgado con un triángulo encima para cada intervalo). Pensemos en una estrella con un período de 30 días. Eso es ~40,000 minutos, así que probablemente usaría intervalos de 5 u 8 minutos... comenzando en el minuto 0 ¿cuál es el área bajo la curva de t=0 a t=8 minutos? luego de 8 a 16, hasta el minuto 43,200. La suma de esas áreas es el área total. Una hoja de cálculo de Excel puede hacer eso fácilmente.

En realidad, si también usa intervalos de 5 minutos y compara los resultados, verá cuánto mejor (con suerte NO mucho) es un intervalo más pequeño... la integración analítica es esto llevado al límite del intervalo cercano a cero, pero en la vida real , la integración numérica es necesaria para casi todos los problemas 'interesantes'. Aquí hay 2 enlaces, fwiw:

Demasiado complicado.