zoukankan      html  css  js  c++  java
  • 离散事件模拟(银行业务模拟。实现算法3.6、3.7的程序)

    // func3-3.cpp、algo3-12.cpp和algo3-13.cpp用到的函数及变量等
    #include"c1.h"
    typedef struct // 定义ElemType为结构体类型
    {
    	int OccurTime; // 事件发生时刻
    	int NType; // 事件类型,Qu表示到达事件,0至Qu-1表示Qu个窗口的离开事件
    }Event,ElemType; // 事件类型,有序链表LinkList的数据元素类型(见图3.40)
    #include"c2-5.h" // 从实际应用角度出发重新定义的线性链表结构
    typedef LinkList EventList; // 事件链表指针类型,定义为有序链表
    #include"bo2-6.cpp" // 基于c2-5.h存储结构的基本操作
    typedef struct
    {
    	int ArrivalTime; // 到达时刻
    	int Duration; // 办理事务所需时间
    }QElemType; // 定义队列的数据元素类型QElemType为结构体类型(见图3.41)
    #include"c3-2.h" // 链队列存储结构
    #include"bo3-2.cpp" // 链队列基本操作
    // 程序中用到的主要变量(全局)
    EventList ev; // 事件表头指针
    Event en,et; // 事件,临时变量
    //FILE *fp; // 文件型指针,用于指向b.txt或d.txt文件
    long int TotalTime=0; // 累计客户逗留时间(初值为0)
    int CloseTime,CustomerNum=0; // 银行营业时间(单位是分),客户数(初值为0)
    int cmp(Event a,Event b)
    { // 依事件a的发生时刻<、=或>事件b的发生时刻分别返回-1、0或1
    	if(a.OccurTime==b.OccurTime)
    		return 0;
    	else
    		return (a.OccurTime-b.OccurTime)/abs(a.OccurTime-b.OccurTime);
    }
    void Random(int &d,int &i)
    { // 生成两个随机数
    	d=rand()%Blsj+1; // 1到Blsj之间的随机数(办理业务的时间)
    	i=rand()%(Khjg+1); // 0到Khjg之间的随机数(客户到达的时间间隔)
    }
    void OpenForDay();
    void CustomerArrived();
    void CustomerDeparture();
    void Bank_Simulation()
    { // 银行业务模拟函数
    	Link p;
    	OpenForDay(); // 初始化事件表ev且插入第1个到达事件,初始化队列
    	while(!ListEmpty(ev)) // 事件表ev不空
    	{
    		DelFirst(ev,ev.head,p); // 删除事件表ev的第1个结点,并由p返回其指针,在bo2-6.cpp中
    		// if(p->data.OccurTime<50) // 输出前50分钟内发生的事件到文件d.txt中
    		// fprintf(fp,"p->data.OccurTime=%3d p->data.NType=%d
    ",p->data.OccurTime,p->data.NType);
    		en.OccurTime=GetCurElem(p).OccurTime;
    		// GetCurElem()在bo2-6.cpp中,返回p->data(ElemType类型)
    		en.NType=GetCurElem(p).NType;
    		if(en.NType==Qu) // 到达事件
    			CustomerArrived(); // 处理客户到达事件
    		else // 由某窗口离开的事件
    			CustomerDeparture(); // 处理客户离开事件
    	} // 计算并输出平均逗留时间
    	printf("窗口数=%d 两相邻到达的客户的时间间隔=0~%d分钟每个客户办理业务的时间=1~%d分钟
    ",Qu,Khjg,Blsj);
    	printf("客户总数:%d, 所有客户共耗时:%ld分钟,平均每人耗时:%d分钟,",CustomerNum,TotalTime,TotalTime/CustomerNum);
    	printf("最后一个客户离开的时间:%d分
    ",en.OccurTime);
    }
    

    // algo3-12.cpp 银行业务模拟。实现算法3.6、3.7的程序
    #define Qu 4 // 客户队列数
    #define Khjg 5 // 两相邻到达的客户的时间间隔最大值
    #define Blsj 30 // 每个客户办理业务的时间最大值
    #include"func3-3.cpp" // 包含algo3-12.cpp和algo3-13.cpp共同用到的函数和变量等
    LinkQueue q[Qu]; // Qu个客户队列
    QElemType customer; // 客户记录,临时变量
    //FILE *fq; // 文件型指针,用于指向a.txt文件
    void OpenForDay()
    { // 初始化事件链表ev且插入第1个到达事件,初始化Qu个队列
    	int i;
    	InitList(ev); // 初始化事件链表ev为空(见图3.42)
    	en.OccurTime=0; // 设定第1位客户到达时间为0
    	// (银行一开门,就有客户来)
    	//fprintf(fq,"首位客户到达时刻=%3d,",en.OccurTime);
    	en.NType=Qu; // 到达
    	OrderInsert(ev,en,cmp); // 将第1个到达事件en有序插入事件表ev中,在bo2-6.cpp中(见图3.43)
    	for(i=0;i<Qu;++i) // 初始化Qu个队列(见图3.44)
    		InitQueue(q[i]);
    }
    int Minimum(LinkQueue Q[])
    { // 返回最短队列的序号,若有并列值,返回队列序号最小的
    	int l[Qu];
    	int i,k=0;
    	for(i=0;i<Qu;i++)
    		l[i]=QueueLength(Q[i]);
    	for(i=1;i<Qu;i++)
    		if(l[i]<l[0])
    		{
    			l[0]=l[i];
    			k=i;
    		}
    		return k;
    }
    void CustomerArrived()
    { // 处理客户到达事件en(en.NType=Qu)
    	QElemType f;
    	int durtime,intertime,i;
    	++CustomerNum; // 客户数加1
    	Random(durtime,intertime);//生成当前客户办理业务的时间和下一个客户到达的时间间隔2个随机数
    	et.OccurTime=en.OccurTime+intertime; // 下一客户et到达时刻等于当前客户en的到达
    	// 时间加时间间隔
    	et.NType=Qu; // 下一客户到达事件
    	i=Minimum(q); // 求长度最短队列的序号,等长为最小的序号(到达事件将入该队)
    	//if(CustomerNum<=20) // 输出前20个客户到达信息到文件a.txt中
    	// fprintf(fq,"办理业务的时间=%2d,所排队列=%d
    下一客户到达时刻=%3d,",durtime,i,et.OccurTime);
    	if(et.OccurTime<CloseTime) // 下一客户到达时银行尚未关门
    		OrderInsert(ev,et,cmp); // 按升序将下一客户到达事件et插入事件表ev中,在bo2-6.cpp中
    	f.ArrivalTime=en.OccurTime; // 将当前客户到达事件en赋给队列元素f
    	f.Duration=durtime;
    	EnQueue(q[i],f); // 将f入队到第i队列(i=0~Qu-1)
    	if(QueueLength(q[i])==1) // 该元素为队头元素
    	{
    		et.OccurTime=en.OccurTime+durtime; // 设定一个离开事件et
    		et.NType=i;
    		OrderInsert(ev,et,cmp); // 将此离开事件et按升序插入事件表ev中
    	}
    }
    void CustomerDeparture()
    { // 处理客户离开事件en(en.NType<Qu)
    	int i;
    	i=en.NType; // 确定离开事件en发生的队列序号i
    	DeQueue(q[i],customer); // 删除第i队列的排头客户
    	TotalTime+=en.OccurTime-customer.ArrivalTime;
    	// 客户逗留时间=离开事件en的发生时刻-该客户的到达时间
    	if(!QueueEmpty(q[i]))
    	{ // 删除第i队列的排头客户后,第i队列仍不空
    		GetHead(q[i],customer); // 将第i队列新的排头客户赋给customer
    		et.OccurTime=en.OccurTime+customer.Duration;
    		// 设定离开事件et,新排头的离开时间等于原排头的离开时间加新排头办理业务的时间
    		et.NType=i; // 第i个队列的离开事件
    		OrderInsert(ev,et,cmp); // 将此离开事件et按升序插入事件表ev中
    	}
    }
    void main()
    {
    	//fq=fopen("a.txt","w"); // 打开a.txt文件,用于写入客户到达信息
    	//fp=fopen("b.txt","w"); // 打开b.txt文件,用于写入有序事件表的历史记录
    	printf("请输入银行营业时间长度(单位:分): ");
    	scanf("%d",&CloseTime);
    	//srand(time(0));
    	// 设置随机数种子,以使每次运行程序产生的随机数不同(time(0)是长整型数,与调用时间有关)
    	Bank_Simulation();
    	//fclose(fq); // 关闭a.txt
    	//fclose(fp); // 关闭b.txt
    }
    代码的运行结果如下:

    /*
    请输入银行营业时间长度(单位:分): 480
    窗口数=4 两相邻到达的客户的时间间隔=0~5分钟每个客户办理业务的时间=1~30分钟
    客户总数:184, 所有客户共耗时:22684分钟,平均每人耗时:123分钟,最后一个客户离开的时
    间:753分
    Press any key to continue
    */


    文件a.txt 的内容是每位客户的到达时刻和他办理业务所需的时间(二者根据随机函数
    依次生成)以及他所排的队列等信息。排队的原则是排在人数最少的队列,如果有不止一
    队的人数都同为最少,则排在序号最小的队列。由文件a.txt 的第3 条记录可知,第8 分
    钟来了一个办理业务需要17 分钟的客户,这时除0 队有1 人外(1 队在第4 分钟来的人,
    已在第7 分钟离去),1~3 队均无人,客户排在了序号最小的1 队。
    文件b.txt 的内容是事件表ev 的历史记录。ev 就象一个安装在门口的监视器,按照时
    间顺序,记录客户的到达和离去。对于离去的客户,还要记录是由哪个窗口离去。文件
    b.txt 的第3 条记录记下了第7 分钟在1 队(1 号窗口)发生了一个离去事件。这和由文件

    a.txt 推算的结果(见图345)相吻合。

    主程序main()的第5 行语句的作用是根据程序运行当前时间的不同,产生不同的随机
    函数。否则,随机数总是一样的。在调试程序时,希望每次产生相同的随机数,以便分析
    程序运行结果。而在实际应用时,则希望每次产生不同的随机数。
    在运行algo3-12.cpp 时不会产生a.txt 和b.txt 文件,因为有关产生这两个文件的语句
    已标为注释。在调试程序时,可加一些输出语句帮助分析,最后提交的程序应只输出需要
    的结果。实际上,迷宫求解问题中,图311(到达终点时栈S 的内容)的数据也是根据增
    加的输出语句得出的。
    以上是银行开4 个服务窗口,两相邻客户到达的时间间隔为0~5 分钟,每个客户办
    理业务所需的时间为1~30 分钟的模拟结果。模拟结果显示,客户等候耗时太多。最后一
    个客户离开时,已14 个半小时了。如果开6 个窗口(修改程序algo3-12.cpp 的第1 行,定
    义Qu 为6,6 个队列),程序运行结果如下


    平均每位客户等候耗时大大缩短。
    目前银行大多使用排队机,即1 个队列,多个窗口。algo3-13.cpp 是根据这种情况编
    制的程序。其中OpenForDay()、CustomerArrived()和CustomerDeparture()等3 个函数与
    algo3-12.cpp 中的相应函数是同名函数,但内容不同。而Bank_Simulation(),甚至main()
    函数(除了打开的文件名不同之外)都是相同的。

    // algo3-13.cpp 使用排队机的银行业务模拟
    #define Qu 4 // 窗口数
    #define Khjg 5 // 两相邻到达的客户的时间间隔最大值
    #define Blsj 30 // 每个客户办理业务的时间最大值
    #include"func3-3.cpp" // 包含algo3-12.cpp和algo3-13.cpp共同用到的函数和变量等
    LinkQueue q; // 排队机队列q
    QElemType customer[Qu]; // Qu个客户队列元素, 存放正在窗口办理业务的客户的信息
    //FILE *fq; // 文件型指针,用于指向c.txt文件
    //int j=0; // 计数器,产生c.txt文件用到
    Boolean chk[Qu]; // 窗口状态,1为闲,0为忙
    void OpenForDay()
    { // 初始化事件链表ev且插入第1个到达事件,初始化排队机q,初始化Qu个窗口为1(空闲)
    	int i;
    	InitList(ev); // 初始化事件链表ev为空(见图3.42)
    	en.OccurTime=0; // 设定第1位客户到达时间为0(银行一开门,就有客户来)
    	en.NType=Qu; // 到达
    	OrderInsert(ev,en,cmp); // 将第1个到达事件en有序插入事件表ev中,在bo2-6.cpp中(见图3.43)
    	InitQueue(q); // 初始化排队机队列q(见图3.47)
    	for(i=0;i<Qu;i++)
    		chk[i]=1; // 初始化Qu个窗口为1(空闲) (见图3.48)
    }
    int ChuangKou()
    { // 返回空闲窗口的序号(0~Qu-1),若有多个窗口空闲,返回序号最小的。若无空闲窗口,返回Qu
    	int i;
    	for(i=0;i<Qu;i++)
    		if(chk[i])
    			return i;
    		return i;
    }
    void CustomerArrived()
    { // 处理客户到达事件en(en.NType=Qu),与algo3-12.cpp不同
    	QElemType f;
    	int durtime,intertime,i;
    	++CustomerNum; // 客户数加1
    	Random(durtime,intertime);//生成当前客户办理业务的时间和下一个客户到达的时间间隔2个随机数
    	et.OccurTime=en.OccurTime+intertime; // 下一客户et到达时刻=当前客户en的到达时间+时间间隔
    	et.NType=Qu; // 下一客户到达事件
    	if(et.OccurTime<CloseTime) // 下一客户到达时银行尚未关门
    		OrderInsert(ev,et,cmp); // 按升序将下一客户到达事件et插入事件表ev中,在bo2-6.cpp中
    	f.ArrivalTime=en.OccurTime; // 将当前客户到达事件en赋给队列元素f
    	f.Duration=durtime;
    	EnQueue(q,f); // 将当前客户f入队到排队机
    	i=ChuangKou(); // 求空闲窗口的序号
    	if(i<Qu) // 有空闲窗口
    	{
    		DeQueue(q,customer[i]); // 删去排队机的排头客户(也就是刚入队的f由排队机到i号窗口)
    		// if(j++<20) // 输出前20个客户到达信息及处理业务窗口信息到文件c.txt中
    		// fprintf(fq,"客户到达时刻=%3d,办理业务的时间=%2d,所在窗口=%d
    ",customer[i].ArrivalTime,customer[i].Duration,i);
    		et.OccurTime=en.OccurTime+customer[i].Duration; // 设定一个i号窗口的离开事件et
    		et.NType=i; // 第i号窗口的离开事件
    		OrderInsert(ev,et,cmp); // 将此离开事件et按升序插入事件表ev中
    		chk[i]=0; // i号窗口状态变忙
    	}
    }
    void CustomerDeparture()
    { // 处理客户离开事件en(en.NType<Qu),与algo3-12.cpp不同
    	int i;
    	i=en.NType; // 确定离开事件en发生的窗口序号i
    	chk[i]=1; // i号窗口状态变闲
    	TotalTime+=en.OccurTime-customer[i].ArrivalTime;
    	// 客户逗留时间=离开事件en的发生时刻-该客户的到达时间
    	if(!QueueEmpty(q))
    	{ // 第i窗口的客户离开后,排队机仍不空
    		DeQueue(q,customer[i]);
    		// 删去排队机的排头客户并将其赋给customer[i](排队机的排头客户去i窗口办理业务)
    		// if(j++<20) // 输出前20个客户到达信息及处理业务窗口信息到文件c.txt中
    		// fprintf(fq,"客户到达时刻=%3d,办理业务的时间=%2d,所在窗口=%d
    ",customer[i].ArrivalTime,customer[i].Duration,i);
    		chk[i]=0; // i号窗口状态变忙
    		et.OccurTime=en.OccurTime+customer[i].Duration;
    		// 设定customer[i]的离开事件et,客户的离开时间等于原客户的
    		// 离开时间加当前客户办理业务的时间
    		et.NType=i; // 第i号窗口的离开事件
    		OrderInsert(ev,et,cmp); // 将此离开事件et按升序插入事件表ev中
    	}
    }
    void main()
    {
    	//fq=fopen("c.txt","w"); // 打开c.txt文件,用于写入客户到达信息
    	//fp=fopen("d.txt","w"); // 打开d.txt文件,用于写入有序事件表的历史记录
    	printf("请输入银行营业时间长度(单位:分): ");
    	scanf("%d",&CloseTime);
    	//srand(time(0));
    	// 设置随机数种子,以使每次运行程序产生的随机数不同(time(0)是长整型数,与调用时间有关)
    	Bank_Simulation();
    	//fclose(fq); // 关闭c.txt
    	//fclose(fp); // 关闭d.txt
    }

    代码的运行结果如下:

    /*
    请输入银行营业时间长度(单位:分): 480
    窗口数=4 两相邻到达的客户的时间间隔=0~5分钟每个客户办理业务的时间=1~30分钟
    客户总数:184, 所有客户共耗时:22323分钟,平均每人耗时:121分钟,最后一个客户离开的时
    间:724分
    Press any key to continue
    
    */



    和algo3-12.cpp 的运行结果相比,在营业时间、营业窗口数和客户总数相同条件下,
    algo3-13.cpp 的平均每人耗时和最后一个客户离开的时间比algo3-12.cpp 要略小。因为没
    有调用srand()函数,两程序产生的随机数是一样的(通过比较文件a.txt 和c.txt 的内容可
    看出)。不同的是,客户所排队列或窗口不完全相同。algo3-12.cpp 在关门后(不再有客户
    进入)可能出现某队已空,而其他队还有若干人的情况。这就延长了最后一个客户离开的
    时间,也增大了平均每人耗时。

    修改程序algo3-13.cpp 的第1 行,定义Qu 为6,6 个队列,程序运行结果如下:


  • 相关阅读:
    Winform中设置ZedGraph曲线图的字体样式是避免出现边框
    50本精品前端开发书籍免费下载
    一切“归零”的开始
    小公司大企业
    给“精点们”的一封信
    一个程序员的自我救赎
    一个微服务框架的情节
    《深入理解Java虚拟机》虚拟机类加载机制
    一个微服务框架的故事
    一个程序员的自我剖析
  • 原文地址:https://www.cnblogs.com/KongkOngL/p/3945957.html
Copyright © 2011-2022 走看看