ex3055.m

来自「the file contain many matlab signal sour」· M 代码 · 共 21 行

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);
ep=sqrt(10^(Rp/10)-1);
Ripple=sqrt(1/(1+ep*ep));
Attn=1/(10^(As/20));
[cs,ds]=afd_butt(OmegaP,OmegaS,Rp,As);
[b,a]=bilinear(cs,ds,T);
[C,B,A]=dir2cas(b,a)
% Plots
[db,mag,pha,grd,w]=freqz_m(b,a);
subplot(221)
plot(w/pi,mag);axis([0,1,0,1.1]);grid
subplot(222)
plot(w/pi,db);axis([0,1,-40,5]);grid
subplot(223)
plot(w/pi,pha/pi);axis([0,1,-1,1]);grid
subplot(224)
plot(w/pi,grd);axis([0,1,0,10]);grid