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

lxsoft110

我编这段程序的目的是为了模拟火灾中，一个房间里面的人员的撤离情况，运行出来，蓝色的点是障碍物、黑色是可以空地、绿色代表人。
+ z( k- q# c7 G. g我用的算法是：1.建立一个60*60的网格，并且定义边界和出口
                    2.给每个格子赋一个状态。（0，0，0）代表空地。（0，1，0）代表人，（0，0，1）代表障碍物。
                     3.规定每隔一段时间，计算人的状态（看前后左右，哪个格子离门口近，而且无障碍物），就向那个格子移动。
+ a% Y0 j- ]* s" Z9 M+ o' j6 g                      4。看所有人离开房间需要多少时间。
下面附上MATLAB程序：盼大牛解答！！不胜感激！！！
               clc
" A. N  M( H( ^7 }8 W; [# Nclear
! T, D& L' [1 k+ Y' y& `nx=60;
ny=60;
%定义固定障碍物的状态
3 b. r; W9 t. o5 u! f& q5 tz=zeros(nx,ny);
peo=z;
%peoNew=z;
/ N8 P' q9 V+ W' ^1 t% |8 e* Twall=z;
) P4 j& q! P6 [2 F# M9 A0 lwall(1,1:ny)=1;
0 U) V# ]$ u* }4 rwall(1:nx,1)=1;
6 d2 U; x/ T$ q8 b3 w: H# s! ~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;
( y3 m+ ], |: Z2 u) K  R: i: pwall(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;
7 x  Z; g0 Y9 G% A- L4 p* c2 G%定义人的状态
4 _8 D5 c# |. V9 {1 M! A  Ffor x=1:1:12
( c, w- x8 }, J9 P3 M: r; f! U( f  J3 T/ ?    for y=1:1:12
  e" N; [* b4 i. H" `$ k3 L, M        x1=round(rand(1)*59+1);
1 t/ g% [: h) ^: h6 C+ Q        x2=round(rand(1)*59+1);
; A  q0 h: \% T6 Z        if wall(x1,x2)~=1
           peo(x1,x2)=1;
       end
* ~( S- Q" B! ~# X, w  K9 R    end
end
%定义安全矩阵，指引人群运动
safe=zeros(60,60);
1 T( w& ^0 b7 Q7 d/ psafe(nx,ny/2+1)=0;
& N* ^  f5 O% b* hsafe(nx,ny/2+2)=0;
8 D* v! U1 ~) e. g, E' M* zsafe(nx,ny/2)=0;
safe(nx,ny/2-1)=0;
safe(nx,ny/2-2)=0;
1 z$ c$ {; j' }: K/ Y% \9 A! N* kfor i=nx-1:-1:2
    for j=ny/2+3:-1:2
( @& O5 k7 r& H* y8 @: E        safe(i,ny/2)=safe(i+1,ny/2)+1;
) q) T  s  d$ S        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;
: ]- M/ E% o9 R) H# Q+ [0 n3 H    end
( i! c; {6 _' R# Nend
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;
( @) u9 u  u; }# ~4 n  m        safe(i,j)=safe(i,j-1)+1;
1 s2 R2 b+ a: q2 }8 r1 L    end
end
        
imh=image(cat(3,z',peo',wall'));
set(imh,'erasemode', 'none')
axis equal
3 h; S$ g# Q2 A' D' A8 Caxis tight
%peoNew=peo;
/ x) e+ b" {# Uxind=2:59;
. J. c3 D0 x+ y0 C  q. b5 j1 hyind=2:59;
" [, r$ S8 I5 a0 m5 C, f4 Yfor m=1:1000
    if peo(xind,yind)==1
       if safe(xind,yind-1)<safe(xind,yind)&wall(xind,yind-1)~=1
0 l0 m% t% M' r' ]1 X           peo(xind,yind-1)=1;
+ g) M4 ^: d$ m. t5 [0 D           peo(xind,yind)=0;
       elseif safe(xind,yind+1)<safe(xind,yind)&wall(xind,yind+1)~=1
! j7 |% Q! I# O           peo(xind,yind+1)=1;
) T& ~6 s8 D" k4 \           peo(xind,yind)=0;
: d; N! m! d+ N) D1 i       elseif safe(xind-1,yind)<safe(xind,yind)&wall(xind-1,yind)~=1
, y& Y% y# K4 h           peo(xind-1,yind)=1;
& b! H, L7 K1 I           peo(xind,yind)=0;
" N/ R9 ?3 t3 V- O9 n1 W( D       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') )
" @. f7 {+ J) }- [; I( i         drawnow
: P$ c. a' O) @4 C& V         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

实际上不用那么麻烦，
首先这个问题必然是满足 贪心 前提的，所以可以首先计算出每个人移动到门口的时间，然后排序，从小到大出门，考虑相等时的拥塞，就可以了。
- D( Y3 M; c# |当然，也有一点不足，这里没有考虑一个人挡在另一个人的路上的情况。