TP LVfinal
TP LVfinal
TP LVfinal
Laurent Dietrich
TP – Équations de Lotka-Volterra (2)
Il est recommandé de créer un nouveau répertoire dans lequel seront placés les fi-
chiers relatifs à ce TP afin de ne pas mélanger, ou pire, écraser, d’anciens fichiers. Il
est chaudement recommandé de taper help quelquechose dans l’invite de commande de
Matlab, lorsque vous voulez des renseignements sur la commande quelquechose et/ou de
m’interpeller. Il est aussi chaudement recommandé de ne pas aller recopier du code d’un
ancien TP afin d’éviter de recopier au pire des erreurs, au mieux quelque chose dont
le code n’aura pas la même structure qu’ici, ce qui vous fera perdre du temps, et de
l’entraînement.
1
Champ de vecteurs associé à a=2, b=c=d=1
4.5
3.5
2.5
2
y
1.5
0.5
−0.5
−0.5 0 0.5 1 1.5 2 2.5 3 3.5
x
3 Données
En écologie, la collecte de donnée est un problème compliqué. Néanmoins, et avec
toutes les précautions que l’on peut prendre, on prendra en considération le graphe
suivant ainsi que ce texte explicatif de Douglas R. Hundley
(http ://people.whitman.edu/ hundledr/)
NOTE about the data : I have not been able to verify this data,but this is the
data (or rather the graph) that is always cited. This particular set of data
came from scanning in the graph from Odum’s ”Fundamentals of Ecology”, p.
191 which is often cited. Odum says that his graph is taken from MacLulich’s
”Fluctuations in the numbers of varying hare”, 1937, which is not widely
available. Some authors caution that this data is actually a composition of
several time series, and should probably not be analyzed as a whole, and
that some of the lynx data was actually missing. It is said that the data was
collected from Hudson’s Bay historical records, and does not reflect animal
2
populations, but rather the number of pelts turned in for trading (a large
number of which came from Native Americans-mentioned because there were
some medical outbreaks during these years which could account for skewed
data). The data is presented here with these cautions.
Figure 2 – Relevés
On pourra mettre en rapport les résultats obtenus et ces relevés, ainsi qu’avec les
propriétés globales des solutions qu’on a vues en cours.
H(x, y) = x2 + y 2
Le calcul de H le long des trajectoires approchées peut alors se faire à la main et
illustrer précisément les observations précédentes.
3
3. Dans ce cadre, la méthode de Crank-Nicholson est dite symplectique. Quelle pro-
priété particulière a-t-elle par rapport aux autres ?
5 Ordres de convergence
Comme on ne connait pas d’expression analytique des solutions des équations de
Lotka-Volterra, on ne peut pas mettre en valeur l’ordre de convergence des méthodes
numériques précédemment employées. On peut cependant le faire sur des problèmes plus
simples. Proposez donc un protocole simple pour rendre compte visuellement des divers
ordres de convergences (1, 2, 4) des méthodes employées.