33t.m

来自「包含50个数字信号处理方面的matlab程序」· M 代码 · 共 8 行

Rp=0.8;Rs=25;T=0.001;fp=300;fs=200;
wp=2*pi*fp*T;ws=2*pi*fs*T;
wp1=(2/T)*tan(wp/2);ws1=(2/T)*tan(ws/2);
[n,wn]=buttord(wp1,ws1,Rp,Rs,'s');
[b,a]=butter(n,wn,'high','s');
[bz,az]=bilinear(b,a,1/T);
[db,mag,pha,grd,w]=freqz_m(bz,az);plot(0.5*w/(T*pi),db);axis([0,500,-300,5]);grid