单服务台排队系统仿真报告
单服务台排队系统仿真报告
一、模型准备
1、 顾客到达特性
在该系统中，顾客的到达规模（成批到达还是单个到达）是单个到达，假设顾客到达率Ai 服从均值为 的指数分布，即
2、 顾客服务时间
顾客服务时间为Si ，服从指数分布，假设均值为 ，即
二、 仿真模型设计
1、 元素定义（Define ）
本系统的元素定义如表1所示。
2、 元素可视化设置（Display ）
本系统中各个元素的显示特征定义设置如图2所示：
min 5=A βA
s A
e A
f ββ/)(-=
)
0(≥A min 4=s βS
A S
e S
f ββ/)(-=
)
0(≥S
图2 各元素的显示特征
（1）Part元素可视化设置
在元素选择窗口选择customer元素，鼠标右键点击Display，跳出Display 对话框（图3），设置它的Text（图4）、Icon（图5）。
图3 Display对话框
图4 Display Text对话框
图5 Display Icon对话框
（2）Buffer元素可视化设置
在元素选择窗口选择paidui元素，鼠标右键点击Display，跳出Display对话框（图3），设置它的Text、Icon、Rectangle（图6）。
图6 Display Rectangle对话框
（3）Machine元素可视化设置
在元素选择窗口选择Fuwuyuan元素，鼠标右键点击Display，跳出Display 对话框（图3），设置它的Text、Icon、Part Queue（图7）。
图7 Display Part Queue对话框
（4）Variable元素可视化设置
在元素选择窗口选择Jifen0元素，鼠标右键点击Display，跳出Display对话框（图3），设置它的Text 、Value（图8）。
图8 Display Value对话框
（5）Timeseries元素可视化设置
在元素选择窗口选择duichang元素，鼠标右键点击Display，跳出Display 对话框（图3），设置它的Text、Timeseries（图9）。
图9 Timeseries对话框
3、元素细节（Detail）设计
（1）对Part元素Guke细节设计
Type：Active
Input to Model. Inter Arrival：-5 * LN（Random（1））Lot Size：1
Input to Model. To…：Push to paidui
Actions on Create：Icon = 58
Detail Part 对话框如图10所示。
图10 Detail Part对话框
（2）对Part元素Jifen细节设计
Type：Active
Input to Model. Inter Arrival：1.0
Lot Size：1
Input to Model. To…：Push to Ship
Actions on Create：Jifen0 = Jifen0 + Nparts（paidui）
（3）对Buffer元素Paidui细节设计
Capacity：100
（4）对Machine元素Fuwuyuan细节设计
Type：Single
Input. From：Pull from Paidui
Duration. Cycle Time：-4 * LN（Random（2））
Output. To…：Push to Ship
Detail machine 对话框如图11所示。
图11 Detail Machine 对话框
（5）对Timeserises元素Duichang细节设计
●Recording：5.0
●在Plot Expressions下第一个Plot中将“Undefined”改为Nparts（paidui）●在Plot Expressions下第二个Plot中将“Undefined”改为Nparts
（Fuwuyuan）
Detail Timeseries对话框如图12所示。
图12 Detail Timeseries对话框
三、模型运行和数据报告
为了模仿邮局的工作时间（约8小时一天），模型仿真时钟取480个时间单位。
选中系统中所有元素，点击Reports菜单下的Statistics选项（图13），得到如图14所示数据统计报告。
图13 Statistics查看统计报告
图14 数据统计报告
以上为各元素的统计数据报告。
通过报告，可以获得该排队系统此次仿真的平均队长、最大队长、最小队长和平均队长以及平均每位顾客的等待时间、平均服务时间等信息。
通过数据报告发现，不同顾客的服务时间和顾客的到达特性，对应的仿真结果有所不同。顾客的到达特性以及顾客的服务时间都影响着排队系统的最大队
长、最小队长和平均队长以及平均每位顾客的等待时间。
实验单服务台单队列排队系统仿真
实验2排队系统仿真 一、学习目的 1．了解仿真的特点 2．学习如何建构模型 3．熟悉eM-Plant基本的对象和操作 4．掌握排队系统的特点与仿真的实现方法 二、问题描述 该银行服务窗口为每个到达的顾客服务的时间是随机的，表2.4是顾客服务时间纪录的统计结果 表2.4 每个顾客服务时间的概率分布 服务时间（min）概率密度累计概率 1 0.1 0.1 2 0.2 0.3 3 0.3 0.6 4 0.2 5 0.85 5 0.1 0.95 6 0.05 1.0 对于上述这样一个单服务待排队系统，仿真分析30天，分析该系统中顾客的到
达、等待和被服务情况，以及银行工作人员的服务和空闲情况。 三、系统建模 3.1 仿真目标 通过对银行排队系统的仿真，研究银行系统的服务水平和改善银行服务水平的方法，为银行提高顾客满意度，优化顾客服务流程服务。 3.2．系统建模 3.2.1 系统调研 1. 系统结构: 银行服务大厅的布局, 涉及的服务设备 2. 系统的工艺参数: 到达-取号-等待-服务-离开 3. 系统的动态参数: 顾客的到达时间间隔, 工作人员的服务时间 4. 逻辑参数: 排队规则, 先到先服务 5. 系统的状态参数: 排队队列是否为空, 如果不为空队长是多少, 服务台是否为空 6. 系统的输入输出变量:输入变量确定其分布和特征值,顾客的到达时间间隔的概率分布表和每个顾客被服务时间的概率分布. 输出变量根据仿真目标设定. 包括队列的平均队长、最大队长、仿真结束时队长、总服务人员、每个顾客的平均服务时间、顾客平均排队等待服务时间、业务员利用率等。 3.2.2系统假设 1．取号机前无排队，取号时间为0 2．顾客排队符合先进先出的排队规则 3．一个服务台一次只能对一个顾客服务 4．所有顾客只有一种单一服务 5．仿真时间为1个工作日（8小时） 6．等候区的长度为无限长 3.2.3系统建模 系统模型： 3.2.4 仿真模型 1．实体：银行系统中的实体是人（主动体）
控制系统仿真课程设计报告.
控制系统仿真课程设计 （2011级） 题目控制系统仿真课程设计学院自动化 专业自动化 班级 学号 学生姓名 指导教师王永忠/刘伟峰 完成日期2014年6月
控制系统仿真课程设计一 ———交流异步电机动态仿真 一 设计目的 1．了解交流异步电机的原理，组成及各主要单元部件的原理。 2. 设计交流异步电机动态结构系统； 3．掌握交流异步电机调速系统的调试步骤，方法及参数的整定。 二 设计及Matlab 仿真过程 异步电机工作在额定电压和额定频率下，仿真异步电机在空载启动和加载过程中的转速和电流变化过程。仿真电动机参数如下： 1.85, 2.658,0.2941,0.2898,0.2838s r s r m R R L H L H L H =Ω=Ω===， 20.1284Nm s ,2,380,50Hz p N N J n U V f =?===，此外，中间需要计算的参数如下： 21m s r L L L σ=-，r r r L T R =，22 2 s r r m t r R L R L R L +=，10N m TL =?。αβ坐标系状态方程： 其中，状态变量： 输入变量： 电磁转矩： 2p m p s r s L r d ()d n L n i i T t JL J βααωψψβ=--r m r r s r r d 1d L i t T T ααβαψψωψ=--+r m r r s r r d 1d L i t T T ββαβψψωψ=-++22s s r r m m m s r r s s 2r r r r d d i R L R L L L L i u t L T L L ααβαα σψωψ+=+-+22 s s r r m m m s r r s s 2 r r r r d d i R L R L L L L i u t L T L L ββαββ σψωψ+=--+[ ] T r r s s X i i αβαβωψψ=[ ] T s s L U u u T αβ=()p m e s s s s r n L T i i L βααβ ψψ=-
排队系统仿真matlab实验报告
M/M/1排队系统实验报告 一、实验目的 本次实验要求实现M/M/1单窗口无限排队系统的系统仿真，利用事件调度法实现离散事件系统仿真，并统计平均队列长度以及平均等待时间等值，以与理论分析结果进行对比。 二、实验原理 根据排队论的知识我们知道，排队系统的分类是根据该系统中的顾客到达模式、服务模式、服务员数量以及服务规则等因素决定的。 1、 顾客到达模式 设到达过程是一个参数为λ的Poisson 过程，则长度为t 的时间内到达k 个呼叫的概 率 服从Poisson 分布，即e t k k k t t p λλ-=!)()(，?????????=,2,1,0k ，其中λ>0为一常数，表示了 平均到达率或Poisson 呼叫流的强度。 2、 服务模式 设每个呼叫的持续时间为i τ，服从参数为μ的负指数分布，即其分布函数为{}1,0t P X t e t μ-<=-≥ 3、 服务规则 先进先服务的规则（FIFO ） 4、 理论分析结果 在该M/M/1系统中，设 λρμ=，则稳态时的平均等待队长为1Q ρλρ=-，顾客的平均等待时间为T ρ μλ=-。 三、实验内容 M/M/1排队系统：实现了当顾客到达分布服从负指数分布，系统服务时间也服从负指数分布，单服务台系统，单队排队，按FIFO （先入先出队列）方式服务。 四、采用的语言 MatLab 语言 源代码： clear; clc;
%M/M/1排队系统仿真 SimTotal=input('请输入仿真顾客总数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:SimTotal t_Arrive(i)=t_Arrive(i-1)+Interval_Arrive(i); ArriveNum(i)=i; end t_Leave(1)=t_Arrive(1)+Interval_Serve(1);%顾客离开时间LeaveNum(1)=1; for i=2:SimTotal if t_Leave(i-1) 实验2---单服务台排队系统的仿真 姓名：学号： 一、目标任务 ①模拟路由器缓存区M|M|1|m实验。 ②设定：λ＝8/s，μ＝10/s，ρ＝0.8，m=10。 ③模拟系统106s，求系统报文的丢失率及报文在路由器中停留时间的均值。 ④模拟100次，图展示每次的模拟结果，并与理论值0.0184比较。 二、编程语言 Matlab 三、关键代码 lamda = 8; %报文到达强度 u = 10; %路由器处理强度 m = 10; %路由器缓冲区长度 T = 1000000; %模拟时间 a = []; %模拟运行时丢失率的运行结果 mean_a = 0; %模拟运行时丢失率的平均运行结果 ref_value = 0.0184; %丢失率理论值大小 b = []; %模拟运行时报文在路由器中的停留时间 mean_b = 0; %模拟运行时报文在路由器中停留时间的均值 %模拟运行一百次 for i=1:100 time = 0; %绝对时钟 t = 0; %路由器的下一空闲时刻 N = 0; %到达报文数 NI = 0; %丢失报文数 q = 0; %队长 stay_time = 0; %报文在路由器中的停留时间 %按指数分布产生随机到达时间和服务时间 while 1 CRTime = exprnd(1/lamda); %按指数分布产生下一报文的到达随机时间间隔 time = CRTime + time; %下一个报文到达的时间 if time > T break; end N = N + 1; q = q + 1; while q > 0 & t < time q = q - 1; ServeTime = exprnd(1/u);%按指数分布产生报文的随机服务时间 if q == 0 t = time + ServeTime; else t = t + ServeTime; end stay_time = stay_time + ServeTime * (q + 1); end if q == m + 1 %如果超过缓冲区长，则丢失报文数加1，队长减1 NI = NI + 1; q = q - 1; end end a = [a, NI/N]; b = [b, stay_time/(N-NI)]; end %计算结果 mean_a = mean(a); mean_b = mean(b); %绘图 x = 1:100; plot(x, a, x, mean_a); %绘制模拟运行时丢包率变化图以及均值线 scatter(x, a, '.'); %绘制模拟运行时丢包率变化散点图 scatter(x, b, '.'); %绘制模拟运行时平均停留时间变化散点图 fprintf('平均丢包率%6.5f\n', mean_a); % 打印平均丢包率 fprintf('平均停留时间%6.5f\n', mean_b); % 打印平均停留时间 四、实验结果与分析 课程设计报告 题目PID控制器应用 课程名称控制系统仿真院部名称机电工程学院专业 班级 学生姓名 学号 课程设计地点 课程设计学时 指导教师 金陵科技学院教务处制成绩 一、课程设计应达到的目的 应用所学的自动控制基本知识与工程设计方法，结合生产实际，确定系统的性能指标与实现方案，进行控制系统的初步设计。 应用计算机仿真技术，通过在MATLAB软件上建立控制系统的数学模型，对控制系统进行性能仿真研究，掌握系统参数对系统性能的影响。 二、课程设计题目及要求 1．单回路控制系统的设计及仿真。 2．串级控制系统的设计及仿真。 3．反馈前馈控制系统的设计及仿真。 4．采用Smith 补偿器克服纯滞后的控制系统的设计及仿真。 三、课程设计的内容与步骤 (1)．单回路控制系统的设计及仿真。 （a）已知被控对象传函W(s) = 1 / (s2 +20s + 1)。 （b）画出单回路控制系统的方框图。 （c）用MatLab的Simulink画出该系统。 （d）选PID调节器的参数使系统的控制性能较好，并画出相应的单位阶约响应曲线。注明所用PID调节器公式。PID调节器公式Wc（s）=50(5s+1)/(3s+1) 给定值为单位阶跃响应幅值为3。 有积分作用单回路控制系统 无积分作用单回路控制系统 大比例作用单回路控制系统 （e）修改调节器的参数，观察系统的稳定性或单位阶约响应曲线，理解控制器参数对系统的稳定性及控制性能的影响？ 答：由上图分别可以看出无积分作用和大比例积分作用下的系数响应曲线，这两个PID调节的响应曲线均不如前面的理想。增大比例系数将加快系统的响 matlab 单服务台排队系统实验报告 一、实验目的 本次实验要求实现M/M/1单窗口无限排队系统的系统仿真，利用事件调度法实现离散事件系统仿真，并统计平均队列长度以及平均等待时间等值，以与理论分析结果进行对比。 二、实验原理 根据排队论的知识我们知道，排队系统的分类是根据该系统中的顾客到达模式、服务模式、服务员数量以及服务规则等因素决定的。 1、 顾客到达模式 设到达过程是一个参数为λ的Poisson 过程，则长度为t 的时间内到达k 个呼 叫的概率 服从Poisson 分布，即 e t k k k t t p λλ-= !)()(，?????????=,2,1,0k ，其中λ>0为一 常数，表示了平均到达率或Poisson 呼叫流的强度。 2、 服务模式 设每个呼叫的持续时间为i τ，服从参数为μ的负指数分布，即其分布函数为 {}1,0t P X t e t μ-<=-≥ 3、 服务规则 先进先服务的规则（FIFO ） 4、 理论分析结果 在该M/M/1系统中，设λρμ= ，则稳态时的平均等待队长为1Q ρλ ρ= -，顾客 的平均等待时间为 T ρμλ= -。 三、实验内容 M/M/1排队系统：实现了当顾客到达分布服从负指数分布，系统服务时间也服 从负指数分布，单服务台系统，单队排队，按FIFO 方式服务。 四、采用的语言 MatLab 语言 源代码： clear; clc; %M/M/1排队系统仿真 SimTotal=input('请输入仿真顾客总数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:SimTotal t_Arrive(i)=t_Arrive(i-1)+Interval_Arrive(i); ArriveNum(i)=i; end t_Leave(1)=t_Arrive(1)+Interval_Serve(1);%顾客离开时间LeaveNum(1)=1; for i=2:SimTotal if t_Leave(i-1) 物流系统建模与仿真 09级自动化学院物流工程1班 20085435 詹乐思 20095277 安静 20095278 陈红玲 20095289 陈均剑 20095290 翟瑞 20095291 胡旺 单服务台排队系统仿真研究报告 ——选重庆大学A区门口中国银行分行某一服务窗口为单服务台排队系统研究 对象 一、系统基本背景 社会的进步越来越快，人们的生活节奏也随之越来越快。在科技的发展，新技术的普及下, 我国的银行业以计算机和信息技术、互联网技术为前提, 通过大量资金和科技的投入, 不断地开发出新产品和新业务。另外有网上银行、支付宝等新业务的出现, 大大提高了工作效率。然而现代的金融服务并不是都可以靠刷卡来解决, 许多技术还不完善, 这些新技术也并不适合所有顾客群,去银行办理业务的顾客仍然经常性地出现排队现象。顾客等待时间过长, 造成顾客满意度下降, 矛盾较为突出, 因此本报告试利用单服务台排队论的方法, 定性定量地对具有排队等候现象的银行服务系统进行统计调查与分析研究,希望能帮助改进银行工作效率, 优化系统的运营。 本报告研究对象为中国银行重庆大学处分行某一服务窗口，数据取自银行内唯一非现金业务柜台。研究对象的选取虽然不是最典型的，但是综合考虑了研究地域范围和小组成员作业时间有限，另有其他方案由于各种原因无法进行，故选择离学校较近的有代表性的中国银行中的服务窗口作为最终方案。 中国银行简介：中国银行是中国历史最为悠久的银行之一，在大家对银行的概念中有着一定地位。中国银行主营传统商业银行业务，包括公司金融业务、个人金融业务和金融市场业务。公司业务以信贷产品为基础，致力于为客户提供个性化、创新的金融服务和融资、财务解决方案。个人金融业务主要针对个人客户的金融需求，提供包括储蓄存款、消费信贷和银行卡在内的服务。作为中国金融行业的百年品牌，中国银行在稳健经营的同时，积极进取，不断创新，创造了国内银行业的许多第一，在国际结算、外汇资金和贸易融资等领域得到业界和客户的广泛认可和赞誉。 二、系统描述 该银行工作时间为上午8：30至下午16：30（周一至周日），另周末不办理对公业务，属于每天8小时工作制。系统调查对象为银行内唯一非现金业务柜台，可知到达的顾客中，需要办理非现金业务的顾客在正常现金业务柜台忙碌的情况下可以选择该服务台。在队列中，等待服务的顾客和服务台构成了一个排队系统。由于银行前台出纳员逐个接待顾客，当顾客较多的时候就会出现排队等待的现象。其中，顾客的到达是随机的，每两个先后到达的顾客的到达间隔时间是不确定的。 本排队系统用顾客的数目、到达模式、服务模式、系统容量和排队规则来描述。 为探求此排队系统的规律, 首先需确定顾客流在一定时间内到达的概率分布函数。抵达本银行服务窗口的顾客流量大体上服从Poisson 分布, 顾客流抵达银行便按先后顺序排队, 进入单服务窗口，即排队论中的M/M/1系统。所谓M/M/1排队系统是指这样的一种排队模型: 顾客的到达为Poisson 流, 银行对 控制系统仿真课程设计 （2010级） 题目控制系统仿真课程设计学院自动化 专业自动化 班级 学号 学生姓名 指导教师王永忠/刘伟峰 完成日期2013年7月 控制系统仿真课程设计（一） ——锅炉汽包水位三冲量控制系统仿真1.1 设计目的 本课程设计的目的是通过对锅炉水位控制系统的Matlab仿真，掌握过程控制系统设计及仿真的一般方法，深入了解反馈控制、前馈-反馈控制、前馈-串级控制系统的性能及优缺点，实验分析控制系统参数与系统调节性能之间的关系，掌握过程控制系统参数整定的方法。 1.2 设计原理 锅炉汽包水位控制的操作变量是给水流量，目的是使汽包水位维持在给定的范围内。汽包液位过高会影响汽水分离效果，使蒸汽带水过多，若用此蒸汽推动汽轮机，会使汽轮机的喷嘴、叶片结垢，严重时可能使汽轮机发生水冲击而损坏叶片。汽包液位过低，水循环就会被破坏，引起水冷壁管的破裂，严重时会造成干锅，甚至爆炸。 常见的锅炉汽水系统如图1-1所示，锅炉汽包水位受汽包中储水量及水位下汽包容积的影响，而水位下汽包容积与蒸汽负荷、蒸汽压力、炉膛热负荷等有关。影响水位变化的因素主要是锅炉蒸发量（蒸汽流量）和给水流量，锅炉汽包水位控制就是通过调节给水量，使得汽包水位在蒸汽负荷及给水流量变化的情况下能够达到稳定状态。 图1-1 锅炉汽水系统图 在给水流量及蒸汽负荷发生变化时，锅炉汽包水位会发生相应的变化，其分别对应的传递函数如下所示： （1）汽包水位在给水流量作用下的动态特性 汽包和给水可以看做单容无自衡对象，当给水增加时，一方面会使得汽包水位升高，另一方面由于给水温度比汽包内饱和水的温度低，又会使得汽包中气泡减少，导致水位降低，两方面的因素结合，在加上给水系统中省煤器等设备带来延迟，使得汽包水位的变化具有一定的滞后。因此，汽包水位在给水流量作用下，近似于一个积分环节和惯性环节相串联的无自衡系统，系统特性可以表示为 ()111()()(1)K H S G S W S s T s ==+ （1.1） （2）汽包水位在蒸汽流量扰动下的动态特性 在给水流量及炉膛热负荷不变的情况下，当蒸汽流量突然增加时，瞬间会导致汽包压力的降低，使得汽包内水的沸腾突然加剧，水中气泡迅速增加，将整个水位抬高；而当蒸汽流量突然减小时，汽包内压力会瞬间增加，使得水面下汽包的容积变小，出现水位先下降后上升的现象，上述现象称为“虚假水位”。虚假水位在大中型中高压锅炉中比较显著，会严重影响锅炉的安全运行。“虚假水位”现象属于反向特性，变化速度很快，变化幅值与蒸汽量扰动大小成正比，也与压力变化速度成正比，系统特性可以表示为 222()()()1f K K H s G s D s T s s ==-+ （1.2） 常用的锅炉水位控制方法有：单冲量控制、双冲量控制及三冲量控制。单冲量方法仅是根据汽包水位来控制进水量，显然无法克服“虚假水位”的影响。而双冲量是将蒸汽流量作为前馈量用于汽包水位的调节，构成前馈-反馈符合控制系统，可以克服“虚假水位”影响。但双冲量控制系统要求调节阀具有好的线性特性，并且不能迅速消除给水压力等扰动的影响。为此，可将给水流量信号引入，构成三冲量调节系统，如图1-2所示。图中LC 表示水位控制器（主回路），FC 表示给水流量控制器（副回路），二者构成一个串级调节系统，在实现锅炉水位控制的同时，可以快速消除给水系统扰动影响；而蒸汽流量作为前馈量用于消除“虚假水位”的影响。 排队系统仿真实验报告 1 实验题目 考虑一个排队系统，顾客到达后若无空闲服务台则等待，否则接受服务后离开，仿真目的是统计服务台前的平均队长和服务台的忙闲程度，要求能利用离散事件仿真的方法对该系统进行仿真，系统的规则和参数自定。 2 实验原理 离散事件驱动的模型只在一些离散点上由事件改变其状态，故离散事件模型是由事件驱动的。驱动某一模型的所有事件按其发生的时间先后构成一个序列，通常要求按时间先后顺序处理事件，而不能颠倒。离散事件系统仿真的关键是按时间顺序确定这一序列。除初始事件，事件序列中的事件不能在仿真前事先确定，而是在仿真进行中产生，离散事件系统仿真中一般采用事先策划事件的方式，即在仿真系统处理任何事件之前该事件必须已被策划。确定事件发生的时间的方式有三类：直接方式、间接方式、混合方式。从事件、活动、进程三个层次来组织事件即构成了处理离散事件模型的三种典型处理方法：事件调度法、活动扫描法、进程交互法。本次实验针对排队系统采用事件调度法，仿真策略如图1所示： 图1事件调度法仿真排队系统流程图 3 仿真程序 #include 物流系统建模与仿真 单服务台排队系统仿真研究报告 ——选大学A区门口中国银行分行某一服务窗口为单服务台排队系统研究对象一、系统基本背景 社会的进步越来越快，人们的生活节奏也随之越来越快。在科技的发展，新技术的普及下, 我国的银行业以计算机和信息技术、互联网技术为前提, 通过大量资金和科技的投入, 不断地开发出新产品和新业务。另外有网上银行、支付宝等新业务的出现, 大大提高了工作效率。然而现代的金融服务并不是都可以靠刷卡来解决, 许多技术还不完善, 这些新技术也并不适合所有顾客群,去银行办理业务的顾客仍然经常性地出现排队现象。顾客等待时间过长, 造成顾客满意度下降, 矛盾较为突出, 因此本报告试利用单服务台排队论的方法, 定性定量地对具有排队等候现象的银行服务系统进行统计调查与分析研究,希望能帮助改进银行工作效率, 优化系统的运营。 本报告研究对象为中国银行大学处分行某一服务窗口，数据取自银行唯一非现金业务柜台。研究对象的选取虽然不是最典型的，但是综合考虑了研究地域围和小组成员作业时间有限，另有其他方案由于各种原因无法进行，故选择离学校 较近的有代表性的中国银行中的服务窗口作为最终方案。 中国银行简介：中国银行是中国历史最为悠久的银行之一，在大家对银行的概念中有着一定地位。中国银行主营传统商业银行业务，包括公司金融业务、个人金融业务和金融市场业务。公司业务以信贷产品为基础，致力于为客户提供个性化、创新的金融服务和融资、财务解决方案。个人金融业务主要针对个人客户的金融需求，提供包括储蓄存款、消费信贷和银行卡在的服务。作为中国金融行业的百年品牌，中国银行在稳健经营的同时，积极进取，不断创新，创造了国银行业的许多第一，在国际结算、外汇资金和贸易融资等领域得到业界和客户的广泛认可和赞誉。 二、系统描述 该银行工作时间为上午8：30至下午16：30（周一至周日），另周末不办理对公业务，属于每天8小时工作制。系统调查对象为银行唯一非现金业务柜台，可知到达的顾客中，需要办理非现金业务的顾客在正常现金业务柜台忙碌的情况下可以选择该服务台。在队列中，等待服务的顾客和服务台构成了一个排队系统。由于银行前台出纳员逐个接待顾客，当顾客较多的时候就会出现排队等待的现象。其中，顾客的到达是随机的，每两个先后到达的顾客的到达间隔时间是不确定的。 本排队系统用顾客的数目、到达模式、服务模式、系统容量和排队规则来描述。 为探求此排队系统的规律, 首先需确定顾客流在一定时间到达的概率分布 单服务台排队系统仿真报告 一、模型准备 1、 顾客到达特性 在该系统中，顾客的到达规模（成批到达还是单个到达）是单个到达，假设顾客到达率Ai 服从均值为 的指数分布，即 2、 顾客服务时间 顾客服务时间为Si ，服从指数分布，假设均值为 ，即 二、 仿真模型设计 1、 元素定义（Define ） 本系统的元素定义如表1所示。 2、 元素可视化设置（Display ） 本系统中各个元素的显示特征定义设置如图2所示： min 5=A βA s A e A f ββ/)(-= ) 0(≥A min 4=s βS A S e S f ββ/)(-= ) 0(≥S 图2 各元素的显示特征 （1）Part元素可视化设置 在元素选择窗口选择customer元素，鼠标右键点击Display，跳出Display 对话框（图3），设置它的Text（图4）、Icon（图5）。 图3 Display对话框 图4 Display Text对话框 图5 Display Icon对话框 （2）Buffer元素可视化设置 在元素选择窗口选择paidui元素，鼠标右键点击Display，跳出Display对话框（图3），设置它的Text、Icon、Rectangle（图6）。 图6 Display Rectangle对话框 （3）Machine元素可视化设置 在元素选择窗口选择Fuwuyuan元素，鼠标右键点击Display，跳出Display 对话框（图3），设置它的Text、Icon、Part Queue（图7）。 图7 Display Part Queue对话框 （4）Variable元素可视化设置 在元素选择窗口选择Jifen0元素，鼠标右键点击Display，跳出Display对话框（图3），设置它的Text 、Value（图8）。 图8 Display Value对话框 控制系统数字仿真课程设计 1．课程设计应达到的目的 1、通过Matlab仿真熟悉课程设计的基本流程； 2、掌握控制系统的数学建模及传递函数的构造； 3、掌握控制系统性能的根轨迹分析； 4、学会分析系统的性能指标； 2．课程设计题目及要求 设计要求 1、进行系统总体设计，画出原理框图。（按给出的形式，自行构造数学模型，构造成1 个零点，三个极点的三阶系统，主导极点是一对共轭复根） G（s）=10(s+2)/(s+1)(s2+2s+6) 2、构造系统传递函数，利用MATLAB绘画系统的开环和闭环零极点图；(分别得 到闭环和开环的零极点图)参考课本P149页例题4-30 clear; num = [10,20]; den =[1 3 8 6]; pzmap(num,den) 3、利用MATLAB绘画根轨迹图，分析系统随着根轨迹增益变化的性能。并估算超 调量=16.3%时的K值（计算得到）。参考课本P149页例题4-31 clear num=[10,20]; den=[1 3 8 6]; sys=tf(num,den); rlocus(sys) hold on jjx(sys); s=jjx(sys); [k,Wcg]=imwk(sys) set(findobj('marker','x'),'markersize',8,'linewidth',1.5,'Color','k'); set(findobj('marker','o'),'markersize',8,'linewidth',1.5,'Color','k'); function s=jjx(sys) sys=tf(sys); num=sys.num{1}; den=sys.den{1}; p=roots(den); z=roots(num); n=length(p); m=length(z); if n>m s=(sum(p)-sum(z))/(n-m) sd=[]; if nargout<1 for i=1:n-m sd=[sd,s] end sysa=zpk([],sd,1); hold on; [r,k]=rlocus(sysa); for i=1:n-m plot(real(r(i,:)),imag(r(i,:)),'k:'); end end else disp; s=[]; end function [k,wcg]=imwk(sys) sys=tf(sys) num=sys.num{1} den=sys.den{1}; asys=allmargin(sys); wcg=asys.GMFrequency; k=asys. GainMargin; 排队论实验报告 《排队现象的建模、解析与模拟》 课程设计 姓名： 学号： 班级： 题目描述：排队系统的稳定性与什么有关？与系统的一步概率转移矩阵有什么关系？收敛速度快慢与什么有关？ 解答过程： （1）初始设定： 设初始状态X=(P1 P2 P3 … Pn),一步状态概率转移矩阵为P ，最终系统趋于稳定的状态为Y=(Y1 Y2 Y3 … Yn)，可知X 和Y 是一个固定不变的行向量，且P1+P2+P3+…+Pn=1，Y1+Y2+Y3+…+Yn=1。 （2）描述模型： 对排队系统最终趋于稳定的描述为：Y=X*P n ，n>N(N 是一个足够大的数)。 （3）提出假想： 由（2）中对于系统最终趋于稳定状态的描述，因为X 和Y 都是固定的向量，所以，若系统趋于稳定，则P n 收敛。假设P 最终收敛为 P σ=(a1 a2 ?an ???x1x2?xn ) ， 由概率转移矩阵的性质可知各行概率之和为1，即a1+a2+…+an=1。 因为Y* P σ= (Y1 Y2 Y3 … Yn)* (a1 a2 ?an ???x1x2?xn )=Y=(Y1 Y2 Y3 … Yn)，故提出猜测：概率转移矩阵收敛后各列的元素值相等。 (4)MATLAB 验证猜想： ① 当n ≥73时收敛： ② 当n≥38时收敛 ③ 当n≥11时收敛 ④ 当n≥3时收敛 ⑤ P本身就是收敛后的结果 (5)结论： 经过一系列验证，得出系统的稳定性只与一步转移概率矩阵P 有关，若P 收敛，则系统趋于稳定，反之系统不稳定。并且P 收敛后行和为1，每列元素值相同。 因为Y* P σ= (Y1 Y2 Y3 …… Yn)* (a1 a2 ?an ???a1a2?an ) =((Y1+Y2+Y3+…Yn)*a1 (Y1+Y2+Y3+…Yn)*a2 … (Y1+Y2+Y3+…Yn)*an) =(a1 a2 … an) 所以最终的概率分布的结果是矩阵收敛后的一行。 收敛速度快慢与一步概率转移矩阵每列元素值的分布有关，若每列元素值分布比较均匀，则收敛速度较快，反之收敛速度较慢。每列元素值相等的矩阵，本身就是收敛后的结果。单位阵是一个特例，它每列元素值不相等，但是单位阵收敛。与单位阵类似的一类矩阵，即 每列有且仅有一个1出现的矩阵，这类矩阵不会收敛。 实验3--- 多服务台排队系统的仿真 姓名：学号： 一、目标任务 已知一个系统有N 个服务员，能力相等，服务时间服从指数分布。顾客的到达时间间隔服从指数分布。用Monte-Carlo 仿真，分别求按下列方案的总体平均排队时间： ①M|M|N 。 ②N 个单通道系统并列，按1/N 概率分裂到达流。 ③N 个单通道并列，挑选最短的队。 要求： ①给出程序设计的过程。 ②如果采用固定的N,则要求N>2。 ③至少取p二和p二两种强度运行程序。 ④对结果进行分析。 二、编程语言 Matlab 三、关键代码 N = 3; % 服务员人数 r = 6; % 顾客到达流强度 u = 20; % 服务员服务强度 T = 1000000; % 仿真运行时间 avg_wait_time = []; % 平均等待时间 for i=1:100 % 模拟排队函数 server_time = [, , ]; % 用来保存服务员下一空闲时间 time = 0; % 绝对时钟，初始为 0 client_num = 0; % 顾客总数，初始为 0 CRTime = 0; % 顾客到达时间间隔 ServeTime = 0; % 顾客服务时间 server_id = 0 ; % 当前进入排队窗口的服务员编号 total_wait_time = 0;% 系统中到达顾客的总等待时间 while 1 按 1..N 的顺序循环排入服务 员窗口 if server_id ==0 server_id = N; end if server_time(1, server_id) <= time % 如果当前 server_id 号 服务员空闲， 则直接接收服务 server_time(1, server_id) = time + ServeTime; % 服务员下 一空闲时间为当 前绝对时钟加上当前服务时间 else % 否则所有服务员都在忙碌，顾客要排队等候 total_wait_time = total_wait_time + server_time(1, server_id) - time; % 顾客排队等候时间为当前服务员下一空闲时间减去绝对时 钟 server_time(1, server_id) = server_time(1, server_id) + ServeTime; end end avg_wait_time = [avg_wait_time, total_wait_time/client_num]; end % 计算平均等待时间 mean_avg_wait_time = mean(avg_wait_time); CRTime = exprnd(1/r); % 按指数分布产生顾客到达时间间隔 time = time + CRTime; % 更新系统的绝对时钟 if time > T break; end client_num = client_num + 1; % 顾客数加 1 ServeTime = exprnd(1/u); % 按指数分布产生顾客服务间隔 server_id = mod(client_num, N); % 课程名称：数学建模与数学实验学院： 专业： 姓名： 学号: 指导老师： 利用Monte Carlo方法模拟单服务台排队系统和多服务台排队系统 摘要 蒙特卡罗方法（Monte Carlo）又称统计模拟法随机抽样技术，是一种随机模拟方法，以概率和统计理论方法为基础的一种计算方法，是使用随机数（或更常见的伪随机数）来解决很多计算问题的方法。将所求解的问题同一定的概率模型相联系，用电子计算机实现统计模拟或抽样，以获得问题的近似解。本文通过两个具体的服务机构为例，分别说明如何利用蒙特卡洛方法模拟单服务台排队系统和多服务台排队系统。 单服务台排队系统（排队模型之港口系统）：通过排队论和蒙特卡洛方法解决了生产系统的效率问题，通过对工具到达时间和服务时间的计算机拟合，将基本模型确定在//1 M M排队模型，通过对此基本模型的分析和改进，在概率论相关理论的基础之上使用计算机模拟仿真（蒙特卡洛法）对生产系统的整个运行过程进行模拟，得出最后的结论。 多服务台排队系统（开水供应模型）：为了解决水房打水时的拥挤问题。根据相关数据和假设推导，最终建立了多服务窗排队M/G/n模型，用极大似然估计和排队论等方法对其进行了求解，并用Matlab软件对数据进行了处理和绘图。用灵敏度分析对结果进行了验证。本模型比较完美地解决了水房排队拥挤问题，而且经过简单的修改，它可以用于很多类似的排队问题。 关键词：蒙特卡洛方法，排队论，拟合优度，泊松流，灵敏度分析。 一、问题重述 港口排队系统：一个带有船只卸货设备的小港口，任何时间仅能为一艘船只卸货。船只进港是为了卸货，响铃两艘船到达的时间间隔在15分钟到145分钟变化。一艘船只卸货的时间有所卸货物的类型决定，在15分钟到90分钟之间变化。 开水供应系统：学院开水房的供水时间有限，水房面积有限，水管易受水垢堵塞。根据调查数据可知：通畅时几乎无人排队，堵塞时水房十分拥挤。由此可以看出水房设计存在问题，我们可以把开水房看成是一个随即服务系统，应用排队论的方法对系统运行状态做定量的描述。 二、基本假设 港口排队系统：通过对问题的重述，那么，每艘船只在港口的平均时间和最长时间是多少？ 若一艘船只的等待时间是从到达到开始卸货的时间，每艘船只的平均等待时间和最长等待时间是多少？ 卸货设备空闲时间的百分比是多少？ 船只排队最长的长度是多少？ 开水供应系统： 假设Ⅰ、顾客流满足参数为λ的Poisson分布，其中λ为单位时间到达的顾客平均数。每个顾客所需的服务时间相互独立，顾客流是无限的，在观测期间平稳。 假设Ⅱ、排队方式为单一队列的等候制，先到先服务。虽然水房内有多个服务台，每个服务台都有自己的队列，但同时顾客总是自由转移到最短的队列上，不可能出现有顾客排队而服务器空闲的情况。本文最后对两种排队方式的比较也表明这一假设是合理的。 假设Ⅲ、水房共有20个并联的服务台（水龙头），设每个服务台的服务时间服从某个相同的分布，t和σ分别是服务时间的均值和均方差，γ=σ/ t为偏离系数。由于锅炉及输水管容量的限制，使t依赖于正在进行服务的水龙头个数m，设此时平均服务时间t(m)。且存在一临界值当m<= m0 时，t(m)为常数 《通信网理论基础》 实验二：二次排队问题——M/M/1排队系统的级联一、实验目的 M/M/1是最简单的排队系统，其假设到达过程是一个参数为λ的Poisson过程，服务时间是参数为μ的负指数分布，只有一个服务窗口，等待的位置有无穷多个，排队的方式是FIFO。 M/M/1排队系统的稳态分布、平均队列长度，等待时间的分布以及平均等待时间，可通过泊松过程、负指数分布、生灭过程以及Little公式等进行理论上的分析与求解。 本次实验的目标有两个： 实现M/M/1单窗口无限排队系统的系统仿真，利用事件调度法实现离散事件系统仿真，并统计平均队列长度以及平均等待时间等值，以与理论分析结果进行对 比。 仿真两个M/M/1级联所组成的排队网络，统计各个队列的平均队列长度与平均系统时间等值，验证Kleinrock有关数据包在从一个交换机出来后，进入下一个 交换机时，随机按负指数分布取一个新的长度的假设的合理性。 二、实验原理 1、M/M/1排队系统 根据排队论的知识我们知道，排队系统的分类是根据该系统中的顾客到达模式、服务模式、服务员数量以及服务规则等因素决定的。 设到达过程是一个参数为λ的Poisson过程，则长度为t的时间内到达k个呼叫的概率 )(t P k 服从Poisson 分布，即()()! k t k t P t k e λλ-= ，?????????=,2,1,0k ，其中λ>0为 一常数，表示了平均到达率或Poisson 呼叫流的强度。设每个呼叫的持续时间为 i τ，服从 参数为μ的负指数分布，即其分布函数为 {}1,0t P X t e t μ-<=-≥.服务规则采用先进 先服务的规则（FIFO ）。 在该M/M/1系统中，设 λρμ= ，则稳态时的平均队长为[]1E N ρ ρ= -，顾客的平均 等待时间为 1 T μλ= -。 2、 二次排队网络 由两个M/M/1排队系统所组成的级联网络，顾客以参数为λ的泊松过程到达第一个排队系统A ，服务时间为参数为1μ的负指数分布；从A 出来后直接进入第二个排队系统B ，B 的服务时间为参数为 2μ的负指数分布，且与A 的服务时间相互独立。 在该级联网络中，如稳态存在，即 1λμ<且2λμ<，则两个排队系统相互独立，顾客 穿过网络的总时延为各个排队系统的时延之和，即 1211 T μλμλ= + --。 如将该模型应用于数据包穿越网络的平均时延的计算，假设数据包的包长服从负指数分布，平均包长为b ；排队系统A 的信道速率为 1C ，B 的信道速率为2C 。为保证两次排队的 独立性，Kleinrock 假设数据包在从一个交换机出来后，进入下一个交换机时，随机按负指数分布取一个新的长度。 三、 实验内容单服务台排队系统的仿真
matlab控制系统仿真课程设计
matlab单服务台排队系统实验报告
单服务台排队系统仿真研究报告
控制系统仿真课程设计
排队系统实验报告
单服务台排队系统建模与仿真研究报告
单服务台排队系统仿真报告
控制系统仿真课程设计
排队论
多服务台排队系统的仿真
数学建模论文(蒙特卡罗的多服务台和单服务台排队系统)
北邮通信网性能分析实验二MM1排队系统实验报告