港口系统仿真过程11.doc1.docVIP

港口系统仿真 课程实验报告 姓名 杨学云 班级 交运专132班 学号 201120611042 院系 交通运输学院 指导老师 罗莉华 港口系统模拟仿真过程 交运专132班 杨学云201120611042 1） 用线性同余法产生1000个[0,1]独立均匀分布的随机数，要求按照以下规则尝试两组参数，产生两组1000个随机数 a) m为2的幂，即m为2的b次方 （比如b取20）并且 c不等于0 ，c是相对于m的素数（两者最大公约数为1），且（k=0,1…） 具体步骤如下： a取1+4K为4949，c取一个奇数11579，m为2的32次方为4294967296，根据线性同余法公式Xi+1=(aXi+c)modm，计算得出X1为5415887，以及相应的X1000,再将得到的X1-X1000每一个都去除以m得到1000个0-1的随机数，平均值为0.487。具体操作见（excel蓝色部分）。 b）m为2的幂，即m为2的b次方,并且 c等于0 ，种子X0为奇数，且乘子a 满足 a=3+8k 或者a=5+8k（k=0,1…） 具体步骤如下： a取3+8K为83，c取一个奇数0，m为2的35次方根据线性同余法公式Xi+1=(aXi+c)modm，计算得出X1以及相应的X1000,再将得到的X1-X1000每一个都去除以m得到1000个0-1的随机数，平均值为0.4924。具体操作见（excel黄色部分）。 （2）假设在某港口装卸服务系统中，通过统计，有以下数据： 船舶到港过程：服从每天平均3.2艘船的泊松到达过程 （A）以第(1)问产生的第一组随机数为基础，按照上述分布特点产生1000艘船舶的到港时间间隔（以min为单位），并计算你所产生的这1000个到达时间间隔的平均值。 具体步骤如下： （1）由于船舶到港过程，服从每天平均3.2艘船的泊松到达过程，则船舶到达时间间隔服从指数分布，则平均到达时间间隔为1440分钟除以3.2分钟为450分钟，根据反变换技术法，计算所要求的随机变量X的累积分布函数(cdf) F(x)，在X的范围内令F(X)=R(R服从[0,1]上的均匀分布)，求解F(X)=R，以得到X ，计算出平均值为434.19903，具体操作见（excel紫色部分）。 （B）以100min为划分间隔，画出产生到达时间间隔在各区间的频率分布图。具体步骤如下： 求出0-100，100-200,200-300，300-400,400-500.，，，1900-2000各区间船舶到达艘数，求出各区间的频率，得到频率分布图，具体操作见（excel紫色部分）。 （3）统计了200艘船舶的装卸服务时间，如下表 区间（小时 3到5 5到7 7到9 9到11 11到13 频数 35 60 45 30 15 （A）根据该数据拟合出装卸服务时间这个随机变量的累积分布函数. 具体步骤如下：总频数为180，的累积分布频率表如下 区间（小时 3到5 5到7 7到9 9到11 11到13 频数 35 60 45 30 15 185 频率 0.189189 0.324324 0.243243 0.162162 0.081081 累积分布点 1 3 5 7 9 11 13 累积频率 0 0 0.189 0.513 0.756 0.918 1 制作累积分布图, 具体操作见（excel橘黄色部分）。 （B）以第(1)问产生的第二组随机数为基础，按照上述统计规律模拟产生1000艘船舶的装卸服务时间（单位min），并计算你所产生的1000艘船舶装卸服务时间的平均值。具体步骤如下： （1）根据表格搜集的数据，估计船舶装卸时间这个随机变量的累积分布函数；（2）利用反变换技术法，用（0~1）均匀分布的随机数反变换得到装卸时间。 第一步求斜率 求斜率 斜率 方程为 0=y0.189189189 ⊿x/⊿y 10.58201 X=3+10.582*(Y-0) 0.189=y0.0.5135 ⊿x/⊿y 6.17284 X=5+6.17*(Y-0.189) 0.51351=y0.0.75675 ⊿x/⊿y 8.230453 X=7+8.230*(Y-0.513) 0.75675=y0.918 ⊿x/⊿y 12.34568 X=9+12.34*(Y-0.757) 0.918=y1 ⊿x/⊿y 24.39024 X=11+24.39*(Y-0.918) 第二步根据第一题第二问得到的1000个随机数得到1000个船舶装卸时间，平均值为439.3821分钟，具体操作见（excel橘黄色部分）。 （C）以100min为区间间隔，画出产生数据在各区间的频率分布图。 根据得到的1000个随机数，按不同区间得出频率分布如下表： 根据上表制作下图：