题目：根据下列参数完成IIR和FIR数字滤波器设计：通带范围300HZ~500HZ.pdfVIP

题目：根据下列参数完成 IIR 和 FIR 数字滤波器设计： 通 带范围 300HZ~500HZ 题目 : 根据下列参数完成 IIR 和 FIR 数字滤波器设计 : 通带范围 300Hz~500Hz 带内最大衰减 Rp=-3dB 阻带范围 250Hz550Hz带内最小衰减 Rs=-40dB 采样频率 Fs=2000Hz 要求 : 1、分别完成 IIR 和 FIR 滤波器的设计 2、IIR 设计不可使用 butter 、cheby1、cheby2 和 ellip 这四个完全设计函数 3、谈谈自己对两种滤波器设计的感受 一、 IIR 数字滤波器的设计 : ( 一) 设计 Butterworth 型带通滤波器 1. 冲激响应法 clear; close all; Fs=2000; Wp(1)=300*2*pi; Wp(2)=500*2*pi; Ws(1)=250*2*pi; Ws(2)=550*2*pi; Rp=3; Rs=40; [n,Wn]=buttord(Wp,Ws,Rp,Rs,s); [z,p,k]=buttap(n); Wo=sqrt(Wn(1)*Wn(2)); [b0,a0]=zp2tf(z,p,k); [b,a]=lp2bp(b0,a0,Wo,Wn(2)-Wn(1)); [bz,az]=impinvar(b,a,Fs); freqz(bz,az) 阶数 n=15 2. 双线性变换法 clear; close all; Fs=2000; Wp(1)=2*Fs*tan(300*pi/Fs); Wp(2)=2*Fs*tan(500*pi/Fs); Ws(1)=2*Fs*tan(250*pi/Fs); Ws(2)=2*Fs*tan(550*pi/Fs); Rp=3; Rs=40; [n,Wn]=buttord(Wp,Ws,Rp,Rs,s); [z,p,k]=buttap(n); Wo=sqrt(Wn(1)*Wn(2)); [b0,a0]=zp2tf(z,p,k); [b,a]=lp2bp(b0,a0,Wo,Wn(2)-Wn(1)); [bz,az]=bilinear(b,a,Fs); freqz(bz,az) n =12 Wn =1.0e+003 * 2.0237 4.0285 ( 二)Chebyshev I 型带通滤波器 clear; close all; Fs=2000; Wp(1)=300*2*pi; Wp(2)=500*2*pi; Ws(1)=250*2*pi; Ws(2)=550*2*pi; Rp=3; Rs=40; [n,Wn]=cheb1ord(Wp,Ws,Rp,Rs,s); Wo=sqrt(Wn(1)*Wn(2)); % 求带通滤波器中心频率 [z,p,k]=cheb1ap(n,Rp); [b0,a0]=zp2tf(z,p,k); [b,a]=lp2bp(b0,a0,Wo, Wn(2)-Wn(1)); % 带通滤波 器带宽 [bz,az]=impinvar(b,a,Fs); freqz(bz,az) n =7 Wn = 1.0e+003 * 1.8850 3.1416 Wo=2.4335e+003 ( 三) 完全设计 椭圆带通滤波器 clear; close all; clc; Fs=2000; Wp(1)=300*2/Fs; % 数字归一化频率， 300*2*pi/(Fs*pi) Wp(2)=500*2/Fs;% 与 后面数字滤波器设计对应 Ws(1)=250*2/Fs; Ws(2)=550*2/Fs; Rp=3; Rs=40; [n,Wn]=ellipord(Wp,Ws,Rp,Rs); [b,a]=ellip(n,Rp,Rs,Wn,bandpass) freqz(b,a,Fs) n =4 Wn =0.3000 0.5000 二、 F