ВУЗ:
Составители:
Рубрика:
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
Страницы
- « первая
- ‹ предыдущая
- …
- 252
- 253
- 254
- 255
- 256
- …
- следующая ›
- последняя »
