离散事件模拟(银行业务模拟。实现算法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\n",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分钟\n",Qu,Khjg,Blsj);
	printf("客户总数:%d, 所有客户共耗时:%ld分钟,平均每人耗时:%d分钟,",CustomerNum,TotalTime,TotalTime/CustomerNum);
	printf("最后一个客户离开的时间:%d分\n",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\n下一客户到达时刻=%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\n",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\n",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 个队列,程序运行结果如下:


转载于:https://www.cnblogs.com/KongkOngL/p/3945957.html

这是这次期末课程设计的一道题目<br>刚写出来<br>离散事件模拟<br>①银行业务模拟<br>【问题描述】客户业务分为两种,第一种是申请从银行得到一笔资金,即取款或者借款。第二种是向银行投入一笔资金,即存款或者还款。银行有两个服务窗口,相应的有两个队列。客户到达银行后先排第一个队。处理每个客户业务时,如果属于第一种,且申请额超出银行现存资金总额而得不到满足,则立即排入第二个队等候,直至满足时才离开银行;否则业务处理完后立即离开银行。每接待完一个第二种业务的客户,则顺序检查和处理(如果可能)第二个队列中的客户,对能满足的申请者予以满足,不能满足的者重新排到第二个队列的队尾。<br>注意:在此检查过程中,一旦银行的资金总额少于或者等于刚才第一个队列中最后一个客户(第二种业务)被接待之前的数额,或者本次已将第二个队列检查或处理了一遍,就停止检查(因为此时已不可能还有能满足者)转而继续接待第一个队列客户。任何时刻都只开一个窗口。假设检查不需要时间。营业时间结束时所有客户立即离开银行。<br>【基本要求】利用动态存储结构实现模拟。<br>【测试数据】一天营业开始时银行拥有的款额为10000(元),营业时间为600(分钟)。其他模拟参量自定,注意测定两种极端的情况:一是两个到达事件之间的间隔时间很短,而客户的交易时间很长,另一个恰好相反,设置两个到达事件的间隔时间很长,而客户的交易时间很短。<br>【实现提示】事件有两类:到达银行的和离开银行。初始时银行现存资金总额为total。开始营业后的第一个事件是客户到达,营业时间从0到closetime。到达事件发生时随机地设置此客户的交易时间和距下一到达事件之间的时间间隔。每一个客户要办理的款额也是随机确定的,用负值和正值分别表示第一类和第二类业务。变量total,closetime以及上述两个随机量的上下界均交互地从终端读入,作为模拟参数。<br>两个队列和一个事件表均要用动态存储结构实现。注意弄清应该在什么条件下设置离开事件,以及第二个队列用怎么样的存储结构实现时可以获得较高的效率。注意:事件表是按时间顺序有序的。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值