Als «ode» getaggte Fragen

Gewöhnliche Differentialgleichungen (ODEs) enthalten Funktionen nur einer unabhängigen Variablen und einer oder mehrerer ihrer Ableitungen in Bezug auf diese Variable. Dieses Tag ist für Fragen zur Modellierung von Phänomenen mit ODEs, zur Lösung von ODEs und zu anderen verwandten Aspekten gedacht.

17
Gibt es einen hochwertigen nichtlinearen Programmierlöser für Python?
Ich habe mehrere herausfordernde nicht konvexe globale Optimierungsprobleme zu lösen. Derzeit verwende ich die Optimization Toolbox von MATLAB (speziell fmincon()mit algorithm = 'sqp'), was sehr effektiv ist . Der größte Teil meines Codes ist jedoch in Python, und ich würde die Optimierung gerne auch in Python durchführen. Gibt es einen …
3
Wie testet man eine Implementierung eines numerischen ODE-Lösers?
Ich beginne mit der Arbeit an einer Softwarebibliothek numerischer ODE-Löser und habe Probleme damit, Tests für die Löserimplementierungen zu formulieren. Mein Ziel ist es, dass die Bibliothek schließlich Löser für nicht steife und steife Probleme sowie mindestens einen impliziten Löser (mehr oder weniger den Fähigkeiten der odeRoutinen in Matlab gleichwertig …
26 ode 
2
Was bedeutet "symplektisch" in Bezug auf numerische Integratoren und werden sie von SciPy Odeint verwendet?
In diesem Kommentar schrieb ich: ... Standard-SciPy-Integrator, von dem ich annehme, dass er nur symplektische Methoden verwendet. in dem ich mich auf SciPy's beziehe odeint, das entweder eine "nicht steife (Adams) Methode" oder eine "steife (BDF) Methode" verwendet. Nach der Quelle : def odeint(func, y0, t, args=(), Dfun=None, col_deriv=0, full_output=0, …
4
Die Definition eines steifen ODE-Systems
Betrachten Sie ein IVP für das ODE-System , y ( x 0 ) = y 0 . Am häufigsten wird dieses Problem als steif angesehen, wenn die Jacobi-Matrix ∂ f isty′= f( x , y)y′=f(x,y)y'=f(x,y)y( x0) = y0y(x0)=y0y(x_0)=y_0hatsowohlEigenwerte mit sehr großem negativen Realteil als auch Eigenwerte mit sehr kleinem negativen …
17 ode  stiffness 
5
Warum verlässt die numerische Lösung einer ODE ein instabiles Gleichgewicht?
Ich möchte das Verhalten eines Doppelpendelsystems simulieren. Das System ist ein Roboter-Manipulator mit 2 Freiheitsgraden, der nicht betätigt wird und sich daher hauptsächlich wie ein durch die Schwerkraft beeinflusstes Doppelpendel verhält. Der einzige Hauptunterschied bei einem Doppelpendel besteht darin, dass es sich aus zwei starren Körpern zusammensetzt, deren Massen- und …
3
Numerische Methoden für diskontinuierliche rs ODEs
Was sind die neuesten Methoden zur numerischen Lösung von ODEs mit diskontinuierlicher rechter Seite? Mich interessieren vor allem stückweise glatte rechte Seitenfunktionen, zB sign. Ich versuche die Gleichung von folgendem Typ zu lösen: x˙v˙= v= { ( | Fextern| - | FReibung| )unterschreiben( Fextern)0: | Fextern| < | FReibung|: sonstx˙=vv˙={(|Fextern|-|FReibung|)Zeichen⁡(Fextern):|Fextern|<|FReibung|0:Andernfalls\begin{align*} …
15 ode 
3
Best Practices zur Beschreibung von agentenbasierten Modellen
Ich arbeite ziemlich stark in der mathematischen Biologie / Epidemiologie, wo der Großteil der Modellierungs- / Computerwissenschaftsarbeit immer noch von Sätzen von ODEs dominiert wird, zugegebenermaßen manchmal ziemlich ausgefeilten Sätzen von ihnen. Ein Pluspunkt dieser Modelle ist, dass sie sich recht einfach beschreiben und nachbilden lassen. Eine Tabelle mit Parameterwerten …
2
Optionen zum Lösen von ODE-Systemen auf GPUs?
Ich möchte das Lösen von ODE-Systemen auf GPUs in einer "trivial parallelisierbaren" Umgebung durchführen. Führen Sie beispielsweise eine Sensitivitätsanalyse mit 512 verschiedenen Parametersätzen durch. Idealerweise möchte ich ODE-Lösungen mit einem intelligenten adaptiven Zeitschrittlöser wie CVODE durchführen, anstatt mit einem festen Zeitschritt wie Forward Euler, aber auf einer NVIDIA-GPU anstelle einer …
14 ode  gpu 
3
Wahl der Schrittgröße mit ODEs in matlab
Hey, und danke, dass du dir Zeit genommen hast, meine Frage anzuschauen. Dies ist eine aktualisierte Version meiner Frage, die ich zuvor unter physics.stackexchange.com gepostet habe Ich studiere derzeit ein 2D-Excitonen-Spinor-Bose-Einstein-Kondensat und bin gespannt auf den Grundzustand dieses Systems. Die mathematische Methode zum Erreichen des Grundzustands wird als imaginäre Zeitmethode …
1
Algorithmen für das lineare System von ODEs
Ich frage mich: Was ist der beste Algorithmus zu lösen dudt= A ududt=EINu\begin{equation} \frac{du}{dt} = Au \end{equation} WobeiEINEINAeine reellen × nn×nn\times nMatrix ist. A ist nicht explizit zeitabhängig, in der Regel spärlich, aber nicht unbedingt gebändert. Ihre Eigenwerte haben nicht positive Realteile. A ist ebenfalls diagonalisierbar, kann jedoch zu groß …
Durch die Nutzung unserer Website bestätigen Sie, dass Sie unsere Cookie-Richtlinie und Datenschutzrichtlinie gelesen und verstanden haben.
Licensed under cc by-sa 3.0 with attribution required.