Implizites Eulerverfahren
Nachdem ich gerade meinen Osterurlaub verlebe habe ich begonnen, ein Thema anzufassen, das mich schon länger beschäftigt - aber für Mathematik brauch ich Ruhe... Ich wollte zunächst hinter die Magie der impliziten Verfahren steigen und schauen, ob es eien Möglichkeit gibt, diese ebenfalls so generisch zu betrachten, dass man auf einfache Art und Weise ein Framework dafür erstellen kann
Ich stellte dazu zunächst die Verfahren anhand einer einfachen Differentialgleichung ersten Grades gegenüber:
Dazu benutzte ich ein kleines Python-Script:
#!/usr/bin/python
'''
Test for implicit and explicit Euler method for
numerically solving differential equations using
y'=i-my as an example and showing instability as
well as slow convergence as disadvantages
in either method
'''
import matplotlib.pyplot as plt
import numpy as np
def f(t):
return 1-np.exp(-t)
t2 = np.arange(0.0, 11.0, 0.02)
h=1.3
m=-1
I=1
x = 0
y= 0
plt.figure(1)
plt.subplot(211)
plt.plot([x],[y],'ro')
while x<10:
ynew=(y+h*I)/(1-m*h)
print " %f " % (ynew)
y=ynew
x=x+h
plt.plot([x],[y],'ro')
plt.plot(t2, f(t2), 'k')
plt.axis([0, 11, 0, 1.5])
h=1.3
m=-1
I=1
x = 0
y= 0
plt.subplot(212)
plt.plot([x],[y],'ro')
while x<10:
ynew=y+h*(I+m*y)
print " %f " % (ynew)
y=ynew
x=x+h
plt.plot([x],[y],'ro')
plt.plot(t2, f(t2), 'k')
plt.axis([0, 11, 0, 1.5])
plt.show()
Dieses Skript liefert folgendes Ergebnis (oben das implizite Verfahren, unten das explizite):
Man sieht deutlich die sehr langsame Konvergenz des impliziten Verfahrens - aber auch das Überschwingen und die Instabilität des expliziten Verfahrens bei zu groß gewählter Schrittweite...
Leider ist mir noch keine Möglichkeit eingefallen, wie man das implizite Verfahren so generisch gestalten könnte, dass man einfach die Differentialgleichung in ein Framework fallen lässt und die Berechnung startet: Dafür müsste man sich mit symbolischer Algebra und Reihenentwicklung auseinandersetzen (jedenfalls ist das meine aktuelle Meinung nach Studium des Internet).
Artikel, die hierher verlinken
Roessler Attractor und implizites Eulerverfahren
26.02.2020
Nachdem ich so viel Spaß bei der Untersuchung der Behandlung eines chaotischen Systems mittels impliziten Eulerverfahrens hatte, habe ich mir gleich noch eines vorgenommen...
Strange Attractors mit BigDecimal
26.02.2020
Nachdem ich mir in meinem früheren Artikel zum Thema der Mathematik mit beliebiger Präzision in Java noch nicht sicher war ob es Sinn machen würde, mein Framework für die numerische Behandlung von Differentialgleichungssystemen auf die Option der Benutzung von BigDecimal als unterliegenden Datentyp zu erweitern, habe ich nun doch damit experimentiert:
Mathematik mit beliebiger Präzision in Java
11.11.2018
Wann immer ich hier über Experimente mit numerischen Lösungsverfahren für Differentialgleichungssysteme berichte, habe ich im Hinterkopf, dass diese Verfahren eigentlich völlig ungeeignet dafür sind, solche Systeme zu analysieren, da heutige Digitalcomputer bereits rationale Zahlen nicht exakt darstellen können - von irrationalen ganz zu schweigen...
Chaotische Systeme und implizites Eulerverfahren
01.10.2018
Nachdem ich erfolgreich einige chaotische Systeme mittels numerischer Verfahren untersucht hatte, reifte in mir der Entschluss, für diese Systeme implizite und explizite numerische Verfahren gegenüberzustellen.
Vor 5 Jahren hier im Blog
-
Multi-User-WebDAV, Docker, GitHub
17.11.2019
Nachdem ich mich in letzter Zeit verstärkt mit Docker und dem zugehörigen Ökosystem beschäftige, habe ich begonnen, verschiedenste Dienste in Containern zu testen um zu sehen, ob in manchen Fällen LXC oder KVM nicht doch die bessere Wahl wäre...
Weiterlesen...
Tags
Android Basteln C und C++ Chaos Datenbanken Docker dWb+ ESP Wifi Garten Geo Go GUI Gui Hardware Java Jupyter Komponenten Links Linux Markdown Markup Music Numerik OpenSource PKI-X.509-CA Python QBrowser Rants Raspi Revisited Security Software-Test sQLshell TeleGrafana Verschiedenes Video Virtualisierung Windows Upcoming...
Neueste Artikel
- Migration der Webseite und aller OpenSource Projekte
In eigener Sache...
Weiterlesen... - AutoHideToolbar für Java Swing
Ich habe eine neue Java Swing Komponente erstellt: Es handelt sich um einen Wrapper für von JToolBar abgeleitete Klassen, die die Werkzeugleiste minimieren und sie nur dann einblenden, wenn der Mauszeiger über ihnen schwebt.
Weiterlesen... - Integration von EBMap4D in die sQLshell
Ich habe bereits in einem früheren Artikel über meine ersten Erfolge berichtet, der sQLshell auf Basis des bestehenden Codes aus dem Projekt EBMap4D eine bessere Integration für Geo-Daten zu spendieren und entsprechende Abfragen, bzw. deren Ergebnisse auf einer Kartenansicht zu visualisieren.
Weiterlesen...
Manche nennen es Blog, manche Web-Seite - ich schreibe hier hin und wieder über meine Erlebnisse, Rückschläge und Erleuchtungen bei meinen Hobbies.
Wer daran teilhaben und eventuell sogar davon profitieren möchte, muss damit leben, daß ich hin und wieder kleine Ausflüge in Bereiche mache, die nichts mit IT, Administration oder Softwareentwicklung zu tun haben.
Ich wünsche allen Lesern viel Spaß und hin und wieder einen kleinen AHA!-Effekt...
PS: Meine öffentlichen Codeberg-Repositories findet man hier.