Problem solving ODE system in Simulink

1 次查看（过去 30 天）
Oscar Espinosa
Oscar Espinosa 2020-11-7
评论： Oscar Espinosa 2020-11-7
Hi, I'm having some troubles when trying to solve a system of ODE.
This is the equation that I'm trying to solve: with initial conditions .
When solving the equation numerically I obtained the next plot in the interval from to .
I tried to do the same with Simulink, but the resulting graph from the scope its the next one:
The Simulink diagram is the next:
Thanks in advance.
Oscar Espinosa

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

采纳的回答

Stephan
Stephan 2020-11-7
I think something went wrong, when you solved numeric:
syms y(t)
eq = diff(y,t,4) + 3 * diff(y,t,2) - sin(t) * diff(y,t,1) + 8 * y == t^2;
[V,S] = odeToVectorField(eq);
fun = matlabFunction(V,'Vars',{'t','Y'});
[tsol, ysol] = ode45(fun,[0 20], [1 2 3 4]);
plot(tsol, ysol(:,1))
result is:
The equation in Simulink:
gives:
Which is the same result. I think you have to switch the x0-Inputs with the signal inputs of your 4 integrators to solve the issue.
  1 个评论
Oscar Espinosa
Oscar Espinosa 2020-11-7
Hi, Stephan
That helped to me, thank you. By the way I'll like to show how I tried to solve the ODE like a system of ODE, so:
xconds = [1;2;3;4];
function dxdt = ode_4th_order(t,x)
dxdt(1,1) = x(2,1);
dxdt(2,1) = x(3,1);
dxdt(3,1) = x(4,1);
dxdt(4,1) = -3*x(3,1)+sin(t)*x(2,1)-x(1,1)+t^2;
end
[t,x] = ode45('ode_4th_order',tspan,xconds);
plot(t,x(:,1));
But this give a completely different graph, the one that you see in the original post. Does this approach work? I mean, there is something that I did wrong or something missing?

请先登录，再进行评论。

更多回答（0 个）

类别

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

标签

产品


版本

R2016a

Community Treasure Hunt

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

Start Hunting!

Translated by