排队论经典程序MM1代码Word文件下载.docx
《排队论经典程序MM1代码Word文件下载.docx》由会员分享，可在线阅读，更多相关《排队论经典程序MM1代码Word文件下载.docx（14页珍藏版）》请在冰豆网上搜索。
area_3_in_q+=time_since_last_event;
3 //店内至少有一个顾客的概率
if（server_status==BUSY）//服务台忙,则店内至少有一个顾客
abv_1+=time_since_last_event;
4 //在店内顾客的平均数
if（server_status==BUSY）//服务台忙,总的顾客数为排队顾客数加一
area_num_in_h+=（num_in_q+1）*time_since_last_event;
5 total_of_server+=time_next_event[2]-sim_time;
//总的服务时间加一个服务时间为新的服务总时间
delay=sim_time-time_arrival[1];
//排队时间=当前时间-这个人来的时间
total_of_delays+=delay;
6 //离开时总的消耗时间大于15,必须在店内消耗15分钟以上的顾客数加一
if（（delay+time_next_event[2]-sim_time）>
15）
abv_15++;
//到达时总的服务时间大于15,必须在店内消耗15分钟以上的顾客数加一
if（（time_next_event[2]-sim_time）>
abv_15++;
程序代码：
/*Externaldefinitionsforsingle-serverqueueingsystem.*/
#include<
stdio.h>
math.h>
/*#include"
lcgrand.h"
Headerfileforrandom-numbergenerator.*/
#defineQ_LIMIT100/*Limitonqueuelength.队伍最长100人*/
#defineBUSY1/*Mnemonicsforserver'
sbeingbusy忙碌状态*/
#defineIDLE0/*andidle.空闲状态*/
intnext_event_type,//下一个事件类型
num_custs_delayed,//已模拟的顾客数
num_delays_required,//模拟的顾客数
num_events,//事件数
num_in_q,//队列中的顾客数
server_status;
//服务状态
floatarea_num_in_q,//有顾客的时间
area_server_status,//总的服务时间
mean_interarrival,//平均顾客到达时间间隔
mean_service,//平均服务时间
sim_time,//模拟时间
time_arrival[Q_LIMIT+1],//到来的时间
time_last_event,//上一个事件的时间
time_next_event[3],//下一个事件的时间
total_of_delays;
//总的排队时间
////////////////////////////////////////////////////////////////////////////////////
//添加的变量
floatabv_15,//15分钟以上的顾客数量
total_of_server,//所有顾客的总的服务时间
area_3_in_q,//有3个顾客的时间
abv_1,//至少有一个顾客的时间
area_num_in_h;
//顾客总数
FILE*infile,*outfile;
/*Thefollowing3declarationsareforuseoftherandom-numbergenerator
lcgrandandtheassociatedfunctionslcgrandstandlcgrandgtforseed
management.Thisfile（namedlcgrand.h）shouldbeincludedinanyprogram
usingthesefunctionsbyexecuting
#include"
beforereferencingthefunctions.*/
floatlcgrand（intstream）;
voidlcgrandst（longzset,intstream）;
longlcgrandgt（intstream）;
voidinitialize（void）;
voidtiming（void）;
voidarrive（void）;
voiddepart（void）;
voidreport（void）;
voidupdate_time_avg_stats（void）;
floatexpon（floatmean）;
main（）/*Mainfunction.*/
{
/*Openinputandoutputfiles.*/
infile=fopen（"
mm1.in"
"
r"
）;
outfile=fopen（"
mm1.out"
w"
/*Specifythenumberofeventsforthetimingfunction.*/
num_events=2;
//两种事件
/*Readinputparameters.*/
fscanf（infile,"
%f%f%d"
&
mean_interarrival,&
mean_service,
&
num_delays_required）;
/*Writereportheadingandinputparameters.输出*/
fprintf（outfile,"
Single-serverqueueingsystem\n\n"
Meaninterarrivaltime%11.3fminutes\n\n"
mean_interarrival）;
Meanservicetime%16.3fminutes\n\n"
mean_service）;
Numberofcustomers%14d\n\n"
num_delays_required）;
/*Initializethesimulation.初始化仿真*/
initialize（）;
//初始化
/*Runthesimulationwhilemoredelaysarestillneeded.没服务完，仿真继续*/
while（num_custs_delayed<
num_delays_required）{//当已服务顾客数小于1000时
/*Determinethenextevent.确定下一事件*/
timing（）;
/*Updatetime-averagestatisticalaccumulators.时间记录更新*/
update_time_avg_stats（）;
/*Invoketheappropriateeventfunction.根据事件的不同，调用不同的函数*/
switch（next_event_type）{
case1:
arrive（）;
//到达
break;
case2:
depart（）;
//离开
}
/*Invokethereportgeneratorandendthesimulation.*/
report（）;
fclose（infile）;
fclose（outfile）;
return0;
}
voidinitialize（void）/*Initializationfunction.*/
/*Initializethesimulationclock.仿真时间置为0*/
sim_time=0.0;
/*Initializethestatevariables.最开始状态初始化*/
server_status=IDLE;
//服务空闲
num_in_q=0;
//队伍里无人排队
time_last_event=0.0;
//上一个事件的时间，最开始肯定是0开始
/*Initializethestatisticalcounters.*/
num_custs_delayed=0;
//已经服务的人数
total_of_delays=0.0;
//总的排队时间
area_num_in_q=0.0;
//有顾客的时间
area_server_status=0.0;
//总的服务时间
/////////////////////////////////////////////////////////////////////////////////
//添加的变量的初始化
area_3_in_q=0.0;
//有3个顾客的时间
abv_1=0.0;
area_num_in_h=0.0;
//顾客的总数
total_of_server=0.0;
//所有顾客的所有的服务的时间
abv_15=0.0;
//消耗15分钟以上的顾客数
/*Initializeeventlist.初始化事件列表Sincenocustomersarepresent,thedeparture
（servicecompletion）eventiseliminatedfromconsideration.无顾客存在和离开*/
time_next_event[1]=sim_time+expon（mean_interarrival）;
//下一事件是来的时间
time_next_event[2]=1.0e+30;
//下一事件是离开的时间
voidtiming（void）/*Timingfunction.*/
inti;
floatmin_time_next_event=1.0e+29;
//像指针一样的对于当前服务的人来说下一个事件的时间
next_event_type=0;
/*Determinetheeventtypeofthenexteventtooccur.接下来将要发生的事件的类型*/
for（i=1;
i<
=num_events;
++i）
if（time_next_event[i]<
min_time_next_event）{//下一事件是来的时间跟离开时间比较
min_time_next_event=time_next_event[i];
next_event_type=i;
/*Checktoseewhethertheeventlistisempty.*/
if（next_event_type==0）{
/*Theeventlistisempty,sostopthesimulation.无事件，停止仿真过程*/
\nEventlistemptyattime%f"
sim_time）;
exit
（1）;
/*Theeventlistisnotempty,soadvancethesimulationclock.有事件，进行仿真过程*/
sim_time=min_time_next_event;
//仿真的时间就是当前事件的时间
voidarrive（void）/*Arrivaleventfunction.*/
floatdelay;
/*Schedulenextarrival.计划下一次的到来*/
/*Checktoseewhetherserverisbusy.检测是否在服务状态*/
if（server_status==BUSY）{
/*Serverisbusy,soincrementnumberofcustomersinqueue.在服务则排队多一人*/
++num_in_q;
/*Checktoseewhetheranoverflowconditionexists.检测人数是否超出*/
if（num_in_q>
Q_LIMIT）{
/*Thequeuehasoverflowed,sostopthesimulation.*/
\nOverflowofthearraytime_arrivalat"
time%f"
exit
（2）;
/*Thereisstillroominthequeue,sostorethetimeofarrivalofthe
arrivingcustomeratthe（new）endoftime_arrival.队列中仍有空间时，记录新到达的时间*/
time_arrival[num_in_q]=sim_time;
//在这个时间的时候有这么多的排队人数，用于计算3顾客的问题
else{//服务空闲的状况
/*Serverisidle,soarrivingcustomerhasadelayofzero.（The
followingtwostatementsareforprogramclarityanddonotaffect
theresultsofthesimulation.）*/
delay=0.0;
/*Incrementthenumberofcustomersdelayed,andmakeserverbusy.*/
++num_custs_delayed;
//已经模拟的顾客数加1
server_status=BUSY;
//人到来，服务开始
/*Scheduleadeparture（servicecompletion）.服务完成*/
time_next_event[2]=sim_time+expon（mean_service）;
//这个人离开的时间为现在时间+服务时间
//总的服务时间加上当前服务时间,更新总的服务时间
total_of_server+=time_next_event[2]-sim_time;
//总的服务时间大于15,必须在店内消耗15分钟以上的顾客数加一
15）//如果这个人服务时间超过15分钟则耗费15分钟人数加1
voiddepart（void）/*Departureeventfunction.讨论离开事件*/
/*Checktoseewhetherthequeueisempty.检测队列是否为空*/
if（num_in_q==0）{
/*Thequeueisemptysomaketheserveridleandeliminatethe
departure（servicecompletion）eventfromconsideration.队列空，服务空闲*/
//离开的时间无限大（无人离开）
else{
/*Thequeueisnonempty,sodecrementthenumberofcustomersin
queue.有人离开，队列人数减少*/
--num_in_q;
/*Computethedelayofthecustomerwhoisbeginningserviceandupdate
thetotaldelayaccumulator.*/
/*Incrementthenumberofcustomersdelayed,andscheduledeparture.已经服务人数+1*/
//服务人数加1
//当前接受服务的人的离开时间
//总的消耗时间大于15,必须在店内消耗15分钟以上的顾客数加一
/*Moveeachcustomerinqueue（ifany）uponeplace.有人离开，队列前移*/
=num_in_q;
++i）
time_arrival[i]=time_arrival[i+1];
//人的到达时间也前移
voidreport（void）/*Reportgeneratorfunction.*/
/*Computeandwriteestimatesofdesiredmeasuresofperformance.*/
\n\nAveragedelayinqueue%11.3fminutes\n\n"
total_of_delays/num_custs_delayed）;
Averagenumberinqueue%10.3f\n\n"
area_num_in_q/sim_time）;
Serverutilization%15.3f\n\n"
area_server_status/sim_time）;
Timesimulationended%12.3fminutes"
printf（"
统计量:
\n"
////////////////////////////////////////////////////////////////////////
//总时间减去服务台忙的时间除以总时间,得到服务台空闲的概率
%24.3f\n"
%20.3f\n"
%16.3f\n"
④在
升级会员 