mm1但服务台排队模型

Simtotal=input('请输入仿真顾客总数SimTotal=');%仿真顾客总数； Lambda=14;%到达率Lambda; Mu=0.078;%服务率Mu; t_Arrive=zeros(1,SimTotal); t_Leave=zeros(1,SimToal); ArriveNum=zeros(1,simToal); LeaveNum=zeros(1,SimToal); 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)<t_Arrive(i) t_Leave(i)=t_Arrive(i)+Interval_Serve(i); else t_Leave(i)=t_Leave(i-1)+Interval_Serve(i); end LeaveNum(i)=i; end t_Wait=t_Leave-t_Arrive; %各顾客在系统中的等待时间 t_Wait_avg=mean(t_Wait); t_Queue=t_Wait-Interval_Serve;%各顾客在系统中的排队时间 t_Queue_avg=mean(t_Queue); Timepoint=[t_Arrive,t_Leave];%系统中顾客数随时间的变化 Timepoint=sort(Timepoint); ArriveFlag=zeros(size(Timepoint));%到达时间标志 CusNum=zeros(size(Timepoint)); temp=2; CusNum(1)=1; for i=2:length(Timepoint) if (temp<=length(t_Arrive))&&(Timepoint(i)==t_Arrive(temp)) CusNum(i)=CusNum(i-1)+1; temp=temp+1; ArriveFlag(i)=1; else CusNum(i)=CusNum(i-1)-1; end end %系统中平均顾客数计算 Time_interval=zeros(size(Timepoint)); Time_interval(1)=t_Arrive(1); for i=2:length(Timepoint) Time_interval(i)=Timepoint(i)-Timepoint(i-1); end CusNum_fromStart=[0 CusNum]; CusNum_avg=sum(CusNum_fromStart.*[Time_interval 0] )/Timepoint(end); QueLength=zeros(size(CusNum)); for i=1:length(CusNum) if CusNum(i)>=2 QueLength(i)=CusNum(i)-1; else QueLength(i)=0; end end QueLength_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)]);