数字信号处理第五次实验报告.docVIP

电子科技大学生命科学与技术学院 标 准 实 验 报 告 （实验）课程名称 数字信号处理 2017-2018-第2学期 电子科技大学教务处制表 电 子 科 技 大 学 实 验 报 告 学生姓名 学 号 20160102008 指导教师 实验时间 2018年 4月 26 日 一、实验室名称：清水河校区，基础实验大楼 506 机房 二、实验名称： 三、实验学时：2学时 四、实验原理： 五、实验目的： 六、实验内容： clear clc; data=load(EEGdata.txt); EEG=fft(data); data1=data+0.2*rand(length(data),1)-0.1; EEG1=fft(data1); x=1:fix(length(data)/2); figure (1); subplot(2,1,1);plot(1:length(data),data,b,1:length(data),data1,r); title(subplot 1:时域信号图); ylabel(振幅),xlabel(样本点序号); subplot(2,1,2);plot(x/fix(length(data)/2),abs(EEG(x)),b,x/fix(length(data)/2),abs(EEG1(x)),r); legend(原信号,噪声信号) title(subplot 信号幅度谱); ylabel(振幅),xlabel(归一化频率); %低通滤波器，截止频率为归一化频率wp=0.092，wc=0.1； Wp=tan(0.092*pi/2);Ws=tan(0.178*pi/2); Rp=1;Rs=40; [N,Wn]=buttord(Wp,Ws,Rp,Rs,s); N [b,a]=butter(N,Wn,s); [b1,a1]=bilinear(b,a,0.5); [h1,omega]=freqz(b1,a1); db=20*log10(abs(h1)); figure (2); subplot(4,3,1);plot((1:512)/512,abs(h1)); title(subplot 巴特沃兹滤波器幅度谱); ylabel(振幅),xlabel(归一化频率); subplot(4,3,2);plot((1:512)/512,db); title(subplot 波特图); ylabel(振幅),xlabel(归一化频率); [N,Wn]=cheb1ord(Wp,Ws,Rp,Rs,s); N [b,a]=cheby1(N,Rp,Wn,s); [b2,a2]=bilinear(b,a,0.5); [h2,omega]=freqz(b2,a2); db=20*log10(abs(h2)); subplot(4,3,4);plot((1:512)/512,abs(h2)); title(subplot 切比雪夫1型滤波器幅度谱); ylabel(振幅),xlabel(归一化频率); subplot(4,3,5);plot((1:512)/512,db); title(subplot 波特图); ylabel(振幅),xlabel(归一化频率); [N,Wn]=cheb2ord(Wp,Ws,Rp,Rs,s); N [b,a]=cheby2(N,Rs,Wn,s); [b3,a3]=bilinear(b,a,0.5); [h3,omega]=freqz(b3,a3); db=20*log10(abs(h3)); subplot(4,3,7);plot((1:512)/512,abs(h3)); subplot(4,3,8);plot((1:512)/512,db); title(subplot 波特图); ylabel(振幅),xlabel(归一化频率); [N,Wn]=ellipord(Wp,Ws,Rp,Rs,s); N [b,a]=ellip(N,Rp,Rs,Wn,s); [b4,a4]=bilinear(b,a,0.5); [h4,omega]=freqz(b4,a4); subplot(4,3,10);plot((1:512)/512,abs(h4)); title(subplot 椭圆滤波器幅度谱); ylabel(振幅),xlabel(归一化频率); subplot(4,3,11);plot((1:512)/512,db); title(subplot 波特图); ylabel(振幅),xlabel(归一化频率); data11=filtfilt(b1,a1,data1); subplot(4,3,3) plot(1:length(data1),data11,b,