Solution to ODE in matlab or simulink

Question

Soheil 2019-10-4

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/483655-solution-to-ode-in-matlab-or-simulink

回答： 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.

2 个评论
显示无隐藏无

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 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.

请先登录，再进行评论。

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

Answer 1

darova 2019-10-5

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/483655-solution-to-ode-in-matlab-or-simulink#answer_394924

在 MATLAB Online 中打开

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