查看: 5561|回复: 9

谁能真正解决这个问题，送20个矩阵币

字体大小: 正常放大

18 主题	3 听众	193 积分

升级 46.5%

该用户从未签到

电梯直达

1^#

发表于 2009-2-5 21:21 |只看该作者 |倒序浏览

|招呼Ta 关注Ta

我编这段程序的目的是为了模拟火灾中，一个房间里面的人员的撤离情况，运行出来，蓝色的点是障碍物、黑色是可以空地、绿色代表人。
我用的算法是：1.建立一个60*60的网格，并且定义边界和出口
                  2.给每个格子赋一个状态。（0，0，0）代表空地。（0，1，0）代表人，（0，0，1）代表障碍物。
                  3.规定每隔一段时间，计算人的状态（看前后左右，哪个格子离门口近，而且无障碍物），就向那个格子移动。
                  4。看所有人离开房间需要多少时间。
下面附上MATLAB程序：盼大牛解答！！不胜感激！！！
            clc
clear
nx=60;
ny=60;
%定义固定障碍物的状态
z=zeros(nx,ny);
peo=z;
%peoNew=z;
wall=z;
wall(1,1:ny)=1;
wall(1:nx,1)=1;
wall(1:nx,ny)=1;
wall(nx,1:ny/2-2)=1;
wall(nx,ny/2+2:ny)=1;
wall(nx/4:nx/4+1,ny/4:ny/4+1)=1;
wall(3*nx/4:3*nx/4+1,ny/4:ny/4+1)=1;
wall(nx/4:nx/4+1,3*ny/4:3*ny/4+1)=1;
wall(3*nx/4:3*nx/4+1,3*ny/4:3*ny/4+1)=1;
%定义人的状态
for x=1:1:12
for y=1:1:12
      x1=round(rand(1)*59+1);
      x2=round(rand(1)*59+1);
      if wall(x1,x2)~=1
         peo(x1,x2)=1;
   end
end
end
%定义安全矩阵，指引人群运动
safe=zeros(60,60);
safe(nx,ny/2+1)=0;
safe(nx,ny/2+2)=0;
safe(nx,ny/2)=0;
safe(nx,ny/2-1)=0;
safe(nx,ny/2-2)=0;
for i=nx-1:-1:2
for j=ny/2+3:-1:2
      safe(i,ny/2)=safe(i+1,ny/2)+1;
      safe(i,ny/2-1)=safe(i+1,ny/2-1)+1;
      safe(i,ny/2-2)=safe(i+1,ny/2-2)+1;
      safe(i,j)=safe(i,j+1)+1;
end
end
for i=nx-1:-1:2
for j=ny/2+3:1:ny-1
      safe(i,ny/2+1)=safe(i+1,ny/2+1)+1;
      safe(i,ny/2+2)=safe(i+1,ny/2+2)+1;
      safe(i,j)=safe(i,j-1)+1;
end
end

imh=image(cat(3,z',peo',wall'));
set(imh,'erasemode', 'none')
axis equal
axis tight
%peoNew=peo;
xind=2:59;
yind=2:59;
for m=1:1000
if peo(xind,yind)==1
   if safe(xind,yind-1)<safe(xind,yind)&wall(xind,yind-1)~=1
         peo(xind,yind-1)=1;
         peo(xind,yind)=0;
   elseif safe(xind,yind+1)<safe(xind,yind)&wall(xind,yind+1)~=1
         peo(xind,yind+1)=1;
         peo(xind,yind)=0;
   elseif safe(xind-1,yind)<safe(xind,yind)&wall(xind-1,yind)~=1
         peo(xind-1,yind)=1;
         peo(xind,yind)=0;
   elseif  safe(xind+1,yind)<safe(xind,yind)&wall(xind+1,yind)~=1
         peo(xind+1,yind)=1;
         peo(xind,yind)=0;
   end
      peo=peoNew;
      set(imh, 'cdata', cat(3,z',peo',wall') )
      drawnow
      end
end

zan