% Aproximacion por polos dominantes

clear
close all
s=tf('s');
format long
t=0:0.001:25;

Gos=136/(s^4 + 18*s^3 + 87*s^2 + 70*s + 136);

[Zos,Pos,Kos]=zpkdata(Gos,'v');

% Ceros=0;
% Polos:
p1=-9.154278363562772;
p2=-8.336925311837442;
p3=-0.254398162299890 + 1.310452758155082i;
p4=-0.254398162299890 - 1.310452758155082i;

figure; step(Gos)
figure; pzmap(Gos)

wn=abs(p3); xita=abs(real(p3))/wn;

Gor1=wn^2/((s - p3)*(s - p4));
Gor2=wn^2/(s^2 + 2*xita*wn*s + wn^2);

figure; step(Gos,Gor)
legend('G orden superior','G orden reducido')

figure; pzmap(Gos,Gor)
legend('G orden superior','G orden reducido')