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

lxsoft110

我编这段程序的目的是为了模拟火灾中，一个房间里面的人员的撤离情况，运行出来，蓝色的点是障碍物、黑色是可以空地、绿色代表人。
我用的算法是：1.建立一个60*60的网格，并且定义边界和出口
                    2.给每个格子赋一个状态。（0，0，0）代表空地。（0，1，0）代表人，（0，0，1）代表障碍物。
$ y3 s) F( h' G6 n+ K; a" N# N                     3.规定每隔一段时间，计算人的状态（看前后左右，哪个格子离门口近，而且无障碍物），就向那个格子移动。
                      4。看所有人离开房间需要多少时间。
/ g' A( a* V# N7 K& Q, W# \, Q下面附上MATLAB程序：盼大牛解答！！不胜感激！！！
               clc
clear
9 {/ l- T7 _( Nnx=60;
( V- |/ `/ [8 q1 ^0 Zny=60;
%定义固定障碍物的状态
- Z& U* i4 d3 C# @( Bz=zeros(nx,ny);
: ^2 {% S0 M3 t. I* j6 C/ l' Ypeo=z;
, x1 @: \, E3 U0 I%peoNew=z;
- b+ E' ~( G$ X* g0 ]wall=z;
wall(1,1:ny)=1;
( |$ b9 m& A8 n) L+ N" Zwall(1:nx,1)=1;
5 L  z$ T- o8 O$ i* 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;
wall(3*nx/4:3*nx/4+1,ny/4:ny/4+1)=1;
! I# K! F5 t# e# `3 U# owall(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;
" q3 f0 e; e* r, m8 N%定义人的状态
for x=1:1:12
" D, h1 w7 B' N5 N7 v; p& j    for y=1:1:12
        x1=round(rand(1)*59+1);
% f8 N! J" b5 l7 I% q        x2=round(rand(1)*59+1);
7 O: E! g' Z) `7 p+ R        if wall(x1,x2)~=1
( N( c4 `3 P+ r           peo(x1,x2)=1;
       end
9 n& m  H6 i% `+ s    end
end
%定义安全矩阵，指引人群运动
safe=zeros(60,60);
safe(nx,ny/2+1)=0;
& F' W& Y  b/ X* k% osafe(nx,ny/2+2)=0;
safe(nx,ny/2)=0;
safe(nx,ny/2-1)=0;
safe(nx,ny/2-2)=0;
3 P2 o: S6 [# Y" ]* `1 Lfor i=nx-1:-1:2
    for j=ny/2+3:-1:2
+ `/ R' P$ ^6 F6 o( b        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;
/ ], O+ g! H  {  P    end
# p: q( E4 F& h: g0 ]" F8 |end
for i=nx-1:-1:2
    for j=ny/2+3:1:ny-1
5 w/ ]- t" K: L. n9 f+ K        safe(i,ny/2+1)=safe(i+1,ny/2+1)+1;
        safe(i,ny/2+2)=safe(i+1,ny/2+2)+1;
* X2 r7 y' Y2 K5 }: e: x        safe(i,j)=safe(i,j-1)+1;
; S2 a5 }7 @( \/ C/ v/ P    end
& M! S. s: F2 x7 i8 i% G9 S- u8 Dend
( o* u) u: w" u1 C( L- n        
imh=image(cat(3,z',peo',wall'));
set(imh,'erasemode', 'none')
axis equal
# d% P# b1 n" o7 k0 Laxis tight
%peoNew=peo;
xind=2:59;
yind=2:59;
for 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;
           peo(xind,yind)=0;
       elseif safe(xind,yind+1)<safe(xind,yind)&wall(xind,yind+1)~=1
, U2 K" J( t# _1 g           peo(xind,yind+1)=1;
  S4 k' n( w# P  K; i           peo(xind,yind)=0;
& q6 ?# O: y1 g/ T5 D* K; n' E1 Z       elseif safe(xind-1,yind)<safe(xind,yind)&wall(xind-1,yind)~=1
" ?+ r5 y4 F  I5 O" r! M           peo(xind-1,yind)=1;
           peo(xind,yind)=0;
       elseif  safe(xind+1,yind)<safe(xind,yind)&wall(xind+1,yind)~=1
! X) ?2 }% D4 K. }           peo(xind+1,yind)=1;
           peo(xind,yind)=0;
       end
6 g4 {% B: T& H: L5 S" t7 x         peo=peoNew;
# c+ w* u( V: L, B) I         set(imh, 'cdata', cat(3,z',peo',wall') )
  a- b$ u1 Q3 s8 p' _* Q         drawnow
         end
end

etzhu

实际上不用那么麻烦，
* J- u1 h& R& s; O5 d/ |- j首先这个问题必然是满足 贪心 前提的，所以可以首先计算出每个人移动到门口的时间，然后排序，从小到大出门，考虑相等时的拥塞，就可以了。
3 N; U8 x5 t: 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，友情帮顶。