人员疏散-元胞自动机-算法程序.pdfVIP

人员疏散-元胞自动机 clear;clc %定义 button plotbutton=uicontrol(style,pushbutton,... string,Run,... fontsize,12,... position,[100,400,50,20],... callback,run=1;); erasebutton=uicontrol(style,pushbutton,... string,Stop,... fontsize,12,... position,[300,400,50,20],... callback,freeze=1;); number=uicontrol(style,text,... string,1,... fontsize,12,... position,[20,400,50,20]); z=zeros(16,24); cells=z;%元胞矩阵 weith=16;length=24; x0=1;y0=12;%最佳出逃位置 peoplemidu=0.4;%人员密度 %危险度 for i=1:weith; for j=1:length; if(j10 j13 i1) pd(i,j)=j-1; end if(j10 j13 i==1) pd(i,j)=0; else pd(i,j)=sqrt((i-x0)^2+(j-y0)^2); end end end %%%%初始化选择矩阵 for i=1:weith; for j=1:length; choice{i,j}=[0 0]; end end %%%%初始化人群 for i=1:weith; for j=1:length; if(rand=peoplemidu) cells(i,j)=1; pv(i,j)=1; choice{i,j}=[i j]; end end end imh = imshow(cat(3,cells,cells,cells)); axis equal axis tight stop= 0; run = 0; freeze = 0; %%%%开始逃生 while (stop==0) if(run==1) for i=1:weith; for j=1:length; if(cells(i,j)~=0) if (j=10j13i==1) choice{i,j}=[0 0]; cells(i,j)=0; else [newi,newj]=findnew(i,j,pd); %找到选择位置 if(pd(i,j)pd(newi,newj)) if(same(newi,newj,choice)~=0) choice{i,j}=[newi newj]; else choice{i,j}=[i j]; end