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

lxsoft110

我编这段程序的目的是为了模拟火灾中，一个房间里面的人员的撤离情况，运行出来，蓝色的点是障碍物、黑色是可以空地、绿色代表人。
我用的算法是：1.建立一个60*60的网格，并且定义边界和出口
                    2.给每个格子赋一个状态。（0，0，0）代表空地。（0，1，0）代表人，（0，0，1）代表障碍物。
                     3.规定每隔一段时间，计算人的状态（看前后左右，哪个格子离门口近，而且无障碍物），就向那个格子移动。
  [" C+ i- }8 t' P4 ]6 }% _                      4。看所有人离开房间需要多少时间。
下面附上MATLAB程序：盼大牛解答！！不胜感激！！！
1 L# }7 Y# w. R& b7 m- d& _               clc
+ E, n  X8 [! A, c) w# fclear
9 l( @0 F$ X; u6 ?5 X9 Y) p5 R6 Onx=60;
ny=60;
( J% R' P! t" M. A%定义固定障碍物的状态
z=zeros(nx,ny);
! T. {" F. h2 X. s! W7 F, B* L* Xpeo=z;
%peoNew=z;
5 }0 u! o; ~* Ewall=z;
6 c0 X# X, a2 r4 f- [1 Rwall(1,1:ny)=1;
. }( P& S5 M6 x3 u$ h4 Awall(1:nx,1)=1;
wall(1:nx,ny)=1;
6 f  W& h; E% }7 {1 e5 h5 Y: dwall(nx,1:ny/2-2)=1;
0 y* p  v/ j3 j! iwall(nx,ny/2+2:ny)=1;
/ j$ h" Y) L& y# F: A5 ^2 x/ zwall(nx/4:nx/4+1,ny/4:ny/4+1)=1;
* o2 v% o) c) Gwall(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
( r0 u, @9 {2 T! e; a+ ?    for y=1:1:12
  _; g$ A+ e6 t0 d        x1=round(rand(1)*59+1);
2 Q, ]2 E9 }4 ~        x2=round(rand(1)*59+1);
  ~; z" w3 A! A( C  l3 b) J        if wall(x1,x2)~=1
3 r! M% |& z+ I+ f6 D& w, C) k$ T$ K0 i% u           peo(x1,x2)=1;
. ~3 @. }# ]- s       end
    end
: R4 |. e) K" N1 K+ }! u0 ?# H2 eend
3 w/ p8 B, {) Q$ _# \%定义安全矩阵，指引人群运动
safe=zeros(60,60);
safe(nx,ny/2+1)=0;
safe(nx,ny/2+2)=0;
safe(nx,ny/2)=0;
1 h! X' h$ I9 m  s/ x. J  lsafe(nx,ny/2-1)=0;
safe(nx,ny/2-2)=0;
+ ?9 d" T# f' L! L6 Lfor i=nx-1:-1:2
/ B! o, y  W/ d5 k    for j=ny/2+3:-1:2
5 ?- a2 w# Z! O7 x. T( G. p        safe(i,ny/2)=safe(i+1,ny/2)+1;
        safe(i,ny/2-1)=safe(i+1,ny/2-1)+1;
1 k! Z% H& e( H        safe(i,ny/2-2)=safe(i+1,ny/2-2)+1;
4 o& q2 w* ?2 B5 o# w1 e        safe(i,j)=safe(i,j+1)+1;
0 J2 r* `# b+ l    end
end
7 |1 x, d0 A0 K& Afor i=nx-1:-1:2
    for j=ny/2+3:1:ny-1
* ], J4 ^  B$ C2 U% b, `1 Y        safe(i,ny/2+1)=safe(i+1,ny/2+1)+1;
: |  h0 p) J" X- K! H( f        safe(i,ny/2+2)=safe(i+1,ny/2+2)+1;
5 j/ \  Q% H8 c        safe(i,j)=safe(i,j-1)+1;
    end
  l- R* w4 W: ]8 P/ N1 z% yend
. o+ o5 F2 ^+ Z' c9 d( a" \        
imh=image(cat(3,z',peo',wall'));
  c  e  H2 ^6 Y( g! mset(imh,'erasemode', 'none')
axis equal
axis tight
%peoNew=peo;
4 I% A/ _( K4 _. g0 vxind=2:59;
4 Q: J: N3 Q8 B7 ]8 Cyind=2:59;
" R* t7 L' M* W7 A6 |) Ofor 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;
  O( O: _: g" B% r# Z+ o; r           peo(xind,yind)=0;
' Z, z: E0 ?7 G: W# A" [       elseif safe(xind,yind+1)<safe(xind,yind)&wall(xind,yind+1)~=1
. t( h- q- B! w% [           peo(xind,yind+1)=1;
' h% e& m. O2 u) `* |, x           peo(xind,yind)=0;
       elseif safe(xind-1,yind)<safe(xind,yind)&wall(xind-1,yind)~=1
9 C0 F$ X% i1 X% ~           peo(xind-1,yind)=1;
           peo(xind,yind)=0;
7 h* N  k; i4 I6 ?- T       elseif  safe(xind+1,yind)<safe(xind,yind)&wall(xind+1,yind)~=1
, H0 ?- M3 a% Z; f0 q% m           peo(xind+1,yind)=1;
           peo(xind,yind)=0;
       end
6 T+ _! E% K- ^6 l         peo=peoNew;
         set(imh, 'cdata', cat(3,z',peo',wall') )
# ?5 |  u* d1 S% ]" D8 D' o! H! G/ L         drawnow
! _5 Q+ y" H/ p2 _         end
end

etzhu

实际上不用那么麻烦，
首先这个问题必然是满足 贪心 前提的，所以可以首先计算出每个人移动到门口的时间，然后排序，从小到大出门，考虑相等时的拥塞，就可以了。
; h4 u7 i8 F9 I- O/ l! i2 Z当然，也有一点不足，这里没有考虑一个人挡在另一个人的路上的情况。

suxpert

up，友情帮顶。

sphereballball

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

lxsoft110

2# etzhu
呵呵~兄弟。你还是没有说明我这个算法错误在哪里啊~

qsoft_cn

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

wxl1988

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

cldndx

bless明天遇到顺手的题目！

victoriajh

注意休息啊 楼主

terrance_ho

up，友情帮顶。