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

UptoLike

150
d
dt
p
1
t( ) = -7.51 p
1
t( ) + 2.01 p
2
t( ) + 3.67 p
3
t()
> 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);
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()
> 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);
d
dt
p
3
t( ) = -8.51 p
3
t( ) + 4.17 p
1
t( ) + 3.49 p
4
t()
> 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);
d
dt
p
4
t( ) = -8.15 p
4
t( ) + 2.84 p
2
t()
> ODESystem:={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) };
ODESystem :=
?
?
?
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()
?
?
?
3. Решение уравнений Колмогорова А.Н. для
вероятностей состояний pi(t) в общем виде
PSi:=dsolve(ODESystem):
                                  d
                                    p (t) = -7.51 p1(t) + 2.01 p2(t) + 3.67 p3(t)
                                  dt 1

>  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);
                            d
                              p (t) = -4.85 p2(t) + 3.34 p1(t) + 4.84 p3(t) + 4.66 p4(t)
                            dt 2

> 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);
                                  d
                                    p (t) = -8.51 p3(t) + 4.17 p1(t) + 3.49 p4(t)
                                  dt 3

> 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);
                                         d
                                           p (t) = -8.15 p4(t) + 2.84 p2(t)
                                         dt 4

> ODESystem:={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) };
             ? d                                    d
ODESystem := ?    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                                                d
         p (t) = -7.51 p1(t) + 2.01 p2(t) + 3.67 p3(t),   p (t) = -4.85 p2(t) + 3.34 p1(t) + 4.84 p3(t) + 4.66 p4(t)
       dt 1                                             dt 2
      ?
      ?
      ?


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

PSi:=dsolve(ODESystem):


                                                       150