新建 文本文档.txt

EDA,VHDL文档

Transfer function:
       3 s^2 + 4
-----------------------
s^3 + 2.5 s^2 + 2 s + 2
 
 
Transfer function:
0.6667 z^3 + 2.776e-016 z^2 - 4.996e-016 z + 0.6667
---------------------------------------------------
         z^3 - z^2 + 0.6667 z + 5.583e-017
 
Sampling time: 1
 
Transfer function:
      3 z^2 - 2.789 z + 2.344
-----------------------------------
z^3 - 1.018 z^2 + 0.726 z - 0.08208
 
Sampling time: 1
>> 

shiyan1
clc
clear all
b=[3 0 4 ];a=[1 2.5 2 2];
T=1;Fs=1;
sys=tf(b,a)
[bz1,az1]=bilinear(b,a,Fs);
[bz2,az2]=impinvar(b,a,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);
figure(1);
plot(w/pi,20*log10(abs(h1+eps)/max(abs(h1))));axis([-1,1,-50,0]);grid
xlabel('Normalized freqency(rad/s');
ylabel('Normalized Magnitude(dB)');
figure(2);
plot(w/pi,20*log10(abs(h2+eps)/max(abs(h2))));axis([-1,1,-10,0]);grid
xlabel('Normalized freqency(rad/s)');
ylabel('Normalized Magnitude(dB)');





22222
clc
clear all
close all
omegac2=0.3*pi;Apmax=1;
Armin=60;
omegar1=0.35*pi;
T=1;
wp=[omegac2/pi];wr=[omegar1/pi];
[N,Wc]=buttord(wp,wr,Apmax,Armin);
[b,a]=butter(N,Wc);
w=linspace(0,pi,500);
H=freqz(b,a,w);
mag=20*log10(abs(H+eps));
plot(w/pi,mag);axis([0,2,-110,5]);grid;
set(gca,'ytick',[(-100:20:-20),-3,0])
set(gca,'xtick',[(0:0.2:0.4),0.50,(0.6:0.1:2)]);
sys=tf(b,a,T)




clc
clear all
close all
omegac2=0.3*pi;Apmax=1;
Armin=60;
omegar1=0.35*pi;
T=1;
wp=[omegac2/pi];wr=[omegar1/pi];
[N,Wc]=cheb1ord(wp,wr,Apmax,Armin);
[b,a]=cheby1(N,Apmax,Wc);
w=linspace(0,pi,500);
H=freqz(b,a,w);
mag=20*log10(abs(H+eps));
plot(w/pi,mag);axis([0,2,-110,5]);grid;
set(gca,'ytick',[(-100:20:-20),-3,0])
set(gca,'xtick',[(0:0.2:0.4),0.50,(0.6:0.1:2)]);
sys=tf(b,a,T)