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

lxsoft110

我编这段程序的目的是为了模拟火灾中，一个房间里面的人员的撤离情况，运行出来，蓝色的点是障碍物、黑色是可以空地、绿色代表人。
我用的算法是：1.建立一个60*60的网格，并且定义边界和出口
) g" g8 b- c! r$ D/ Q( x                    2.给每个格子赋一个状态。（0，0，0）代表空地。（0，1，0）代表人，（0，0，1）代表障碍物。
2 {! K$ s1 X2 N+ k; ?2 r  t5 Z                     3.规定每隔一段时间，计算人的状态（看前后左右，哪个格子离门口近，而且无障碍物），就向那个格子移动。
                      4。看所有人离开房间需要多少时间。
下面附上MATLAB程序：盼大牛解答！！不胜感激！！！
               clc
. E0 W1 F4 ~) Rclear
nx=60;
* v% a# F; N& M5 j% T) K5 zny=60;
' H8 @" K. M. ]3 x%定义固定障碍物的状态
0 U" U( q6 k* [5 u2 kz=zeros(nx,ny);
' p- L7 l6 D" F, E1 |peo=z;
%peoNew=z;
wall=z;
wall(1,1:ny)=1;
& U  F9 ~4 r9 cwall(1:nx,1)=1;
. |( ?6 u# |2 |wall(1:nx,ny)=1;
4 s$ E/ s; X7 W: p5 J* Z4 I) Qwall(nx,1:ny/2-2)=1;
# ~6 N( |$ o( U0 Iwall(nx,ny/2+2:ny)=1;
8 Q+ A3 K4 A6 y: [) {$ B( Mwall(nx/4:nx/4+1,ny/4:ny/4+1)=1;
3 S. K# U& b! dwall(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;
( Y# [4 B/ K. V& H4 T/ q%定义人的状态
& p  b& K! ~5 O9 E( a- I  ^3 Efor x=1:1:12
    for y=1:1:12
7 I; W' F) H) d( X# W1 j        x1=round(rand(1)*59+1);
- l1 {2 |/ n8 u        x2=round(rand(1)*59+1);
# V6 C% v. F+ K9 c! W        if wall(x1,x2)~=1
& z6 }  L, L% v% C+ @* s           peo(x1,x2)=1;
       end
) ?) W9 x. g- N% l* Z    end
: K+ N3 v. L" z% Kend
/ z1 J! I; Z9 u* P& B5 w4 a, ^  l%定义安全矩阵，指引人群运动
safe=zeros(60,60);
safe(nx,ny/2+1)=0;
safe(nx,ny/2+2)=0;
  U6 c, Y$ I( E9 m7 [/ Ksafe(nx,ny/2)=0;
( V# F+ K' u6 ^/ {7 l8 Msafe(nx,ny/2-1)=0;
2 Y% v  D' I- e5 t- E/ R0 f3 usafe(nx,ny/2-2)=0;
% T/ s4 v: H- ~: Xfor i=nx-1:-1:2
3 R$ L7 k* `) I    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;
" X( [/ I" x) n    end
end
$ k# v" T5 K0 Q6 f; }for i=nx-1:-1:2
: a" E; ^" G7 E' [! F/ @. X5 U    for j=ny/2+3:1:ny-1
        safe(i,ny/2+1)=safe(i+1,ny/2+1)+1;
  G+ G( v6 ^! [: q        safe(i,ny/2+2)=safe(i+1,ny/2+2)+1;
5 D7 P8 T+ ~5 u. N4 A2 \+ j        safe(i,j)=safe(i,j-1)+1;
$ b6 y. @* {; \    end
6 c& J& @0 r! {4 Q. wend
. J0 n( |& L! l2 z        
imh=image(cat(3,z',peo',wall'));
+ I9 D2 Q+ M. n9 S- ~1 h) vset(imh,'erasemode', 'none')
axis equal
axis tight
+ Y! s* U1 ~) L9 l%peoNew=peo;
xind=2:59;
$ O0 c" k+ b" u$ q* \2 Uyind=2:59;
for m=1:1000
    if peo(xind,yind)==1
! i" o; E1 ^% B0 g8 [       if safe(xind,yind-1)<safe(xind,yind)&wall(xind,yind-1)~=1
2 c- k* u: a6 M; o  L           peo(xind,yind-1)=1;
7 T, P6 w( S* k+ U6 y           peo(xind,yind)=0;
. ~$ `& I1 A1 x       elseif safe(xind,yind+1)<safe(xind,yind)&wall(xind,yind+1)~=1
           peo(xind,yind+1)=1;
: }* c1 g4 \1 a           peo(xind,yind)=0;
       elseif safe(xind-1,yind)<safe(xind,yind)&wall(xind-1,yind)~=1
           peo(xind-1,yind)=1;
9 v% V8 w5 {6 n8 B! m           peo(xind,yind)=0;
8 W- r6 ~0 J2 n( `9 r1 Y6 V       elseif  safe(xind+1,yind)<safe(xind,yind)&wall(xind+1,yind)~=1
0 e! z" Q" d0 V1 n3 L           peo(xind+1,yind)=1;
           peo(xind,yind)=0;
       end
  n( w  U5 b6 d# u4 f; D' _. e         peo=peoNew;
         set(imh, 'cdata', cat(3,z',peo',wall') )
- @& t7 i' m$ a& a* t, }         drawnow
         end
end

terrance_ho

up，友情帮顶。

victoriajh

注意休息啊 楼主

cldndx

bless明天遇到顺手的题目！

wxl1988

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

qsoft_cn

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

lxsoft110

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

sphereballball

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

suxpert

up，友情帮顶。

etzhu

实际上不用那么麻烦，
首先这个问题必然是满足 贪心 前提的，所以可以首先计算出每个人移动到门口的时间，然后排序，从小到大出门，考虑相等时的拥塞，就可以了。
" ?. }+ Y0 D7 Z' e) m- Y当然，也有一点不足，这里没有考虑一个人挡在另一个人的路上的情况。