Solution to ODE in matlab or simulink

1 visualización (últimos 30 días)
Soheil
Soheil el 4 de Oct. de 2019
Respondida: darova el 5 de Oct. de 2019
Hi,
I would be thankful if someone gives a comment on how to solve the following equation in MATLAB/Simulink with ODE45 or ODE15s...It's equation number 3 in this publication [Link]. (Picture of the equation is attached).
It's basically a material balance equation of solid particles in ractor.
Thanks in advance.
matlabequation.jpg
  2 comentarios
Sulaymon Eshkabilov
Sulaymon Eshkabilov el 5 de Oct. de 2019
Hi
Your equation as it does not appear to be correct. Could you please post the correct one?
Just a guess as it looks like a 1st order ODE BVP that means you'd need to employ bvp solver
Soheil
Soheil el 5 de Oct. de 2019
Thanks for the comment; I will check the BVP solver.
I post a picture of equations from a paper below.
They've mentioned at the end that ODE45 is used.
A general guideline toward solving this would be super helpful.
matlabequation.jpg

Iniciar sesión para comentar.

Respuestas (1)

darova
darova el 5 de Oct. de 2019
I tried to solve this equation first
What are constants values?
% opengl software
G = 2;
V = 2.4;
Q = 5;
J = 1;
n0 = 1.1;
num1 = 30;
num2 = 40;
n = zeros(num1,num2);
L = linspace(0,30,num1);
t = linspace(0,100,num2);
dL = L(2) - L(1);
dt = t(2) - t(1);
for i = 1:100
n(end,:) = J/G;
n(:,1) = n0;
[dnt, dnL] = gradient(n,dt,dL);
n = -V/Q*(dnt+G*dnL);
surf(t,L,n)
pause(0.1)
end
The result
gif_animation.gif

Categorías

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

Productos


Versión

R2019b

Community Treasure Hunt

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

Start Hunting!

Translated by