Forum Sceptique

// Un jour, l’environnement de la population des pikmins change : alors que jusqu'à présent, leur nourriture, le pikmiam, se trouvait au sol, un accident magnetico-telluro-quantique le fait maintenant léviter. Comment les pikmins vont-ils s'adapter à leur nouvel environnement ?

// effacer tous les graphs des simulations précédentes
xdel(winsid());
// effacer toutes les variables des simulations précédentes
clear;

//------------------------------------------------------------------------------
//paramètres qu'on peut bidouiller

// nombre de "pas de temps" (jours par ex) par génération
T=10;
// nombre total de générations
G=10;
// nombre d'individus de la première génération
N=1000;
// hauteur minimale de la nourriture ie difficulté à obtenir de la nourriture (entre 0 et 1)
difF=0;
// difficulté à faire un bébé (entre 0 (lapin) et 1 (panda) :-) )
difB=0.3;
// stock initial de nourriture de chaque individu
stockI=0;

//------------------------------------------------------------------------------

// la moitié de la population (les h premiers sont des filles, les autres des garçons ou vice-versa, l'individu n est en couple avec l'individu n+h)
h=floor(N/2);

//taille des individus de la première génération définie aléatoirement entre 0.1 et 1
taille=0.9*rand(1,N)+0.1;

// définition du vecteur des abscisses pour l'histogramme
[x]=linspace(0.1,1,90);

// calcul et enregistrement des tailles min, max et moy du nombre d'individus de la population
m(1,1)=min(taille);
M(1,1)=max(taille);
mm(1,1)=mean(taille);
Nn(1,1)=N;

// histogramme des tailles
figure;
histplot(x,taille);
xtitle('Histogramme des tailles', 'Tailles', 'Nombre d''individus');

// à chaque génération
for g=1:G

// chaque individu reçoit un stock de nourriture égal à stockI
stock=stockI*ones(1,N);

j=1;

// à chaque pas de temps
for t=1:T

// chaque individu trouve de la nourriture à une certaine hauteur déterminée aléatoirement entre 0 et 1
testfood=rand(1,N);

// pour chaque individu
for n=1:N

// si l'individu est aussi grand ou plus grand que la hauteur où se situe la nourriture qu'il a trouvé, il la prend
if taille(1,n)>=testfood(1,n)+difF

// dans ce cas son stock augmente de 1
stock(1,n)=stock(1,n)+1;

end

// si le stock de l'individu atteint 0, il meurt de faim :'-(
if stock(1,n)<=0;

// dans ce cas, il est "effacé" de la liste des tailles
taille(1,n)=0;

end

end

// chaque couple fait un lancé de dé pour avoir un bébé (si, si, c'est comme ça qu'on fait les bébés...:-) )
testbaby=rand(1,h);

// pour chaque couple
for n=1:h

// si le résultat du lancer de dé est supérieur à la valeur de la difficulté à avoir un bébé, et qu'aucun des deux parents n'est déjà mort, la cigogne leur amène un bébé
if taille(1,n)>0 & taille(1,n+h)>0 & testbaby(1,n)>difB

// la taille du bébé est une valeur aléatoire entre les tailles des deux parents
taillebaby(1,j)=rand()*abs(taille(1,n)-taille(1,n+h))+min(taille(1,n),taille(1,n+h));

// il y a un bébé de plus !
j=j+1;

end

end //(boucle sur les individus)

// chaque individu mange et perd 1 de nourriture (à chaque pas de temps)
stock=stock-1;

end //(boucle sur le temps)

// le nombre d'individus de la nouvelle génération est égal aux nombres de bébés
N=size(taillebaby,2);

// la moitié (h)sont des filles, l'autre des garçons
h=floor(N/2);

// la liste des tailles de la nouvelle génération est la liste des tailles des bébés
taille=taillebaby;

// histogramme des tailles
figure;
histplot(x,taille);
xtitle('Histogramme des tailles', 'Tailles', 'Nombre d''individus');

// calcul et enregistrement des tailles min, max et moy du nombre d'individus de la population
m(g+1,1)=min(taille);
M(g+1,1)=max(taille);
mm(g+1,1)=mean(taille);
Nn(g+1,1)=N;

end // (boucle sur les générations)

// définition du vecteur des abscisses pour les graphs
xx=[1:1:G+1];

// graph tailles min, max, moy de la population au cours des générations
figure;
plot2d(xx,[m M mm],leg="min@max@moy");
xtitle('Evolution des tailles', 'Génération', 'Taille');

//graph nombre d'individus de chaque génération
figure;
plot2d(xx,Nn);
xtitle('Evolution de la population', 'Génération', 'Nombre d''individus');

