% FPS week08_ex01

C=0.0; D=0.64; NT=30; NF=512;
while(C <= 2.0)
   % a)  h[n]
   subplot(2,2,1)
   b=[1 0 C]; a=[1 0 D];
   h=impz(b,a,NT);
   stem([0:NT-1],h)

   title('Impulse response')

   % b)   Z-P diagram
   subplot(2,2,3)
   zplane(b,a)

   title('Zero-Pole diagram')

   % c)   Magnitude frequency response
   subplot(2,2,4)
   [H W]=freqz(b,a,NF,'whole');
   plot(W/pi,abs(H))

   title('Magnitude frequency response')
   C=C+1/50; pause;
end