脉冲响应不变法设计iir数字滤波器matlab讲课讲稿.pdf

实验三 脉冲响应不变法设计 IIR 数字滤波器 一、 实验目的 1. 掌握利用脉冲响应不变法设计 IIR数字滤波器的原理和具体方法。 2. 加深理解数字滤波器与连续时间滤波器之间的技术指标转化。 3. 掌握脉冲响应不变法设计 IIR数字滤波器的优缺点及使用范围。 二、 实验内容 1. 利用巴特沃思模拟滤波器 ,通过脉冲响应不变法设计巴特沃思数字滤波器 ,数 字滤波器的技术指标为 采样周期为 T=2 程序代码 T=2; %设置采样周期为 2 fs=1/T; %采样频率为周期倒数 Wp=0.25*pi/T; Ws=0.35*pi/T; %设置归一化通带和阻带截止频率 Ap=20*log10(1/ 0.9); As=20*log10(1/ 0.18); %设置通带最大和最小衰减 [N,Wc]=buttord(Wp,Ws,Ap,As,s); %调用 butter 函数确定巴特沃斯滤波器阶数 [B,A]=butter(N,Wc,s); %调用 butter 函数设计巴特沃斯滤波器 W=linspace(0,pi,400*pi); %指定一段频率值 hf=freqs(B,A,W); %计算模拟滤波器的幅频响应 subplot(2,1,1); plot(W/pi,abs(hf)/abs(hf(1))); %绘出巴特沃斯模拟滤波器的幅频特性曲线 grid on; title( 巴特沃斯模拟滤波器 ); xlabel(Frequency/Hz); ylabel(Magnitude); [D,C]=impinvar(B,A,fs); %调用脉冲响应不变法 Hz=freqz(D,C,W); %返回频率响应 subplot(2,1,2); plot(W/pi,abs(Hz)/abs(Hz(1))); %绘出巴特沃斯数字低通滤波器的幅频特性曲线 grid on; title( 巴特沃斯数字滤波器 ); xlabel(Frequency/Hz); ylabel(Magnitude); 进行试验 输出图像： 巴 特 沃 斯 模 拟 滤 波 器 1 e d u t i n 0.5 g a M 0 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 Frequency/Hz 巴 特 沃 斯 数 字 滤 波 器 1 e d u t i n 0.5 g a M 0 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1