Gewöhnliche Differentialgleichung

Eine gewöhnliche Differential- bzw. Differenzialgleichung ist eine Differentialgleichung, die nur Ableitungen nach einer reellen Variablen enthält. Ihre Lösung ist somit eine Funktion, die von einer Variablen abhängt.

Inhaltsverzeichnis

Motivation

Differentialgleichungen werden oft benötigt, um Vorgänge zu beschreiben, bei denen die Veränderung einer Größe durch sie selbst bestimmt wird.

Die historisch ersten Differentialgleichungen waren die der gleichmäßigen und ungleichmäßig beschleunigten Bewegung, welche von Galileo Galilei noch mit geometrischen Methoden bearbeitet werden konnten. Als Isaac Newton jedoch auch Bewegungen unter zum Betrag oder Quadrat der Geschwindigkeit proportionaler Reibung betrachtete, war er genötigt, den Differentialkalkül und die heute geläufige Form einer Differentialgleichung einzuführen.

Das Zerfallsgesetz in der Physik etwa besagt, dass die Anzahl der pro Zeiteinheit zerfallenden Atome einer Menge instabiler Atome von der gesamten Anzahl N der vorhandenen Atome abhängt. Insofern ist die Abnahme der Anzahl der Atome proportional zur Anzahl aller Atome:

\frac{\mathrm{d}}{\mathrm{d}t}N(t) = c\; N(t).

Durch Berechnen der Funktion N(t)\! aus dieser Differentialgleichung kann die Anzahl der Atome zu jedem Zeitpunkt bestimmt werden.

Ein anderes einfaches Beispiel ist der ungedämpfte harmonische Oszillator mit der Differentialgleichung

m\;a = m\; \frac{\mathrm{d}^2}{\mathrm{d}t^2}x(t) = -k \; x(t).

Die gesuchte Funktion ist hier die Funktion x(t)\!, deren zweite zeitliche Ableitung als Beschleunigung aus den Bewegungsgesetzen stammt.

Ordnung einer Differentialgleichung und Ordnungsreduktion

Die Ordnung einer Differentialgleichung ist durch die Ordnung der höchsten vorkommenden Ableitung gegeben. Gewöhnliche Differentialgleichungen beliebiger Ordnung lassen sich immer in ein System von gewöhnlichen Differentialgleichungen erster Ordnung umwandeln:

Sei

