clear all; close all; clc;
pkg load control % paquete donde estan las funciones lsim, step, impulse, ss y tf
% Tanque agitado V = 5m3, v = 1 m3/min 

% para trabajar en el dominio de Laplace primero hay que definir la variable s 
% como un polinomio y = s , esto es con raiz en 0
s = tf('s');

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

% para las simulaciones en el dominio del tiempo podemos usar  lsim  pero 
% previamente tenemos que haber definido el sistema lineal

t = 0:0.1:20; 
Du = 0.5* ones(length(t),1); % altura del escalón definida para cada tiempo
y=lsim(g,Du,t);

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

% otra forma de simular un escalón es con el comando step
Du = 0.5; % sobreescribimos altura del escalón
y=Du* step(g,t);

pause

y = impulse(g,t);

plot(t,y)
title('sistema de primer orden - respuesta a un pulso unitario')
xlabel('t')
ylabel('y')