I have a system of 2 second order ODEs, which i converted to a system of fist order ODEs, which ode45 (or ode15s) needs as an input
my system input looks like this:
as an output i get the vector , but I need the derivatives and . How can i get those?
I tried using the gradient function as suggested by Star Strider in this question: How to output derivative from ODE 45 but unfortunately that did not help.
After integrating your ODE-system you will have for some , and you have the function for your ODE-equations, f(t,X), one way to calculate the accelerations at times will be to simply call:
dxd2xdyd2ydt(i,:) = f(t_i,X(i,:));
Typically I use a regularly space array for the time in the call to ode45 and use the gradient-function, but the above should work - that is after all what we expect odeNN: to give us a solution X where that is true...
ty for your answer! I dont quite understand what you mean by "function for your ODE- equations" or what the f would be. Could you please elaborate give me simple example?
Please show us the expressions for ODE1 and ODE2. I still don't understand why you are having trouble plugging x and xdot into these expressions to get xdotdot. Why can't you just call ode_exo_ho( ) with the appropriate inputs to get what you want?
If the ode_exo_ho( ) function is vectorized (meaning it will accept vectors and properly calculate a vector result), then yes it is just this easy. If ode_exo_ho( ) is not vectorized, then you would have to call it in a loop, but still pretty easy.
The ease of vectorizing your function depends on what it is doing. Sometimes it is just a simple matter of replacing the operators * and / and ^ with their element-wise counterparts .* and ./ and .^
Impossibile completare l'azione a causa delle modifiche apportate alla pagina. Ricarica la pagina per vedere lo stato aggiornato.
Translated by
Seleziona un sito web
Seleziona un sito web per visualizzare contenuto tradotto dove disponibile e vedere eventi e offerte locali. In base alla tua area geografica, ti consigliamo di selezionare: .
Puoi anche selezionare un sito web dal seguente elenco:
Come ottenere le migliori prestazioni del sito
Per ottenere le migliori prestazioni del sito, seleziona il sito cinese (in cinese o in inglese). I siti MathWorks per gli altri paesi non sono ottimizzati per essere visitati dalla tua area geografica.