ex3052.m

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

wp=0.2*pi;ws=0.3*pi;Rp=1;As=15;
T=1;
OmegaP=wp/T;
OmegaS=ws/T;
ep=sqrt(10^(Rp/10)-1);
Ripple=sqrt(1/(1+ep*ep));
Attn=1/(10^(As/20));
% Analog Butterworth Prototype Filter Calculation:
[cs,ds]=afd_butt(OmegaP,OmegaS,Rp,As);
[b,a]=imp_invr(cs,ds,T)
% 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