Kurs:Analysis II/Kapitel VI: Gewöhnliche Differentialgleichungen/Differenzierbare Abhängigkeit von den Anfangswerten (§6)
Aus Wikiversity
Inhaltsverzeichnis |
[Bearbeiten] Voraussetzung (c):
Seien die stetigen Funktionen
mit
in Rfür
auf dem Rechteck R gegeben. Deren folgende partielle Ableitungen
für 
existieren und dort stetig sind.
[Bearbeiten] Bemerkung
Die Menge R ist kompakt und konvex und somit liefert der Mittelwertsatz der Differentialrechnung eine Lipschitzbedingung gemäß Voraussetzung (b) mit

als Lipschitzkonstante.
Wir wählen nun die Anfangswerte
![\begin{matrix} \eta = (\eta_1, \ldots, \eta_m) \text{ und zu festem Index } j \in \{1, \ldots, m\} \\ \overline{\eta} = (\overline{\eta}_1, \ldots, \overline{\eta}_m) := (\eta_1, \ldots, \eta_j + \lambda, \ldots, \eta_m) = \Bigl( \eta_k + \lambda \cdot \delta_{kj} \Bigr)_{k = 1, \ldots, m} \\ \text{mit dem Parameter } \lambda \in [- \varepsilon, 0) \cup (0, \varepsilon]\ \mathrm{f\ddot ur\ ein\ gegebenes}\ \varepsilon > 0. \end{matrix}](http://upload.wikimedia.org/math/b/8/8/b8837c9ffbc8cd918b0ee082c4218e2e.png)
Wir betrachten die zugehörigen Lösungen (2) und (3) der Anfangswertprobleme aus §5 und wir untersuchen die Differenzenquotienten
![\begin{matrix} z_{ij}(x, \lambda) := \frac{\overline{y}_i(x) - y_i(x)}{\lambda}, \quad x \in [\xi - h, \xi + h] \\ \mathrm{f\ddot ur}\ \lambda \in [- \varepsilon, 0) \cup (0, \varepsilon] \text{ und } i, j = 1, \ldots, m. \end{matrix}](http://upload.wikimedia.org/math/8/b/3/8b3eeef5db20cdd0565482dda3b152e0.png)
Wir verwenden die äquivalenten Integralgleichungssysteme (4) sowie (5) aus §5 und wir erhalten die folgende Integralgleichung für die Differenzenquotienten
![\begin{matrix} z_{ij}(x, \lambda) = \frac{\overline{y}_i(x) - y_i(x)}{\lambda} \\ = \delta_{ij} + \int^x_\xi \frac{1}{\lambda} \{f_i(t, \overline{y}_1(t), \ldots, \overline{y}_m(t)) - f_i(t, y_1(t), \ldots, y_m(t))\}\, dt \\ = \delta_{ij} + \int^x_\xi \left\{ \sum\limits^m_{k = 1} \frac{\partial f_i}{\partial y_k} \Bigl( t, y(t) + \tau_i(t) (\overline{y}(t) - y(t)) \Bigr) z_{kj}(t, \lambda) \right\}\, dt \\ = \delta_{ij} + \int^x_\xi \left\{ \sum\limits^m_{k = 1} p_{ik}(t, \lambda) \cdot z_{kj}(t, \lambda) \right\}\, dt \\ \mathrm{f\ddot ur}\ x \in [\xi - h, \xi + h], \lambda \in [- \varepsilon, 0) \cup (0, \varepsilon] \text{ und } i, j = 1, \ldots, m. \end{matrix}](http://upload.wikimedia.org/math/3/f/d/3fd2e12399e02b2c83b24b492c821a7d.png)
Hierbei erfüllen die Zwischenwerte die Bedingung
für
. Die Koeffizienten
![p_{ik}(x, \lambda) := \frac{\partial f_i}{\partial y_k} \Bigl( x, y(x) + \tau_i(x) (\overline{y}(x) - y(x)) \Bigr), x \in [\xi - h, \xi + h], \lambda \in [- \varepsilon, 0) \cup (0, \varepsilon]](http://upload.wikimedia.org/math/9/c/4/9c4d922970213a8d6e7195afd99b1b7e.png)
besitzen den Grenzwert
![\lim_{\lambda \to 0, \lambda \neq 0} p_{ik}(x, \lambda) = \frac{\partial f_i}{\partial y_k} \Bigl( x, y(x) \Bigr), \quad x \in [\xi - h, \xi + h]](http://upload.wikimedia.org/math/9/b/b/9bbd3604b21c34960a081a58d27dbd45.png)
für
. Wir haben also das lineare, parameterabhängige System
![\begin{matrix} z_{ij}(x, \lambda) = \delta_{ij} + \int^x_\xi \left\{ \sum\limits^m_{k = 1} p_{ik}(t, \lambda) \cdot z_{kj}(t, \lambda) \right\}\, dt \\ \mathrm{f\ddot ur}\ x \in [\xi - h, \xi + h], \lambda \in [- \varepsilon, 0) \cup (0, \varepsilon] \text{ und } i, j = 1, \ldots, m \end{matrix}](http://upload.wikimedia.org/math/4/c/8/4c886afa15877aa0c650475dd32ead84.png)
mit den Koeffizientenfunktionen (4) unter der asymptotischen Bedingung (5) zu studieren.
Zunächst zeigen wir, dass die Lösungen der Integralgleichung (6) gleichmäßig beschränkt sind. Hierzu formen wir diese um in die Gestalt
![\begin{matrix} z_{ij}(x, \lambda) - \delta_{ij} = \int^x_\xi \left\{ \sum\limits^m_{k = 1} p_{ik}(t, \lambda) \cdot \Bigl( z_{kj}(t, \lambda) - \delta_{kj} \Bigr) + p_{ij}(t, \lambda) \right\}\, dt \\ \mathrm{f\ddot ur}\ x \in [\xi - h, \xi + h], \lambda \in [- \varepsilon, 0) \cup (0, \varepsilon] \text{ und } i, j = 1, \ldots, m. \end{matrix}](http://upload.wikimedia.org/math/5/d/5/5d564348168c5e7be825f7b3bd906f5d.png)
Dann betrachten wir für
die Hilfsfunktion
![\Psi_j(x) := \sum^m_{i = 1} |z_{ij}(t, \lambda) - \delta_{ij}|, x \in [\xi - h, \xi + h], \lambda \in [- \varepsilon, 0) \cup (0, \varepsilon]](http://upload.wikimedia.org/math/e/2/6/e26b39e661879c24ec35e6fe443daaeb.png)
und entnehmen (7) die folgenden Integralungleichungen:
für alle
.Das Gronwallsche Lemma impliziert die Abschätzung
für alle
und
.Somit sind die Lösungen von (6) gleichmäßig beschränkt.
Wir vergleichen jetzt die Lösungen von (6) zu zwei verschiedenen Parametern
und erhalten für
die Identität
![\begin{matrix} z_{ij}(x, \lambda) - z_{ij}(x, \mu) = \int^x_\xi \left\{ \sum\limits^m_{k = 1} p_{ik}(t, \lambda) \cdot \Bigl( z_{kj}(t, \lambda) - z_{kj}(t, \mu) \Bigr) \right. \\ \left. + \sum\limits^m_{k = 1} \Bigl( p_{ik}(t, \lambda) - p_{ik}(t, \mu) \Bigr) \cdot z_{kj}(t, \mu) \right\}\, dt, \quad x \in [\xi - h, \xi + h]. \end{matrix}](http://upload.wikimedia.org/math/1/8/a/18a9a844b94949d66761a152b9c84349.png)
Nun genügen für
die Funktionen
![\Phi_j(x) := \sum^m_{i = 1} |z_{ij}(t, \lambda) - z_{ij}(t, \mu)|, \quad x \in [\xi - h, \xi + h]](http://upload.wikimedia.org/math/c/c/9/cc9c99e7a9e64e7f10a1040018d7ac5a.png)
wegen (11) Integralungleichungen der Form
.Nach (10) sind nämlich die Funktionen zkj(t,μ) gleichmäßig beschränkt und zusammen mit der asymptotischen Bedingung (5) erhalten wir
.Wenden wir nun das Gronwallsche Lemma auf die Ungleichung (13) an, so existieren für
die Limites
.Hierbei fassen wir y = y(x,η) als Funktion ihrer Anfangswerte auf. Wir erklären nun die Matrix
.Dann erhalten wir beim Grenzübergang
in (6) das folgende Integralgleichungssystem
für
und
,worin der Anfangswert η als Parameter auftritt. Da nach Satz 2 aus §5 die Lösung y(x,η) bereits stetig von den Anfangswerten η abhängt, sind die Koeffizienten Q = Q(x,η) in (16) und (17) stetig von diesen Daten abhängig. Wie wir oben für den Differenzenquotienten gezeigt haben, so beweisen wir genauso für die Funktionen yij(x,η) als Lösungen der Integralgleichung (17) mit dem Gronwallschen Lemma, dass sie gleichmäßig beschränkt auf ihrem Definitionsintervall sind und stetig von den Anfangswerten abhängen.
[Bearbeiten] Satz 1 (Differenzierbare Abhängigkeit von den Anfangswerten)
- Unter der Voraussetzung (c) hängt die Lösung y(x,η) des Differentialgleichungssystems (2) aus §5 einmal stetig differenzierbar von den Anfangswerten η ab.
[Bearbeiten] Voraussetzung (d):
Seien die stetigen Funktionen
mit
in Rfür
auf dem Rechteck R gegeben. Deren folgende zweite partielle Ableitungen
für 
existieren und sind dort stetig.
Wir vergleichen jetzt die Lösung des Systems (17) zum Parameter η mit derjenigen zum verschobenen Parametervektor
und ![\lambda \in [- \varepsilon, 0) \cup (0, \varepsilon]](http://upload.wikimedia.org/math/6/0/a/60acf1638d5def7e076e9d9ea3e54221.png)
für ein beliebiges
. Dann definieren wir – parallel zu (2) – für
die Differenzenquotienten
mit
.Aus (17) erhalten wir für
die Integralgleichungen
![\begin{matrix} z_{ijl}(x, \lambda) = \int^x_\xi \left\{ \sum\limits^m_{k = 1} q_{ik}(t, \eta) \cdot z_{ijl}(x, \lambda) + \sum\limits^m_{k = 1} \frac{q_{ik}(t, \overline{\eta}) - q_{ik}(t, \eta)}{\lambda} \cdot y_{kj}(t, \overline{\eta}) \right\}\, dt, \\ x \in [\xi - h, \xi + h] \end{matrix}](http://upload.wikimedia.org/math/7/8/2/782e69c150093adf4b318572229bbebb.png)
mit dem Parameter
. Wegen der Voraussetzung (d) und der Setzung (16) können wir die in (20) auftretenden Differenzenquotienten über die Kettenregel mit gewissen Zwischenwerten
wie folgt angeben:
![\begin{matrix} \frac{q_{ik}(x, \overline{\eta}) - q_{ik}(x, \eta)}{\lambda} = \frac{1}{\lambda} \cdot \left( \frac{\partial f_i}{\partial y_k} (x, y(x, \overline{\eta})) - \frac{\partial f_i}{\partial y_k} (x, y(x, \eta)) \right) \\ = \sum\limits^m_{n = 1} \frac{\partial^2 f_i}{\partial y_k \partial y_n} (x, y(x, \eta) + \tau_{ik}(x) \cdot [y(x, \overline{\eta}) - y(x, \eta)]) \cdot z_{nl}(x, \lambda), \\ x \in [\xi - h, \xi + h]. \end{matrix}](http://upload.wikimedia.org/math/f/5/f/f5f404fedb5bd9679770bd81ba54452f.png)
Somit sind diese Differenzenquotienten beschränkt auf ihrem Definitionsintervall und wir berechnen für alle
ihre Grenzwerte
.Wir definieren nun für
die beschränkten Funktionen
![r_{ijl}(x, \lambda) := \sum^m_{k = 1} \frac{q_{ik}(x, \overline{\eta}) - q_{ik}(x, \eta)}{\lambda} \cdot y_{kj}(x, \overline{\eta}), \quad x \in [\xi - h, \xi + h]](http://upload.wikimedia.org/math/f/e/0/fe0474a1dd27218bc11979b4bbd74489.png)
und die Integralgleichungen (20) verwandeln sich in
![z_{ijl}(x, \lambda) = \int^x_\xi \left\{ \sum^m_{k = 1} q_{ik}(t, \eta) \cdot z_{kjl}(t, \lambda) + r_{ijl}(t, \lambda) \right\}\, dt, \quad x \in [\xi - h, \xi + h]](http://upload.wikimedia.org/math/3/1/d/31dc941df1ee836ecbabbb0df8d4933b.png)
Wie oben sehen wir mit dem Gronwallschen Lemma leicht ein, dass die Lösungen der Integralgleichungen (24) gleichmäßig beschränkt sind und für
die folgenden Grenzwerte existieren:
.Wir können jetzt in den Integralgleichungen (24) mit Hilfe von (22) sowie (23) den Grenzübergang
vollziehen. Mit
erhalten wir die Integralgleichungen des Tensors der zweiten Ableitungen
![\begin{matrix} y_{ijl}(x, \eta) = \int^x_\xi \left\{ \sum\limits^m_{k = 1} q_{ik}(t, \eta) \cdot y_{kjl}(t, \eta) + s_{ijl}(t, \eta) \right\}\, dt \\ \mathrm{f\ddot ur\ alle}\ x \in [\xi - h, \xi + h] \text{ mit den Funktionen} \\ s_{ijl}(x, \eta) = \sum\limits^m_{k, n = 1} \frac{\partial^2 f_i}{\partial y_k \partial y_n} \Bigl( x, y(x, \eta) \Bigr) \cdot y_{nl}(x, \eta) y_{kj}(x, \eta). \end{matrix}](http://upload.wikimedia.org/math/5/5/7/5573c1640d8cb4bec29f3f01909d6fc1.png)
Wir haben in (26) ein inhomogenes, lineares System mit 0 Anfangsbedingungen vor uns, worin die rechte Seite stetig vom Parameter η abhängt. Wir ziehen die Integralgleichungen für einen festen Anfangswert
ab von derjenigen für einen benachbarten Anfangswert η und erhalten
![\begin{matrix} y_{ijl}(x, \eta) - y_{ijl}(x, \stackrel{\circ}{\eta}) = \int^x_\xi \left\{ \sum\limits^m_{k = 1} q_{ik}(t, \eta) \cdot \left( y_{kjl}(t, \eta) - y_{kjl}(t, \stackrel{\circ}{\eta}) \right) \right\}\, dt \\ + \int^x_\xi \left\{ \sum\limits^m_{k = 1} \left( q_{ik}(t, \eta) - q_{ik}(t, \stackrel{\circ}{\eta}) \right) \cdot y_{kjl}(t, \stackrel{\circ}{\eta}) + \left( s_{ijl}(t, \eta) - s_{ijl}(t, \stackrel{\circ}{\eta}) \right) \right\}\, dt \\ \mathrm{f\ddot ur\ alle}\ x \in [\xi - h, \xi + h]. \end{matrix}](http://upload.wikimedia.org/math/c/e/5/ce5fae24d4495a64847d98a4fbf48608.png)
Für
betrachten wir die Hilfsfunktion
.Dann leiten wir aus (27) eine Differentialungleichung der Form
![\Theta_{jl}(x) \le A \cdot \int^x_\xi \left\{ \Theta_{jl}(t) + \varepsilon(\eta, \stackrel{\circ}{\eta}) \right\}\, dt, \quad x \in [\xi - h, \xi + h]](http://upload.wikimedia.org/math/e/4/3/e438c716d4313bc480ac3af7288e3faa.png)
her; dabei ist
eine Konstante und die Größe
besitzt die asymptotische Eigenschaft
für
. Mit dem Gronwallschen Lemma zeigen wir nun, dass
gleichmäßig auf dem Intervall [ξ − h,ξ + h] für
erfüllt ist. Somit ist der Tensor der zweiten Ableitungen
stetig von den Anfangswerten η abhängig.
[Bearbeiten] Satz 2 (C2-Abhängigkeit von den Anfangswerten)
- Unter der Voraussetzung (d) hängt die Lösung y(x,η) des Differentialgleichungssystems (2) aus §5 zweimal stetig differenzierbar von den Anfangswerten η ab.
[Bearbeiten] Satz 3 (Existenz des integrierenden Faktors)
- Auf der Umgebung
des Punktes
schreiben wir die Funktionen p = p(x,y) und q = q(x,y) der Klasse
mit der Eigenschaft

- vor. Dann gibt es eine Umgebung
von (x0,y0) und Funktionen

- sowie
,- so dass folgendes gilt:
und
in
.- Wir erhalten also mit M einen Eulerschen Multiplikator der regulären, ebenen Differentialgleichung
in
.[Bearbeiten] Beweis
Wir können ohne Einschränkung
in
annehmen. Dann transformieren wir

in die explizite Differentialgleichung
.Nun gehören p,q,f zur Klasse C2 und f erfüllt insbesondere eine Lipschitzbedingung. Für ein hinreichen klein gewähltes
lösen wir für alle Parameter
das parameterabhängige Anfangswertproblem
.Auf dem Rechteck

betrachten wir die Transformation
.Wegen der Eigenschaft
gehört sie zur Klasse C2 und erfüllt die Bedingung
.Nach dem Fundamentalsatz über die inverse Abbildung existiert die Umkehrabbildung auf einer gewissen Umgebung
von (x0,y0), nämlich
.Wir setzen
mit
.Nach Konstruktion ist klar, dass die Niveaulinien F = const die Differentialgleichung
lösen. Da
senkrecht auf den Niveaulinien steht, folgt die Identität
in 
mit einer Funktion

q.e.d.