画零极点 b=[1 2 1]; a=[1 -0.5 -0.005 0.3]; figure(1); zplane(b, a); 画脉冲响应 num=[0 1 2 1]; den=[1 -0.5 -0.005 0.3]; h=impz(num, den);%脉冲响应 figure(2); stem(h, '.');%画脉冲响应 画频率响应 [H, w]=freqz(num, den);%频率响应 figure(3); plot(w/pi, abs(H));