27.11.2014 Views

Quantitative analysis of EEG signals: Time-frequency methods and ...

Quantitative analysis of EEG signals: Time-frequency methods and ...

Quantitative analysis of EEG signals: Time-frequency methods and ...

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

(ml)<br />

j ( n ) =<br />

2 m X<br />

d=1<br />

<strong>and</strong> the corresponding coecients are easily computed as<br />

where 1 d 2 m .<br />

D (ml)<br />

j ( d ) =<br />

D (ml)<br />

j ( d ) (ml)<br />

jd<br />

( n ) : (40)<br />

l+2 m ;1<br />

X<br />

k=l<br />

M (m)<br />

dk<br />

C j ( k ) (41)<br />

The trigonometric wavelet packets (ml)<br />

jd<br />

(n) have zero mean, oscillate on the interval<br />

2 ;j l n 2 ;j (l +2 m ) <strong>and</strong> decay with exponential ratio. Moreover, their wave-forms<br />

resemble modulated sines or cosines. In fact, it can be demonstrated that each Fourier<br />

(ml)<br />

transform ^<br />

jd<br />

(!) is centered at the fundamental <strong>frequency</strong> ! mh , when d = 2h or<br />

(ml)<br />

d =2h +1. Moreover, ^ jd (!) =0on the other fundamental frequencies.<br />

In other words, the coecients fD (ml)<br />

j<br />

spectrum<br />

for the local signal r (ml)<br />

j<br />

<strong>of</strong> coecients fC j (k)D (ml)<br />

j<br />

(d)g can be considered as the discrete Fourier<br />

(n). Summing up, we can resume in the double set<br />

(d)g the time-scale-<strong>frequency</strong> information <strong>of</strong> the local signal<br />

r (ml)<br />

j (n).<br />

Finally, to analyze the complete function r j (n), that is, the details at level j, we<br />

choose some partition in local components r (m il i )<br />

j (n), according the structure <strong>of</strong> the<br />

signal,<br />

r j ( n ) = X m i<br />

r (m il i )<br />

j ( n ) (42)<br />

where the sequence <strong>of</strong> index l i veries l i+1 = l i +2 m i<br />

. Then, we implement the above<br />

refereed time-scale-<strong>frequency</strong> technique for each local signal.<br />

4.3 Short review <strong>of</strong> wavelets applied to the study <strong>of</strong> <strong>EEG</strong> <strong>signals</strong><br />

Several works applied the Wavelet Tranform to the study <strong>of</strong> <strong>EEG</strong>s <strong>and</strong> ERPs (see a<br />

review in Unser <strong>and</strong> Aldroubi, 1996 or in Samar et al., 1995). One rst line <strong>of</strong> applications<br />

is for pattern recognition in the <strong>EEG</strong>. This is achieved by correlating dierent<br />

transients <strong>of</strong> the <strong>EEG</strong> with wavelet coecients <strong>of</strong> dierent scales. Schi et al. (1994a)<br />

used amultiresolution decomposition implemented with B-Splines mother functions for<br />

extracting features <strong>of</strong> <strong>EEG</strong> seizure recordings. They showed a better performance <strong>of</strong><br />

wavelets in comparison with Gabor Transform, <strong>and</strong> a similar resolution <strong>of</strong> the multiresolution<br />

decomposition compared with the continuous Wavelet Transform but with a<br />

44

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!