Grundlagen der anwendungsbezogenen Hochschulmathematik
7. Gewöhnliche Differentialgleichungen
7.3. Gewöhnliche Differentialgleichungen zweiter und höherer Ordnung
Für Differentialgleichungen höherer Ordnung ist es meistens unmöglich analytische Lösungen zu finden. In diesem Abschnitt stellen wir ein paar Sonderfälle vor, für die analytische Lösungen gefunden werden können. Die meisten dieser Fälle sind lineare Differentialgleichungen. Unser Fokus liegt auf Differentialgleichungen zweiter Ordnung, da diese häufiger praktische Anwendung finden und es wahrscheinlicher ist analytische Lösungen zu finden im Vergleich zu Differentialgleichungen dritter oder höherer Ordnung.

Lösen einer gewöhnlichen homogenen Differentialgleichung
Für lineare Gleichungen zweiter Ordnung gibt es keine einfache Weise eine allgemeine Lösung zu finden. Wir beginnen zuerst mit der Untersuchung einer homogenen Gleichung
\( y’’ + p(x)y’ + q(x)y = 0,\)
wobei \(p\) und \(q\) stetige Funktionen in ihrem Definitionsbereich sind. Dann gilt, dass
1) die Gleichung linear unabhängige Lösungen \(y_1\) und \(y_2\), genannt Fundamental-Lösungen, besitzt. Grob gesagt bedeutet lineare Unabhängigkeit, dass das Verhältnis \(y_2(x)/y_1(x)\) nicht konstant ist, so dass die Lösungen im Grund voneinander verschieden sind.
2) die allgemeine Lösung durch ein beliebiges linear unabhängiges Lösungspaar der Form \(y(x) = C_1y_1(x) + C_2y_2(x) \), wobei \( C_1\) und \( C_2\) Konstanten sind, ausgedrückt werden kann.
3) wenn die Anfangswerte \(y(x_0) = a, y'(x_0) = b\) fix sind, dann ist die Lösung eindeutig.
Eine allgemeine Methode die Fundamentallösungen \(y_1(x)\) und \(y_2(x)\) explizit zu finden existiert nicht. Um die Lösung zu finden, ist ein typischer Ansatz eine gute Vermutung über die Form der Lösung abzugeben und die Details durch Einsetzen in die Gleichung zu überprüfen.
Die obigen Ergebnisse können auch auf homogene Gleichungen höherer Ordnung verallgemeinert werden, aber dann nimmt die Anzahl an benötigten Fundamentallösungen und Anfangswerten in Anbetracht der Ordnung der Gleichung zu.
Beispiel 1.
Die Gleichung \( y’’-y= 0\) hat die Lösungen \( y = e^x\) und \( y = e^{-x}.\) Diese Lösungen sind linear unabhängig, so dass die allgemeine Lösung der Form \( y(x) = C_1e^x + C_2e^{-x}\) ist.

Gleichungen mit konstanten Koeffizienten
Betrachten wir einen relativ simplen Sonderfall einer Gleichung zweiter Ordnung
\( y’’ + py’ + qy = 0.\)
Um die Gleichung zu lösen, nutzen wir folgende Schätzung \( y(x) = e^{\lambda x}\), wobei \( \lambda\) eine unbekannte Konstante ist. Einsetzen der Vermutung in die Gleichung liefert
\( \lambda^2 e^{\lambda x} + p\lambda e^{\lambda x} + qe^{\lambda x} = 0.\)
\( \lambda^2 + p\lambda + q = 0.\)
Die letzte Gleichung heißt charakteristische Gleichung der gewöhnlichen Differentialgleichung. Das Lösen der charakteristischen Gleichung erlaubt es uns die Lösungen für die eigentliche gewöhnliche Differentialgleichung zu finden. Die Wurzel der charakteristischen Gleichung kann in drei Fälle aufgeteilt werden:
1) Die charakteristische Gleichung hat zwei unterscheidbare Wurzelergebnisse. Dann hat die gewöhnliche Differentialgleichung die Lösungen \(y_1(x) = e^{\lambda_1x} \) und \(y_2(x) = e^{\lambda_2x}. \)
2) Die charakteristische Gleichung hat ein Wurzelergebnis. Dann hat die gewöhnliche Differentialgleichung die Lösungen \(y_1(x) = e^{\lambda x} \) und \(y_2(x) = xe^{\lambda x}. \)
3) Die Wurzel der charakteristischen Gleichung hat die Form \(\lambda = a \pm bi.\) Dann hat die gewöhnliche Differentialgleichung die Lösungen \(y_1(x) = e^{ax}\cos(bx) \) und \(y_2(x) = e^{ax}\sin(bx). \)
Der zweite Fall kann durch Einsetzen in die ursprüngliche gewöhnliche Differentialgleichung erhalten werden und der dritte Fall kann durch Anwendung der Euler Formel \( e^{ix}=\cos x+i\sin x\) erhalten werden. Mit kleinen Änderungen können diese Ergebnisse auch für Differentialgleichungen höherer Ordnung verallgemeinert werden.
Dadurch dass die charakteristische Gleichung die exakt gleichen Koeffizienten wie die ursprüngliche gewöhnliche Differentialgleichung besitzt, ist es nicht notwendig diese nochmals für konkrete Beispiel zu erörtern: schreib sie einfach hin, indem du die gewöhnliche Differentialgleichung anschaust!
Beispiel 2.
Lösen wir das Anfangswertproblem
\( \left\{\begin{align}y'' -y' +2y=0 \\ y(0) = 1, y(1)=0 \end{align} \right. \)
Die charakteristische Gleichung ist \(\lambda^2 -\lambda -2 = 0,\) welche die Wurzeln \( \lambda_1 = 2\) und \( \lambda_2 = -1\) hat. Somit ist die allgemeine Lösung \( y(x) = C_1e^{2x} + C_2e^{-x}.\) Die Konstanten können durch Verwendung der Anfangswerte ermittelt werden:
\( \left\{\begin{align}C_1 + C_2=1 \\ e^2C_1 + e^{-1}C_2 = 0 \end{align} \right. \)
\( \left\{\begin{align}C_1 = -\frac{1}{e^3-1} \\ C_2 = \frac{e^3}{e^3-1} \end{align} \right. \)
Somit ist die allgemeine Lösung \( y(x) = \frac{1}{e^3-1} (-e^{2x} + e^{3-x}).\)
Beispiel 3.
Werfen wir einen Blick darauf wie die obigen Ergebnisse sich auf Gleichungen höherer Ordnung auswirken, indem wir Folgendes lösen
\( y^{(4)} - 4y''' +14y'' -20y' +25y = 0.\)
Nun ist die charakteristische Gleichung \( \lambda^4 - 4\lambda^3 +14\lambda^2 -20\lambda +25 = 0,\) welches die Wurzeln \( \lambda_1 = \lambda_2 = 1 + 2i\) und \( \lambda_3 = \lambda_4 = 1 - 2i\) liefert. Somit sind die Fundamentallösungen der gewöhnlichen Differentialgleichung \(e^x\sin(2x)\), \(e^x\cos(2x)\), \(xe^x\sin(2x)\) und \(xe^x\cos(2x)\). Die allgemeine Lösung ist
\( y = C_1e^x\sin(2x) + C_2e^x\cos(2x) + C_3xe^x\sin(2x) + C_4xe^x\cos(2x).\)

Eulers Differentialgleichung
Ein weiterer häufiger Typ von Differentialgleichungen zweiter Ordnung ist Eulers Differentialgleichung
\( x^2y'' + axy' + by = 0,\)
wobei \(a\) und \(b\) Konstanten sind. Eine Gleichung dieser Form wird anhand der Schätzung \(y(x)= x^r\) gelöst. Einsetzen liefert:
\( r^2 + (a-1)r + b = 0.\)
Durch Wurzelziehen dieser Gleichung erhalten wir die Lösungen für die gewöhnliche Differentialgleichung auf folgende Art und Weise:
1) Wenn die Wurzel eindeutig und reell ist, dann \( y_1(x)= |x|^{r_1}\) und \( y_2(x)= |x|^{r_2}\).
2) Wenn die Gleichung nur ein Wurzelergebnis liefert, dann \( y_1(x)= |x|^{r}\) und \( y_2(x)= |x|^{r}\ln |x|\).
3) Wenn das Ergebnis der Wurzel die Form \(r = a \pm bi\) hat, dann \( y_1(x)= |x|^{a}\cos(b\ln |x|)\) und \( y_2(x)= |x|^{a}\sin(b\ln |x|)\).
Beispiel 4.
Lösen wir die Gleichung \( x^2y'' - 3xy' + y = 0.\) Da wir merken, dass die Gleichung Eulers Differentialgleichung ist, gehen wir mit der Schätzung \(y= x^r\) vor. Durch Einsetzen der Schätzung in die Gleichung erhalten wir \( r(r-1)x^r - 3rx^r + x^r = 0 \Rightarrow r^2 - 4r + 1 = 0,\) welches \( r = 2 \pm \sqrt{3}\) liefert. Damit ist die allgemeine Lösung der gewöhnlichen Differentialgleichung
\(y = C_1 x^{2+\sqrt{3}} + C_2x^{2-\sqrt{3}}\).

Inhomogene lineare Differentialgleichungen
Die allgemeine Lösung einer inhomogenen Gleichung
\(y'' + p(x)y' + q(x)y = r(x)\)
ist die allgemeine Lösung der entsprechenden homogenen Gleichung \(+\) partikuläre Lösung zur inhomogenen Gleichung, d. h.
\(y(x) = C_1y_1(x) + C_2y_2(x) + y_0(x)\).
Die partikuläre Lösung \(y_0\) lässt sich normalerweise durch eine Schätzung, die der gleichen Form wie \(r(x)\) mit allgemeinen Koeffizienten entspricht. Durch Einsetzen dieser Schätzung in die gewöhnliche Differentialgleichung können wir die Koeffizienten auflösen, aber nur wenn die Schätzung von der korrekten Form ist.
In der unteren Tabelle haben wir eine Liste an möglichen Schätzungen für Differentialgleichungen zweiter Ordnung mit konstanten Koeffizienten zusammengestellt. Die Form der Schätzung hängt davon ab aus welcher Art von elementaren Funktionen \(r(x)\) besteht. Wenn \(r(x)\) eine Kombination aus mehreren verschiedenen elementaren Funktionen ist, dann müssen wir die entsprechenden Elemente für all diese Funktionen in unsere Vermutung einbeziehen. Die charakteristische Gleichung der entsprechenden homogenen Differentialgleichung ist \(P(\lambda)=\lambda^2+p\lambda+q=0\).
| \(r(x)\) enthält | die Schätzung besteht aus |
|---|---|
| Polynom von Grad \(n\) | \(A_0+A_1x+\dots +A_nx^n +A_{n+1}x^{n+1}\), wenn \(q=P(0)=0\)) |
| \(\sin kx,\ \cos kx\) | \(A\cos kx+B\sin kx\), wenn \(P(ik)\neq 0\) |
| \(\sin kx,\ \cos kx\) | \(Ax\cos kx+Bx\sin kx\), wenn \(P(ik)=0\) |
| \(e^{cx}\sin kx,\ e^{cx}\cos kx\) | \(Ae^{cx}\cos kx+Be^{cx}\sin kx\), wenn \(P(c+ik)\neq 0\) |
| \(e^{kx}\) | \(Ae^{kx}\), wenn \(P(k)\neq 0\) |
| \(e^{kx}\) | \(Axe^{kx}\), wenn \(P(k)=0\) and \(P'(k)\neq 0\) |
| \(e^{kx}\) | \(Ax^2e^{kx}\), wenn \(P(k)=P'(k)=0\) |
Anmerkung: Für Wurzeln von Polynomen zweiten Grades müssen wir beachten, dass
-
\(P(k)=0\) und \(P'(k)\neq 0\) \(\Leftrightarrow\) \(k\in\mathbb{R}\) ist eine Wurzel von \(P\).
-
\(P(k)=P'(k)= 0\) \(\Leftrightarrow\) \(k\in\mathbb{R}\) ist die einzige Wurzel von \(P\).
-
\(P(ik)\neq 0\) \(\Leftrightarrow\) \(ik\in\mathbb{C}\) ist kein Wurzelergebnis \(P\); d.h. \(\sin kx\) und \(\cos kx\) sind keine Lösungen der homogenen Gleichung.
Beispiel 5.
Suchen wir die allgemeine Lösung der gewöhnlichen Differentialgleichung \(y''+y'-6y=r(x)\), wenn
a) \(r(x)=12e^{-x}\)
b) \(r(x)=20e^{2x}\).
Die Lösungen sind von der Form \(y(x)=C_1e^{-3x}+C_2e^{2x}+y_0(x)\).
a) Durch Substitution der Vermutung \(y_0(x)=Ae^{-x},\) erhalten wir \((A -A -6A)e^{-x} =12e^{-x}\), welches für \(A=-2\) lösbar ist.
b) In diesem Fall ist die Schätzung von der Form \(Be^{2x}\) nutzlos, da es ein Teil der Lösung der entsprechenden homogenen Gleichung ist und liefert nur Null, wenn es in die linke Seite der gewöhnlichen Differentialgleichung eingesetzt wird. Hier ist die richtige Schätzung von der Form \(y_0(x)=Bxe^{2x}\). Einsetzen liefert
\[ (4B+2B-6B)xe^{2x}+(4B+B)e^{2x} = 20e^{2x}, \]welches für \(B=4\) lösbar ist.
Verwenden wir diese Werte für \(A\) und \(B\) können wir die allgemeinen Lösungen für die gegebenen Differentialgleichungen niederschreiben.
Beispiel 6.
Finden wir die Lösung der gewöhnlichen Differentialgleichung \(y''+y'-6y=12e^{-x}\) mit den Anfangswerten \(y(0)=0\), \(y'(0)=6\).
Basierend auf den vorherigen Beispielen ist die allgemeine Lösung von der Form \(y(x)=C_1e^{-3x}+C_2e^{2x}-2e^{-x}\). Differenzieren liefert \(y'(x)=-3C_1e^{-3x}+2C_2e^{2x}+2e^{-x}\). Durch die Anfangswerte erhalten wir das folgende Paar an Gleichungen:
\[ \begin{cases} 0=y(0)=C_1+C_2-2 &\\ 6=y'(0)=-3C_1+2C_2+2, &\\ \end{cases} \]welches für \(C_1=0\) und \(C_2=2\) lösbar ist. Damit ist die Lösung des Anfangswertproblems \(y(x)=2e^{2x}-2e^{-x}\).
Beispiel 7.
Eine typische Anwendung einer inhomogenen gewöhnlichen Differentialgleichung zweiter Ordnung ist ein RLC Stromkreislauf, der aus einem Widerstand (mit Widerstand \( R\)), einer Spule (mit Induktion \( L \)), einem Kondensator (mit Kapazität \( C \)) ) und einer zeitabhängigen elektromotorischen Kraft \( E(t)\) besteht Der elektrische Strom \( y(t)\) im Kreislauf erfüllt die gewöhnliche Differentialgleichung \[ Ly''+Ry'+\frac{1}{C}y=E'(t).\] Lösen wir diese gewöhnliche Differentialgleichung mit künstlich gewählten numerischen Werten von der Form \[ y''+10y'+61y=370\sin t.\]
Der homogene Teil hat eine charakteristische Gleichung von der Form \( \lambda^2+10\lambda +61=0\) mit Lösungen \( \lambda = -5\pm 6i\). Dies liefert die Lösungen \( y_1(t)=e^{-5t}\cos(6t)\) und \( y_2(t)=e^{-5t}\sin(6t) \) für die homogene Gleichung. Für eine partikuläre Lösung probieren wir \(y_0(t)=A\cos t +B\sin t\). Durch Einsetzen dieses in die inhomogene gewöhnliche Differentialgleichung und sammeln ähnlicher Terme liefert \[ (60A+10B)\cos t +(60B-10A)\sin t = 370\sin t. \] Die Gleichung ist für alle \( t\) erfüllt, (nur) wenn
\[ \begin{cases} 60A+10B=0 &\\ -10A+60B=370. &\\ \end{cases} \]welches für \( A=-1\) und \( B=6\) lösbar ist. Somit ist die allgemeine Lösung \[ y(t)=e^{-5t}(C_1\cos(6t)+C_2\sin(6t)) -\cos t+6\sin t .\] Merke, dass die exponentiellen Terme sehr schnell gegen Null gehen und schließlich der Stromfluss von der Form \[ y(t)\approx -\cos t+6\sin t\] oszilliert.
