用Matlab模拟一个缓冲区大小为10的随机排队系统

function p=queue(mu,lamda)
customer=zeros(1000,6);                         %每一行代表是一位顾客,每一行的每一列所代表的意思如下所示
NO=1;                                           %/*编号*/
intime=2;                                       %/*进入队列时间*/
intertime=3;                                    %/*他的下一个人到来的时间间隔*/
starttime=4;                                    %/*开始接受服务时间*/
durtime=5;                                      %/*服务所需时间*/
leavetime=6;                                    %/*离开队列的时间*/
flag=1;                                         %初始化flag=1,当flag=1时,采用泊松和指数分布,当flag=0时,采用平均和瑞利分布

%*********************************InitQueue();%初始化队列***********************

queue_customernum=0;                            %/*队列中接受服务的顾客数*/
curtime=1;                                      %/*当前时间*/
first=0;                                        %/*当前接受服务的顾客的编号*/
last=0;                                         %/*当前队列最后一位顾客的编号*/
totalnum=0;                                     %/*进入队列的总的顾客数*/
cusnum=zeros(1,100);                            %/*存放观察到的队列中接受服务的顾客数*/
p=zeros(1,11);                                  %/*存放观察到的队列中接受服务的某个顾客数的出现的次数*/

%********************************customer_in();%顾客进入队列*******************

if queue_customernum~=10
    totalnum=totalnum+1;
    queue_customernum=queue_customernum+1;
    customer(totalnum,NO)=totalnum;
    customer(totalnum,intime)=curtime;                                  %/*记录顾客进入时间*/
    
    if flag==1
        customer(totalnum,durtime)=poissrnd(mu,1,1);                       %/*产生随机数记录顾客服务所需时间*/ 泊松分布
        customer(totalnum,intertime)=round(exprnd(lamda,1,1)*10);          %/*此顾客的下一个顾客来的时间间隔*/ 指数分布
    else
        customer(totalnum,durtime)=unidrnd(mu,1,1);                         %/*产生随机数记录顾客服务所需时间*/ 均匀分布
        customer(totalnum,intertime)=round(raylrnd(lamda,1,1)*10);          %/*此顾客的下一个顾客来的时间间隔*/ 瑞利分布
    end
    last=last+1;
else
    customer(totalnum,intertime)=customer(totalnum,intertime)+1;
end

%*****************************customer_serve();%顾客接受服务********************

first=first+1;
customer(first,starttime)=curtime;                                      %/*得出开始服务时间和离开时间*/
customer(first,leavetime)=customer(first,durtime)+curtime;

%**************************************************************************

while curtime<=1000                                                     %/*当前时间属于营业时间,允许顾客进入*/
     curtime=curtime+1;
     if customer(first,leavetime)==curtime
        %customer_leave();%顾客离开队列
        queue_customernum=queue_customernum-1; 
        
        %customer_serve();%顾客接受服务
        first=first+1;
        customer(first,starttime)=curtime;                              %/*得出开始服务时间和离开时间*/
        customer(first,leavetime)=customer(first,durtime)+curtime;
        
     end
     
     if (customer(last,intime)+customer(last,intertime))==curtime
        %customer_in(); %顾客进入队列                                               %/*判断是否有人符合要进的条件*/
        if queue_customernum~=10
                totalnum=totalnum+1;
                queue_customernum=queue_customernum+1;
                customer(totalnum,NO)=totalnum;
                customer(totalnum,intime)=curtime;                              %/*记录顾客进入时间*/
                if flag==1
                    customer(totalnum,durtime)=poissrnd(mu,1,1);                   %/*产生随机数记录顾客服务所需时间*/ 泊松分布
                    customer(totalnum,intertime)=round(exprnd(lamda,1,1)*10);      %/*此顾客的下一个顾客来的时间间隔*/ 指数分布
                else
                    customer(totalnum,durtime)=unidrnd(mu,1,1);                     %/*产生随机数记录顾客服务所需时间*/ 均匀分布
                    customer(totalnum,intertime)=round(raylrnd(lamda,1,1)*10);      %/*此顾客的下一个顾客来的时间间隔*/ 瑞利分布
                end
                last=last+1;
        else
            customer(totalnum,intertime)=customer(totalnum,intertime)+1;
        end
        
     end
     
     if rem(curtime,10)==0
         cusnum(1,curtime/10)=queue_customernum;
     end
end

%************************probability();%计算队列中顾客数为(0,1,。。。,10)时的概率**********

  for i=1:1:100
    if cusnum(1,i)==0
        p(1,1)=p(1,1)+1;
    elseif cusnum(1,i)==1
        p(1,2)=p(1,2)+1;    
    elseif cusnum(1,i)==2
        p(1,3)=p(1,3)+1;
    elseif cusnum(1,i)==3
        p(1,4)=p(1,4)+1;
    elseif cusnum(1,i)==4
        p(1,5)=p(1,5)+1;
    elseif cusnum(1,i)==5
        p(1,6)=p(1,6)+1;
    elseif cusnum(1,i)==6
        p(1,7)=p(1,7)+1;
    elseif cusnum(1,i)==7
        p(1,8)=p(1,8)+1;
    elseif cusnum(1,i)==8
        p(1,9)=p(1,9)+1;
    elseif cusnum(1,i)==9
        p(1,10)=p(1,10)+1;
    elseif cusnum(1,i)==10
        p(1,11)=p(1,11)+1;
    end
  end
  
%***************************list();%打印结果****************************
  for j=1:1:11
    p(1,j)=p(1,j)/100;
  end
  用Matlab模拟一个缓冲区大小为10的随机排队系统_第1张图片
  D118

 

你可能感兴趣的:(MATLAB,其他,matlab,算法,随机排队系统)