Este fin de semana toca, una vez más, práctica numérica de Fundamentos de Astrofísica. Una de las tareas que hay que hacer es, a partir de la curva de velocidad radial de una estrella binaria espectroscópica, hallar la velocidad del centro de masas.
Según el método de Lehmann-Filhes, lo que hay que hacer es trazar una linea horizontal en la gráfica de velocidad radial, de modo que las áreas encerradas por la curva por encima y por debajo de esta linea horizontal sean iguales. Parece trivial, pero en una licenciatura como la nuestra (CC. Físicas en la UCM) no nos enseñan a manejar programas de cálculo numérico lo suficiente como para resolver este tipo de problemas.
Yo he seguido un proceso iterativo. En primer lugar, he puesto los datos de la curva en una tabla, los he representado, y he usado la función de integración de Qtiplot (un programa clon de Origin) para realizar la integral de la curva por el método del trapecio.
Como se puede ver, el valle de la curva de velocidad radial es mucho más extenso que la cresta, y es por lo que al final las áreas suman negativo. La estrategia consiste en desplazar la curva hacia arriba hasta que la integral de cero, y la distancia que lo hayamos desplazado, cambiada de signo, equivaldrá a la velocidad del centro de masas en km/s. Es decir, si para que la integral se anule tenemos que desplazar la gráfica 8 km/s hacia arriba, la velocidad del centro de masas será -8 km/s.
En este caso, el desplazamiento adecuado fue de -13.3 km/s.
Son bienvenidos comentarios y sugerencias sobre otras formas de hacer el cálculo.
Una preguntilla…eres de la Complutense?Porque si es así, imagino que vas a clase conmigo xD
Andaba mirando por internet más información del método de Lehmann-Filhes y me he encontrado con tu blog. Un placer conocer a alguien que anda liado con las mismas cosas que un servidor!
Sí claro, la práctica que yo estoy haciendo es la misma que estás haciendo tú.
Puse esto en el blog precisamente con la esperanza de que alguien, de ésta o futuras generaciones, lo encontrara y le fuese útil.
Me parece que tambien yo voy a la misma clase. Yo fui a preguntarle al susodicho profesor cómo podía hallar el área y me dijo que representase los puntos en papel milimetrado y luego contase cuadritos…no comment
Al final me lo hizo mi novio con Autocad ;D