WICHTIG: Der Betrieb von goMatlab.de wird privat finanziert fortgesetzt. - Mehr Infos...

Mein MATLAB Forum - goMatlab.de

Mein MATLAB Forum

  
Gast > Registrieren       Autologin?   

Partner:




Forum
      Option
[Erweitert]
  • Diese Seite per Mail weiterempfehlen
     


Gehe zu:  
Neues Thema eröffnen Neue Antwort erstellen

Wie löse ich eine mit Octave lsode gelöste ODE mit ode45?
 

Monk3y
Forum-Newbie

Forum-Newbie

Beiträge: 1
Anmeldedatum: 25.04.17
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 25.04.2017, 14:05     Titel: Wie löse ich eine mit Octave lsode gelöste ODE mit ode45?
  Antworten mit Zitat      
Hallo Zusammen

Ganz neu im Forum und gleich die erste Frage Very Happy
Ich schreibe gerade ein Octave Skript in Matlab um. Alles klappt tip top. Nur mit dem Differentialgleichungs-Löser habe ich noch ein wenig mühe.

Ich habe eine gegebene Funktion (bereits in Matlab code)

Code:
function [ y ] = dgeodesic( x,s )
n = size(x,1)/2;
y = zeros(2*n,1);
p = position(x);
u = velocity(x);
y(1:n,1) = u;
for alpha = 1:n
    Gamma = christoffel(p, alpha);
    y(n+alpha,1) = -u' * Gamma * u;
end
end


In Octave wurde das mit diesem Code gelöst:

Code:
function solution = geodesic(x0, s)
   [solution, istate, msg] = lsode("dgeodesic", x0, s);
   if istate != 2
      printf("integration failed. %s\n", msg);
   endif
   solution = cat(2, s', solution);
endfunction


Wie löst man sowas jetzt mit dem ode45 von Matlab? Ich weiss das einmal der Initialisierung Vektor x0 und der Zeitparameter vertauscht sind. Des weiteren will ode45 ein Zeitintervall während lsode einfach das "s" Schluckte. (s ist bloß N = 800; s = 1 * (0:N)). Ich brauche also irgendwie eine for schleife oder so denke ich.

Hatt da jemand einen Tipp wie ich das machen könnte?
Besten Dank, Grüsse
Monk3y

eidt: der initial-vektor x0 enthält 8 Anfangsbedingungen, dgeodesic liefert 8 gleichungen
Private Nachricht senden Benutzer-Profile anzeigen


Neues Thema eröffnen Neue Antwort erstellen



Einstellungen und Berechtigungen
Beiträge der letzten Zeit anzeigen:

 Du kannst Beiträge in dieses Forum schreiben.
Du kannst auf Beiträge in diesem Forum antworten.
Du kannst deine Beiträge in diesem Forum nicht bearbeiten.
Du kannst deine Beiträge in diesem Forum nicht löschen.
Du kannst an Umfragen in diesem Forum nicht mitmachen.
Du kannst Dateien in diesem Forum posten
Du kannst Dateien in diesem Forum herunterladen
.





 Impressum  | Nutzungsbedingungen  | Datenschutz | FAQ | goMatlab RSS Button RSS

Hosted by:


Copyright © 2007 - 2025 goMatlab.de | Dies ist keine offizielle Website der Firma The Mathworks

MATLAB, Simulink, Stateflow, Handle Graphics, Real-Time Workshop, SimBiology, SimHydraulics, SimEvents, and xPC TargetBox are registered trademarks and The MathWorks, the L-shaped membrane logo, and Embedded MATLAB are trademarks of The MathWorks, Inc.