Estoy tratando de usar Kepler Data for Eclipsing Binaries para estimar el período de tiempo y luego otros parámetros como la masa, la excentricidad, el semieje mayor, la distancia, etc. del sistema. Quiero escribir código en MATLAB que usará FFT. Los datos disponibles tienen las siguientes columnas:
bjd(fecha) | fase | raw_flux | raw_err | corr_flujo | corr_err | flujo_dtr | dtr_err
¿Los datos requieren alguna limpieza? Si quiero tomar la FFT de los datos para encontrar el período de tiempo, ¿qué datos debo usar?
Para esto, debe crear una curva de luz, un gráfico de brillo a lo largo del tiempo, para ver los datos. Para los datos de Kepler, la columna bjd(fecha) es la hora en BJD . dtr_flux significa flujo sin tendencia, lo que significa que debe contener datos "limpios".
Primero, trazaría la columna dtr_flux sobre bjd (tiempo) para obtener una curva de luz. Dependiendo del tipo de binario eclipsante, debería ver una caída mínima y máxima que se produce periódicamente. Los datos también contienen una columna de "fase", lo que significa que los datos se han "doblado" para permitirle trazar cada ciclo del binario eclipsante uno encima del otro. Consulte http://www.southampton.ac.uk/~sdc1g08/BinningFolding.html
Después de eso, puede tomar la FFT si lo desea. Recuerde que dado que este es un binario eclipsante con caídas de brillo máximas y mínimas, la FFT tendrá un pico tanto en P como en P/2.
céfiro