a812.m

很多matlab的源代码

wp = 0.2*pi;
ws = 0.3*pi;
Rp = 1;
As = 15;
T = 1;
OmegaP = wp/T;
OmegaS = ws/T;
[cs,ds] = afd_chb1(OmegaP,OmegaS,Rp,As);
[b,a] = imp_invr(cs,ds,T);
[C,B,A] = dir2par(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
xlabel('w');
ylabel('样本');