fork download
copy 
  1. const 
  2. lambda=3.333;                                    
  3. mu=0.57;                                              
  4. eta=0.667; 
  5. time_t=1000000;                                    
  6. n_kanalov=4;                                         
  7. l=2;                                                         
  8. var 
  9. n_svob,p:integer;                                    
  10. vypolneno:real;                                       
  11. kolvo,otkaz:longint;                               
  12. i,ok:integer; 
  13. t:longint;                                                 
  14. t_z0,t_z,t_obs:real;                                 
  15. n:array[1..n_kanalov] of real;                 
  16. t_zan:array[1..n_kanalov] of real;          
  17.  
  18. n_vyp:array[1..n_kanalov] of integer;                                                                   
  19. n_l_z:array[1..n_kanalov] of integer; 
  20. Pobs,Pzk,Ksr:real;                                         
  21. tpk:real;                                                   
  22. n_zero:array[1..l] of integer;                   
  23. begin 
  24. randomize; 
  25. t_z0:=0; 
  26. kolvo:=0; 
  27. otkaz:=0; 
  28. vypolneno:=0; 
  29. for i:=0 to n_kanalov do                                    
  30.     begin 
  31.         n_vyp[i]:=0; 
  32.         t_zan[i]:=0; 
  33.     end; 
  34. for t:=0 to time_t do                                           
  35.     begin 
  36.       for i:=1 to n_kanalov do                               
  37.          begin 
  38.              if ((t>=n[i])AND(n[i]<>0)) then           
  39.                begin                                                    
  40.                      n[i]:=0; 
  41.                      vypolneno:=vypolneno+1/n_l_z[i]; 
  42.                      n_l_z[i]:=0; 
  43.                 end; 
  44.          end; 
  45.       if (t>=t_z0) then                                            
  46.          begin                                                           
  47.              t_z:=-100/lambda*ln((random(1000)+1)/1000); 
  48.              t_z0:=t_z0+t_z; 
  49.              inc(kolvo);                                               
  50.              n_svob:=0; 
  51.              for i:=1 to l do n_zero[i]:=0;                    
  52.              p:=1; 
  53.              for i:=1 to n_kanalov do                           
  54.                  begin 
  55.                     if (n[i]=0) then 
  56.                        begin 
  57.                          inc(n_svob); 
  58.                          n_zero[p]:=i; 
  59.                          inc(p); 
  60.                        end; 
  61.                  end; 
  62.              if (n_svob>=l) then begin                                                                                                              
  63.  
  64.               t_obs:=-100/(l*mu+eta)*ln((random(1000)+1)/1000); 
  65.                 for i:=1 to l do                                            
  66.                   begin 
  67.                           n[n_zero[i]]:=t+t_obs;                     
  68.                            n_l_z[n_zero[i]]:=l; 
  69.          t_zan[n_zero[i]]:=t_zan[n_zero[i]]+t_obs;              
  70.                           n_vyp[n_zero[i]]:=n_vyp[n_zero[i]]+1;    
  71.                           ok:=1;                                                         
  72.                   end; 
  73.               end 
  74.              else 
  75.              if (n_svob>0)AND(n_svob<l) then begin    
  76.  
  77.              t_obs:=-100/(n_svob*mu+eta)*ln((random(1000)+1)/1000); 
  78.  
  79.                 for i:=1 to n_svob do                                 
  80.  
  81.                   begin 
  82.                         n[n_zero[i]]:=t+t_obs;                     
  83.                         n_l_z[n_zero[i]]:=n_svob; 
  84. 						t_zan[n_zero[i]]:=t_zan[n_zero[i]]+t_obs;
  85. 						n_vyp[n_zero[i]]:=n_vyp[n_zero[i]]+1;    
  86.                         ok:=1;                                                         
  87.                   end; 
  88.               end 
  89.               else ok:=0;                                                              
  90.               if (ok<>1) then inc(otkaz);                                     
  91.             end;                                                                           
  92.  
  93.     end;                                                                          
  94.  
  95.          {Вычисление интересующих параметров по результатам моделирования} 
  96. pzk:=0; 
  97. for i:=1 to n_kanalov do pzk:=pzk+t_zan[i]/time_t; 
  98. 	pzk:=pzk/n_kanalov;                                                 
  99. Pobs:=1-(otkaz)/kolvo;  
  100. Ksr:=Pzk*n_kanalov;
  101.  
  102. writeln('Вероятность обслуживания ',Pobs:5:3);             
  103. writeln('Вероятность занятости канала ',pzk:5:3); 
  104. writeln('Ksr  ',Ksr:5:3); 
  105. end. 
  106.  
Success #stdin #stdout 0.02s 5288KB
comments (?)
stdin
copy 
Standard input is empty
stdout
copy 
Вероятность обслуживания 0.629
Вероятность занятости канала 0.576
Ksr  2.306
https://ideone.com/1SgK4m
language:
Pascal (fpc 3.0.4)
created:
3 days ago
visibility:
public

Share or Embed source code

Discover > Sphere Engine API

The brand new service which powers Ideone!

Discover > IDE Widget

Widget for compiling and running the source code in a web browser!