shishi_dongtai.m

经典的排队论程序

clear all
%av1=0.5:0.01:0.6;
%ad3=1;
%av1=60:20:240;
%ad3=40;%:20:120;
av1=0.01:0.001:0.02;
ad3=0.1:0.01:0.2;
[M,N]=meshgrid(av1,ad3); 
%bv=0.03;
%bd=10;
bv=1/180;
bd=1;
c1=15;
c2=17;
c3=19;
d1=10;
d2=15;
d3=19;
c=25;
s=19;
 for i=1:size(M,1)
     for j=1:size(N,2)
    [Pv1,Pv2,Pv3,Pd1,Pd2,Pd3,PL]=sanwei_vip1(M(i,j),N(i,j),bv,bd,c1,c2,c3,d1,d2,d3,c,s);
    P11(i,j)=Pv1;
    P22(i,j)=Pv2;
    P33(i,j)=Pv3;
    Pd11(i,j)=Pd1;
    Pd22(i,j)=Pd2;
    Pd33(i,j)=Pd3;
     ll(i,j)=PL;
     end
 end
 for i=1:size(M,1)
     for j=1:size(N,2)
    [Pv1,Pv2,Pv3,Pd1,Pd2,Pd3,PL]=sanwei_vip1(M(i,j),N(i,j),bv,bd,c1,c2,c3,d1,d2,d3,c,s);
    P1(i,j)=Pv1;
    P2(i,j)=Pv2;
    P3(i,j)=Pv3;
    PD1(i,j)=Pd1;
    PD2(i,j)=Pd2;
    PD3(i,j)=Pd3;
    l(i,j)=PL;
    %G1(i)=G;
    %[Pn,Ph,Pd0,G0,PL0,T0]=normal(a1n(i),0.18);
    %P5(i)=Pn;
    %P6(i)=Ph;
    %G2(i)=G0;
if mod(i,3)==0
cc3=c3;
cc2=c2;
cc1=c1;
DD3=d3;
DD2=d2;
DD1=d1;
%u6=4e-40;d6=2e-40;
%u5=2e-29;d5=1e-29;
%u4=3e-27;d4=1e-27;
uv3=3e-9;dv3=1e-9;
%uv3=1e-9;dv3=8e-10;
uv2=7e-7;dv2=2e-7;
%uv2=5e-7;dv2=4e-7;
uv1=12e-5;dv1=2e-5;
%uv1=10e-4;dv1=3e-5;
ud3=1e-10;dd3=1e-11;
ud2=3e-4;dd2=1e-6;
ud1=0.05;dd1=0.025;
if P3(i)>uv3&c3>cc3&c2<=19&c2>=15
    c2=cc2;
    elseif P3(i)>uv3&c3==cc3&c2<=19&c2>15
    c2=cc2-1;
elseif P3(i)>uv3&c3<cc3&c2<=19&c2>16
    c2=cc2-2;
elseif P3(i)<dv3&c3>cc3&c2<=17&c2>=15
    c2=cc2+2; 
    c3=cc3+2;
elseif P3(i)<dv3&c3==cc3&c2<=18&c2>=15
    c2=cc2+1; 
    c3=cc3+1;
elseif P3(i)<dv3&c3<cc3&c2<=19&c2>=15
    c2=cc2; 
end
if P2(i)>uv2&c2>cc2&c1<=17&c1>=13
    c1=cc1;
elseif P2(i)>uv2&c2==cc2&c1<=17&c1>13
    c1=cc1-1;
elseif P2(i)>uv2&c2<cc2&c1<=17&c1>14
    c1=cc1-2;
elseif P2(i)<dv2&c2>cc2&c1<=15&c1>=13
    c1=cc1+2; 
elseif P2(i)<dv2&c2==cc2&c1<=16&c1>=13
    c1=cc1+1; 
elseif P2(i)<dv2&c2<cc2&c1<=17&c1>=13
    c1=cc1; 
