% comparacao entre bilinear e impinvar

Ordem=4;
Ripple=3; %3dB
Fs=2;  %taxa de amostragem = 2Hz

% Frequencia de Corte
CorteHz= 0.5;
CorteRds= 2*pi*CorteHz;
CorteNorm= CorteHz*2/Fs;  % CorteNormalizada igual CorteHz pois Fs=2

[bs as]=cheby1(Ordem,Ripple,CorteRds,'s');
[Hs ws]=freqs(bs,as,0:.01:(2*pi));

[bz az]=cheby1(Ordem,Ripple,CorteNorm);
[Hz wz]=freqz(bz,az);

[Nz Dz]=bilinear(bs,as,pi/2);
[Hzb wzb]=freqz(Nz,Dz);

[bzi azi]=impinvar(bs,as,Fs);
[Hzi wzi]=freqz(bzi,azi);

plot(ws/(2*pi),abs(Hs),wz/pi,abs(Hz),wzb/pi,abs(Hzb)+0.01,wzi/pi,abs(Hzi),[0 0.5 0.5 1],[1 1 0 0])
legend('analog','digital','bilin','impinvar','ideal')