Kurs:Mathematik (Osnabrück 2009-2011)/Teil II/Vorlesung 37

Vergleichskriterien mit Reihen

Lemma

Es sei ${}I=[1,\infty ]$ ein rechtsseitig unbeschränktes Intervall und sei

$f\colon I\longrightarrow \mathbb {R}$

eine stetige fallende Funktion mit ${}f(x)\geq 0$ für alle ${}x\in I$.

Dann existiert das uneigentliche Integral

$\int _{1}^{\infty }f(t)\,dt$

genau dann, wenn die Reihe

$\sum _{n=1}^{\infty }f(n)$

konvergiert.

Beweis

Wenn das uneigentliche Integral existiert, so betrachten wir die Abschätzung

${}\sum _{n=2}^{k}f(n)\leq \int _{1}^{k}f(t)\,dt\,,$

die darauf beruht, dass die linke Seite das Treppenintegral zu einer unteren Treppenfunktion für ${}f$ auf ${}[1,k]$ ist. Da die rechte Seite beschränkt ist, gilt dies auch für die linke Seite, sodass wegen ${}f(n)\geq 0$ die Reihe konvergiert.
Ist umgekehrt die Reihe konvergent, so betrachten wir die Abschätzung

${}\int _{1}^{k}f(t)\,dt\leq \sum _{n=1}^{k-1}f(n)\,,$

die gilt, da die rechte Seite das Treppenintegral zu einer oberen Treppenfunktion ist. Wegen ${}f(t)\geq 0$ ist die Integralfunktion ${}x\mapsto \int _{1}^{x}f(t)\,dt$ wachsend und beschränkt, da die rechte Seite wegen der Konvergenz der Reihe beschränkt ist. Daher besitzt die Integralfunktion für ${}x\mapsto \infty$ einen Grenzwert und das uneigentliche Integral existiert.

$\Box$

Beispiel

Die Funktion

$f\colon [1,\infty ]\longrightarrow \mathbb {R} ,\,t\longmapsto {\frac {1}{t}},$

ist streng fallend. Daher ist die Funktion ${}g$, die für ${}x$ mit ${}k\leq x<k+1$ (${}k\in \mathbb {N}$) durch ${}{\frac {1}{k}}$ definiert ist, eine „Majorante“ für ${}f$, also ${}g(t)\geq f(t)$. Auf jedem Intervall ${}[1,n]$ liefert ${}g$ eine obere Treppenfunktion zu ${}f$. Ebenso liefert die durch ${}{\frac {1}{k+1}}$ bei ${}k\leq x<k+1$ definierte Funktion ${}h$ eine untere Treppenfunktion für ${}f$. Daher gelten die Abschätzungen

${}\sum _{k=1}^{n-1}{\frac {1}{k}}\geq \int _{1}^{n}{\frac {1}{t}}\,dt\geq \sum _{k=1}^{n-1}{\frac {1}{k+1}}\,.$

Das Integral in der Mitte besitzt den Wert ${}\ln n$. Daraus ergibt sich mit Lemma 37.1 ein neuer Beweis, dass die harmonische Reihe divergiert.

Die blaue Fläche stellt die Eulersche Konstante dar, die Darstellung ist überhöht.

Die Differenz zwischen der linken und der rechten Summe ist ${}1-{\frac {1}{n}}$. Daher ist die Differenz

$\sum _{k=1}^{n-1}{\frac {1}{k}}-\ln n$

für jedes ${}n$ positiv, mit ${}n$ wachsend und nach oben beschränkt. Daher existiert für ${}n\rightarrow \infty$ der Limes, und dieser Limes ändert sich nicht, wenn man vorne in der Summe bis ${}n$ aufsummiert anstatt bis ${}n-1$. Wir setzen

${}\gamma :=\lim _{n\rightarrow \infty }{\left(\sum _{k=1}^{n}{\frac {1}{k}}-\ln n\right)}\,$

und nennen sie die eulersche Konstante (oder Mascheronische Konstante). Ihr numerischer Wert ist ungefähr

