Löser für gewöhnliche Differentialgleichungen erster Ordnung

Löser für gewöhnliche Differentialgleichungen erster Ordnung

Der Löser für gewöhnliche Differentialgleichungen erster Ordnung nimmt eine gewöhnliche Differentialgleichung in der Form dy/dx = f(x, y) entgegen, klassifiziert automatisch deren Struktur (separierbar, linear, autonom, exakt oder allgemein) und erstellt sowohl eine symbolische geschlossene Lösung (sofern möglich) als auch eine hochpräzise numerische Lösung. Eine Live-Visualisierung des Richtungsfelds mit überlagerter Lösungskurve macht die geometrische Bedeutung der Gleichung sofort ersichtlich — Lösungen sind genau die Kurven, die tangential zu jedem Pfeil verlaufen.

Was ist eine gewöhnliche Differentialgleichung erster Ordnung?

Eine gewöhnliche Differentialgleichung erster Ordnung enthält eine unbekannte Funktion y(x) und nur deren erste Ableitung y'(x). Die explizite Standardform lautet:

Zusammen mit einer Anfangsbedingung y(x₀) = y₀ definiert dies ein Anfangswertproblem (AWP). Der Satz von Picard-Lindelöf garantiert eine eindeutige Lösung in einer Umgebung von x₀, sofern f in der Nähe von (x₀, y₀) Lipschitz-stetig bezüglich y ist. Geometrisch sucht das AWP nach der eindeutigen Kurve, die durch (x₀, y₀) verläuft und deren Steigung an jedem Punkt mit f an diesem Punkt übereinstimmt — also genau die Kurve, die das Richtungsfeld tangiert.

Sechs Klassen, die der Löser erkennt

Methode der geschlossenen Form: Linear mit konstanten Koeffizienten

Wenn sich die rechte Seite zu dy/dx = a·y + b mit den Konstanten a und b vereinfacht, liefert der integrierende Faktor μ(x) = e^(-a·x) eine exakte Lösung. Die allgemeine Lösung lautet:

Durch Anwendung der Anfangsbeding y(x₀) = y₀ wird die Konstante C festgelegt und man erhält die eindeutige spezielle Lösung. Diese eine Klasse deckt eine enorme Anzahl von Lehrbuchproblemen ab:

• Exponentielles Wachstum — dy/dx = k·y, spezielle Lösung y(t) = y₀·e^(k·t).
• Exponentieller Zerfall — dy/dx = -k·y, Halbwertszeit ln 2 / k.
• Newtonsches Abkühlungsgesetz — dy/dx = -k·(y - T_amb), die Körpertemperatur nähert sich exponentiell der Umgebungstemperatur an.
• Laden eines RC-Glieds — dV/dt = (1/RC)·(V_in - V), die Kondensatorspannung nähert sich der Quelle an.
• Arzneimittel-Clearance — Pharmakokinetik erster Ordnung mit Eliminationsrate k.

Ein Richtungsfeld lesen

An jedem Gitterpunkt (x, y) zeichnet das Tool ein kurzes Liniensegment, dessen Steigung f(x, y) entspricht. Drei nützliche Beobachtungen:

• Gleichgewichtspunkte sind Stellen, an denen f(x, y) = 0 gilt — das Richtungsfeld ist dort horizontal. Bei autonomen Gleichungen sind dies Fixpunkte y*, die f(y*) = 0 erfüllen; nahegelegene Trajektorien nähern sich y* entweder an (stabil) oder entfernen sich davon (instabil).
• Isoklinen sind Kurven, auf denen f(x, y) gleich einer Konstanten c ist, sodass alle Pfeile entlang der Kurve die gleiche Steigung c haben.
• Lösungskurven kreuzen sich nie (wenn f Lipschitz-stetig ist) — visuell offensichtlich, da zwei sich kreuzende Kurven am Schnittpunkt unterschiedliche Steigungen haben müssten.

Numerisches Verfahren: Klassisches Runge-Kutta (RK4)

Ausgehend von (x_n, y_n) wird der nächste Wert durch Mittelung von vier Steigungsschätzungen berechnet:

RK4 hat einen lokalen Abbruchfehler von O(h⁵) und einen globalen Fehler von O(h⁴), was bei der Standard-Schrittzahl für nicht-steife Gleichungen eine Genauigkeit von etwa sechs Stellen liefert. Der Löser integriert vom Anfangspunkt aus in beide x-Richtungen und stoppt sauber, wenn der Betrag von y den Wert 10¹⁵ überschreitet — typisch für Lösungen, die in endlicher Zeit explodieren, wie dy/dx = y².

So verwenden Sie diesen Rechner

1. Geben Sie die rechte Seite in das Feld dy/dx = ... ein. Verwenden Sie x und y als Variablen, * für Multiplikation, ^ oder ** für Potenzen und Standardfunktionen wie sin, cos, exp, log, sqrt. Die Konstanten pi und e werden erkannt.
2. Geben Sie die Anfangsbedingung (x₀, y₀) an — die eindeutige Lösungskurve wird durch diesen Punkt verlaufen.
3. Wählen Sie den x-Bereich, über den das Richtungsfeld und die Lösungskurve geplottet werden sollen. Der y-Bereich wird automatisch an die integrierte Lösung angepasst.
4. Klicken Sie auf Lösen & Visualisieren. Zuerst läuft die Klassifizierung; wenn Ihre Gleichung einem Muster für eine geschlossene Form entspricht (linear mit konstanten Koeffizienten), erhalten Sie die symbolische Antwort. Das Richtungsfeld und die Lösungskurve werden immer gerendert.
5. Schalten Sie das Richtungsfeld ein oder aus, um sich auf die Lösungskurve zu konzentrieren, oder wiederholen Sie die Animation der Kurvenzeichnung, um zu sehen, wie die Integration vom Anfangspunkt aus fortschreitet.

