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

fd=1/T;

fn=fd/2;

f=-fd/2:fd/(length(y)-1):fd/2;

subplot(5,1,2);

stem(f,fftshift(A));

title('spectrum');

AKF=xcorr(y);

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

subplot(5,1,3);

plot(tau,AKF);

title('AKF');

[P,f]=pmtm(y,[],[],fd);

subplot(5,1,4);

plot(f,P);

title('periodogram');

subplot(5,1,5);

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

title('specgram');

%filt param

n=170;

a=1;

w=hamming(length(y)+1);

fn=fd/2;

%alpha Fp1

figure(2);

fx=[8 13];

b=fir1(n,fx/fn);

[h,wn]=freqz(b,a);

subplot(6,1,1);