matlab編程 頻率采樣法設(shè)計(jì)濾波器

wp=0.2*pi;



ws=0.35*pi;



wc=wp;



N=34;



alpha=(N-1)/2;



k=0:N-1;



wk=(2*pi)/N*k;



Hk=[ones(1,4),0.5925,0.1099,zeros(1,23),-0.1099,-0.5925,-ones(1,3)];



angH=-alpha*(2*pi)/N*k;



H=Hk.*exp(i*angH);



h=real(ifft(H,N));



[db,mag,pha,grd,w]=freqz_m(h,1);



[Hr,ww,b,l]=Hr_Type2(h);



plot(w/pi,db,'r');



axis([0,1,-100,2]);