数字高通带通带阻滤波器设计绪论.pptVIP

数字高通、带通、带阻滤波器的设计——第三次翻转课堂 目录 数字高通、带通、带阻滤波器设计简介 数字带通滤波器设计实例 直接使用Matlab自带函数求解 数字高通、带通、带阻滤波器设计简介 1 对于IIR数字滤波器的设计具体步骤如下: (1)按照一定的规则将给出的数字滤波器的技术指标转换为模拟滤波器的技术指标。 (2)根据转换后的技术指标设计模拟低通滤波器G(p)(低通滤波器的归一化原型系统函数)。 (3)利用频率变换公式将G(p)转换成H(s)(模拟高通滤波器的传递函数)。 (4)利用双线性变换法将H(s)转换成H(z)(数字高通滤波器的传递函数)。 1 设计 步骤 数字高通、带通、带阻滤波器设计简介 2 技术 指标 数字技术指标： 通带边界频率：wp(低通) wph (高通) wpt(带通) wpz (带阻) 阻带截止频率：ws(低通) wsh (高通) wst(低通) wsz (高通) 模拟技术指标： 通带边界频率：Ωp (低通) Ωph (高通) Ωpt (低通) Ωpz (高通) 阻带截止频率：Ωs (低通) Ωsh (高通) Ωst (低通) Ωsz (高通) 通带内允许的最大衰减：αp 阻带内允许的最大衰减：αs 数字高通、带通、带阻滤波器设计简介 3 技术 指标 转换 脉冲响应不变法： Ωp =wp/T Ωs =ws/T 双线性变换法：(预畸变校正) Ωp =2tan(wp/2)/T Ωs =2tan(ws/2)/T 高通到低通： Ωp=Ωc=1 Ωs =Ωph/Ωsh 带通到低通： Ωp=Ωc=1 Ωs = 带阻到低通： Ωp=Ωc=1 Ωs = 注： 数字高通、带通、带阻滤波器设计简介 4 阶数N计算： λsp=Ωs/Ωp ksp = N=lgksp /lgλsp 由N查表得到归一化原型系统函数G(p) 计算 过程 数字高通、带通、带阻滤波器设计简介 5 频率 变换 公式 低通到高通： 低通到带通： 低通到带阻： 数字带通滤波器设计实例 1 确定技术指标 数字带通技术指标: 通带范围： 阻带范围： 通带最大衰减：αp=3db 阻带最大衰减：αs=40db 数字低通技术指标: 通带边界频率：wp=1 阻带截止频率：ws=2 通带最大衰减：αp=3db 阻带最大衰减：αs=40db 数字带通滤波器设计实例 2 Matlab求解： 模拟低通滤波器程序： wp=1;ws=2;rp=3;rs=40;T=1; OmegaP=(2/T)*tan(wp/2);OmegaS=(2/T)*tan(ws/2); [cs,ds]=afd_butt(OmegaP,OmegaS,rp,rs); [h,w]=freqs(cs,ds); subplot(3,1,1);plot(w/pi,abs(h)); title(low anolog filter Magnitude Response) 数字带通滤波器设计实例 2 Matlab求解： 模拟带通滤波器程序： wp1=2*pi*[0.2,0.3];ws1=2*pi*[0.1,0.4]; [N1,wc1]=buttord(wp1,ws1,rp,rs,s); [B1,A1]=butter(N1,wc1,s); [h1,w1]=freqs(B1,A1); subplot(3,1,2);plot(w1/pi,abs(h1)); title(high anolog filter Magnitude Response) 数字带通滤波器设计实例 2 Matlab求解： 数字带通滤波器程序： wp2=[0.4,0.6];ws2=[0.2,0.8];rp=3;rs=40; [N2,wc2]=buttord(wp2,ws2,rp,rs,s); [B2,A2]=butter(N2,wc2,s); [h2,w2]=freqs(B2,A2); subplot(3,1,3);plot(w2/pi,abs(h2)); title(digital filter Magnitude Response) 数字带通滤波器设计实例 3 Matlab结果 模拟低通滤波器： Butter