${}\gamma =0{,}5772156649\dotso \,.$

Es ist ein offenes mathematisches Problem, ob diese Zahl rational ist oder nicht.

Die Riemannsche Zeta-Funktion im Reellen

Nach Beispiel 36.8 existiert für ${}c<-1$ das uneigentliche Integral ${}\int _{1}^{\infty }t^{c}\,dt$, so dass aufgrund von Lemma 37.1, auch die Reihen ${}\sum _{n=1}^{\infty }n^{c}=\sum _{n=1}^{\infty }{\frac {1}{n^{-c}}}$ konvergieren. Daher ist die folgende Funktion wohldefiniert.

Definition

Die Riemannsche ${}\zeta$-Funktion ist für ${}s\in \mathbb {R}$ mit ${}s>1$ durch

${}\zeta (s)=\sum _{n=1}^{\infty }{\frac {1}{n^{s}}}\,$

definiert.

Diese Funktion lässt sich komplex fortsetzen und spielt eine wichtige Rolle in der Zahlentheorie.

Die Fakultätsfunktion

Die Fakultät einer natürlichen Zahl ${}n$ ist ${}n!=\cdot (n-1)\cdot (n-2)\cdots 3\cdot 2\cdot 1$. Dabei gilt die rekursive Beziehung ${}n!=n((n-1)!)$. Gibt es eine Möglichkeit, diese für die natürlichen Zahlen definierte Funktion auf ${}\mathbb {R} _{\geq 0}$ durch eine differenzierbare Funktion ${}f$ fortzusetzen? Ist es sogar möglich, dass dabei die Beziehung ${}f(x)=xf(x-1)$ für jedes ${}x$ gilt? Wir werden mit Hilfe von uneigentlichen Integralen zeigen, dass dies in der Tat möglich ist.

Beispiel

Es sei ${}x>-1$. Wir betrachten die Funktion

$\mathbb {R} _{+}\longrightarrow \mathbb {R} ,\,t\longmapsto t^{x}e^{-t}.$

Wir behaupten, dass das uneigentliche Integral

$\int _{0}^{\infty }t^{x}e^{-t}\,dt$

existiert. Für den rechten Rand (also ${}\infty$) betrachten wir eine natürliche Zahl ${}n\geq x$. Da die Exponentialfunktion schneller wächst als jede Polynomfunktion (siehe Aufgabe 25.22), gibt es ein ${}a\in \mathbb {R} _{+}$ derart, dass ${}t^{n}e^{-{\frac {t}{2}}}\leq 1$ für alle ${}t\geq a$ gilt. Daher ist

${}{\begin{aligned}\int _{a}^{b}t^{x}e^{-t}\,dt&\leq \int _{a}^{b}t^{n}e^{-t}\,dt\\&=\int _{a}^{b}t^{n}e^{-{\frac {t}{2}}}e^{-{\frac {t}{2}}}\,dt\\&\leq \int _{a}^{b}e^{-{\frac {t}{2}}}\,dt\\&=2{\left(e^{-{\frac {a}{2}}}-e^{-{\frac {b}{2}}}\right)}\\&\leq 2e^{-{\frac {a}{2}}}.\end{aligned}}$

Für ${}b\rightarrow \infty$ wächst das linke Integral und ist durch ${}2e^{-{\frac {a}{2}}}$ beschränkt, sodass der Grenzwert existiert. Für das Verhalten am linken Rand (das nur bei ${}-1<x\leq 0$ problematisch ist) müssen wir wegen ${}e^{-t}\leq 1$ nach Lemma 36.6 nur ${}\int _{0}^{1}t^{x}\,dt$ betrachten. Eine Stammfunktion davon ist ${}{\frac {1}{x+1}}t^{x+1}$, deren Exponent positiv ist, sodass der Limes für ${}t\rightarrow 0$ existiert.

Das uneigentliche Integral

$\int _{0}^{\infty }t^{x}e^{-t}\,dt$

