ex3060.m

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

wplp=0.2*pi;
wslp=0.3*pi;
Rp=1;
As=15;
T=1;Fs=1/T;
OmegaP=(2/T)*tan(wplp/2);
OmegaS=(2/T)*tan(wslp/2);
ep=sqrt(10^(Rp/10)-1);
Ripple=sqrt(1/(1+ep*ep));
Attn=1/(10^(As/20));
[cs,ds]=afd_chb1(OmegaP,OmegaS,Rp,As);
[blp,alp]=bilinear(cs,ds,T);
wphp=0.6*pi;
% LP - HP frequency-band transformation:
alpha=-(cos((wplp+wphp)/2))/(cos((wplp-wphp)/2))
Nz=-[alpha,1];Dz=[1,alpha];
[bhp,ahp]=zmapping(blp,alp,Nz,Dz);
[C,B,A]=dir2cas(bhp,ahp)
% Plots
[db1,mag1,pha1,grd1,w]=freqz_m(blp,alp);
subplot(221)
plot(w/pi,mag1);axis([0,1,0,1]);grid
subplot(222)
plot(w/pi,db1);axis([0 1 -30 0]);grid
[dbh,magh,phah,grdh,w]=freqz_m(bhp,ahp);
subplot(223)
plot(w/pi,magh);axis([0,1,0,1]);grid
subplot(224)
plot(w/pi,dbh);axis([0 1 -30 0]);grid