04仿真模型设计与实现.ppt

* * * function sqssMoreReport() %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % calculate the systems performance when 100 customers served % edit by jiannywang@163.com at 2014/2/10 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% sqss=initData(); runTimes=0; %preventing the unlimited while loop while(sqss.servedCustomerssqss.totalCustomers) sqss=timeUpdate(sqss); if sqss.nextEvent==1 sqss=arriveProcess(sqss); else sqss=departProcess(sqss); end runTimes=runTimes+1; if runTimes1000 break; end end report(sqss); end function sqss=initData() sqss=struct(... simuClock,0,totalCustomers,100,nextEvent,0,... arriveMean,5,serveTimeMin,2,serveTimeMax,8,... event1Time,0,event2Time,inf,... queueLength,0,serverStatus,0,... comingCustomers,0,servedCustomers,0,... eachArriveTime,[],eachLeaveTime,[],... eachServeTime,[],eachServeStartTime,[],... eachQueueLength,[],eachQueueStartTime,[]); sqss.event1Time=arriveInterval(sqss.arriveMean); end function sOut=timeUpdate(sIn) sOut=sIn; if sIn.event2TimesIn.event1Time sOut.nextEvent=2; sOut.simuClock=sIn.event2Time; else sOut.nextEvent=1; sOut.simuClock=sIn.event1Time; end end function sOut=arriveProcess(sIn) sOut=sIn; sOut.event1Time=sIn.simuClock+arriveInterval(sIn.arriveMean); sOingCustomers=sIingCustomers+1; % record this customer arrive time sOut.eachArriveTime=[sOut.eachArriveTime;sIn.simuClock]; if sIn.serverStatus==1 sOut.queueLength=sIn.queueLength+1; else sOut.serverStatus=1; midServeTime=serveTime(sIn); sOut.event2Time=sIn.simuClock+midServeTime; % update server time data sOut.eachServeTime=[sOut.eachServeTime;midServeTime]; sOut.eachServeStartTime=[sOut.eachServeStartTime;sIn.simuClock]; end % update