function growth
%×÷Õߣºdbb627
%2011.10.27
%======É趨³õÖµ=============
clc;clear;
rand('state',0);
set(gcf,'DoubleBuffer','on');%Ïû³ýÕñ¶¯
Cd=0.01;%Éú³¤ºË·Ö²¼¸ÅÂÊÊý
Vr=0.3;%Éú³¤ÏàµÄÌå»ý·ÖÊý
global kn
kn=500;
p=[1/8 1/8 1/8 1/8 1/4 1/4 1/4 1/4];%·½ÏòÔö³¤¸ÅÂÊÊý
A=rand(kn,kn);%³õʼÍø¸ñkn¡Ákn
A(find(A<=Cd))=0;%ºÚÉ«£¬Éú³¤ºË
A(find(A>Cd))=1;%°×É«
imshow(A,'InitialMagnification','fit')%³õʼͼÏñ
k=0;%µü´ú²½Öè
N=kn*kn;%¸ñµãÊý
[x,y]=find(A==0);%Éú³¤ºË½Å±ê
Ixy=sub2ind([kn,kn],x,y); % ת»¯½Å±êΪË÷Òý
gn=length(Ixy);%Éú³¤ºË¼ÆÊý
Vs=gn/N;       %Ìå»ý·ÖÊý
ti=title(['time = 0 Ìå»ý·ÖÊý=',num2str(Vs)],'Fontsize',14,'Fontname','Times new roman'); % ÔÚͼÌ⴦ʵʱÏÔʾʱ¼ä
%===========================
while Vs     xL=x-1;xR=x+1;
    xL(xL==0)=1;% ¶Ô½üÁÚ½øÐб߽çÌõ¼þ´¦Àí
    xR(xR==kn+1)=kn;% ¶Ô½üÁÚ½øÐб߽çÌõ¼þ´¦Àí
    yD=y-1;yU=y+1;
    yD(yD==0)=1;% ¶Ô½üÁÚ½øÐб߽çÌõ¼þ´¦Àí
    yU(yU==kn+1)=kn;% ¶Ô½üÁÚ½øÐб߽çÌõ¼þ´¦Àí
   
    InE=sub2ind([kn,kn],xR,y); % ÓÒ·½×ª»¯½Å±êΪË÷Òý
    InS=sub2ind([kn,kn],x,yD); % Ï·½×ª»¯½Å±êΪË÷Òý
    InW=sub2ind([kn,kn],xL,y); % ×ó·½×ª»¯½Å±êΪË÷Òý
    InN=sub2ind([kn,kn],x,yU); % ÉÏ·½×ª»¯½Å±êΪË÷Òý
   
    InEN=sub2ind([kn,kn],xR,yU); % ÓÒÉÏ·½×ª»¯½Å±êΪË÷Òý
    InES=sub2ind([kn,kn],xR,yD); % ÓÒÏ·½×ª»¯½Å±êΪË÷Òý
    InWS=sub2ind([kn,kn],xL,yD); % ×óÏ·½×ª»¯½Å±êΪË÷Òý
    InWN=sub2ind([kn,kn],xL,yU); % ×óÉÏ·½×ª»¯½Å±êΪË÷Òý
     %===Óú¯ÊýµÄ·½·¨Ð´°Ë¸ö·½ÏòµÄÉú³¤µü´ú=====
          [A,Ixy]=growthmod(InE,Ixy,p(5),A);  
          [A,Ixy]=growthmod(InS,Ixy,p(6),A);
          [A,Ixy]=growthmod(InW,Ixy,p(7),A);
          [A,Ixy]=growthmod(InN,Ixy,p(8),A);
          [A,Ixy]=growthmod(InEN,Ixy,p(1),A);
          [A,Ixy]=growthmod(InES,Ixy,p(2),A);
          [A,Ixy]=growthmod(InWS,Ixy,p(3),A);
          [A,Ixy]=growthmod(InWN,Ixy,p(4),A);
  %==========×öͼ========   
    [x,y]=find(A==0);
    Ixy=sub2ind([kn,kn],x,y);
    gn=length(Ixy);
    Vs=gn/N %Ìå»ý·ÖÊý
    k=k+1
    pause(0.1)
    imshow(A,'InitialMagnification','fit')%¸üÐÂͼÏñ
    title(['time = ',num2str(k),'Ìå»ý·ÖÊý=',num2str(Vs)],'Fontsize',14,'Fontname','Times new roman'); % ¸üÐÂʱ¼ä²ÎÊý
   
end
   
function [A,Ixy]=growthmod(Ind1,Indxy,p,A)
%Ind1 µ±Ç°·½ÏòË÷Òý
%,Indxy µ±Ç°Éú³¤ºËË÷ÒýµÄ
%p µ±Ç°·½ÏòµÄ¸ÅÂÊÖµ
%A µ±Ç°µÄÉú³¤¾ØÕó
global kn
    CInE=setdiff(Ind1,Indxy);%ÌáÈ¡·ÇÉú³¤ºË±ß½ç¸ñµãË÷Òý
    En=rand(1,length(CInE));%¸÷¸öµ¥Ôª¾ù»á±»¸³ÖµÒ»¸öËæ»úÊý
    [i,j]=find(En<=p);%´óÓڸ÷½ÏòÉϵķ½ÏòÔö³¤¸ÅÂÊÊý£¬Ôò´Ëµ¥ÔªÍ¬Ñù³ÉΪһ¸öеÄÉú³¤Ïà
    pE=sub2ind([1,length(CInE)],i,j);%еÄÉú³¤ÏàË÷ÒýµÄϱê
    A(CInE(pE))=0; %¸üÐÂÉú³¤¾ØÕó
    [x,y]=find(A==0);%¸üÐÂÉú³¤ºË×ø±ê
    Ixy=sub2ind([kn,kn],x,y);%¸üÐÂÉú³¤ºËË÷Òý

1Â¥: Originally posted by dbb627 at 2011-10-27 09:20:06:
Õâ¸ö°å¿éÌû×ÓÌ«ÉÙ£¬·¢Ò»¸öÔö¼ÓµãÈËÆø£¬Ô­Ê¼ÎÊÌâ¼û
http://muchong.com/bbs/viewthread.php?tid=3723739
Ëã·¨ÈçÏ£º
Ëæ»úÉú³¤µÄ²½ÖèÈçÏ£º
1) ½¨Á¢²ÄÁÏÉú³¤µÄ³õʼÍø¸ñ£¬²¢É趨ËùÓÐÍø¸ñ²ÄÁϽÔΪÖ÷Ïࣻ ...

4Â¥: Originally posted by cxguo at 2011-10-28 12:41:57:
Â¥Ö÷¸ßÊÖ°¡£¬Çë½ÌÒÔÏ¿ÉÒÔÓÃMatalab»­Ô­×ÓÉú³¤Ä£Äâͼ²»£¿Èýλµº×´Éú³¤µÄ¡£