Ich versuche, Kepler-Daten für Eclipsing Binaries zu verwenden, um den Zeitraum und dann andere Parameter wie Masse, Exzentrizität, große Halbachse, Entfernung usw. des Systems zu schätzen. Ich möchte Code in MATLAB schreiben, der FFT verwendet. Die verfügbaren Daten haben die folgenden Spalten:
bjd(Datum) | Phase | raw_flux | raw_err | corr_flux | corr_err | dtr_flux | dtr_err
Müssen die Daten gereinigt werden? Wenn ich die FFT der Daten nehmen möchte, um den Zeitraum zu finden, welche Daten sollte ich verwenden?
Dazu sollten Sie eine Lichtkurve erstellen, ein Diagramm der Helligkeit über der Zeit, um die Daten anzuzeigen. Bei Kepler-Daten ist die Spalte bjd(date) die Uhrzeit in BJD . dtr_flux steht für detrended flux, was bedeutet, dass es "bereinigte" Daten enthalten sollte.
Zuerst würde ich die dtr_flux-Spalte über bjd(time) zeichnen, um eine Lichtkurve zu erhalten. Abhängig von der Art der Eclipse-Binärdatei sollten Sie periodisch einen minimalen und maximalen Einbruch sehen. Die Daten enthalten auch eine „Phasen“-Spalte, was bedeutet, dass die Daten „gefaltet“ wurden, damit Sie jeden Zyklus der Finsternis-Binärdatei übereinander darstellen können. Siehe http://www.southampton.ac.uk/~sdc1g08/BinningFolding.html
Danach können Sie die FFT durchführen, wenn Sie möchten. Denken Sie daran, dass die FFT sowohl bei P als auch bei P / 2 eine Spitze aufweist, da dies eine verdunkelnde Binärdatei mit maximalen und minimalen Helligkeitseinbrüchen ist.
Zephyr