From 698aec2e06b917a365e7483935d071fac5a6e142 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Andreas=20M=C3=BCller?= Date: Sat, 27 Nov 2021 21:25:46 +0100 Subject: Verallgemeinerte Potenzreihen, Bessel-Funktionen --- buch/chapters/050-differential/bessel.tex | 330 +++++++++++++++++++- .../050-differential/potenzreihenmethode.tex | 343 +++++++++++++++++++++ 2 files changed, 671 insertions(+), 2 deletions(-) (limited to 'buch/chapters/050-differential') diff --git a/buch/chapters/050-differential/bessel.tex b/buch/chapters/050-differential/bessel.tex index 25092ca..19460f5 100644 --- a/buch/chapters/050-differential/bessel.tex +++ b/buch/chapters/050-differential/bessel.tex @@ -5,11 +5,337 @@ % \section{Bessel-Funktionen \label{buch:differntialgleichungen:section:bessel}} +Die Besselsche Differentialgleichung +erlaubt Wellen mit zylindrischer +Symmetrie und die Strömung in einem zylindrischen Rohr zu beschreiben. +Die Auflösung eines optischen Systems wird durch die Beugung limitiert, +die Helligkeitskverteilung des Bildes einer Punktquelle ist +zylindersymmetrisch und kann mit Hilfe von Lösungen der Besselschen +Differentialgleichung beschrieben werden. +Die Besselsche Differentialgleichung hat im Allgemeinen keine Lösung, +die sich durch bekannte Funktionen ausdrücken lassen, es ist also +nötig, eine neue Familie von speziellen Funktionen zu definieren, +die Bessel-Funktionen. \subsection{Die Besselsche Differentialgleichung} % XXX Wo taucht diese Gleichung auf +Die Besselsche Differentialgleichung ist die Differentialgleichung +\[ +x^2\frac{d^2y}{dx^2} + x\frac{dy}{dx} + (x^2-\alpha^2)y = 0 +\] +zweiter Ordnung +für eine auf dem Interval $[0,\infty)$ definierte Funktion $y(x)$. +Der Parameter $\alpha$ ist eine beliebige komplexe Zahl $\alpha\in \mathbb{C}$, +die Lösungsfunktionen hängen daher von $\alpha$ ab. + +Die Besselsche Differentialgleichung ist eine Differentialgleichung +der Art~\eqref{buch:differentialgleichungen:eqn:dglverallg} mit +\[ +p(x) = 1 +\qquad\text{und}\qquad +q(x) = x^2-\alpha^2. +\] +Nach den Ausführungen von +Abschnitt~\ref{buch:differentialgleichungen:subsection:verallgemeinrt}, +muss die Lösung in der Form einer verallgemeinerten Potenzreihe +gesucht werden. +Dazu muss zunächst die Indexgleichung +\[ +0 += +X(X-1) + Xp_0 + q_0 += +X(X-1) + X - \alpha^2 += +X^2-\alpha^2 += +(X-\alpha)(X+\alpha) +\] +gelöst werden. +Die Nullstellen sind offenbar $\varrho_1=\alpha$ und $\varrho_2=-\alpha$. + +Die beiden Vorzeichen der Nullstellen der Indexgleichung führen +auf die gleiche Differentialgleichung. +Der Lösungsraum der Differentialgleichung ist natürlich trotzdem +zweidimensional, so dass es immer noch möglich ist, den +beiden Nullstellen der Indexgleichung verschiedene Lösungen +zuzuordnen. +Die Diskussion in +Abschnitt~\ref{buch:differentialgleichungen:subsection:verallgemeinrt} +hat Kriterien ergeben, unter denen zwei linear unabhängige Lösungen +mit Hilfe einer verallgemeinerten Potenzreihe gefunden werden können. +Falls nur eine solche Lösung gefunden werden kann, wird sie der grösseren +der beiden Zahlen $\alpha$ und $-\alpha$ zugeordnet +(oder $0$, falls $\alpha=-\alpha=0$). +Eine weitere Lösung kann mit Hilfe analytischer Fortsetzung gefunden werden, +wie später gezeigt wird. + +Für nicht reelles $\alpha$ kann $\varrho_1-\varrho_2=2\alpha$ keine +Ganzzahl sein, es ist also garantiert, dass zwei linear unabhängig +Lösungen der Form +\begin{equation} +y_1(x) = x^\alpha\sum_{k=0}^\infty a_kx^k +\qquad\text{und}\qquad +y_2(x) = x^{-\alpha}\sum_{k=0}^\infty b_kx^k. +\label{buch:differentialgleichungen:eqn:besselloesungen} +\end{equation} +existieren. + +Für reelles $\alpha\in\mathbb{R}$ gibt es zwei Lösungen der +Form~\eqref{buch:differentialgleichungen:eqn:besselloesungen} +genau dann, wenn $\varrho_1-\varrho_2$ keine Ganzzahl ist. +Nur eine Lösung kann man finden, wenn +\[ +\alpha-(-\alpha)=2\alpha \in \mathbb{Z} +\qquad\Rightarrow\qquad +\alpha = \frac{k}{2},\quad k\in\mathbb{Z} +\] +ist. + + \subsection{Bessel-Funktionen erster Art} -% XXX Result der Potenzreihenentwicklung +Für $\alpha \ge 0$ gibt es immer mindestens eine Lösung der Besselgleichung +als verallgemeinerte Potenzreihe mit $\varrho=\alpha$. +Die Funktion $q(x)=x^2-\alpha^2$ ist ein Polynom, die einzigen +von $0$ verschiedenen Koeffizienten sind $q_0=-\alpha^2$ +und $q_2=1$. +Für den ersten Koeffizienten $a_0$ gibt es keine Einschränkungen, +wir wählen $a_0=1$. + +Die Rekursionsformel für $n=1$ ist +\[ +F(\varrho+1) a_1 = (\varrho p_1+q_1)a_0, +\] +aber die Koeffizienten $p_1$ und $q_1$ verschwinden beide und damit +die ganze rechte Seite. +Da $F(\varrho+1)\ne 0$ ist, folgt dass $a_1=0$ sein muss. + +% Fall n=1 gesondert behandeln + +\subsubsection{Der allgemeine Fall} +Für die höheren Potenzen $n>1$ wird die Rekursionsformel für die +Koeffizienten $a_n$ der verallgemeinerten Potenzreihe +\[ +a_{n} = +-\frac{ q_2 a_{n-2} }{F(\varrho+n)} += +-\frac{a_{n-2}}{(\varrho+n)^2-\alpha^2} += +-\frac{a_{n-2}}{\varrho^2 + 2\varrho n+n^2-\alpha^2} += +-\frac{a_{n-2}}{n(n+2\varrho)}. +\] +Im letzten Schritt haben wir verwendet, dass $\varrho=\pm\alpha$ +und damit $\varrho^2=\alpha^2$ ist. +Daraus folgt wegen $a_1=0$, dass auch $a_{2k+1}=0$ für alle $k$. +Damit können wir jetzt die Reihe hinschreiben: +\begin{align*} +y(x) +&= +x^{\varrho}\biggl( +1 +- +\frac{1}{2(2+2\varrho)} x^2 ++ +\frac{1}{2(2+2\varrho)4(4+2\varrho)} x^4 +- +\frac{1}{2(2+2\varrho)4(4+2\varrho)6(6+2\varrho)} x^6 ++ +\dots +\biggr) +\\ +&= +x^{\varrho} +\biggl( +1 ++ +\frac{(-x^2/4)}{1\cdot (1+\varrho)} ++ +\frac{(-x^2/4)^2}{1\cdot 2\cdot (1+\varrho)\cdot(2-\varrho)} ++ +\frac{(-x^2/4)^3}{1\cdot 2\cdot 3\cdot (1+\varrho)\cdot(2+\varrho)\cdot(3+\varrho)} ++ +\dots +\biggr) +\\ +&= +x^\varrho\biggl( +1 ++ +\frac{1}{(\varrho+1)}\frac{(-x^2/4)}{1!} ++ +\frac{1}{(\varrho+1)(\varrho+2)}\frac{(-x^2/4)^2}{2!} ++ +\frac{1}{(\varrho+1)(\varrho+2)(\varrho+3)}\frac{(-x^2/4)^3}{3!} ++ +\dots +\biggr) +\\ +&= +x^\varrho \sum_{k=0}^\infty +\frac{1}{(\varrho+1)_k} \frac{(-x^2/4)}{k!} += +\mathstrut_0F_1\biggl(;\varrho+1;-\frac{x^2}{4}\biggr) +\end{align*} +Falls also $\alpha$ kein ganzzahliges Vielfaches von $\frac12$ ist, finden +wir zwei Lösungsfunktionen +\begin{align} +J_\alpha(x) +&= +x^{\alpha\phantom{-}} +\sum_{k=0}^\infty +\frac{1}{(\alpha+1)_k} +\frac{(-x^2/4)^k}{k!} += +\mathstrut_0F_1\biggl(;\alpha+1;-\frac{x^2}{4}\biggr), +\label{buch:differentialgleichunge:bessel:erste} +\\ +J_{-\alpha}(x) +&= +x^{-\alpha} \sum_{k=0}^\infty +\frac{1}{(-\alpha+1)_k} \frac{(-x^2/4)^k}{k!} += +\mathstrut_0F_1\biggl(;-\alpha+1;-\frac{x^2}{4}\biggr). +\label{buch:differentialgleichunge:bessel:zweite} +\end{align} +Die Funktionen $J_{\pm\alpha}(x)$ heissen {\em Bessel-Funktionen +der Ordnung $\alpha$}. + +\subsubsection{Der Fall $\alpha=0$} +Im Fall $\alpha=0$ hat das Indexpolynom eine doppelte Nullstelle, wir +können daher nur eine Lösung erwarten. +Im Fall $\alpha=0$ wird das Produkt im Nenner zu $n!$, so dass die +Lösungsfunktion +\[ +J_0(x) += +\sum_{k=0}^\infty +\frac{(-1)^k}{(k!)^2} +\biggl(\frac{x}{2}\biggr)^{2k} +\] +geschrieben werden kann. + +% XXX Zweite Lösung explizit durchrechnen + +\subsubsection{Der Fall $\alpha=p$, $p\in\mathbb{N}, p > 0$} +In diesem Fall kann nur die erste +Lösung~\eqref{buch:differentialgleichunge:bessel:erste} +verwendet werden. +Die Pochhammer-Symbole im Nenner können ebenfalls als +Quotient +\[ +\frac{1}{(p+1)_k} += +\frac{1}{(p+1)\cdot(p+k)} += +\frac{p!}{(p+k)!} +\] +von Fakultäten geschrieben werden. +Damit erhält die Lösungsfunktion die Form +\[ +J_p(x) += +\sum_{k=0}^\infty +\frac{(-1)^k}{k!(p+k)!}\biggl(\frac{x}{2}\biggr)^{p+2k}. +\] + +\subsubsection{Der Fall $\alpha=n+\frac12$, $n\in\mathbb{N}$} +Obwohl $2\alpha$ eine Ganzzahl ist, sind die beiden Lösungen +\label{buch:differentialgleichunge:bessel:erste} +und +\label{buch:differentialgleichunge:bessel:zweite} +linear unabhängig. + +Man kann zeigen, dass sich die Lösungsfunktionen in diesem Fall +durch bereits bekannte elementare Funktionen ausdrücken lassen. +Wir rechnen dies für $n=0$ nach. +Zunächst drücken wir die Pochhammer-Symbole im Nenner anders aus. +Es ist +\begin{align*} +\biggl(\frac12 + 1\biggr)_k +&= +\biggl(\frac12 + 1\biggr) +\biggl(\frac12 + 2\biggr) +\cdots +\biggl(\frac12 + k\biggr) += +\frac{1}{2^k}\bigl(3\cdot 5\cdot\ldots\cdot (2k+1)\bigr) += +\frac{(2k+1)!}{2^{2k+1}\cdot k!} +\\ +\biggl(-\frac12 + 1\biggr)_k +&= +\biggl(-\frac12 + 1\biggr) +\biggl(-\frac12 + 2\biggr) +\cdots +\biggl(-\frac12 + k\biggr) +\\ +&= +\biggl(\frac12 + 0\biggr) +\biggl(\frac12 + 1\biggr) +\cdots +\biggl(\frac12 + k-1\biggr) += +\frac{1}{2^k}\bigl(1\cdot 3 \cdot\ldots\cdot (2(k-1)+1)\bigr) += +\frac{(2k-1)!}{2^{2k}\cdot (k-1)!} +\end{align*} +Damit können jetzt die Reihenentwicklungen der Lösung wie folgt +umgeformt werden +\begin{align*} +y_1(x) +&= +\sqrt{x} +\sum_{k=0}^\infty +\frac{1}{(\alpha+1)_k} +\frac{(-x^2/4)^k}{k!} += +\sqrt{x} +\sum_{k=0}^\infty +\frac{2^{2k+1}k!}{(2k+1)!} +\frac{(-x^2/4)^k}{k!} += +\sqrt{x} +\sum_{k=0}^\infty +(-1)^k +\frac{2\cdot x^{2k}}{(2k+1)!} +\\ +&= +\frac{1}{2\sqrt{x}} +\sum_{k=0}^\infty +(-1)^k +\frac{x^{2k+1}}{(2k+1)!} += +\frac{1}{2\sqrt{x}} \sin x +\\ +y_2(x) +&= +\frac{1}{\sqrt{x}} +\sum_{k=0}^\infty +\frac{2^{2k}\cdot (k-1)!}{(2k-1)!} +\frac{(-x^2/4)^k}{k!} += +\frac{1}{\sqrt{x}} +\sum_{k=0}^\infty +(-1)^k +\frac{x^{2k}}{(2k-1)!\cdot k} +\\ +&= +\frac{2}{\sqrt{x}} +\sum_{k=0}^\infty +(-1)^k +\frac{x^{2k}}{(2k-1)!\cdot 2k} += +\frac{2}{\sqrt{x}} \cos x. +\end{align*} + +% XXX Nachrechnen, dass diese Funktionen +% XXX Lösungen der Differentialgleichung sind + +\subsection{Analytische Fortsetzung und Bessel-Funktionen zweiter Art} + + + + -\subsection{Analytische Fortsetzung und Bessel-Funktioenn zweiter Art} diff --git a/buch/chapters/050-differential/potenzreihenmethode.tex b/buch/chapters/050-differential/potenzreihenmethode.tex index 6d30129..247b962 100644 --- a/buch/chapters/050-differential/potenzreihenmethode.tex +++ b/buch/chapters/050-differential/potenzreihenmethode.tex @@ -238,4 +238,347 @@ Die Newtonsche Reihe für $(1-t)^\alpha$ ist der Wert der hypergeometrischen Funktion $\mathstrut_1F_0$. \end{satz} +% +% Verallgemeinerte Potenzreihen +% +\subsection{Lösung mit verallgemeinerten Potenzreihen +\label{buch:differentialgleichungen:subsection:verallgemeinrt}} +In vielen Anwendungsfällen hat die Differentialgleichung die Form +\begin{equation} +x^2y'' + p(x)xy' + q(x)y = 0, +\label{buch:differentialgleichungen:eqn:dglverallg} +\end{equation} +gesucht ist eine Lösung $y(x)$ auf dem Intervall $[0,\infty)$. +Für die folgende Diskussion nehmen wir an, dass sich die Funktionen +$p(x)$ und $q(x)$ in konvergente Potenzreihen +\begin{align*} +p(x)&=\sum_{k=0}^\infty p_kx^k = p_0+p_1x+p_2x^2+p_3x^3+\dots +\\ +q(x)&=\sum_{k=0}^\infty q_kx^k = q_0+q_1x+q_2x^2+q_3x^3+\dots +\end{align*} +entwickeln lassen. + +\subsubsection{Die Potenzreihenmethode funktioniert nicht} +Für die Differentialgleichung +\eqref{buch:differentialgleichungen:eqn:dglverallg} +funktioniert die Potenzreihenmethod oft nicht. +Sind die Funktionen $p(x)$ und $q(x)$ zum Beispiel Konstante +$p(x)=p_0$ und $q(x)=q_0$, dann führt der Potenzreihenansatz +\[ +y(x) = \sum_{k=0}^\infty a_kx^k +\] +auf die Gleichung +\begin{align*} +x^2\sum_{k=0}^\infty a_kk(k-1)x^{k-2} ++ +p_0x\sum_{k=0}^\infty a_kkx^{k-1} ++ +q_0\sum_{k=0}^\infty a_kx^k +&= +0 +\\ +\Rightarrow\qquad +\sum_{k=0}^\infty\bigl( +k(k-1) ++ +p_0k ++ +q_0 +\bigr)a_kx^k +&= +0. +\end{align*} +Durch Koeffizientenvergleich folgt dann, dass für jedes $k$ mindestens +eine der Gleichungen +\[ +k(k-1) +p_0k +q_0 = k^2 + (p_0-1)k +q_0 = 0 +\qquad\text{und}\qquad +a_k=0 +\] +erfüllt sein muss. +Die erste Gleichung ist eine quadratische Gleichung in $k$, es gibt also +höchstens zwei Koeffizienten, für die die erste Gleichung erfüllt sein +kann, für die also auch die Koeffizienten $a_k\ne 0$ sein können. +Sind die Lösungen nicht ganzzahlig, dann müssen alle Koeffizienten +$a_k=0$ sein, die einzige Potenzreihe ist die triviale Funktion $y(x)=0$. + +\subsubsection{Verallgemeinerte Potenzreihe} +Für Differentialgleichungen der Art +\eqref{buch:differentialgleichungen:eqn:dglverallg} +ist also ein anderer Ansatz nötig. +Die Schwierigkeit bestand darin, dass die Gleichungen für die einzelnen +Koeffizienten $a_k$ voneinander unabhängig waren. +Mit einem zusätzlichen Potenzfaktor $x^\varrho$ mit nicht +notwendigerweise ganzzahligen Wert kann die nötige Flexibilität +erreicht werden. +Wir verwenden daher den Ansatz +\[ +y(x) += +x^\varrho \sum_{k=0}^\infty a_kx^k += +\sum_{k=0}^\infty a_k x^{\varrho+k} +\] +und versuchen nicht nur die Koeffizienten $a_k$ sondern auch den +Exponenten $\varrho$ zu bestimmen. +Durch Modifikation von $\varrho$ können wir immer erreichen, dass +$a_0\ne 0$ ist. + +Die Ableitungen von $y(x)$ mit der zugehörigen Potenz von x sind +\begin{align*} +xy'(x) +&= +\sum_{k=0}^\infty +(\varrho+k)a_kx^{\varrho+k} += +\sum_{k=1}^\infty +(\varrho+k-1)a_{k-1}x^{\varrho+k} +\\ +x^2y''(x) +&= +\sum_{k=0}^\infty +(\varrho+k)(\varrho+k-1)a_kx^{\varrho+k}. +\end{align*} +Diese Ableitungen setzen wir jetzt in die Differentialgleichung ein, +die dadurch zu +\begin{equation} +\sum_{k=0}^\infty (\varrho+k)(\varrho+k-1) a_k x^{\varrho+k} ++ +\sum_{k=0}^\infty \sum_{l=0}^\infty (\varrho+k) p_l a_kx^{\varrho+k+l} ++ +\sum_{k=0}^\infty \sum_{l=0}^\infty q_l a_k x^{\varrho+k+l} += +0 +\label{buch:differentialgleichungen:eqn:veralgpotenzsumme} +\end{equation} +wird. + +Ausgeschrieben geben die einzelnen Terme +\begin{align*} +0 +&= +\varrho(\varrho-1)a_0x^\varrho ++ +(\varrho+1)\varrho a_1x^{\varrho+1} ++ +(\varrho+2)(\varrho+1)a_2x^{\varrho+2} ++ +(\varrho+3)(\varrho+2)a_3x^{\varrho+3} ++ +\dots +\\ +&+ +\varrho p_0 a_0 x^{\varrho} ++ +\bigl((\varrho +1)a_1p_0 + \varrho a_0 p_1\bigr) x^{\varrho+1} ++ +\bigl((\varrho +2)a_2p_0 + (\varrho+1)a_1p_1 + \varrho a_0 p_2\bigr) x^{\varrho+2} ++ +\dots +\\ +&+ +q_0a_0x^{\varrho} ++ +(q_0a_1+q_1a_0) x^{\varrho+1} ++ +(q_0a_2+q_1a_1+q_2a_0) x^{\varrho+2} ++ +(q_0a_3+q_1a_2+q_2a_1+q_3a_0) x^{\varrho+3} ++ +\dots +\end{align*} +Fasst man die Terme mit gleichem Exponenten zusammen, findet man +\begin{align*} +0 +&= +\bigl( +\varrho(\varrho-1) + \varrho p_0 + q_0 +\bigr)a_0 x^{\varrho} +\\ +&+ +\bigl( +((\varrho+1)\varrho ++ +(\varrho+1) p_0 ++ +q_0) a_1 ++ +( \varrho p_1 + q_1)a_0 +\bigr)x^{\varrho+1} +\\ +&+ +\bigl( +( +(\varrho+2)(\varrho+1) ++ +(\varrho+2)p_0 ++ +q_0)a_2 ++ +(\varrho+1)p_1 a_1 ++ +\varrho p_2 a_0 ++q_1a_1+q_2a_0 +\bigr)x^{\varrho+2} +\\ +&+\dots +\end{align*} +Der Koeffizientenvergleich ergibt dann +\[ +\renewcommand{\arraycolsep}{0pt} +\begin{array}{rcrlcrlcrl} +0&\mathstrut=\mathstrut&(\varrho(\varrho-1) + \varrho p_0 + q_0)&a_0 + & & & + && +\\ +0&\mathstrut=\mathstrut&((\varrho+1)\varrho + \varrho p_0 + q_0)&a_1 + &\mathstrut+\mathstrut&(\varrho p_1+q_1)&a_0 + & & +\\ +0&\mathstrut=\mathstrut&((\varrho+2)(\varrho+1) + \varrho p_0 + q_0)&a_2 + &\mathstrut+\mathstrut&((\varrho+1)p_1+q_1)&a_1 + &\mathstrut+\mathstrut&(\varrho p_2+q_0)&a_0 +\end{array} +\] + +Diese Rechnung kann man auch allgemein durchführen. +Für den Koeffizientenvergleich müssen die Terme in +\eqref{buch:differentialgleichungen:eqn:veralgpotenzsumme} +mit gleicher +Potenz $x^{\varrho+n}$ zusammengefasst werden. +Dazu schreiben wir zunächst die Summen alle so, dass die Potenz von $x$ +in der Form $x^{\varrho+n}$ auftritt. +So entsteht die Gleichung +\begin{align*} +\sum_{n=0}^\infty +(\varrho+n)(\varrho+n-1) a_n x^{\varrho+n} ++ +\sum_{n=0}^\infty +\biggl( +\sum_{l=0}^n +(\varrho+n-l) p_{n-l} a_{l} +\biggr) +x^{\varrho+n} ++ +\sum_{n=0}^\infty +\biggl(\sum_{l=0}^n q_{n-l} a_{l}\biggr) +x^{\varrho+n} +&= +0 +\end{align*} +Jetzt kann der Koeffizientenvergleich durchgeführt werden. +Der Koeffizient von $x^{\varrho+n}$ ist +\[ +(\varrho+n)(\varrho+n-1) a_n x^{\varrho+n} ++ +\sum_{l=0}^n +(\varrho+n-l) p_{n-l} a_{l} ++ +\sum_{l=0}^n q_{n-l} a_{l}. +\] +Alle diese Koeffizienten müssen verschwinden. +Indem wir die Terme in den beiden Summen über $l$ zusammenfassen, +erhalten wir die Gleichungen +\begin{equation} +\bigl( +(\varrho+n)(\varrho + n-1) ++ +\varrho p_0 ++ +q_0 +\bigr)a_n ++ +\sum_{l=0}^{n-1} +\bigl( +(\varrho+n-l) p_{n-l} ++ +q_{n-l} +\bigr) a_{l} += 0, +\label{buch:differentialgleichungen:eqn:verallgkoefgl} +\end{equation} +die für jedes $n$ erfüllt sein müssen. + +\subsubsection{Indexgleichung} +Die Gleichungen~\eqref{buch:differentialgleichungen:eqn:verallgkoefgl} +müssen erfüllt sein, wenn eine Lösung in Form einer verallgemeinerten +Potenzreihe existieren soll. +Der Koeffizient $a_n$ mit dem grössten $n$ in jeder Gleichung hat +den gemeinsamen Faktor $F(\varrho+n)$ für das Polynom +\[ +F(X) = X(X+1) +Xp_0 + q_0. +\] +Da wir in der Definition einer verallgemeinerten Potenzreihe vorausgesetzt +haben, dass $a_0\ne 0$ sein muss, ist der Ansatz überhaupt nur dann +erfolgreich, wenn \begin{equation} +F(\varrho) = \varrho(\varrho-1) + \varrho p_0 + q_0 = 0 +\label{buch:differentialgleichungen:eqn:indexgleichung} +\end{equation} +gilt. +Die Gleichung~\eqref{buch:differentialgleichungen:eqn:indexgleichung} +heisst die {\em Indexgleichung}. +Der Exponent $\varrho$ muss also eine Nullstelle der Indexgleichung sein. + +Die Indexgleichung ist eine quadratische Gleichung und hat daher +im allgemeinen zwei Lösungen. +Wir bezeichnen die beiden Nullstellen mit $\varrho_1$ und $\varrho_2$. +Wenn $p_0$ und $q_0$ reell sind, sind die Nullstellen entweder reell +oder konjugiert komplex. + +\subsubsection{Rekursive Bestimmung der $a_n$} +Der Koeffizient $a_{n}$ kann nur dann aus den vorangegangene +Koeffizienten $a_{n-1},a_{n-2},\dots$ bestimmt werden, wenn +$F(\varrho+n)\ne 0$ ist. +In diesem Fall gilt +\begin{equation} +a_n += +\frac{1}{F(\varrho+n)} +\sum_{l=0}^{n-1}\bigl( (\varrho+n-l)p_{n-l} + q_{n-l}\bigr)a_l. +\label{buch:differentialgleichungen:eqn:anrekursion} +\end{equation} +Dies funktioniert aber nur, wenn $F(\varrho+n)\ne 0$ für alle +natürlichen $n > 0$ gilt. +Dies ist gleichbedeutend damit, dass die Differenz $\varrho_1-\varrho_2$ +keine ganze Zahl ist. + +\begin{itemize} +\item +Fall 1: $\varrho_1-\varrho_2$ ist keine ganze Zahl. +In diesem Fall lassen sich zwei Lösungen +\begin{align*} +y_1(x) &= x^{\varrho_1}\sum_{k=0}^\infty a_k x^k +\\ +y_2(x) &= x^{\varrho_2}\sum_{k=0}^\infty b_k x^k +\end{align*} +bestimmen, wobei die Koeffizienten $a_k$ und $b_k$ für $k>0$ durch +die Rekursionformel~\eqref{buch:differentialgleichungen:eqn:anrekursion} +aus $a_0$ und $b_0$ bestimmt werden müssen. + +\item +Fall 2: $\varrho$ ist eine doppelte Nullstelle ($\varrho_1-\varrho_2=0$). +In diesem Fall kann nur eine Lösung als verallgemeinerte Potenzreihe +gefunden werden. +Um eine zweite Lösung zu finden, muss die Technik der analytischen +Fortsetzung verwendet werden, die in +Kapitel~\ref{buch:chapter:funktionentheorie} +dargestellt werden. + +\item +Fall 3: $\varrho_1-\varrho-2$ ist eine positive ganze Zahl. +In diesem Fall ist im Allgemeinen nur eine Lösung in Form einer +verallgemeinerten Potenzreihe möglich. +Auch hier müssen Techniken der Funktionentheorie aus +Kapitel~\ref{buch:chapter:funktionentheorie} +verwendet werden, um eine zweite Lösung zu finden. + +\end{itemize} + +Wenn $\varrho_1-\varrho_2$ eine negative ganze Zahl ist, kann man die +beiden Nullstellen vertauschen. +Es folgt dann, dass es eine + + + + -- cgit v1.2.1