shiyan3.m

EDA,VHDL文档

oc1=15
T=0.1;Fs=1/T;
wc1=oc1*T;
oc2=c2*tam(wc1/2);
N=3;

[b1,a1]=butter(N,oc1,'s');
[b2,a2]=butter(N,oc2,'s');
[bz1,az1]=bilinear(b1,a1,Fs);
[bz2,az2]=impinvar(b2,a2,Fs);
sys1=tf(bz1,az1,T)
sys2=tf(bz2,az2,T)
w=linspace(-pi,pi,512);
h1=freqz(bz1,az1,w);
h2=freqz(bz2,az2,w);
mag1=20*log10(abs(h1+eps));
mag2=20*log10(abs(h2+eps));
subplot[2,1,1];plot(w,mag1);grid;
axis([0,2*pi,-100,5]);
set(gca,'ytick',[(-100:20:-20),-3,5]);
subplot[2,1,1];plot(w,mag2);grid;
axis([0,2*pi,-100,5]);
set(gca,'ytick',[(-100:20:-20),-3,5]);