第6讲蒙特卡罗方法与计算机模拟(第2次课)解说.ppt

离散系统的计算机模拟实例一 范例 海港系统的卸载货物问题(p110-111,p119-129) 海港系统的卸载货物-问题1 问题1 海港系统的卸载货物 (p122) 程序片段(船只到港时间均匀分布,船只卸货时间均匀分布) ShipBetweenTime(1)=unifrnd(15,145,1,1); %船只到港间隔时间随机化(均匀分布) ShipUnloadTime(1)=unifrnd(45,90,1,1); %船只卸货时间随机化(均匀分布) 通用程序haibor1.m可实现多次模拟，结果保存到H1.txt delete H1.txt %清除历史数据 harbor1(100,15,145,45,90) load H1.txt; Hmean1=mean(H1); %导入H1并按列取平均值 海港系统的卸载货物-问题2 问题2 海港系统的卸载货物 (p124) 程序片段(船只到港时间指数分布,船只卸货时间均匀分布) ShipBetweenTime(1)=exprnd(60,1,1); %船只到港间隔时间随机化(指数分布) ShipUnloadTime(1)=unifrnd(45,90,1,1); %船只卸货时间随机化(均匀分布) 通用程序haibor2.m可实现多次模拟，结果保存到H2.txt delete H2.txt %清除历史数据 harbor2(100,60,45,90) load H2.txt; Hmean2=mean(H2); %导入H2并按列取平均值 海港系统的卸载货物-问题3-[1] 问题3 海港系统的卸载货物 (p124-129) 程序片段(船只到港时间离散分布,船只卸货时间离散分布) [1] 编写船只到港间隔离散累积分布函数并作阶梯图： xs=15:10:145; for i=1:length(xs)-1,x(i)=(xs(i)+xs(i+1))/2;end px=[0.009,0.029,0.035,0.051,0.090,0.161,0.200,0.172,0.125,0.071,0.037,0.017,0.003]; Fx=0; for i=1:length(px), Fx=[Fx,sum(px(1:i))]; end plot([10,x],Fx,-rs); hold on; stairs([0,x-5,145],[Fx,1]); set(gca,xtick,0:5:145); set(gca,xgrid,on); axis tight; 海港系统的卸载货物-问题3-[2] 问题3 海港系统的卸载货物 (p124-129) 程序片段(船只到港时间离散分布,船只卸货时间离散分布) [2] 编写船只到港间隔离散累积分布反函数并作线性插值： Fxi=0:0.001:1-eps; xi=interp1(Fx,[0,x],Fxi,linear); index=[find(xixs(1)),find(xixs(end))]; xi(index)=[]; Fxi(index)=[]; r=rand(1,n); rnd=[]; for i=1:n index=find(r(i)=Fxi); rnd=[rnd,xi(index(1))]; end %以上程序已编写通用M函数文件 harborrnd(xs,px,n) %即给出n个满足离散分布(xs,px)的船只到港间隔随机数 海港系统的卸载货物-问题3-[3] 问题3 海港系统的卸载货物 (p124-129) 程序片段(船只到港时间离散分布,船只卸货时间离散分布) [3] 编写船只卸货时间离散累积分布函数并作阶梯图： xs=45:5:90; for i=1:length(xs)-1,x(i)=(xs(i)+xs(i+1))/2;end px=[0.017,0.045,0.095,0.086,0.130,0.185,0.208,0.143,0.091]; Fx=0; for i=1:length(px), Fx=[Fx,sum(px(1:i))]; end plot([40,x],Fx,-rs); hold on; stairs([40,x-2.5,90],[Fx,1]); set(gca,xtick,40:2.5:90); set(gca,xgrid,on); axis tight; 海港系统的卸载货物-问题3-[4] 问题3 海港系统的卸载货物 (p124-129) 程序片段(船只到港时间离散分布,船只卸货时间离散分布) [4] 编写船只卸货时间离散累积分布反函数并作线性插值： Fxi=0:0.001:1-eps; xi=interp1(Fx,[0,x],Fxi,linear);