Ableitung der Phase eines pulsierenden Sterns aus der Liste der julianischen Datumswerte der Beobachtungshelligkeit

INT bezieht sich auf die Bodenfunktion . So$\mathrm{INT}(x)$ist die größte ganze Zahl, die nicht überschreitet$x$.

In diesem Fall wird die Funktion verwendet, um den Bruchteil von zu erhalten$\frac{\mathrm{HJD}-T_0}{P}$, also der Teil nach dem Komma. Also wenn zB$\frac{\mathrm{HJD}-T_0}{P}=3.4$, dann$\frac{\mathrm{HJD}-T_0}{P}-\mathrm{INT}\left(\frac{\mathrm{HJD}-T_0}{P}\right)=0.4$.

Sie können Ihren Referenzpunkt frei wählen, wie Sie möchten. Du könntest dir eine aussuchen$\text{HJD}$bei dem der Stern am hellsten oder am dunkelsten oder genau im Durchschnitt ist, oder Sie können den auswählen$\text{HJD}$der ersten Beobachtung, oder Sie könnten sogar festlegen$T_0 = 0$. Wenn zwei Leute unterschiedlich auswählen$T_0$, dann unterscheiden sich die gefundenen Phasen um einen konstanten Wert ( modulo 1).

Es ist bequem einzustellen$T_0$bei oder vor dem kleinsten$\text{HJD}$in Ihren Beobachtungen, damit alle$\text{HJD} - \text{T_0}$sind nichtnegativ, weil die$\text{INT}$-ähnliche Funktionen, die von Taschenrechnern, Computersprachen und Softwarebibliotheken bereitgestellt werden, runden nicht alle negativen Zahlen in die gleiche Richtung. Viele runden alle Zahlen gegen Null, damit$+3.4$Erträge$+3$(wie wir wollen) und$-3.4$Erträge$-3$, aber hier wollen wir$-3.4$nachzugeben$-4$, sonst bekommt man negative Phasen vorher$T_0$und positive Phasen danach$T_0$.

Wenn Ihre Helligkeitsbeobachtungen zeitlich gleichmäßig verteilt sind, können Sie eine schnelle Fourier-Transformation verwenden , um die ungefähre Periode und Phase zu ermitteln.