shiyan1.m

来自「EDA,VHDL文档」· M 代码 · 共 20 行

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)');