existiert also für ${}x\in \mathbb {R} ,\ x>-1$. Dies ist der Ausgangspunkt für die Definition der Fakultätsfunktion.

Definition

Für ${}x\in \mathbb {R}$, ${}x>-1$, heißt die Funktion

$x\longmapsto \operatorname {Fak} \,(x):=\int _{0}^{\infty }t^{x}e^{-t}\,dt$

die Fakultätsfunktion.

Die für ${}x>0$ durch

${}\Gamma (x):=\operatorname {Fak} \,(x-1)=\int _{0}^{\infty }t^{x-1}e^{-t}\,dt\,$

definierte Funktion heißt Gammafunktion, mit der häufiger gearbeitet wird. Mit der Fakultätsfunktion werden aber die Formeln etwas schöner und insbesondere wird der Zusammenhang zur Fakultät, der in der folgenden Aussage aufgezeigt wird, deutlicher.

Satz

Die Fakultätsfunktion besitzt die folgenden Eigenschaften.

1. Es ist ${}\operatorname {Fak} \,(x)=x\cdot \operatorname {Fak} \,(x-1)$ für ${}x>0$.
2. Es ist ${}\operatorname {Fak} \,(0)=1$.
3. Es ist ${}\operatorname {Fak} \,(n)=n!$ für natürliche Zahlen ${}n\in \mathbb {N}$.
4. Es ist ${}\operatorname {Fak} \,{\left(-{\frac {1}{2}}\right)}={\sqrt {\pi }}$.

Beweis

(1) Mittels partieller Integration ergibt sich (für reelle Zahlen ${}b\geq a>0$ bei fixiertem ${}x>0$)

${}{\begin{aligned}\int _{a}^{b}t^{x}e^{-t}\,dt&=-t^{x}e^{-t}|_{a}^{b}+\int _{a}^{b}xt^{x-1}e^{-t}\,dt\\&=-b^{x}e^{-b}+a^{x}e^{-a}+x\cdot \int _{a}^{b}t^{x-1}e^{-t}\,dt.\end{aligned}}$

Für ${}b\rightarrow \infty$ geht ${}b^{x}e^{-b}\rightarrow 0$ und für ${}a\rightarrow 0$ geht ${}a^{x}e^{-a}\rightarrow 0$ (da ${}x$ positiv ist). Wendet man auf beide Seiten diese Grenzwertprozesse an, so erhält man ${}\operatorname {Fak} \,(x)=x\operatorname {Fak} \,(x-1)$.
(2). Es ist

${}\operatorname {Fak} \,(0)=\int _{0}^{\infty }e^{-t}\,dt=-e^{-t}|_{0}^{\infty }=1\,.$

(3) folgt aus (1) und (2) durch Induktion.
(4). Es ist

${}\operatorname {Fak} \,{\left(-{\frac {1}{2}}\right)}=\int _{0}^{\infty }t^{-{\frac {1}{2}}}e^{-t}\,dt=2\int _{0}^{\infty }e^{-s^{2}}\,ds=\int _{-\infty }^{\infty }e^{-s^{2}}\,ds={\sqrt {\pi }}\,.$

Dies ergibt sich mit der Substitution ${}t=s^{2}$ und dem sogenannten Fehlerintegral.

$\Box$

Gewöhnliche Differentialgleichungen

Welche Bewegung vollzieht ein Löwenzahnfallschirmchen? Das Fallschirmchen lässt sich zu jedem Zeitpunkt von dem Wind tragen, der an der Stelle herrscht, wo es sich gerade befindet. Der Wind, seine Stärke und seine Richtung, hängt sowohl von der Zeit als auch vom Ort ab. Das bedeutet, dass hier ein gewisser „Rückkopplungsprozess“ vorliegt: Die bisherige Bewegung (also die Vergangenheit) bestimmt, wo sich das Fallschirmchen befindet und damit auch, welcher Wind auf es einwirkt und damit den weiteren Bewegungsablauf. Solche Bewegungsprozesse werden durch Differentialgleichungen beschrieben.

