Symbolic Vector ODE numerical solution

13 vues (au cours des 30 derniers jours)
Leonardo Dipilato
Leonardo Dipilato le 14 Juin 2020
Commenté : Walter Roberson le 15 Juin 2020
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 commentaire
Walter Roberson
Walter Roberson le 15 Juin 2020
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.

Connectez-vous pour commenter.

Connectez-vous pour répondre à cette question.

Réponses (1)

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

Connectez-vous pour commenter.

Catégories

En savoir plus sur Calculus dans Help Center et File Exchange

Tags

Community Treasure Hunt

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

Start Hunting!

Translated by