Matlab - Sistema EDO

   
Vista:

Sistema EDO

Publicado por Felipe (8 intervenciones) el 23/09/2011 02:20:41
Hola, agradecería mucho su ayuda para resolver este problema, bastante simple quizas, pero mi poca experiencia en matlab me ha perjudicado:

Tengo el siguiente sistema de ecuaciones diferenciales:

dX/dt = -6*A
dY/dt = 6*A -8*B
dZ/dt = 8B

Y debo resolverlo mediante Ode23 o Ode45.

De antemano muchas gracias
Valora esta pregunta
Me gusta: Está pregunta es útil y esta claraNo me gusta: Está pregunta no esta clara o no es útil
0
Responder

Sistema EDO

Publicado por Felipe (8 intervenciones) el 23/09/2011 06:27:57
Corrijo:

dX/dt = -6*X
dY/dt = 6*X -8*Y
dZ/dt = 8Y
Valora esta respuesta
Me gusta: Está respuesta es útil y esta claraNo me gusta: Está respuesta no esta clara o no es útil
0
Comentar
Imágen de perfil de JOSE JEREMIAS CABALLERO

Sistema EDO

Publicado por JOSE JEREMIAS CABALLERO jjcc94@hotmail.com (3440 intervenciones) el 23/09/2011 14:09:36
Hola Felipe.
Sigue corregiendo, aun faltan los datos iniciales. Por eso no he habia contestado.

saludos
JOSE JEREMIAS CABALLERO
ASESOR DE PROYECTOS CON MATLAB
PROFESOR DE METODOS NUMERICOS CON MATLAB
PROGRAMADOR MATLAB
jjcc94@hotmail.com
Valora esta respuesta
Me gusta: Está respuesta es útil y esta claraNo me gusta: Está respuesta no esta clara o no es útil
0
Comentar

Sistema EDO

Publicado por Felipe (8 intervenciones) el 23/09/2011 15:21:59
dX/dt = -6*X
dY/dt = 6*X -8*Y
dZ/dt = 8Y

Para t = 0 , X = 1
t = 0 , Y = 0
t = 0 , Z = 0

Esa es la información con la que cuento

Muchas gracias
Valora esta respuesta
Me gusta: Está respuesta es útil y esta claraNo me gusta: Está respuesta no esta clara o no es útil
0
Comentar
Imágen de perfil de JOSE JEREMIAS CABALLERO

Sistema EDO 23

Publicado por JOSE JEREMIAS CABALLERO jjcc94@hotmail.com (3440 intervenciones) el 23/09/2011 19:10:11
Hola Ivan.
%=====================
% dX/dt=-6*X X(0)=1
% dY/dt=6*X -8*Y Y(0)=0
% dZ/dt=8Y Z(0)=0

function S=ode23_caballero
opciones=odeset('RelTol',1e-5,'AbsTol',[1e-5 1e-5 1e-5]);
[T,V]=ode45(@ecua3,[0 1],[1 0 0 ],opciones);
S=[T V];
figure(gcf)
plot(T,V(:,1),T,V(:,2),T,V(:,3));
[x,y,z]=dsolve('Dx =-6*x', 'Dy = 6*x-8*y', 'Dz=8*y','x(0)=1', 'y(0)=0', 'z(0)=1');
legend(['x(t)=',char(x)],['y(t)=',char(y)],['z(t)=',char(z)])
grid

function dv=ecua3(t,v)
dv=zeros(3,1);
dv(1)=-6*v(1);
dv(2)=6*v(1)-8*v(2);
dv(3)=8*v(2);
%============================


EJECUCION.
>> S=ode23_caballero

S =

