用蒙特卡洛法实现对排队等待问题的计算机模拟

用蒙特卡洛法实现

对排队等待问题的计算机模拟蒙特卡洛(MonteCarlo)法，或称统计试验法、计算机随机模拟方法，起源于美国在第一次世界大战进研制原子弹的“曼哈顿计划”。该计划的主持人之一、数学家冯•诺伊曼用驰名世界的赌城一摩纳哥的MonteCarlo—来命名这种方法，为它蒙上了一层神秘色彩。一、蒙特卡洛法的基本思想及其应用MonteCarlo方法是一种基于“随机数”，采用统计抽样方法，近似求解数学问题或物理问题的过程。把统计模拟法用于数值计算已有200多年的历史，最早是法国数学家蒲丰(1707-1788)。他进行了著名的“蒲丰投针实验”,早以此来求圆周率丸的近似值。本世纪40年代，随着电子计算机的出现，特别是近年来高速电子计算机的出现，使得用数学方法在计算机上大量、快速地模拟这样的试验成为可能。统计试验法通常用来研究概率过程，研究问题时常涉及下列一些与随机因素有关的概率，如各类概率等，一般来说，建立描述过程的复杂的概率模型是不成问题的，但用数学方法研究与分析这些模型是却很困难，问题的维数(即变量的个数)可能高达数百甚至数千。对这类问题，难度随维数的增加呈指数增长，这就是所谓的维数的灾难”(CourseDimensionality)o传统的数值方法难以对付(即使使用速度最快的计算机)，甚至达到了无法进行的地步。因此，唯一可取的研究方法是统计实验法。统计模拟(蒙特卡洛法)，在系统工程中的应用日益广泛，据国夕卜有关文献报道其应用领域大致有：.航空运输排队，机场设计等；港口设计，泊位研究等；消防车或救护车的布局和调派；城市公共汽车作业调度；出租汽车调度计划；铁路货运调度计划；加油站、停车场等设计；售票所布局；.存储模拟，仓库布局等；.设备维修计划；11.生产过程的安排；.工厂的单件、小批生产的作业计划；.销售预测；二、排队或等待问题的分析在日常生活中，我们每天都会遇到各种各样的排队。比如：银行中取款要排队，火车站买票要排队，超市、商场中购物付款要排队，预订旅馆或机票时也要排队，人们仿佛置身于一个排队的社会。所谓排队，就是等候消费服务的顾客在进入点前排队（意大利G•佩里切利）。排队的类型一般来说有以下几种：一个服务点。即只有一个点可供顾客选择等候服务，顾客只能按顺序一个接一个等侯，这样排队，顾客多的时候最混乱。现在这种排队类型已不多见。多个服务点。即有两个或两个以上的点可供顾客选择，每一点都可排成一条队。多个服务点可减轻一个服务点的负担，增加顾客选择服务点的灵活性，同时顾客还可选择自己喜欢的服务员。专门服务点。即专门为某些特殊的顾客开辟的服务点。专门的服务点可以保证某些特殊顾客的特殊权益，如：减少顾客等待服务时间，同时使一些顾客分离出来，减轻其服务口的压力。我们拿一个理发店的实际例子来分析：某理发馆每天早晨开业后，顾客的到来总是陆续不断的，他们到来的间隔时间，给统计如下：到达时间间隔（分钟）频率40.0750.1060.52

