ex3051.m

the file contain many matlab signal source code and many example .

Wp=0.2*pi;Ws=0.3*pi;Rp=1;As=16;
Ripple=10^(-Rp/20);Attn=10^(-As/20);
[b,a]=afd_elip(Wp,Ws,Rp,As);
[C,B,A]=sdir2cas(b,a)
[db,mag,pha,w]=freqs_m(b,a,0.5*pi);
[ha,x,t]=impulse(b,a);
% plots
subplot(221)
plot(w/pi,mag);axis([0,0.5,0,1.1]);grid
subplot(222)
plot(w/pi,db);axis([0,0.5,-30,1]);grid
subplot(223)
plot(w/pi,pha/pi);axis([0,0.5,-1,1]);grid
subplot(224)
plot(t,ha,[0,max(t)],[0,0])
axis([0,max(t),min(ha),max(ha)])