Welche Daten sollte ich verwenden, um die schnelle Fourier-Transformation (FFT) durchzuführen, um den Zeitraum für die Verfinsterung von Binärdateien zu finden?

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?

Sofern Sie nicht versuchen, mit FFT oder etwas Ähnlichem zu üben, ist Ihre vorgeschlagene Methode zum Ermitteln der Periode übertrieben. Das bekommt man direkt aus dem BJD und den Phasendaten. Eine FFT der Flussdaten ist nicht erforderlich.

Antworten (1)

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.

Also habe ich die Handlung von bjd vs dtr_flux als imgur.com/gallery/3ZFNW bekommen . Ich habe versucht, die FFT davon zu nehmen, aber alle Werte kommen imaginär, außer dem Maximalwert, der als 6,0836e + 04 kommt. Was ist hier das Problem?
@ruskin23 Das ist zu erwarten . Es klingt für mich so, als sollten Sie zurückgehen und genau studieren, was eine FFT ist und wie sie funktioniert, um Ihre Ergebnisse besser zu verstehen. Darüber hinaus sollten Sie wahrscheinlich darauf verzichten, die Lücke in Ihren Daten auszugleichen, da dies Ihre Ergebnisse bis zu einem gewissen Grad verfälschen wird.
Die FFT geht davon aus, dass die verwendeten Samples in regelmäßigen Abständen abgetastet werden. Es funktioniert nicht mit unregelmäßig abgetasteten Daten.