Matlab - Sistema EDO

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

Corrijo:

dX/dt = -6*X
dY/dt = 6*X -8*Y
dZ/dt = 8Y

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
[email protected]

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

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
[email protected]

Muchas gracias, me ha sido de gran ayuda nuevamente

Problema resuelto, y he logrado aplicarlo a otros sistemas.

Reitero mis agradecimientos.