Gewöhnliche Differentialgleichung

Eine gewöhnliche Differentialgleichung (oft abgekürzt mit GDGL oder ODE, englisch ordinary differential equation) ist eine Differentialgleichung, bei der zu einer gesuchten Funktion nur Ableitungen nach genau einer Variablen auftreten.

Viele physikalische, chemische und biologische Vorgänge in der Natur lassen sich mit solchen Gleichungen mathematisch beschreiben, z.B. der radioaktive Zerfall, Bewegungsvorgänge von Körpern, viele Arten von Schwingungsvorgängen oder das Wachstumsverhalten von Tier-Populationen. In naturwissenschaftlichen Modellen werden gewöhnliche Differentialgleichungen daher häufig eingesetzt, um solche Vorgänge zu analysieren, zu simulieren oder um Vorhersagen abgeben zu können. In vielen Fällen kann die Differentialgleichung nicht analytisch gelöst werden. Man ist daher auf numerische Verfahren angewiesen.

Historische Entwicklung

Historisch gesehen wurden die ersten Differentialgleichungen verwendet, um die Bewegung von Objekten zu modellieren. Besonders hervorzuheben sind dabei die Gleichungen für die Bewegung mit konstanter Geschwindigkeit bzw. konstanter Beschleunigung. Im Jahr 1590 erkannte Galileo Galilei den Zusammenhang zwischen der Fallzeit eines Körpers und seiner Fallgeschwindigkeit sowie dem Fallweg und formulierte (noch) mit geometrischen Mitteln das Gesetz des freien Falles.

Als Isaac Newton auch Bewegungen mit Reibungen betrachtete, die zum Betrag oder zum Quadrat der Geschwindigkeit proportional sind, war er genötigt, die Differentialrechnung und den heute geläufigen Formalismus der Differentialgleichungen einzuführen.

Durch die exakte Formulierung des Grenzwertbegriffes, der Ableitung und des Integrals stellte schließlich Augustin Louis Cauchy im 19. Jahrhundert die Theorie der gewöhnlichen Differentialgleichungen auf ein festes Fundament und machte sie somit vielen Wissenschaften zugänglich.

Das wissenschaftliche Interesse an Differentialgleichungen ist im Wesentlichen darin begründet, dass mit ihnen auf Grund vergleichsweise einfacher Beobachtungen und Experimente vollständige Modelle geschaffen werden können.

Nur wenige Typen von Differentialgleichungen lassen sich analytisch lösen. Trotzdem lassen sich qualitative Aussagen wie Stabilität, Periodizität oder Bifurkation auch dann treffen, wenn die Differentialgleichung nicht explizit gelöst werden kann. Eines der wichtigsten Hilfsmittel für skalare Differentialgleichungen sind Argumente mittels eines Vergleichssatzes.

Allgemeine Definition

Seien $\Omega \subseteq \mathbb {R} \times \left(\mathbb {R} ^{m}\right)^{n+1},n\in \mathbb {N}$ und $f\colon \Omega \to \mathbb {R} ^{m}$ eine stetige Funktion. Dann heißt

