Electronica

UNIVERSIDAD NACIONAL
PEDRO RUÍZ GALLO












TRABAJO DE PROGRAMACION


ALUMNO : Campos Silva Franck



DOCENTE :
Ing. Victor Jara Sandoval

CICLO :
2013-II


Lambayeque, Junio del 2014


CONTROL II


En la figura mostrada la válvula automática cv-01 control el ingreso del fluido qi la válvula manual cm-01 regula el flujo de descarga deltanque cónico q0





Datos de entrada:

Flujo de entrada:
Flujo de salida:
Transmisor de nivel: LT-01
Volumen del tanque:
Controlador de nivel: LIC-01
Válvula de control: cv-01Válvula normal: cm-01
Balance de materiales: *A=Área de la base

Realizando la programación el software MATLAB, y mediante parámetros arbitrarios se analizara y hallara las siguientes ecuaciones yvalores:
Ecuación de estado en forma simbólica.
Determinar la matriz de realimentación k.
Los gráficos y diagramas de bloques del sistema se podrán observar por medio del plot del software.%CONTROL 2
clear all
close all
clc
syms V H a1 cv c1 t u x1 x2 s
%INGRESA FUNCION PARA DETERMINANTE JACOBIANO
f=[x2;(cv*u*(-2)*(x1^(-3))*x2)/a1 - (c1*(-3/2)*(x1^(-5/2))*x2)/a1];
V=[x1,x2];U=[u];
subs(jacobian(f,V));
subs(jacobian(f,U));
x1=5.9
x2=0.1
u=0.68
cv= 137
c1=51.46 ; R=1 ; H=3;
a1=(pi*(R/H)^2)
A=subs(jacobian(f,V))
B=subs(jacobian(f,U))
A1=vpa(A,3)
B1=vpa(B,3)
A1=[0, 1.0 ; 0.0213, 0.0164]
B1= [ 0 ; -0.382]
C1=[1,0]
D1=[0]
[num,den]=ss2tf(A1,B1,C1,D1)
p=[1.0000 -0.0164 -0.0213];
r=roots(p)
%SCORE ROOTS = 0.1544 AND -0.318

%EVALUANDO CONTROLM=[B1 A1*B1]
rank(M)
J=[-2,-3];
W=[-0.0164,1;1,0];
T=M*W
%CANONICA PARA CONTROL
AA=inv(T)*A1*T
BB=inv(T)*B1
CC=[1 0]
DD=[0]

K1=acker(AA,BB,J)
AA2=AA-BB*K1
BB2=BB*K1
I=eye(2)EC2=det(s*I-AA2)
r2=roots([1 5 6])%s^2 + 5*s + 6

x0=[5.9,0.1] %x1=5.9//x2=0.1
sys=ss(AA2, I, I, I)
t=0:0.01:10;
x=initial(sys,x0,t);
x1=[1,0]*x';
x2=[0,1]*x';
plot(t,x1,'-',t,x2,'*');
grid;...