数据说明。
本实验共分为4组子实验，分别将四组数据命名为***2.1~**2.4，刺激频率7hz～30hz，采样率为1000hz。刺激序列：
arbrcrdrerfr
8s4s8s4s8s4s8s4s8s4s8s30s
每个子实验按照上述序列闪烁6次，即每个频率刺激8s×6=48s。刺激频率：
刺激序列。a/(hz)b/(hz)c/(hz)d/(hz)e/(hz)f/(hz)
子实验1子实验2子实验3子实验4
数据格式为。mat格式，是一个68(channel)×n(sample)的数组形式。读入后有两个数组，data是数据，event是标记点。
题目：1.提取12hz的数据，共6个trials，采样率为1000hz，每次任务持续8秒，休息4秒，按照event提取12hz数据，并用detrend函数去基线。
最后得到channels×samples×trials的数据格式。
2.绘制第62导，第3个trial的前5000点脑电数据。
3.计算第62导，第5个trial数据的功率谱，并使用方法一，对所有trial进行叠加平均后再求功率谱，绘制功率谱曲线进行比较。（提示：
对于随机信号，可以直接使用fft做傅立叶变换，再求幅频响应以表示功率谱；注意：绘制功率谱曲线注意频率轴的定标。）
4.设置带通滤波器，要求带通范围为8～20hz，通带衰减和阻带衰减分别为3，20db，左右各设置3hz宽度的边缘。参看matlab中cheb1ord例子使用cheb1方法设计滤波器，并绘制频率响应。
5.采用方法二，对第62导，第5个trial数据进行滤波，求功率谱，绘制原始数据、平均后估计数据以及滤波后数据得到的功率谱曲线进行比较分析。
1、提取12hz的数据，共6个trials，采样率为1000hz，每次任务持续8秒，休息4秒，按照event提取12hz数据，并用detrend函数去基线。最后得到channels×samples×trials的数据格式。
load 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点脑电数据。figure
plot(datapick(62,1:5000,3));
title('第62导第3个tril的前5000点脑电数据');
3、计算第62导，第5个trial数据的功率谱，并使用方法一，对所有trial进。
行叠加平均后再求功率谱，绘制功率谱曲线进行比较。（提示：对于随机信号，可以直接使用fft做傅立叶变换，再求幅频响应以表示功率谱；注意：
绘制功率谱曲线注意频率轴的定标。）figure
y=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|')
figure
x=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、设置带通滤波器，要求带通范围为8～20hz，通带衰减和阻带衰减分别为。
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数据进行滤波，求功率谱，绘制原。
始数据、平均后估计数据以及滤波后数据得到的功率谱曲线进行比较分析。figure
df=filter(bz1,az1,datapick(62,:,5));yf=fft(df,fs*t)/(fs*t);
ff=fs/2*linspace(0,1,(fs*t)/2+1);plot(ff,2*abs(yf(1:(fs*t)/2+1)))
title('channel62,tril5滤波后的单边幅度谱');axis([0,100,0,4]);xlabel('frequence(hz)')ylabel('|yf|')
数字信号处理大作业
题目一 利用matlab设计模拟带通巴特沃斯滤波器，要求通带下 上截止频率分别为fhpl 4khz,fhpu 6khz，阻带下 上截止频率分别为fhsl 2khz，fhsu 9khz，通带最大衰减rp 1db，阻带最小衰减rs 20db。问题分析 1 该题目中模拟带通滤波器指标 通带下截止频率 wp...
数字信号处理大作业
设计步骤 这次选择设计的是巴特沃斯数字高通滤波器，对输入信号xn cos 8 pi t cos 20 pi t 最高频率f 20 pi 2 pi 10hz，所以采样频率fs应大于等于20hz，在这里我选择的采样频率fs 80hz。设计思想 从归一化模拟低通原型出发，先在模拟域内经频率变换成为所需类型...
数字信号处理大作业
数字信号处理。大作业。第二次大作业。n 0 10 x 0.9 n h 400 400 w pi 100 h xjw1 x exp j pi 100 n h subplot 2,2,1 plot w pi,abs xjw1 grid on xlabel w1 pi ylabel x jw1 title...