From b95aef80656a0d476fcc94f2e344fc349cb5f087 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Andreas=20M=C3=BCller?= Date: Thu, 25 Nov 2021 08:38:34 +0100 Subject: add new structure --- buch/chapters/050-differential/beispiele.tex | 105 +++++++++- buch/chapters/050-differential/chapter.tex | 49 +++++ .../050-differential/potenzreihenmethode.tex | 233 +++++++++++++++++++++ 3 files changed, 382 insertions(+), 5 deletions(-) (limited to 'buch/chapters/050-differential') diff --git a/buch/chapters/050-differential/beispiele.tex b/buch/chapters/050-differential/beispiele.tex index 1182e7c..c0a57d8 100644 --- a/buch/chapters/050-differential/beispiele.tex +++ b/buch/chapters/050-differential/beispiele.tex @@ -6,12 +6,107 @@ \section{Beispiele \label{buch:differentialgleichungen:section:beispiele}} \rhead{Beispiele} +Viele der bisher betrachteten speziellen Funktionen können +durch gewöhnliche Differentialgleichungen charakterisiert werden, +als deren Lösungen sie auftreten. -\subsection{Exponentialfunktion +\subsection{Potenzen und Wurzeln +\label{buch:differentialgleichungen:subsection:potenzen-und-wurzeln}} +Die Potenzfunktionen und die zugehörigen Wurzeln als die ältesten +speziellen Funktionen bieten bereits eine erste kleine Schwierigkeit. +Die Differentialgleichung, die man aus einem naiven Ansatz ableitet, +ist singulär. + +\subsubsection{Differentialgleichung in $(0,\infty)$} +Die Ableitung einer Potenzfunktion $x\mapsto y(x)=x^\alpha$ ist +\[ +y'(x) = +\begin{cases} +\alpha x^{\alpha-1} &\qquad \alpha\ne -1\\ +\log x&\qquad\text{sonst} +\end{cases} +\] +Im Folgenden wollen wir uns auf den Fall $\alpha\ne -1$ konzentrieren. +Die Ableitungsoperation läuft in diesem Fall darauf hinaus, dass der +Grad um $1$ reduziert wird. +Dies könnte man mit einem Faktor $x$ komponsieren. +Wir fragen daher nach der allgmeinen Lösung der linearen +Differentialgleichung der Form +\begin{equation} +xy' = \alpha y. +\label{buch:differentialgleichungen:eqn:wurzeldgl} +\end{equation} +Diese Gleichung ist separierbar, die Separation von $x$ und $y$ liefert +die Integrale +\[ +\int \frac{dy}{y} = \alpha \int \frac{dx}{x} + C. +\] +Die Durchführunge der Integration liefert +\[ +\log |y| = \alpha \log|x| + C. +\] +Wendet man die Exponentialfunktion an, erhält man wieder +\[ +y = Dx^\alpha,\quad D=\exp C. +\] + +Die Differentialgleichung~\eqref{buch:differentialgleichungen:eqn:wurzeldgl} +hat aber eine schwerwiegenden Mangel. +Ihre explizite Form lautet +\begin{equation} +y' = \frac{\alpha}{x}\cdot y. +\label{buch:differentialgleichungen:eqn:wurzelsing} +\end{equation} +Dies ist zwar durchaus eine lineare Differentialgleichung erster Ordnung, +aber der Koeffiziente $\alpha/x$ wächst für $x\to 0$ über alle Grenzen. +Man kann daher den Wert der Potenzfunktion im Nullpunkt gar nicht aus der +Differentialgleichung erhalten, es ist dazu mindestens noch ein Grenzübergang +$x\to 0+$ nötig. + +\subsubsection{Differentialgleichung in der Nähe von $x=1$} +Um dem Problem des singulären Koeffizienten der +Differentialgleichung~\eqref{buch:differentialgleichungen:eqn:wurzelsing} +aus dem Weg zu gehen, verwenden wir die Variable $t$ mit $x=1+t$ und +versuchen eine Differentialgleichung für die Potenzfunktion +$(1+t)^\alpha$ zu finden. +Es gilt natürlich +\begin{equation} +\frac{d}{dt} (1+t)^\alpha += +\alpha (1+t)^{\alpha-1} +\qquad\Rightarrow\qquad +(1+t) \dot{y} = \alpha y. +\label{buch:differentialgleichungen:eqn:wurzeldgl1} +\end{equation} +Diese Differentialgleichung kann natürlich auch wieder mit Separation +gelöst werden, es ist +\begin{equation} +\int +\frac{dy}{y} += +\alpha +\int +\frac{dt}{1+t} ++ +C +\qquad\Rightarrow\qquad +\log|y| = \alpha \log|1+t| + C +\label{buch:differentialgleichungen:eqn:wurzeldgl1loesung} +\end{equation} +und daraus die Potenzfunktion +\[ +y=D(1+t)^\alpha +\] +wie vorhin. +Der Vorteil der +Form~\eqref{buch:differentialgleichungen:eqn:wurzeldgl1} +wird sich später bei dem Versuch zeigen, die Fuktion $y(t)$ +direkt als Potenzreihenlösung der Differentialgleichung zu finden. + + +\subsection{Exponentialfunktion und ihre Varianten \label{buch:differentialgleichungen:subsection:exponentialfunktion}} -\subsection{Trigonometrische Funktionen -\label{buch:differentialgleichungen:subsection:trigonometrisch}} +\subsubsection{Lineare Differentialgleichung erster Ordnung mit konstanten Koeffizienten} -\subsection{Hyperbelfunktionen -\label{buch:differentialgleichungen:subsection:hyperbelfunktionen}} +\subsubsection{Lineare Differentialgleichungen zweiter Ordnung mit konstanten Koeffizienten} diff --git a/buch/chapters/050-differential/chapter.tex b/buch/chapters/050-differential/chapter.tex index 1cfc1dd..f8b0dc3 100644 --- a/buch/chapters/050-differential/chapter.tex +++ b/buch/chapters/050-differential/chapter.tex @@ -9,6 +9,55 @@ \label{buch:chapter:differential}} \lhead{Differentialgleichungen} \rhead{} +Allgemeine Sätze über die Existenz und Eindeutigkeit der Lösungen +gewöhnlicher Differentialgleichungen garantieren für fast jeder +einigermassen vernünftige Gleichung mindestens für kurze Zeit +eine eindeutige Lösung für fast jede Anfangsbedingung. +Die Konstruktion solcher Lösungen stellt sich jedoch als deutlich +schwieriger heraus. + +Für einzelne Kategorien von Differentialgleichungen sind +gut funktionierende Lösungsverfahren gefunden worden, zum Beispiel +für lineare Differentialgleichungen mit konstanten Koeffizienten. +Damit konnten auch Gleichungen gelöst werden, die sich zum Beispiel +durch eine Variablentransformation auf eine lineare Differentialgleichung +mit konstanten Koeffizienten reduzieren lassen, wie die Eulersche +Differentialgleichung. + +Die Methode der Separation der Variablen liefert führt die Lösung +einer Differentialgleichung erster Ordnung auf die Bestimmung +zweier Stammfunktionen und deren Invertierung zurück. +Dieses Verfahren ist jedoch nicht auf Vektordifferentialgleichungen +oder auf Differentialgleichungen höherer Ordnung verallgemeinerungsfähig. + +Daneben gibt es eine Reihe von ``Spezialfällen'' wie die +Clairaut-Differentialgleichung oder die damit verwandte +Lagrangesche Differentialgleichung, deren Lösung eine sehr +spezielle Form haben. + +Sehr viele Differentialgleichungen in den Anwendungen können aber +mit keinem der genannten Verfahren gelöst werden. +Hier bleibt nichts anderes übrig, als neue spezielle Funktionen +zu definieren, die Lösungen dieser Differentialgleichungen sind. +Dabei ist man bestrebt, möglichst universell einsetzbare Funktionen +zu definieren, die ein breites Anwendungsfeld haben. + +In den folgenden Abschnitten wird zunächst gezeigt, dass viele +der bereits bekannten speziellen Funktionen ebenfalls als Lösungen +gewöhnlicher Differentialgleichungen erhalten werden können. +Die numerische Lösung gewöhnlicher Differentialgleichungen ist +oft keine effizientes Vorgehen zur Bestimmung von einzelnen Werten, +daher wird in +Abschnitt~\ref{buch:differentialgleichungen:section:potenzreihenmethode} +eine universelle Methode vorgestellt, mit der eine Potenzreihenentwicklung +gefunden werden kann. +Eine Potenzreihendarstellung ermöglicht nicht nur die Berechnung +einzelner Werte, sondern auch beliebiger Ableitungen und die +analytische Untersuchung der Funktion mit den Methoden der +komplexen Analysis. +Als Beispiel für dieses Verfahren werden in +Abschnitt~\ref{buch:differntialgleichungen:section:bessel} +die Bessel-Funktionen erster Art vorgestellt. \input{chapters/050-differential/beispiele.tex} \input{chapters/050-differential/potenzreihenmethode.tex} diff --git a/buch/chapters/050-differential/potenzreihenmethode.tex b/buch/chapters/050-differential/potenzreihenmethode.tex index feafacb..6d30129 100644 --- a/buch/chapters/050-differential/potenzreihenmethode.tex +++ b/buch/chapters/050-differential/potenzreihenmethode.tex @@ -5,4 +5,237 @@ % \section{Potenzreihenmethode \label{buch:differentialgleichungen:section:potenzreihenmethode}} +Die Potenzreihenmethode versucht die Lösung einer gewöhnlichen +Differentialgleichung als Potenzreihe um die Anfangsbedingung zu +entwickeln. +Wir gehen in diesem Abschnitt von einer Differentialgleichung der +Form +\begin{equation} +a_n(x)y^{(n)}(x) ++ +a_{n-1}(x)y^{(n-1)}(x) ++ +\dots ++ +a_1(x)y'(x) ++ +a_0(x)y(x) += +f(x) +\label{buch:differentialgleichungen:eqn:potenzreihendgl} +\end{equation} +mit der Randbedingung $y(0)=y_0$ aus. +Schon im einfachsten Fall einer homogenen Differentialgleichung erster +Ordnung ergibt sich die Beziehung +\[ +a_1(x) y'(x) = a_0(x)y(x), +\] +wobei wir uns $y(x)$ und damit auch $y'(x)$ als Potenzreihe vorstellen. +Insbesondere ist +\[ +\frac{a_1(x)}{a_0(x)} = \frac{y(x)}{y'(x)} +\] +ein Quotient von Potenzreihen, den man natürlich wieder als +Potenzreihe schreiben kann. +Da es nur auf den Quotienten ankommt, kann man sich auf den Fall +beschränken, dass die Koeffizienten Potenzreihen sind. +Tatsächlich gilt der folgende sehr viel allgemeinere Satz von +Cauchy und Kowalevskaja: + +\begin{satz}[Cauchy-Kowalevskaja] +Eine partielle Differentialgleichung der Ordnung $k$ für eine +Funktion $u(x_1,\dots,x_n,t)=u(x,t)$ +in expliziter Form +\[ +\frac{\partial^k}{\partial t^k} += +G\biggl(x,t, +\frac{\partial^j\partial^\alpha}{\partial t^j\,\partial x^k} +\biggr) +\quad\text{mit $j\alpha$. +Für nicht ganzzahliges $\alpha$ sind alle Koeffizienten $a_k\ne 0$. + +Die Lösung der +Differentialgleichung~\eqref{buch:differentialgleichungen:eqn:wurzeldgl1} +ist daher die Reihe +\begin{equation} +(1+t)^\alpha += +\sum_{k=0}^\infty +\frac{\alpha(\alpha-1)\dots(\alpha-k+1)}{k!}\, t^k. +\label{buch:differentialgleichungen:eqn:newtonreihe} +\end{equation} +Für ganzzahliges $\alpha$ wird daraus die binomische Formel +\[ +(1+t)^\alpha += +\sum_{k=0}^\infty +\frac{\alpha(\alpha-1)\dots(\alpha-k+1)}{k!}\, t^k += +\sum_{k=0}^\alpha \binom{\alpha}{k} t^k. +\] + +% +% Lösung als hypergeometrische Riehe +% +\subsubsection{Lösung als hypergeometrische Funktion} +Die Newtonreihe verwendet ein absteigendes Produkt im Zähler. +Man kann sie aber in eine Form bringen, die besser zu den aufsteigenden +Produkten bringen, die wir im Zusammenhang mit der Gamma-Funktion +angetroffen und als Pochhammer-Symbole formalisiert haben. + +Eine hypergeometrische Funktion zeichnet sich dadurch aus, dass +die Quotienten aufeinanderfolgender Koeffizienten der Reihe rationale +Funktionen von $k$ sind. +Der Quotient ist +nach~\eqref{buch:differentialgleichungen:eqn:newtonreiherekursion} +\[ +\frac{a_{k+1}}{a_k} += +\frac{\alpha-k}{k+1}. +\] +Der Nenner wird nie $0$, aber das Zählerpolynom hat genau die Nullstelle +$-\alpha$. +Die Newtonsche Reihe muss sich daher als Wert der hypergeometrischen +Funktion $\mathstrut_1F_0$ schreiben lassen. + +Das Produkt im Zähler von $a_k$ hat $k$ Faktoren, indem wir jeden Faktor +mit $-1$ multiplizieren, erhalten wir +\begin{align*} +\alpha(\alpha-1)(\alpha-2)\dots(\alpha-k+1) +&= +(-\alpha)(-\alpha+1)(-\alpha+2)\dots(-\alpha+k-1) (-1)^k +\\ +&= +(-\alpha)_k (-1)^k. +\end{align*} +Indem wir den Faktor $-1$ in der Variablen absorbieren, erhalten +wir die Darstellung +\[ +(1+t)^\alpha += +\sum_{k=0}^\infty +(-\alpha)_k\frac{(-t)^k}{k!}. +\] +Damit haben wir den folgenden Satz gezeigt. + +\begin{satz} +Die Newtonsche Reihe für $(1-t)^\alpha$ ist der Wert +\[ +(1-t)^\alpha += +\sum_{k=0}^\infty (-\alpha)_k \frac{t^k}{k!} += +\mathstrut_1F_0(-\alpha;t) +\] +der hypergeometrischen Funktion $\mathstrut_1F_0$. +\end{satz} + -- cgit v1.2.1