SISTEMA DE TANQUES INTERCONECTADOS

Código en Matlab para hallar los puntos de operación

A1=10;

As=1;

g=10;

k=As*sqrt(2*g);

syms h1 h2 qi

definimos las ecuaciones en una variable (sistema no lineal)

sis(1,1)= (1/A1)*(qi-k*sqrt(h1-h2));

sis(1,2)= (1/A1)*(k*sqrt(h1-h2)-k*sqrt(h2));

PUNTOS DE OPERACION DEL SISTEMA

Ph2=2;

solución de la ecuación para obtener los puntos de operación de h1 y qi

s=solve(sis,h1,qi);

punto de operación de h1 y qi sustituyendo en s el valor de ph2

Ph1= double(subs(s.h1,h2,Ph2));

Pqi= double(subs(s.qi,h2,Ph2));

Código en Matlab para hallar  las matrices A, B, C y D

declaración de estradas, salidas y estados del sistema

entradas= qi;

salidas=h2;

estados=[h1 h2];

MATRIZ JACOBIANA

A=jacobian(sis,estados);

B=jacobian(sis,entradas);

C=jacobian(salidas,estados);

D=jacobian(salidas,entradas);

Sustitución de h1 por Ph1 ..... Pqi para obtener las matrices numéricas

A=double(subs(A,[h1 h2 qi],[Ph1 Ph2 Pqi]));

B=double(subs(B,[h1 h2 qi],[Ph1 Ph2 Pqi]));

C=double(subs(C,[h1 h2 qi],[Ph1 Ph2 Pqi]));

D=double(subs(D,[h1 h2 qi],[Ph1 Ph2 Pqi]));

obtención de la función de transferencia con respecto h1 con las matrices A B C y D

[num den]=ss2tf (A,B,C,D);

G=tf(num,den);