Symbolic Vector ODE numerical solution

1 次查看(过去 30 天)
I'm trying to solve an ODE using vectorized symbolic calculus.
N = 4;
syms theta(t) [1 N]
Q = tril(ones(N));
theta = theta*Q';
eq = diff(theta, t, t) - cos(theta);
Dt = diff(theta, t);
How can I use this with ode45 or other solvers?
  1 个评论
Walter Roberson
Walter Roberson 2020-6-15
To use that with ode45 or similar, I recommend reading the first example for odeFunction() as it shows the steps to use to prepare symbolic differential equations for use with ode45 or similar.

请先登录,再进行评论。

回答(1 个)

Paresh yeole
Paresh yeole 2020-6-14
in your equation for 'eq', if you are looking for double differentiation then it should be
eq = diff(theta, t, 2) - cos(theta);
  1 个评论
Walter Roberson
Walter Roberson 2020-6-14
diff(theta, t, 2) and diff(theta, t, t) are treated the same, both become diff(theta, t, t) on output.

请先登录,再进行评论。

类别

Help CenterFile Exchange 中查找有关 Symbolic Math Toolbox 的更多信息

Community Treasure Hunt

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

Start Hunting!

Translated by