We analyse the preservation of physical properties of numerical approximations tosolutions of the Lotka-Volterra system: its positivity and the conservation of theHamiltonian. We focus on two numerical methods : the symplectic Euler method andan explicit variant of it. We first state under which conditions they are symplectic andwe prove they are both Poisson integrators for the Lotka-Volterra system. Then, westudy under which conditions they stay positive. For the symplectic Euler method,we derive a simple condition under which the numerical approximation always stayspositive. For the explicit variant, there is no such simple condition. Using propertiesof Poisson integrators and backward error analysis, we prove that for initial conditionsin a given set in the positive quadrant, there exists a bound on the step size, such thatnumerical approximations with step sizes smaller than the bound stay positive overexponentially long time intervals. We also show how this bound can be estimated.We illustrate all our results by numerical experiments.
| Identifer | oai:union.ndltd.org:LACETR/oai:collectionscanada.gc.ca:QMM.19583 |
| Date | January 2003 |
| Creators | Beck, Mélanie |
| Publisher | McGill University |
| Source Sets | Library and Archives Canada ETDs Repository / Centre d'archives des thèses électroniques de Bibliothèque et Archives Canada |
| Language | English |
| Detected Language | English |
| Type | Electronic Thesis or Dissertation |
| Format | application/pdf |
| Coverage | Master of Science (Department of Mathematics and Statistics) |
| Rights | All items in eScholarship@McGill are protected by copyright with all rights reserved unless otherwise indicated. |
| Relation | alephsysno: 002022446, Theses scanned by McGill Library. |
Page generated in 0.0016 seconds