wp=(22/59)*pi;wr=(26/59)*pi;
tr_width=wr-wp;
N=ceil(10*pi/tr_width);
n=0:1:N-1;
wc=(wr+wp)/2;
hd=ideal_lp(wc,N);
w_ham=(hamming(N))';
h=hd.*w_ham;
[db,mag,pha,w]=freqz_m(h,[1]);
delta_w=2*pi/59000;
subplot(221);stem(n,hd);
subplot(222);stem(n,w_ham);
subplot(223);stem(n,h);
subplot(224);plot(w/pi,db);
ap=-(min(db(1:1:wp/delta_w+1)));
ar=-(max(db(wr/delta_w:1:59000/2)));