河北大学工商数字信号处理实验四.docVIP

实验4 数字滤波器的设计 一、实验目的 （1）加深对脉冲响应不变法、双线性变换法、窗函数法和频率采样法的了解； （2）掌握MATLAB进行滤波器设计的子函数。 二、实验内容 1. 利用脉冲响应不变法，用巴特沃斯滤波器原型设计一个低通滤波器，满足：，采样频率为10000Hz。 fs=10000; T=1/fs; Wp=0.2*pi/T; Ws=0.3*pi/T; Ap=1;%通带的衰减值为1 As=15;%阻带的衰减值为15 [N,Wc]=buttord(Wp,Ws,Ap,As,s);%计算ButterWorth滤波器的阶数N和截止频率Wc [B,A]=butter(N,Wc,s);%根据滤波器的阶数N和截止频率Wc计算ButterWorth滤波器分子和分母系数 W=linspace(0,pi,400*pi);%将产生从0到400pi步长为pi的数组 [D,C]=impinvar(B,A,fs);%用冲激响应不变法将模拟滤波器变换成数字滤波器 Hz=freqz(D,C,W); plot(W/pi,abs(Hz)/abs(Hz(1))); grid on; title(巴特沃斯数字滤波); xlabel(Frequency/Hz); ylabel(Magnitude); 2. 设计巴特沃斯高通数字滤波器，满足：通带边界频率为400Hz，阻带边界频率为200Hz，通带衰减小于3dB