ВУЗ:
Составители:
Рубрика:
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
Страницы
- « первая
- ‹ предыдущая
- …
- 149
- 150
- 151
- 152
- 153
- …
- следующая ›
- последняя »
