clear;

 clc;

 Alpha=0.1; %置信水平

 Gama=0.1; %相对精度

 Beta=0.1;

 Lambda=0.2; %到达率 Lambda

 Mu=0.25; %服务率 Mu

 time=; %单回模拟次数

 %序贯法实现

 All_vector=Func( Lambda,Mu );%函数返回的向量

 for i=:time

 All_vector=[All_vector;Func( Lambda,Mu )];

 end

 vect=sum(All_vector)/time;%各列求和取平均

 %下面系统中平均等待队长作为每次模拟的 X，用来评估运行结果

 S=sum((All_vector()-vect()).*(All_vector()-vect()))/(ti

 me-);%后面假设 S 不变

 Betan=sqrt(S/time)*tinv(-Alpha/,time-);

 Gaman=Betan/vect();

 while(Gaman>=Gama) %Betan>=Beta||

 time=time+;

 All_vector=[All_vector;Func( Lambda,Mu )];

 vect=sum(All_vector)/time;

 S=sum((All_vector()-vect()).*(All_vector()-vect()))/(ti

 me-);%后面假设 S 不变

 Betan=sqrt(S/time)*tinv(-Alpha/,time-);

 Gaman=Betan/vect();

 end

 time

 disp(['系统中平均等待队长的置信区间下界

 =',num2str(vect(3)-Betan)]);

 disp(['系统中平均等待队长的置信区间上界

 =',num2str(vect(3)+Betan)]);

 function [ vector ] = Func( Lambda,Mu )

 %单次的排队模拟，样本数 CusTotal

 CusTotal=; %仿真顾客总数%=input('请输入仿真顾客

 总数 CusTotal=');

 Cus_Arrive=zeros(,CusTotal);%到达时刻

 Cus_Leave=zeros(,CusTotal);%离开时刻

 IntervaCus_Arrive=-log(rand(,CusTotal))/Lambda;%到达时间间

 隔

 Cus_Arrive=cumsum(IntervaCus_Arrive); %每列

 累加，形成到达初始时刻；如果只有一行，该行向后叠加

 Interval_Serve=-log(rand(,CusTotal))/Mu; %服务时间间隔

 %为事件调度法做准备

 Cus_Leave()=Cus_Arrive()+Interval_Serve();%顾客离开时间

 for i=:CusTotal

 if Cus_Leave(i-)<Cus_Arrive(i)

 Cus_Leave(i)=Cus_Arrive(i)+Interval_Serve(i);

 else

 Cus_Leave(i)=Cus_Leave(i-)+Interval_Serve(i);

 end

 end

 Cus_Wait=Cus_Leave-Cus_Arrive; %各顾客在系统中的等待时间

 %mean：如果是 n*m 的矩阵，mean 对各列分别进行求平均；当 n= 时，

 对一行求平均

 Cus_Wait_avg=mean(Cus_Wait); %平均等待时间

 Cus_Queue=Cus_Wait-Interval_Serve;%各顾客在系统中的排队时间

 Cus_Queue_avg=mean(Cus_Queue); %平均排队时间

 %TimePoint 系统调度时间

 TimePoint=[Cus_Arrive,Cus_Leave];%系统中顾客数随时间的变化

 TimePoint=sort(TimePoint);

 CusNum=zeros(size(TimePoint));

 temp=; %指向事件表

 CusNum()=;

 %统计 dt 时间内的排队人数——事件调度法

 %截止到 i，还有多少个事件

 for i=:length(TimePoint)

 %后一事件的结束时间晚于前一事件的结束时间

 %所以只要第 temp- 个事件没有结束， 后面的事件到了发生时间

 (事件仿真钟<=系统仿真钟)，都要加入 CusNum 中计数

 if

 (temp<=length(Cus_Arrive))&&(TimePoint(i)==Cus_Arrive(temp)

 )

 CusNum(i)=CusNum(i-)+;

 temp=temp+;

 else

 CusNum(i)=CusNum(i-)-;

 end

 end

 %系统中平均顾客数计算

 Time_interval=zeros(size(TimePoint));

 Time_interval()=Cus_Arrive();

 for i=:length(TimePoint)

 Time_interval(i)=TimePoint(i)-TimePoint(i-);

 end

 %这里画一下图即可.i 时间段*(i-)时刻统计的事件数量。类似于向

 后积分

 CusNum_fromStart=[ CusNum];

 Cus_Wait_CusNum_avg=sum(CusNum_fromStart.*[Time_interval

 ] )/TimePoint(end);

 %系统平均等待队长

 QueLength=zeros(size(CusNum));

 for i=:length(CusNum)

 if CusNum(i)>=

 %还有等待事件（满足事件仿真钟<=系统仿真钟，但没有发

 生）

 QueLength(i)=CusNum(i)-;

 else

 QueLength(i)=;

 end

 end

 Cus_Wait_Queue_avg=sum([ QueLength].*[Time_interval

 ] )/TimePoint(end);

 %仿真值与理论值比较

 row=Lambda/Mu;

 QueLength_avg=row*row/(-row);%Q

 CusNum_avg=row/(-row);%L

 Queue_avg=QueLength_avg/Lambda;%d

 Wait_avg=CusNum_avg/Lambda;%w

 %返回的向量

 vector=[Cus_Queue_avg,Cus_Wait_avg,Cus_Wait_Queue_avg,Cus_W

 ait_CusNum_avg];

 end