sy21.m

数字信号处理工具箱

Rp=1.5;Rs=40;T=0.001;fp=0.3;fs=0.2;
wp=2*pi*fp;ws=2*pi*fs;
wp1=(2/T)*tan(wp/2);
ws1=(2/T)*tan(ws/2);
[n,wn]=cheb1ord(wp1,ws1,Rp,Rs,'s');      % 求阶数n ,3dB截止频率wn
[z,p,k]=cheb1ap(n,Rp);                   % 模拟低通原型
[c,d]=zp2tf(z,p,k);                      % 转化为多项式形式
[c1,d1]=lp2hp(c,d,wn);                   % 模拟低通到模拟高通
[b,a]=bilinear(c1,d1,1/T);               % 双线性变换法
[db,mag,pha,grd,w]=freqz_m(b,a);         % 求频域响应
subplot(2,1,1);plot(w/pi,db); axis([0,1,-60,5]);
subplot(2,1,2);plot(w/pi,pha);