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

lxsoft110

我编这段程序的目的是为了模拟火灾中，一个房间里面的人员的撤离情况，运行出来，蓝色的点是障碍物、黑色是可以空地、绿色代表人。
$ r, b  c0 t1 W' ]3 r我用的算法是：1.建立一个60*60的网格，并且定义边界和出口
                    2.给每个格子赋一个状态。（0，0，0）代表空地。（0，1，0）代表人，（0，0，1）代表障碍物。
, ]- m) @; e3 R                     3.规定每隔一段时间，计算人的状态（看前后左右，哪个格子离门口近，而且无障碍物），就向那个格子移动。
                      4。看所有人离开房间需要多少时间。
下面附上MATLAB程序：盼大牛解答！！不胜感激！！！
$ p2 }) r( R5 F7 w( w) _, Q               clc
" V* `& V2 I4 \' v! pclear
9 H7 O2 y( q$ R$ [: z! M3 Qnx=60;
# C3 T" q  i) T- x5 B. `ny=60;
%定义固定障碍物的状态
. u) T  B6 \- Jz=zeros(nx,ny);
- C! v+ v. R- k, ypeo=z;
! o$ c6 }, o+ C0 L  b. m; ], N%peoNew=z;
8 p  m/ q# r; b6 F' ]. O0 Zwall=z;
3 N  j4 F  O9 nwall(1,1:ny)=1;
wall(1:nx,1)=1;
- I# T. D0 F+ Rwall(1:nx,ny)=1;
wall(nx,1:ny/2-2)=1;
5 M: U: B  s+ G" X( Gwall(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;
5 g/ E' _% Q( x1 u+ jwall(3*nx/4:3*nx/4+1,3*ny/4:3*ny/4+1)=1;
%定义人的状态
. R" ]. O( `7 Z( P. ifor 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;
# N4 E1 ^) V2 c       end
4 S- j0 U7 L. D    end
end
8 E/ f! B: {& l; g%定义安全矩阵，指引人群运动
, J) X- e' h0 ]- t/ Msafe=zeros(60,60);
safe(nx,ny/2+1)=0;
5 u* r% x4 r4 {/ Z2 ]' Usafe(nx,ny/2+2)=0;
+ _! d5 d0 I; U% [0 Usafe(nx,ny/2)=0;
safe(nx,ny/2-1)=0;
safe(nx,ny/2-2)=0;
& `6 b  x( K' H5 n+ U7 {/ Ifor 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;
( ?9 f2 I( h$ R5 p+ }0 U        safe(i,ny/2-2)=safe(i+1,ny/2-2)+1;
        safe(i,j)=safe(i,j+1)+1;
3 ^9 l6 ]7 R2 Q) L( h) X: C    end
end
8 c& R, r1 M- `for i=nx-1:-1:2
    for j=ny/2+3:1:ny-1
: I5 m! u7 Z: @) C6 k& G8 o& ^        safe(i,ny/2+1)=safe(i+1,ny/2+1)+1;
/ k9 Z! S: ]) F  c        safe(i,ny/2+2)=safe(i+1,ny/2+2)+1;
& W9 Q6 W; c' ~7 k, I" x$ L        safe(i,j)=safe(i,j-1)+1;
    end
# J  x3 V' V7 M3 G/ M3 {/ n% H8 ?end
        
" v% l6 ?  v& Limh=image(cat(3,z',peo',wall'));
# h2 R. ?' u, L8 l, N! B- Bset(imh,'erasemode', 'none')
* `1 A2 G* T! B+ s; ^  b5 `! ]0 L0 [axis equal
axis tight
& s/ T: \; K1 l, @/ s8 M%peoNew=peo;
/ q: I5 R: {5 @$ R4 Gxind=2:59;
0 W% ]9 J! z; r5 q. K! m3 Ryind=2:59;
for m=1:1000
    if peo(xind,yind)==1
  G- J$ x! [9 `. K7 _% h       if safe(xind,yind-1)<safe(xind,yind)&wall(xind,yind-1)~=1
           peo(xind,yind-1)=1;
7 c( {7 j. S7 t9 b           peo(xind,yind)=0;
       elseif safe(xind,yind+1)<safe(xind,yind)&wall(xind,yind+1)~=1
( V9 e# ]9 C5 V9 B* T6 ~           peo(xind,yind+1)=1;
1 ]! e; u- X* G6 b; i           peo(xind,yind)=0;
( M/ D* ]0 W6 @3 Y- X1 @1 [. u& W       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;
6 f0 f& ?& A  V1 _         set(imh, 'cdata', cat(3,z',peo',wall') )
# m) J4 w. f+ H& q1 m- L         drawnow
         end
! r$ q" U) R# R2 J3 jend

terrance_ho

up，友情帮顶。

victoriajh

注意休息啊 楼主

cldndx

bless明天遇到顺手的题目！

wxl1988

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

qsoft_cn

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

lxsoft110

2# etzhu
( }: |8 f5 ^) M+ j7 Q# s, {7 X2 J7 o呵呵~兄弟。你还是没有说明我这个算法错误在哪里啊~

sphereballball

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

suxpert

up，友情帮顶。

etzhu

实际上不用那么麻烦，
3 a/ r! V% M) c9 i% t* L7 S- s首先这个问题必然是满足 贪心 前提的，所以可以首先计算出每个人移动到门口的时间，然后排序，从小到大出门，考虑相等时的拥塞，就可以了。
+ _6 A& s0 ]% D' i9 z当然，也有一点不足，这里没有考虑一个人挡在另一个人的路上的情况。