Курсовая работа: Обработка электроэнцефалограмм в частотной области

title('filter');

y1=filtfilt(b,a,y);

subplot(6,1,2);

plot(t,y1);

title('alpha eeg_ Fp2_1.txt');

AKF=xcorr(y1);

tau=-t(length(t)):2*t(length(t))/(length(AKF)-1):t(length(t));

subplot(6,1,3);

plot(tau,AKF);

title('alpha AKF');

A1=abs(fft(y1));

f1=-fd/2:fd/(length(y1)-1):fd/2;

subplot(6,1,4);

stem(f1,fftshift(A1));

title('alpha spectrum');

[P,f_1]=pmtm(y1,[],[],fd);

subplot(6,1,5);

plot(f_1,P);

title('periodogram');

subplot(6,1,6);

specgram(y1,[],fd,100);

title('specgram');

Приложение В

РЕЗУЛЬТАТЫ АНАЛИЗА ЭЭГ В ПРОГРАММНОЙ СРЕДЕ M AT L AB

Рисунок B.1 – Входной ЭЭГ сигнал ('eeg_Fp1.txt'), его АЧХ, АКФ, периодограмма и спектрограмма

Рисунок B.2 – Фильтр альфа-ритма сигнала 'eeg_Fp1.txt', альфа-ритм, его АКФ, АЧХ, периодограмма и спектрограмма

Рисунок B.3 – Входной ЭЭГ сигнал ('eeg_Fp2.txt'), его АЧХ, АКФ, периодограмма и спектрограмма