Kurs:Analysis (Osnabrück 2013-2015)/Teil II/Vorlesung 40

Wir haben schon im ersten Semester gewöhnliche Differentialgleichungen samt einiger Lösungsverfahren besprochen. Dort ging es um die Bewegungen auf einer Geraden, die durch ein von der Zeit und dem Ort (der Lage auf der Geraden) abhängigen Vektorfeld bestimmt wurden. Eine physikalische Bewegung spielt sich aber häufig höherdimensional (im $ℝ^{2}$ oder im $ℝ^{3}$ ) ab, so dass wir jetzt gewöhnliche Differentialgleichungen allgemein besprechen. Die Zeitkomponente wird sich nach wie vor in einem reellen Intervall bewegen, die Orts- und die Richtungskomponente wird ein Element in einem beliebigen endlichdimensionalen reellen Vektorraum sein. Diesen statten wir mit einem Skalarprodukt aus, so dass wir eine Norm, eine Metrik, offene Mengen, stetige Abbildungen, Limiten, etc. zur Verfügung haben.

Gewöhnliche Differentialgleichungen

Definition

Es sei $V$ ein endlichdimensionaler reeller Vektorraum, $I \subseteq ℝ$ ein reelles Intervall, $U \subseteq V$ eine offene Menge und

f : I \times U ⟶ V, (t, v) ⟼ f (t, v),

ein Vektorfeld auf $U$ . Dann nennt man

