港口系统仿真课程设计交运172王星宇2017106111581.在计算机上利用线性同余法或乘同余法产生[0,1]区间独立均匀分布的随机数序列。(A)要求产生2组随机数序列，每组1000个。(B)计算每组随机数的均值，并且以0.1为间隔，(0,0.1)，(0.1,0.2),(0.2,0.3),...画出随机数落在每个子区间内的频率分布曲线。解：⑴用a=9,c=7,m=2M5,x0=7生成一个1000行1列的矩阵，用循环语句生成1000个0到1之间的随机数，用mean函数求平均值并用histfit进行作图，如图1以及图2。图35975►Desktap►1Untitled.lm^=9;e=7;m.=2'15;±0=7;B-zeras(1OOLL1)P感-C;\UItimefun.m图图22）以下是第二组随机数的生成，流程与第一小组随机数相同，如图3和图4。图3a.=LL;e=O：tii=2'15：D^eroaflQQO,1);n-1:Iwhilen<-1000:n.=n*l;x=r™((3l*K0+c),nJ;■endmeantn)histritfriLdli?5»■Destttap►4110&8必砂整・山titled歹5=s/m;D(n-1.l)-s:图42.为模拟船舶的装卸过程，首先需要模拟船舶的到港过程，假设泊位充足，且船舶的到达服从每天平均3.6艘船的泊松分(A)以T1产生的第一组随机数为基础，按照上述分布特点产生1000艘船舶的到港时间间隔(以min为单位),并计算你所产生的这1000个到达时间间隔的平均值。(B)以100min为划分间隔，画出产生到达时间间隔在各区间的频率分布图解：首先将每天3.6艘船的单位进行转换，然后建立循环语句，因为船舶达港服从泊松到达过程且符合指数分布，所以到达间隔时间可以用R=1-e"a*x),然后求出平均数并作图，如图5和图6。图6pinaju=y.eJlan职=3,6/14401titervBl=zer&a(1(XX)・]J.D=l；rtillerK=lWD；mterval?：hl£t!3t(ititervai.:K)：i图5ans-图73.为仿真船舶的装卸过程，需模拟船舶的装卸过程。为此，研究人员统计了200艘船舶的装卸所需时间，如下表所示：区间(小时)2~44~66~88~1010~11频数4050603020(A)根据该数据拟合出装卸服务时间这个随机变量的累积分布函数(B)以T1产生的第二组随机数为基础，按照上述统计规律模拟产生1000艘船舶的装卸服务时间(单位:min),并计算所产生的1000艘船舶装卸服务时间的平均值。(C)以100min为区间间隔，画出产生数据在各区间的频率分布图。解：可以套用累积分布函数得到函数tan(i)=(x(i+1)-x(i))/(y(i+1)-y(i))并利用switch语句生成装卸时间。如图7与图8。timefun.mUntrtled.1.mUntiTled^^Untitled4feUntitled5*+7=10,0.2,0.45,0.75,0.9,11:^.an=zer05(.5,1j;fori=l:5tfa(i+1)-jc(i)>.;-endstevedare=zeros(1000.1).*60;fori-l:1000switchceil(D(I.)4=100)casenum2cellID:20)stevedore(.i.1=tmnU*D⑴:casenun2cell(21:45)stevedoreU)=tan⑵.*(D(i)-0.2)+x⑵;casenun2cell(46:75)stevedore(i)=tan(3).*(D(i)-0.45)七(3);casemm2匚ell:g。)stevedore(i)=tan〔4).*(D(i)-0.75)+玄(4);casenum2cell(91,100)stevedore(1)^tan(5).*(D(1)-0.9)十x⑸;end-endmean(stevedore.*60)令行奮口324.4588图9sie^doreFn1000x1double112.0222.2G34.6743.2S55,6365.97710.M80.00■&Q.W107.4611S.48124.2713€.831斗5.0115£.92166.09172.94IE2.3«195.583D4.87215.62225.80237,4S245.H255.40餌3.19?7图8利用histfit作图可得出如下装卸时间的分布，如图9。＞）sums匕点甲亡日皿己/门点葢七闪口匕!;图10>>(suma-ffait4suiristevedorc,/1000ans-0.0139图11图124.采用T2产生的1000个船舶到港时间间隔数据，T3产生的1000艘船舶装卸服务时间,假设只有1台桥吊，对1000艘船舶的装卸排队服务过程进行仿真：统计桥吊利用率、每艘船舶平均在港总时间、每艘船舶平均等待时间解：吊桥利用率=总服务时间/总时间平均在港时间=（总服务时间+总装卸时间）/船舶数量平均等待时间=总等待时间/船舶数量其中总的服务时间就是第一艘船开始到最后一艘船的结束，如图10，11，12，13。yiire_fun.m[Untitled/Lm[Urrtitfed旷[UnBBSt瓷[llrrtitfed!5畫Untitled6*V|+n^ztwork=0;2arTival=0;3sumstevedore=sum〔石tEvedore)45umawait=0;5Efor1=1:10006ifnext.work>arrival”sumawa1t~simawa11+(nest.wark-arrIvalJ8nextwork=arr1val+(nestwork-arrivalj9else0tiextwork=arrival+stevedore(i)Lend2arriva1=arrivaL+i...
