Solution to ODE in matlab or simulink

1 次查看（过去 30 天）
Soheil
Soheil 2019-10-4
回答： darova 2019-10-5
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 个评论
Sulaymon Eshkabilov
Sulaymon Eshkabilov 2019-10-5
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 2019-10-5
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

请先登录，再进行评论。

请先登录，再回答此问题。

回答（1 个）

darova
darova 2019-10-5
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

类别

在 Help Center 和 File Exchange 中查找有关 Ordinary Differential Equations 的更多信息

标签

产品


版本

R2019b

Community Treasure Hunt

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

Start Hunting!

Translated by