第PAGE页码2页/总页数NUMPAGES总页数2页排队论(matlab代码可运行)排队论(matlab代码可运行)篇一%/M/PH/1(k);%p=stationary_prob(4,beta,S,order,lambda);function[p,p_minus,p_add,mean,time]=stationary_prob(k,beta,S,lambda)order=length(beta);e=ones(order,1);R=lambda*inv(lambda*eye(size(S))lambda*e*betaS);R_k=eye(size(S));fori=1:kR_k=R_k+R^i;endp0=inv(beta*(R_klambda*R^k*inv(S))*e);p_t=p0;fori=1:kp_t=[p_t,p0*beta*R^i];endp_t=[p_t,p0*beta*(R^k)*(lambda*inv(S))];%thequeuelengthdistributioninanytime%addeveryphaseinthesamelevelp=p0;fori=1:k+1sum=0;forj=1:ordersum=sum+p_t(1+(i1)*order+j);endp=[p,sum];endp_minus=p;%thequeuelengthdistributionatdeparture%p_add=p_minus(i)/(1p_minus(k+1))p_add=[0];fori=0:kp_add=[p(ki+1)/(1p(k+2)),p_add];end%themeanqueuelengthmean=0;fori=1:kmean=mean+i*p(i+1);end%themeanwaitingtimetime=mean/lambda/(1p(k+2));%%application%%M/E2/1(k)lambda=2;mu=4;beta=[1,0];S=[mu,mu;0,mu];[p,p_minus,p_add,mean,time]=stationary_prob(4,beta,S,lambda);pp_minusp_addmeantime%M/H2/1(k)lambda=2;mu=2;beta=[0.5,0.5];S=[3,0;0,1.5];[p,p_minus,p_add,mean,time]=stationary_prob(4,beta,S,lambda);pp_minusp_addmeantime%M/M/1(k)lambda=2;mu=2;beta=[1];S=[mu];[p,p_minus,p_add,mean,time]=stationary_prob(4,beta,S,lambda);pp_minusp_addmeantime