Beispiel - SAM - ETH ZÃ¼rich

Weitere Magazine

Empfehlungen

Info

2.2.3.1 Konsistenzordnung . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 172 2.2.3.2 Spektrale Konvergenz . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 185 2.3 Runge-Kutta-Verfahren . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 222 2.3.1 Konstruktion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 222 2.3.2 Konvergenz . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 237 2.4 Extrapolationsverfahren [8, Sect. 4.3] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 249 2.4.1 Der Kombinationstrick . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 249 2.4.2 Extrapolationsidee . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 252 2.4.3 Extrapolation von Einschrittverfahren . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 259 2.4.4 Lokale Extrapolations-Einschrittverfahren . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 265 2.4.5 Ordnungssteuerung . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 271 2.4.6 Extrapolation reversibler Einschrittverfahren . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 274 2.5 Splittingverfahren [16, Sect. 2.5] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 278 2.6 Schrittweitensteuerung [8, Kap. 5], [19, Sect. 2.8] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 287 3 Stabilität [8, Kap. 6] 320 3.1 Modellproblemanalyse . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 322 3.2 Vererbung asymptotischer Stabilität . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 339 3.2.1 Attraktive Fixpunkte . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 339 3.2.2 Attraktive Fixpunkte von Einschrittverfahren . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 345 3.3 Nichtexpansivität [8, Abschn. 6.3.3] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 352 3.4 Gleichmässige Stabilität . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 363 3.5 Steifheit . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 375 3.6 Linear-implizite Runge-Kutta-Verfahren [8, Sect. 6.4] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 386 3.7 Exponentielle Integratoren [24, 28, 25] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 396 3.8 Differentiell-Algebraische Anfangswertprobleme . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 402 3.8.1 Grundbegriffe . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 402 3.8.2 Runge-Kutta-Verfahren für Index-1-DAEs . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 408 3.8.3 DAEs mit höherem Index . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 419 Numerische Mathemtik R. Hiptmair rev 35327, 13. Mai 2011 0.0 p. 4
4 Strukturerhaltende numerische Integration 433 4.1 Polynomiale Invarianten . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 434 Numerische Mathemtik 4.2 Volumenerhaltung . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 446 4.3 Verallgemeinerte Reversibilität . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 456 4.4 Symplektizität . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 466 4.4.1 Symplektische Evolutionen Hamiltonscher Differentialgleichungen . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 466 4.4.2 Symplektische Integratoren . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 483 4.4.3 Rückwärtsanalyse . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 507 4.4.4 Modifizierte Gleichungen: Fehleranalyse . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 519 4.4.5 Strukturerhaltende modifizierte Gleichungen . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 545 4.5 Methoden für oszillatorische Differentialgleichungen [23] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 556 Verzeichnisse 570 R. Hiptmair rev 35327, 13. Mai 2011 Stichwortverzeichnis . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 570 Verzeichnis der <strong>Beispiel</strong>e und Bemerkungen . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 581 Verzeichnis der Definitionen und Konzepte . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 585 Verzeichnis der MATLAB-CODE-Fragmente . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 587 Symbolverzeichnis . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 588 0.0 p. 5
Seite 1 und 2: Vorlesung 401-2654-00L, Numerische
Seite 3: 1.3.3.1 Grundbegriffe . . . . . . .
Seite 7 und 8: • MATLAB-basierte Programmieraufg
Seite 9 und 10: Hinweise auf Fehler in den Vorlesun
Seite 11 und 12: Numerische Mathemtik 1 Einleitung V
Seite 13 und 14: ✎ Notation (Newton): Punkt ˙ ˆ=
Seite 15 und 16: 1.5 1.5 Numerische Mathemtik 1 1 y
Seite 17 und 18: 5 6 % plot tangent field 7 f i g u
Seite 19 und 20: Verallgemeinerung: Eine gewöhnlich
Seite 21 und 22: Daher sind Darstellungsformeln wie
Seite 23 und 24: 1.2.1 Ökologie Numerische Mathemti
Seite 25 und 26: Bemerkung 1.2.4 (AWP-Löser in MATL
Seite 27 und 28: (1.2.6) ⇒ 0 =(δ − γ u )˙u−
Seite 29 und 30: Definition 1.2.7 (Erstes Integral).
Seite 31 und 32: Beispiel 1.2.12 (Oregonator-Reaktio
Seite 33 und 34: Grössen: l = l(t) ˆ= Länge der H
Seite 35 und 36: 5 6 f u n c t i o n beat(alpha,file
Seite 37 und 38: 9 axis([tspan -3 3]); legend(’l(t
Seite 39 und 40: Definition 1.2.20 (Hamiltonsche Dif
Seite 41 und 42: ←→ Hamiltonsches System (→ De
Seite 43 und 44: himpuls exakt, jedoch nicht die Ene
Seite 45 und 46: y (i) Numerische Mathemtik (iii) :
Seite 47 und 48: Ein einfaches Kriterium für lokale
Seite 49 und 50: Definition 1.3.7 (Evolutionsoperato
Seite 51 und 52: Lösung: ⎧ ⎨ y(t) = 1 y0 −1
Seite 53 und 54: Vorbereitung: Basiswechsel im Zusta
Seite 55 und 56:
Ansatz: y(t) = exp(At)z(t) mitz ∈
Seite 57 und 58:
1.3.3 Sensitivität [8, Sect. 3.1]
Seite 59 und 60:
wobei‖·‖ ˆ= Matrixnorm induzi
Seite 61 und 62:
Intervallweise Kondition: in[0,T]:
Seite 63 und 64:
Bemerkung 1.3.32 („Gronwall-Schra
Seite 65 und 66:
zum AWP (1.1.13) erfüllt Anfangswe
Seite 67 und 68:
Listing 1.3: Numerische Integration
Seite 69 und 70:
σ = 10, ρ = 28, β = 2.666667e+00
Seite 71 und 72:
m 1 = 2, m 2 = 1, l 1 = 1, l 2 = 1.
Seite 73 und 74:
1.4.1 Das explizite Euler-Verfahren
Seite 75 und 76:
Bemerkung 1.4.3 (Explizites Eulerve
Seite 77 und 78:
7 err = [err;erri]; % assemble matr
Seite 79 und 80:
Definition 1.4.5 (Arten der Konverg
Seite 81 und 82:
0.7 0.6 β = 0.5, γ = 0.5 β = 1.0
Seite 83 und 84:
10 0 λ = 1.000000 λ = 3.000000 λ
Seite 85 und 86:
Rekursion des expliziten Eulerverfa
Seite 87 und 88:
Bemerkung 1.4.14 (Implizites Eulerv
Seite 89 und 90:
( ) 1 k y k = y 0 ⇒ |y 1−λh k
Seite 91 und 92:
5 % Implicit Euler 6 y_imp = y0; y
Seite 93 und 94:
1 y = y_expl; 2 3 E_kin = 0.5*(y(2,
Seite 95 und 96:
6 E_tot = E_kin + E_pot; 7 8 f i g
Seite 97 und 98:
9 8 Energies for explicit Euler dis
Seite 99 und 100:
✸ Numerische Mathemtik 1.4.3 Impl
Seite 101 und 102:
Beispiel 1.4.21 (Implizite Mittelpu
Seite 103 und 104:
Beispiel 1.4.24 (Implizite Mittelpu
Seite 105 und 106:
Gegebeny k−1 ≈ y(t k−1 ),y k
Seite 107 und 108:
Berechney 1 aus (1.4.30) & (1.4.31)
Seite 109 und 110:
%f,\\alpha(0)=%f,p(0)=%f’,g,l,y0,
Seite 111 und 112:
(1.4.27) angewandt auf (1.2.18) 5 4
Seite 113 und 114:
10 9 Energies for Stoermer−Verlet
Seite 115 und 116:
y/v v k+ 1/2 Numerische Mathemtik P
Seite 117 und 118:
Numerische Mathemtik 2 Einschrittve
Seite 119 und 120:
Numerische Mathemtik Baustein: Verf
Seite 121 und 122:
Definition 2.1.5 (Explizite und imp
Seite 123 und 124:
Beachte: Konvergenz gemäss Def. 2.
Seite 125 und 126:
Definition 2.1.8 (Konsistenz einer
Seite 127 und 128:
Interpretation: D Ψ t,t+h y τ(t,y
Seite 129 und 130:
D(y) := x -> f(y(x)); D(y) := x ↦
Seite 131 und 132:
Annahme: rechte Seitef : Ω ↦→
Seite 133 und 134:
✬ ✩ Lemma 2.1.20 (Diskretes Gro
Seite 135 und 136:
➁ Annahme A1:(y k ) N k=0 existie
Seite 137 und 138:
Merkregel: (Nur) für Einschrittver
Seite 139 und 140:
Fehlerrekursion e k = y(t k )−y k
Seite 141 und 142:
Wir haben gesehen: Eine approximati
Seite 143 und 144:
Der Beweis verwendet folgendes Hilf
Seite 145 und 146:
✎ Notation: P s ˆ= Raum der univ
Seite 147 und 148:
➤ (2.2.3) ˆ= (Nichtlineares) Gle
Seite 149 und 150:
„Beweis” (von Lemma 2.2.7, unte
Seite 151 und 152:
✎ Übliche Notation für Koeffizi
Seite 153 und 154:
∥ ≤ |h|·‖A‖ ∞ max ∥f(g
Seite 155 und 156:
Beweis auf der Grundlage des Satzes
Seite 157 und 158:
Mit g y := (g y 1 ,...,gy s), g z :
Seite 159 und 160:
Verschobene äquidistante Kollokati
Seite 161:
• Falls = 1 & c 1 = 1/2 (↔ einf
Seite 164 und 165:
△ Numerische Mathemtik Bekannt au
Seite 166 und 167:
Bemerkung 2.2.25 (Fixpunktform von
Seite 168 und 169:
mit KompaktumK ⊂ D, für das (rü
Seite 170 und 171:
✬ ✩ Lemma 2.2.31 (Lipschitz-Ste
Seite 172 und 173:
2.2.3 Konvergenz von Kollokationsve
Seite 174 und 175:
Hilfsmittel beim Beweis: Restglieda
Seite 176 und 177:
Aus der Lösungseigenschaft vont
Seite 178 und 179:
10 0 h 1 Numerische Mathemtik y 0.8
Seite 180 und 181:
Bsp. 2.2.49 legt die Vermutung nahe
Seite 182 und 183:
Wegenẏ(t) = f(y(t)) folgt für di
Seite 184 und 185:
Geniale Idee: Abschätzung von ∫
Seite 186 und 187:
Dieses Beispiel studiert den Einsch
Seite 188 und 189:
Polynominterpolationsfehlerabschät
Seite 190 und 191:
1 f u n c t i o n errinf = gaussint
Seite 192 und 193:
Interpolationsfehleranalyse mit Hil
Seite 194 und 195:
✬ Lemma 2.2.62 (Residuenformel f
Seite 196 und 197:
f(t) = − s∑ f(τ j ) j=1 P(t) (
Seite 198 und 199:
Vermutung: |P n (t)| ≤ 1 für all
Seite 200 und 201:
Dabei wurde im Sinne der komplexen
Seite 202 und 203:
Wir betrachten nun die Joukowski-Tr
Seite 204 und 205:
ρ n /10 10 5 10 4 10 3 10 2 10 1 1
Seite 206 und 207:
Erinnerung: Formel von Cauchy-Hadam
Seite 208 und 209:
✬ ✩ Theorem 2.2.81 (Fehlerabsch
Seite 210 und 211:
Wenden wir ein Gauss-Kollokations-E
Seite 212 und 213:
Im 0.8 0.6 0.4 0.2 0 y 0 = 10,λ =
Seite 214 und 215:
λ = 5.5: Fehler zum Zeitpunkt T=1,
Seite 216 und 217:
Füry 0 < 1: Pole in−1− 2 λ ln
Seite 218 und 219:
400 350 300 y0=0.80111 λ=1 y0=0.80
Seite 220 und 221:
λ = 10: Fehler zum Zeitpunkt T=1,
Seite 222 und 223:
2.3 Runge-Kutta-Verfahren Numerisch
Seite 224 und 225:
Definition 2.3.5 (Runge-Kutta-Verfa
Seite 226 und 227:
Interpretation: Runge-Kutta-Verfahr
Seite 228 und 229:
1 Numerische Mathemtik grün: Expli
Seite 230 und 231:
Kuttas 3/8-Regel 1 Numerische Mathe
Seite 232 und 233:
Bemerkung 2.3.15 (Autonomisierungsi
Seite 234 und 235:
Bemerkung 2.3.18 (“Dense output
Seite 236 und 237:
h „klein” ➣ Natürliche Anfan
Seite 238 und 239:
1 0.9 0.8 0.7 y(t) Explicit Euler E
Seite 240 und 241:
Konsistenzfehler: τ(t,y,h) := (Φ
Seite 242 und 243:
k i =f(y 0 )+hDf(y 0 ) ⎛ 1 2 h 2
Seite 244 und 245:
s∑ b i c i = 1 2 , (2.3.30) i=1 s
Seite 246 und 247:
Numerische Mathemtik Bemerkung 2.3.
Seite 248 und 249:
Ziel: vorgegebene Reduktion des Feh
Seite 250 und 251:
(Annahme: Äquidistante Zeitschritt
Seite 252 und 253:
2.4.2 Extrapolationsidee Numerische
Seite 254 und 255:
Bemerkung 2.4.4 (Skalierungsinvaria
Seite 256 und 257:
☞ Extrapolation „funktioniert
Seite 258 und 259:
k∑ i=1 L i (0)h j i = ⎧ ⎪⎨
Seite 260 und 261:
✬ ✩ Theorem 2.4.11 (Asymptotisc
Seite 262 und 263:
Beweis von (2.4.15) durch Induktion
Seite 264 und 265:
= y(t+h)−y(t)−hψ(t,y(t)+e(t)h
Seite 266 und 267:
Fixiere Sequenz(n l ) k+1 l=1 ,n l
Seite 268 und 269:
AWP für logistische Dgl. (→ Bsp.
Seite 270 und 271:
wobeiC > 0 nur von den Verhältniss
Seite 272 und 273:
MATLAB-CODE : Adaptives Euler-Extra
Seite 274 und 275:
0.6 0.4 5 y (t ) (Naeherung) 1 k y
Seite 276 und 277:
Beobachtung aus Bsp. 2.4.20 einfach
Seite 278 und 279:
Bemerkung 2.4.23 (DIFEX). Numerisch
Seite 280 und 281:
Φ t f y = 1 1+(y −1 −1)e−λt
Seite 282 und 283:
Kutta-Inkremente in Bsp. 2.3.24, (2
Seite 284 und 285:
( ) 0 Splitting: F(y) = + a(r) } {{
Seite 286 und 287:
10 −2 Zeitschrittweite h |y(T)−
Seite 288 und 289:
Gedankenexperiment: wie wird sich w
Seite 290 und 291:
Concentration of HBrO 2 Fig. 89 4 3
Seite 292 und 293:
Effizienz ✬ Ziel: N so klein wie
Seite 294 und 295:
Ψ↔Explizites Euler-Verfahren (1.
Seite 296 und 297:
3 while ((t(end) < T) (h > hmin)) %
Seite 298 und 299:
! Gemäss unserer Heuristik, siehe
Seite 300 und 301:
☞ Adaptive Zeitschrittweitensteue
Seite 302 und 303:
Effizient ? Genaueres (teureres) Ve
Seite 304 und 305:
Anfangswertproblem für skalare log
Seite 306 und 307:
Beispiel 2.6.9 (“Versagen” adap
Seite 308 und 309:
Beispiel 2.6.10 (Schrittweitensteue
Seite 310 und 311:
Beispiel 2.6.12 (Schrittweitensteue
Seite 312 und 313:
0 1 3 1 3 1 2 1 3 1 1 6 6 1 8 0 3 8
Seite 314 und 315:
Beispiel 2.6.15 (Adaptive RK-ESV zu
Seite 316 und 317:
4 3 abstol = 0.001000, reltol = 0.0
Seite 318 und 319:
AWP aus Bsp. 2.4.19 ode45 mit versc
Seite 320 und 321:
Numerische Mathemtik 3 Stabilität
Seite 322 und 323:
Überlegung: Linearisierung um Fixp
Seite 324 und 325:
Dies ist eine Frage nach Strukturer
Seite 326 und 327:
Was sagt uns diese Stabilitätsfunk
Seite 328 und 329:
✬ ✩ Korollar 3.1.8. ✫ Explizi
Seite 330 und 331:
Dieses hat nur die triviale Lösung
Seite 332 und 333:
60 50 40 30 exp(z) Explicit Euler C
Seite 334 und 335:
3 2 * λ Stabilitätsbedingte Schri
Seite 336 und 337:
Hilfsmittel: Schur-Zerlegung Numeri
Seite 338 und 339:
Klar ist p(A) = s∑ c j A j ∑ f
Seite 340 und 341:
Die Begriffsbildung ist klar: eine
Seite 342 und 343:
✬ ✩ Theorem 3.2.4 (Hinreichende
Seite 344 und 345:
Linearisierung, siehe Bem. 1.3.19:
Seite 346 und 347:
Wir betrachten weiterhin ein autono
Seite 348 und 349:
✬ Numerische Mathemtik ✩ Lemma
Seite 350 und 351:
Für welche Verfahren entfällt Sch
Seite 352 und 353:
3 2 implizite Mittelpunktsregel (1.
Seite 354 und 355:
Erinnerung: Eine AbbildungV : R d
Seite 356 und 357:
Nichtexpansivität einer Evolution:
Seite 358 und 359:
Aus Kollokationsbedingungen (2.2.1)
Seite 360 und 361:
0.9 0.9 1 1 1 △ Numerische Mathem
Seite 362 und 363:
{ −y 3 füry < 0 , angewandt auf
Seite 364 und 365:
ODE mit stark attraktivem Fixpunkty
Seite 366 und 367:
➣ Ungenügende Dämpfung der Anfa
Seite 368 und 369:
Bemerkung 3.4.6 (Invertierbarkeit d
Seite 370 und 371:
order of quadrature rule 23 21 19 1
Seite 372 und 373:
1.5 1.5 0.7 0.4 10 20 30 8 6 4 2 0.
Seite 374 und 375:
✬ ✩ Theorem 3.4.8 (Radau-ESV ni
Seite 376 und 377:
Beispiel 3.5.2 (Adaptive explizite
Seite 378 und 379:
Fallsc A (0) > c B (0) ➢ 2. Reakt
Seite 380 und 381:
Schaltkreisanalyse im Zeitbereich:
Seite 382 und 383:
☞ Falls‖y 0 ‖ = 1 ⇒‖y(t)
Seite 384 und 385:
1 ode45 for rigid motion 0.2 Numeri
Seite 386 und 387:
3.6 Linear-implizite Runge-Kutta-Ve
Seite 388 und 389:
Logistic ODE, y 0 = 0.100000, λ =
Seite 390 und 391:
2-stufiges Radau-ESV, Butcher Schem
Seite 392 und 393:
Ein Newton-Schritt mit Startwertk (
Seite 394 und 395:
(3.6.9) ⇒ y 1 = y 0 +h ( ∑ s )
Seite 396 und 397:
3.7 Exponentielle Integratoren [24,
Seite 398 und 399:
exponentielles Euler-Verfahren (auf
Seite 400 und 401:
• Steifes AWP → Bsp. 3.5.2, 3.4
Seite 402 und 403:
Bemerkung 3.7.8. Herausforderung: e
Seite 404 und 405:
Numerische Mathemtik Beachte: ( 1 0
Seite 406 und 407:
In Bsp. 3.8.1: Transformationen ➣
Seite 408 und 409:
Definition 3.8.7 (DAE vom Index 1).
Seite 410 und 411:
( ) C+Cp −C )(˙u1 −C C ˙u2 =
Seite 412 und 413:
Formale Rechnung: Singuläre Störu
Seite 414 und 415:
⇒ v 1 = g v s ⇒ c(u 1 ,v 1 ) =
Seite 416 und 417:
Erwünscht: S(−∞) = 0 für Stab
Seite 418 und 419:
MATLAB-CODE : Lösung einer Index-1
Seite 420 und 421:
x 2 0 x 1 Mathematisches Pendel (Au
Seite 422 und 423:
Bemerkung 3.8.23 (Hamiltonsche Bewe
Seite 424 und 425:
FürH der Form (3.8.25): ∂ 2 H
Seite 426 und 427:
✸ Numerische Mathemtik Beispiel 3
Seite 428 und 429:
Man beobachtet algebraische Konverg
Seite 430 und 431:
➣ Stufengleichungen (→ Bem. 2.3
Seite 432 und 433:
△ Numerische Mathemtik Konvergenz
Seite 434 und 435:
Beachte: In diesem Kapitel beschrä
Seite 436 und 437:
Beweis. (für den autonomen Fall) L
Seite 438 und 439:
➣ Keine Erhaltung von‖y(t)‖
Seite 440 und 441:
y h (t) ∈ P s ˆ= Gauss-Kollokati
Seite 442 und 443:
g i = y 0 +h s ∑ j=1 a ij f(g j )
Seite 444 und 445:
☞ Bekannt aus der linearen Algebr
Seite 446 und 447:
S(A) = diag(S(µ),S(ν),S(−(µ+ν
Seite 448 und 449:
Beweis. (basierend auf Lemma 4.2.2,
Seite 450 und 451:
Volumenerhaltende numerische ODE-L
Seite 452 und 453:
so gilt detW(t) ≡ 1, also is detW
Seite 454 und 455:
⎛ ⎞ f 1 (y) f 2 (y) f 3 (y) f(y
Seite 456 und 457:
4.3 Verallgemeinerte Reversibilitä
Seite 458 und 459:
⇒ ⎧ s∑ k ⎪⎨ i = f(y 0 +h
Seite 460 und 461:
Numerische Mathemtik Nachtrag zu Se
Seite 462 und 463:
Annahme: Für alley 0 ∈ D existie
Seite 464 und 465:
Alternative Perspektive: Hamiltonsc
Seite 466 und 467:
➁ direkte Verifikation von Def. 4
Seite 468 und 469:
q 9 8 7 6 5 4 (p(t),q(t)) RK4 Metho
Seite 470 und 471:
11 10 9 t=0 t=0.5 t=1 t=2 t=3 t=5 N
Seite 472 und 473:
Verallgemeinerung auf höhere Dimen
Seite 474 und 475:
Beweis. L = −L T ⇒ unitär diag
Seite 476 und 477:
Istψ : U ↦→ R d eine Parametri
Seite 478 und 479:
✎ Notation: ∇ 2 H ˆ= (symmetri
Seite 480 und 481:
Einschrittverfahren fürẏ = f(y)
Seite 482 und 483:
Beweis (von Thm. 4.4.15) Numerische
Seite 484 und 485:
Bemerkung 4.4.19 (Einfache symplekt
Seite 486 und 487:
✬ ✩ Theorem 4.4.23 (Symplektisc
Seite 488 und 489:
(4.4.26) = explizite symplektische
Seite 490 und 491:
in Stufenform, vgl. Bem. 2.3.7: ⎧
Seite 492 und 493:
11 10 9 t=0 t=0.5 t=1 t=2 t=3 t=5 1
Seite 494 und 495:
Symplektische ESV: Symplektisches p
Seite 496 und 497:
1 0.5 Stoermer−Verlet, h = 0.2000
Seite 498 und 499:
➥ Hamiltonsche Differentialgleich
Seite 500 und 501:
100 Trajektorien der Atome, Verlet,
Seite 502 und 503:
2D konservatives Vielteilchensystem
Seite 504 und 505:
1.4 1.2 Stoermer−Verlet: 10 x 10
Seite 506 und 507:
5 Stoermer−Verlet Explicit TR Pro
Seite 508 und 509:
Konkrete Anwendung dieser Philosoph
Seite 510 und 511:
? Modifizierte Gleichung im allgeme
Seite 512 und 513:
Durchwegs “stillschweigende Annah
Seite 514 und 515:
➣ Beobachtung: Taylorentwicklung
Seite 516 und 517:
Explizites Euler-Verfahren (1.4.2):
Seite 518 und 519:
Explicit Euler h=0.050000, logistic
Seite 520 und 521:
Im Sinne der Rückwärtsanalyse (
Seite 522 und 523:
JEDOCH: γ = 1, L = 1 Numerische Ma
Seite 524 und 525:
Wir sind frei in der Wahl der Absch
Seite 526 und 527:
✬ ✩ Theorem 4.4.66 (Konsistenzf
Seite 528 und 529:
Ziel, vgl. (4.4.67): Abschätzung d
Seite 530 und 531:
Wir haben nun gesehen, wie man unte
Seite 532 und 533:
Idee: „∀α” in (4.4.73) ➣ n
Seite 534 und 535:
Wir sind „frei” in der Wahl von
Seite 536 und 537:
Dann weiter wie zuvor skizziert, si
Seite 538 und 539:
10 0 1/C 2 h 10 −2 40 35 Numerisc
Seite 540 und 541:
✎ Notationen: ỹ ˆ= Lösung des
Seite 542 und 543:
✬ ✩ Lemma 4.4.83 (Störungsabsc
Seite 544 und 545:
Also: |˜f h,l (y)−f(y)| ≤ Ch f
Seite 546 und 547:
➀ Explizites Euler-Verfahren für
Seite 548 und 549:
Konkret: mathematisches Pendel →
Seite 550 und 551:
✬ ✩ Theorem 4.4.88 (Symplektizi
Seite 552 und 553:
Erklärung der “Langzeitenergieer
Seite 554 und 555:
✬ ✩ Theorem 4.4.91 (Langzeitene
Seite 556 und 557:
4.5 Methoden für oszillatorische D
Seite 558 und 559:
ω ↑ ⇒ Oszillationen iny(t) ↑
Seite 560 und 561:
1.5 1 Gautschi−Verfahren: h = 0.1
Seite 562 und 563:
10 1 Gautschi−Verfahren: Konverge
Seite 564 und 565:
15 Gautschi−Verfahren: Konvergenz
Seite 566 und 567:
1 0.8 Numerische Mathemtik sinc(x)
Seite 568 und 569:
0.14 0.12 Filterndes Gautschi−Ver
Seite 570 und 571:
Numerische Mathemtik Verzeichnisse
Seite 572 und 573:
Newtonsche, 37 bimolekulare Reaktio
Seite 574 und 575:
Gauss-Radau-Quadratur, 369 Gaussqua
Seite 576 und 577:
lokale, 46 Logistische Differential
Seite 578 und 579:
Runge-Kutta 3/8-Regel, 230 Affin-Ko
Seite 580 und 581:
Zwangskraft, 420 Zweischrittverfahr
Seite 582 und 583:
Einfache reversible Einschrittverfa
Seite 584 und 585:
Schrittweitensteuerung für explizi
Seite 586 und 587:
Stetiger linearer Operator, 164 Sym
Seite 588 und 589:
Numerische Mathemtik Notationen C l
Seite 590 und 591:
• Beispiel 1.4.15 • Beispiel 1.
Seite 592 und 593:
• Beispiel 3.7.6 • Beispiel 3.8
Seite 594 und 595:
[8] P. DEUFLHARD AND F. BORNEMANN,
Seite 596:
[27] E. LORENZ, Deterministic non-p
Alle anzeigen

Beispiel - SAM - ETH ZÃ¼rich

Erfolgreiche ePaper selbst erstellen

Template löschen?

Als Template speichern?