$f\left(x,y,y',y'',\dotsc,y^{(n)}\right)=0$

ein gewöhnliches Differentialgleichungssystem -ter Ordnung von Gleichungen ( ist hier die unabhängige Variable, $y'={\tfrac {dy}{dx}}$ usw.). Im Fall m=1 nennt man dies eine gewöhnliche Differentialgleichung -ter Ordnung.

Ihre Lösungen sind -mal differenzierbare Funktionen $y\colon I\to \mathbb {R} ^{m}$ , welche die Differentialgleichung auf einem zu bestimmenden Intervall $I\subset \mathbb{R}$ erfüllen. Sucht man eine spezielle Lösung, welche zu gegebenen $x_0 \in I$ und $y_{0},\dotsc ,y_{n-1}\in \mathbb {R} ^{m}$ zusätzlich

$y(x_0) = y_0,\;y'(x_0) = y_1, \dotsc, y^{(n-1)}(x_0) = y_{n-1}$

erfüllt, so bezeichnet man dies als Anfangswertproblem.

Kann die Differentialgleichung nach der höchsten vorkommenden Ableitung aufgelöst werden und hat somit die Form

$y^{(n)}=f\left(x,y,y',y'',\dotsc ,y^{(n-1)}\right)$ ,

so heißt sie explizit, andernfalls implizit; siehe auch Satz von der impliziten Funktion.

Zur Notation

In der Literatur zu gewöhnlichen Differentialgleichungen werden standardmäßig zwei unterschiedliche Notationen verwendet. In der einen Variante wird die unabhängige Variable mit bezeichnet und die Ableitungen der Funktion nach mit $y',y''$ usw. Die andere Schule verwendet eine auf Newton zurückgehende Notation. Dabei ist die unabhängige Variable bereits mit einem Sinn versehen; ist die Zeit. Lösungen werden dann oft mit bezeichnet und die Ableitungen nach der Zeit werden als ${\dot {x}},{\ddot {x}}$ notiert. Da dieser Artikel von Vertretern beider Schulen bearbeitet wurde, finden sich beide Notationen wieder.

Existenz und Eindeutigkeit

Ob überhaupt eine Lösung existiert, lässt sich anhand einiger Kriterien erkennen. Die Differentialgleichung selbst reicht im Allgemeinen nicht aus, um die Lösung eindeutig zu bestimmen.

Beispielsweise ist der grundsätzliche Bewegungsablauf aller schwingenden Pendel gleich und kann durch eine einzige Differentialgleichung beschrieben werden. Der konkrete Bewegungsablauf ist jedoch durch die Rand- oder Anfangsbedingung(en) (wann wurde das Pendel angestoßen, und wie groß ist die Anfangsauslenkung) bestimmt.

Die lokale Lösbarkeit von Anfangswertproblemen bei gewöhnlichen Differentialgleichungen erster Ordnung wird durch den Satz von Picard-Lindelöf und den Satz von Peano beschrieben. Aus der Existenz einer lokalen Lösung kann man in einem zweiten Schritt auf die Existenz einer nicht-fortsetzbaren Lösung schließen. Mit Hilfe des Satzes vom maximalen Existenzintervall kann man darauf aufbauend von dieser nicht-fortsetzbaren Lösung dann gelegentlich Globalität nachweisen. Die Eindeutigkeit bekommt man als Anwendung der gronwallschen Ungleichung.

Reduktion von Gleichungen höherer Ordnung auf Systeme erster Ordnung

Gewöhnliche Differentialgleichungen beliebiger Ordnung lassen sich immer auf ein System von gewöhnlichen Differentialgleichungen erster Ordnung zurückführen. Hat eine gewöhnliche Differentialgleichung die Ordnung , so führt man dazu die voneinander abhängigen Funktionen $y_{1},y_{2},\dotsc ,y_{n}$ ein:

${\begin{aligned}y_{1}&:=y\\y_{2}&:=y'_{1}\\y_{3}&:=y'_{2}\\&\vdots \\y_{n}&:=y'_{n-1}\\\end{aligned}}$

Aus der expliziten Differentialgleichung -ter Ordnung für wird dabei:

$y'_{n}=f(x,y_{1},\ldots ,y_{n})$

Man erhält also ein System von gewöhnlichen Differentialgleichungen erster Ordnung:

$(y_{1},\dotsc ,y_{n-1},y_{n})'=(y_{2},\dotsc ,y_{n},f(x,y_{1},\ldots ,y_{n}))$ .

Umgekehrt kann man aus manchen Differentialgleichungssystemen eine einzige Differentialgleichung höherer Ordnung ableiten.

Beispiele

Ein einfaches Beispiel aus der Physik ist das Zerfallsgesetz:

$\dot N \sim N$

Dieses besagt, dass bei einer Menge instabiler Atome die Anzahl der zerfallenden Atome von der gesamten Anzahl der vorhandenen Atome proportional abhängt.

Eine wichtige Klasse weiterer Differentialgleichungen bilden die newtonschen Bewegungsgleichungen:

$m \cdot \ddot{\vec{r}}(t) = \vec{F} \left(\vec{r}(t),t\right)$

Durch die Kenntnis der von der Zeit

und der Position eines Teilchens abhängenden Kraft treffen diese Gleichungen Aussagen über die Bewegung des Teilchens selbst.

Neben einfachen Zusammenhängen der Änderungen einer einzelnen Größe lassen sich aber auch Vorhersagen über mehrere Größen in einem System treffen. In etwa die Lotka-Volterra-Gleichungen der Ökologie:

$\dot r = Z_r r b - M_r r$

$\dot b = Z_b b - M_b r b$

Dieses System beschreibt die zeitliche Veränderung der Räuberpopulation und der Beutepopulation

bei konstanten natürlichen Geburtenraten und Sterberaten . Einige wichtige Eigenschaften dieses Modells lassen sich in Form der sogenannten Lotka-Volterra-Regeln zusammenfassen. Dieses und ähnliche Systeme finden in der theoretischen Biologie auch zur Beschreibung von Ausbreitungsprozessen und in Epidemiemodellen breite Anwendung.

Spezielle Typen von Differentialgleichungen

Den bekanntesten Typ der gewöhnlichen Differentialgleichungen bildet die lineare Differentialgleichung -ter Ordnung mit:

$\sum_{i=0}^n a_i(x) y^{(i)}(x) = b(x)$ für stetige $a_i\colon \mathbb{R} \rightarrow \mathbb{R}$ .

Weitere wichtige Typen von gewöhnlichen Differentialgleichungen sind die folgenden:

d’Alembertsche Differentialgleichung

$y(x) = x g(y'(x)) + f(y'(x)) \,$ .

Bernoullische Differentialgleichung

$y'(x)=f(x)y(x) + g(x)y^\alpha(x) \,$ mit $\alpha\neq 1$ .

Exakte Differentialgleichung

$p\left(x,y(x)\right)+q(x,y(x))y'(x)=0 \,$ , worin das Vektorfeld (p,q)

eine Potentialfunktion besitzt.

Jacobische Differentialgleichung

$y'(x)= f\left(\frac{ax + by(x) + c}{\alpha x + \beta y(x) + \gamma}\right)$ .

Lineares Differentialgleichungssystem erster Ordnung von Gleichungen

$\ y'(x) = A(x)y(x) + b(x)$ für stetige $A\colon \mathbb{R} \rightarrow \mathbb{R}^{m \times m}$ und $b\colon \mathbb{R} \rightarrow \mathbb{R}^m \,$ .

Riccatische Differentialgleichung

$y'(x) = f(x)y^2(x)+ g(x)y(x) + h(x) \,$ .

Separierbare Differentialgleichung

$y'(x) = f\left(y(x)\right)g(x) \,$ .

Autonome Systeme

Ein Differentialgleichungssystem heißt autonom oder zeitinvariant, falls die beschreibende Gleichung nicht von der unabhängigen Variable abhängt. D.h., wenn das System von der Form

$f\left(y,y',y'',\dotsc,y^{(n)}\right)=0$

ist.

Ein Differentialgleichungssystem

$y'=f(y)\ ,\ y(0)=y_{0}\in \mathbb {R} ^{m}$

heißt vollständig, wenn zu jedem Anfangswert $y_{0}$ die globale Lösung auf ganz $\mathbb {R}$ definiert und eindeutig ist. Dies ist z.B. der Fall, wenn $f\colon\mathbb{R}^m\rightarrow\mathbb{R}^m$ linear beschränkt und Lipschitz-stetig ist. Es bezeichne $\varphi(\cdot,y_0)$ diese (eindeutig bestimmte globale) Lösung. Dann nennt man $\varphi :\mathbb {R} \times \mathbb {R} ^{m}\to \mathbb {R} ^{m}$ den Fluss der Differentialgleichung y'=f(y) , und $(\mathbb{R}, \mathbb{R}^m, \varphi)$ bildet dann ein dynamisches System.

Besonders einfach zu analysieren ist der Fall $n=1,\ m=2$ der ebenen autonomen Systeme. Mit Hilfe des Satzes von Poincaré-Bendixson kann man oft die Existenz periodischer Lösungen nachweisen. Ein wichtiges ebenes autonomes System bildet das Lotka-Volterra-Modell.

Da die Poincaré-Bendixson-Theorie zentral auf den jordanschen Kurvensatz aufbaut, sind höherdimensionale Analoga falsch. Insbesondere ist es sehr schwierig, periodische Lösungen höherdimensionaler autonomer Systeme zu finden.

Lösungsverfahren für lineare gewöhnliche Differentialgleichungen mit konstanten Koeffizienten

Durch gewöhnliche Differentialgleichungen lassen sich viele dynamische Systeme aus der Technik, Natur und Gesellschaft beschreiben. Viele auf den ersten Blick sehr verschiedene physikalische Probleme lassen sich mit der GDGL jedoch formal identisch darstellen.

Ein dynamisches System ist eine Funktionseinheit zur Verarbeitung und Übertragung von Signalen, wobei die Eingangsgröße u(t) als Ursache und die Ausgangsgröße y(t) als Folge des zeitlichen Übertragungsverhaltens des Systems definiert ist. Ist die Eingangsgröße u(t)=0 , so handelt es sich um eine homogene GDGL, anderenfalls um eine inhomogene GDGL.

Ein dynamisches System verhält sich linear, wenn die Wirkungen zweier linear überlagerter Eingangssignale sich am Ausgang des Systems in gleicher Weise linear überlagern. Eine lineare GDGL enthält die gesuchte Funktion und deren Ableitungen nur in der ersten Potenz. Es dürfen keine Produkte der gesuchten Funktion und ihren Ableitungen auftreten. Die gesuchte Funktion darf auch nicht in Argumenten von Winkelfunktionen, Logarithmen usw. erscheinen.

Ein bekanntes Beispiel aus der Mechanik ist die lineare GDGL zweiter Ordnung eines gedämpften Federpendels mit der Federsteifigkeit , Masse und Dämpfungskonstante . Dabei ist die Eingangsgröße: die Kraft , die Ausgangsgröße der Weg .

${\ddot {x}}(t)+{\frac {d}{m}}\cdot {\dot {x}}(t)+{\frac {c}{m}}\cdot x(t)={\frac {1}{m}}\cdot F(t)$ .

Linear zeitinvariante Systeme können durch die nachfolgenden Verfahren berechnet werden:

Klassisch mit Hilfe des Exponentialansatzes,
Laplace-Transformation,
Numerisch.

Lösung mit Hilfe des Exponentialansatzes

Die Lösung einer inhomogenen GDGL besteht aus der allgemeinen Lösung der homogenen GDGL und einer speziellen Lösung (partikuläre Lösung) der inhomogenen GDGL. Deshalb erfolgt das Lösungsverfahren der inhomogenen GDGL, unabhängig von der Ordnung, in zwei Stufen. Die Gesamtlösung ist die Summe der beiden Lösungen:

Die homogene Lösung der GDGL beschreibt das Systemverhalten mit Anfangswerten der Systemspeicher zum Zeitpunkt und dem Eingangssignal . Dies bedeutet für das dynamische System, es ist sich selbst überlassen und hat nur ein Ausgangssignal. Die homogene Lösung der GDGL ist Null, wenn alle Anfangsbedingungen von $y_{0}=0$ und deren Ableitungen Null sind.

Die partikuläre Lösung der GDGL beschreibt das Übertragungsverhalten von für $u(t)\neq 0$ als erzwungene Bewegung. Je nach Systemordnung müssen alle Anfangsbedingungen $y_{0}=0$ und deren Ableitungen Null sein.

Ist die Übertragungsfunktion G(s)

als Laplace-transformierte GDGL gegeben, so ist die Berechnung des System-Ausgangssignals Y(s)

für ein gegebenes Eingangssignal U(s)

bei Anwendung der inversen Laplace-Transformation immer eine partikuläre Lösung. Die partikuläre Lösung der GDGL ist in der Regelungstechnik meist von hauptsächlichem Interesse.

Mit Hilfe des Exponentialansatzes und der sich daraus ergebenden charakteristischen Gleichung lassen sich auch GDGL höherer Ordnung lösen. Dieser Exponentialansatz gilt als universelles Lösungsverfahren für homogene GDGL beliebiger Ordnungen mit konstanten Koeffizienten.

Hat eine GDGL die Ordnung n, so hat ihre Lösung n Integrationskonstanten. Dazu müssen n Anfangsbedingungen gegeben sein.

Der Exponentialansatz $y(t)=e^{\lambda \cdot t}$ liefert Ableitungen der Form: $y^{(n)}(t)=\lambda ^{n}\cdot e^{\lambda \cdot t}$ .

Werden diese Beziehungen in die homogene GDGL eingesetzt, entsteht die charakteristische Gleichung als Polynom n-ter Ordnung für $\lambda$ . Die homogene Lösung einer inhomogenen Differenzialgleichung lautet damit allgemein für den Fall reeller ungleicher Nullstellen $\lambda _{i}$ des charakteristischen Polynoms:

$y_{H}(t)=C_{1}\cdot e^{{\lambda }_{1}\cdot t}+C_{2}\cdot e^{{\lambda }_{2}\cdot t}+\cdots +C_{n}\cdot e^{{\lambda }_{n}\cdot t}$

Die Lösung einer GDGL erfolgt durch Integration. Jede Integration ergibt Integrationskonstanten C_i , deren Anzahl durch die Ordnung der GDGL bestimmt ist. Die Lösung einer GDGL n-ter Ordnung enthält voneinander unabhängige Integrationskonstanten. Diese sind für eine spezielle (partikuläre) Lösung der GDGL abhängig von den Eigenwerten und gegebenen Anfangsbedingungen des Übertragungssystems zu bestimmen.

Die Bestimmung der Integrationskonstanten C_i bei Systemen höherer Ordnung (> 2) ist sehr umständlich. Weitere Informationen liefert die Fachliteratur.

Anfangswertproblem und Integrationskonstanten für eine homogene GDGL 2. Ordnung

Eine homogene GDGL n-ter Ordnung hat n Anfangswerte. Für die homogene GDGL zweiter Ordnung mit zwei vorzugebenden Anfangswerten $y_{0}$ und ${\dot {y}}_{0}$ können die Koeffizienten $C_{1}$ und $C_{2}$ errechnet werden, wenn die Nullstellen des charakteristischen Polynoms bekannt sind.

Für jede Anfangsbedingung ergibt sich eine Gleichung ( $y_{H}(0)=y_{0}$ , ${\dot {y}}_{H}(0)={\dot {y}}_{0}$ ).

Beispiel für eine homogene GDGL mit zwei reellen Nullstellen $\lambda _{1}=-0{,}5$ und $\lambda _{2}=-1$ und Anfangswerten ${\dot {y}}_{0}=1$ ; $y_{0}=1$ :

Lösung $y_{H}(t)$ der homogenen DGL 2. Ordnung:

$y_{H}(t)=C_{1}\cdot e^{-0{,}5\cdot t}+C_{2}\cdot e^{-1\cdot t}$

Berechnung der Koeffizienten:

$y_{0}\ {\xrightarrow {t=0}}1=C_{1}+C_{2}$

${\dot {y}}_{0}{\xrightarrow {d/dt}}1=C_{1}\cdot -0{,}5\cdot e^{-0{,}5\cdot t}+C_{2}\cdot -1\cdot e^{-1\cdot t}$

${\dot {y}}_{0}{\xrightarrow {t=0}}1=C_{1}\cdot -0{,}5+C_{2}\cdot -1=-0{,}5\cdot C_{1}-C_{2}$

Aus den beiden Gleichungen von $y_{0}$ für t=0 und ${\dot {y}}_{0}$ für t=0 lassen sich die Koeffizienten $C_{1}$ und $C_{2}$ bestimmen.

Anmerkung: Die Ableitung $d/dt$ von $e^{-a\cdot t}=-a\cdot e^{-a\cdot t}$

Tabelle: Durch die verschiedenen Arten der Lösungen der quadratischen Gleichung, bedingt durch die Größe der Diskriminante, ergeben sich drei unterschiedliche Fälle der Eigenwerte λ der GDGL wie:

Lösung der homogenen linearen Differenzialgleichung 2. Ordnung mit konstanten Koeffizienten	Nullstellen	Anfangswertproblem Bestimmung C1, C2
Radikand > 0: 2 reelle Nullstellen $y_{H}(t)=C_{1}\cdot e^{{\lambda }_{1}\cdot t}+C_{2}\cdot e^{{\lambda }_{2}\cdot t}$	${\lambda }_{1;2}=-{\frac {a_{1}}{2}}\pm {\sqrt {{\frac {{a_{1}}^{2}}{4}}-a_{0}}}$	$C_{1}={\frac {{\dot {y}}_{0}-y_{0}\cdot \lambda _{2}}{\lambda _{1}-\lambda _{2}}}\,$ $C_{2}=y_{0}-C_{1}\,$
Radikand = 0: 2 gleiche Nullstellen $y_{H}(t)=C_{1}\cdot e^{{\lambda }_{1}\cdot t}+C_{2}\cdot t\cdot e^{{\lambda }_{1}\cdot t}$	$\lambda =\lambda _{1;2}=-{\frac {a_{1}}{2}}$	$C_{1}=y_{0}\,$ $C_{2}={\dot {y}}_{0}-y_{0}\cdot \lambda \,$
Radikand < 0: konjugiert komplexen Nullstellen $y_{H}(t)=e^{\alpha \cdot t}\cdot [C_{1}\cdot \cos(\beta \cdot t)+C_{2}\cdot \sin(\beta \cdot t)]$	${\lambda }_{1;2}=\alpha \pm i\cdot \beta$ $\alpha =-a_{1}/2\,$ $\beta ={\sqrt {a_{0}-{\frac {{a_{1}}^{2}}{4}}}}$	$C_{1}={\frac {{\dot {y}}_{0}-\alpha \cdot y_{0}}{\beta }}$ $C_{2}=y_{0}\,$

Berechnungsbeispiel der Lösung einer GDGL 2. Ordnung mit reellen Nullstellen

* Homogene Lösung der GDGL einer Reihenschaltung von zwei PT1-Gliedern mit Anfangswerten.
* Partikuläre Lösung der GDGL für einen Eingangssprung.

Übertragungsfunktion eines dynamischen Systems bestehend aus zwei PT1-Gliedern: $G(s)={\frac {Y(s)}{U(s)}}={\frac {1}{(2\cdot s+1)\cdot (s+1)}}={\frac {1}{2\cdot s^{2}+3\cdot s+1}}$

Zugehörige systembeschreibende GDGL:

$2\cdot {\ddot {y}}(t)+3\cdot {\dot {y}}(t)+y(t)=u(t)$

Die höchste Ableitung freigestellt:

${\ddot {y}}(t)+1{,}5\cdot {\dot {y}}(t)+0{,}5\cdot y(t)=0{,}5\cdot u(t)\qquad {\text{Koeffizienten:}}\ a_{1}=1{,}5\quad a_{0}=0{,}5\quad b_{0}=0{,}5$

Vorgegeben: Willkürlich gewählte Anfangswerte der Energiespeicher (Integratoren): $y'_{0}(t)=1$ ; $y_{0}(t)=1$ ;
Vorgegeben: Eingangsgröße $u(t)=1$ ist eine normierte Sprungfunktion für .
Gesucht: Homogene Lösung der GDGL $y_{H}(t)$ und partikuläre Lösung $y_{p}(t)$ :

Für die homogene Lösung wird u(t)=0

gesetzt.

Errechnet laut der oben dargestellten Tabelle der homogenen Lösung:

Es ergeben sich zwei reelle Nullstellen: ${\lambda }_{1}=-0{,}5;\quad {\lambda }_{2}=-1$

Errechnet: Die Integrationskonstanten errechnen sich laut Tabelle mit $C_{1}=4$ ; $C_{2}=-3$ .
Analytische homogene Lösung laut Tabelle für zwei reelle Nullstellen:

$y_{H}(t)=C_{1}\cdot e^{{\lambda }_{1}\cdot t}+C_{2}\cdot e^{{\lambda }_{2}\cdot t}\qquad \to$ daraus folgt:

Mit den eingesetzten Zahlenwerten lautet die analytische Lösung der homogenen GDGL:

${\underline {\underline {y_{H}(t)=4\cdot e^{-0{,}5\cdot t}-3\cdot e^{-t}}}}$

Partikuläre Lösung:

Die Berechnung der Systemantwort $y_{p}(t)$ des Eingangs-Ausgangsverhaltens über das Faltungsintegral ist aufwendig.

Einfacher ist die Lösung – wie nachfolgend dargestellt – durch die Anwendung der Laplace-Transformation.

Lösung mittels der Übertragungsfunktion

Die allgemeine Form einer Differentialgleichung mit konstanten Koeffizienten $a_{i}$ der Ausgangsgröße y(t) und mit $b_{i}$ der Eingangsgröße u(t) lautet:

$a_{n}\cdot y^{(n)}(t)+\ \cdots \ +a_{2}\cdot {\ddot {y}}(t)+a_{1}\cdot {\dot {y}}(t)+a_{0}\cdot y(t)=b_{0}\cdot u(t)+b_{1}\cdot {\dot {u}}(t)+b_{2}\cdot {\ddot {u}}(t)+\ \cdots \ +b_{m}\cdot u^{(m)}(t)$ .

Durch Anwendung des Laplace-Differentiationssatzes einer GDGL entstehen algebraische Gleichungen mit sogenannten Zähler- und Nennerpolynomen. $s=\delta +j\cdot \omega$ ist die komplexe Laplace-Variable, die mit einem Exponenten anstelle der Ordnung einer Ableitung steht. Die Übertragungsfunktion G(s) ist definiert als das Verhältnis des Ausgangssignals Y(s) zum Eingangssignal U(s) , wobei die Anfangswerte des Systems gleich Null sind.

$G(s)={\frac {Y(s)}{U(s)}}={\frac {b_{m}s^{m}+\dotsb +b_{2}s^{2}+b_{1}s+b_{0}}{a_{n}s^{n}+\dotsb +a_{2}s^{2}+a_{1}s+a_{0}}}$ .

Die Berechnung des Zeitverhaltens eines Übertragungssystems aus der Übertragungsfunktion G(s) wird üblicherweise für normierte Eingangssignale U(s) durchgeführt. Zur Berechnung der Sprungantwort mit dem Eingangssignal $u(t)=1(t):=U(s):={\frac {1}{s}}$ wird der Übertragungsfunktion der Term ${\frac {1}{s}}$ multiplikativ angehängt. Wird letzteres nicht durchgeführt, erhält man an Stelle der Sprungantwort die Impulsantwort.

Übertragungsfunktion in Polynomdarstellung, Pol- Nullstellendarstellung und Zeitkonstantendarstellung:

Die Pole und Nullstellen der Übertragungsfunktion sind die wichtigsten Kenngrößen des Systemverhaltens. Die Pole (Nullstellen des Nennerpolynoms) sind gleichzeitig die Lösung des Systems und bestimmen das System-Zeitverhalten. Die Nullstellen des Zählerpolynoms haben nur Einfluss auf die Amplituden der Systemantwort.

${\begin{aligned}G(s)&={\frac {Y(s)}{U(s)}}={\frac {b_{m}s^{m}+\dotsb +b_{1}s+b_{0}}{a_{n}s^{n}+a_{n-1}s^{n-1}+\dotsb +a_{1}s+a_{0}}}\\\\&=k\cdot {\frac {(s-s_{n1})(s-s_{n2})\dotsm (s-s_{nm})}{(s-s_{p1})(s-s_{p2})\dotsm (s-s_{pn})}}=K\cdot {\frac {(T_{v1}\cdot s+1)(T_{v2}\cdot s+1)\cdots }{(T_{1}\cdot s+1)(T_{2}\cdot s+1)\cdots }}\end{aligned}}$

Die Lösung erfolgt durch Partialbruch-Zerlegung der Produktdarstellung in einfache additive Terme, die sich leicht in den Zeitbereich transformieren lassen. Die Partialbruch-Zerlegung von Übertragungsfunktionen höherer Ordnung ist nicht immer einfach, insbesondere wenn konjugiert komplexe Nullstellen vorliegen.

Alternativ können Laplace-Transformationstabellen benutzt werden, welche die häufigsten korrespondierenden Gleichungen im Zeitbereich enthalten.

Partikuläre Lösung der GDGL 2. Ordnung mit Hilfe der Laplace-Transformation

Die partikuläre Lösung beschreibt das Übertragungsverhalten des Systems als Funktion des Eingangssignals u(t) und ist meist von hauptsächlichem Interesse. Die Anfangsbedingungen $y_{0}$ und ${\dot {y}}_{0}$ haben dabei den Wert 0.

Lösung der gegebenen GDGL 2. Ordnung:

${\ddot {y}}(t)+a_{1}\cdot {\dot {y}}(t)+a_{0}\cdot y(t)=b_{0}\cdot u(t)$ .

Die Übertragungsfunktion eines Systems entsteht nach dem Differentiationssatz durch Austausch der zeitabhängigen Terme einer GDGL mit den Laplace-Transformierten. Voraussetzung ist, dass die Anfangsbedingung des Systems Null ist. Je nach Grad der Ableitungen einer Funktion y(t) entstehen nach der Transformation folgende Laplace-Transformierte Y(s):

$s^{2}\cdot Y(s)+a_{1}\cdot s\cdot Y(s)+a_{0}\cdot Y(s)=b_{0}\cdot U(s)$

Mit den transformierten Termen kann die Übertragungsfunktion des dynamischen Systems G(s) aufgestellt werden:

$G(s)={\frac {Y}{U}}(s)={\frac {b_{0}}{s^{2}+a_{1}\cdot s+a_{0}}}$

Polynome einer Übertragungsfunktion werden durch Nullstellenbestimmungen in Linearfaktoren (Grundpolynome: Monom, Binom und Trinom) zerlegt. Liegen Zahlenwerte der Koeffizienten einer Übertragungsfunktion 2. Ordnung vor, können die Pole (= Nullstellen im Nenner der Übertragungsfunktion) durch die bekannte Formel zur Lösung einer gemischt-quadratischen Gleichung ermittelt werden.

Durch die verschiedenen Arten der Lösungen der Pole bedingt durch die Größe des Radikanden der quadratischen Gleichung ergeben sich drei unterschiedliche Fälle der Eigenwerte $s_{i}$ (der Pole $s_{pi}$ ) der Übertragungsfunktion. Nachfolgend ist eine Korrespondenztabelle des s-Bereichs mit $Y(s)=U(s)\cdot G(s)$ und des Zeitbereichs für y(t) für einen transformierten Eingangssprung $u(t)=1(t)\to \ U(s):=1/s$ .

Folgende Grundpolynome (Binome und Trinome bei konjugiert komplexen Polen) entstehen in Abhängigkeit von den Nullstellen. Die Lösungen der Übertragungsfunktionen als Sprungantwort im Zeitbereich sind einer Laplace-Transformationstabelle entnommen worden:

Die Laplace-Transformationstabellen können in zwei Formen der Produkt-Darstellung aufgeführt sein, wobei unterschiedliche Faktoren a₀ und K berücksichtigt werden müssen. Die Umrechnung der Pole- Nullstellen-Darstellung in Zeitkonstanten-Darstellung ist einfach, sie sind algebraisch identisch. $T_{i}=1/{s_{pi}}$ .

Pol-Nullstellen-Darstellung (Stabiles System) und Zeitkonstanten-Darstellung:

$G(s)={\frac {a_{0}}{(s+s_{p1})\cdot (s+s_{p2})}}={\frac {K}{(T_{1}\cdot s+1)\cdot (T_{2}\cdot s+1)}}$

Tabelle: Berechnung der Sprungantworten eines Übertragungssystems 2. Ordnung in Abhängigkeit von den Polstellenarten:

Sprungantworten PT2-Glied:
* 1): 2 reelle Polstellen,
* 2): 2 konjugiert komplexe Polstellen.

f(s) Übertragungsfunktion 2. Ordnung Eingangssprung u(t) = 1 := Multiplikation mit 1/s	f(t) Partikuläre Lösung Sprungantwort im Zeitbereich	Bestimmung der Pole $s_{p1}$ und $s_{p2}$ aus der Polynom-Darstellung
2 reelle Polstellen: $Y(s)={\frac {1}{s\cdot (T_{1}\cdot s+1)\cdot (T_{2}\cdot s+1)}}$	$y(t)=1-{\frac {1}{T_{1}-T_{2}}}\cdot$ $\cdot \left(T_{1}\cdot e^{-{\frac {t}{T_{1}}}}-{T_{2}\cdot e^{-{\frac {t}{T_{2}}}}}\right)$	$s_{p1;2}=-{\frac {a_{1}}{2}}\pm {\sqrt {{\frac {{a_{1}}^{2}}{4}}-a_{0}}}$ $T_{1}=-{\frac {1}{s_{p1}}}\quad T_{2}=-{\frac {1}{s_{p2}}}$
2 gleiche Polstellen: $Y(s)={\frac {1}{s\cdot (T\cdot s+1)^{2}}}$	$y(t)=1-{\frac {T+t}{T}}\cdot e^{-{\frac {t}{T}}}$	$s_{p1}=s_{p2}=-{\frac {a_{1}}{2}}$ $T=-{\frac {1}{s_{p}}}$
Konjugiert komplexe Polstellen: $Y(s)={\frac {{\omega _{0}}^{2}}{s\cdot (s^{2}+2\cdot D\cdot \omega _{0}\cdot s+{\omega _{0}}^{2})}}\quad {\text{oder:}}$ $Y(s)={\frac {1}{s\cdot (T^{2}\cdot s^{2}+2\cdot D\cdot T\cdot s+1)}}$ $\omega _{0}={\text{ungedämpfte Kreisfrequenz}},\quad T=1/\omega _{0}$	$y(t)=1-{\frac {1}{\sqrt {1-D^{2}}}}\cdot$ $\cdot e^{-D\cdot \omega _{0}\cdot t}\cdot \sin(\omega _{d}\cdot t+\phi )$ $\omega _{d}=\omega _{0}\cdot {\sqrt {1-D^{2}}}$ $\phi =\arccos(D)$	$s_{p1;2}=-D\cdot \omega _{0}\pm j\omega _{0}\cdot {\sqrt {1-D^{2}}}$ Dämpfung D: $\quad -1<D<1$ $D={\frac {a_{1}}{2\cdot {\sqrt {a_{0}}}}}$

Wird für den Fall der zwei reellen Nullstellen in die Gleichung für $f(t)\to T_{1}=T_{2}$ eingesetzt, entsteht eine Division durch Null $1/(T_{1}-T_{1})$ , was nicht zulässig ist. Als „verschiedene“ Nullstellen gelten bereits Nullstellen, wenn sie sich in einer theoretisch unendlichen Dezimalstelle eines Wertes unterscheiden.

Die Gesamtlösung einer GDGL ergibt sich aus der Überlagerung der Systemantworten auf die Anfangsbedingungen und auf das Eingangssignal:

$y(t)=y_{H}(t)+y_{P}(t)$

Die partikuläre Lösung der GDGL bezieht sich darauf, dass die Anfangswerte $y_{0};\ {\dot {y}}_{0};\ {\ddot {y}}_{0}\dots$ gleich Null sind und das Eingangssignal $u(t)\neq 0$ ist. Sie lässt sich aus der Übertragungsfunktion G(s) bestimmen, indem die Differentialgleichung einer Laplace-Transformation unterzogen wird.

Berechnungsbeispiel der partikulären Lösung einer GDGL 2. Ordnung mit der Laplace-Transformationstabelle

Vorgegeben:

Eingangssignal: Sprungfunktion $U(s)=1/s$ .

Übertragungsfunktion des Systems:

$G(s)={\frac {Y(s)}{U(s)}}={\frac {1}{(2\cdot s+1)\cdot (s+1)}}\qquad T_{1}=2;\,T_{2}=1$

Gesucht: Partikuläre Lösung $y_{P}(t)$ für die gegebene Übertragungsfunktion:

Suchbegriff für die Laplace-Transformationstabelle:

${Y(s)}={\frac {U(s)}{(2\cdot s+1)\cdot (s+1)}}={\frac {1}{s\cdot (2\cdot s+1)\cdot (s+1)}}\qquad T_{1}=2;\,T_{2}=1$

Errechnet:

Die gefundene analytische Gleichung f(t)

der partikulären Lösung laut Transformationstabelle durch Eingabe der Koeffizienten lautet:

$y_{p}(t)=1-{\frac {1}{T_{1}-T_{2}}}\cdot [T_{1}\cdot e^{-{\frac {t}{T_{1}}}}-T_{2}\cdot e^{-{\frac {t}{T_{2}}}}]$ .

Zahlenwerte der Zeitkonstanten eingesetzt:

${\underline {\underline {y_{p}(t)=\ 1-[2\cdot {e^{-{\frac {t}{2}}}-e^{-t}}]}}}$ .

Grafische Darstellung der partikulären Lösung siehe letztes Bild.

Anmerkung: Enthält die Ausgangsgröße eines Übertragungssystems Schwingungsanteile, ergeben sich laut Transformationstabellen aufwendige trigonometrische Gleichungen.

Lösung einer linearen gewöhnlichen Differentialgleichung mittels der numerischen Berechnung

Eine lineare gewöhnliche Differentialgleichung mit konstanten Koeffizienten, die ein dynamisches System mit einem Eingangssignal und einem Ausgangssignal beschreibt, wird in eine Differenzengleichung umgeformt, indem die Differentialquotienten der GDGL durch Differenzenquotienten ausgetauscht werden. Die rekursive Lösung einer Differenzengleichung erster Ordnung erfolgt von einer Anfangsbedingung durch nummerierte Folgegleichungen. Jede Folgegleichung bezieht sich auf das Ergebnis einer zurückliegenden Folgegleichung.

Bei allen numerischen Verfahren wird der Wert der gesuchten Funktion für diskrete Zeitpunkte ermittelt.

Die numerische Gesamtlösung des Systems erfolgt – bei einfachen Differenzengleichungen – rekursiv (sich selbst aufrufend) über viele Berechnungsfolgen $k=(0,1,2,3,\dots ,k_{{\mathrm {max}}})$ in meist je kleinen konstanten Zeitstufen. Die Form der Gesamtlösung ist damit tabellarisch für die gesuchten Werte $y_{(k)}$ (Stützpunkt, Knoten) eines Funktionsverlaufs im zeitlichen Abstand $\Delta t$ .

Das einfachste numerische Verfahren ist das Euler-Streckenzugverfahren. Zu den Einschritt-Verfahren gehören Euler-Vorwärtsverfahren (Explizites Euler-Verfahren), Euler-Rückwärtsverfahren (Implizites Eulerverfahren), Runge-Kutta-Verfahren, Heun-Verfahren. Bei den Mehrschrittverfahren wird die Information aus den zuvor bereits errechneten Stützpunkten gebildet.

Zum weiteren Verständnis der numerischen Berechnung linearer dynamischer Systeme in Form von linearen gewöhnlichen Differentialgleichungen siehe:

Anwendung der numerischen Berechnung nach dem Euler-Streckenzugverfahren:

$\to$ siehe Hauptartikel Differenzengleichung.

Tabelle der Differenzengleichungen dynamischer Systeme (Euler-Rückwärts) und diskreten Signalen:

$\to$ siehe Hauptartikel „Digitaler Regler“: Tabelle Differenzengleichungen.

Beispiel der numerischen Berechnung der Sprungantwort eines Schwingungsgliedes ( $PT2_{{kk}}$ -Glied) als Differenzengleichung 2. Ordnung:

$\to$ Siehe Hauptartikel Differenzengleichung#Differenzengleichungen höherer Ordnung.

Beispiel für die numerische Berechnung eines senkrechten Federpendels mit Anfangswerten für ein System 2. Ordnung nach dem Signalflussplan der Regelungsnormalform:

$\to$ Siehe Hauptartikel Dynamisches System (Systemtheorie)#Federpendel.

Software

Einige CAS können Differentialgleichungen lösen, z.B.:

ExpressionsinBar: desolve(y'=b*y-a,y)
Maple: dsolve
SageMath: desolve
Xcas: desolve(y'=k*y,y)

Siehe auch

Literatur

Herbert Amann: Gewöhnliche Differentialgleichungen, 2. Auflage, Gruyter – de Gruyter Lehrbücher, Berlin New York, 1995, ISBN 3-11-014582-0
Martin Hermann: Numerik gewöhnlicher Differentialgleichungen, Anfangs- und Randwertprobleme, Oldenbourg Verlag, München und Wien, 2004, ISBN 3-486-27606-9
Harro Heuser: Gewöhnliche Differentialgleichungen, Teubner, März 2004, ISBN 3-519-32227-7

Basierend auf einem Artikel in:

Wikipedia.de