70.2080.11顾客到来后，由管理人员引导至A、B、C三个服务椅（理发员固定在服务椅上），各服务椅的服务内容相同，但服务时间有多有少，根据以往统计资料，他们的服务时间分布如下：ABC时间（分）频率时间（分）频率时间（分）频率80.18100.18120.1590.22110.19130.22100.33120.35140.36110.27130.28150.271.001.001.00当有一个以上服务椅空闲时，管理人员按字母排列次序引导顾客至服务椅去进行服务。顾客中有10%属“特殊顾客”理发很费时间，每人需多4分钟。现模拟10位顾客的到来和服务情况，并分析三个服务椅的忙闲情况。顾客到达时间间隔频率累计频率对应的随机数A服务时间频率累计频率对应的随机数40.070.070.00〜0.0680.180.180.00〜0.1750.100.170.07〜0.1690.220.400.18〜0.3960.520.690.17〜0.68100.330.730.40〜0.7270.200.890.69〜0.88110.271.001.0080.111.000.89〜1.00B服务时间频率累计频率对应的随机数C服务时间频率累计频率对应的随机数100.180.180.00〜0.17120.150.150.00〜0.14110.190.370.18〜0.36130.220.370.15〜0.36120.350.720.37〜0.71140.360.730.37〜0.72130.281.000.72〜1.00150.371.000.72〜1.00我们采用二位随机数来模拟顾客到来间隔和服务情况，并选定碰到个位数是5，就作为“特殊顾”（占10%）现从随机数中抽取随机数如下：

顾客顾客到来时间间隔服务随机数随机数时间间隔（分）12345678910（14）656686678009946201119663260409020336501（特殊顾客）72689287设第一个顾客到达的时间是在3点正，随机数是14,此数本不需要，主要是看一下是不是“特殊顾客”，现此数尾数不是5说明第一位不是“特殊顾客”。第一个顾客到达，按字母的顺序应到A服务椅，随机数00，表示要服务8分钟，第二个顾客到达与第一个顾客的间隔时间，按随机数46应为6分钟，即在8.06时，服务椅A正在使用，因此到服务椅B，服务随机数99，应服务13分钟，如此编成表如下：顾客到达时间服务椅等待时间服务时间ABC到达离开到达离开到达离开18:00A88:008:0828:06B138:068:1938:11C138:118:2448:17A98:178:2658:23B118:238:3468:31A108:318:4178:37B108:378:4788:43A128:438:5598:50B128:509:02108:58A118:589:09共计109*为“特殊顾客”服务时间8（分）+4（分）=12（分）10个顾客在系统内的总时间为63分。分析：.10个顾客在系统(理发馆内)的总时间为69分，平均每人69/10=6.9分.10个顾客总的服务时间为109分。平均每人服务109/10=10.9分.统内无等待时间，即没出现顾客排队等待问题.4•椅的服务情况如下：ABC服务顾客数总的服务时间(分)541504613为每个顾客的平均服务时间1011.513192356空闲时间空闲占总的时间％27.3%33.3%81.2%三、排队等待问题的计算机模拟在此基础上，可用计算机编写程序，在计算机上实现仿真模拟，用以模拟更多时间和更多顾客时的情况。以下我我组用Matlab编写的模拟程序的算法部分：function[num,pass]=computing(tim0)seat=[000];%服务员属性pass=rand(1,4);%顾客信息：序号、到达时间、特殊要求时间、正常理发时间pass(5)=0;%理发员pass(6)=0;%离开时间pass(7)=0;%等待时间num=1;%服务人数初始化tim=0;%时间计数器temp=0;%whiletim<=tim0pass(num,1)=num;%装入序号pass(num,2)=rand;pass(num,3)=rand;pass(num,4)=rand;%计算顾客到达时间ifpass(num,2)<=0.07temp=4;elseifpass(num,2)<=0.17

