Instytut Automatyki Politechniki Łódzkiej - Metody Numeryczne wykład 11
W11-1
Obszary stabilności absolutnej
Adamsa-Bashfortha
l. kroków
1
2
3
4
rząd 1
2
3
4
q
*
-2 -1 -6/11
-3/10
Adamsa-Moultona
l. kroków
1
2
3
4
rząd 2
3
4
5
q
*
-
∞
-6 -3 -90/49
Instytut Automatyki Politechniki Łódzkiej - Metody Numeryczne wykład 11
W11-2
Wstecznego różniczkowania
k 1
2
3
4
5
6
a 0
0
-0.1
-0.7
-2.4 -6.1
α
90 90 88 73 51 18
Instytut Automatyki Politechniki Łódzkiej - Metody Numeryczne wykład 11
W11-3
Instytut Automatyki Politechniki Łódzkiej - Metody Numeryczne wykład 11
W11-4
Rungego-Kutty
p=m 1
2
3 4
Lewa
granica
obszaru
stabilności
-2 -2 -2.51 -2.78
Instytut Automatyki Politechniki Łódzkiej - Metody Numeryczne wykład 11
W11-5
Przykłady
Rozwiązanie wolnozmienne
10
2
1
3
)
(
,
50
0
,
1
)
0
(
3
)
(
1
)
(
1
.
0
)
(
t
e
t
f
t
f
t
f
t
f
dt
t
df
−
+
=
<
<
=
−
=
Instytut Automatyki Politechniki Łódzkiej - Metody Numeryczne wykład 11
W11-6
10
-12
10
-10
10
-8
10
-6
10
-4
10
-2
10
-14
10
-12
10
-10
10
-8
10
-6
10
-4
10
-2
Błąd w funkcji RelTol
RKF23
RKF45
Adam s PC
Gear W R
Instytut Automatyki Politechniki Łódzkiej - Metody Numeryczne wykład 11
W11-7
10
-12
10
-10
10
-8
10
-6
10
-4
10
-2
10
1
10
2
10
3
10
4
Kroki w funkcji RelTol
RKF23
RKF45
Adam s PC
Gear W R
Instytut Automatyki Politechniki Łódzkiej - Metody Numeryczne wykład 11
W11-8
10
-12
10
-10
10
-8
10
-6
10
-4
10
-2
10
1
10
2
10
3
10
4
10
5
LW F w funkcji RelTol
RKF23
RKF45
Adam s PC
Gear W R
Instytut Automatyki Politechniki Łódzkiej - Metody Numeryczne wykład 11
W11-9
Równanie sztywne
(
)
)
(
)
(
,
2
1
3
)
(
,
50
0
500
,
2
)
0
(
)
(
)
(
)
(
)
(
10
t
f
e
t
y
e
t
f
t
a
y
dt
t
df
t
f
t
y
a
dt
t
dy
at
t
+
=
+
=
<
<
=
=
+
−
−
=
−
−
Instytut Automatyki Politechniki Łódzkiej - Metody Numeryczne wykład 11
W11-10
0
5
10
15
20
25
30
35
40
45
50
1
1.2
1.4
1.6
1.8
2
2.2
2.4
2.6
2.8
3
Instytut Automatyki Politechniki Łódzkiej - Metody Numeryczne wykład 11
W11-11
10
-6
10
-5
10
-4
10
-3
10
-2
10
1
10
2
10
3
10
4
10
5
Udane kroki w funkcji RelTol
RKF23
RKF45
Adam s PC
Gear W R
Instytut Automatyki Politechniki Łódzkiej - Metody Numeryczne wykład 11
W11-12
10
-6
10
-5
10
-4
10
-3
10
-2
10
0
10
1
10
2
10
3
10
4
Kroki odrzucone w funkcji RelTol
RKF23
RKF45
Adam s PC
Gear W R
Instytut Automatyki Politechniki Łódzkiej - Metody Numeryczne wykład 11
W11-13
10
-6
10
-5
10
-4
10
-3
10
-2
10
1
10
2
10
3
10
4
10
5
LW F w funkcji RelTol
RKF23
RKF45
Adams PC
Gear W R