xunh2.m

经典的排队论程序

function  [J4]=xunh2(x,av1,ad3,bv,bd,c1,c2,c3,d1,d2,d3,c,s)
% sutiful for the numberof living channels  more than c1  and less than c2
%av1=0.07;%:0.005:0.03;
%avv=av1;
%avv1=av1;
%avv2=av1;
av2=0.25*av1;
av3=0.1*av1;
ad2=0.25*ad3;
ad1=0.1*ad3;
avv=av1+av2+av3;
avv1=av2+av3;
avv2=av3;
add=ad1+ad2+ad3;
add1=ad2+ad3;
add2=ad3;
%add=ad1;
%add1=add;
%add2=add;
    A4=zeros(s+1,(x-1)*(s+1)); 
    I4=zeros(s+1,(c3-1-x)*(s+1));
   for i=1:s+1
     B4(i,1)= -avv1;
   end
     B4=spdiags([B4],[0],s+1,s+1);
     for i=1:s
if  i>=1&i<=d1;  
    C4(i,1)=-add;
elseif i>=d1+1&i<=d2; 
    C4(i,1)=-add1;
elseif  i>=d2+1&i<=d3; 
    C4(i,1)=-add2;    
end 
end
C4=[C4;0];
 for i=1:s
if  i>=1&i<=d1;  
    D4(i,1)=add+avv1+x*bv+min(c-x,i-1)*bd;
elseif i>=d1+1&i<=d2; 
    D4(i,1)=add1+avv1+x*bv+min(c-x,i-1)*bd;
elseif i>=d2+1&i<=d3; 
   D4(i,1)=add2+avv1+x*bv+min(c-x,i-1)*bd;    
end 
end
 D4(s+1,1)=avv1+x*bv+min(c-x,s)*bd;   
  for i=1:s
     E4(i,1)=-min(c-x,i)*bd;
  end
  E4=[0;E4];
  F4=spdiags([C4,D4,E4],[-1,0,1],s+1,s+1);
  for i=1:s+1
  G4(i,1)=-(x+1)*bv;
  end
  H4=spdiags([G4],[0],s+1,s+1);
  J4=[A4,B4,F4,H4,I4];