matlab排队系统仿真,MM1排队系统仿真matlab实验报告
MM1排队系统仿真matlab实验报告
M/M/1 排队系统实验报告一、实验目的本次实验要求实现 M/M/1 单窗口无限排队系统的系统仿真,利用事件调度法实现离散事件系统仿真,并统计平均队列长度以及平均等待时间等值,以与理论分析结果进行对比。二、实验原理根据排队论的知识我们知道,排队系统的分类是根据该系统中的顾客到达模式、服务模式、服务员数量以及服务规则等因素决定的。1、 顾客到达模式设到达过程是一个参数为 的 Poisson 过程,则长度为 的时间内到达 个tk呼叫的概率 服从 Poisson 分布,即 , ,其中 >0etkktp!)(,210为一常数,表示了平均到达率或 Poisson 呼叫流的强度。2、 服务模式设每个呼叫的持续时间为 ,服从参数为 的负指数分布,即其分布函数i为 {}1,0tPXte3、 服务规则先进先服务的规则(FIFO)4、 理论分析结果在该 M/M/1 系统中,设 ,则稳态时的平均等待队长为 ,顾客1Q的平均等待时间为 。T三、实验内容M/M/1 排队系统:实现了当顾客到达分布服从负指数分布,系统服务时间也服从负指数分布,单服务台系统,单队排队,按 FIFO(先入先出队列)方式服务。四、采用的语言MatLab 语言源代码:clear;clc;%M/M/1排队系统仿真SimTotal=( 请输入仿真顾客总数SimTotal= ); %仿真顾客总数;Lambda=0.4; %到达率Lambda;Mu=0.9; %服务率Mu;t_Arrive=zeros(1,SimTotal); t_Leave=zeros(1,SimTotal);ArriveNum=zeros(1,SimTotal);LeaveNum=zeros(1,SimTotal);Interval_Arrive=-log(rand(1,SimTotal))/Lambda;%到达时间间隔Interval_Serve=-log(rand(1,SimTotal))/Mu;%服务时间t_Arrive(1)=Interval_Arrive(1);%顾客到达时间ArriveNum(1)=1;for i=2:SimTotalt_Arrive(i)=t_Arrive(i-1)+Interval_Arrive(i);ArriveNum(i)=i;endt_Leave(1)=t_Arrive(1)+Interval_Serve(1);%顾客离开时间LeaveNum(1)=1;for i=2:SimTotalif t_Leave(i-1)=2QueLength(i)=CusNum(i)-1;elseQueLength(i)=0;endendQueLength_avg=sum([0 QueLength].*[Time_interval 0] )/Timepoint(end);%系统平均等待队长%仿真图figure(1);set(1, position ,[0,0,1000,700]);subplot(2,2,1);title( 各顾客到达时间和离去时间 );stairs([0 ArriveNum],[0 t_Arrive], b );hold on;stairs([0 LeaveNum],[0 t_Leave], y );legend( 到达时间 , 离去时间 );hold off;subplot(2,2,2);stairs(Timepoint,CusNum, b )title( 系统等待队长分布 );xlabel( 时间 );ylabel( 队长 );subplot(2,2,3);title( 各顾客在系统中的排队时间和等待时间 );stairs([0 ArriveNum],[0 t_Queue], b );hold on;stairs([0 LeaveNum],[0 t_Wait], y );hold off;legend( 排队时间 , 等待时间 );%仿真值与理论值比较disp([ 理论平均等待时间t_Wait_avg= ,num2str(1/(Mu-Lambda))]);disp([ 理论平均排队时间t_Wait_avg= ,num2str(Lambda/(Mu*(Mu-Lambda)))]);disp([ 理论系统中平均顾客数= ,num2str(Lambda/(Mu-Lambda))]);disp([ 理论系统中平均等待队长= ,num2str(Lambda*Lambda/(Mu*(Mu-Lambda)))]);disp([ 仿真平均等待时间t_Wait_avg= ,num2str(t_Wait_avg)])disp([ 仿真平均排队时间t_Queue_avg= ,num2str(t_Queue_avg)])disp([ 仿真系统中平均顾客数= ,num2str(CusNum_avg)]);disp([ 仿真系统中平均等待队长= ,num2str(QueLength_avg)]);五、数据结构1.仿真设计算法(主要函数)利用负指数分布与泊松过程的关系,产生符合泊松过程的顾客流,产生符合负指数分布的随机变量作为每个顾客的服务时间:Interval_Arrive=-log(rand(1,SimTotal))/Lambda;%到达时间间隔,结果与调用exprnd(1/Lambda,m)函数产生的结果相同Interval_Serve=-log(rand(1,SimTotal))/Mu;%服务时间间隔t_Arrive(1)=Interval_Arrive(1);%顾客到达时间时间计算t_Wait=t_Leave-t_Arrive; %各顾客在系统中的等待时间t_Queue=t_Wait-Interval_Serve; %各顾客在系统中的排队时间由事件来触发仿真时钟的不断推进。每发生一次事件,记录下两次事件间隔的时间以及在该时间段内排队的人数:Timepoint=[t_Arrive,t_Leave]; %系统中顾客数变化CusNum=zeros(size(Timepoint));CusNum_avg=sum(CusNum_fromStart.*[Time_interval 0] )/Timepoint(end); %系统中平均顾客数计算QueLength_avg=sum([0 QueLength].*[Time_interval 0] )/Timepoint(end); %系统平均等待队长2.算法的流程图开始计算第 1 个顾客的离开时间: i-2输入仿真人数计算第 i 个顾客的等待时间、离开时间、标示位: i+1标志位置 0:i=i+1系统是否接纳第i 个顾客?仿真时间是否越界?结束输出结果六、仿真结果分析顾客的平均等待时间与顾客的平均等待队长,计算其方差如下:仿 真 顾 客 总 数=100000 1 2 3 4 5 平 均 值 方 差平 均 等 待 时 间 2.023 1.9971 1.9945 1.9961 2.0043 2.003 0.000556360平 均 排 队 时 间 0.91147 0.8865 0.88293 0.88404 0.89495 0.89198 0.000
matlab排队系统仿真,MM1排队系统仿真matlab实验报告相关推荐
- matlab 滑动平均滤波,滑动平均滤波器实验报告
滑动平均滤波器实验报告 所属分类:matlab例程 开发工具:matlab 文件大小:798KB 下载次数:19 上传日期:2018-01-27 16:12:36 上 传 者:玉玲珑 说明: 给出一 ...
- matlab摆线等时性程序,摆的等时性实验报告.doc
摆的等时性实验报告 摆的等时性实验报告 篇一:科学实验报告单 科学实验报告单1 实验名称 物体的沉浮 实验目的 观察物体的沉浮 实验材料 水槽.水.塑料.小刀.泡沫.橡皮.萝卜.曲别针等各种材料 实验 ...
- led的伏安特性曲线 matlab实现_小灯泡伏安特性曲线实验报告
2016 全新精品资料 - 全新公文范文 - 全程指导写作 – 独家原创 1 / 5 小灯泡伏安特性曲线实验报告 篇一: <描绘小灯泡的伏安特性曲线>的实验报告 一.实验目的 描绘小灯泡的 ...
- matlab计算鲈鱼质量,鲈鱼数学建模实验报告.pdf
数学建模实验报告 姓名:胡斌 学号: 一. 摘 要 题目提供了哈德逊河鲈鱼的年龄分组. 成年鱼的年龄. 允许捕捞 鱼的年龄段.各年龄段的鱼的存活率以及各组成年雌性鱼每年能产雌 性后代的个数.题目初始数 ...
- 数值分析方程求根实验matlab,基于matlab的数值分析( 非线性方程求根)上机实验报告1...
数值分析试验报告 非线性方程求根 二分法解方程 1. 题目:用二分法求方程0133 =--x x 的所有根 2. 方法:二分法 3. 程序 function x=erfenfa(a,b) if (a* ...
- cgns matlab,MATLABSimulink系统建模与仿真实验报告(详细解析)(word文档良心出品).docx...
MATLAB Simuli nk 电力系统建模与仿真 ******实验报告 ****** 姓名: 专业?电气工程及其自动化 班级: ******************* 学号: ********** ...
- matlab 高斯迭代法求解,高斯迭代法matlab算例
Matlab 线性方程组的迭代解法 Gauss-Seidel 迭代法 Matlab 线性方程组的迭代解法 Gauss-Seidel 迭代法实验报告 1.熟悉 Gauss-Seidel 迭代法,并编写 ...
- mm1排队系统仿真matlab实验报告,MM1排队系统仿真matlab实验报告.doc
M/M/1排队系统实验报告 一.实验目的 本次实验要求实现M/M/1单窗口无限排队系统的系统仿真,利用事件调度法实现离散事件系统仿真,并统计平均队列长度以及平均等待时间等值,以与理论分析结果进行对比. ...
- mm1排队系统仿真matlab实验报告,matlab仿真实验报告_mm1排队系统仿真matlab实验报告资料...
M/M/1排队系统实验报告 一.实验目的 本次实验要求实现M/M/1单窗口无限排队系统的系统仿真,利用事件调度法实现离散事件系统仿真,并统计平均队列长度以及平均等待时间等值,以与理论分析结果进行对比. ...
- mm1 matlab,(高质)MM1排队系统仿真matlab实验报告.doc
M/M/1排队系统实验报告 一.实验目的 本次实验要求实现M/M/1单窗口无限排队系统的系统仿真,利用事件调度法实现离散事件系统仿真,并统计平均队列长度以及平均等待时间等值,以与理论分析结果进行对比. ...
最新文章
- Geany 权限问题:Error opening file ... : permission denied.,原因及解决办法。
- leveldb 学习记录(四)Log文件
- AJAX Control Toolkit——DropShadow(阴影和圆角框)
- 使用复合设计模式扩展持久化的CURD,Select能力
- CSMA/CD协议(一目了然,看过都说好)
- nginx 后端获取真实ip
- 2020跨境电商独立站将喷发式增长?
- paip.一千 常用汉字 高频汉字 覆盖率90%
- 计算机科学导论:第六章 计算机网络和因特网
- 领导力有哪些方面?什么是领导力?
- java外卖项目介绍_JavaWeb网上订餐系统项目
- h5 app跳转客服咨询 临时会话 (没有开通在线咨询、无法会话)
- java工程师js_Ajax, 励志做一名Java全干工程师
- 敬业签云便签提示“基础连接已关闭,发送时发生错误!”如何解决
- 举个栗子!Tableau 技巧(93):用 热图标记 地图或散点图
- word-breakbreak-all和word-wrapbreak-word的区别
- 始于2015年11月23日,终于2018年08月08日
- (十二) 灰度共生矩阵
- Android tint颜色渲染
- RLException: Invalid roslaunch XML syntax: mismatched tag: line 3, column 2报错