17t.m

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

Rp=0.5;
T=0.001;
ws=200*2*pi*T;
ws1=(2/T)*tan(ws/2);
[b,a]=cheby1(9,Rp,ws1,'high','s');
[bz,az]=bilinear(b,a,1/T);
[db,mag,pha,grd,w]=freqz_m(bz,az);
subplot(2,1,1);
plot(w/pi,db);
grid;
axis([0,1,-400,100]);
title('系统的幅频响应');
subplot(2,1,2);
plot(w/pi,pha);
title('系统的相频响应');