v^{'} = f (t, v)

die gewöhnliche Differentialgleichung (oder gewöhnliches Differentialgleichungssystem) zum Vektorfeld $f$ .

(Zeitabhängige) Vektorfelder und gewöhnliche Differentialgleichungssysteme sind im Wesentlichen äquivalente Objekte. Man spricht auch von einem dynamischen System. Von Differentialgleichungen spricht man insbesondere dann, wenn man sich für die Lösungen im Sinne der folgenden Definition interessiert.

Definition

Es sei $V$ ein endlichdimensionaler reeller Vektorraum, $I \subseteq ℝ$ ein reelles Intervall, $U \subseteq V$ eine offene Menge und

f : I \times U ⟶ V, (t, v) ⟼ f (t, v),

ein Vektorfeld auf $U$ . Zur gewöhnlichen Differentialgleichung

v^{'} = f (t, v)

heißt eine Abbildung

v : J ⟶ V, t ⟼ v (t),

auf einem offenen (Teil)Intervall^[1] $J \subseteq I$ eine Lösung der Differentialgleichung, wenn folgende Eigenschaften erfüllt sind.

Es ist $v (t) \in U$ für alle $t \in J$ .
Die Abbildung $v$ ist differenzierbar.
Es ist $v^{'} (t) = f (t, v (t))$ für alle $t \in J$ .

Eine Lösung ist also eine differenzierbare Kurve, d.h. eine (orts-)vektorwertige Abbildung

v : J ⟶ V, t ⟼ v (t) .

Wenn $V = ℝ^{n}$ ist, so wird eine solche Abbildung durch ihre Komponenten

(v_{1} (t), \dots, v_{n} (t))

beschrieben. Ebenso wird das Vektorfeld durch $n$ , von $t$ und $v = (v_{1}, \dots, v_{n})$ abhängige Funktionen $(f_{1}, \dots, f_{n})$ beschrieben. Die Differentialgleichung lautet dann ausgeschrieben

(\begin{matrix} {v_{1}}^{'} \\ ⋮ \\ {v_{n}}^{'} \end{matrix}) = (\begin{matrix} f_{1} (t, v_{1}, \dots, v_{n}) \\ ⋮ \\ f_{n} (t, v_{1}, \dots, v_{n}) \end{matrix}) .

Daher spricht man auch von einem Differentialgleichungssystem.

Häufig soll eine Kurve nicht nur eine Differentialgleichung erfüllen, sondern sich zusätzlich zu einem bestimmten Zeitpunkt an einem bestimmten Ort befinden. Dies führt zum Begriff des Anfangswertproblems.

Definition

Es sei $V$ ein endlichdimensionaler reeller Vektorraum, $I \subseteq ℝ$ ein reelles Intervall, $U \subseteq V$ eine offene Menge und

f : I \times U ⟶ V, (t, v) ⟼ f (t, v),

ein Vektorfeld auf $U$ . Es sei $(t_{0}, w) \in I \times U$ gegeben. Dann nennt man

v^{'} = f (t, v) und v (t_{0}) = w

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

Definition

Es sei $V$ ein endlichdimensionaler reeller Vektorraum, $I \subseteq ℝ$ ein reelles Intervall, $U \subseteq V$ eine offene Menge und

f : I \times U ⟶ V, (t, v) ⟼ f (t, v),

ein Vektorfeld auf $U$ . Es sei $(t_{0}, w) \in I \times U$ vorgegeben. Dann nennt man eine Abbildung

v : J ⟶ V, t ⟼ v (t),

auf einem Intervall $J \subseteq I$ mit $t_{0} \in J$ eine Lösung des Anfangswertproblems

v^{'} = f (t, v) und v (t_{0}) = w,

wenn $v$ eine Lösung der Differentialgleichung $v^{'} = f (t, v)$ ist und wenn zusätzlich

v (t_{0}) = w

gilt.

Eine zu einem Vektorfeld, einer gewöhnlichen Differentialgleichung und einem Anfangswertproblem passende Vorstellung ist das Windmodell. Das Vektorfeld

F : I \times U ⟶ V

beschreibt zu einem jeden Zeitpunkt $t \in I$ und einem Ortspunkt $P \in U$ die in diesem Punkt herrschende Windrichtung (oder Windgeschwindigkeit). Die Lösung einer Differentialgleichung ist die Bewegung eines Teilchens, das (beschleunigungsfrei und verzögerungsfrei) vom Wind getragen wird, dessen Momentangeschwindigkeit also zu jedem Zeitpunkt gleich der Windgeschwindigkeit an dem Ort ist, an dem sich das Teilchen gerade befindet. Die Lösung eines Anfangswertproblems beschreibt die Bewegung, wenn das Teilchen an einem bestimmten Punkt losgelassen wird.

Die Vorstellung, dass eine Differentialgleichung die Bewegung in einem Kraftfeld^[2] beschreibt, kann irreführend sein. Ein Kraftfeld ist ein Beschleunigungsfeld und kein Geschwindigkeitsfeld. Allerdings führt ein Kraftfeld zu einer Differentialgleichung zweiter Ordnung, die in eine Differentialgleichung erster Ordnung (unter Hinzunahme neuer Variablen) übersetzt werden kann.

Erste Beispiele

Beispiel

Wir betrachten ein konstantes Vektorfeld auf dem $ℝ^{n}$ , also eine Abbildung

ℝ \times ℝ^{n} ⟶ ℝ^{n}, (t, x_{1}, \dots, x_{n}) ⟼ w,

wobei $w \in ℝ^{n}$ ein fixierter Vektor ist. Im „Windmodell“ bedeutet dies, dass überall und zu jeder Zeit eine konstante Windgeschwindigkeit herrscht. Die Bewegung eines (durch den Wind getragenen) Teilchens muss sich also auf der durch einen Startpunkt und den Richtungsvektor $w$ gegebenen Geraden vollziehen. In der Tat besitzt das Anfangswertproblem

v^{'} = w und v (t_{0}) = P

die eindeutige^[3] (affin-lineare) Lösung

v : ℝ ⟶ ℝ^{n}, t ⟼ v (t) = P + (t - t_{0}) w,

wie man durch Ableiten bestätigt.

Beispiel

Wir betrachten ein stetiges ortsunabhängiges Vektorfeld auf dem $ℝ^{n}$ , d.h. es sei eine stetige Abbildung

g : I ⟶ ℝ^{n}

auf einem reellen Intervall $I$ gegeben, die wir als Vektorfeld

F : I \times ℝ^{n} ⟶ ℝ^{n}, (t, x_{1}, \dots, x_{n}) ⟼ g (t),

auffassen. Im „Windmodell“ bedeutet dies, dass zu einem festen Zeitpunkt überall die gleiche Windgeschwindigkeit herrscht, diese sich aber mit der Zeit ändert. Die Bewegungskurven der (durch den Wind getragenen) Teilchen müssen also parallel zueinander sein, also durch eine Ortsverschiebung auseinander hervorgehen. Der Differenzvektor zwischen den Positionen von zwei Teilchen bleibt während des Bewegungsvorgangs erhalten. Die Lösungskurven zu einem Anfangswertproblem

v^{'} = F (t, v) = g (t) und v (t_{0}) = P

lassen sich einfach berechnen: Die eindeutige Lösung ist die Integralkurve

v : ℝ ⟶ ℝ^{n}, t ⟼ v (t) = P + (\int_{t_{0}}^{t} g_{1} (s) d s, \int_{t_{0}}^{t} g_{2} (s) d s, \dots, \int_{t_{0}}^{t} g_{n} (s) d s),

wobei die $g_{i}$ die Komponentenfunktionen von $g$ sind.

Beispiel

Es seien $n$ reellwertige Funktionen $f_{i} (t, x)$ in zwei Variablen gegeben. Diese kann man zu einem Vektorfeld

F : ℝ \times ℝ^{n} ⟶ ℝ^{n}, (t, x_{1}, \dots, x_{n}) ⟼ (f_{1} (t, x_{1}), f_{2} (t, x_{2}), \dots, f_{n} (t, x_{n}))

zusammenfassen. Dabei hängt die $i$ -te Koordinatenfunktion des Vektorfeldes nur von $t$ und der $i$ -ten Ortskoordinaten $x_{i}$ ab. Eine Lösungskurve $x (t) = (x_{1} (t), \dots, x_{n} (t))$ muss die Bedingungen

x'_{i} (t) = F_{i} (t, x_{1}, \dots, x_{n}) = f_{i} (t, x_{i})

(für $i = 1, \dots, n$ ) erfüllen. Diese $n$ Bedingungen sind unabhängig voneinander, d.h. man kann die $n$ Komponentenfunktionen $x_{i} (t)$ getrennt mit einem eindimensionalen Ansatz bestimmen. Daher spricht man von einem entkoppelten Differentialgleichungssystem.

Manchmal ist ein Differentialgleichungssystem in den ursprünglich gegebenen Koordinaten nicht entkoppelt, lässt sich aber durch einen Koordinatenwechsel entkoppeln und dann lösen. Dies ist vor allem für lineare Differentialgleichungssysteme mit konstanten Koeffizienten wichtig, die mit Mitteln der linearen Algebra entkoppelt werden können.

Beispiel

Wir betrachten das (zeitunabhängige) Vektorfeld

F : ℝ \times ℝ^{2} ⟶ ℝ^{2}, (t, x, y) ⟼ (- y, x) .

Hier steht also der Richtungsvektor $F (t, x, y) = (- y, x)$ stets senkrecht auf dem Ortsvektor $(x, y)$ , und ihre Normen stimmen überein. Man erwartet kreisförmige Bewegungen. In der Tat ist zur Anfangsbedingung $v (0) = (r, 0)$ die Kurve

v : ℝ ⟶ ℝ^{2}, t ⟼ (r \cos t, r \sin t),

die eindeutige Lösung.

Definition

Es sei $U \subseteq V$ eine offene Teilmenge in einem endlichdimensionalen reellen Vektorraum $V$ , $I \subseteq ℝ$ ein Intervall und es sei

g : I \times U ⟶ ℝ, (t, v) ⟼ g (t, v)

eine Funktion. Dann heißt das Vektorfeld

F : I \times U ⟶ V, (t, v) ⟼ F (t, v) = g (t, v) \cdot v,

ein Zentralfeld.

Bei einem Zentralfeld sind also der Ortsvektor und der Richtungsvektor linear abhängig, d.h. der Richtungsvektor weist in Richtung des Ortsvektors. Daher findet die durch ein Zentralfeld definierte Bewegung allein auf der durch einen Ortspunkt und den Nullpunkt (dem Zentrum) festgelegten Geraden statt. Es handelt sich also im Grunde um einen eindimensional festgelegten Bewegungsvorgang, was auch im folgenden Lemma zum Ausdruck kommt.

Lemma

Es sei $U \subseteq V$ eine offene Teilmenge in einem endlichdimensionalen reellen Vektorraum $V$ . Es sei

F : I \times U ⟶ V, (t, v) ⟼ F (t, v) = g (t, v) \cdot v,

ein stetiges Zentralfeld zur stetigen Funktion

g : I \times U ⟶ ℝ, (t, v) ⟼ g (t, v) .

Es sei $w \in U$ und es sei

α : J ⟶ ℝ

eine Lösung der eindimensionalen Differentialgleichung

z^{'} = h (t, z) := g (t, z w) \cdot z mit α (t_{0}) = 1.

Dann ist

$v (t) = α (t) \cdot w$

eine Lösung des Anfangswertproblems

$v^{'} = F (t, v) mit v (t_{0}) = w .$

Beweis

Es ist

\begin{matrix} v^{'} (t) & = (α (t) \cdot w)^{'} \\ = α^{'} (t) \cdot w \\ = g (t, α (t) \cdot w) \cdot α (t) \cdot w \\ = F (t, α (t) \cdot w) \\ = F (t, v (t)) \end{matrix}

und

v (t_{0}) = α (t_{0}) \cdot w = w,

sodass eine Lösung des Anfangswertproblems vorliegt.

◻

Beispiel

Wir betrachten das Zentralfeld zur Funktion

g : ℝ \times ℝ \times ℝ ∖ {0} ⟶ ℝ, (t, x, y) ⟼ g (t, x, y) = \frac{t^{2} x^{2}}{y},

also das Vektorfeld

F : ℝ \times ℝ \times ℝ ∖ {0} ⟶ ℝ^{2}, (t, x, y) ⟼ \frac{t^{2} x^{2}}{y} \cdot (x, y) = (\frac{t^{2} x^{3}}{y}, t^{2} x^{2}),

und die Anfangsbedingung $v (0) = (4, - 3)$ . Um dieses Anfangswertproblem zu lösen, müssen wir gemäß Lemma 40.10 die eindimensionale gewöhnliche Differentialgleichung

z^{'} = g (t, 4 z, - 3 z) \cdot z = \frac{t^{2} 16 z^{2}}{- 3 z} \cdot z = - \frac{16}{3} t^{2} z^{2}

mit der Anfangsbedingung $z (0) = 1$ lösen. Dies ist eine Differentialgleichung mit getrennten Variablen, nach Korollar 30.6 ist

z (t) = \frac{1}{\frac{16}{9} t^{3} + 1}

die Lösung mit $z (0) = 1$ . Daher ist

v (t) = \frac{1}{\frac{16}{9} t^{3} + 1} (4, - 3)

die Lösung des Anfangswertproblems zum Zentralfeld.

Fußnoten

↑ Rein formal gesehen ist hier auch das leere Intervall zugelassen, wobei diese „leere Lösung“ natürlich uninteressant ist. Bei einem Anfangswertproblem sichert bereits die Anfangsbedingung, dass die Lösung nicht leer ist.
↑ Die physikalische Interpretation eines Vektorfeldes als Kraftfeld ist hingegen bei Wegintegralen (nämlich als Arbeitsintegral) richtig.
↑ Ob die Lösung einer Differentialgleichung (existiert und) eindeutig ist, ist ein wichtiges Problem. Der wichtigste Satz zu dieser Fragestellung ist der Satz von Picard-Lindelöf, den wir später besprechen werden. In vielen der hier besprochenen Beispiele ist die Eindeutigkeit der Lösung direkt klar oder folgt aus den Eindeutigkeitsaussagen aus dem ersten Semester.

<< | Kurs:Analysis (Osnabrück 2013-2015)/Teil II | >>

PDF-Version dieser Vorlesung

Arbeitsblatt zur Vorlesung (PDF)

[1] Rein formal gesehen ist hier auch das leere Intervall zugelassen, wobei diese „leere Lösung“ natürlich uninteressant ist. Bei einem Anfangswertproblem sichert bereits die Anfangsbedingung, dass die Lösung nicht leer ist.

[2] Die physikalische Interpretation eines Vektorfeldes als Kraftfeld ist hingegen bei Wegintegralen (nämlich als Arbeitsintegral) richtig.

[3] Ob die Lösung einer Differentialgleichung (existiert und) eindeutig ist, ist ein wichtiges Problem. Der wichtigste Satz zu dieser Fragestellung ist der Satz von Picard-Lindelöf, den wir später besprechen werden. In vielen der hier besprochenen Beispiele ist die Eindeutigkeit der Lösung direkt klar oder folgt aus den Eindeutigkeitsaussagen aus dem ersten Semester.

[1]

[2]

[3]