temp=5;elseifpass(num,2)<=0.69temp=6;elseifpass(num,2)<=0.89temp=7;elsetemp=8;endendendendtim=tim+temp; %装入顾客到达时间pass(num,2)=tim;ifpass(num,3)<=0.1pass(num,3)=4;%装入需要特殊服务的时间elsepass(num,3)=0;endnum=num+1;endnum=num-1;fori=1:num%计算顾客的理发席位ifseat(1)<=pass(i,2)+pass(i,7)pass(i,5)=1;%由1号服务员理发temp1=timinge1(1,pass(i,4));seat(1)=pass(i,2)+pass(i,3)+temp1;pass(i,4)=temp1;%装入正常理发所需时间elseifseat(2)<=pass(i,2)+pass(i,7)pass(i,5)=2;%由2号服务员理发temp1=timinge1(2,pass(i,4));seat(2)=pass(i,2)+pass(i,3)+temp1;pass(i,4)=temp1;%装入正常理发所需时间elseifseat(3)<=pass(i,2)+pass(i,7)pass(i,5)=3;%由2号服务员理发temp1=timinge1(3,pass(i,4));seat(3)=pass(i,2)+pass(i,3)+temp1;pass(i,4)=temp1;else %计算等待时间x=seat(1);y=1;ifx>seat(2)x=seat(2);y=2;endifx>seat(3)x=seat(3);y=3；endpass(i,5)=y;temp1=timinge1(y,pass(i,4));pass(i,7)=seat(y)-pass(i,2);seat(y)=seat(y)+temp1+pass(i,3);pass(i,4)=temp1;endendendpass(i,6)=pass(i,2)+pass(i,3)+pass(i,4);endfunctionxxxx=timinge1(vect1,vect)switchvect1case1ifvect<=0.18xxxx=8;elseifvect<=0.4xxxx=9;elseifvect<=0.77xxxx=10;elsexxxx=11;endendendcase2ifvect<=0.18xxxx=10;elseifvect<=0.37xxxx=11;elseifvect<=0.72xxxx=12;elsexxxx=13;endendendotherwiseifvect<=0.15xxxx=12;

elseifvect<=0.37xxxx=13;elseifvect<=0.74xxxx=14;elsexxxx=15;endendendend下面是在一次模拟过程中产生的顾客用户信息表:理发店顾客数据统计序号到达时间正常理发所需时间特殊要求时间等待时间离开时间理发员序号到达时间正常理发所需时间特殊要求时间等待时间离开时间理发员18100018158369140038332121200242593758003831318110029160381130039424251200372613871100398153310004316239315004083641114056263398004071

97488406016440110041528541400683654180041919608006816641900428116713008026742120043721731100841684380043911791400933694313004502186104010017044110045411921200104271441500464319913001123724580046311101000115173461200472211110001212744610004781

718111110012917547100048421121200137276481100491121310001431774812004982214110015127849100050212141000156179501140515221513001662805010005141215100016918151140052632161400179382511000528121711001821835212405402217100018728453110054212189001941855314005523

85821910002022865411005561319900205187551100562232011002132885590056613201140223189561200575232112002262905710005801322130023339157120058823221000236192581000592132312002442935813006022323100024919459100060513241000255295601200615232510402651966010006191

91942510002672976112006272426150027939862110063314271000281199631140645242711002882106310006461428100029511064140065434291200304210641000658143090030911065130066724301000317210651100670143110003231106615006803431120033121067100068115329403361106712406952

