result.m

通信网络的混合轮询系统仿真和等待时间解析求解程序。每个站有两类数据

Gagan的仿真结果
数据 clear;
r=16;%共有r个  假 站     r/2个工作站
l=zeros(1,r);u=zeros(1,r);
s_load=0.5;
s=0.18618*ones(1,r);%平均服务长度
ss=0.0693259848*ones(1,r);
%ss=0.0661152*ones(1,r);
u=0.021/r*ones(1,r);%
uu=zeros(1,r);
                      rule=ones(1,r); 
 %                                 rule=zeros(1,r);
l(1)=s_load/s(1)*.90909;
for i=2:r,
    l(i)=s_load/s(1)*(1-.90909)/15;
end
门限 
仿真时间1e5
w_avrg =   0.2215    0.1501    0.1550    0.1505    0.1606    0.1529    0.1576    0.1609    0.1483
   0.1488    0.1561    0.1452    0.1482    0.1522    0.1416    0.1555
ans = 0.0128    0.0130    0.0130    0.0130    0.0130    0.0129    0.0129    0.0129    0.0129
    0.0129    0.0129    0.0129    0.0129    0.0128    0.0128    0.0128
解析
 w 0.2220    0.1555    0.1554    0.1552    0.1551    0.1549    0.1547    0.1546    0.1544
    0.1543    0.1541    0.1539    0.1538    0.1536    0.1534    0.1533
 ddn  0.0128    0.0130    0.0130    0.0130    0.0130    0.0130    0.0130    0.0129    0.0129
   0.0129    0.0129    0.0129    0.0129    0.0129    0.0128    0.0128
   
穷尽
仿真1e5
w_avrg =0.1814    0.3458    0.3235    0.3369    0.3531    0.3386    0.3431    0.3358    0.3562
 0.3702    0.2976    0.3509    0.3641    0.3274    0.3376    0.3666
 ans = 0.0013    0.0277    0.0278    0.0278    0.0279    0.0279    0.0280    0.0281    0.0281
 0.0282    0.0283    0.0283    0.0284    0.0285    0.0285    0.0286
 解析
0.1841    0.3372    0.3379    0.3387    0.3394    0.3402    0.3409    0.3417    0.3424
   0.3432    0.3440    0.3448    0.3456    0.3464    0.3472    0.3480
ddn = 0.0013    0.0282    0.0283    0.0283    0.0284    0.0284    0.0285    0.0286    0.0286
  0.0287    0.0288    0.0288    0.0289    0.0290    0.0290    0.0291

  
  
%小工况  debug   不对称  用于调试门穷
clear;
r=10;%共有r个  假 站     r/2个工作站
s=1*ones(1,r);%平均服务长度
ss=2*ones(1,r);
s_load=0.5;
percent=0.5;
for i=1:r/2,
   if i==1,
       l(2*i-1)=s_load/s(1)*percent/2;        l(2*i)=l(2*i-1);
   else 
       l(2*i-1)=s_load/s(1)*(1-percent)/(r/2-1)/2;  l(2*i)=l(2*i-1);
   end
       rule(2*i-1)=1;                    rule(2*i)=0;
       u(2*i-1)=0;                          u(2*i)=0.1;
end


 rule(2*i-1)=1;                    rule(2*i)=0;
解析   x1 = 1.3045    1.4600    1.4612    1.4626    1.4643
       x2 = 1.6772    1.5542    1.5554    1.5569    1.5587
  仿真 
  
 rule(2*i-1)=0;                    rule(2*i)=1;
解析   x1 = 
       x2 =
仿真   1.6865    1.3131    1.5484    1.4536    1.5495    1.4549    1.5468    1.4517    1.5470    1.4537
  
 


table II of 文章   (不对称负载 )
  clear;
r=10;%共有r个  假 站     r/2个工作站
s=1*ones(1,r);%平均服务长度
ss=2*ones(1,r);
s_load=0.5;
for i=1:r/2,
   if i==1,
       l(2*i-1)=s_load/s(1)*.5/2;        l(2*i)=l(2*i-1);
   else 
       l(2*i-1)=s_load/s(1)*(1-.5)/4/2;  l(2*i)=l(2*i-1);
   end
       rule(2*i-1)=1;                    rule(2*i)=0;
       u(2*i-1)=0;                          u(2*i)=0.1;
end
仿真 1.3030    1.6759    1.4641    1.5570    1.4668    1.5546    1.4613    1.5559    1.4622    1.5532
解析w =  1.3045    1.6772    1.4600    1.5542    1.4612    1.5554    1.4626    1.5569    1.4643    1.5587

不对称 
rule(2*i-1)=1;                    rule(2*i)=1;
 解析 1.3896    1.3801    1.5296    1.5314    1.5333    1.5354    1.5377    1.5402    1.5429    1.5459
 仿真 1.3894    1.3816    1.5316    1.5304    1.5369    1.5428    1.5396    1.5387    1.5458    1.5470

 rule(2*i-1)=0;                    rule(2*i)=0;
 解析  1.6006    1.6063    1.4794    1.4782    1.4770    1.4757    1.4742    1.4727    1.4710    1.4692
 仿真  1.6020    1.5998    1.4687    1.5154    1.4820    1.4813    1.4795    1.4856    1.4706    1.4665 1e6秒

   1e7秒 1.5944    1.6042    1.4838    1.4697    1.4740    1.4782    1.4703    1.4715    1.4644    1.4619
  
   1e7秒  1.6050    1.6137    1.4803    1.4808    1.4844    1.4799    1.4837    1.4824    1.4748    1.4709

    
    
 table II of 文章   （ 对称  ）
  clear;
r=10;%共有r个  假 站     r/2个工作站
s=1*ones(1,r);%平均服务长度
ss=2*ones(1,r);
s_load=0.5;
for i=1:r/2,
       l(2*i-1)=s_load/s(1)/r;  l(2*i)=l(2*i-1);
       rule(2*i-1)=1;                    rule(2*i)=0;
       u(2*i-1)=0;                          u(2*i)=0.1;
end
解析   1.4250    1.5750    1.4250    1.5750    1.4250    1.5750    1.4250    1.5750    1.4250    1.5750
仿真   1.4228    1.5701    1.4223    1.5764    1.4283    1.5809    1.4291    1.5707    1.4189    1.5731


全穷
rule(2*i-1)=1;                    rule(2*i)=1;
解析  1.4750    1.4750    1.4750    1.4750    1.4750    1.4750    1.4750    1.4750    1.4750    1.4750
仿真  1.4833    1.4777    1.4843    1.4743    1.4742    1.4789    1.4775    1.4807    1.4726    1.4826


全门
rule(2*i-1)=0;                    rule(2*i)=0;
解析  1.5250    1.5250    1.5250    1.5250    1.5250    1.5250    1.5250    1.5250    1.5250    1.5250
仿真  1.5258    1.5295    1.5220    1.5245    1.5292    1.5257    1.5270    1.5212    1.5260    1.5316

 table II of 文章   24个工作站（ 不  对称  ）
 clear;
r=24;%共有r个  假 站     r/2个工作站
s=1*ones(1,r);%平均服务长度
ss=2*ones(1,r);
s_load=0.5;
percent=0.5;
for i=1:r/2,
   if i==1,
       l(2*i-1)=s_load/s(1)*percent/2;        l(2*i)=l(2*i-1);
   else 
       l(2*i-1)=s_load/s(1)*(1-percent)/(r/2-1)/2;  l(2*i)=l(2*i-1);
   end
       rule(2*i-1)=1;                    rule(2*i)=0;
       u(2*i-1)=0;                          u(2*i)=0.1;
end