Differentialgleichungen sind ein fundamentaler Bestandteil der Mathematik und der Naturwissenschaften. Sie drücken eine Beziehung zwischen einer abhängigen Größe (häufig $y(t)$) und der Änderung dieser Größe ($y'(t)$) aus. Viele Gesetzmäßigkeiten in der Natur wie Bewegungsprozesse, Ablauf von chemischen Reaktionen, Wachstumsverhalten von Populationen werden durch Differentialgleichungen beschrieben. Hier besprechen wir nur solche Differentialgleichungen, die durch Integration gelöst werden können.

Die Vektorfelder werden im Folgenden nicht immer auf ganz ${}\mathbb {R} ^{2}$ definiert sein, sondern auf einer offenen Menge ${}U\subseteq \mathbb {R} ^{2}$. Dabei heißt ${}U$ offen, wenn es zu jedem Punkt ${}P\in U$ eine offene Ballumgebung ${}U(P,r)$ gibt, die ganz in ${}U$ liegt.

Definition

Es sei ${}U\subseteq \mathbb {R} ^{2}$ eine offene Teilmenge und es sei

$f\colon U\longrightarrow \mathbb {R} ,\,(t,y)\longmapsto f(t,y),$

eine Funktion. Dann nennt man

${}y'=f(t,y)\,$

die (gewöhnliche) Differentialgleichung zu ${}f$ (oder zum Vektorfeld oder zum Richtungsfeld ${}f$).

Dabei ist ${}y'=f(t,y)$ erstmal nur ein formaler Ausdruck, dem wir aber sofort eine inhaltliche Interpretation geben. Das ${}y$ soll eine Funktion in einer Variablen repräsentieren und ${}y'$ ihre Ableitung. Dies wird präzisiert durch den Begriff der Lösung einer Differentialgleichung.

Definition

Es sei ${}U\subseteq \mathbb {R} ^{2}$ eine offene Teilmenge und es sei

$f\colon U\longrightarrow \mathbb {R} ,\,(t,y)\longmapsto f(t,y),$

eine Funktion. Zur gewöhnlichen Differentialgleichung

${}y'=f(t,y)\,$

heißt eine Funktion

$y\colon I\longrightarrow \mathbb {R} ,\,t\longmapsto y(t),$

auf einem (mehrpunktigen) Intervall ${}I\subseteq \mathbb {R}$ eine Lösung der Differentialgleichung, wenn folgende Eigenschaften erfüllt sind.

1. Es ist ${}(t,y(t))\in U$ für alle ${}t\in I$.
2. Die Funktion ${}y$ ist differenzierbar.
3. Es ist ${}y'(t)=f(t,y(t))$ für alle ${}t\in I$.

Statt Lösung sagt man auch Lösungsfunktion oder Lösungskurve. Es sei betont, dass anders als bei vielen Gleichungen wie quadratische Gleichungen oder lineare Gleichungssysteme, wo die Lösung eine Zahl oder ein Vektor ist, die Lösungen von Differentialgleichungen Funktionen sind.

Differentialgleichungen beschreiben häufig physikalische Prozesse, insbesondere Bewegungsprozesse. Daran soll auch die Notation erinnern, es steht ${}t$ für die Zeit und ${}y$ für den Ort. Dabei ist hier der Ort eindimensional, d.h. die Bewegung findet nur auf einer Geraden statt. Den Wert ${}f(t,y)$ sollte man sich als eine zu einem Zeit- und Ortspunkt vorgegebene Richtung auf der Ortsgeraden vorstellen. Eine Lösung ist dann eine Funktion

$y\colon I\longrightarrow \mathbb {R} ,\,t\longmapsto y(t),$

die differenzierbar ist und deren Ableitung, vorgestellt als Momentangeschwindigkeit, zu jedem Zeitpunkt ${}t$ mit dem durch ${}f(t,y(t))$ gegebenen Richtungsvektor übereinstimmt. Später werden wir auch Bewegungen betrachten, die sich in der Ebene oder im Raum abspielen, und die durch ein entsprechendes Richtungsfeld gesteuert werden.

Die Lösung einer Differentialgleichung ist im Allgemeinen nicht eindeutig, man muss noch Anfangsbedingungen festlegen.

Definition

Es sei ${}U\subseteq \mathbb {R} ^{2}$ eine offene Teilmenge und es sei

$f\colon U\longrightarrow \mathbb {R} ,\,(t,y)\longmapsto f(t,y),$

eine Funktion. Es sei ${}(t_{0},y_{0})\in U$ vorgegeben. Dann nennt man

$y'=f(t,y){\text{ und }}y(t_{0})=y_{0}$

das Anfangswertproblem zur gewöhnlichen Differentialgleichung ${}y'=f(t,y)$ mit der Anfangsbedingung ${}y(t_{0})=y_{0}$.

Definition

Es sei ${}U\subseteq \mathbb {R} ^{2}$ eine offene Teilmenge und es sei

$f\colon U\longrightarrow \mathbb {R} ,\,(t,y)\longmapsto f(t,y),$

eine Funktion. Es sei ${}(t_{0},y_{0})\in U$ vorgegeben. Dann nennt man eine Funktion

$y\colon I\longrightarrow \mathbb {R} ,\,t\longmapsto y(t),$

auf einem Intervall ${}I\subseteq \mathbb {R}$ eine Lösung des Anfangswertproblems

$y'=f(t,y){\text{ und }}y(t_{0})=y_{0},$

wenn ${}y$ eine Lösung der Differentialgleichung ist und wenn zusätzlich

${}y(t_{0})=y_{0}\,$

gilt.

Es gibt kein allgemeines Verfahren eine Differentialgleichung bzw. ein Anfangswertproblem explizit zu lösen. Die Lösbarkeit hängt wesentlich von der gegebenen Funktion ${}f(t,y)$ ab.

Das eine Differentialgleichung beschreibende Vektorfeld ${}f(t,y)$ hängt im Allgemeinen von beiden Variablen ${}t$ und ${}y$ ab. Einfache, aber keineswegs triviale Spezialfälle von Differentialgleichungen liegen vor, wenn das Vektorfeld nur von einer der beiden Variablen abhängt.

Definition

Eine gewöhnliche Differentialgleichung

${}y'=f(t,y)\,$

heißt ortsunabhängig, wenn die Funktion ${}f$ nicht von ${}y$ abhängt, wenn also ${}f(t,y)=g(t)$ mit einer Funktion ${}g$ in der einen Variablen ${}t$ gilt.

Eine ortsunabhängige gewöhnliche Differentialgleichung

$y'=g(t)$

zu einer stetigen Funktion ${}g$ ist nichts anderes als das Problem, eine Stammfunktion ${}G(t)$ von ${}g$ zu finden; eine Lösung ${}y$ der Differentialgleichung ist ja genau durch die Bedingung ausgezeichnet, dass ${}y'(t)=g(t)$ ist. Da eine Stammfunktion nur bis auf die Integrationskonstante bestimmt ist, besitzt ein ortsunabhängiges Anfangswertproblem eine eindeutige Lösung.

Definition

Eine gewöhnliche Differentialgleichung

${}y'=f(t,y)\,$

heißt zeitunabhängig, wenn die Funktion ${}f$ nicht von ${}t$ abhängt, wenn also ${}f(t,y)=h(y)$ mit einer Funktion ${}h$ in der einen Variablen ${}y$ gilt.

Bei einer zeitunabhängigen Differentialgleichung hängt nur das zugrunde liegende Vektorfeld nicht von der Zeit ab, die Lösungskurven sind hingegen im Allgemeinen zeitabhängig.

<< | Kurs:Mathematik (Osnabrück 2009-2011)/Teil II | >> PDF-Version dieser Vorlesung Arbeitsblatt zur Vorlesung (PDF)