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

lxsoft110

我编这段程序的目的是为了模拟火灾中，一个房间里面的人员的撤离情况，运行出来，蓝色的点是障碍物、黑色是可以空地、绿色代表人。
我用的算法是：1.建立一个60*60的网格，并且定义边界和出口
                    2.给每个格子赋一个状态。（0，0，0）代表空地。（0，1，0）代表人，（0，0，1）代表障碍物。
                     3.规定每隔一段时间，计算人的状态（看前后左右，哪个格子离门口近，而且无障碍物），就向那个格子移动。
                      4。看所有人离开房间需要多少时间。
下面附上MATLAB程序：盼大牛解答！！不胜感激！！！
               clc
2 `- V- w9 k$ [  a) p) J1 Dclear
1 e& X  V/ W3 o2 dnx=60;
ny=60;
+ n! K* E" Q7 A& b%定义固定障碍物的状态
z=zeros(nx,ny);
peo=z;
%peoNew=z;
; l' w+ C$ }9 b- Q5 c6 D1 p  k+ Owall=z;
, f! h& ]7 B9 }, F! rwall(1,1:ny)=1;
wall(1:nx,1)=1;
9 u* Q0 h- i! {9 k! wwall(1:nx,ny)=1;
wall(nx,1:ny/2-2)=1;
wall(nx,ny/2+2:ny)=1;
2 S" y/ B' K" rwall(nx/4:nx/4+1,ny/4:ny/4+1)=1;
% N8 z- G# e, h) r4 y; ?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;
  s! K8 e+ z% u! x' [wall(3*nx/4:3*nx/4+1,3*ny/4:3*ny/4+1)=1;
%定义人的状态
/ v; V' w7 p2 r+ dfor x=1:1:12
! Z- \  a* Z$ R) d    for y=1:1:12
5 z, c4 f/ {. D* [        x1=round(rand(1)*59+1);
        x2=round(rand(1)*59+1);
& e" A4 {" p' m" F: b5 |. r        if wall(x1,x2)~=1
           peo(x1,x2)=1;
: W4 d2 R+ F. Z; S" ]/ u& n, M       end
) }6 }2 T4 W5 J( ?    end
end
%定义安全矩阵，指引人群运动
safe=zeros(60,60);
9 l% ~5 |4 s3 h& a" h) hsafe(nx,ny/2+1)=0;
safe(nx,ny/2+2)=0;
1 [4 D( }( i' M0 z% T! Lsafe(nx,ny/2)=0;
safe(nx,ny/2-1)=0;
safe(nx,ny/2-2)=0;
1 P' S+ m) l6 N5 Y* }( kfor i=nx-1:-1:2
    for j=ny/2+3:-1:2
, V6 @4 T( f0 L        safe(i,ny/2)=safe(i+1,ny/2)+1;
( c5 k. b7 ]# N" E2 b% h! L- j        safe(i,ny/2-1)=safe(i+1,ny/2-1)+1;
5 l4 @9 a* ^# P& F        safe(i,ny/2-2)=safe(i+1,ny/2-2)+1;
        safe(i,j)=safe(i,j+1)+1;
    end
end
0 F: ]/ L: T8 d1 q: W8 Y3 xfor 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;
8 A* ?; [$ m% N( X    end
end
        
( S& l! X8 g' b$ f4 Cimh=image(cat(3,z',peo',wall'));
5 ~* o7 D/ D+ ?* x" Aset(imh,'erasemode', 'none')
9 |2 d# A+ [% Z2 F2 W7 p+ I0 w9 gaxis equal
axis tight
%peoNew=peo;
# |& ~/ ^, [4 I) A2 j6 Rxind=2:59;
% o3 z) c' }- E( w, _, cyind=2:59;
for m=1:1000
    if peo(xind,yind)==1
$ v5 A5 N9 `+ |1 \       if safe(xind,yind-1)<safe(xind,yind)&wall(xind,yind-1)~=1
           peo(xind,yind-1)=1;
+ r: ]2 K3 I5 p1 U0 M           peo(xind,yind)=0;
4 C8 Q/ m+ G5 _  C# _/ P# W* f9 n       elseif safe(xind,yind+1)<safe(xind,yind)&wall(xind,yind+1)~=1
           peo(xind,yind+1)=1;
+ `& r- R) S4 W" u           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;
" @; Z$ Q5 M& @7 @9 E/ u. ]           peo(xind,yind)=0;
       end
         peo=peoNew;
         set(imh, 'cdata', cat(3,z',peo',wall') )
. B4 ?/ R6 d5 |$ Y% y9 A8 I         drawnow
         end
end

etzhu

实际上不用那么麻烦，
* ^& P  x+ I' b+ V. d! h5 p& t) X& V首先这个问题必然是满足 贪心 前提的，所以可以首先计算出每个人移动到门口的时间，然后排序，从小到大出门，考虑相等时的拥塞，就可以了。
* c& ~. u- Q) f. z当然，也有一点不足，这里没有考虑一个人挡在另一个人的路上的情况。

suxpert

up，友情帮顶。

sphereballball

现在是休息重要，别再动脑了

lxsoft110

2# etzhu
/ J5 ]: }! f  t  L. h7 U. a; k呵呵~兄弟。你还是没有说明我这个算法错误在哪里啊~

qsoft_cn

你调试一下可以发现，，if peo(xind,yind)==1 这个有问题，for循环里面的一次都没有执行。。。我正在想，应该有if find(peo(xind,yind)==1),，，厄，现在我该休息了。没时间再想了。元胞刚刚学习。。。。

wxl1988

我也要休息了~不看了~~好成绩·

cldndx

bless明天遇到顺手的题目！

victoriajh

注意休息啊 楼主

terrance_ho

up，友情帮顶。