26t.m

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

Rp=1;
Rs=15;
wp1=0.2*pi;
ws1=0.3*pi;
T=0.001;
wp=(2/T)*tan(wp1/2);
ws=(2/T)*tan(ws1/2);
[n,wn]=ellipord(wp,ws,Rp,Rs,'s')  %从此处可以计算阶数n
[b,a]=ellip(n,Rp,Rs,wn,'low','s')  %由b,a的值可以得到系统函数
[bz,az]=bilinear(b,a,1/T);
[b0,B,A]=dir2cas(bz,az)
[db,mag,pha,grd,w]=freqz_m(bz,az);
subplot(2,1,1);
plot(w/pi,db);
axis([0,1,-65,10]);
xlabel('');
ylabel('幅频');
title('幅频相应(相对幅度)')
subplot(2,1,2);
plot(w/pi,pha);
xlabel('频率（单位：pi）');
ylabel('相位');
title('相频相应')