基于蒙特卡罗模拟的银行排队业务的建模与仿真

基于蒙特卡罗模拟的银行排队业务的建模与仿真

1蒙特卡罗模拟的实际应用随着国内银行市场化的全面转型，银行代扣代缴的中间业务，如水电、电费、煤气费、工资和社会保障基金等，这增加了银行的业务领域的压力。最近，由于样本额的折扣，贷款的折扣和提前贷款的突然增加。与此同时，随着股市和基金的升高，银行推出的基金和产品也出现了大量的董事会业务，这些业务相对复杂，需要很长时间。多种因素使银行团队的矛盾集中。银行网点排队问题成为社会舆论和新闻媒体关注的焦点。本文结合运筹学排队论问题的基本特点,利用蒙特卡罗模拟在刻画排队系统动态性方面的优势,重点介绍了利用Excel强大的函数、公式编辑功能以及加载宏CrystaBall软件建立模型,实现蒙特卡罗模拟的方法在处理实际排队问题中的应用。本文第3部分重点介绍应用Excel实现对某银行排队系统各运行指标的动态模拟,并对结果进行了敏感性分析,取得了良好的效果。本文最后介绍了该方法在实际应用中的意义,并给出全文的总结。2队列问题的方法论2.1长时间和保证时间解决排队问题的目的,是研究排队系统的运行效率,估计服务质量,确定系统参数的最优值,以确定系统结构是否合理、研究设计改进措施等。所以必须确定用以判断系统运行优劣的数量指标,求出这些数量指标的概率分布。这些数量指标通常是,队长L:系统中顾客数的期望值;排队长Lq:等待排队的顾客数期望值;逗留时间W:顾客在系统中的停留时间期望值;等待时间Wq:排队等待的期望值;空闲率I:服务台空闲的概率。这些参数反映的都是稳定状态下的期望值,没有体现系统的动态特性。在一个实际的等待系统中,这些值会围绕着期望值有较大的波动。当系统开始为空值时,可能需要很长的时间才能达到这种期望水平。这些数量指标反映了系统长期运行下的静态参数,而仿真则可以模拟排队的动态行为信息等一系列优于一般静态指标分析的特点。相比其他使用设计复杂的程序和烦琐的计算机编程语言,如Matlab6.5、CISMAN、VisualBasic,用Excel建模具有直观简便、操作性更强、适用性更广等优点。2.2不确定性问题模拟蒙特卡罗方法又名随机模拟法或统计试验法,属于试验数学的一个分支。它利用随机数进行统计试验,以求得的统计特征值(如均值、概率等)作为待解问题的数值解,所作的统计试验称为蒙特卡罗模拟。蒙特卡罗模拟的本质是从概率分布中重复抽样以建立输出变量的分布。其原理是:在所研究的系统中,采用某种特定方法产生随机数和随机变量,仿真随机事件,对结果进行统计处理,从而得到问题的解。在解决风险分析、项目管理等不易处理的不确定性问题时,该方法特别有用。在应用蒙特卡罗方法解决实际问题的时候主要有两部分工作:(1)用蒙特卡罗方法模拟某一过程时,确定需要产生各种概率分布的随机变量。(2)用统计方法把模型的数字特征估计出来,从而得到实际问题的数值解。CrystalBall软件是由美国Decisioneering公司开发的,为Excel电子表格提供的功能强大的加载宏。在CrystalBall中,将随机的输入变量规定为假设单元(assumptioncell),将感兴趣的随机输出变量规定为预测单元(forecastcell)。CrystalBall的DistributionGallery为随机输入变量提供了多种可供选择的分布类型,能够适应常见的分布需要。以上这些功能使得在Excel中运用CrystalBall进行蒙特卡罗仿真可以非常方便地实现。排队系统涉及的多是关于时间先后关系的变量,只需用公式表示出各顾客到达时间变量的逻辑关系,将顾客的到达率和业务员的服务率所在单元格设置为假设单元,为其设定一个适当的分布,将需要重点考量的顾客平均等待时间、平均排队长、平均空闲率等所在的单元格设置为预测单元,最后设定重复抽样次数即可实现模拟。要得到一个未知总体参数的一个优良估计,必须模拟仿真很多次。3fcfs排他系统的分析例如,某银行营业网点现只开放一个窗口,在其正常的营业时间内,顾客的到达率λ服从U(8,10)的均匀分布,业务员的服务率μ服从U(10,12)的均匀分布。求顾客等待时间大于1小时的概率。现对该系统做如下假设:(1)顾客到达时间间隔服从指数分布,且各顾客到达间隔时间相互独立;(2)服务时间服从指数分布,且各顾客的服务时间相互独立;(3)顾客按照FCFS规则排队接受服务;如果系统忙,顾客愿意等待;队伍长短无限制。这样该问题正好满足M/M/1型的排队系统的条件。解决这种复杂的动态问题,模拟方法体现出较大的优越性。3.1基于党组织有限元的模型模拟变量定义如下:AT(i)=顾客i的到达时刻;BT(i)=顾客i的开始服务时刻;FT(i)=顾客i的完成服务时刻;WT(i)=顾客i的等待时间;IT(i)=顾客i-1和顾客i之间的空闲时间;NIQ(i)=队列中等待的顾客数;ST(i)=顾客i的服务时间;TBA(i)=顾客i-1和顾客i的到达时间差;i=1,2,…,n。其中,顾客到达时间间隔TBA(i)、顾客i的服务时间ST(i)是假设变量,其他变量都可以由这两个变量表达出来。模拟的逻辑流程图如图1所示。第一步,将CrystalBall宏文件加载至Excel,Excel界面菜单栏会增加Cell,Run和CBTools3个新按钮条。利用Excel的公式编辑功能将图1所示的逻辑关系表示出来,如图2所示。使用“Ctrl+~”组合键,可将公式图切换至仿真结果图,如图3所示。第二步,使用Cell菜单下的DefineAssumption功能分别将顾客平均到达率、业务员平均服务率所在单元格B3、B4设置为假设变量,并设置其服从的分布。因为M1003、O1003、Q1003分别表示第1000名顾客到达时队列平均顾客数、平均等待时间、平均空闲率的值,可以将其作为预测变量显示在B12:B14所在的结果中,并使用Cell菜单下的DefineForecasts将其定义为预测变量。第三步,使用Run菜单下的Run功能,就可以实现对模型的仿真模拟。关于建模公式图(图2)的两点说明:(1)D列是由CrystalBall随机数生成功能产生的顾客到达间隔TBA(i),如产生以模型中顾客到达率λ为参数的指数分布的随机数,只需在D4中输入公式=CB.Exponential($B$3)”。同样H列需要产生服务率为μ的指数分布的随机顾客服务时间ST(i),同样也可在H4中输入公式“=CB.Exponential($B$4)”。(2)在编辑完第2位顾客到达的各变量的公式后,选中“$D$5:$Q$5”单元格所在区域,单击鼠标拖动区域右下方的黑点(拖动柄)并一直拖动到第1003行,即可实现对1000名到达的顾客的所有变量的模拟。3.2顾客等待时间研究结果表明:客户排队等候超过10分钟,情绪开始急躁;超过20分钟,情绪表现为厌烦;超过40分钟,常因恼火而拂袖离去。顾客对于服务质量的感知与等待时间等因素有关,顾客普遍不喜欢等待。本文用顾客的平均等待时间来衡量服务质量和客户满意度。用CrystalBall对模型进行10000次模拟的结果如图4所示。可知顾客的平均等待时间为0.3171小时,等待时间超过0.6667小时(40分钟)的概率为8.58%,超过1小时的概率为3.52%。随着银行业竞争的日益加剧,处于长期等待状态的这部分顾客会对银行的服务质量增加不满和抱怨,银行将面临损失这部分顾客的可能。3.3平均到达率和平均等待时间的关系灵敏度分析是分析一个系统或模型的状态或输出变化对系统参数或周围条件变化的敏感程度的方法。通过灵敏度分析可以确定哪些参数对系统或模型有较大的影响。CrystalBall对顾客平均等待时间的灵敏度分析,见图5。图5显示的灵敏度分析是由秩相关系数度量的,秩相关系数提供了各种假设和预测变量共同变化程度的度量。图5显示,顾客平均到达率λ、业务员平均服务率μ和平均等待时间的相关度为0.84、-0.36。表明顾客平均到达率λ和增加平均等待时间之间存在显著正相关,业务员平均服务率μ和减少平均等待时间存在一定负相关。可知,业务员平均服务率μ的提升在一定程度上可以缩短顾客等待时间。本文利用Excel的加载宏软件CrystalBall对一个简化的银行排队问题进行模拟和分析。对于目前的银行排队难这一问题,希望银行短期内通过快速开设新的网点加以改善的做法难切实际,而通过增加现有网点服务窗口的个数则是可能被银行接受的方案。增加窗口的投入和提高的顾客满意度之间存在什么关系,很难做出定量的回答。而本文提供的方法,避开了设计复杂的程序和烦琐的计算机编程语言的的障碍,利用一般管理人员熟悉的Excel实现蒙特卡罗模拟,为银行的设计者和管理者方便地做出切合实际的解决方案提供科学的根据和启示。但由于搜集的实际数据