Discussion :

[Rorchar]

Bonjour,

Voici ce que j'aimerai coder :
- Un cell array de dimension n*n contient dans chacune de ses cellules un vecteurs X[x;y;z] avec x,y,z des réels positifs (des densité de population)
- Une fonction de transition doit s'appliquer sur le cell array afin de calculer les densités à l'instant t+1 à partir des densités de l'instant t. Dans cette fonction de transition la partie s'appliquant au calcul de X(3){t+1} nécessite de faire la somme d'une proportion de X(3) des huit cellules voisines (ce qui correspond dans le contexte à un phénomène de dispersion entre les cellules).

Voici la fonction de projection dans le cas simple de seulement deux cellules. En rouge les termes qui décrivent l'immigration

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
% Projection function
% X=(F1,B1,P1,F2,B2,P2)
% X_n+1=G2(X_n)
G2 = @(X) [v*phib(R(X(1),X(3))).*X(1) + r(rho(X(1),X(2))).*X(2);...
          mub(rho(X(1),X(2))).*X(2).*(1+gamma(rho(X(1),X(2)))).*(1-X(2)./B_max);...
          f.*X(3).*S(R(X(1),X(3)))*sigma*deltab(X(3),X(1))+delta(X(6),X(4))*sigma*S(R(X(4),X(6)))*X(6)*f;...
          v*phib(R(X(4),X(6))).*X(4) + r(rho(X(4),X(5))).*X(5);...
          mub(rho(X(4),X(5))).*X(5).*(1+gamma(rho(X(4),X(5)))).*(1-X(5)./B_max);...
          f.*X(6).*S(R(X(4),X(6)))*sigma*deltab(X(6),X(4))+delta(X(3),X(1))*sigma*S(R(X(1),X(3)))*X(3)*f];

Par rapport à ce code l'idée serait donc maintenant que l'immigration proviennent des huit cellules voisines.

Donc au final mes problèmes sont :
Comment modifier la fonction de transition
Comment appliquer la fonction sur l'ensemble du cell array en même temps

Merci

[le fab]

salut

pour la deuxième question, regarde la fonction cellfun

[Rorchar]

Bjr ,

Merci pour la réponse, fonction très utile ! Mais dans mon cas elle me pose un petit problème par rapport au dispersion qui risque d'être comptées plusieurs fois, enfin bref je me suis débrouillé différemment avec une boucle.
Mais cela m'amène à une nouvelle question plutôt graphique, comment visualiser mon automate graphiquement.

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
for t=1:tmax
    A = [] ;
    for i=1:n
        for j=1:n
            A = [A;M{i,j}(3)*1/8*delta(M{i,j}(1),M{i,j}(3))*sigma*f] ; 
        end
    end
    B = reshape(A,[n,n]) ; 
    C = conv2(B, F, 'same') ; 
 
    for i=1:n
        for j=1:n
        X = [M{i,j}(1) ; M{i,j}(2) ; M{i,j}(3)] ;
        M{i,j} = Gca(X) ;
        end
    end
end

L'idée étant la suivante :
Visualiser l'évolution des valeurs des cellules de M au fur et à mesure de la boucle avec un code couleur, par exemple :
- Vert si M{i,j}(1) < 0
- Rouge si M{i,j}(3) < 0
- Marron si M{i,j}(1) < 0 et M{i,j}(2) > 0
- Noir si M{i,j}(2) < 0

J'avais une piste avec pcolor() mais j'ai l'impression que ça ne prête pas trop à l'utilisation d'un cell array...