clear all;
% d=wavread('007.wav');
fg=100;
fp=44100;
fp2=fp*fp;
fg=44100;
w=2*pi*fg;
w2=w*w;

b0=w2/(4*fp2)+(sqrt(2)*w)/(2*fp)+1;
b1= 2*(w/(4*fp2))-2;
b2= w2/(4*fp2)-(sqrt(2)*w)/(2*fp)+1;

k=10001;

for n=1:k
   W(n)=(pi/k)*(n-1);
   H(n)=(1-2*exp(-1i*W(n))+exp(-2*1i*W(n)))/(b0+b1*exp(-1i*W(n))+b2*exp(-2*1i*W(n)));
end
plot(W/(2*pi)*fp,abs(H));