排队随机服务系统的模拟与R.docVIP

统计计算课程设计 题 目 摘 要 1 绪论 2方法介绍 一两个服务台排队系统，在这个系统中，顾客按照指数分布，每个到达的顾客首先接受第一个服务台的服务，其服务完成之后再接受第二个服务员的服务。这样的系统称为串联排队系统。如果服务台1是空闲的，来到的顾客就接受服务员1?的服务；否则，来到的顾客就进入到等候接受服务员1服务的队伍中排队。类似地，当顾客接受完服务员1的服务后，如果服务员2是空闲的，此顾客就接受服务员2?的服务；否则，此顾客进入到等候接受服务员2服务的队伍中排队。在接受完服务员2的服务后此顾客就离开系统。假设每个顾客接受服务员i的服务时间有分布函数Gi,i=1，2.假设我们的目的是研究顾客在两个服务台的平均逗留时间和排队中的顾客平均数。我们采用模拟的方法实现这个目的。在模拟中我们将用到下述变量。 队长（Ls) 队长是指在系统中的顾客的平均数 顾客的平均等待时间与平均逗留时间（Ws) 顾客的平均等待时间是指从顾客进入系统的时刻起直到开始接受服务止的平均时间，平均逗留时间是指系统中的平均等待时间与平均服务时间之和，平均等待时间与平均服务时间是顾客最关心的数量指标。 系统的忙期与闲期 从顾客到达空闲的系统，服务立即开始，到系统再次变为空闲，这段时间是系统连续繁忙的时间，我们称为系统的忙期，他反应了系统中服务机构的工作强度，是衡量服务机构利用效率的指标。 时间变量t? 系统状态变量 ss.(n1,n2):有n1个顾客在服务员1（包括正在接受服务的）前排队，有n2个顾客在服务员2（包括正在接受服务的）前排队。 计数变量?：到时刻t时来到的顾客数。???????????? ：到时刻t时离开的顾客数。? 输出变量 (n): 顾客n?到达的时刻， (n): 顾客n?(n): 顾客n?到达的时刻， D（n)：顾客n离开的时刻， 事件清单 其中是下一个顾客到达的时间，是正在接受服务员i，i=1，2服务完成的时间。如果在服务员i处没有顾客接受服务，那么=∞，i=假设是具有分布，i=1，2的随机变量，C为系统中顾客数。 3 程序及实现结果 3.1 R程序 js-function(lambda, mu, T){ k-0; wt-0; wn-0; ws-0; tp-0; nA-0; n-0; t-0 r-runif(1); tA--1/lambda*log(r); tD-Inf repeat{ k-k+1; wt[k]-t; wn[k]-n if (tA T){ ws[k]-min(tA, tD)-t if (tA tD){ t-tA; n-n+1; nA-nA+1 r-runif(1); tA-t-1/lambda*log(r) if (n==1){ r-runif(1); tD-t-1/mu*log(r) } }else{ t-tD; n-n-1 if (n==0){ tD-Inf }else{ r-runif(1); tD-t-1/mu*log(r) } } }else{ ws[k]-if(tD==Inf) 0 else tD-t if (n0){ t-tD; n-n-1 if (n0){ r-runif(1); tD-t-1/mu*log(r) } }else tp-1 } if (tp==1) break } data.frame(Ls=sum(ws*wn)/t, Ws=sum(ws*wn)/nA, Pwait=sum(ws[wn=1])/t) } js(0.6,0.8,600) js(0.6,0.8,1000) 3.2 运行结果 5总结 当T=600时，系统中顾客的平均数，即队长数为2,.6441，平均逗留时间为4.3357，平均等待时间为0,.7611。 当T=1000时，系统中顾客的平均数，即队长数为3.7637，平均逗留时间为6.4381，平均等待时间为0。7901。 其理论值为队长数为5，平均逗留时间为5，平均等待时间为0,.75。 从模拟600分钟和模拟1000分钟的排队服从系统的运行情况来看还是很接近的