22.01.2015 Aufrufe

INSTITUT FÜR INFORMATIK Entwurf und Implementierung - Ludwig ...

INSTITUT FÜR INFORMATIK Entwurf und Implementierung - Ludwig ...

INSTITUT FÜR INFORMATIK Entwurf und Implementierung - Ludwig ...

MEHR ANZEIGEN
WENIGER ANZEIGEN

Sie wollen auch ein ePaper? Erhöhen Sie die Reichweite Ihrer Titel.

YUMPU macht aus Druck-PDFs automatisch weboptimierte ePaper, die Google liebt.

A.2 Dynamical Similarity Index<br />

A.2 Dynamical Similarity Index<br />

Beim Dynamical Similarity Index wird die Dynamik eines gleitenden Datenfensters<br />

S t mit der eines Referenzfensters S ref verglichen [20, 37], wie in Abschnitt 2.3.2<br />

angeführt. Zum Zweck der Vorhersage epileptischer Anfälle wird dabei aus den EEG-<br />

Daten ein Referenzfenster gewählt, das fern von Anfällen liegt, <strong>und</strong> möglichst frei von<br />

Anfallsaktivität ist.<br />

Um die den Fenstern zugr<strong>und</strong>eliegende Dynamik zu rekonstruieren, werden die<br />

Zeitpunkte der positiven Nulldurchgänge T n der Daten ermittelt, <strong>und</strong> die Zeitintervalle<br />

zwischen den Nulldurchgängen betrachtet, was zu einer neuen Zeitreihe<br />

I n = T n+1 −T n führt. Aus dieser werden über eine Einbettung der Dimension m = 15<br />

die Einbettungsvektoren A n = (I n , I n−1 , ..., I n−m−1 ) bestimmt [38]. Zur Reduktion<br />

des Rauschens wird eine Singulärwertzerlegung der Trajektorienmatrizen A(S ref )<br />

bzw. A(S t ) angewendet. Diese werden anschlieÿend auf die vier gröÿten Hauptachsen<br />

des Referenzfensters projiziert, wodurch man X(S t ) bzw. X(S ref ) erhält.<br />

Da das Referenzfenster normalerweise mit einer Länge von fünf Minuten gewählt<br />

wird, während für das gleitende Fenster eine Länge von 25 Sek<strong>und</strong>en gewählt wurde,<br />

wird eine zufällig eine Teilmenge Y(S ref ) von X(S ref ) mit N ref ≈ N t Punkten für<br />

den Vergleich der Fenster bestimmt. Als ein Maÿ der Ähnlichkeit beider Fenster wird<br />

das Kreuzkorrelationsintegral<br />

C(S ref , S t ) =<br />

1<br />

N ref N t<br />

∑ ∑<br />

i=1,N ref<br />

j=1,N t<br />

Θ (r − ||Y i (S ref ) − X j (S t )||)<br />

(A.7)<br />

berechnet [39], wobei || · || die Euklidische Norm ist. Die Heaviside-Funktion Θ liefert<br />

hier den Wert 1 zurück, falls die Punkte Y i (S ref ) <strong>und</strong> X j (S t ) in einem Abstand<br />

geringerer als r liegen, ansonsten 0. Standardmäÿig wird der Abstand so gewählt,<br />

dass sich im Mittel von jedem beliebigen Punkt des Referenzfensters mindestens 30%<br />

aller Punkte des Referenzfensters in einem Abstand von höchstens r benden. Zur<br />

Normierung auf Werte im Bereich [0, 1] ist der Dynamical Similarity Index deniert<br />

als<br />

γ(S t ) =<br />

C(S ref , S t )<br />

√<br />

C(Sref , S ref )C(S t , S r ) .<br />

(A.8) v

Hurra! Ihre Datei wurde hochgeladen und ist bereit für die Veröffentlichung.

Erfolgreich gespeichert!

Leider ist etwas schief gelaufen!