180 likes | 347 Views
EEG analysis during hypnagogium. Petr Svoboda Laboratory of System Reliability Faculty of Transportation Czech Technical University e-mail: svobodap@spel.cz. Presented methods. Traditional methods Fourier transform Parametrical methods Autoregressive estimator
E N D
EEG analysis during hypnagogium Petr Svoboda Laboratory of System Reliability Faculty of Transportation Czech Technical University e-mail: svobodap@spel.cz
Presented methods • Traditional methods Fourier transform • Parametrical methods Autoregressive estimator • Nonlinear methods - Chaos theory Delay-time embedding, Correlation dimension, Takens estimator, State Space dimension, Lyapunov exponents
EEG activity • electric potential of brain‘s neural activity • registered on the skelet • four basic frequencies
Traditional methods Estimate of a periodogram using the Fourier transform Potencial problems • Signal‘s stationarity • frequency resolution • leakage of frequency spectra • quality of the spectral estimate • phase of the signal is lost
Parametric model Approximation of an EEG signal by adequate parametric model Autoregressive (AR) model: Approximation of an EEG signal bylinear time invariant filter with transfer function H(z)=1/A(z) Whitening of signal by AR filter:
Analyzed signal Pole placement Autocovariance function Spectral estimate
Comparison of traditional and parametric methods Traditional methods: + low noise sensitivity - frequency resolution Parametric methods: + frequency resolution + parametric description of analyzed signal - estimate of AR model order - high noise sensitivity
Microsleep classification Traditional methods + Alpha, deltha and theta activity of spectral estimate Parametrical methods + Alpha, deltha and theta activity of spectral estimate - estimate of AR model order - placement of poles in a complex plane
Classification by spectral estimate • Classification into 2 states • RELAXATION • DROWSINESS • Classification based on neural network (back propagation) • Classical methods: accuracy of about 87% • Parametrical methods: accuracy of about 90%
Relaxation Drowsiness
Chaos theory • analysis of dynamic deterministic systems • high sensitivity on initial conditions • known dynamics and phase of the system • detecting nonlinearity by surrogate data testing • delay-time embedding • state-space dimension estimate • estimate of delay time • estimate of fractal dimension D2 • Takens estimator for D2 dimension • largest Lyapunov exponents
Delay-time embedding Si=[x(i),x(i+L),… x(i+(m-1)L)] L… time delay Si… state-space vector m… state dimension x… analyzed signal
Selection of Delay Time L Time delay should be set so, x(i),x(i+L),… are independent • autocorrelation method • method of Mutual Information (MI)
Fractal dimension & Takens estimator Fractal dimension is a measure of complexity of the analyzed signal D2 = log C(r) / log r Where C(r) is correlation integral D2 computed by maximum likelihood estimator is known as Takens estimator
Microsleep classification Chaos theory + matematical description of state-space trajectory reconstruction nonlinearity detection - correlation dimension D2 estimate + Takens estimator + largest Lyapunov exponents
Largest Lyapunov Exponents Sensitive dependence on initial conditions