eg10_10.m

来自「MATLAB7.0从入门到精通(光盘).对MATLAB 7.0进行了比较全面的介」· M 代码 · 共 11 行

fsamp=8000;
fcuts=[1000 1300 2210 2410];
mags=[0 1 0];
devs=[0.01 0.05 0.01];
[n,wn,beta,ftype]=kaiserord(fcuts,mags,devs,fsamp);
n=n+rem(n,2);
hh=fir1(n,wn,ftype,kaiser(n+1,beta),'noscale');
[H,f]=freqz(hh,1,1024,fsamp);
plot(f,abs(H));
grid on