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

UptoLike

178
>
5. Решение уравнений Эрланга для вероятностей
состояний pi(t) в общем виде
PSi:=dsolve(ODESystem):
> PSi[1]:
> diff((PSi[1]), t):
> PSi[2]:
> diff((PSi[2]), t):
> PSi[3]:
> diff((PSi[3]), t):
> PSi[4]:
> diff((PSi[4]), t):
> PSi[5]:
> diff((PSi[5]), t):
6. Решение уравнений Эрланга для вероятностей
состояний pi(t) при заданных начальных условиях
> ODESystemInit:={diff(p0(t), t)=-
Lambda*p0(t)+Mu*p1(t),
diff(p1(t),
t)=Lambda*p0(t)-(Lambda+1*Mu)*p1(t)+2*Mu*p2(t),
diff(p2(t),
t)=Lambda*p1(t)-(Lambda+2*Mu)*p2(t)+3*Mu*p3(t),
# diff(p3(t),
t)=Lambda*p2(t)-(Lambda+3*Mu)*p3(t)+4*Mu*p4(t),
diff(p3(t),
t)=Lambda*p2(t)-3*Mu*p3(t),
p0(0)=1.0,
p1(0)=0.0, p2(0)=0.0, p3(0)=0.0
#, p4(0)=0.0
};
ODESystemInit :=
?
?
?
d
dt
p3 t( ) = 2.17 p2 t( ) - 10.77 p3 t(),
d
dt
p0 t( ) = -2.17 p0 t( ) + 3.59 p1 t(),
d
dt
p1 t( ) = 2.17 p0 t( ) - 5.76 p1 t( ) + 7.18 p2 t(),
d
dt
p2 t( ) = 2.17 p1 t( ) - 9.35 p2 t( ) + 10.77 p3 t(), p1 0( ) = 0.,
>

5. Решение уравнений Эрланга для вероятностей
состояний pi(t) в общем виде

PSi:=dsolve(ODESystem):
> PSi[1]:
> diff((PSi[1]), t):
> PSi[2]:
> diff((PSi[2]), t):
> PSi[3]:
> diff((PSi[3]), t):
> PSi[4]:
> diff((PSi[4]), t):
> PSi[5]:
> diff((PSi[5]), t):

6. Решение уравнений Эрланга для вероятностей
состояний pi(t) при заданных начальных условиях

> ODESystemInit:={diff(p0(t), t)=-
Lambda*p0(t)+Mu*p1(t),
                                  diff(p1(t),
t)=Lambda*p0(t)-(Lambda+1*Mu)*p1(t)+2*Mu*p2(t),
                                  diff(p2(t),
t)=Lambda*p1(t)-(Lambda+2*Mu)*p2(t)+3*Mu*p3(t),
#                                  diff(p3(t),
t)=Lambda*p2(t)-(Lambda+3*Mu)*p3(t)+4*Mu*p4(t),
                                  diff(p3(t),
t)=Lambda*p2(t)-3*Mu*p3(t),
                                  p0(0)=1.0,
p1(0)=0.0, p2(0)=0.0, p3(0)=0.0
#, p4(0)=0.0
                               };

                 ? d                                    d
ODESystemInit := ?    p3(t) = 2.17 p2(t) - 10.77 p3(t),    p0(t) = -2.17 p0(t) + 3.59 p1(t),
                 ? dt                                   dt

       d                                                d
          p1(t) = 2.17 p0(t) - 5.76 p1(t) + 7.18 p2(t),    p2(t) = 2.17 p1(t) - 9.35 p2(t) + 10.77 p3(t), p1(0) = 0.,
       dt                                               dt

                                                        178