Prof. Dr. Guido Sweers WS 20/21 Inka Schnieders, M.Sc.
Gewöhnliche Differentialgleichungen
Übungsblatt 5
Die Lösungen müssen eingescannt über Ilias eingereicht werden. Sollten dabei Probleme auf- treten melden Sie sich bei Inka Schnieders. Abgabeschluss ist am Dienstag, den 08.12.2020, um 12 Uhr.
Aufgabe 1: Gegeben sei das Anfangswertproblem
y
0(x) =
1 2 0 x
2y(x) + 3x
0
, y(0) = 0
1
.
a) Berechnen Sie die Picard-Iterationen y
1, y
2und y
3, wenn Sie mit der konstanten Funktion
y
0(x) = 0
1
beginnen.
b) Verwenden Sie das Euler-Vorwärts-Verfahren, um eine Näherung für die exakte Lösung in x =
32zu bekommen. Dabei soll die Schrittweite h =
12verwendet werden.
Lösung 1: a) Wir haben die Funktion
y
0(x) = 0
1
gegeben. Wir finden
y
n+1(x) = y
0(x) + Z
x0
1 2 0 s
2y
n(s) + 3s
0
ds.
Also gilt
y
1(x) = 0
1
+ Z
x0
1 2 0 s
20 1
+ 3s
0
ds
= 0
1
+ Z
x0
2 + 3s s
2ds =
0 1
+
2x +
32x
21 3
x
3=
2x +
32x
21 3
x
3+ 1
. Und analog finden wir auch
y
2(x) = 0
1
+ Z
x0
1 2 0 s
22s +
32s
21 3
s
3+ 1
+
3s 0
ds
= 0
1
+ Z
x0
5s +
32s
2+
23s
3+ 2
1
3
s
5+ s
2ds =
52
x
2+
12x
3+
16x
4+ 2x
1
18
x
6+
13x
3+ 1
.
1
Genauso folgt auch
y
3(x) = 0
1
+ Z
x0
1 2 0 s
25
2
s
2+
12s
3+
16s
4+ 2s
1
18
s
6+
13s
3+ 1
+ 3s
0
ds
= 0
1
+ Z
x0
52
s
2+
76s
3+
16s
4+ 5s +
19s
6+ 2
1
18
s
8+
13s
5+ s
2ds
=
56
x
3+
247x
4+
301x
5+
52x
2+
631x
7+ 2x
1
162
x
9+
181x
6+
13x
3+ 1
.
b) Wir berechnen einige diskreten Stellen der Approximation. Wir finden
y
n+1= y
n+ h
1 2 0 x
2ny
n+ h 3x
n0
und damit
x
0= 0 y
0= 0
1
, x
1=
12y
1=
1 1
, x
2= 1 y
2=
1349 8
, x
3=
32y
3=
15272 16
.
Aufgabe 2 (6+4 Punkte): Wir betrachten das Anfangswertproblem ( (
13+ x
2)y
00(x) − 4xy
0(x) + 6y(x) = 0,
y(0) = 1, y
0(0) = 2.
a) Nehmen Sie an, dass sich die Lösung als Potenzreihe schreiben lässt und finden Sie so eine Lösung des Anfangswertproblems.
b) Zeigen Sie, dass das Problem keine andere Lösung haben kann.
Hinweis: Schreiben Sie das Problem in ein System erster Ordnung um und verwenden Sie ein Ergebnis zur Eindeutigkeit von Lösungen.
Lösung 2: a) Wir nehmen an, dass man die Lösung als Potenzreihe schreiben kann. Es gilt also
y(x) =
∞
X
k=0
b
kx
k.
Wegen der Anfangsbedingungen muss b
0= 1 und b
1= 2 gelten. Außerdem finden wir y
0(x) =
∞
X
k=0
b
k+1(k + 1)x
k,
y
00(x) =
∞
X
k=0
b
k+2(k + 2)(k + 1)x
k. Mit der Differentialgleichung folgt
(
13+ x
2)
∞
X
k=0
b
k+2(k + 2)(k + 1)x
k− 4x
∞
X
k=0
b
k+1(k + 1)x
k+ 6
∞
X
k=0
b
kx
k= 0
∞
X
k=0
x
k 13(k + 2)(k + 1)b
k+2+ 6b
k+ x
k+1(−4(k + 1)b
k+1) + x
k+2((k + 2)(k + 1)b
k+2)
= 0.
Dies können wir umschreiben und finden
0 = (
23b
2+ 6b
0) + x(2b
1+ 2b
3) +
∞
X
k=2
x
k 13(k + 2)(k + 1)b
k+2+ b
k(k
2− 5k + 6) Wir können b
0= 1 und b
1= 2 einsetzen und erhalten
0 = (
23b
2+ 6) + x(4 + 2b
3) +
∞
X
k=2
x
k 13(k + 2)(k + 1)b
k+2+ b
k(k
2− 5k + 6) Nun können wir die Koeffizienten berechnen. Wir finden
b
2= −9, b
3= −2, b
4= 0, b
5= 0,
und damit folgt auch b
k= 0 für alle k ≥ 6. Wir erhalten also y(x) = 1 + 2x − 9x
2− 2x
3.
3
b) Wir können das Problem in ein System erster Ordnung umschreiben. Mit z
1(x) = y(x) und z
2(x) = y
0(x) finden wir
z
10(x) z
20(x)
= z
2(x)
4x
1
3+x2
z
2(x) −
1 63+x2
z
1(x)
!
= 0 1
−
1 6 3+x24x
1 3+x2
! z
1(x) z
2(x)
Das bedeutet, dass wir das folgende Anfangswertproblem betrachten:
z
10(x) z
20(x)
!
= 0 1
−
1 6 3+x24x
1 3+x2
! z
1(x) z
2(x)
! , z
1(0)
z
2(0)
!
= 1 2
! .
Nach Lemma 3.5 hat das Problem höchstens eine Lösung, da die Matrix beschränkt ist.
Aufgabe 3 (10 Punkte): Sei {y
n}
n∈Ndie Picard-Iteration zu dem folgenden Anfangswertpro- blem:
( y
0(x) = xy(x) + 1, y(0) = 0.
Berechnen Sie genügend viele Iterationsschritte, um eine Formel für y
nzu erkennen. Beweisen Sie diese dann durch vollständige Induktion.
Lösung 3: Wir setzen
y
0(x) = 0.
Dann finden wir im ersten Iterationsschritt
y
1(x) = 0 + Z
x0
s0 + 1ds = Z
x0
1ds = x.
Weiter folgt
y
2(x) = Z
x0
(sy
1(s) + 1)ds = Z
x0
(s
2+ 1)ds =
13x
3+ x, y
3(x) =
Z
x 0(sy
2(s) + 1)ds = Z
x0
(
13s
4+ s
2+ 1)ds =
5·31x
5+
13x
3+ x, y
4(x) =
Z
x 0(sy
3(s) + 1)ds = Z
x0
(
151s
6+
13s
4+ s
2+ 1)ds =
7·5·31x
7+
5·31x
5+
13x
3+ x.
Wir vermuten, dass
y
n(x) =
n−1
X
k=0
x
2k+11 Q
kj=0
(2j + 1) für n ∈ N
+.
Für n = 1 stimmt die Gleichung. Wir nehmen also an, dass die Formel für ein n ∈ N
+stimmt.
Wir zeigen, dass sie dann auch für n + 1 korrekt ist. Es gilt y
n+1(x) =
Z
x 0(sy
n(s) + 1)ds = Z
x0
1 +
n−1
X
k=0
s
2(k+1)1
Q
kj=0
(2j + 1)
! ds
= x +
n−1
X
k=0
1
2(k + 1) + 1 x
2(k+1)+11 Q
kj=0
(2j + 1)
= x +
n−1
X
k=0
x
2(k+1)+11
Q
k+1j=0
(2j + 1)
= x +
n
X
k=1
x
2k+11 Q
kj=0
(2j + 1)
=
n
X
k=0
x
2k+11 Q
kj=0
(2j + 1) .
Damit haben wir gezeigt, dass die Formel für alle n ∈ N
+gültig ist.
5
Bemerkung: Mathematica bekommt als Lösung des Anfangswertproblems y(x) =
r π
2 e
x2/2erf
√x 2
heraus.
-3 -2 -1 1 2 3
-100 -50 50 100
In der Abbildung ist die Lösung in blau und die Approximationen y
1(orange), y
2(grün), y
3(rot) und y
4(lila) dargestellt.
Aufgabe 4: Wir betrachten die folgenden Funktionen f
imit i ∈ {1, 2}:
f
1: ( R , | · |) → ( R , | · |), f
1(x) = x + e
x2,
f
2: (C([0, 1]), k · k
∞) → (C([0, 1]), k · k
∞), (f
2(x))(t) =
13x(t
2) − 1.
Beantworten Sie mit Begründung die folgenden Fragen:
a) Gilt kf
i(x) − f
i(y)k < kx − yk für alle x, y aus den angegebenen Definitionsgebieten mit x 6= y?
b) Gibt es eine Konstante L > 0, sodass kf
i(x) − f
i(y)k ≤ Lkx − yk für alle x, y aus den angegebenen Definitionsgebieten?
c) Wie viele Fixpunkte hat f
i?
Lösung 4: a) 1. Für f
1ist die Ungleichung nicht erfüllt. Beispielsweise kann x = 0 und y = 1 betrachtet werden. Es gilt
|f
1(0) − f
1(1)| = |1 − 1 − e
1| = e
1> 1 = |x − y|.
2. Da t 7→ t
2auf [0, 1] bijektiv ist, finden wir für beliebige x, y ∈ C([0, 1]) kf
2(x) − f
2(y)k
∞=
13sup
t∈[0,1]
|x(t
2) − y(t
2)| =
13sup
t∈[0,1]
|x(t) − y(t)| =
13kx − yk
∞. Die Bedingung ist also erfüllt.
b) 1. Mit dem Mittelwertsatz finden wir für jedes Paar x < y ein ξ
x,y∈ (x, y), sodass
|f
1(x) − f
1(y)| = |f
10(ξ
x,y)||x − y| = |1 + 2ξ
x,ye
ξx,y2||x − y|.
Wählt man jedoch x
n= n und y
n= n + 1, finden wir, dass |f
10(ξ
xn,yn)| → ∞ für n → ∞. Also gilt die Bedingung nicht.
2. Ja, mit L =
13folgt die Aussage wie in Aufgabenteil a).
c) 1. Für einen Fixpunkt x ∈ R muss gelten f
1(x) = x, also x = f
1(x) = x + e
x2⇔ 0 = e
x2. Da jedoch e
x2> 0 für alle x ∈ R , hat f
1keinen Fixpunkt.
2. Für einen Fixpunkt x ∈ C([0, 1]) muss iterativ für n ∈ N gelten:
x(t) =
13x(t
2) − 1 =
312x(t
4) −
13− 1
=
313x(t
8) −
312−
13− 1
=
3n+11x(t
2n+1) −
n
X
k=0 1 3k
. Da x beschränkt ist, folgt für n → ∞, dass
x(t) = −
∞
X
k=0 1
3k