num=[1];
den=[1 2 2 1];
sys=tf(num,den);
figure(1); pzmap(sys); title('零极点图');
figure(2); impulse(sys);
w=0:0.1:10*pi;
figure(3); freqs(num,den,w)