(海港系统卸货模型).docVIP

§8—4 海港系统卸载货物的模拟 8.4.1 关于原问题的模拟 读书P119-----P122之“流程框图” 关于下面主程序（syp123）的说明： g 是100个数，服从“15到145的均匀分布”，代表100个船到来的间隔时间； G 是100个数，代表100个船的到来时刻； f 是100个数，服从“45到90的均匀分布”，代表100个船的卸载时间； F 是100个数，代表100个船的离开时刻 clear g=15+130*rand(1,100); G(1)=g(1); for i=2:100 G(i)=G(i-1)+g(i); end f=45+45*rand(1,100); F(1)=g(1)+f(1); for i=2:100 F(i)=max(F(i-1),G(i))+f(i); end ZXZSJ=sum(f);ZKXSJ=F(100)-ZXZSJ; KXL=ZKXSJ/F(100) dd(1)=0;tl(1)=dd(1)+f(1); for i=2:100 dd(i)=max(F(i-1)-G(i),0); tl(i)=dd(i)+f(i); end pjdd=mean(dd),pjtl=mean(tl) zcdd=max(dd),zctl=max(tl) 由于随机因素，每次运行，机器会产生不同的随机数，故计算结果也会不同。下表给出了十次计算得到的结果： 次 平均停留 最长停留 平均等待 最长等待 空闲率 1 101 322 36 267 0.22 2 88 200 20 146 0.20 3 112 275 44 189 0.16 4 134 296 65 246 0.11 5 123 269 56 181 0.16 6 117 271 50 218 0.13 7 88 202 21 146 0.22 8 92 193 25 127 0.19 9 103 256 37 189 0.15 10 106 291 37 223 0.13 书P124之“想(1)”：时间间隔改为“10到90的均匀分布”，则结果如下： 次 平均停留 最长停留 平均等待 最长等待 空闲率 1 781 1389 713 1323 0.0180 2 752 1358 686 1307 0.0053 3 822 1524 753 1448 0.0079 8.4.2 根据统计数据确定间隔时间、卸载时间 看书P124之“表8.9”，如何根据统计规律来产生间隔时间、卸栽时间随机数呢？ 法一：书P125----P127之方法。 法二：命令g=rand(1,100)产生0到1的100个均匀数，逐个检验每个g(i)，若它介于0---0.009则令g(i)=20就是第i个间隔时间，若它介于0.009---0.038则令g(i)=30就是第i个间隔时间，若它介于0.038---0.073则令g(i)=40就是第i个间隔时间，等等。以同样方法产生卸载时间。 下面两个函数文件就是用“法二”的方法，根据“表8.9”的数据，产生100个船的到来间隔g、卸载时间f ： function g=syp125hswj1 g=rand(1,100); for i=1:100 if g(i)0.0009 g(i)=20; elseif g(i)0.038 g(i)=30; elseif g(i)0.073 g(i)=40; elseif g(i)0.124 g(i)=50; elseif g(i)0.214 g(i)=60; elseif g(i)0.375 g(i)=70; elseif g(i)0.575 g(i)=80; elseif g(i)0.747 g(i)=90; elseif g(i)0.872 g(i)=100; elseif g(i)0.943 g(i)=110; elseif g(i)0.980 g(i)=120; elseif g(i)0.997 g(i)=130; else g(i)=140; end end 第二个： function f=syp125hswj2 f=rand(1,100); for i=1:100 if f(i)0.017 f(i)=47; elseif f(i)0.062 f(i)=52;