34t.m

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

Rp=0.3;Rs=15;T=1/20000;
wp=2*pi*1400*T;ws=2*pi*1000*T;
wp1=(2/T)*tan(wp/2);ws1=(2/T)*tan(ws/2);
[n,wn]=cheb2ord(wp1,ws1,Rp,Rs,'s');
[b,a]=cheby2(n,Rs,wn,'high','s');
[bz,az]=bilinear(b,a,1/T);
[db,mag,pha,grd,w]=freqz_m(bz,az);plot((w*20000)/(2*pi),db);axis([0,1800,-65,5]);grid