Beispiel - SAM - ETH Zürich

1
0.5
Stoermer−Verlet, h = 0.200000
1
0.5
Explicit trapezoidal rule, h = 0.200000
Numerische
Mathemtik
0
0
−0.5
−0.5
−1
−1
q 2
q 2
−1.5
−1.5
−2
−2
−2.5
−2.5
−3
0 < t < 50
1000 < t < 1050
−3.5
−1.5 −1 −0.5 0 0.5 1 1.5
q 1
Fig. 170
Animation
✄
Störmer-Verlet ESV
−3
0 < t < 50
1000 < t < 1050
−3.5
−1.5 −1 −0.5 0 0.5 1 1.5
q 1
Fig. 171
Explizite Trapezregel
R. Hiptmair
rev 35327,
25. April
2011
Symplektischer Integrator: Positionen im “zulässigen Bereich” auch bei Langzeitintegration
Explizite Trapezregel: Trajektorien verlassen bei Langzeitintegration den “zulässigen Bereich” (Energiedrift
!)
4.4
p. 496