Kurs:Analysis (Osnabrück 2021-2023)/Teil I/Vorlesung 28/latex
\setcounter{section}{28}
\zwischenueberschrift{Gewöhnliche Differentialgleichungen}
Wir beginnen mit zwei Beispielen, die beide zu gewöhnlichen Differentialgleichungen führen.
\inputbeispiel{}
{
\bild{ \begin{center}
\includegraphics[width=5.5cm]{\bildeinlesung {COVID-19-Germany.svg} }
\end{center}
\bildtext {} }
\bildlizenz { COVID-19-Germany.svg } {} {Hbf878} {Commons} {CC-by-sa 4.0} {}
\bild{ \begin{center}
\includegraphics[width=5.5cm]{\bildeinlesung {COVID-19-Sweden.svg} }
\end{center}
\bildtext {} }
\bildlizenz { COVID-19-Sweden.svg } {} {Hbf878} {Commons} {CC-by-sa 4.0} {}
Wir versuchen, die Ausbreitung einer Virusinfektion wie bei den Wellen der Corona-Pandemie seit 2020 zu modellieren. Die Ausbreitung wird durch eine Funktion
\maabbeledisp {y} {\R} {\R
} {t} {y(t)
} {,}
beschrieben, wobei
\mavergleichskette
{\vergleichskette
{t
}
{ \in }{\R
}
{ }{
}
{ }{
}
{ }{
}
}
{}{}{}
für die Zeit und
\mathl{y(t)}{} für die Gesamtanzahl der bis zum Zeitpunkt $t$ Infizierten
\zusatzklammer {einschließlich der Genesenen} {} {}
angibt. Dies ist zunächst eine empirische Funktion, die man aus verschiedenen Gründen auch gar nicht genau kennt, insbesondere, da nicht jeder getestet wird. Man kann stattdessen auch die Entwicklung der bestätigt Infizierten betrachten. Diese empirische Funktion wird durch die Daten, die jeden Tag das Robert-Koch-Institut übermittelt, beschrieben, und ist so gesehen zunächst
eine Abbildung von einer Anfangsmenge der natürlichen Zahlen
\zusatzklammer {die ersten $100$ Tage seit Ausbruch} {} {}
in die natürlichen Zahlen, wobei jedem Tag die Anzahl der bis dahin Infizierten zugeordnet wird.
Wenn man zu den Daten aus verschiedenen Ländern
\zusatzklammer {oder verschiedenen Wellen} {} {}
den Verlauf skizziert, ergibt sich jeweils ein ähnliches Bild. Die Ausbreitung scheint einer Gesetzmäßigkeit zu folgen, die man in der \stichwort {mathematischen Modellierung} {} verstehen möchte. Das bedeutet
\zusatzklammer {in einem ersten Schritt} {} {,}
dass man die empirische Funktion, also das vorliegende Datenmaterial, durch eine mathematische Funktion, also einen funktionalen Ausdruck, annähern möchte, um so den qualitativen und den quantitativen Verlauf der Ausbreitung zu verstehen und auch Extrapolationen
\zusatzklammer {Prognosen} {} {}
formulieren zu können. Hierbei wird man den Definitionsbereich und den Wertebereich als die reellen Zahlen
\zusatzklammer {oder Intervalle davon} {} {}
und die Funktion als stetig oder differenzierbar ansetzen. Man kann mit verschiedenen Zeiteinheiten arbeiten und auch die Gesamtzahl absolut oder aber prozentual
\zusatzklammer {bezogen auf die Erdbevölkerung, ein Land, ...} {} {}
angeben. So oder so ergibt sich, dass der Verlauf gut durch eine
\definitionsverweis {Exponentialfunktion}{}{}
beschrieben werden kann, also von der Bauart
\maabbeledisp {} {\R} {\R
} {t} {b^t
} {,}
mit einer Basis
\mavergleichskette
{\vergleichskette
{b
}
{ > }{1
}
{ }{
}
{ }{
}
{ }{
}
}
{}{}{}
ist. Welche Basis $b$ zu nehmen ist, hängt von der Skalierung und auch von länderspezifischen Gegebenheiten ab.
Diese Basis ist äquivalent zum Verdoppelungszeitraum der Ausbreitung, man kann das eine aus dem andern berechnen, siehe
Aufgabe 28.2.
Diese Modellierung ist bisher aber nur die Beobachtung einer Übereinstimmung einer mathematischen Funktionsklasse mit empirischen Funktionen. In einem zweiten Schritt kann man sich fragen, ob es \anfuehrung{in der Natur der Sache liegt}{,} dass die Ausbreitung eines Virus exponentiell verläuft. Gibt es einen mathematischen Grund dafür, eine innere Dynamik, eine zu jedem Zeitpunkt gültige Gesetzmäßigkeit, die den Verlauf erklären kann? Die Antwort zu dieser Frage erfolgt im Rahmen der Theorie der gewöhnlichen Differentialgleichungen, und beruht auf einer einfachen Beobachtung. Wir nehmen die Funktion $y(t)$ als differenzierbar an. Die Ableitung $y'(t)$ beschreibt dann den momentanen Zuwachs zu jedem Zeitpunkt, ist also ein Maß für die Neuansteckungen. Der naheliegende Ansatz ist nun zu sagen, dass zu jedem Zeitpunkt die Anzahl der Infizierten, also $y(t)$, proportional zur Anzahl der Begegnungen zwischen Infizierten und Nichtinfizierten ist und damit proportional zur Anzahl der Neuinfektionen, also zu $y'(t)$
\zusatzklammer {für Einschränkungen zu dieser Überlegung siehe weiter unten} {} {.}
Dies führt zur Beziehung
\mavergleichskettedisp
{\vergleichskette
{ y'(t)
}
{ =} { c y(t)
}
{ } {
}
{ } {
}
{ } {
}
}
{}{}{}
mit einem konstanten Proportionalitätsfaktor $c$, der ein Maß für die Ansteckungswahrscheinlichkeit ist und vom Virus, der Saison, vom Abstandsverhalten der Bevölkerung u. Ä. abhängt. Wir haben also eine Beziehung zwischen der gesuchten Funktion und ihrer Ableitung, die in jedem Moment gilt und für die Ausbreitung eines Virus charakteristisch sein sollte. Ein solcher Zusammenhang zwischen einer Funktion und ihrer Ableitung heißt eine \stichwort {gewöhnliche Differentialgleichung} {.} Wenn eine solche Differentialgleichung vorliegt, fragt man sich, welche Funktionen $y(t)$ diese Gleichung erfüllen. Dies ist im Allgemeinen schwierig. Im vorliegenden Fall lässt sich direkt durch Ableiten bestätigen, dass die Funktionen
\mavergleichskettedisp
{\vergleichskette
{ y(t)
}
{ =} { a e^{ct}
}
{ } {
}
{ } {
}
{ } {
}
}
{}{}{}
mit
\mavergleichskette
{\vergleichskette
{ a
}
{ \in }{\R
}
{ }{
}
{ }{
}
{ }{
}
}
{}{}{}
Lösungen sind. Der Vorfaktor $a$ ist dabei durch
\mavergleichskettedisp
{\vergleichskette
{y(0)
}
{ =} {a
}
{ } {
}
{ } {
}
{ } {
}
}
{}{}{}
festgelegt, also durch den Wert der Funktion zum Zeitpunkt $0$, und das $c$ im Exponenten ist direkt der Proportionalitätsfaktor aus der Differentialgleichung. Wegen
\mavergleichskettedisp
{\vergleichskette
{ a e^{ct}
}
{ =} { e^{ \ln a } e^{ct}
}
{ =} { e^{ \ln a + ct}
}
{ } {
}
{ } {
}
}
{}{}{}
ist der Vorfaktor $a$ im Wesentlichen eine Verschiebung im Zeitargument, und $c$ kann man durch eine Umskalierung der Zeit zu $1$ normieren. Man kann nun sogar zeigen, dass die Exponentialfunktionen die einzigen Funktionen sind, die diese Differentialgleichung erfüllen, siehe
Aufgabe 20.37
bzw.
Aufgabe 28.5.
Dies bedeutet, dass eine Virusausbreitung durch den Faktor $c$ und dem Wert an einem einzigen Zeitpunkt eindeutig bestimmt ist. Dies ist ein Spezialfall des Satzes, dass ein Anfangswertproblem eine eindeutige Lösung besitzt, von dem wir verschiedene Varianten kennenlernen werden.
Kommen wir nun zu einigen Einschränkungen der oben formulierten Modellierung. Zunächst ist klar, dass die Exponentialfunktion zu jeder Basis
\mavergleichskette
{\vergleichskette
{b
}
{ > }{1
}
{ }{
}
{ }{
}
{ }{
}
}
{}{}{}
gegen unendlich geht, es aber nur endlich viele Menschen gibt. Also kann irgendwas nicht stimmen. Der Punkt ist, dass in unserer Modellierung die Anzahl der Infizierten zur Anzahl der Begegnungen von Infizierten mit der Gesamtbevölkerung proportional ist, aber nicht mit der Anzahl der Begegnungen mit den Nichtinfizierten. Dieser Unterschied ist zu Beginn der Ausbreitung unerheblich, da zu Beginn die Gesamtbevölkerung nahezu vollständig nicht infiziert ist. Im Verlauf der Epidemie, wenn sich der Durchseuchungsgrad erhöht, wird es zunehmend wahrscheinlicher, dass sich Infizierte und Infizierte begegnen, was zu keiner Neuansteckung führt.
Ferner haben wir ignoriert, dass die Genesenen nicht mehr andere Leute anstecken können. Hier muss man den Unterschied zwischen infiziert und akut infiziert berücksichtigen. Dieser Unterschied ist für den Anfangsverlauf der Ausbreitung ebenfalls unerheblich, spielt aber im späteren Verlauf eine wichtige Rolle. Die Neuansteckung ist also proportional zur Anzahl der akut Infizierten, dies ist die Differenz zwischen der Gesamtinfiziertenzahl und der Gesamtinfiziertenzahl vor einem gewissen Genesungszeitraum $d$
\zusatzklammer {bei Corona ca. $2-3$ Wochen} {} {.}
Dies führt auf die Bedingung
\mavergleichskettedisp
{\vergleichskette
{y'(t)
}
{ =} { c (y(t) - y(t-d))
}
{ } {
}
{ } {
}
{ } {
}
}
{}{}{,}
man spricht von einer \stichwort {Differentialgleichung mit Verzögerung} {,} was wir nicht behandeln werden. Für den ersten Zeitraum der Länge $d$ nach Ausbruch spielt der Korrekturterm aber keine Rolle.
Schließlich ist der Faktor $c$ keine Konstante, sondern wird durch politische Maßnahmen und Verhaltensregeln beeinflusst.
}
\inputbeispiel{}
{
\bild{ \begin{center}
\includegraphics[width=5.5cm]{\bildeinlesung {Taraxacum_sect_Ruderalia13_ies.jpg} }
\end{center}
\bildtext {} }
\bildlizenz { Taraxacum sect Ruderalia13 ies.jpg } {Frank Vincentz} {} {Commons} {CC-by-sa 3.0} {}
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 \anfuehrung{Rückkopplungsprozess}{} vorliegt: Die bisherige Bewegung \zusatzklammer {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 \zusatzklammer {häufig \mathlk{y(t)}{}} {} {} und der Änderung dieser Größe \zusatzklammer {\mathlk{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 $\R^2$ definiert sein, sondern auf einer \stichwort {offenen Menge} {}
\mavergleichskette
{\vergleichskette
{ U
}
{ \subseteq }{ \R^2
}
{ }{
}
{ }{
}
{ }{
}
}
{}{}{.}
Dabei heißt $U$ offen, wenn es zu jedem Punkt
\mavergleichskette
{\vergleichskette
{ P
}
{ \in }{ U
}
{ }{
}
{ }{
}
{ }{
}
}
{}{}{}
eine offene Ballumgebung
\mathl{U(P,r)}{} gibt, die ganz in $U$ liegt.
\inputdefinition
{}
{
Es sei
\mavergleichskette
{\vergleichskette
{ U
}
{ \subseteq }{ \R^2
}
{ }{
}
{ }{
}
{ }{
}
}
{}{}{}
eine offene Teilmenge und es sei
\maabbeledisp {f} {U} {\R
} {(t,y)} {f(t,y)
} {,}
eine
\definitionsverweis {Funktion}{}{.} Dann nennt man
\mavergleichskettedisp
{\vergleichskette
{y'
}
{ =} {f(t,y)
}
{ } {
}
{ } {
}
{ } {
}
}
{}{}{}
die
\zusatzklammer {gewöhnliche} {} {}
\definitionswort {Differentialgleichung}{} zu $f$
\zusatzklammer {oder zum \stichwort {Vektorfeld} {} oder zum \stichwort {Richtungsfeld} {} $f$} {} {.}
}
Dabei ist
\mavergleichskette
{\vergleichskette
{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 \stichwort {Lösung einer Differentialgleichung} {.}
\inputdefinition
{}
{
Es sei
\mavergleichskette
{\vergleichskette
{ U
}
{ \subseteq }{ \R^2
}
{ }{
}
{ }{
}
{ }{
}
}
{}{}{}
eine offene Teilmenge und es sei
\maabbeledisp {f} {U} {\R
} {(t,y)} {f(t,y)
} {,}
eine
\definitionsverweis {Funktion}{}{.} Zur
\definitionsverweis {gewöhnlichen Differentialgleichung}{}{}
\mavergleichskettedisp
{\vergleichskette
{y'
}
{ =} {f(t,y)
}
{ } {
}
{ } {
}
{ } {
}
}
{}{}{}
heißt eine
\definitionsverweis {Funktion}{}{}
\maabbeledisp {y} {I} {\R
} {t} {y(t)
} {,}
auf einem
\zusatzklammer {mehrpunktigen} {} {}
\definitionsverweis {Intervall}{}{}
\mavergleichskette
{\vergleichskette
{I
}
{ \subseteq }{\R
}
{ }{
}
{ }{
}
{ }{
}
}
{}{}{}
eine \definitionswort {Lösung der Differentialgleichung}{,} wenn folgende Eigenschaften erfüllt sind.
\aufzaehlungdrei{Es ist
\mavergleichskette
{\vergleichskette
{ (t,y(t))
}
{ \in }{U
}
{ }{
}
{ }{
}
{ }{
}
}
{}{}{}
für alle
\mavergleichskette
{\vergleichskette
{t
}
{ \in }{I
}
{ }{
}
{ }{
}
{ }{
}
}
{}{}{.}
}{Die Funktion $y$ ist
\definitionsverweis {differenzierbar}{}{.}
}{Es ist
\mavergleichskette
{\vergleichskette
{y'(t)
}
{ = }{f(t,y(t))
}
{ }{
}
{ }{
}
{ }{
}
}
{}{}{}
für alle
\mavergleichskette
{\vergleichskette
{t
}
{ \in }{I
}
{ }{
}
{ }{
}
{ }{
}
}
{}{}{.}
}
} Statt Lösung sagt man auch \stichwort {Lösungsfunktion} {} oder \stichwort {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
\mathl{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
\maabbeledisp {y} {I} { \R
} {t} { y(t)
} {,}
die differenzierbar ist und deren Ableitung, vorgestellt als Momentangeschwindigkeit, zu jedem Zeitpunkt $t$ mit dem durch
\mathl{f(t,y(t))}{} gegebenen Richtungsvektor übereinstimmt. In Analysis II 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.
\inputdefinition
{}
{
Es sei
\mavergleichskette
{\vergleichskette
{ U
}
{ \subseteq }{ \R^2
}
{ }{
}
{ }{
}
{ }{
}
}
{}{}{}
eine offene Teilmenge und es sei
\maabbeledisp {f} {U} {\R
} {(t,y)} {f(t,y)
} {,}
eine
\definitionsverweis {Funktion}{}{.} Es sei
\mavergleichskette
{\vergleichskette
{ (t_0,y_0)
}
{ \in }{U
}
{ }{
}
{ }{
}
{ }{
}
}
{}{}{}
vorgegeben. Dann nennt man
\mathdisp {y'=f(t,y) \text{ und } y(t_0)=y_0} { }
das \definitionswort {Anfangswertproblem}{}
zur
\definitionsverweis {gewöhnlichen Differentialgleichung}{}{}
\mavergleichskette
{\vergleichskette
{y'
}
{ = }{f(t,y)
}
{ }{
}
{ }{
}
{ }{
}
}
{}{}{}
mit der \definitionswort {Anfangsbedingung}{}
\mavergleichskette
{\vergleichskette
{ y(t_0)
}
{ = }{y_0
}
{ }{
}
{ }{
}
{ }{
}
}
{}{}{.}
}
\inputdefinition
{}
{
Es sei
\mavergleichskette
{\vergleichskette
{ U
}
{ \subseteq }{ \R^2
}
{ }{
}
{ }{
}
{ }{
}
}
{}{}{}
eine offene Teilmenge und es sei
\maabbeledisp {f} {U} {\R
} {(t,y)} {f(t,y)
} {,}
eine
\definitionsverweis {Funktion}{}{.} Es sei
\mavergleichskette
{\vergleichskette
{ (t_0,y_0)
}
{ \in }{ U
}
{ }{
}
{ }{
}
{ }{
}
}
{}{}{}
vorgegeben. Dann nennt man eine
\definitionsverweis {Funktion}{}{}
\maabbeledisp {y} {I} {\R
} {t} {y(t)
} {,}
auf einem
\definitionsverweis {Intervall}{}{}
\mavergleichskette
{\vergleichskette
{I
}
{ \subseteq }{ \R
}
{ }{
}
{ }{
}
{ }{
}
}
{}{}{}
eine \definitionswort {Lösung des Anfangswertproblems}{}
\mathdisp {y'=f(t,y) \text{ und } y(t_0)=y_0} { , }
wenn $y$ eine
\definitionsverweis {Lösung der Differentialgleichung}{}{}
ist und wenn zusätzlich
\mavergleichskettedisp
{\vergleichskette
{ 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
\mathl{f(t,y)}{} ab.
Das eine Differentialgleichung beschreibende Vektorfeld
\mathl{f(t,y)}{} hängt im Allgemeinen von beiden Variablen
\mathkor {} {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.
\zwischenueberschrift{Ortsunabhängige Differentialgleichungen}
\inputdefinition
{}
{
Eine
\definitionsverweis {gewöhnliche Differentialgleichung}{}{}
\mavergleichskettedisp
{\vergleichskette
{ y'
}
{ =} { f(t,y)
}
{ } {
}
{ } {
}
{ } {
}
}
{}{}{}
heißt \definitionswort {ortsunabhängig}{,} wenn die Funktion $f$ nicht von $y$ abhängt, wenn also
\mavergleichskette
{\vergleichskette
{ f(t,y)
}
{ = }{ g(t)
}
{ }{
}
{ }{
}
{ }{
}
}
{}{}{}
mit einer Funktion $g$ in der einen Variablen $t$ gilt.
}
Eine ortsunabhängige gewöhnliche Differentialgleichung
\mavergleichskettedisp
{\vergleichskette
{ y'
}
{ =} { g(t)
}
{ } {
}
{ } {
}
{ } {
}
}
{}{}{}
zu einer stetigen Funktion $g$ ist nichts anderes als das Problem, eine Stammfunktion
\mathl{G(t)}{}
von $g$ zu finden; eine Lösung $y$ der Differentialgleichung ist ja genau durch die Bedingung ausgezeichnet, dass
\mavergleichskette
{\vergleichskette
{ y'(t)
}
{ = }{ g(t)
}
{ }{
}
{ }{
}
{ }{
}
}
{}{}{}
ist. Da eine Stammfunktion nur bis auf die Integrationskonstante bestimmt ist, besitzt ein ortsunabhängiges Anfangswertproblem eine eindeutige Lösung.
\inputbeispiel{}
{
Wir betrachten das
\definitionsverweis {ortsunabhängige}{}{}
\definitionsverweis {Anfangswertproblem}{}{}
\mathdisp {y'= { \frac{ 1 }{ t^2-1 } } \text{ mit der Anfangsbedingung } y(5)= 3} { . }
Die Funktion
\mathl{{ \frac{ 1 }{ t^2-1 } }}{} besitzt die
\definitionsverweis {Partialbruchzerlegung}{}{}
\mavergleichskettedisp
{\vergleichskette
{ { \frac{ 1 }{ t^2-1 } }
}
{ =} { { \frac{ 1 }{ 2 } } \cdot { \frac{ 1 }{ t-1 } } - { \frac{ 1 }{ 2 } } \cdot { \frac{ 1 }{ t + 1 } }
}
{ } {
}
{ } {
}
{ } {
}
}
{}{}{,}
daher sind die Stammfunktionen
\zusatzklammer {wir beschränken uns auf
\mavergleichskettek
{\vergleichskettek
{t
}
{ > }{1
}
{ }{
}
{ }{
}
{ }{
}
}
{}{}{}} {} {}
gleich
\mavergleichskettedisp
{\vergleichskette
{y(t)
}
{ =} { { \frac{ 1 }{ 2 } } \cdot \ln (t-1) - { \frac{ 1 }{ 2 } } \cdot \ln ( t +1) +c
}
{ } {
}
{ } {
}
{ } {
}
}
{}{}{.}
Die Anfangsbedingung
\mavergleichskette
{\vergleichskette
{ y(5)
}
{ = }{ 3
}
{ }{
}
{ }{
}
{ }{
}
}
{}{}{}
führt auf
\mavergleichskettedisp
{\vergleichskette
{ { \frac{ 1 }{ 2 } } \cdot \ln 4 - { \frac{ 1 }{ 2 } } \cdot \ln 6 +c
}
{ =} { 3
}
{ } {
}
{ } {
}
{ } {
}
}
{}{}{,}
also ist
\mavergleichskettedisp
{\vergleichskette
{ c
}
{ =} { 3 - { \frac{ 1 }{ 2 } } \cdot \ln 4 + { \frac{ 1 }{ 2 } } \cdot \ln 6
}
{ } {
}
{ } {
}
{ } {
}
}
{}{}{}
und die Lösungsfunktion des Anfangswertproblems ist
\mavergleichskettedisphandlinks
{\vergleichskettedisphandlinks
{ y(t)
}
{ =} { { \frac{ 1 }{ 2 } } \cdot \ln (t-1) - { \frac{ 1 }{ 2 } } \cdot \ln (t+1) + 3 - { \frac{ 1 }{ 2 } } \cdot \ln 4 + { \frac{ 1 }{ 2 } } \cdot \ln 6
}
{ } {
}
{ } {
}
{ } {
}
}
{}{}{.}
}
\inputbeispiel{}
{
Wir betrachten das
\definitionsverweis {ortsunabhängige}{}{}
\definitionsverweis {Anfangswertproblem}{}{}
\mathdisp {y' = { \frac{ 1 }{ \cosh t } } \text{ mit der Anfangsbedingung } y(0)= 5} { . }
Es ist
\mavergleichskettedisp
{\vergleichskette
{ { \frac{ 1 }{ \cosh t } }
}
{ =} { { \frac{ 2 }{ e^t +e^{-t} } }
}
{ =} { { \frac{ 2 e^t }{ e^{2t} + 1 } }
}
{ } {
}
{ } {
}
}
{}{}{,}
sodass eine rationale Funktion in der Exponentialfunktion vorliegt, die wir nach
Lemma 27.1
über die Substitution
\mathl{t= \ln s}{} lösen können. Das transformierte Integral ist dabei
\mathdisp {\int { \frac{ 2s }{ s^2 +1 } } \cdot { \frac{ 1 }{ s } } ds} { . }
Eine Stammfunktion dazu ist
\mathdisp {2 \arctan s} { . }
Somit ist
\mathdisp {2 \arctan { \left( e^t \right) }} { }
eine Stammfunktion von
\mathl{{ \frac{ 1 }{ \cosh t } }}{.} Für das Anfangswertproblem setzen wir
\mavergleichskettedisp
{\vergleichskette
{ 2 \arctan { \left( e^0 \right) } +c
}
{ =} { 5
}
{ } {
}
{ } {
}
{ } {
}
}
{}{}{}
an. Dies führt auf
\mavergleichskettedisp
{\vergleichskette
{c
}
{ =} { 5- 2 \arctan 1
}
{ } {
}
{ } {
}
{ } {}
}
{}{}{,}
also ist
\mavergleichskettedisp
{\vergleichskette
{y(t)
}
{ =} {2 \arctan { \left( e^t \right) } +5 -2 \arctan 1
}
{ } {
}
{ } {
}
{ } {
}
}
{}{}{}
die Lösung des Anfangswertproblems.
}
\zwischenueberschrift{Zeitunabhängige Differentialgleichungen}
\inputdefinition
{}
{
Eine
\definitionsverweis {gewöhnliche Differentialgleichung}{}{}
\mavergleichskettedisp
{\vergleichskette
{y'
}
{ =} {f(t,y)
}
{ } {
}
{ } {
}
{ } {
}
}
{}{}{}
heißt \definitionswort {zeitunabhängig}{,} wenn die Funktion $f$ nicht von $t$ abhängt, wenn also
\mavergleichskette
{\vergleichskette
{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 \anfuehrung{Vektorfeld}{} nicht von der Zeit ab, die Lösungsfunktionen sind hingegen im Allgemeinen zeitabhängig.
\inputbeispiel{}
{
Wir betrachten die zeitliche Entwicklung einer Population, die durch folgende Eigenschaften charakterisiert ist.
\aufzaehlungvier{Die Individuen der Population leben ewig.
}{Alle Individuen beteiligen sich ab ihrer Geburt mit gleichem
\zusatzklammer {durchschnittlichen} {} {}
Engagement und Erfolg an der Fortpflanzung.
}{Zeugung und Geburt finden gleichzeitig statt.
}{Der Fortpflanzungserfolg eines Individuums ist unabhängig von der Größe der Gesamtpopulation.
}
Unter diesen Bedingungen ist die Vermehrung, also der Zuwachs der Population, allein von der momentanen Populationsgröße abhängig und proportional zu dieser. Wenn man die Populationsentwicklung als
\mathl{y(t)}{} ansetzt, so erhält man eine gewöhnliche Differentialgleichung
\mavergleichskettedisp
{\vergleichskette
{y'(t)
}
{ =} { c y(t)
}
{ } {
}
{ } {
}
{ } {
}
}
{}{}{}
\zusatzklammer {oder kurz
\mavergleichskettek
{\vergleichskettek
{ y'
}
{ = }{ cy
}
{ }{
}
{ }{
}
{ }{
}
}
{}{}{}} {} {}
mit einer Konstanten
\mavergleichskette
{\vergleichskette
{c
}
{ \in }{\R_+
}
{ }{
}
{ }{
}
{ }{
}
}
{}{}{.}
Die Lösungsfunktionen sind
\mathdisp {\lambda e^{ct}} { }
\zusatzklammer {wobei im Populationsbeispiel
\mavergleichskettek
{\vergleichskettek
{ \lambda
}
{ \in }{ \R_+
}
{ }{
}
{ }{
}
{ }{
}
}
{}{}{}
ist} {} {.}
Man spricht von \stichwort {exponentiellem Wachstum} {} der Population, und zwar unabhängig davon, ob $c$ groß oder klein ist.
}
\inputbeispiel{}
{
Eine Wüste
\zusatzklammer {oder ein Kornblumenfeld} {} {}
sei kreisrund und breite sich mit der Zeit kontinuierlich aus, indem die Grenze gleichmäßig nach außen geschoben werde, und zwar pro Zeiteinheit um einen gewissen Vortrieb. Die Fläche der Wüste werde durch die Funktion
\mathl{z(t)}{} beschrieben. Die Grenze der Wüste hat somit die Länge
\mathl{2 \sqrt{\pi} \sqrt{z(t)}}{} und diese Länge ist proportional zum Wüstenzuwachs zum Zeitpunkt $t$. Es ergibt sich daher eine Differentialgleichung
\mavergleichskettedisp
{\vergleichskette
{z'(t)
}
{ =} { c \sqrt{z(t)}
}
{ } {
}
{ } {
}
{ } {
}
}
{}{}{}
mit einer Konstanten
\mavergleichskette
{\vergleichskette
{c
}
{ \in }{ \R_+
}
{ }{
}
{ }{
}
{ }{
}
}
{}{}{.}
Die Lösungen haben die Form
\mavergleichskettedisp
{\vergleichskette
{z(t)
}
{ =} { { \frac{ c^2 }{ 4 } } t^2
}
{ } {
}
{ } {
}
{ } {
}
}
{}{}{,}
wie man direkt durch Ableiten bestätigen kann.
}
\zwischenueberschrift{Differentialgleichungen höherer Ordnung}
Viele physikalische Bewegungsprozesse sind nicht dadurch determiniert, dass zu jedem Zeit- und Ortspunkt die Bewegungsrichtung \zusatzklammer {also die gerichtete Geschwindigkeit} {} {} vorgegeben wird, sondern dadurch, dass zu jedem Zeit- und Ortspunkt eine Kraft auf ein Teilchen wirkt, die dieses beschleunigt. In diesem Fall kann die Bewegung also nicht durch die erste Ableitung \zusatzklammer {Geschwindigkeit} {} {} modelliert werden, sondern durch die zweite Ableitung \zusatzklammer {Beschleunigung} {} {.} Typische Beispiele hierzu sind die durch die Gravitation oder eine Federkraft hervorgerufenen Bewegungen.
\bild{ \begin{center}
\includegraphics[width=5.5cm]{\bildeinlesung {Galileo.arp.300pix.jpg} }
\end{center}
\bildtext {Galileo Galilei (1564-1642) entdeckte das Gesetz des freien Falls.} }
\bildlizenz { Galileo.arp.300pix.jpg } {} {} {Commons} {} {}
\inputbeispiel{}
{
Ein Gegenstand der Masse $m$ wird im Vakuum aus einer Höhe $0$ zum Zeitpunkt $0$ losgelassen und fällt unter dem Einfluss der Gravitation zu Boden
\zusatzklammer {freier Fall im Vakuum} {} {.}
Dabei wirkt auf den Körper die Gravitationskraft
\mathl{gm}{}
\zusatzklammer {die Erdbeschleunigung $g$ nehmen wir für diesen Bewegungsvorgang als konstant an} {} {,}
die ihn nach dem Gesetz \anfuehrung{Kraft ist Masse mal Beschleunigung}{} beschleunigt. Die Beschleunigung ist also konstant und unabhängig von der Masse. Dies bedeutet, dass die Geschwindigkeit
\mathl{v(t)}{} des Körpers die Differentialgleichung
\mavergleichskettedisp
{\vergleichskette
{v'(t)
}
{ =} {-g
}
{ } {
}
{ } {
}
{ } {
}
}
{}{}{}
erfüllt
\zusatzklammer {die Wahl des Vorzeichens bewirkt, dass der Körper ins Negative fällt} {} {.}
Die durch die Anfangsbedingung
\zusatzklammer {der Gegenstand ruhe zum Zeitpunkt $0$} {} {}
\mavergleichskette
{\vergleichskette
{v(0)
}
{ = }{ 0
}
{ }{
}
{ }{
}
{ }{
}
}
{}{}{}
festgelegte Lösung für die Geschwindigkeit ist daher
\mavergleichskettedisp
{\vergleichskette
{v(t)
}
{ =} {-gt
}
{ } {
}
{ } {
}
{ } {
}
}
{}{}{.}
Der zurückgelegte Weg
\mathl{y(t)}{} des Körpers ergibt sich wiederum aus der Differentialgleichung
\mavergleichskettedisp
{\vergleichskette
{y'(t)
}
{ =} {v(t)
}
{ =} {-gt
}
{ } {
}
{ } {
}
}
{}{}{,}
die besagt, dass die Ableitung des Weges nach der Zeit die Momentangeschwindigkeit beschreibt. Die Lösung davon ist
\mavergleichskettedisp
{\vergleichskette
{ y(t)
}
{ =} { - { \frac{ g }{ 2 } } t^2
}
{ } {
}
{ } {
}
{ } {
}
}
{}{}{.}
Den Gesamtvorgang kann man durch die Differentialgleichung zweiter Ordnung
\mavergleichskettedisp
{\vergleichskette
{ y^{\prime \prime}
}
{ =} { - g
}
{ } {
}
{ } {
}
{ } {
}
}
{}{}{}
ausdrücken.
}
\inputbeispiel{}
{
Ein Gegenstand der Masse $m$ wird aus der Höhe losgelassen und fällt unter dem Einfluss der Gravitation zu Boden. Dabei wirkt auf den Körper einerseits die Gravitationskraft
\mathl{gm}{}
\zusatzklammer {die Erdbeschleunigung $g$ nehmen wir für diesen Bewegungsvorgang als konstant an} {} {,}
die ihn beschleunigt, andererseits wird diese Beschleunigung durch den Luftwiderstand verringert. Nach einem physikalischen Gesetz ist die Reibung
\zusatzklammer {bei relativ kleinen Geschwindigkeiten} {} {}
proportional und entgegengesetzt zur Geschwindigkeit des Körpers. Es sei $\beta$ der Reibungswiderstand, also dieser Proportionalitätsfaktor. Die auf den Körper
\zusatzklammer {nach unten} {} {}
wirkende Gesamtkraft ist daher
\mavergleichskettedisp
{\vergleichskette
{F(t)
}
{ =} { gm - \beta y'(t)
}
{ } {
}
{ } {
}
{ } {
}
}
{}{}{.}
Wegen
\mavergleichskettedisp
{\vergleichskette
{ y^{\prime \prime} (t)
}
{ =} { { \frac{ F(t) }{ m } }
}
{ } {
}
{ } {
}
{ } {
}
}
{}{}{}
gilt daher für diesen Bewegungsvorgang die Differentialgleichung zweiter Ordnung
\mavergleichskettedisp
{\vergleichskette
{ y^{\prime \prime}
}
{ =} {- { \frac{ \beta }{ m } } y' +g
}
{ } {
}
{ } {
}
{ } {
}
}
{}{}{.}
Wenn wir dies mit der Ableitungsfunktion
\mavergleichskette
{\vergleichskette
{v
}
{ = }{y'
}
{ }{
}
{ }{
}
{ }{
}
}
{}{}{}
schreiben, so erhalten wir die Bedingung
\mavergleichskettedisp
{\vergleichskette
{v'
}
{ =} { - { \frac{ \beta }{ m } } v +g
}
{ } {
}
{ } {
}
{ } {
}
}
{}{}{,}
die nach
Beispiel 29.11
die Lösungen
\mavergleichskettedisp
{\vergleichskette
{v(t)
}
{ =} { c e^{- { \frac{ \beta }{ m } } t } + { \frac{ gm }{ \beta } }
}
{ } {
}
{ } {
}
{ } {
}
}
{}{}{}
besitzt. Durch Intergration erhält man für die Differentialgleichung zweiter Ordnung die Lösungsfunktionen
\mavergleichskettedisp
{\vergleichskette
{ y(t)
}
{ =} { - c { \frac{ m }{ \beta } } e^{- { \frac{ \beta }{ m } } t } + { \frac{ g m }{ \beta } } t +d
}
{ } {
}
{ } {
}
{ } {
}
}
{}{}{}
mit beliebigen Konstanten
\mavergleichskette
{\vergleichskette
{c,d
}
{ \in }{ \R
}
{ }{
}
{ }{
}
{ }{
}
}
{}{}{.}
Siehe auch
Aufgabe 29.21.
}
\bild{ \begin{center}
\includegraphics[width=5.5cm]{\bildeinlesung {Simple_Harmonic_Motion_Orbit.gif} }
\end{center}
\bildtext {} }
\bildlizenz { Simple Harmonic Motion Orbit.gif } {} {Mazemaster} {Commons} {PD} {}
\inputbeispiel{}
{
Wir betrachten die Bewegung eines Punktes auf einer Geraden, wobei die auf den Punkt
\zusatzklammer {in Richtung des Nullpunkts} {} {}
wirkende Kraft
\zusatzklammer {bzw. Beschleunigung} {} {}
proportional zur Lage des Punktes sein soll. Wenn der Punkt sich in $\R_+$ befindet und sich in die positive Richtung bewegt, so wirkt diese Kraft bremsend, wenn er sich in die negative Richtung bewegt, so wirkt die Kraft beschleunigend. Mit der Proportionalitätskonstante $1$ gelangt man zur
\definitionsverweis { Differentialgleichung}{}{}
\zusatzklammer {zweiter Ordnung} {} {}
\mavergleichskettedisp
{\vergleichskette
{ y^{\prime \prime}
}
{ =} { -y
}
{ } {
}
{ } {
}
{ } {
}
}
{}{}{,}
die diesen Bewegungsvorgang beschreibt. Als Anfangsbedingung wählen wir
\mathkor {} {y(0)=0} {und} {y'(0)=v} {,}
zum Zeitpunkt $0$ soll die Bewegung also durch den Nullpunkt gehen und dort die Geschwindigkeit $v$ besitzen. Man kann sofort die Lösung
\mavergleichskettedisp
{\vergleichskette
{ y(t)
}
{ =} { v \cdot \sin t
}
{ } {
}
{ } {
}
{ } {
}
}
{}{}{}
angeben.
}