Anwendungsbeispiel: Newtonsches Abkühlungsgesetz

Eine Tasse Kaffee mit 80 °C kühlt in einem 20 °C warmen Raum ab. Die Wärmeübertragungsrate ist proportional zur Temperaturdifferenz:

Dies ist linear mit konstanten Koeffizienten (a = -0.1, b = 2). Die geschlossene Form lautet:

Nach 30 Minuten: T(30) = 20 + 60·e⁻³ ≈ 22.99 °C. Die Richtungsfeld-Ansicht macht das Grenzverhalten deutlich — jede Lösungskurve nähert sich, unabhängig von der Starttemperatur, asymptotisch der horizontalen Linie T = 20 an.

Häufige Anwendungen

• Populationsdynamik — Modelle für exponentielles, logistisches Wachstum oder den Allee-Effekt.
• Pharmakokinetik — Wirkstoffaufnahme und -elimination, Halbwertszeitberechnungen.
• Wärmeübertragung — Newtonsches Abkühlungsgesetz, Modelle mit konzentrierter Kapazität.
• RC- und RL-Schaltkreise — Lineare elektrische Einschwingvorgänge erster Ordnung.
• Radioaktiver Zerfall — Zerfallsketten einzelner Isotope.
• Mischtanks — Konzentration eines gelösten Stoffes bei Zu- und Abfluss.
• Fallendes Objekt mit Luftwiderstand — Analyse der Endgeschwindigkeit dv/dt = g - kv.

Häufig gestellte Fragen

Was ist eine gewöhnliche Differentialgleichung erster Ordnung?

Eine gewöhnliche Differentialgleichung (ODE) erster Ordnung ist eine Gleichung der Form dy/dx = f(x, y), die die unbekannte Funktion y(x) und ihre erste Ableitung enthält. Das Lösen der Differentialgleichung bedeutet, die Funktion y(x) zu finden, deren Ableitung mit der rechten Seite übereinstimmt. Mit einer Anfangsbedingung y(x₀) = y₀ ist die Lösung unter milden Regularitätsbedingungen (Satz von Picard-Lindelöf) eindeutig.

Was ist ein Richtungsfeld?

Ein Richtungsfeld (oder Vektorfeld) zeichnet an jedem Gitterpunkt (x, y) ein kleines Liniensegment, dessen Steigung f(x, y) entspricht. Lösungskurven der Differentialgleichung sind genau die Kurven, die an jedem Punkt tangential zu diesen Segmenten verlaufen. Das Richtungsfeld liefert eine sofortige visuelle Intuition für das globale Verhalten der Lösungen, ohne die Gleichung symbolisch lösen zu müssen.

Welche Klassen von Differentialgleichungen erster Ordnung löst dieses Tool?

Das Tool klassifiziert die Gleichung automatisch in eine der folgenden Kategorien: integrierbar (hängt nur von x ab, gelöst durch direkte Integration), linear mit konstanten Koeffizienten y' = a·y + b (vollständige geschlossene Form wird bereitgestellt), autonom (hängt nur von y ab), separierbar (faktorisierbar als g(x)·h(y)), linear mit variablen Koeffizienten (P(x)·y + Q(x)) oder allgemein. Für jede Klasse werden eine hochpräzise Runge-Kutta-Numeriklösung und eine Richtungsfeld-Visualisierung erstellt.

Welches numerische Verfahren wird verwendet?

Das klassische Runge-Kutta-Verfahren vierter Ordnung (RK4) wird mit 300 Teilschritten in jede Richtung vom Anfangspunkt aus angewendet. RK4 hat einen lokalen Abbruchfehler von O(h⁵) und ist das Standardwerkzeug für nicht-steife Differentialgleichungen in dieser Größenordnung. Der Löser erkennt Divergenz (Überlauf oder NaN) und stoppt die Integration sauber, damit der Plot gültig bleibt.

Was ist die Methode des integrierenden Faktors für lineare Differentialgleichungen?

Für eine lineare Differentialgleichung erster Ordnung y' + P(x)·y = Q(x) multipliziert man beide Seiten mit dem integrierenden Faktor μ(x) = e^∫P(x) dx. Die linke Seite wird zur exakten Ableitung d/dx[μ·y], sodass y(x) = (1/μ(x)) · (∫μ(x)·Q(x) dx + C) gilt. Wenn P und Q Konstanten sind, vereinfacht sich dies zur geschlossenen Form y = -b/a + C·e^(a·x), die das Tool automatisch zurückgibt.

Kann dieses Tool steife Gleichungen oder Systeme von Differentialgleichungen verarbeiten?

Dieser Löser ist für nicht-steife skalare Differentialgleichungen erster Ordnung gedacht. Sehr steife Probleme (bei denen die Lösung mehrere Zeitskalen hat, die sich um viele Größenordnungen unterscheiden) erfordern möglicherweise ein implizites Verfahren wie den impliziten Euler oder Rosenbrock; gekoppelte Systeme benötigen einen vektorwertigen Löser. Verwenden Sie für diese Fälle ein spezielles Paket wie solve_ivp von SciPy oder spezialisierte Löser für steife Differentialgleichungen.

Weiterführende Literatur

• Gewöhnliche Differentialgleichung — Wikipedia
• Richtungsfeld — Wikipedia
• Runge-Kutta-Verfahren — Wikipedia
• Integrierender Faktor — Wikipedia
• Satz von Picard-Lindelöf — Wikipedia