数字信号处理—基于计算机的方法双线性变换法的设计IIR数字滤波器.docVIP

电子科技大学生命科学与技术学院 标 准 实 验 报 告 （实验）课程名称 数字信号处理 2013-2014-第2学期 电子科技大学教务处制表 电 子 科 技 大 学 实 验 报 告 学生姓名 学 号 指导教师 实验时间 2014年 5月 29日 一、实验室名称：清水河校区，实验大楼 A508 机房 二、实验名称： 三、实验学时：2学时 四、实验原理： 五、实验目的：（详细填写） 借助于模拟滤波器，利用双线性变换法设计IIR数字滤波器： 1）掌握其基本思想和步骤； 2）学会使用有关的MATLAB函数。 六、实验内容：（详细填写） 1)加载信号(EEGdata.txt)； 2)分析信号的幅度谱,确定低通IIR数字滤波器的指标； 3)利用双线性变换法设计Butterworth、切比雪夫I型、II型、椭圆数字低通滤波器； 4)对信号加随机噪声； 5)用设计的滤波器对含噪声的信号进行滤波。 要求： 1）给出程序，画出滤波器幅度谱及其增益图(dB)；分析设计的滤波器是否达到指标。 2）分析不同方式所需的滤波器阶数； 3）利用设计的四种滤波器对加载的信号分别进行滤波，对比分析滤波结果(时间域及频谱)。 七、实验器材（设备、元器件）： 八、实验步骤： 九、实验数据及结果分析：（详细填写）（包括程序、图、结果等） Butterworth滤波器 程序： %输入原始信号并观察特点 data=load(EEGdata.txt); n0=0:length(data)-1; x0=n0/350*2*pi; f=fft(data);amp=abs(f); figure(1); subplot(2,1,1);plot(x0,data);title(原信号); subplot(2,1,2);plot(x0,amp,r);title(原信号的幅度谱)%构造原信号的幅度谱 %巴特沃斯滤波器的设计 wp=0.5095;ws=1.025;% 模拟指标 [n,wn]=buttord(wp,ws,0.5,50,s); [b,a]=butter(n,wn,s); [b,a]=bilinear(b,a,0.5); [h,w0]=freqz(b,a,512); figure(2); subplot(2,1,1);plot(w0/pi,abs(h));title(巴特沃斯滤波器幅度谱); subplot(2,1,2);plot(w0/pi,20*log10(abs(h)));title(巴特沃斯滤波器增益谱) %叠加噪声后的信号 noise=0.2*rand(350,1)-0.1; y=data+noise; amp1=abs(fft(y)); %滤波比较效果 lvbo=filtfilt(b,a,y); amp2=abs(fft(lvbo)); figure(3); subplot(2,1,1);plot(x0,y,r,x0,lvbo,b,x0,data,k-);title(滤波前后信号的对比); subplot(2,1,2);plot(x0,amp1,r,x0,amp2,b,x0,amp,k-);title(滤波前后的信号幅度谱对比); 实验结果 切比雪夫I型： 程序 %% %输入原始信号并观察特点 data=load(EEGdata.txt); n0=0:length(data)-1; x0=n0/350*2*pi; f=fft(data);amp=abs(f); figure(1); subplot(2,1,1);plot(x0,data);title(原信号); subplot(2,1,2);plot(x0,amp,r);title(原信号的幅度谱)%构造原信号的幅度谱 %% %切比雪夫一型滤波器设计 wp=0.5095;ws=1.025;% 模拟指标 [n,wn]=cheb1ord(wp,ws,0.5,50,s); [b,a]=cheby1(n,0.5,wn,s); [b,a]=bilinear(b,a,0.5); [h,w0]=freqz(b,a,512); figure(2); subplot(2,1,1);plot(w0/pi,abs(h));title(切比雪夫一型滤波器幅度谱); subplot(2,1,2);plot(w0/pi,20*log10(abs(h)));title(切比雪夫一型滤波器增益谱) %% %叠加噪声后的信号 noise=0.2*rand(350,1)-0.1; y=data+noise; amp1=abs(fft(y)); %% %滤波比较效果 lvbo=filtfilt(b,a,y); amp2=abs(fft(lvbo)); figure(3