Tarea 3 C
Tarea 3 C
Tarea 3 C
function p21
clear
clc
t0=[0 50];
h0=0;
[t,hs]=ode45(@p2,t0,h0);
plot(t,hs)
end
%estado estacionario
function dhsdt=p2(~,hs)
As=10;
F1s=30;
F2s=10;
Bs=11.619;
dhsdt=(F1s+F2s-Bs*sqrt(hs))/As;
hold on
end
1
function p2v
clear
clc
t0=[36 50];
h0=11.84;
[t,h]=ode45(@p5,t0,h0);
plot(t,h)
end
function dhdt=p5(~,h)
A=10;
F1=31;
F2=10;
B=11.619;
dhdt=(F1+F2-B*sqrt(h))/A;
end
2
2. Considere el caso de tres tanques de nivel colocados en serie.
a) Desarrolle el modelo matematico de este sistema dinamico sabiendo que se ali-
mentan dos corrientes al primer tanque (F y F ) y que las corrientes de salida son
proporcionales al nivel del tanque.
b) Integre las ecuaciones diferenciales utilizando MATLAB usando como condiciones
inciales un valor 50% del valor en el estado estacionario.
function P2T3
tspan=[0,100];
h10=0;
h20=2;
h30=1;
subplot(1,1,1),plot(h(:,1),t)
subplot(1,1,2),plot(h(:,2),t)
subplot(1,1,3),plot(h(:,3),t)
end
function dhdt=P2T(~,h)
x=h(1);
T=h(2);
P=h(3);
%Datos calculados
%Datos
a1=5;
a2=10;
a3=6;
3
F1=10;
F2=10;
A1=10;
A2=15;
A3=12;
%Ecuacion de Diseño
dhdt(1)=(F1+F2-(a1*x))/A1;
dhdt(2)=((a1*x)-(a2*T))/A2;
dhdt(3)=((a2*T)-(a3*P))/A3;
dhdt=dhdt';
end
4
dx 1 2 dx 6 8
a) = x d) = x
dt 3 4 dt 2 2
dx dx 1 1
1 2 = x
b) = x e) dt 5 3
dt 3 4
dx 1 3 dx 2 6
c) = x f) = x
dt 4 2 dt 1 2
.
function p3t3
clear
clc
t0=[0 10];
x10=-30;
x20=0;
x0=[x10 x20] ;
[t,x]=ode45(@p3t,t0,x0);
plot(t,x(:,1),t,x(:,2))
plot(x(:,1),x(:,2))
end
function dxdt=p3t(~,x)
global a11 a12 a21 a22
a11=1;%input('Valor a11= ');
a12=-2;%input('Valor a12= ');
a21=3;%input('Valor a21= ');
a22=-4;%input('Valor a22= ');
dxdt(1,:)=a11*x(1,:)+a12*x(2,:);
dxdt(2,:)=a21*x(1,:)+a22*x(2,:);
end
a)
5
_____________________________________
b)
6
_____________________________________
c)
__________________________________ 7
d)
_______________________________
e)
_______________________________________
8
f)
______________________________________