1 / 68

Le modèle de Lotka-Volterra

Le modèle de Lotka-Volterra. Version Volterra. Croissance exponentielle proies x en absence de prédateurs Décroissance exponentielle des prédateurs y en absence de proies

avak
Download Presentation

Le modèle de Lotka-Volterra

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. Le modèle de Lotka-Volterra

  2. Version Volterra • Croissance exponentielle proies x en absence de prédateurs • Décroissance exponentielle des prédateurs y en absence de proies • L’ « interaction » proies-prédateurs est proportionnelle au nombre de rencontres entre les deux espèces, donc à leurs densités respectives(les coeff. B et D sont différents car la disparition d’une proie n’entraîne pas de facto la naissance d’un prédateur)

  3. Modèle de Lotka-Volterra k1 Z + X 2X k2 Y + X 2Y Version Lotka k3 Y W k2 = B k3 = C Si [Z] >> [X] alors k1Z=cste=A

  4. Etats stationnaires 1ère solution : Premier état stationnaire : u0=0, v0=0 2ème solution : Deuxième état stationnaire : u0=1, v0=1

  5. Analyse des modes normaux Matrice jacobienne Eléments de la jacobienne Jacobienne à l’état stationnaire 1 (u0=0, v0=0) Jacobienne à l’état stationnaire 2 (u0=1, v0=1)

  6. Jacobienne à l’état stationnaire 1 (u0=0, v0=0) Analyse des modes normaux Trace de la jacobienne : Déterminant de la jacobienne : Equation caractéristique : Discriminant de l’équation caractéristique : La constante a étant toujours positive, D est négatif, ce qui signifie que les deux racines R1 et R2 sont de signes contraires. L'état stationnaire est donc instable et de type point de selle. Les trajectoires s'écartent de cet état. 5 6

  7. Jacobienne à l’état stationnaire 2 (u0=1, v0=1) Analyse des modes normaux Trace de la jacobienne : Déterminant de la jacobienne : Equation caractéristique : Discriminant de l’équation caractéristique : Le discriminant étant toujours négatif, l'équation caractéristique admet deux racines toujours imaginaires. La constante a étant toujours positive, D est positif. La trace est toujours nulle. L'état stationnaire est un centre (cas très particulier). Les trajectoires sont (dans le système linéarisé) des cercles centrés sur l'état stationnaire. 3 2

  8. Sensibilité du modèle de Lotka-Volterra aux conditions initiales

  9. Oscillations glycolytiques

  10. Pye, 1971 Extrait S. carlbergensis Milieu de culture : tréhalose

  11. Influence de la vitesse d’injection du substrat

  12. + S P vPFK (a) (g) a

  13. + S P vin vPFK vout

  14. Etat(s) stationnaire(s) du système soit avec dont l’unique racine positive est :

  15. Analyse des modes normaux

  16. Analyse des modes normaux : déterminant de la jacobienne L’état stationnaire (a*, g*) ne peut donc être en aucun cas du type point de selle.

  17. Analyse des modes normaux : trace de la jacobienne Condition d’instabilité : Tr* > 0 soit

  18. Analyse des modes normaux : discriminant La valeur stationnaire du discriminant peut être positive ou négative. L’état stationnaire peut donc être, a priori, stable ou instable, de type « nœud » ou « foyer », selon la valeur des paramètres.

  19. 1. 0 < µ < µ1(d) = 0.032144 tr* < 0, d* > 0 nœud stable 2. µ1(d) = 0.032144 < µ < µ1(tr) = 0.03598 tr* < 0, d* < 0 foyer stable 3. µ1(tr) = 0.03598 < µ < µ2(d) = 0.04098 tr* > 0, d* < 0 foyer instable 4. µ2(d) = 0.04098 < µ < µ3(d) = 0.3574 tr* > 0, d* > 0 nœud instable 5. µ3(d) = 0.3574 < µ < µ2(tr) = 0.4572 tr* > 0, d* < 0 foyer instable 6. µ2(tr) = 0.4572 < µ < µ4(d) = 0.5543 tr* < 0, d* < 0 foyer stable 7. µ4(d) = 0.5543 < µ tr* < 0, d* > 0 nœud stable µ = 0.36 3

  20. 1. 0 < µ < µ1(d) = 0.032144 tr* < 0, d* > 0 nœud stable 2. µ1(d) = 0.032144 < µ < µ1(tr) = 0.03598 tr* < 0, d* < 0 foyer stable 3. µ1(tr) = 0.03598 < µ < µ2(d) = 0.04098 tr* > 0, d* < 0 foyer instable 4. µ2(d) = 0.04098 < µ < µ3(d) = 0.3574 tr* > 0, d* > 0 nœud instable 5. µ3(d) = 0.3574 < µ < µ2(tr) = 0.4572 tr* > 0, d* < 0 foyer instable 6. µ2(tr) = 0.4572 < µ < µ4(d) = 0.5543 tr* < 0, d* < 0 foyer stable 7. µ4(d) = 0.5543 < µ tr* < 0, d* > 0 nœud stable µ = 0.45 3

  21. 1. 0 < µ < µ1(d) = 0.032144 tr* < 0, d* > 0 nœud stable 2. µ1(d) = 0.032144 < µ < µ1(tr) = 0.03598 tr* < 0, d* < 0 foyer stable 3. µ1(tr) = 0.03598 < µ < µ2(d) = 0.04098 tr* > 0, d* < 0 foyer instable 4. µ2(d) = 0.04098 < µ < µ3(d) = 0.3574 tr* > 0, d* > 0 nœud instable 5. µ3(d) = 0.3574 < µ < µ2(tr) = 0.4572 tr* > 0, d* < 0 foyer instable 6. µ2(tr) = 0.4572 < µ < µ4(d) = 0.5543 tr* < 0, d* < 0 foyer stable 7. µ4(d) = 0.5543 < µ tr* < 0, d* > 0 nœud stable µ = 0.55 2

  22. Méthode des isoclines nulles

  23. Diffusion Première loi de Fick : le flux matériel J (cellules, concentration, nombre d’individus appartenant à une population, etc.) est proportionnel au gradient de « concentration » de ce matériel, ce qui, dans un espace à une dimension, s’exprime mathématiquement par : D ( coefficient de diffusion) a pour dimension [longueur]2[temps]-1 (cm2 sec-1 par exemple). Deuxième loi de Fick : Equation générale d’un processus de réaction-diffusion : (diffusion à une dimension)

  24. Diffusion Réaction – diffusion à une dimension : Formalisme plus général : avec, dans le cas d’une diffusion monodirectionnelle : Cas général d’une diffusion à 3 dimensions :

  25. Scilab cd ..\g

  26. Qu‘est-ce que le phénomène d’excitabilité ? • Dans quelles conditions paramétriques peut-il être observé ? • Pourquoi n’observe-t-on jamais qu’une seule vague (et non, par exemple, un train d’oscillations) ? • Pourquoi lorsqu’ils se rencontrent, les deux fronts de vague s’annihilent-ils (ce qui a pour conséquence de faire que chaque point du cortex est soumis une et une seule fois au passage de la vague) ? • Quel est le mécanisme expliquant la transmission d’une unité à l’autre de la fluctuation initiale (phénomène de relais) ? La transmission est-elle soumise à des contraintes paramétriques particulières ?

  27. L’excitabilité est une propriété d’un état stationnaire stable d’un système qui posséde au moins une isocline nulle ayant une branche de pente négative (forme en S). Dans la plupart des zones de pente positive, l’état stationnaire est très stable : le système ne peut amplifier un pulse de concentration. Cependant, à mesure que l’état stationnaire s’approche de l’extremum de la branche stable passant par l’origine (0, 0), le système devient excitable, c’est-à-dire capable d’amplifier un pulse de concentration (ici g). La zone d’excitabilité est en général très étroite.

  28. Condition à l’excitation effective d’un système excitable : Pour que la fluctuation soit amplifiée, il faut que son amplitude dépasse un seuil. Le seuil diminue à mesure que l’état stationnaire soumis à la fluctuation approche le domaine instable. L’excitabilité se définit ainsi comme la capacité d’un système à amplifier, au-delà d’un seuil, un pulse d’une de ses variables, à partir d’un état stationnaire stable.

  29. Trajectoire dans le plan de phase, suite à une excitation supraliminaire

  30. A t= 0, Dg=0.4

  31. Perturbations et période réfractaire Dg = 2.0 Dg = 0.4 Dg = 0.2 g seuil

  32. Période réfractaire annihilation de la vague

  33. Le phénomène de relais • Transmission de l’excitation d’une unité à l’autre, via la diffusion du calcium cytosolique Programme VagueCaRelais (nb_cycles = 100) Existence d’une limite supérieure et d’une limite inférieure à la diffusion entre lesquelles le signal d’excitation amplifié est transmis entre unités corticales. Détermination des valeurs limites de d pour lesquelles le relais se produit

  34. De l'excitabilité aux oscillations µ = 0.2 µ = 0.0333

  35. Relais d'oscillations Tous les compartiments basculent dans un état oscillant : synchronisation des oscillations, en phase dans tous les compartiments Un seul compartiment bascule dans un état oscillant : ondes concentriques en provenance d’un émetteur Une rangée de compartiments basculent dans un état oscillant : vagues successives Dans les 2 derniers cas, comment sont relayées les oscillations par les unités restées dans l’état excitable ? Programme VagueCaRelais

More Related