a818.m

很多matlab的源代码

wp = 0.2*pi;
ws = 0.3*pi;
Rp = 1;
As = 15;
T = 1;
Fs = 1/T;
OmegaP = (2/T)*tan(wp/2);
OmegaS = (2/T)*tan(ws/2);
[cs,ds] = afd_chb1(OmegaP,OmegaS,Rp,As);
[b,a] = bilinear(cs,ds,Fs);
[C,B,A] = dir2cas(b,a)
[db,mage,pha,grd,w] = freqz_m(b,a);
subplot(2,3,4);
plot(w/pi,db);grid
xlabel('w');
ylabel('dB');
subplot(2,3,1);
plot(w/pi,mage);grid
xlabel('w');
ylabel('|H|');
axis([0 1 0 1.1])
subplot(2,3,2);
plot(w/pi,pha);grid
xlabel('w');
ylabel('相位');
axis([0 1 min(pha) max(pha)]);
subplot(2,3,5);
plot(w/pi,grd);grid
axis([0 1 0 15]);
xlabel('w');
ylabel('样本');