How to solve ODE's and find m(t) using matlab? Urgent!!

1 visualización (últimos 30 días)
Shubham Maurya
Shubham Maurya el 25 de Jul. de 2014
Comentada: Star Strider el 26 de Jul. de 2014
Here is the given ODE for which I need a solution:
Please help me on how to do this!!!
  2 comentarios
Ben11
Ben11 el 25 de Jul. de 2014
Do you have the Symbolic Math Toolbox?
Shubham Maurya
Shubham Maurya el 26 de Jul. de 2014
No, where it is?

Iniciar sesión para comentar.

Respuesta aceptada

Star Strider
Star Strider el 25 de Jul. de 2014
Editada: Star Strider el 25 de Jul. de 2014
Use ode45. It’s five lines of code, including the plot:
mdot = @(t,m) sqrt(2E+5.*((1+(0.4-m)./0.5).^-1.4 - 1));
[t,m] = ode45(mdot, [0 1], 0.001);
figure(1)
semilogy(t, real(m), '-b', t, imag(m), '-g')
grid
  3 comentarios
Shubham Maurya
Shubham Maurya el 26 de Jul. de 2014
My initial value is m(0)=0.4, and domain of t is [0,0.5]
Star Strider
Star Strider el 26 de Jul. de 2014
Since I did’t have other information, I chose m(0)=0.001 and ran it from [0 1] in time. With your m(0)=0.4 and time span [0 0.5] the call to ode45 is:
[t,m] = ode45(mdot, [0 0.5], 0.4);
and the entire code is now:
mdot = @(t,m) sqrt(2E+5.*((1+(0.4-m)./0.5).^-1.4 - 1));
[t,m] = ode45(mdot, [0 0.5], 0.4);
figure(1)
plot(t, m, '-*b')
grid
axis([xlim 0 0.5])
With the square root, I’m somewhat surprised that there aren’t two solutions, for instance ±0.4 but the routine produces only one. Yours is a nonlinear equation, and the Symbolic Math Toolbox cannot solve it. (I have it, and I tried that to see what it would do. There is no analytic solution.)
The output doesn’t change from the initial conditions, and is uniformly 0.4 from 0 to 0.5.

Iniciar sesión para comentar.

Más respuestas (0)

Categorías

Más información sobre Ordinary Differential Equations en Help Center y File Exchange.

Productos

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by