// Un jour, l'environnement de la population des pikmins change : alors que jusqu'à présent, leur nourriture, le pikmiam, se trouvait au sol, un accident magnetico-telluro-quantique le fait maintenant léviter. Comment les pikmins vont-ils s'adapter à leur nouvel environnement ?

// effacer tous les graphs des simulations précédentes
xdel(winsid());
// effacer toutes les variables des simulations précédentes
clear;

//------------------------------------------------------------------------------
//paramètres qu'on peut bidouiller

// nombre de "pas de temps" (jours par ex) par génération
T=10;
// nombre total de générations
G=10;
// nombre d'individus de la première génération
N=1000;
// hauteur minimale de la nourriture ie difficulté à obtenir de la nourriture (entre 0 et 1)
difF=0;
// difficulté à faire un bébé (entre 0 (lapin) et 1 (panda) :-) )
difB=0.3;
// stock initial de nourriture de chaque individu
stockI=0;
// nombre max de pikmins, au delà, le programme quitte
limite=50000;

//------------------------------------------------------------------------------

// la moitié de la population (les h premiers sont des filles, les autres des garçons ou vice-versa, l'individu n est en couple avec l'individu n+h)
h=floor(N/2);

//taille des individus de la première génération définie aléatoirement entre 0.1 et 1
taille=0.9*rand(1,N)+0.1;

// définition du vecteur des abscisses pour l'histogramme
[x]=linspace(0.1,1,90);

// calcul et enregistrement des tailles min, max et moy du nombre d'individus de la population
m(1,1)=min(taille);
M(1,1)=max(taille);
mm(1,1)=mean(taille);
Nn(1,1)=N;

// histogramme des tailles
figure;
histplot(x,taille);
xtitle('Histogramme des tailles', 'Tailles', 'Nombre d''individus');

// à chaque génération
for g=1:G

// chaque individu reçoit un stock de nourriture égal à stockI
stock=stockI*ones(1,N);

j=1;

// à chaque pas de temps
for t=1:T

// chaque individu trouve de la nourriture à une certaine hauteur déterminée aléatoirement entre 0 et 1
testfood=rand(1,N);

// pour chaque individu
for n=1:N

// si l'individu est aussi grand ou plus grand que la hauteur où se situe la nourriture qu'il a trouvé, il la prend
if taille(1,n)>=testfood(1,n)+difF

// dans ce cas son stock augmente de 1
stock(1,n)=stock(1,n)+1;

end

// si le stock de l'individu atteint 0, il meurt de faim :'-(
if stock(1,n)<=0;

// dans ce cas, il est "effacé" de la liste des tailles
taille(1,n)=0;

end

end

// chaque couple fait un lancé de dé pour avoir un bébé (si, si, c'est comme ça qu'on fait les bébés...:-) )
testbaby=rand(1,h);

// pour chaque couple
for n=1:h

// si le résultat du lancer de dé est supérieur à la valeur de la difficulté à avoir un bébé, et qu'aucun des deux parents n'est déjà mort, la cigogne leur amène un bébé
if taille(1,n)>0 & taille(1,n+h)>0 & testbaby(1,n)>difB

// la taille du bébé est une valeur aléatoire entre les tailles des deux parents
taillebaby(1,j)=rand()*abs(taille(1,n)-taille(1,n+h))+min(taille(1,n),taille(1,n+h));

// il y a un bébé de plus !
j=j+1;

if j>limite

disp("Désolé, une météorite a anéanti votre population de pikmins (nombre max de pikmins atteint)");

resume;

end


end

end //(boucle sur les individus)

// chaque individu mange et perd 1 de nourriture (à chaque pas de temps)
stock=stock-1;

end //(boucle sur le temps)

// le nombre d'individus de la nouvelle génération est égal aux nombres de bébés
N=size(taillebaby,2);

// la moitié (h)sont des filles, l'autre des garçons
h=floor(N/2);

// la liste des tailles de la nouvelle génération est la liste des tailles des bébés
taille=taillebaby;

// histogramme des tailles
figure;
histplot(x,taille);
xtitle('Histogramme des tailles', 'Tailles', 'Nombre d''individus');

// calcul et enregistrement des tailles min, max et moy du nombre d'individus de la population
m(g+1,1)=min(taille);
M(g+1,1)=max(taille);
mm(g+1,1)=mean(taille);
Nn(g+1,1)=N;

end // (boucle sur les générations)

// définition du vecteur des abscisses pour les graphs
xx=[1:1:G+1];

// graph tailles min, max, moy de la population au cours des générations
figure;
plot2d(xx,[m M mm],leg="min@max@moy");
xtitle('Evolution des tailles', 'Génération', 'Taille');

//graph nombre d'individus de chaque génération
figure;
plot2d(xx,Nn);
xtitle('Evolution de la population', 'Génération', 'Nombre d''individus');