[n1,wn1]=buttord(wp,ws,Ap,As);%估计巴特沃什滤波器的阶数 %wn1为归一化截止频率
[num,den]=butter(n1,wn1);%设计巴特沃什滤波器 %计算单位为dB的gain函数
%在单位圆上将函数转化为256个等份 w=0:pi/255:pi;%计算相位 h=freqz(num,den,w);
g=20*log10(abs(h));%计算幅度 pha=angle(h); %画图
subplot(1,2,1);plot(w/pi,g);grid%画带通滤波器的幅度响应曲线图 axis([0 1 -60 5]);
xlabel('\omega/\pi');ylabel('gain in db');%横坐标表示频率，纵坐标表示分贝数 title('巴特沃什带通滤波器幅频响应');%写标题
subplot(1,2,2);plot(w/pi,pha);grid%画带通滤波器的相位响应曲线图 xlabel('\omega/\pi');ylabel('phase'); title('巴特沃什带通滤波器相位响应');
3.1.3IIR数字滤波器MATLAB仿真结果与分析
可以得出：
IIR数字滤波器的阶数为：n1=7
3dB截止频率为：wn1=0.4284π，0.6693π 系统函数H(z)的分子系数为：
num=0.0003，0，-0.0019，0，0.0057，0，-0.0095，0.0095，0，-0.0057，0，0.0019，
0，-0.0003
系统函数H(z)的分母系数为：
den=1.0000，1.7451，4.9282，6.1195，9.8134，9.2245，10.4323，7.5154，6.4091，
3.4595，2.2601，0.8470，0.4167，0.0856，0.0299 根据运行结果可写出IIR滤波器的系统函数为