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

UptoLike

151
> diff(rhs(PSi[1]), t):
> rhs(PSi[1]):
> diff(rhs(PSi[1]), t):
4. Решение уравнений Колмогорова А.Н. для
вероятностей состояний pi(t) при заданных
начальных условиях
> ODESystemInit:={diff(p[1](t), t)=-add(Lij[1,j],
j=1..n)*p[1](t)+add(Lij[i,1]*p[i](t),i=1..n),
diff(p[2](t),
t)=-add(Lij[2,j],
j=1..n)*p[2](t)+add(Lij[i,2]*p[i](t),i=1..n),
diff(p[3](t),
t)=-add(Lij[3,j],
j=1..n)*p[3](t)+add(Lij[i,3]*p[i](t),i=1..n),
diff(p[4](t),
t)=-add(Lij[4,j],
j=1..n)*p[4](t)+add(Lij[i,4]*p[i](t),i=1..n),
p[1](0)=0.2,
p[2](0)=0.25, p[3](0)=0.3, p[4](0)=0.25
};
PSiInit:=dsolve(ODESystemInit):
ODESystemInit :=
?
?
?
d
dt
p
4
t( ) = -8.15 p
4
t( ) + 2.84 p
2
t(),
d
dt
p
3
t( ) = -8.51 p
3
t( ) + 4.17 p
1
t( ) + 3.49 p
4
t(),
d
dt
p
1
t( ) = -7.51 p
1
t( ) + 2.01 p
2
t( ) + 3.67 p
3
t(),
d
dt
p
2
t( ) = -4.85 p
2
t( ) + 3.34 p
1
t( ) + 4.84 p
3
t( ) + 4.66 p
4
t(), p
1
0( ) = 0.2, p
2
0( ) = 0.25, p
3
0( ) = 0.3,
p
4
0( ) = 0.25
?
?
?
5. Визуализация решений уравнений Колмогорова
А.Н. для вероятностей состояний pi(t)
при заданных начальных
условиях
> p1(t):=rhs(PSiInit[2]):
>   diff(rhs(PSi[1]), t):
>   rhs(PSi[1]):
>   diff(rhs(PSi[1]), t):
4. Решение уравнений Колмогорова А.Н. для
вероятностей состояний pi(t) при заданных
начальных условиях

> ODESystemInit:={diff(p[1](t), t)=-add(Lij[1,j],
j=1..n)*p[1](t)+add(Lij[i,1]*p[i](t),i=1..n),
                                 diff(p[2](t),
t)=-add(Lij[2,j],
j=1..n)*p[2](t)+add(Lij[i,2]*p[i](t),i=1..n),
                                 diff(p[3](t),
t)=-add(Lij[3,j],
j=1..n)*p[3](t)+add(Lij[i,3]*p[i](t),i=1..n),
                                 diff(p[4](t),
t)=-add(Lij[4,j],
j=1..n)*p[4](t)+add(Lij[i,4]*p[i](t),i=1..n),
                                 p[1](0)=0.2,
p[2](0)=0.25, p[3](0)=0.3, p[4](0)=0.25
                               };
PSiInit:=dsolve(ODESystemInit):
                     ? d                                    d
    ODESystemInit := ?    p4(t) = -8.15 p4(t) + 2.84 p2(t),    p3(t) = -8.51 p3(t) + 4.17 p1(t) + 3.49 p4(t),
                     ? dt                                   dt

           d
             p (t) = -7.51 p1(t) + 2.01 p2(t) + 3.67 p3(t),
           dt 1
           d
             p (t) = -4.85 p2(t) + 3.34 p1(t) + 4.84 p3(t) + 4.66 p4(t), p1(0) = 0.2, p2(0) = 0.25, p3(0) = 0.3,
           dt 2
                      ?
          p4(0) = 0.25?
                      ?


5. Визуализация решений уравнений Колмогорова
А.Н. для вероятностей состояний pi(t)
                                           при заданных начальных
условиях

>   p1(t):=rhs(PSiInit[2]):
                                                        151