Hochpassfilterung von Sensordaten
Verfasst: Mi Mär 29, 2017 5:13 pm
Hallo zusammen,
Im Zuge meiner Masterarbeit analysiere ich Daten von Drohnenflügen (z.B. Beschleunigungswerte in X, Y oder Z-Richtung einer Inertial Measurement Unit). Die Daten sollen zunächst hochpass-gefiltert werden.
Dafür werden die Daten als erstes mit einem Gaußfilter (smth.gaussian, alpha = ?, window = ?, tail =?) tiefpass-gefiltert. Die dabei erzeugten Werte sollen dann von den Ausgangswerten subtrahiert werden, was als Ergebnis hochpassgefilterte Werte liefert. Im Anschluss soll ein "moving windows" mit einer bestimmten Breite (in Sekunden) über die ganze Datenreihe geschoben werden. Dabei soll mit einer Schleife jede Zeile i bis i+1 usw. des gesamten Datensatzes erfasst und für jedes Zeitfenster Mittelwert und Standardabweichung berechnet werden. Eine genaue Beschreibung findet sich im Anhang.
Ich wäre sehr dankbar, wenn ich Ratschläge bekommen könnte, wie ich dieses Problem in einem R-Skript umsetzen könnte!
-welche Werte für alpha, window, tail
-implementierung eines moving windows
Beste Grüße
Benedikt
Im Zuge meiner Masterarbeit analysiere ich Daten von Drohnenflügen (z.B. Beschleunigungswerte in X, Y oder Z-Richtung einer Inertial Measurement Unit). Die Daten sollen zunächst hochpass-gefiltert werden.
Dafür werden die Daten als erstes mit einem Gaußfilter (smth.gaussian, alpha = ?, window = ?, tail =?) tiefpass-gefiltert. Die dabei erzeugten Werte sollen dann von den Ausgangswerten subtrahiert werden, was als Ergebnis hochpassgefilterte Werte liefert. Im Anschluss soll ein "moving windows" mit einer bestimmten Breite (in Sekunden) über die ganze Datenreihe geschoben werden. Dabei soll mit einer Schleife jede Zeile i bis i+1 usw. des gesamten Datensatzes erfasst und für jedes Zeitfenster Mittelwert und Standardabweichung berechnet werden. Eine genaue Beschreibung findet sich im Anhang.
Ich wäre sehr dankbar, wenn ich Ratschläge bekommen könnte, wie ich dieses Problem in einem R-Skript umsetzen könnte!
-welche Werte für alpha, window, tail
-implementierung eines moving windows
Beste Grüße
Benedikt