clc; clear all; close all; 
pkg load control;
% Sistema de segundo orden
s = tf('s');

k = 1; % ganancia (adimensional en este caso)
tau = 5; % constante de tiempo (min)
z = input('factor de dumping = ');

g = k/(tau^2*s^2+2*z*tau*s+1);

% para las simulaciones en el dominio del tiempo podemos usar lsim o step

t = 0:0.1:100; Du = 0.5; % altura del escalón
y=Du*step(g,t);

clf
plot(t,y)
title('sistema de segundo orden - respuesta a un escalón de altura 0.5')
xlabel('t')
ylabel('y')

pause
y = impulse(g,t);
clf
plot(t,y)
title('sistema de segundo orden - respuesta a un pulso unitario')
xlabel('t')
ylabel('y')