y^{(n)}(t)=f(t,y(t),y'(t),y''(t),\dots,y^{(n-1)}(t))

eine gewöhnliche Differentialgleichung n-ter Ordnung. Dann werden folgende Hilfsfunktionen eingeführt:

u_k(t)=y^{(k)}(t) \quad k \in {0 \dots n-1}

Wir erhalten dadurch ein System von n gewöhnlichen Differentialgleichungen 1.Ordnung:

u'_0(t)=u_1(t)\!

u'_1(t)=u''_0(t)=u_2(t)\!

\dots

u'_{n-1}=u''_{n-2}= \dots = u^{(n)}_0=y^{(n)}=f(t,y,u_0,u_1,\dots\,u_{n-1}).

Umgekehrt kann man natürlich auch aus einem Differentialgleichungssystem eine einzige Differentialgleichung höherer Ordnung ableiten.

Lineare Differentialgleichungen

Die Theorie der linearen Differentialgleichungen ist ein wichtiger Teilbereich der gewöhnlichen Differentialgleichungen. Lineare Differentialgleichungen lassen sich oft mit Standardmethoden lösen.

Eine lineare Differentialgleichung ist linear in der unbekannten Funktion und ihren Ableitungen. Die allgemeine Form für eine gewöhnliche lineare Differentialgleichung für die Funktion y(t) ist

\sum_{i=0}^n c_i(t) y^{(i)}(t) = s(t)

c_i(t)\! und s(t)\! sind hierbei bekannte Funktionen, y(t)\! wird gesucht, y^{(i)}(t)\! ist die i-te Ableitung von y nach t. Man unterscheidet lineare Differentialgleichungen mit variablen oder konstanten (von t unabhängigen) Koeffizienten c_i(t)\! und hat in jedem dieser beiden Fälle homogene (mit s(t) = 0\!) und inhomogene (mit s(t) \not= 0\!) Problemstellungen.

Das Lösen von Differentialgleichungen

Die Lösung einer Differentialgleichung ist immer eine Funktion (oder im Falle eines Systems von Differentialgleichungen mehrere Funktionen). Es ist jedoch nicht jede Differentialgleichung lösbar, es gibt allerdings einige Kriterien, anhand derer man Lösbarkeit erkennen kann. Ferner reicht die Differentialgleichung allein im Allgemeinen nicht aus, um die Funktion eindeutig zu bestimmen. Die allgemeine Lösung einer gewöhnlichen Differentialgleichung n-ter Ordnung hat im Allgemeinen n freie Parameter. Die allgemeine Lösung einer partiellen Differentialgleichungen von n Unbekannten enthält im Allgemeinen eine frei wählbare (aber hinreichend oft differenzierbare) Funktion von n-1 Variablen, die selbst Funktionen der n Unbekannten sind (diese Funktionen sind aber natürlich nicht frei wählbar sondern werden durch die Lösung bestimmt).

Beispielsweise werden alle schwingenden Pendel durch eine Differentialgleichung beschrieben, und der generelle Bewegungsablauf folgt immer dem gleichen Prinzip. Der konkrete Bewegungsablauf ist jedoch durch die Rand- oder Anfangsbedingung(en) (wann wurde das Pendel angestoßen, und wie weit) bestimmt. Die Lösbarkeit von Anfangswertproblemen bei gewöhnlichen Differentialgleichungen 1. Ordnung wird durch den Satz von Picard-Lindelöf beschrieben. Nur wenige Typen von Differentialgleichungen lassen sich analytisch lösen. Insbesondere partielle Differentialgleichungen können oft nur mit numerischen Methoden approximiert werden.

Die Menge aller Lösungen einer gewöhnlichen Differentialgleichung bildet ein dynamisches System, auch Fluss der Differentialgleichung genannt.

Lineare Differentialgleichungen

Für Lösungen homogener linearer Differentialgleichung gilt das Superpositionsprinzip: Eine Linearkombination mehrerer Lösungen ist wieder eine Lösung. n unabhängige Lösungen einer Differentialgleichung n-ter Ordnung bilden ein Fundamentalsystem der Differentialgleichung. Es gibt verschiedene (unendlich viele) Fundamentalsysteme für eine gegebene Gleichung, die Anzahl der Funktionen des Fundamentalsystems ist aber immer gleich.

Bei inhomogenen Differentialgleichungen löst man meist zuerst die zugehörige homogene Differentialgleichung. Die allgemeine Lösung der inhomogenen Gleichung ist dann

y(t) = y_h(t) + y_p(t)\!

wobei y_h(t)\! die allgemeine Lösung der homogenen Gleichung und y_p(t)\! eine beliebige Lösung der inhomogenen Gleichung ist: Addiert man zu einer Lösung der inhomogenen Gleichung eine Lösung der homogenen Gleichung dazu, ist das Ergebnis wieder eine Lösung der inhomogenen Gleichung.

Lineare homogene Differentialgleichungen mit konstanten Koeffizienten

Hier führt der Ansatz y(t) = e^{\lambda \cdot t} mit zunächst unbekanntem λ zum Ziel. Dadurch erhält man ein Polynom, dessen Grad gleich der Ordnung der Differentialgleichung ist:

\sum_{i=0}^n c_i \cdot \lambda^i = 0

Dieses Polynom hat im Allgemeinen n komplexe Lösungen, daraus erhält man n Funktionen der Form e^{\lambda_j \cdot t}, die ein Fundamentalsystem bilden. Bei Mehrfachlösungen sind auch die zusätzlichen Funktionen t^k \cdot e^{\lambda_j \cdot t} mit natürlichen Zahlen k zwischen (ausschließlich) 0 und der Vielfachheit der Lösung des Polynoms Lösungen der Differentialgleichung. Sind alle Koeffizienten ci reell, so erhält man ein rein reelles Fundamentalsystem indem man bei den Funktionenpaaren mit nichtreellen, komplex-konjugierten λs e^{\lambda_j \cdot t} einmal durch e^{Re(\lambda_j) \cdot t} \cdot cos(Im(\lambda_j) \cdot t) und einmal durch e^{Re(\lambda_j) \cdot t} \cdot sin(Im(\lambda_j) \cdot t) ersetzt.

Spezielle Lösungsmethoden

Neben linearen Systemen lassen sich Differentialgleichungen, die separierbar sind, durch direkte Integration lösen.

Manche Typen von Differentialgleichungen lassen sich durch Potenzreihen lösen.

Das Lösen von DGL mit getrennten Variablen

Ist eine DGL der Form \dot c(t) = g( c(t) ) f(t)\! gegeben, kann diese mit getrennten Variablen berechnet werden. Als Startwert kann c(t_0) = c_0 \! angenommen werden.

Beispiel

\dot c(t) = 2t c(t) - t = t (2c(t) - 1)

Dann setzt man g( c(t) ) = 2c(t) - 1 \! und f(t) = t \!

Mit dem obigen Verfahren ergibt sich dann: F(t)= \int_{t_o}^{t} s ds = { \frac{s^2}{2} } \mid_{t_0}^{t} = {t^2 \over 2} - \frac{t_0^2}{2}

H(v) = \int_{c_o}^{v} \frac{1}{ 2c(t) - 1 } dx = \ln (2 c(v) - 1) \mid_{c_0}^{v} = \ln (2 c(v) - 1) - \ln (2 c(c_0) - 1)

F(t) = H( c(t) ) \Rightarrow \frac{t^2}{2} - \frac{t_0^2}{2} = \ln (2 c(t) - 1) - \ln (2 c(c_0) - 1)

e^{ \frac{t^2 - {t_0}^2}{2} } = e^{ \ln (2 c(t) - 1) - \ln (2 c(c_0) - 1) }

e^{ \frac{t^2 - {t_0}^2}{2} } = \frac{ e^{\ln (2 c(t) - 1)} }{e^{\ln (2 c(c_0) - 1)} }

e^{ \frac{t^2 - {t_0}^2}{2} } = \frac{2 c(t) - 1}{2 c(c_0) - 1}

c(t) = \frac{e^ { \frac{t^2  - t_0^2}{2} } (2c(c_0) - 1) + 1}{2}

Laplace-Transformation zur Lösung von Differentialgleichungen

Die Laplace-Transformation eignet sich aufgrund ihres Differentiationssatzes u.a. dazu Differentialgleichungen zu lösen. Dazu transformiert man die Differentialgleichung in den Spektralbereich, löst die so erhaltene algebraische Gleichung und transformiert die Lösung in den Zeitbereich zurück.

Besonderen Wert hat hier die Laplace-Transformation bei Differentialgleichungssystemen: Da Ableitungen im Bildbereich als Produkt aus Originalfunktion und Laplace-Faktor s entstehen, werden DGL Systeme zu einfacheren normalen Gleichungssystemen.

Spezielle Differentialgleichungen

d'Alembert-Differentialgleichung
Bernoulli-Gleichung
Clairaut-Gleichung
Exponentialfunktion y'=cy \Rightarrow y=ce^{\alpha x}
Eulersche Differentialgleichung
Riccati-Gleichung

Numerische Verfahren für gewöhnliche Differentialgleichungen

Da sich gewöhnliche Differentialgleichungen höherer Ordnung immer auf Systeme erster Ordnung reduzieren lassen, geht man bei der Konstruktion von numerischen Lösungsverfahren im Normalfall von einem System erster Ordnung aus:

\dot{x}(t)=f(t,x(t)), \; x:\mathbb{R}^{n} \mapsto \mathbb{R}^{n}, f:\mathbb{R}^{n+1} \mapsto \mathbb{R}^{n}

Es gibt zwei wichtige Klassen von numerischen Verfahren zur Lösung gewöhnlicher Differentialgleichungen , die Einschrittverfahren (insbesondere die Runge-Kutta-Verfahren) und die linearen Mehrschrittverfahren. Eine Verallgemeinerung von beiden Klassen stellen die allgemeinen linearen Verfahren (General linear Methods (GLM)) dar.

Siehe auch

  1. partielle Differentialgleichung
  2. Integralgleichung
  3. Anfangswertproblem
  4. Randwertproblem

See also: Gewöhnliche Differentialgleichung, Anfangsbedingung, Anfangswertproblem, Atom, Bernoulli-Differentialgleichung, Beschleunigung, Bewegungsgesetze, Clairaut-Gleichung, D'Alembert-Differentialgleichung, Differentialgleichung