数字信号处理作业
1.两个导联C3,C4位置的脑电信号（已预处理），实验采样频率为
250Hz，每次实验采集8秒数据，总共做了36次实验。

依次求出C3,C4位置第1秒~第8秒数据的功率谱。

clc
clear
load('C:\Users\刘冰\Desktop\数字信号处理\matlab\C3C4.mat')
a(1:8,1:512)=zeros();
for j=1:8
for 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; %求其功率谱end
a(j,:)=p(j,:)./36;%求平均值
end
p(1:8,1:512)=zeros();
for j=1:8
for 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;
end
p(j,:)=q(j,:)./36;
end
for 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
0.20.40.60.81 1.2 1.4 1.6 1.82
0100
200
300
400
500
600
700
第1秒 C3、C4脑电功率谱对照
0.20.40.60.81 1.2 1.4 1.6 1.82
0100
200
300
400
500
600
700
0.20.40.60.81 1.2 1.4 1.6 1.82
0100200300400500600700800第3秒 C3、C4脑电功率谱对照
0.20.40.60.81 1.2 1.4 1.6 1.82
0100
200
300
400
500
600
0.20.40.60.81 1.2 1.4 1.6 1.82
050100150200250300350400第5秒 C3、C4脑电功率谱对照
0.20.40.60.81 1.2 1.4 1.6 1.82
050100150200250300350400
0.2
0.4
0.6
0.8
1
1.2
1.4
1.6
1.8
2050100150200250300350400450500第7秒 C3、C4脑电功率谱对
照
0.20.40.60.81 1.2 1.4 1.6 1.82
050100150200250300350400。