examp1.m

来自「滤波器的设计原理与实现」· M 代码 · 共 11 行

N=65;k=0:(N-1)/2;
Wm=2*pi*k./N;
Ad(1:(N+1)/2)=1;
Ad(18)=0.5886;Ad(19)=0.1065;Ad(20:33)=0;
Hd=Ad.*exp(-j*0.5*(N-1)*Wm);
Hd=[Hd conj(fliplr( Hd(2:(N+1)/2) ) )];
h=real(ifft(Hd));
w=linspace(0,pi-0.1,1000);
H=freqz(h,[1],w);
plot(w/pi,20*log10(abs(H)));grid;