//===============================================================
// Scilab como herramienta en Sistemas y Control
//===============================================================

//
// Definiendo un sistema sabiendo su representacion de estado
//

A = [ -0.7 -0.05; 1 0]
B = [  1 0]'
C = [  0  1]
D = [ 0]

[sist_1]=syslin('c',A,B,C,D)   

pause


//
// Verificamos la correcta entrada de las matrices de estado
//

sist_1("A")
sist_1("B")
sist_1("C")
sist_1("D")
sist_1("X0")   // Valor inicial

pause

//
// Hallando la funcion de transferencia
//

[h]=ss2tf(sist_1) 

Num = numer(h)

Den = denom(h)


//
// Caso contrario de querer pasas de funcion de transferencia a 
// representacion de estado usar tf2ss
//

pause

//
// Autovalores (polos) del sistema 
//

spec(A)

// Verificamos polos

roots(Den)

pause

//
// Simulacion en el tiempo
//

t=0:0.05:100

// Respuesta a una funcion escalon

U = ones(1,length(t));

[y1 X]=csim(U,t,sist_1);

plot2d(t, [U', y1'])
//legend(['u(t)';'Salida(t)'],a=2) Corregir esto


pause

xbasc

// Respuesta a una funcion escalon

[y2 X]=csim("step",t,sist_1);
plot2d(t, y2)

pause

xbasc


// Respuesta a una funcion senoidal

U = sin(t);
[y3 X]=csim(U,t,sist_1);
plot2d(t, [y3' X(1,:)' X(2,:)'], style=[color("red"),color("green"),color("blue")])