03795321400341310689006941533124034821069134070835338403501116910007071534120035631170130071625351000360111711100721153512403732117111007272536800371111728007301从该表中，我们可以做出如下分析:.共计114个顾客中，在系统（理发馆内）的总时间为1319分，平均每人1319/114=11.6分.系统内无等待时间，即没出现顾客排队等待问题.3•理发员的工作情况如下：ABC服务顾客数554415565534211总的服务时间（分）10.312.114.1155186509为每个顾客的平均服务时间21.5%25.8%70.7%空闲时间空闲占总的时间％四、小结在排队等待问题中，存在着在随机的某种不确定需要，而这种不确定的需要，还具备了一定的统计规律，用蒙特卡洛法对此类问题求解，可节省时间和劳力，大大提高工作效率同时确保计算的高效、准确，对可靠性设计提供了保证。使我们对解决问题的决策有可能更加明智，从而能够有效地解决问题。缺点是为提高计算精度必须进行大量模拟计算，需要较长的计算时间和编程时间。附：用matlab7.0实现该项模拟的源程序：clearallclosealltim0=720;%需要模拟的时间[num,pass]=computing(tim0);%计算模拟数据%动画制作envirment %场景title1=annotation('textbox','Position',[0.33770.87120.33480.07885],...'EdgeColor','none','FitHeightToText','off,'FontName','Arial','FontSize',16,...'FontWeight','bold','String',{理发店忙闲情况分析'});time1=annotation('textbox','Position',[0.027540.10190.09420.06538],...'EdgeColor','none','FitHeightToText','off,'FontSize',14,'FontWeight','bold',...'String',{'时间'});hour1=annotation('textbox','Position',[0.10720.10380.075360.06346],...'EdgeColor','none','FitHeightToText','off,'FontSize',14,...'FontWeight','bold','String',{'08:'});minute1=annotation('textbox','Position',[0.14930.10380.072460.06538],...'EdgeColor','none','FitHeightToText','off,'FontSize',14,...'FontWeight','bold','String',{'10'});counter1=1;%时间计数temp1=1; %顾客计数counter=1;counterxy=zeros(7,2);tempa=0;tempb=0;tempc=0;tempd=0;tempe=0;tempf=0;tempg=0;globalman1mana1mana2mana3mana4mana5globalmanb1manb2manb3manb4manb5globalmanc1manc2manc3manc4manc5globalmand1mand2mand3mand4mand5globalmane1mane2mane3mane4mane5man1(1,:)=[0.46250.17860.09820.1381];%[x0y0x1y1]man1(2,:)=[0.48750.26670.50180.2667];%[x0y0x1y1]man1(3,:)=[0.52140.26670.53390.2667];%[x0y0x1y1]man1(4,:)=[0.51250.250.51250.23]; %[x0y0x1y1]man1(5,:)=[0.49680.20810.5250.2072]; %[x0y0x1y1]axisoffwhilecounter1<=10%时间计数%显示时间minutex=rem(counter1,60);hourx=8+(counter1-minutex)/60;set(hour1,'String',{hourx})set(minute1,'String',{minutex})iftemp1<=num%离开与删除人脸对象iftempa〜=0ifpass(tempa,6)==counter1leaveflash(1,counterxy(1,1),counterxy(1,2))tempa=0;delete(mana1,mana2,mana3,mana4,mana5);endendiftempb〜=0ifpass(tempb,6)==counter1leaveflash(2,counterxy(2,1),counterxy(2,2))tempb=0;delete(manb1,manb2,manb3,manb4,manb5);endendiftempc〜=0ifpass(tempc,6)==counter1leaveflash(3,counterxy(3,1),counterxy(3,2))tempc=0;delete(manc1,manc2,manc3,manc4,manc5);endendiftempd〜=0ifpass(tempd,6)==counterlleaveflash(4,counterxy(4,1),counterxy(4,2))tempd=0;delete(mand1,mand2,mand3,mand4,mand5);endendiftempe〜=0ifpass(tempe,6)==counterlleaveflash(5,counterxy(5,1),counterxy(5,2))tempe=0;delete(mane1,mane2,mane3,mane4,mane5);endendiftempf〜=0ifpass(tempf,6)==counter1leaveflash(6,counterxy(6,1),counterxy(6,2))tempf=0;delete(manf1,manf2,manf3,manf4,manf5);endendiftempg〜=0ifpass(tempg,6)==counter1leaveflash(7,counterxy(7,1),counterxy(7,2))tempg=0;delete(mang1,mang2,mang3,mang4,mang5);endend%产生人脸ifpass(temp1,2)==counter1iftempa==0[mana1,mana2,mana3,mana4,mana5]=personcreat(pass(temp1,3))tempa=temp1;ifpass(temp1,5)==1counterxy(1,:)=[-0.0270.055];endifpass(temp1,5)==2counterxy(1,:)=[00.055];endifpass(temp1,5)==3counterxy(1,:)=[0.0270.055];endtemp1=temp1+1;elseiftempb==0[manb1,manb2,manb3,manb4,manb5]=personcreat(pass(temp1,3))tempb=temp1;ifpass(temp1,5)==1counterxy(2,:)=[-0.0270.055];endifpass(temp1,5)==2counterxy(2,:)=[00.055];endifpass(temp1,5)==3counterxy(3,:)=[0.0270.055];endtemp1=temp1+1;elseiftempc==0[manc1,manc2,manc3,manc4,manc5]=personcreat(pass(temp1,3))tempc=temp1;ifpass(temp1,5)==1counterxy(3,:)=[-0.0270.055];endifpass(temp1,5)==2counterxy(3,:)=[00.055];endifpass(temp1,5)==3counterxy(3,:)=[0.0270.055];endtemp1=temp1+1;elseiftempd==0[mand1,mand2,mand3,mand4,mand5]=personcreat(pass(temp1,3))tempd=temp1;ifpass(temp1,5)==1counterxy(4,:)=[-0.0270.055];endifpass(temp1,5)==2counterxy(4,:)=[00.055];endifpass(temp1,5)==3counterxy(4,:)=[0.0270.055];endtemp1=temp1+1;elseiftempe==0[mane1,mane2,mane3,mane4,mane5]=personcreat(pass(temp1,3))tempe=temp1;ifpass(temp1,5)==1counterxy(5,:)=[-0.0270.055];endifpass(temp1,5)==2counterxy(5,:)=[00.055];endifpass(temp1,5)==3counterxy(5,:)=[0.0270.055];endtemp1=temp1+1;elseiftempf==0[manf1,manf2,manf3,manf4,manf5]=personcreat(pass(temp1,3))tempf=temp1;ifpass(temp1,5)==1counterxy(6,:)=[-0.0270.055];endifpass(temp1,5)==2counterxy(6,:)=[00.055];endifpass(temp1,5)==3counterxy(6,:)=[0.0270.055];endtemp1=temp1+1;elseiftempg==0[mang1,mang2,mang3,mang4,mang5]=personcreat(pass(temp1,3))tempg=temp1;ifpass(temp1,5)==1counterxy(7,:)=[-0.0270.055];endifpass(temp1,5)==2counterxy(7,:)=[00.055];endifpass(temp1,5)==3counterxy(7,:)=[0.0270.055];endtemp1=temp1+1;endendend

endendendendend%开始服务iftempa〜=0ifpass(tempa,2)+pass(tempa,7)==counterlserveflash(1,counterxy(1,1),counterxy(1,2))endendiftempb〜=0ifpass(tempb,2)+pass(tempb,7)==counter1serveflash(2,counterxy(2,1),counterxy(2,2))endendiftempc〜=0ifpass(tempc,2)+pass(tempc,7)==counter1serveflash(3,counterxy(3,1),counterxy(3,2))endendiftempd〜=0ifpass(tempd,2)+pass(tempd,7)==counter1serveflash(4,counterxy(4,1),counterxy(4,2))endendiftempe〜=0ifpass(tempe,2)+pass(tempe,7)==counter1serveflash(5,counterxy(5,1),counterxy(5,2))endendiftempf〜=0ifpass(tempf,2)+pass(tempf,7)==counterlserveflash(6,counterxy(6,1),counterxy(6,2))endendiftempg〜=0ifpass(tempg,2)+pass(tempg,7)==counterlserveflash(7,counterxy(7,1),counterxy(7,2))endendendcounter1=counter1+1;endiiii=1;sumtime=0;sumtimea=0;sumtimeb=0;sumtimec=0;numa=0;numb=0;numc=0;whileiiii<numsumtime=sumtime+pass(iiii,3)+pass(iiii,4);ifpass(iiii,5)==1sumtimea=sumtimea+pass(iiii,3)+pass(iiii,4);numa=numa+1;endifpass(iiii,5)==2sumtimeb=sumtimeb+pass(iiii,3)+pass(iiii,4);numb=numb+1;endifpass(iiii,5)==3sumtimec=sumtimec+pass(iiii,3)+pass(iiii,4);numc=numc+1;endiiii=iiii+1;endnum=num-1;freetimea=720-sumtimea;freetimeb=720-sumtimeb;freetimec=720-sumtimec;freetime=freetimea+freetimeb+freetimec;avetimea=round(1000*sumtimea/720)/1000;avetimeb=round(1000*sumtimeb/720)/1000;avetimec=round(1000*sumtimec/720)/1000;avetime=round(1000*sumtime/(3*720))/1000;scaleaa=freetimea/720;scaleab=freetimeb/720;scaleac=freetimec/720;dispdata(1,:)=[numanumbnumc];dispdata(2,:)=[sumtimeasumtimebsumtimec];dispdata(3,:)=[avetimeaavetimebavetimec];dispda芟^n^oaimoafreaimobfreaimoc」；dispda_a(5o''freaimea/720freaimeb/720freaimec\720rfigure(2)subplor(221)息usebar(dispdaMlo)H-A漏澈殍DMW)subplor(222，pauso多夏^澈耳画-)bar(dispdaA2」))subplor(223，pauso多夏长格蜀澈耳画-)pAdispdaA3o)subplor(224，pausoH-lAH)}_ai民w)pie(dispda_a(5」))functionInum9passllcompudng(dmo)sgvo0。「％蜀澈9D1DWpassgnd(L4)」％驯血，四院耳画，茹海m好耳画，用Mffl涎耳画pass(5np%^^aDpass(6np%®Jfs画pass(7np%4^^s画dmno-%耳画牛肄|#Empnp%whileHmAHu.mo，25,pass(num,1)=num;%装入序号pass(num,2)=rand;pass(num,3)=rand;pass(num,4)=rand;%计算顾客到达时间ifpass(num,2)<=0.07temp=4;elseifpass(num,2)<=0.17temp=5;elseifpass(num,2)<=0.69temp=6;elseifpass(num,2)<=0.89temp=7;elsetemp=8;endendendendtim=tim+temp; %装入顾客到达时间pass(num,2)=tim;ifpass(num,3)<=0.1pass(num,3)=4;%装入需要特殊服务的时间elsepass(num,3)=0;endnum=num+1;endnum=num-1;fori=1:num%计算顾客的理发席位ifseat(1)<=pass(i,2)+pass(i,7)pass(i,5)=1;%由1号服务员理发temp1=timinge1(1,pass(i,4));seat(1)=pass(i,2)+pass(i,3)+temp1;pass(i,4)=temp1;%装入正常理发所需时间elseifseat(2)<=pass(i,2)+pass(i,7)pass(i,5)=2;%由2号服务员理发temp1=timinge1(2,pass(i,4));seat(2)=pass(i,2)+pass(i,3)+temp1;pass(i,4)=temp1;%装入正常理发所需时间elseifseat(3)<=pass(i,2)+pass(i,7)pass(i,5)=3;%由2号服务员理发temp1=timinge1(3,pass(i,4));seat(3)=pass(i,2)+pass(i,3)+temp1;pass(i,4)=temp1;else %计算等待时间x=seat(1);y=1;ifx>seat(2)x=seat(2);y=2;endifx>seat(3)x=seat(3);y=3;endpass(i,5)=y;temp1=timinge1(y,pass(i,4));pass(i,7)=seat(y)-pass(i,2);seat(y)=seat(y)+temp1+pass(i,3);pass(i,4)=temp1;endendendpass(i,6)=pass(i,2)+pass(i,3)+pass(i,4);endfunctionenvirment%场景设置pict=figure('color',[0.750.750.75],'position',[5050690520]);door=annotation('rectangle',[0.42140.1405 0.19460.03333],'FaceColor',[010],'EdgeColor',[010]);word1=annotation('textbox','Position',[0.33930.1143 0.08750.07619],'EdgeColor','none','FitHeightToText','off,'FontSize',12,'String',{门口'});seata1=annotation('rectangle',[0.20540.8214 0.11070.02857],'FaceColor',[0.84710.16080],'EdgeColor',[0.84710.16080]);seata2 = annotation('rectangle',[0.2071 0.7476 0.014290.1024],'FaceColor',[0.84710.16080],'EdgeColor',[0.84710.16080]);seata3 = annotation('rectangle',[0.3021 0.7419 0.014290.1024],'FaceColor',[0.84710.16080],'EdgeColor',[0.84710.16080]);word2=annotation('textbox','Position',[0.16070.65 0.11070.06667],'EdgeColor','none','FitHeightToText','off,'String',{'SeatA'});seatbl=annotation('rectangle',[0.45790.8171 0.11070.02857],'FaceColor',[0.84710.16080],'EdgeColor',[0.84710.16080]);seatb2 = annotation('rectangle',[0.4579 0.7386 0.014290.1024],'FaceColor',[0.84710.16080],'EdgeColor',[0.84710.16080]);seatb3 = annotation('rectangle',[0.5564 0.7448 0.014290.1024],'FaceColor',[0.84710.16080],'EdgeColor',[0.84710.16080]);word3=annotation('textbox','Position',[0.37680.6381 0.10360.07381],'EdgeColor','none','FitHeightToText','off,'String',{'SeatB'});seatc1=annotation('rectangle',[0.69860.8157 0.11070.02857],'FaceColor',[0.84710.16080],'EdgeColor',[0.84710.16080]);seatc2 = annotation('rectangle',[0.6996 0.7343 0.014290.1024],'FaceColor',[0.84710.16080],'EdgeColor',[0.84710.16080]);seatc3 = annotation('rectangle',[0.7964 0.7333 0.014290.1024],'FaceColor',[0.84710.16080],'EdgeColor',[0.84710.16080]);word4=annotation('textbox','Position',[0.64290.6286 0.10180.07619],'EdgeColor','none','FitHeightToText','off,'String',{'SeatC'});watierplace=annotation('rectangle',[0.82140.1762 0.019640.3857],'FaceColor',[011],'EdgeColor',[011]);word5=annotation('textbox','Position',[0.71070.2 0.11790.0619],'EdgeColor','none','FitHeightToText','off,'String',{等待席'});functionleaveflash(xx,yy,zz)globalman1mana1mana2mana3mana4mana5globalmanb1manb2manb3manb4manb5globalmanc1manc2manc3manc4manc5globalmand1mand2mand3mand4mand5globalmanelmane2mane3mane4mane5switchxxcase1forii=10:-1:1face(:,1)=man1(:,1)+yy*(ii-1);face(:,3)=man1(:,3)+yy*(ii-1);face(:,2)=man1(:,2)+zz*(ii-1);face(:,4)=man1(:,4)+zz*(ii-1);face(1,3)=0.0982;face(1,4)=0.1381;set(mana1,'position',[face(1,1)face(1,2) face(1,3)face(1,4)]);set(mana2,'X',[face(2,1);face(2,3)],'Y',[face(2,2);face(2,4)]);set(mana3,'X',[face(3,1);face(3,3)],'Y',[face(3,2);face(3,4)]);set(mana4,'X',[face(4,1);face(4,3)],'Y',[face(4,2);face(4,4)]);set(mana5,'X',[face(5,1);face(5,3)],'Y',[face(5,2);face(5,4)]);MM(ii)=getframe;endcase2forii=10:-1:1face(:,1)=man1(:,1)+yy*(ii-1);face(:,3)=man1(:,3)+yy*(ii-1);face(:,2)=man1(:,2)+zz*(ii-1);face(:,4)=man1(:,4)+zz*(ii-1);face(1,3)=0.0982;face(1,4)=0.1381;set(manb1,'position',[face(1,1)face(1,2) face(1,3)face(1,4)]);set(manb2,'X',[face(2,1);face(2,3)],'Y',[face(2,2);face(2,4)]);set(manb3,'X',[face(3,1);face(3,3)],'Y',[face(3,2);face(3,4)]);set(manb4,'X',[face(4,1);face(4,3)],'Y',[face(4,2);face(4,4)]);set(manb5,'X',[face(5,1);face(5,3)],'Y',[face(5,2);face(5,4)]);MM(ii)=getframe;endcase3forii=10:-1:1face(:,1)=man1(:,1)+yy*(ii-1);face(:,3)=man1(:,3)+yy*(ii-1);face(:,2)=man1(:,2)+zz*(ii-1);face(:,4)=man1(:,4)+zz*(ii-1);face(1,3)=0.0982;face(1,4)=0.1381;set(manc1,'position',[face(1,1)face(1,2) face(1,3)face(1,4)]);set(manc2,'X',[face(2,1);face(2,3)],'Y',[face(2,2);face(2,4)]);set(manc3,'X',[face(3,1);face(3,3)],'Y',[face(3,2);face(3,4)]);set(manc4,'X',[face(4,1);face(4,3)],'Y',[face(4,2);face(4,4)]);set(manc5,'X',[face(5,1);face(5,3)],'Y',[face(5,2);face(5,4)]);MM(ii)=getframe;endcase4forii=10:-1:1face(:,1)=man1(:,1)+yy*(ii-1);face(:,3)=man1(:,3)+yy*(ii-1);face(:,2)=man1(:,2)+zz*(ii-1);face(:,4)=man1(:,4)+zz*(ii-1);face(1,3)=0.0982;face(1,4)=0.1381;set(mand1,'position',[face(1,1)face(1,2) face(1,3)face(1,4)]);set(mand2,'X',[face(2,1);face(2,3)],'Y',[face(2,2);face(2,4)]);set(mand3,'X',[face(3,1);face(3,3)],'Y',[face(3,2);face(3,4)]);set(mand4,'X',[face(4,1);face(4,3)],'Y',[face(4,2);face(4,4)]);set(mand5,'X',[face(5,1);face(5,3)],'Y',[face(5,2);face(5,4)]);MM(ii)=getframe;endcase5forii=10:-1:1face(:,1)=man1(:,1)+yy*(ii-1);face(:,3)=man1(:,3)+yy*(ii-1);-32-face(:,2)=man1(:,2)+zz*(ii-1);face(:,4)=man1(:,4)+zz*(ii-1);face(1,3)=0.0982;face(1,4)=0.1381;set(mane1,'position',[face(1,1)face(1,2) face(1,3)face(1,4)]);set(mane2,'X',[face(2,1);face(2,3)],'Y',[face(2,2);face(2,4)]);set(mane3,'X',[face(3,1);face(3,3)],'Y',[face(3,2);face(3,4)]);set(mane4,'X',[face(4,1);face(4,3)],'Y',[face(4,2);face(4,4)]);set(mane5,'X',[face(5,1);face(5,3)],'Y',[face(5,2);face(5,4)]);MM(ii)=getframe;endcase6forii=10:-1:1face(:,1)=man1(:,1)+yy*(ii-1);face(:,3)=man1(:,3)+yy*(ii-1);face(:,2)=man1(:,2)+zz*(ii-1);face(:,4)=man1(:,4)+zz*(ii-1);face(1,3)=0.0982;face(1,4)=0.1381;set(manf1,'position',[face(1,1)face(1,2) face(1,3)face(1,4)]);set(manf2,'X',[face(2,1);face(2,3)],'Y',[face(2,2);face(2,4)]);set(manf3,'X',[face(3,1);face(3,3)],'Y',[face(3,2);face(3,4)]);set(manf4,'X',[face(4,1);face(4,3)],'Y',[face(4,2);face(4,4)]);set(manf5,'X',[face(5,1);face(5,3)],'Y',[face(5,2);face(5,4)]);MM(ii)=getframe;endcase7forii=10:-1:1face(:,1)=man1(:,1)+yy*(ii-1);face(:,3)=man1(:,3)+yy*(ii-1);face(:,2)=man1(:,2)+zz*(ii-1);face(:,4)=man1(:,4)+zz*(ii-1);face(1,3)=0.0982;face(1,4)=0.1381;set(mang1,'position',[face(1,1)face(1,2) face(1,3)face(1,4)]);set(mang2,'X',[face(2,1);face(2,3)],'Y',[face(2,2);face(2,4)]);set(mang3,'X',[face(3,1);face(3,3)],'Y',[face(3,2);face(3,4)]);set(mang4,'X',[face(4,1);face(4,3)],'Y',[face(4,2);face(4,4)]);set(mang5,'X',[face(5,1);face(5,3)],'Y',[face(5,2);face(5,4)]);MM(ii)=getframe;endendfunction[face1,face2,face3,face4,face5]=personcreat(vect)ifvect==0facel=annotation('el