From a08ceeb804e40e13673339aace0c554eb50e921a Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Andreas=20M=C3=BCller?= Date: Fri, 31 Dec 2021 22:55:49 +0100 Subject: complete stuff on hypergeometric differential equations --- buch/chapters/050-differential/arcsin.maxima | 5 + buch/chapters/050-differential/bessel.tex | 1 + .../chapters/050-differential/hypergeometrisch.tex | 626 +++++++++++++++++++-- .../050-differential/potenzreihenmethode.tex | 1 + buch/chapters/050-differential/w.maxima | 101 ++++ 5 files changed, 681 insertions(+), 53 deletions(-) create mode 100644 buch/chapters/050-differential/arcsin.maxima create mode 100644 buch/chapters/050-differential/w.maxima (limited to 'buch/chapters/050-differential') diff --git a/buch/chapters/050-differential/arcsin.maxima b/buch/chapters/050-differential/arcsin.maxima new file mode 100644 index 0000000..260d8a9 --- /dev/null +++ b/buch/chapters/050-differential/arcsin.maxima @@ -0,0 +1,5 @@ +u: asin(x); +up: diff(u,x); +upp: diff(up,x); + +(1-x^2)*upp-x*up ; diff --git a/buch/chapters/050-differential/bessel.tex b/buch/chapters/050-differential/bessel.tex index 16527ad..55eb78a 100644 --- a/buch/chapters/050-differential/bessel.tex +++ b/buch/chapters/050-differential/bessel.tex @@ -5,6 +5,7 @@ % \section{Bessel-Funktionen \label{buch:differntialgleichungen:section:bessel}} +\rhead{Bessel-Funktionen} Die Besselsche Differentialgleichung erlaubt Wellen mit zylindrischer Symmetrie und die Strömung in einem zylindrischen Rohr zu beschreiben. diff --git a/buch/chapters/050-differential/hypergeometrisch.tex b/buch/chapters/050-differential/hypergeometrisch.tex index de6752f..b93bc6e 100644 --- a/buch/chapters/050-differential/hypergeometrisch.tex +++ b/buch/chapters/050-differential/hypergeometrisch.tex @@ -5,6 +5,7 @@ % \section{Hypergeometrische Differentialgleichung \label{buch:differentialgleichungen:section:hypergeometrisch}} +\rhead{Hypergeometrische Differentialgleichung} Die hypergeometrische Funktion $\mathstrut_2F1(a,b;c;x)$ wurde in Abschnitt~\ref{buch:rekursion:section:hypergeometrische-funktion} als Potenzreihe mit sehr speziellen Koeffizienten, die sich aus @@ -377,6 +378,7 @@ x(1-x)\frac{d^2y}{dx^2} -ab y = 0 +\label{buch:differentialgleichungen:eqn:eulerhyper} \end{equation} hat die Lösung \[ @@ -628,20 +630,20 @@ Dies ist natürlich die Differentialgleichung der Exponentialfunktion. \subsubsection{Die Differentialgleichungen für $\mathstrut_1F_0$} In diesen Fälle gibt es nur jeweils einen einzigen Operator anzuwenden. -Wir betrachten zunächst den Fall $w(z) = \mathstrut_1F_0(\alpha; z)$ +Wir betrachten zunächst den Fall $w(z) = \mathstrut_1F_0(a; z)$ und finden direkt die Differentialgleichung \begin{align*} -\biggl(z\frac{d}{dz}+\alpha\biggr)w +\biggl(z\frac{d}{dz}+a\biggr)w &= \frac{d}{dz}w \\ -zw'+\alpha w +zw'+a w &= w' \\ (1-z)w' &= -\alpha w. +a w. \end{align*} \begin{beispiel} @@ -684,29 +686,33 @@ sie erfüllt also die genannte Differentialgleichung. % \subsubsection{Die Differentialgleichungen für $\mathstrut_0F_1$} Für die Funktion $\mathstrut_0F_1$ setzen wir -$w(z)=\mathstrut_0F_1(;\beta;z)$. +$w(z)=\mathstrut_0F_1(;b;z)$. In diesem Fall gibt es keine Operatoren $D_{a_i}$ anzuwenden, die linke Seite der Differentialgleichung ist also einfach die Funktion $w$. -Für die rechte Seite ist der Operator $D_{\beta-1}$ anzuwenden, was auf +Für die rechte Seite ist der Operator $D_{b-1}$ anzuwenden, was auf die Differentialgleichung -\begin{align*} +\begin{align} w &= \frac{d}{dz} -\biggl(z\frac{d}{dz}+\beta -1\biggr)w +\biggl(z\frac{d}{dz}+b -1\biggr)w +\notag \\ w &= -\frac{d}{dz}(zw'+\beta w - w) +\frac{d}{dz}(zw'+b w - w) +\notag \\ w &= -zw''+w'+\beta w' -w' +zw''+w'+b w' -w' +\notag \\ 0 &= -zw''+\beta w' - w -\end{align*} +zw''+b w' - w +\label{buch:differentialgleichungen:0f1:dgl} +\end{align} führt. \begin{beispiel} @@ -787,90 +793,100 @@ tatsächlich die Kosinus-Funktion als Lösung hat. % Die Differentialgleichung für 1F1 % \subsubsection{Die Differentialgleichung für $\mathstrut_1F_1$} -Wir setzen wieder $w(z) = \mathstrut_1F_1(\alpha;\beta;z)$. -Es sind die Operatoren $D_\alpha$ und $D_{\beta-1}$ anzuwenden. +Wir setzen wieder $w(z) = \mathstrut_1F_1(a;b;z)$. +Es sind die Operatoren $D_a$ und $D_{b-1}$ anzuwenden. Es ergibt sich die Differentialgleichung \begin{align*} -\biggl(z\frac{d}{dz}+\alpha\biggr)w +\biggl(z\frac{d}{dz}+a\biggr)w &= -\frac{d}{dz}\biggl(z\frac{d}{dz} +\beta-1\biggr)w +\frac{d}{dz}\biggl(z\frac{d}{dz} +b-1\biggr)w \\ -zw'+\alpha w +zw'+a w &= \frac{d}{dz} -(zw'+\beta w - w) +(zw'+b w - w) \\ -zw'+\alpha w +zw'+a w &= -zw'' +w'+\beta w' - w' +zw'' +w'+b w' - w' \\ 0 &= -zw'' + (\beta - z)w' - \alpha w. +zw'' + (b - z)w' - a w. \end{align*} % % Die hypergeometrische Differentialgleichung für 2F1 % \subsubsection{Die Differentialgleichung für $\mathstrut_2F_1$} -Für die hypergeometrische Funktion $\mathstrut_2F_1(\alpha,\beta;\gamma;z)$ +Für die hypergeometrische Funktion $\mathstrut_2F_1(a,b;c;z)$ ist die Difrentialgleichung von der Form \[ -\biggl(z\frac{d}{dz} + \alpha\biggr) -\biggl(z\frac{d}{dz} + \beta\biggr)w +\biggl(z\frac{d}{dz} + a\biggr) +\biggl(z\frac{d}{dz} + b\biggr)w = \frac{d}{dz} -\biggl(z\frac{d}{dz}+\gamma -1\biggr) +\biggl(z\frac{d}{dz}+c -1\biggr) w. \] Durchführen der Ableitungen auf beiden Seiten ergibt für die linke Seite \begin{align*} -\biggl(z\frac{d}{dz} + \alpha\biggr) -\biggl(z\frac{d}{dz} + \beta\biggr)w +\biggl(z\frac{d}{dz} + a\biggr) +\biggl(z\frac{d}{dz} + b\biggr)w &= -\biggl(z\frac{d}{dz} + \alpha\biggr) -(zw'+\beta w) +\biggl(z\frac{d}{dz} + a\biggr) +(zw'+b w) \\ &= -z^2w'' + zw' + \beta zw' + \alpha(zw'+\beta w) +z^2w'' + zw' + b zw' + a(zw'+b w) \\ &= -z^2w'' + (1+\alpha+\beta )zw' + \beta\alpha w +z^2w'' + (1+a+b )zw' + ba w \intertext{und die rechte Seite} -\frac{d}{dz}\biggl(z\frac{d}{dz}+\gamma-1\biggr)w +\frac{d}{dz}\biggl(z\frac{d}{dz}+c-1\biggr)w &= -\frac{d}{dz}(zw'+\gamma w-w) +\frac{d}{dz}(zw'+c w-w) \\ &= -zw''+w'+\gamma w' - w' +zw''+w'+c w' - w' \\ &= -zw'' +\gamma w'. +zw'' +c w'. \end{align*} Durch Gleichsetzen ergibt sich jetzt \begin{align*} -z^2w'' + (1+\alpha+\beta )zw' + \alpha\beta w +z^2w'' + (1+a+b )zw' + ab w &= -zw'' +\gamma w' +zw'' +c w' \\ 0 &= z(1-z)w'' + -(\gamma-z(1+\alpha+\beta))w' +(c-z(1+a+b))w' - -\alpha\beta +ab w \end{align*} Dies ist die früher definierte hypergeometrische Differentialgleichung. -\subsubsection{Differentialgleichungen für $w(x^2)$ und $xw(x^2)} - - % % Gerade und ungerade Funktionen % -\subsubsection{Gerade und ungerade Funktionen} +\subsection{Differentialgleichung für +$x^\varrho\cdot\mathstrut_pF_q(a_i;b_j;sx^\nu)$} +In verschiedenen Beispielen ist gezeigt worden, wie sich +wohlbekannte Funktionen durch hypergeometrische Funktionen +ausdrücken lassen. +Aus der Differentialgleichung der hypergeometrischen Funktionen +muss sich daher auch eine Differentialgleichung für die +gesuchten Funktionen ergeben. +Zum Beispiel lassen sich die Besselfunktionen durch hypergeometrische +Funktionen des Argumentes $-x^2/4$ schreiben, es muss also auch +möglich sein, die besselsche Differentialgleichung wieder aus +der eulerschen hypergeometrischen Differentialgleichung zu gewinnen. + +\subsubsection{Gerade und ungerade} Hypergeometrische Funktionen, deren Reihe mehr als einen Term enthalten, enthalten immer mindestens eine gerade und eine ungerade Potenz der unabhängigen Variable. @@ -883,6 +899,8 @@ Lösung $\cos x$ und die ungerade Lösung $\sin x$. Auch die Differentialgleichung $y''-y=0$ hat eine gerade Lösung, $\cosh x$, und eine ungerade Lösung, $\sinh x$. +\subsubsection{Symmetrien der eulerschen hypergeometrischen +Differentialgleichung} Hat die hypergeometrische Differentialgleichung gerade und ungerade Lösungen? Wenn es eine gerade Lösung $y(x)$ gibt, dann sollte die Substitution @@ -911,6 +929,7 @@ y(x) &= \frac{1}{-c+1}x^{-c+1} Dies zeigt, dass die hypergeometrische Differentialgleichung im allgemeinen keine geraden oder ungeraden Lösungen hat. +\subsubsection{Zusammengesetzte Funktionen} Die gerade oder ungeraden Funktionen, die in früheren Beispielen als hypergeometrische Funktionen dargestellt wurden, konnten also nicht Lösungen der hypergeometrische Differentialgleichung sein. @@ -927,22 +946,523 @@ g(x) \qquad\text{und}\qquad u(x) = -x\, +x\cdot \mathstrut_pF_q\biggl(\begin{matrix}a_1,\dots,a_p\\b_1,\dots,b_1\end{matrix};x^2\biggr). \] +Viele wohlbekannte Funktionen $f(x)$ können aus einer hypergeometrischen +Funktion $w(z)$ als $g(x)=w(\pm x^2)$ oder $u(x)=xw(\pm x^2)$ erhalten +werden. +Für die hypergeometrische Funktion $w(z)$ ist eine definierende +Differentialgleichung bekannt. +Im Folgenden soll daraus eine Differentialgleichung für $f(x)$ +abgeleitet werden. + +Der Fall $w(-x^2)$ könnte natürlich auch durch Verwendung eines +imaginären Arguments wie in $w((ix)^2)$ auf den Fall $w(x^2)$ +zurückgeführt werden. +Um aber die Differentialgleichungen reell zu belassen, schreiben +wir $g(x)=w(sx^2)$ für eine gerade Funktion beziehungsweise +$u(x)=xw(sx^2)$ für eine ungerade Funktion. +Der Faktor $s$ kann ausserdem dazu verwendet werden, das Argument +zu skalieren, wie es zum Beispiel in der Darstellung der +Kosinus-Funktion als $\cos x = \mathstrut_0F_1(;\frac12;-\frac{x^2}4)$ +nötig ist. + +Um die Differentialgleichung für $g$ oder $u$ zu finden, berechnen wir die +Ableitungen von $g$ und $u$, drücken die Ableitungen von $w$ durch die +Ableitungen $g$ und $u$ aus und setzen sie in die Differentialgleichung +für $w$ ein. +Wir wollen dies im Folgenden nur für ein paar Beispiele niedrigerer +Ordnung tun. + +% +% Differentialgleichungen für g(x)=w(sx^2) +% +\subsubsection{Differentialgleichungen für $g(x)=w(sx^2)$ für $w=\mathstrut_2F_1$} +Die Ableitungen von $g$ sind +\[ +\begin{linsys}{3} + g(x)&=&w(sx^2)& & & & \\ + g'(x)&=& & &2sxw'(sx^2)& & \\ +g''(x)&=& & & 2sw'(sx^2)&+&4s^2x^2w''(sx^2)\\ +\end{linsys} +\] +Dies sind lineare Gleichungssysteme für die Ableitungen +von $w$, die wir nach $w'$ und $w''$ auflösen können. +Dies wird einfacher, wenn wir das Gleichungssysteme in +Matrixschreibweise darstellen mit der Matrix +\begin{equation*} +\begin{pmatrix} + g(x)\\ + g'(x)\\ +g''(x) +\end{pmatrix} += +\underbrace{ +\begin{pmatrix} +1& 0& 0\\ +0& 2sx& 0\\ +0& 2s& 4s^2x^2 +\end{pmatrix} +}_{\displaystyle = A_g} +\begin{pmatrix} +w(sx^2)\\ +w'(sx^2)\\ +w''(sx^2) +\end{pmatrix} +\qquad\Rightarrow\qquad +A_g^{-1} += +\begin{pmatrix} +1& 0 & 0 \\ +0& \frac{1}{2sx} & 0 \\ +0& -\frac{1}{4s^2x^3} & \frac{1}{4s^2x^2} +\end{pmatrix} +\end{equation*} +Damit lassen sich jetzt die Ableitungen von $w$ ausdrücken: +\begin{align*} + w(sx^2) &= g(x) \\ + w'(sx^2) &= \frac{1}{2sx} g'(x) \\ +w''(sx^2) &= \frac{1}{4s^2x^3}g'(x) + \frac{4}{4s^2x^2} g''(x) +\end{align*} +Durch Einsetzen dieser Ausdrücke in die eulersche hypergeometrische +Differentialgleichung~\eqref{buch:differentialgleichungen:eqn:eulerhyper} +%\[ +%z(1-z) \frac{d^2y}{dz^2} + (c+(a+b+1)z)\frac{dy}{dz} - abz = 0 +%\] +finden wir jetzt die Differentialgleichungen für $g(x)$ +wie folgt: +\begin{equation} +x(1-sx^2)g''(x) ++ +(2c-1-(2a+2b+1)sx^2) g'(x) +-4absx g(x) += +0. +\label{buch:differential:hypergeometrisch:geradedgl} +\end{equation} + + % % % -\subsubsection{Hypergeometrische Funktionen von $x^2$ und $x^3$} -Die hypergeometrischen Funktionen $w(z)$ als Lösungen der hypergeometrischen -Differentialgleichungen sind Potenzreihen, in denen kein Koeffizient -verschwindet, sofern die Lösung nicht ein Polynom ist. -Die trigonometrischen Funktionen sind nicht von dieser Art. -Sie lassen sich als Funktionen von $x^2$ schreiben. -Wir untersuchen in diesem Abschnitt, wie sich eine Differentialgleichung -von $y(x) = x^lw(tx^k)$ aus der Differentialgleichung für $w(z)$ gewinnen -lässt. +\subsubsection{Differentialgleichung von $x^\varrho w(sx^\nu)$} +Die Methode der verallgemeinerten Potenzreihen zeigt, dass eine +gewöhnliche Potenzreihe, wie $\mathstrut_pF_q$ eine ist, +manchmal nicht reicht, um eine Lösung einer Differentialgleichung +darzustellen. +Die Besselsche Differentialgleichung ist von dieser Art. +Eine verallgemeinerte Potenzreihe erhält man mit Hilfe eines +zusätzlichen Faktors der Form $x^\varrho$. +Der Fall $\varrho=1$ deckt auch die früher vorgeschlagene +Funktion $u(x)=x\cdot \mathstrut_pF_q(sx^2)$ ab. +In den bisherigen Beispielen haben wir als Argument für eine +hypergeometrische Funktion $\mathstrut_pF_q$ einen Ausdruck +der Form $sx^2$ verwendet, was für die Beispiele gereicht hat, +aber zum Beispiel für die später untersuchten Airy-Funktionen +nicht genügt. +Daher soll jetzt für eine Funktion $f(x) += x^\varrho\cdot \mathstrut_pF_q(sx^\nu)$ eine Differentialgleichung +aus der Differentialgleichung der hypergeometrischen Funktion +abgeleitet werden. +Dabei soll der im vorangegangenen Abschnitt behandelte Fall +$\varrho=0$ und $\nu=2$ als Leitlinie dienen. +Wie vorhin beginnen wir damit, die Ableitungen von $f(x)$ zu +berechnen: +\[ +\begin{linsys}{4} +f(x) &=& x^\varrho w(sx^\nu) + & & + & & +\\ +f'(x) &=& \varrho x^{\varrho-1} w(sx^\nu) + &+& \nu s x^{\varrho+\nu-1} w'(sx^\nu) + & & +\\ +f''(x) &=& (\varrho-1)\varrho x^{\varrho-2} w(sx^\nu) + &+& \nu(2\varrho+\nu-1)sx^{\varrho+\nu-2} w'(sx^\nu) + &+& \nu^2 s^2 x^{\varrho+2\nu-2} w''(sx^\nu) +\end{linsys} +\] +Dies ist ein lineares Gleichungssystem, welches in Matrixform +geschrieben werden kann als +\[ +\begin{pmatrix} +f(x)\\ +f'(x)\\ +f''(x) +\end{pmatrix} += +% [ rho ] +% [ x ] +% [ ] +%(%o40) Col 1 = [ rho - 1 ] +% [ rho x ] +% [ ] +% [ 2 rho - 2 ] +% [ (rho - rho) x ] +% +% [ 0 ] +% [ ] +% [ rho + nu - 1 ] +% Col 2 = [ nu s x ] +% [ ] +% [ 2 rho + nu - 2 ] +% [ (2 nu rho + nu - nu) s x ] +% +% [ 0 ] +% [ ] +% Col 3 = [ 0 ] +% [ ] +% [ 2 2 rho + 2 nu - 2 ] +% [ nu s x ] +\underbrace{ +\begin{pmatrix} +x^\varrho + & 0 + & 0 \\ +\varrho x^{\varrho-1} + & \nu s x^{\varrho+\nu-2} + & 0 \\ +(\varrho-1)\varrho x^{\varrho-2} + & \nu(2\varrho+\nu-1)sx^{\varrho+\nu-2} + & \nu^2 s^2 x^{\varrho+2\nu-2} +\end{pmatrix} +}_{\displaystyle = A_f} += +\begin{pmatrix} +w(sx^\nu)\\ +w'(sx^\nu)\\ +w''(sx^\nu) +\end{pmatrix}. +\] +Die Inverse der Matrix $A_f$ ist +\[ +A_f^{-1} += +% [ 1 ] +% [ ---- ] +% [ rho ] +% [ x ] +% [ ] +% [ (- rho) - nu ] +% [ rho x ] +%(%o42) Col 1 = [ - ----------------- ] +% [ nu s ] +% [ ] +% [ 2 (- rho) - 2 nu ] +% [ (rho + nu rho) x ] +% [ ------------------------------- ] +% [ 2 2 ] +% [ nu s ] +% [ 0 ] +% [ ] +% [ (- rho) - nu + 1 ] +% [ x ] +% [ ----------------- ] +% Col 2 = [ nu s ] +% [ ] +% [ (- rho) - 2 nu + 1 ] +% [ (2 rho + nu - 1) x ] +% [ - ------------------------------------ ] +% [ 2 2 ] +% [ nu s ] +% [ 0 ] +% [ ] +% [ 0 ] +% [ ] +% Col 3 = [ (- rho) - 2 nu + 2 ] +% [ x ] +% [ ------------------- ] +% [ 2 2 ] +% [ nu s ] +\renewcommand{\arraystretch}{1.7} +\frac{1}{x^\varrho} +\begin{pmatrix} +\displaystyle 1 + & 0 + & 0 +\\ +\displaystyle-\frac{\varrho}{\nu s} x^{-\nu} + &\displaystyle \frac{1}{\nu s} x^{-\nu+1} + & 0 +\\ +\displaystyle\frac{\varrho^2+\nu\varrho}{\nu^2s^2}x^{-2\nu} + &\displaystyle -\frac{2\varrho+\nu-1}{\nu^2s^2} x^{-2\nu+1} + &\displaystyle \frac{1}{\nu^2s^2} x^{-2\nu+2} +\end{pmatrix} +\] +Damit kann man jetzt die Funktion $w(sx^\nu)$ und die Ableitungen +$w'(sx^\nu)$ und $w''(sx^\nu)$ durch $f$ und die Ableitungen davon +ausdrücken als +\begin{equation} +\renewcommand{\arraystretch}{2.3} +\begin{linsys}{4} +w(sx^\nu) + &=& \displaystyle \frac{1}{x^\varrho} f(x) + & & + & & +\\ +w'(sx^\nu) + &=& \displaystyle -\frac{\varrho}{\nu s} x^{-\varrho-\nu} f(x) + &+& \displaystyle \frac{1}{\nu s}x^{-\varrho-\nu+1} + & & +\\ +w''(sx^\nu) + &=& \displaystyle \frac{\varrho^2+\nu\varrho}{\nu^2 s^2}x^{-\varrho-2\nu} + &-& \displaystyle \frac{2\varrho+\nu-1}{\nu^2s^2} x^{-\varrho-2\nu+1} + &+& \displaystyle \frac{1}{\nu^2 s^2} x^{-\varrho-2\nu+2}. +\end{linsys} +\label{buch:differentialgleichungen:hypergeometrisch:wsubst} +\end{equation} +Einsetzen in die Differentialgleichung der hypergeometrischen Funktion +$\mathstrut_2F_1$ liefert die Differentialgleichung +\begin{equation} +\begin{aligned} +sx^4(x^2-1) f'' +%(%i50) ratsimp(subst(0,F,subst(0,Fp,ef))/Fpp) +%(%o50) -x^((-rho)-2*nu)*(s*x^6-x^4) +%(%i51) ratsimp(subst(0,F,subst(0,Fpp,ef))/Fp) +%\\ +&-( + ((-2\varrho-\nu+1) s x^5 +x^\nu(b+a+1)\nu s x^3-c\nu x) + + + (2\varrho+\nu-1)x^3 +)f' +\\ +%(%o51) -x^((-rho)-2*nu)*(((-2*rho)-nu+1)*s*x^5+x^nu*((b+a+1)*nu*s*x^3-c*nu*x) +% +(2*rho+nu-1)*x^3) +%(%i52) ratsimp(subst(0,Fp,subst(0,Fpp,ef))/F) +%(%o52) -x^((-rho)-2*nu)*(a*b*nu^2*s*x^(2*nu)+(rho^2+nu*rho)*s*x^4 +% +x^nu +% *(((-b)-a-1)*nu*rho*s*x^2 +% +c*nu*rho) +% +((-rho^2)-nu*rho)*x^2) +&-( + ab\nu^2 sx^{2\nu} + + \varrho(\varrho+\nu)sx^4 + + x^\nu((-b-a-1)\nu\varrho s x^2 + c\nu\varrho) + - \varrho(\varrho+\nu)x^2 +)f +=0 +\end{aligned} +\label{buch:differentialgleichungen:hypergeometrisch:2f1dgl} +\end{equation} +für die Funktion +\[ +f(x)=x^\varrho \cdot\mathstrut_2F_1\biggl( +\begin{matrix}a,b\\c\end{matrix}; sx^\nu +\biggr). +\] +Die Differentialgleichung +\eqref{buch:differentialgleichungen:hypergeometrisch:2f1dgl} +ist etwas unübersichtlich, daher soll sie in einem Beispiel illustriert +werden. +\begin{beispiel} +Früher in Aufgabe \ref{401} auf Seite \pageref{401} +wurde gezeigt, dass +\[ +\arcsin x = x\,\mathstrut_2F_1\biggl( +\begin{matrix}\frac12,\frac12\\\frac32\end{matrix};x^2 +\biggr) += +x^\nu\cdot \mathstrut_2F_1\biggl( +\begin{matrix}a,b\\c\end{matrix}; +sx^\nu +\biggr) +\] +ist. +Die Arkus-Sinus-Funktion ist daher Lösung der Differentialgleichung +\eqref{buch:differentialgleichungen:hypergeometrisch:2f1dgl} +mit +\[ +\varrho=1,\quad +a=\frac12,\quad +b=\frac12,\quad +c=\frac32,\quad +s=1,\quad +\nu=2, +\] +also +\begin{equation} +-\frac{x^2-1}{x}f'' +-f' += +0 +\qquad\Rightarrow\qquad +(1-x^2)f''=xf'. +\label{buch:differentialgleichungen:hypergeometrisch:beispiel:arcsindgl} +\end{equation} +Tatsächlich ist +\[ +\frac{d}{dx}\arcsin x += +\frac{1}{(1-x^2)^{\frac12}} +\qquad\text{und}\qquad +\frac{d^2}{dx^2} \arcsin x += +\frac{x}{(1-x^2)^{\frac32}}, +\] +und nach Einsetzen in die Differentialgleichung +\[ +(1-x^2) +\cdot +\frac{x}{(1-x^2)^{\frac32}} +- +x +\cdot +\frac{1}{(1-x^2)^{\frac12}} += +0. +\] +Die Arkus-Sinus-Funktion ist also tatsächlich eine Lösung der +Differentialgleichung +\eqref{buch:differentialgleichungen:hypergeometrisch:beispiel:arcsindgl}. +\end{beispiel} + +% +% +% +\subsubsection{Differentialgleichung für Funktionen, die aus $\mathstrut_0F_1$ zusammengesetzt sind} +Die Substitutionen +\eqref{buch:differentialgleichungen:hypergeometrisch:wsubst} +angewendet auf die Differentialgleichung +\eqref{buch:differentialgleichungen:0f1:dgl} +der Funktion $\mathstrut_0F_1$ +liefert +% +%(%i60) ratsimp(subst(0,F,subst(0,Fp,e0f1))/Fpp) +%(%o60) x^2 +%(%i61) ratsimp(subst(0,F,subst(0,Fpp,e0f1))/Fp) +%(%o61) ((-2*rho)+(beta-1)*nu+1)*x +%(%i62) ratsimp(subst(0,Fp,subst(0,Fpp,e0f1))/F) +%(%o62) (-nu^2*s*x^nu)+rho^2+(1-beta)*nu*rho +\begin{equation} +x^2f'' ++ +(-2\varrho+(\beta-1)\nu+1)xf' ++ +(-\nu^2sx^\nu + \varrho^2 -(\beta-1)\nu\varrho)f += +0. +\label{buch:differentialgleichungen:0F1:dgl} +\end{equation} +Die nächsten zwei Abschnitte sollen zeigen, wie sich daraus für die +Bessel-Funktionen wie auch die Airy-Funktionen, die sich durch +$\mathstrut_0F_1$ ausdrücken, die Besselsche und die Airysche +Differentialgleichung wiedergewonnen werden kann. + +% +% Besselsche Differentialgleichung +% +\subsubsection{Besselsche Differentialgleichung} +Die Besselfunktionen lassen sich in der Form +\begin{equation} +J_\alpha(x) += +\frac{(x/2)^\alpha}{\Gamma(\alpha+1)} \, +\mathstrut_0F_1\biggl( +\begin{matrix}\text{---}\\\alpha+1\end{matrix};-\frac14x^2 +\biggr) += +\frac{1}{2^\alpha\Gamma(\alpha+1)} +x^\varrho\cdot +\mathstrut_0F_1\biggl( +\begin{matrix}\text{---}\\b\end{matrix};sx^\nu +\biggr) +\label{buch:differentialgleichungen:0f1:besselfunktion} +\end{equation} +schreiben. +Somit sollte sich aus der +Differentialgleichung~\eqref{buch:differentialgleichungen:0f1:dgl} +der Funktion $\mathstrut_0F_1$ die Besselsche Differentialgleichung +\eqref{buch:differentialgleichungen:eqn:bessel} rekonstruieren lassen. +Dazu substituieren wir die aus +\eqref{buch:differentialgleichungen:0f1:besselfunktion} +abgelesenen Parameter +\[ +\varrho=\alpha,\quad\nu=2,\quad s=-\frac14,\quad b=\alpha+1 +\] +in \eqref{buch:differentialgleichungen:0F1:dgl} und erhalten +die Differentialgleichung +\begin{equation} +x^2y'' ++ +%(-2\alpha+2\alpha+1)xy' +xy' ++ +%(-4sx^2 + \alpha^2 -2\alpha^2)y +(x^2 - \alpha^2)y += +0. +\label{buch:differentialgleichungen:0F1:besseldgl} +\end{equation} +Dies ist tatsächlich die Besselsche Differentialgleichung. + +% +% Airy-Differentialgleichung +% +\subsubsection{Die Airy-Differentialgleichung} +Die in Aufgabe \ref{501} untersuchte +Airy-Differentialgleichung $y''-xy=0$ hat die Funktionen +\begin{align*} +y_1(x) +&= +\mathstrut_0F_1\biggl( +\begin{matrix}\text{---}\\\frac23\end{matrix};\frac{x^3}9 +\biggr) +=x^\varrho\cdot \mathstrut_0F_1\biggl( +\begin{matrix}\text{---}\\b\end{matrix};sx^\nu +\biggr) +&&\text{mit $\varrho=0$, $\nu=3$, $s=\frac19$, $b=\frac23$, } +\intertext{und} +y_2(x) +&= +x\cdot +\mathstrut_0F_1\biggl( +\begin{matrix}\text{---}\\\frac23\end{matrix};\frac{x^3}9 +\biggr) +=x^\varrho\cdot \mathstrut_0F_1\biggl( +\begin{matrix}\text{---}\\b\end{matrix};sx^\nu +\biggr) +&&\text{mit $\varrho=1$, $\nu=3$, $s=\frac19$, $b=\frac43$, } +\end{align*} +als Lösungen. +Die Differentialgleichung von $\mathstrut_0F_1$ sollte sich in diesem +Fall also auf die Airy-Differentialgleichung reduzieren lassen. + +Bei der Substition der Parameter in die Differentialgleichung +\eqref{buch:differentialgleichungen:0F1:dgl} beachten wird, dass +die beiden möglichen Werte für $b$ auf $b-1=\pm\frac13$ +führen, mit dem positiven Zeichen für den zweiten Fall, in dem $\varrho=1$ +ist. +So ergibt sich die Differentialgleichung +\begin{align*} +x^2y'' ++ +(-2\varrho\pm\frac13\cdot 3+1)xy' ++ +(-x^\nu + \varrho^2 \mp\frac13\cdot 3\varrho)y +&= +0 +\\ +x^2y'' ++ +(-2\varrho\pm1+1)xy' ++ +(-x^3 + \varrho^2 \mp\varrho)y +&= +0 +\\ +x^2y'' +- +x^3y +&= +0 +\qquad\Rightarrow\qquad y''-xy=0. +\end{align*} +Dies ist wie erwartet die Airy-Differentialgleichung. diff --git a/buch/chapters/050-differential/potenzreihenmethode.tex b/buch/chapters/050-differential/potenzreihenmethode.tex index 127f4d7..2d95fb2 100644 --- a/buch/chapters/050-differential/potenzreihenmethode.tex +++ b/buch/chapters/050-differential/potenzreihenmethode.tex @@ -5,6 +5,7 @@ % \section{Potenzreihenmethode \label{buch:differentialgleichungen:section:potenzreihenmethode}} +\rhead{Potenzreihenmethode} Die Potenzreihenmethode versucht die Lösung einer gewöhnlichen Differentialgleichung als Potenzreihe um die Anfangsbedingung zu entwickeln. diff --git a/buch/chapters/050-differential/w.maxima b/buch/chapters/050-differential/w.maxima new file mode 100644 index 0000000..24dd91d --- /dev/null +++ b/buch/chapters/050-differential/w.maxima @@ -0,0 +1,101 @@ + +z: s*x^2; +equation: z*(1-z)*Ypp + (c-(a+b+1)*z)*Yp-a*b*Y; + +gradef(w(x), wp(x)); +gradef(wp(x), wpp(x)); + +g: w(z); +gp: diff(g, x); +gpp: diff(gp, x); + +Ag: matrix( +[ coeff(g, w(z)), coeff(g, wp(z)), coeff(g, wpp(z)) ], +[ coeff(gp, w(z)), coeff(gp, wp(z)), coeff(gp, wpp(z)) ], +[ coeff(gpp,w(z)), coeff(gpp,wp(z)), coeff(gpp,wpp(z)) ] +); +Agi: invert(Ag); +W: Agi.matrix([G],[Gp],[Gpp]); + +eg: subst(W[1,1], Y, equation); +eg: subst(W[2,1], Yp, eg); +eg: subst(W[3,1], Ypp, eg); +eg: 4*s*x*eg; +eg: expand(ratsimp(eg)); +ratsimp(subst(0, G, subst(0, Gp, eg)) / Gpp); +ratsimp(subst(0, G, subst(0, Gpp, eg)) / Gp); +ratsimp(subst(0, Gp, subst(0, Gpp, eg)) / G); + + +u: x*w(z); +up: diff(u, x); +upp: diff(up, x); + +Au: matrix( +[ coeff(u, w(z)), coeff(u, wp(z)), coeff(u, wpp(z)) ], +[ coeff(up, w(z)), coeff(up, wp(z)), coeff(up, wpp(z)) ], +[ coeff(upp,w(z)), coeff(upp,wp(z)), coeff(upp,wpp(z)) ] +); +Aui: invert(Au); +W: Aui.matrix([U],[Up],[Upp]); + +eu: subst(W[1,1], Y, equation); +eu: subst(W[2,1], Yp, eu); +eu: subst(W[3,1], Ypp, eu); +eu: 4*s*x^3*eu; +eu: expand(ratsimp(eu)); +display2d: false$ +ratsimp(subst(0, U, subst(0, Up, eu)) / Upp); +ratsimp(subst(0, U, subst(0, Upp, eu)) / Up); +ratsimp(subst(0, Up, subst(0, Upp, eu)) / U); + +display2d: true$ + +/* allgemeiner Fall, f(x) = x^nu w(s * x^rho) */ + +z: s*x^nu; +f: x^rho * w(z); +fp: diff(f, x); +fpp: diff(fp, x); + +Af: ratsimp(matrix( +[ coeff(f, w(z)), coeff(f, wp(z)), coeff(f, wpp(z)) ], +[ coeff(fp, w(z)), coeff(fp, wp(z)), coeff(fp, wpp(z)) ], +[ coeff(fpp,w(z)), coeff(fpp,wp(z)), coeff(fpp,wpp(z)) ] +)); +Afi: invert(Af); +Afi: ratsimp(Afi); +W: Afi.matrix([F],[Fp],[Fpp]); + +ef: subst(W[1,1], Y, equation); +ef: subst(W[2,1], Yp, ef); +ef: subst(W[3,1], Ypp, ef); + ef: s * nu^2 * ef; +ef: expand(ratsimp(ef)); +display2d: true$ +c2: ratsimp(subst(0, F, subst(0, Fp, ef)) / Fpp); +c1: ratsimp(subst(0, F, subst(0, Fpp, ef)) / Fp); +c0: ratsimp(subst(0, Fp, subst(0, Fpp, ef)) / F); + +ratsimp(subst(1, s, subst(3/2, c, subst(1/2, b, subst(1/2, a, + subst(2, nu, subst(1, rho, c2))))))); +ratsimp(subst(1, s, subst(3/2, c, subst(1/2, b, subst(1/2, a, + subst(2, nu, subst(1, rho, c1))))))); +ratsimp(subst(1, s, subst(3/2, c, subst(1/2, b, subst(1/2, a, + subst(2, nu, subst(1, rho, c0))))))); + + +/* Differentialgleichung von 0F1 */ + +display2d: true$ +equation0f1: z*Ypp + beta*Yp - Y; +e0f1: subst(W[1,1], Y, equation0f1); +e0f1: subst(W[2,1], Yp, e0f1); +e0f1: subst(W[3,1], Ypp, e0f1); +e0f1: s*nu^2*x^(rho+nu) * e0f1; +e0f1: expand(ratsimp(e0f1)); +display2d: false$ +ratsimp(subst(0, F, subst(0, Fp, e0f1)) / Fpp); +ratsimp(subst(0, F, subst(0, Fpp, e0f1)) / Fp); +ratsimp(subst(0, Fp, subst(0, Fpp, e0f1)) / F); + -- cgit v1.2.1