用双线性变换法设计一个巴特沃斯低通IIR数字滤波器

clear all;
A=[0.090360.09036*20.09036];
B1=[1-1.26860.7051];
B2=[1-1.01060.3583];
B3=[1-0.90440.2155];
[H1w]=freqz（AB150）;
[H2w]=freqz（AB250）;
[H3w]=freqz（AB350）;
h=H1.*H2;
H=h.*H3;
i=1:50;
a=20*log10（abs（H（i））+eps/abs（H（1）））;
b=20*log10（abs（H1（i））+eps/abs（H1（1）））;
c=20*log10（abs（H2（i））+eps/abs（H2（1）））;
d=20*log10（abs（H3（i））+eps/abs（H3（1）））;

figure（1）;
subplot（221）;
plot（w/pia）;
axis（[00.5-5010]）;
title（‘滤波器的幅频响应特性曲线‘）;
xlabel（‘w/\pi‘）;ylabel（‘20lg|H（ej0）/H（ejw）|‘）;
grid on;

subplot（222）;
plot（w/pib）;
axis（[00.5-5010]）;
title（‘第一级二阶滤波器幅频响应特性曲线‘）;
xlabel（‘w/\pi‘）;ylabel（‘20lg|H1（ej0）/H1（ejw）|‘）;
grid on;

subplot（223）;
plot（w/pib）;
axis（[00.5-5010]）;
title（‘第二级二阶滤波器幅频响应特性曲线‘）;
xlabel（‘w/\pi‘）;ylabel（‘20l