clc; clear all; close all; 
pkg load control;
% Tanque agitado V = 5m3, v = 1 m3/min 
s = tf('s');

k = 1 ; % ganancia (adimensional en este caso)
tau = 5 ; % constante de tiempo (min)
g = k/(tau*s+1);

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

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

pause

[A, B, C, D] = tf2ss(g);
sl1 = ss(A, B, C, D); disp(sl1)
pause
y = Du*step(sl1,t);

plot(t,y)
title('idem pero con state-space')
xlabel('t')
ylabel('y')

pause
disp(ss2tf(sl1))