Теория массового обслуживания. Сивохин А.В - 254 стр.

UptoLike

254
# p[4](t):=rhs(PSiInit[3]):
# plot(p[1](t), t=1..10);
# plot(p[2](t), t=1..10);
# plot(p[3](t), t=1..10);
# plot(p[4](t), t=1..10);
8. Решение системы уравнений для финальных
вероятностей переходов pfi
> SymPi:=solve({-Lambda*pff[0]+Mu*pff[1]=0,
seq(Lambda*pff[k-1]-
(Lambda+k*Mu)*pff[k]+(k+1)*Mu*pff[k+1]=0,
k=1..n-1),
Lambda*pff[n-1]-
(Lambda+n*Mu)*pff[n]+n*Mu*pff[n+1]=0,
seq(Lambda*pff[n+si-
1]-(Lambda+n*Mu)*pff[n+si]+n*Mu*pff[n+si+1]=0,
si=1..m-1),
Lambda*pff[n+m-1]-
n*Mu*pff[n+m]=0,
add(pff[i],
i=0..n+m)=1
},
{seq(pff[i],
i=0..n+m)}
);
SymPi := {pff
3
= 0.09099573309, pff
15
= 5.752928072 10
-9
, pff
0
= 0.2912545977, pff
6
= 0.001705770176,
pff
7
= 0.0004206462376, pff
5
= 0.006917099534, pff
13
= 9.460105077 10
-8
, pff
12
= 3.836184342 10
-7
,
pff
11
= 0.000001555618060, pff
10
= 0.000006308214968, pff
9
= 0.00002558055676,
pff
8
= 0.0001037321790, pff
4
= 0.02804965559, pff
2
= 0.2213990671, pff
14
= 2.332880281 10
-8
,
pff
1
= 0.3591197467}
> SumPi:=add(rhs(SymPi[i]), i=1..n+m+1);
SumPi := 0.9999999999
9. Решение дифференциальных уравнений для
#         p[4](t):=rhs(PSiInit[3]):
#         plot(p[1](t), t=1..10);
#         plot(p[2](t), t=1..10);
#         plot(p[3](t), t=1..10);
#         plot(p[4](t), t=1..10);

8. Решение системы уравнений для финальных
вероятностей переходов pfi

> SymPi:=solve({-Lambda*pff[0]+Mu*pff[1]=0,
                           seq(Lambda*pff[k-1]-
(Lambda+k*Mu)*pff[k]+(k+1)*Mu*pff[k+1]=0,
k=1..n-1),
                           Lambda*pff[n-1]-
(Lambda+n*Mu)*pff[n]+n*Mu*pff[n+1]=0,
                           seq(Lambda*pff[n+si-
1]-(Lambda+n*Mu)*pff[n+si]+n*Mu*pff[n+si+1]=0,
si=1..m-1),
                           Lambda*pff[n+m-1]-
n*Mu*pff[n+m]=0,
                           add(pff[i],
i=0..n+m)=1
                           },
                          {seq(pff[i],
i=0..n+m)}
                         );
                                                           -9
    SymPi := {pff3 = 0.09099573309, pff15 = 5.752928072 10 , pff0 = 0.2912545977, pff6 = 0.001705770176,

                                                                               -8                          -7
         pff7 = 0.0004206462376, pff5 = 0.006917099534, pff13 = 9.460105077 10 , pff12 = 3.836184342 10 ,

         pff11 = 0.000001555618060, pff10 = 0.000006308214968, pff9 = 0.00002558055676,

                                                                                                  -8
         pff8 = 0.0001037321790, pff4 = 0.02804965559, pff2 = 0.2213990671, pff14 = 2.332880281 10 ,

         pff1 = 0.3591197467}


>      SumPi:=add(rhs(SymPi[i]), i=1..n+m+1);
                                            SumPi := 0.9999999999


9. Решение дифференциальных уравнений для
                                                     254