Bonjour,
Je suis étudiant un M2 de biologie, avec une orientation en modélisation. Vous comprendrez donc le préfixe "débutant" dans le sens que je ne suis absolument pas passé par une formation du style "maths-info".
Je travaille actuellement, dans le cadre de mon stage de M2, sur un projet de modélisation de dynamique de population.
Un modèle (discret) reliant la dynamique d'un insecte et de sa plante hôte a déjà été réalisé pour une situation avec un unique patch dans lequel les deux populations évoluaient. Modèle qui mène à l'extinction des populations, mais cela est cohérent avec les résultats imaginés.
La problématique du projet est de savoir la dynamique de ces populations dans un modèle spatialisé constituer de multiples patch et donc d'une migration de l'insecte entre ces différents patch (est-ce que la spatialisation permet le maintien des populations ? la disparition d'une seule ?...)
Ainsi, un modèle à deux patchs, avec une dispersion entre les deux a été mis au point.
Pour faire rapide : F représente la biomasse de feuille, B la biomasse de bois (qui est source de feuille), et P la biomasse de l'insecte qui se nourrit des feuilles. La survie du bois, sa croissance, ainsi que la croissance des feuilles sont facteurs de la défoliation par l'insecte. Les indices (1) et (2) signifiant patch1 et patch2.
Je pourrai rentrer plus dans les détails ultérieurement si nécessaire mais je ne suis pas sûr que cela vous sera utile pour mon questionnement. Je rajouterai quand même que le paramètre est le paramètre qui définit la dispersion de l'insecte et qui est fonction de la biomasse de feuille (encore disponible comme ressource) dans le patch. Plus un patch contient beaucoup d'insecte par rapport à sa biomasse de feuille plus une proportion importante d'insecte va se disperser vers un autre patch.
J'en viens donc à ma question.
J'aimerai à présent faire tourner ce modèle non plus seulement avec deux patchs mais sur une grille dont chaque cellule serait un patch, certain vide, certain avec de la biomasse végétale. La population d'insecte serait initialement sur un seul patch et sa dispersion sur la grille serait ensuite visualisée un peu comme le ferait un automate cellulaire.
J'ai effectué quelques recherches sur la programmation d'automate cellulaire sur MATLAB mais rien de bien concluant par rapport à mon problème. En fait (comme le veulent les règles du forum de toute manière) je ne cherche pas de solution toute faite mais une piste de travail pour mettre en place ce modèle spatialisé sur une grille.
Merci pour votre réponse
Partager