32t.m

包含50个数字信号处理方面的matlab程序

clc;
clear all;
Rp=0.5;Rs=60;T=1/200;
wp1=60*2*pi*T;wp2=80*2*pi*T;ws1=55*2*pi*T;ws2=85*2*pi*T;
wp3=(2/T)*tan(wp1/2);wp4=(2/T)*tan(wp2/2);
ws3=(2/T)*tan(ws1/2);ws4=(2/T)*tan(ws2/2);
wp=[wp3,wp4];ws=[ws3,ws4];
[n,wn]=cheb2ord(wp,ws,Rp,Rs,'s');
[b,a]=cheby2(n,Rs,wn,'bandpass','s');
[bz,az]=bilinear(b,a,1/T);
[db,mag,pha,grd,w]=freqz_m(bz,az);plot((w*200)/(2*pi),db);axis([0,120,-100,2]);