Ein Potenzreihenansatz ist ein Lösungsansatz für Differentialgleichungen. Die gesuchte Funktion wird als Potenzreihe mit unbekannten Koeffizienten dargestellt und dann in die Differentialgleichung eingesetzt. Durch Koeffizientenvergleich kann so die Lösung gefunden und in manchen Fällen wieder durch elementare Funktionen ausgedrückt werden.
Im allgemeinen Fall, wenn die Koeffizientenfunktionen meromorph sind wie bei der Fuchsschen Differentialgleichung (zu der die Hypergeometrische Differentialgleichung gehört), muss die Differentialgleichung grundsätzlich im Komplexen (Riemannsche Zahlenkugel) betrachtet werden. Es gibt bei Differentialgleichungen vom Fuchsschen Typ (mit ausschließlich hebbaren Singularitäten auch im Unendlichen) verallgemeinerte Potenzreihenlösungen (siehe Frobenius-Methode) und die lokal als Potenzreihenlösungen gegebenen Fundamentallösungen der Differentialgleichung sind durch Betrachtung von analytischen Fortsetzungen um die singulären Punkte der Koeffizientenfunktionen über Monodromie-Matrizen verbunden.
Die Exponentialfunktion als motivierendes Beispiel
Als einfaches Beispiel betrachten wir folgende Fragestellung: Welche Funktion ergibt abgeleitet ein Vielfaches dieser Funktion? Als Gleichung:
![{\displaystyle {\frac {dy(x)}{dx}}=\lambda \,y(x)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b12ae844ce14d15b2f0504abd0b74b5ad7f94330)
Diese gewöhnliche Differentialgleichung 1. Ordnung ist eindeutig lösbar, wenn noch eine Anfangsbedingung festgelegt wird:
![{\displaystyle y(0)=y_{0}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3b7f7ab4d55beca261df55663bdc6e674d0acef2)
Für
setzen wir nun eine Potenzreihe an:
![{\displaystyle y(x)=\sum _{k=0}^{\infty }a_{k}x^{k}=a_{0}+a_{1}x+a_{2}x^{2}+a_{3}x^{3}+\ldots }](https://wikimedia.org/api/rest_v1/media/math/render/svg/097f3b65978fb2ad38c27fc811e39b3b3bf1071e)
Die Anfangsbedingung übersetzt sich zu
, weil
.
Die Ableitung von
ist folglich:
![{\displaystyle {\frac {dy(x)}{dx}}={\frac {d}{dx}}\sum _{k=0}^{\infty }a_{k}x^{k}=\sum _{k=0}^{\infty }{\frac {d}{dx}}a_{k}x^{k}=\sum _{k=1}^{\infty }a_{k}\,k\,x^{k-1}=a_{1}+2a_{2}x+3a_{3}x^{2}+\ldots }](https://wikimedia.org/api/rest_v1/media/math/render/svg/ec7ded98d9086754dbead46eae7901c8b20a293e)
Eingesetzt in obige Differentialgleichung heißt das:
![{\displaystyle a_{1}+2a_{2}x+3a_{3}x^{2}+\ldots =\lambda a_{0}+\lambda a_{1}x+\lambda a_{2}x^{2}+\lambda a_{3}x^{3}+\ldots }](https://wikimedia.org/api/rest_v1/media/math/render/svg/e1ce96238a6d174a54f309bbc3c687800f545296)
Da dies für alle
gelten soll, müssen die Koeffizienten vor
usw. gleich sein. Folglich ist:
usw. Dies lässt sich umstellen und einsetzen:
,
,
. Allgemein ist:
und somit
für alle
.
Dies ist eine Rekursionsgleichung für die Koeffizienten
und es ergibt sich:
.
Eingesetzt in die Potenzreihe heißt dies:
.
Wenn wir darin die Potenzreihe der Exponentialfunktion wiedererkennen, lässt sich die Lösung noch kompakter schreiben als:
.
Theoretische Begründung
Zur theoretischen Begründung dieses Verfahrens sollte man bereits im Vorfeld wissen, dass es eine holomorphe Lösung gibt, das heißt eine Lösung, die sich in eine Potenzreihe entwickeln lässt.
Natürlich kann man das einfach annehmen, auf Basis dieser Annahme wie im einleitenden Beispiel eine Lösung konstruieren und dann diese durch Einsetzen prüfen. Kann man aber die Rekursion der Koeffizienten nicht auflösen und kann man nur einige Koeffizienten berechnen, so hat man ein Polynom als Approximation einer möglichen Lösung, aber das ist nur sinnvoll, wenn die Existenz einer holomorphen Lösung gesichert ist. Das liefert der folgende Satz:
- Satz: Seien
sowie
gegeben und
holomorph, wobei
und
. Dann existiert genau eine holomorphe Lösung
des Anfangswertproblems
,
- und zwar mindestens auf dem offenen Kreis
.[1]
In obigem Beispiel ist
und
. Für
ist
.
Der durch den Satz zugesicherte Konvergenzradius von
kann also kleiner sein als der tatsächliche Konvergenzradius der Lösung, der im vorliegenden Beispiel bekanntlich unendlich ist. Der Identitätssatz für holomorphe Funktionen zeigt dann, dass die gefundene Lösung auch außerhalb des Konvergenzradius noch das Anfangswertproblem löst, solange man
in einer zusammenhängenden Umgebung des Konvergenzkreises noch bilden kann.
Insbesondere zeigt dieser Satz, dass der Potenzreihenansatz im Falle holomorpher rechter Seite des Anfangswertproblems zum Erfolg führt.
Weiteres Beispiel: Hermitesche Differentialgleichung
Gesucht wird die Lösung der Hermiteschen Differentialgleichung[2]
![{\displaystyle {\frac {d^{2}y(x)}{dx^{2}}}-2\,x\,{\frac {dy(x)}{dx}}+2\,\lambda \,y(x)=0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/417e7d3c7f6e29592741ee542cb25e85c927fc6d)
Man setzt die Lösung als Potenzreihe an:
![{\displaystyle y(x)=\sum _{k=0}^{\infty }{\frac {a_{k}}{k!}}x^{k}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a7dadb219c8b7a7825329926758cf74191f76ea5)
Um die weitere Rechnung einfacher zu gestalten, wurde in diesem Ansatz im Vergleich zum letzten Beispiel ein Faktor
eingeführt.
Folglich ist:
![{\displaystyle 2\,\lambda \,y(x)=2\,\lambda a_{0}+\sum _{k=1}^{\infty }{\frac {2\,\lambda \,a_{k}}{k!}}x^{k}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/16a21ca1e356fb27082317ba2a15a8d692d2bbb0)
![{\displaystyle -2x\,{\frac {dy(x)}{dx}}=-2x\,\sum _{k=1}^{\infty }{\frac {a_{k}}{k!}}k\,x^{k-1}=\sum _{k=1}^{\infty }{\frac {-2a_{k}}{(k-1)!}}\,x^{k}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/dcc393031bee7142f1b0d9382497a848c3d66e1e)
![{\displaystyle {\frac {d^{2}y(x)}{dx^{2}}}=\sum _{k=2}^{\infty }{\frac {a_{k}}{k!}}k\,(k-1)\,x^{k-2}=\sum _{k=2}^{\infty }{\frac {a_{k}}{(k-2)!}}\,x^{k-2}=\sum _{k=0}^{\infty }{\frac {a_{k+2}}{k!}}\,x^{k}=a_{2}+\sum _{k=1}^{\infty }{\frac {a_{k+2}}{k!}}\,x^{k}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/957d52e2f0bcd60a6dce8058eca96e306e3bacfd)
Eingesetzt in die Differentialgleichung heißt das:
![{\displaystyle a_{2}+\sum _{k=1}^{\infty }{\frac {a_{k+2}}{k!}}\,x^{k}\;+\;\sum _{k=1}^{\infty }{\frac {-2a_{k}}{(k-1)!}}\,x^{k}\;+\;2\,\lambda a_{0}+\sum _{k=1}^{\infty }{\frac {2\,\lambda \,a_{k}}{k!}}x^{k}=0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/71ec93b34e809356c6a358f43397642980fed6e4)
![{\displaystyle a_{2}+2\,\lambda a_{0}+\sum _{k=1}^{\infty }\left({\frac {a_{k+2}}{k!}}+{\frac {-2a_{k}}{(k-1)!}}+{\frac {2\,\lambda \,a_{k}}{k!}}\right)x^{k}=0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/dc4712918ed3f49e6ed617664a2ec615a37468b8)
Der Koeffizientenvergleich ergibt für die konstanten Terme (
):
und für alle weiteren (
):
.
Multiplikation mit
ergibt:
, d. h.
.
Sind die Koeffizienten
und
bspw. aus Anfangsbedingungen bekannt, dann lassen sich alle weiteren Koeffizienten
berechnen und ggf. als Reihe zusammenfassen. Die analytische Lösung der Differentialgleichung lautet also:
.
Literatur
- Earl A. Coddington, Norman Levinson: Theory of Ordinary Differential Equations. McGraw–Hill, New York 1955.
- Einar Hille: Ordinary Differential Equations in the Complex Plane, Dover Publications, Mineola, New York, 1976
- Gerald Teschl Ordinary Differential Equations and Dynamical Systems, publisher=American Mathematical Society, Providence, Rhode Island, 2012, ISBN 978-0-8218-8328-0, pdf
Weblinks
- Frobenius Method (MathWorld)
Einzelnachweise
- ↑ W. Walter: Gewöhnliche Differentialgleichungen, Springer-Verlag (1986), ISBN 3-540-16143-0, Kapitel I, §8, Satz II
- ↑ Harro Heuser: Gewöhnliche Differentialgleichungen: Einführung in Lehre und Gebrauch. Teubner, 3. Auflage, 1995, ISBN 3-519-22227-2, S. 262.