数字滤波器设计及应用综合实验样本[参照].pdfVIP

资料内容仅供您学习参考，如有不当或者侵权，请联系改正或者删除。 数字信号处理实验 实验四、 数字滤波器设计及 应用综合实验 学院 : 信息工程学院 班级 : 电子 101 班 姓名 : 学号 : 资料内容仅供您学习参考，如有不当或者侵权，请联系改正或者删除。 一、 实验目的 1．熟悉 IIR 数字滤波器的设计原理及方法。 2 ．熟悉 FIR 数字滤波器的设计原理及方法。 3. 掌握利用 Matlab 实现数字滤波器的方法 4. 掌握利用数字滤波器进行信号处理的方法。 5. 了解基于 Simulink 的动态仿真实现信号滤波的基本方法。 二、 实验内容及要求 实验内容 : 综合运用数字滤波器设计的相关知识 , 根据给定设计方法要求 , 用脉冲响 应不变法和双线性变换法设计 IIR 数字滤波器 ; 利用窗函数设计法设计 FIR数字 滤波器。根据实际信号的频谱特性 , 分析、 确定滤波器设计技术指标 , 实现对 信号的滤波。 1．IIR 数字滤波器设计 ( 1) 用脉冲响应不变法设计巴特沃斯数字滤波器。 该实验所需 M文件如下 : ①、 butterworth 低通滤波器原型设计函数 : function [b,a]=afd_butt(Wp,Ws,Rp,As) N=ceil((log10((10^(Rp/10)-1)/(10^(As/10)-1)))/(2*log10(Wp/Ws))) fprintf(\n Butterworth Filter Order=%2.0f\n,N) OmegaC=Wp/((10^(Rp/10)-1)^(1/(2*N))) [b,a]=u_buttap(N,OmegaC) ②、 非归一化 Butterworth 模拟低通滤波器设计函数 : function [b,a]=u_buttap(N,Omegac); [z,p,k]=buttap(N); p=p*Omegac; k=k*Omegac^N; B=real(poly(z)); 资料内容仅供您学习参考，如有不当或者侵权，请联系改正或者删除。 b=k*B; a=real(poly(p)); ③、 利用脉冲响应不变法从模拟到数字滤波器变换函数 : function [b,a]=imp_invr(c,d,T) [R,p,k]=residue(c,d); p=exp(p*T); [b,a]=residuez(R,p,k); b=real(b); a=real(a); ④、 频率响应函数 freqz 的修正 : function [db,mag,pha,w]=freqz_m(b,a); [H,w]=freqz(b,a,1000,whole); H=(H(1:501)); w=(w(1:501)); mag=abs(H); db=20*log10((mag+eps)/max(mag)); pha=angle(H); 本实验程序如下 : wp=0.2*pi;ws=0.3*pi;Rp=1;As=15;T=1; OmegaP=wp/T;OmegaS=ws/T; [cs,ds]=afd_butt(OmegaP,OmegaS,Rp,As); [b,a]=imp_invr(cs,ds,T) [db,mag,pha,w]=freqz_m(b,a); subplot(2,1,1);plot(w/pi,mag); title(digital filter Magnitude Res