sy23.m

来自「数字信号处理工具箱」· M 代码 · 共 13 行

Rp=3;Rs=15;t=0.001
wp1=0.25*pi/t;wp2=0.45*pi/t;ws1=0.15*pi/t;ws2=0.55*pi/t;
wp=[wp1,wp2];ws=[ws1,ws2];
[n,wn]=buttord(wp,ws,Rp,Rs,'s');         % 求阶数n ,3dB截止频率wn
[z,p,k]=buttap(n);                       % 模拟低通原型
[c,d]=zp2tf(z,p,k);                      % 转化为多项式形式
w0=(wn(1)+wn(2))/2;Bw=wn(1)-wn(2);
[c1,d1]=lp2bp(c,d,w0,Bw);                % 模拟低通到模拟带通
[b,a]=impinvar(c1,d1,1/T);
[db,mag,pha,grd,w]=freqz_m(b,a);         % 求频域响应
subplot(2,1,1);plot(w/pi,db); axis([0.1,0.8,-30,5]);
subplot(2,1,2);plot(w/pi,pha);