脑电信号功率谱.docx

数字信号处理作业两个导联C3,C4位置的脑电信号（已预处理），实验采样频率为250Hz，每次实验采集8秒数据，总共做了36次实验。依次求出C3,C4位置第1秒~第8秒数据的功率谱。clcclearload(C:\Users\刘冰\Desktop\数字信号处理\matlab\C3C4.mat)a(1:8,1:512)=zeros();for j=1:8for k=0:35; z=fft(Left_C3(((j-1)*250+1+2000*k):(2000*k+j*250)),512); %截取特定的一段数据进行傅里叶变换a(j,:)=p(j,:)+z.*conj(z)/512; %求其功率谱enda(j,:)=p(j,:)./36;%求平均值endp(1:8,1:512)=zeros();for j=1:8for k=0:35; z=fft(Left_C4((j-1)*250+1+2000*k:2000*k+j*250),512); 、%截取特定的一段数据进行傅里叶变换p(j,:)=q(j,:)+z.*conj(z)/512;endp(j,:)=q(j,:)./36;endfor i=1:8 w=0:2*pi/255:2*pi; figure plot(w/pi,p(i,1:256),b,w/pi,q(i,1:256),r)%在一幅图里面显示C3C4功率谱，因为其结果是对称的，所以只取前一半结果 legend(C3,C4);%线段标题 title([第,num2str(i),秒C3、C4脑电功率谱对照])end