Berechnen von Periodenänderungen in ungleichmäßig abgetasteten Zeitserien in R (oder Matlab)

Question

Die Überschrift sagt alles, ich versuche verzweifelt, einen Weg zu finden, um die Periode einer Zeitreihe zu berechnen, die ungleichmäßig abgetastet wird. Ich habe versucht, für die Zeiten, in denen es keine Daten gibt, eine gleichmäßig abgetastete Zeitreihe mit NAs zu erstellen, aber es gibt einfach zu viele NAs für irgendeine Imputationsmethode, um einen vernünftigen Job zu machen. Das Hauptproblem besteht darin, dass die Abtastzeiten viel weiter voneinander entfernt sind als die durchschnittliche Periode (SEHR ungefähr 0,5), was nur bei Anwendung der Periodenfaltung offensichtlich ist. Da ich nach einer kleinen Änderung in der Periode suche, kann ich die Sampling-Zeiten nicht runden.

Time-folded period

Hier ist eine Probe der Daten:

HJD(time) Mag err 
2088.91535 18.868 0.078 
2090.87535 19.540 0.165 
2103.92958 18.704 0.040 
2104.94812 19.291 0.098 
2106.84596 18.910 0.066 
... 
4864.56170 18.835 0.061 

Der Datensatz etwa 650 Zeilen.

Ich habe fast eine Woche damit verbracht, mein Problem zu googeln und nichts hat noch geholfen, also würden alle Ideen sehr geschätzt! Ich habe auch etwas Erfahrung mit Matlab, also wenn es möglich ist, es mit Matlab anstatt R zu machen, wäre ich auch damit zufrieden.

Answer 1

Ich glaube nicht, dass es eine Möglichkeit gibt zu tun, was Sie wollen. Die Nyquist-Shannon theorem besagt, dass Ihre durchschnittliche Abtastfrequenz mindestens doppelt so hoch sein muss wie die Häufigkeit der Ereignisse, die Sie erfassen möchten (grob gesagt).

Wenn Sie also Informationen aus Ereignissen mit einer Periode von 0,5 [Einheiten] extrahieren möchten, benötigen Sie eine Probe alle 0,25 [Einheiten].

Beachten Sie, dass dies eine mathematische Begrenzung ist, nicht einer von R.