end

 if PD3(i)>ud3&d3>DD3&d2<=17&d2>=13
    d2=DD2;
    elseif PD3(i)>ud3&d3==dd3&d2<=17&d2>13
    d2=DD2-1;
elseif PD3(i)>ud3&d3<DD3&d2<=17&d2>14
    d2=DD2-2;
elseif PD3(i)<dd3&d3>DD3&d2<=15&d2>=13
    d2=DD2+2;
    %d3=DD3+5;
    %d3=DD3+10;
    %s=s+10;
elseif PD3(i)<dd3&d3==DD3&d2<=16&d2>=13
    d2=DD2+1; 
    %d3=DD3+11;
    %s=s+11;
elseif PD3(i)<dd3&d3<DD3&d2<=17&d2>=13
    d2=DD2; 
end
if PD2(i)>ud2&d2>DD2&d1<=12&d1>=8
    d1=DD1;
elseif PD2(i)>ud2&d2==dd2&d1<=12&d1>8
    d1=dd1-1;
elseif PD2(i)>ud2&d2<DD2&d1<=12&d1>9
    d1=DD1-2;
elseif PD2(i)<DD2&d2>DD2&d1<=10&d1>=8
    d1=DD1+2; 
elseif PD2(i)<DD2&d2==DD2&d1<=11&d1>=8
    d1=DD1+1; 
elseif PD2(i)<DD2&d2<DD2&d1<=12&d1>=8
    d1=DD1; 
end
 end
     end
 end
 %figure(1),plot(a1n,G1,'r*-');
 %hold on;
% plot(a1n,G2,'gd-');
 %   grid on;
%   xlabel('\it\lambda_{v}\rightarrow(calls/s)');
%   ylabel('\itGos'); 
%   title('\fontname{楷书}系统的服务质量');
%   set(gca,'Xtick',linspace(0.07,0.11,5));
av1=0.01:0.001:0.02;
ad3=0.1:0.01:0.2;
[M,N]=meshgrid(av1,ad3);
figure(1)
surf(M,N,P11);
figure(2)
surf(M,N,P1);
 %plot(av1,P11,'gd-',av1,P1);
  %plot(av1,P1,'gd-');
 %hold on;
 figure(3)
surf(M,N,P22);
figure(4)
surf(M,N,P2);
 %plot(av1,P22,'r*-',av1,P2);
 %plot(av1,P2,'r*-');
 %hold on;
 figure(5)
surf(M,N,P33);
figure(6)
surf(M,N,P3);
 %plot(av1,P33,'yo-',av1,2*P3);
  %plot(av1,P3,'yo-');
 grid on;
 % xlabel('\it\lambda_{v}\rightarrow(calls/s)');
 %  ylabel('\itPn'); 
%   title('\fontname{楷书}新呼叫的阻塞率');
   %set(gca,'Xtick',linspace(0.07,0.11,5));
  figure(7)
  %plot(av1,P11,'gd-');
surf(M,N,Pd11);
   figure(8)
surf(M,N,PD1);
  %plot(av1,Pd11,'gd-',av1,PD1);
  %hold on;
  figure(9)
  %plot(av1,P22,'r*-');
surf(M,N,Pd22);
  figure(10)
surf(M,N,PD2);
  %plot(av1,Pd22,'r*-',av1,PD2);
  %hold on;
  figure(11)
  %plot(av1,P33,'yo-');
surf(M,N,Pd33);
  figure(12)
surf(M,N,PD3/3.5);
  %plot(av1,Pd33,'yo-',av1,PD3/3.5);
  %grid on;
  figure(13)
surf(M,N,l);
  %plot(av1,l,'gd-',av1,ll);
  %plot(av1,l,'gd-')
  %grid on
 % xlabel('\it\lambda_{v}\rightarrow(calls/s)');
 %  ylabel('\itPh'); 
 %  title('\fontname{楷书}切换呼叫的阻塞率');
 %  set(gca,'Xtick',linspace(0.07,0.11,5));