数字信号处理大作业.doc

数据说明 本实验共分为4组子实验，分别将四组数据命名为*2.1*2.4，刺激频率7Hz30Hz，采样率为1000Hz。 刺激序列： 每个子实验按照上述序列闪烁6次，即每个频率刺激8s6=48s。 刺激频率：刺激序列A/(HZ)B/(HZ)C/(HZ)D/(HZ)E/(HZ)F/(HZ)子实验1789101112子实验2131415161718子实验3192021222324子实验4252627282930 数据格式为.mat格式，是一个68(channel)n(sample)的数组形式。读入后有两个数组，data是数据，event是标记点。题目：1. 提取12Hz的数据，共6个trials，采样率为1000Hz，每次任务持续8秒，休息4秒，按照event提取12Hz数据，并用detrend函数去基线。最后得到channelssamplestrials的数据格式。2. 绘制第62导，第3个trial的前5000点脑电数据。3. 计算第62导，第5个trial数据的功率谱，并使用方法一，对所有trial进行叠加平均后再求功率谱，绘制功率谱曲线进行比较。（提示：对于随机信号，可以直接使用fft做傅立叶变换，再求幅频响应以表示功率谱；注意：绘制功率谱曲线注意频率轴的定标。）4. 设置带通滤波器，要求带通范围为820Hz，通带衰减和阻带衰减分别为3，20dB，左右各设置3Hz宽度的边缘。参看matlab中cheb1ord例子使用cheb1方法设计滤波器，并绘制频率响应。5. 采用方法二，对第62导，第5个trial数据进行滤波，求功率谱，绘制原始数据、平均后估计数据以及滤波后数据得到的功率谱曲线进行比较分析。1、提取12Hz的数据，共6个trials，采样率为1000Hz，每次任务持续8秒，休息4秒，按照event提取12Hz数据，并用detrend函数去基线。最后得到channelssamplestrials的数据格式。load matlab.mat;fs=1000;t=8;k=1;for i=6:6:length(event); temp=detrend(data(:,event(i):event(i)+fs*8-1); datapick(:,:,k)=temp; k=k+1;end;2、 绘制第62导，第3个trial的前5000点脑电数据。figureplot(datapick(62,1:5000,3);title(第62导第3个tril的前5000点脑电数据);3、 计算第62导，第5个trial数据的功率谱，并使用方法一，对所有trial进行叠加平均后再求功率谱，绘制功率谱曲线进行比较。（提示：对于随机信号，可以直接使用fft做傅立叶变换，再求幅频响应以表示功率谱；注意：绘制功率谱曲线注意频率轴的定标。）figureY=fft(datapick(62,:,5),fs*t)/(fs*t);f=fs/2*linspace(0,1,(fs*t)/2+1);plot(f,2*abs(Y(1:(fs*t)/2+1);title(channel62,tril=5的单边幅度谱);axis(0,100,0,4);xlabel(frequence(hz);ylabel(|Y|);figureX=mean(datapick(62,:,:),3);Y1=fft(X,fs*t)/(fs*t);f1=fs/2*linspace(0,1,(fs*t)/2+1);plot(f1,2*abs(Y1(1:(fs*t)/2+1);title(channel62,tril平均后的单边幅度谱);axis(0,100,0,4);xlabel(frequence(hz);ylabel(|Y1|);4、 设置带通滤波器，要求带通范围为820Hz，通带衰减和阻带衰减分别为3，20dB，左右各设置3Hz宽度的边缘。参看matlab中cheb1ord例子使用cheb1方法设计滤波器，并绘制频率响应。figurefa=8,fb=20;fl=5, fh=23;rp=3, rs=20;wp1=2*pi*fa/fs;wp3=2*pi*fb/fs;wl=2*pi*fl/fs;wh=2*pi*fh/fs;wp=wp1 wp3;ws=wl wh;n,wn=cheb1ord(ws/pi,wp/pi,rp,rs);bz1,az1=cheby1(n,rp,wp/pi)h,w=freqz(bz1,az1,256,fs);h=20*log10(abs(h);plot(w,h);grid on;axis(0,50,-50,1);title(切比雪夫滤波器);5、 采用方法二，对第62导，第5个trial数据进行滤波，求功率谱，绘制原始数据、平均后估计数据以及滤波后数据得到的功率谱曲线进行比较分析。figuredf=filter(bz1,az1,datapick(62,:,5);YF=fft(df,fs*t)/(fs*t);ff=fs/2*linspace(0,1,(fs*t)/2+1);plo