0 1.0000 0 0
0.0033 0.9802 0.0195 0.0003
0.0067 0.9608 0.0382 0.0010
0.0100 0.9418 0.0559 0.0023
0.0133 0.9231 0.0729 0.0040
0.0221 0.8761 0.1134 0.0105
0.0308 0.8314 0.1489 0.0197
0.0395 0.7890 0.1798 0.0312
0.0482 0.7488 0.2065 0.0447
0.0570 0.7103 0.2296 0.0601
0.0658 0.6738 0.2493 0.0769
0.0746 0.6392 0.2658 0.0951
0.0834 0.6063 0.2794 0.1143
0.0927 0.5733 0.2911 0.1356
0.1021 0.5420 0.3003 0.1577
0.1114 0.5125 0.3071 0.1804
0.1208 0.4845 0.3119 0.2036
0.1307 0.4566 0.3150 0.2284
0.1406 0.4302 0.3163 0.2535
0.1505 0.4053 0.3161 0.2786
0.1604 0.3819 0.3145 0.3036
0.1710 0.3584 0.3115 0.3301
0.1816 0.3364 0.3073 0.3563
0.1921 0.3157 0.3022 0.3820
0.2027 0.2963 0.2963 0.4074
0.2140 0.2769 0.2892 0.4339
0.2254 0.2587 0.2816 0.4598
0.2367 0.2417 0.2734 0.4849
0.2480 0.2258 0.2649 0.5093
0.2602 0.2099 0.2555 0.5347
0.2724 0.1951 0.2458 0.5591
0.2846 0.1813 0.2361 0.5826
0.2968 0.1685 0.2263 0.6052
0.3100 0.1557 0.2158 0.6285
0.3232 0.1438 0.2054 0.6508
0.3364 0.1329 0.1952 0.6720
0.3496 0.1227 0.1852 0.6921
0.3641 0.1126 0.1746 0.7128
0.3785 0.1032 0.1644 0.7324
0.3929 0.0947 0.1546 0.7508
0.4073 0.0868 0.1451 0.7681
0.4232 0.0789 0.1352 0.7859
0.4391 0.0717 0.1258 0.8025
0.4550 0.0652 0.1169 0.8179
0.4709 0.0593 0.1085 0.8322
0.4887 0.0533 0.0997 0.8470
0.5064 0.0479 0.0915 0.8606
0.5242 0.0431 0.0839 0.8730
0.5419 0.0387 0.0769 0.8844
0.5620 0.0343 0.0695 0.8962
0.5822 0.0304 0.0627 0.9068
0.6023 0.0269 0.0566 0.9165
0.6225 0.0239 0.0510 0.9251
0.6459 0.0207 0.0451 0.9341
0.6693 0.0180 0.0399 0.9421
0.6928 0.0157 0.0352 0.9491
0.7162 0.0136 0.0311 0.9553
0.7412 0.0117 0.0272 0.9611
0.7662 0.0101 0.0237 0.9662
0.7912 0.0087 0.0207 0.9706
0.8162 0.0075 0.0180 0.9745
0.8412 0.0064 0.0157 0.9779
0.8662 0.0055 0.0137 0.9808
0.8912 0.0048 0.0119 0.9834
0.9162 0.0041 0.0103 0.9856
0.9372 0.0036 0.0092 0.9872
0.9581 0.0032 0.0082 0.9887
0.9791 0.0028 0.0072 0.9899
1.0000 0.0025 0.0064 0.9911

y tambien te va a salir una grafica.


saludos
JOSE JEREMIAS CABALLERO
ASESOR DE PROYECTOS CON MATLAB
PROFESOR DE METODOS NUMERICOS CON MATLAB
PROGRAMADOR MATLAB
jjcc94@hotmail.com
Valora esta respuesta
Me gusta: Está respuesta es útil y esta claraNo me gusta: Está respuesta no esta clara o no es útil
0
Comentar

Sistema EDO 23

Publicado por Felipe (8 intervenciones) el 23/09/2011 21:04:04
Muchas gracias, me ha sido de gran ayuda nuevamente

Problema resuelto, y he logrado aplicarlo a otros sistemas.

Reitero mis agradecimientos.
Valora esta respuesta
Me gusta: Está respuesta es útil y esta claraNo me gusta: Está respuesta no esta clara o no es útil
0
Comentar