Équations de Lotka-Volterra

Équations de Lotka-Volterra

En mathématiques, les équations de Lotka-Volterra, que l'on désigne aussi sous le terme de "modèle proie-prédateur", sont un couple d'équations différentielles non-linéaires du premier ordre, et sont couramment utilisées pour décrire la dynamique de systèmes biologiques dans lesquels un prédateur et sa proie interagissent. Elles ont été proposées indépendamment par Alfred James Lotka en 1925 et Vito Volterra en 1926.

Ce système d'équations est classiquement utilisé comme modèle pour la dynamique du lynx et du lièvre des neiges, pour laquelle de nombreuses données de terrain ont été collectées sur les populations des deux espèces par la Compagnie de la baie d'Hudson au XIXe siècle. Il a aussi été employé par John Allan Hobson pour décrire les relations entre les neurones cholinergiques responsables du sommeil paradoxal et les neurones aminergiques liées à l'état de veille.[réf. nécessaire]

Sommaire

Les équations

Elles s'écrivent fréquemment :

\left\{\begin{array}{ccc}\dfrac{\mathrm{d}x(t)}{\mathrm{d}t} &=& x(t)\ \Big(\alpha - \beta y(t)\Big) \\
\dfrac{\mathrm{d}y(t)}{\mathrm{d}t} &=& -y(t) \ \Big(\delta - \gamma x(t)\Big) 
\end{array}\right.


  • x(t) est l'effectif des proies ;
  • y(t) est l'effectif des prédateurs ;
  • t est le temps ;
  • dx(t) / dt et dy(t) / dt représentent les taux de croissance des populations au cours du temps.

Les paramètres suivants caractérisent les interactions entre les deux espèces :

  • α, taux de reproduction des proies en l’absence de prédateurs ;
  • β, taux de mortalité des proies due aux prédateurs ;
  • δ, taux de mortalité des prédateurs en l’absence de proies.;
  • γ, taux de reproduction des prédateurs en fonction des proies mangées ;

Signification physique des équations

Une fois développées, les équations prennent une forme utile pour une interprétation physique.

Proies

L'équation de la proie devient :

\frac{\mathrm{d}x(t)}{\mathrm{d}t} =\alpha x(t) - \beta x(t) y(t)

Les proies sont supposées avoir une source illimitée de nourriture et se reproduire exponentiellement si elles ne sont soumises à aucune prédation ; cette croissance exponentielle est représentée dans l'équation ci-dessus par le terme αx(t). Le taux de prédation sur les proies est supposé proportionnel à la fréquence de rencontre entre les prédateurs et les proies ; il est représenté ci-dessus par βx(t)y(t). Si l'un des termes x(t) ou y(t) est nul, alors il ne peut y avoir aucune prédation.

Avec ces deux termes, l'équation peut alors être interprétée comme : la variation du nombre de proies est donnée par sa propre croissance moins le taux de prédation qui leur est appliqué.

Prédateurs

L'équation du prédateur devient :

\frac{\mathrm{d}y(t)}{\mathrm{d}t} = \gamma x(t) y(t) - \delta y(t)

Dans cette équation, γx(t)y(t) représente la croissance de la population prédatrice. (Notons la similarité avec le taux de prédation ; cependant, une constante différente est utilisée car la vitesse à laquelle la population des prédateurs augmente n'est pas nécessairement égale à celle à laquelle il consomme la proie). De plus, δy(t) représente la mort naturelle des prédateurs ; c'est une décroissance exponentielle. L'équation représente donc la variation de la population de prédateurs en tant que croissance de cette population, diminuée du nombre de morts naturelles.

Solutions de l'équation

On prouve que pour une condition initiale au temps t0 qui vérifie x(t0)>0 et y(t0)>0, l'unique solution maximale est définie pour tout réel t, et vérifie

 \forall t \in \mathbb{R}, \ x(t)>0, \ y(t)>0


La fonction F : (x,y) \mapsto \beta y + \gamma x - \alpha \ln(y) - \delta \ln(x) est alors une intégrale première du mouvement :  t \mapsto F(x(t),y(t)) est constante.

Les solutions maximales sont alors périodiques, et leur trajectoire est fermée bornée.

Évolution des populations au cours du temps.jpg

Les équations admettent des solutions périodiques qui n'ont pas d'expression simple à l'aide des fonctions trigonométriques habituelles. Néanmoins, une solution approximative linéarisée offre un mouvement harmonique simple, avec la population des prédateurs en retard de 90° (un quart de période) sur celle des proies.

Dynamique du système

Dans le modèle utilisé, les prédateurs prospèrent lorsque les proies sont nombreuses, mais finissent par épuiser leurs ressources et déclinent. Lorsque la population de prédateur a suffisamment diminué, les proies profitant du répit se reproduisent et leur population augmente de nouveau. Cette dynamique se poursuit en un cycle de croissance et déclin.

Equilibres de la population

Un état d'équilibre de la population est observé quand aucune des deux populations en présence n'évolue, c'est-à-dire quand les dérivées correspondantes sont nulles, ce qui se traduit par le système d'équations :

\left\{\begin{array}{lcl} x(t)(\alpha - \beta y(t)) &=& 0 \\-y(t)(\delta - \gamma  x(t)) &=& 0\end{array}\right.

qui a pour solutions :

\left\{y(t) = 0  ,  x(t) = 0\right\}\quad \mathrm{et}\quad \left\{y(t) = \frac{\alpha}{\beta}, x(t) = \frac{\delta}{\gamma}\right\}

La première solution correspond à une extinction définitive des deux espèces, la seconde à des valeurs pour les deux populations, dépendant des quatre paramètres α, β, δ et γ, qui restent stables indéfiniment.


Stabilité des points fixes

La stabilité des points fixes peut être déterminée par une linéarisation du système aux dérivées partielles. La matrice jacobienne du système est

J(x(t),y(t)) = \begin{bmatrix} 
\alpha - \beta y(t) & -\beta x(t) \\
\gamma y(t) & \gamma x(t) - \delta \\
\end{bmatrix}

Premier point fixe

Au premier point fixe (0,0), cette matrice prend la valeur :

J(0,0) = \begin{bmatrix}
\alpha & 0 \\
0 & -\delta \\
\end{bmatrix},

qui a pour valeurs propres :

\lambda_1 = \alpha,\quad \lambda_2 = -\delta.

Ces valeurs propres sont toujours de signes opposés, ce qui montre que ce point fixe est un point selle. Ce n'est donc pas un point fixe stable, ce qui montre en particulier que, suivant ce modèle, l'extinction des deux espèces en jeu est difficile à obtenir.

Second point fixe

En évaluant la matrice jacobienne en le second point fixe, la valeur suivante est obtenue :

J\left(\frac{\delta}{\gamma},\frac{\alpha}{\beta}\right) = \begin{bmatrix}
0 & -\frac{\beta \delta}{\gamma} \\
\frac{\alpha \gamma}{\beta} & 0 \\
\end{bmatrix}

et elle a pour valeurs propres :

\lambda_1 = i \sqrt{\alpha \delta},\quad \lambda_2 = -i \sqrt{\alpha \delta}

Ce point fixe est donc un foyer et plus particulièrement un centre, ce qui signifie que les populations de proies et prédateurs oscillent autour de leurs valeurs en ce point fixe.

Voir également

Bibliographie

  • E. R. Leigh (1968) The ecological role of Volterra's equations, in Some Mathematical Problems in Biology - a modern discussion using Hudson's Bay Company data on lynx and hares in Canada from 1847 to 1903.
  • Understanding Nonlinear Dynamics. Daniel Kaplan et Leon Glass.
  • V. Volterra. Variations and fluctuations of the number of individuals in animal species living together. In Animal Ecology. McGraw-Hill, 1931. Traduit de l'édition de 1928 par R. N. Chapman.

Liens externes


Wikimedia Foundation. 2010.

Contenu soumis à la licence CC-BY-SA. Source : Article Équations de Lotka-Volterra de Wikipédia en français (auteurs)

Игры ⚽ Нужна курсовая?

Regardez d'autres dictionnaires:

  • Equations de Lotka-Volterra — Équations de Lotka Volterra En mathématiques, les équations de Lotka Volterra, que l on désigne aussi sous le terme de modèle proie prédateur , sont un couple d équations différentielles non linéaires du premier ordre, et sont couramment… …   Wikipédia en Français

  • Équations de lotka-volterra — En mathématiques, les équations de Lotka Volterra, que l on désigne aussi sous le terme de modèle proie prédateur , sont un couple d équations différentielles non linéaires du premier ordre, et sont couramment utilisées pour décrire la dynamique… …   Wikipédia en Français

  • Système de Lotka-Volterra — Équations de Lotka Volterra En mathématiques, les équations de Lotka Volterra, que l on désigne aussi sous le terme de modèle proie prédateur , sont un couple d équations différentielles non linéaires du premier ordre, et sont couramment… …   Wikipédia en Français

  • Lotka–Volterra equation — The Lotka–Volterra equations, also known as the predator prey equations , are a pair of first order, non linear, differential equations frequently used to describe the dynamics of biological systems in which two species interact, one a predator… …   Wikipedia

  • Competitive Lotka–Volterra equations — The competitive Lotka–Volterra equations are a simple model of the population dynamics of species competing for some common resource. They can be further generalised to include trophic interactions. Contents 1 Overview 1.1 Two species 1.2 N… …   Wikipedia

  • Équations différentielles — Équation différentielle En mathématiques, une équation différentielle est une relation entre une ou plusieurs fonctions inconnues et leurs dérivées. L ordre d une équation différentielle correspond au degré maximal de différenciation auquel une… …   Wikipédia en Français

  • Équations différentielles ordinaires — Équation différentielle En mathématiques, une équation différentielle est une relation entre une ou plusieurs fonctions inconnues et leurs dérivées. L ordre d une équation différentielle correspond au degré maximal de différenciation auquel une… …   Wikipédia en Français

  • Volterra (homonymie) — Cette page d’homonymie répertorie les différents sujets et articles partageant un même nom. Volterra est une commune d Italie. Elle a donné son nom à : Daniele da Volterra, peintre italien de la Renaissance, au mathématicien Vito Volterra,… …   Wikipédia en Français

  • Volterra (disambiguation) — Volterra may refer to the following:* Volterra a town in Italy * Daniele da Volterra an Italian painter * Francesco da Volterra an Italian painter * Vito Volterra an Italian mathematician * Volterra Semiconductor an American semiconductor… …   Wikipedia

  • Alfred J. Lotka — Alfred James Lotka Alfred James Lotka (1880 1949) est un mathématicien et statisticien américain, théoricien de la dynamique des populations. Biographie Né le 2 mars 1880 à Lemberg en Autriche Hongrie (aujourd hui Lvov en Ukraine) de parents… …   Wikipédia en Français

Share the article and excerpts

Direct link
Do a right-click on the link above
and select “Copy Link”