f=[0:1:2*f_0];
f_0=2000;
f_s=8000;
w=2*pi*f/f_s;
z=exp(w*j);
B=500;

H=((0.16217*(z.*z-1))./(z.*z+0.67565));
h = ifft(H);
n = [1:1:length(h)];

figure
subplot(2,1,1),plot(f,abs(H))
title('Frequency Response')
ylabel('|H[z]|')
xlabel('f [Hz]')


subplot(2,1,2),stem(n(1:25),h(1:25))
title('Impulse Response')
ylabel('h[n]')
xlabel('n')