b=[12];
a=[143];
subplot（511）;
zplane（ba）;
legend（‘零點‘‘極點‘）;
t=-0:0.02:10;
h=impulse（bat）
subplot（513）;
plot（th