Соответствующая этому графу система дифференциальных уравнений имеет вид
^ % L = _X p0(/)
dp]}P~— ~ ~Ь нОР\ (О4- ^Ро(О-
at
Преобразование Лапласа для этих уравнений при отыскании за
кона распределения величины Г0(ро(0) = 1; P i(0)= 0) запишет ся так:
($ |
Ро(5)~~ \**Pi (5) —55 |
(5~Ь ^ ~Ь Iх) Р\ (5) 1Ро(5)~ |
|
Решая эту систему относительно изображения p\(s), получим
Рi(s )= S2 + |
\s |
\s |
5 (2Х + (1) + |
Х2 (S + a) (s + Ь) |
|
откуда |
e - b t_ e-at |
|
|
||
px[t)=\ |
а — b |
|
|
|
|
где |
|
|
а-- |
21 + ц — у |
4Х<л + |
(2.7.4)
^__2Х + [I -Ь \' 4X;J. 4- JJ.2
— 2
корни уравнения х2+ (2Л+ц)х + л2 = 0, взятые с обратным знаком.
Тогда искомая плотность распределения случайной величины Т0 будет
f о(0= ~ ~ ~ ~ —^Pi |
e-bt —e~at |
(*> П). |
|
а — b |
|||
dt |
|
Нетрудно убедиться в том, что условие
j u t ) d t = \
о
выполняется.
Действительно:
81
так как произведение корней равно свободному члену. Найдем математическое ожидание Т0:
- |
>-2 |
^ te~b,dt — ^ te~atdt |
||
M[T0] = t{ |
а — Л |
|
|
|
|
|
м |
|
|
к- |
r j ____i_l=>2 JL+jL— J |
I |
Л |
|
а — Ь |
[ *2 |
a1- j ' a-b- |
>. ' |
>.2 |
При отыскании закона распределения случайной величины Т\
начальными условиями будут ро(0)=0; pi(0) = l. В этом случае преобразование Лапласа для дифференциальных уравнений при мет вид
(s-г'-) AjOO-tJ-Pi (s) = 0;
(s -j- >'+!*) PI (s)— Xp0 (s) = s.
Решая эту систему относительно Pi{s), получим
s(s + l)______
Р\ (s) = - s2 + S (2), + jx) + X2
откуда
a— b
где а и b находятся из формул (2.7.4).
Следовательно:
/ i( 0 = x ? , |
l a — b |
|
+ ± z L e-bi] |
|||
|
|
|
b — a |
J |
||
Математическое ожидание времени |
Tx будет равно |
|||||
М [7^]—tx— ^ t f x (t) dt = X |
|
|
oo |
|
||
a— b |
\ te~atdt 4- |
|||||
|
|
о |
J |
|
||
|
|
|
|
о |
|
|
+ - — - |
^ te~bidt |
[ a- ь op- |
|
b — l |
J _ l_ J_ |
|
b — a |
.) |
|
b — a b2 J |
X 1 12 * |
||
|
о |
|
|
|
|
|
Сравнивая это выражение с выражением для математического ожидания случайной величины Г0, убеждаемся в том, что для любых К>0 и р > 0 значение ?0>Г Ь
82
Это объясняется тем, что система будет блуждать по состоя ниям х0 и Х\ в среднем больше, если она в начальный момент
находится |
в |
состоянии |
х0. На |
рис. 2.7.4в показаны |
графики |
|||||||
плотностей распределения f0(t) |
и f\(t) при Я = р=1. |
|
|
|||||||||
2.7.5. |
|
Через систему ПВО прорывается группа самолетов, |
||||||||||
состоящая из т постановщиков помехи одного бомбардировщика. |
||||||||||||
Огонь сначала ведется по постановщикам помех, а затем по бом |
||||||||||||
бардировщику. |
|
Поток |
пусков |
|
|
|
|
|||||
ракет |
можно |
считать |
пуассо |
|
|
|
|
|||||
новским |
с |
интенсивностью Я; |
|
|
|
|
||||||
вероятность |
|
поражения |
одной |
|
|
|
|
|||||
ракетой одного |
|
самолета (лю |
|
|
|
|
||||||
бого) |
равна |
|
р. |
|
После пораже |
|
|
|
|
|||
ния самолета огонь |
мгновенно |
|
|
|
|
|||||||
переносится на следующий са |
|
|
|
|
||||||||
молет. Время |
нахождения са |
|
|
|
|
|||||||
молетов в зоне обстрела ПВО |
|
|
|
|
||||||||
равно |
т. Определить |
вероят |
|
|
|
|
||||||
ность того, что за время обстре |
|
|
|
|
||||||||
ла т бомбардировщик будет по |
|
|
|
|
||||||||
ражен. |
|
|
|
|
|
|
|
|
|
|
|
|
Р е ш е н и е . |
|
|
|
|
|
|
|
|
||||
Группа |
налетающих само |
|
|
|
|
|||||||
летов может быть в следующих |
|
|
|
|
||||||||
состояниях: Xi — сбито |
i само |
|
|
|
|
|||||||
летов |
(/ = 0, |
|
1, |
2,..., т + 1 ) . |
|
|
|
|
||||
Граф состояний системы по |
|
|
|
|
||||||||
казан |
на |
рис. 2.7.5а. |
|
|
|
|
|
|||||
Величина |
/ч |
определяется |
|
|
|
|
||||||
следующим |
образом: К\ = ЯР. |
|
|
|
|
|||||||
Этот граф с точностью до обозначений совпадает с графом, |
||||||||||||
изображенным на рис. 2.7.36* |
|
|
|
|
||||||||
|
|
|
|
^ 7 |
I |
| ^ 7 |
\ |
1 |
Следовательно: |
|
||
|
|
|
|
'V |
dpm-t-1 (') h (V)"1 |
|
||||||
|
|
|
|
|
|
|
|
** |
|
.-М |
||
|
|
|
|
|
|
|
|
|
|
dt |
т\ |
|
|
|
|
|
|
|
Рис. 2.7.5а |
|
|
||||
|
|
|
|
|
|
|
|
|
|
|||
где pm+\(t) — есть вероятность того, что к моменту t бомбарди
ровщик будет поражен.
Тогда искомую вероятность найдем по формуле
ут-ы |
dPт+1(О |
d t = \ —R(jn, >ux), |
dt |
|
где R(m, а) — табулированная функция (см. приложение).
83
Лг.пример, при Х=1 ------; |
т=16 мин; /7 = 0,5; т = 10 получим |
|||
|
|
мин |
|
|
|
(х) = |
1—/?( 10,8)==1 —(1 —0, 184)= |
0Э184. |
|
2.7.6. |
Граф |
состояний |
системы приведен |
на рис. 2.7.6а. Из |
вестна интенсивность пуассоновского потока событий 'k(t), пере
водящего систему из состояния |
х$. Известно |
также, |
что при |
|||||
|
выходе |
системы |
из со |
|||||
|
стояния х0 она |
с |
веро |
|||||
|
ятностью |
|
оказывает |
|||||
|
ся |
в С О С Т О Я Н И И |
X i |
(/= 1, |
||||
|
2, |
..., |
п). |
|
Определить |
|||
|
а) |
интенсивность |
пуас |
|||||
|
соновского |
потока |
со |
|||||
Рис. 2.7.6а |
бытий, |
|
переводящего |
|||||
систему |
из |
состояния |
||||||
|
х0 |
в состояние |
х\ |
и |
||||
б) вероятность pi(t) того, |
что к моменту |
t |
система |
будет в со |
стоянии x it если в момент t = 0 система была в состоянии х0. |
||||
О т в е т |
|
|
|
|
|
t |
- |
| X (х) dx |
|
а) >-|(0 = Ч 0 Йо |
б)/7,(/)= ( >ч(~)е |
0 |
d-. |
|
|
О |
|
|
|
При л(/) = \ = const получим:
|
a) |
\ l= \ql) б) |
Pl(i) = q,[ 1— е~и ]. |
2.7.7. |
Граф |
состояний |
системы имеет вид, показанный на |
рис. 2.7.7а. Требуется найти закон распределения времени Г0,1 нахождения системы в состояниях х0, хх в стационарном режиме,
если интенсивности всех потоков равны 1----
сек
Ре ш е н и е
Всоответствии с графом состояний система из группы со стояний х0ч Х\ может попасть только в состояние х2. Попав в со
стояние х2, система с вероятностью 7г перейдет в состояние х0,
либо с вероятностью 72 — в состояние х3 (а оттуда в состояние *i). Следовательно, блуждание в группе состояний х0, Х\ с равной
вероятностью 72 может начаться из состояния х0 или из состоя
ния Х\ . Граф состояний подсистемы X, с помощью которого мож но найти искомый закон распределения времени Г0,1 блуждания по состояниям х0, хи представлен на рис. 2.7.76.
84
Этот граф с точностью до обозначений интенсивностей пото ков событий совпадает с графом, изображенным на рис. 2.7.46. Следовательно, и решение с точностью до обозначений будет таким же, как и в примере 2.7.4.
Таким образом, если блуждание начинается из состояния л*0,. то плотность распределения времени Тол будет равна
_ з - УГ t |
_ з + /5~ |
t |
m = e- — -— |
— -— |
( < > о ) . |
V 5 |
|
|
Если блуждание начинается из состояния х\, то плотность рас пределения времени Тол будет равна
3 + V 5 |
3 - / 5 |
К — Ч - 3- # ) |
|
т = - |
{t > 0). |
Уъ
Сучетом того, что блуждание в состояниях х0\ х{ может с веро ятностью 72 начаться из состояния х0 и с вероятностью 7г — из состояния Х\, получим следующее выражение для плотности, распределения времени блуждания в состояниях л:0; Х \ \
|
/о,1 (0 = 4" 1/о (*) + |
f i |
= |
|
||
|
_ |
_ 3 + |
t |
_ |
|
(t> 0). |
= — |
( / 5 - 1 ) 6 |
2 |
- f ( / 5 |
+ |
l)e' |
|
4 У |
5 L |
|
|
|
|
|
8S
График этой функции представлен на рис. 2.7.7в. Математическое ожидание времени пребывания в состояниях
-т0, Xj будет равно
/о,1= \ i/o.i{t)dt=-Y сек.
О
2.7.8. Автомашина при ее эксплуатации может находиться в следующих состояниях:
|
— исправна; |
|
A'I — неисправна, проходит осмотр, ко |
|
торый проводится с целью опре |
|
деления вида ремонта; |
|
а'2 — неисправна, проходит капиталь |
|
ный ремонт; |
|
а'з — неисправна, проходит средний ре |
Рис. 2.7.8а |
монт; |
лг4 — неисправна, проходит текущий ре |
|
|
монт. |
Граф состояний автомашины показан на рис. 2.7.8а. Среднее время межремонтного пробега равно 10. Среднее время осмотра машины равно t\. После осмотра автомашина подвергается капи тальному ремонту с вероятностью q2i среднему ремонту с веро ятностью <7з и текущему ремонту с вероятностью (<72+<7з4-<74= = 1). Среднее время проведения капитального ремонта равно t2,
среднее время проведения среднего ремонта равно /з, среднее
время проведения текущего ремонта равно /4. Определить веро ятность того, что машина будет исправна (для стационарного режима), вероятность того, что время простоя машины будет не более t и среднее время простоя машины.
Р е ш е н и е
Найдем интенсивности потоков, переводящих автомашину из состояния в состояние:
/ OA = |
~ Zr ; |
|
<7? . |
__ |
Яз . |
\ |
Я 4 . |
|
м , 2 = |
__ |
У |
/•1,3 = |
------- I |
Л 1 ,4 = |
— , |
||
|
*0 |
|
и |
|
|
/ 1 |
|
h |
л2,0 = |
1 |
|
1 |
|
|
|
|
|
— ; |
^‘З . о = |
__ |
У |
' • 4 , 0 = |
— |
|
|
|
|
h |
|
h |
|
|
и |
|
|
Считая, что все потоки простейшие с выше указанными интенсивностями, найдем вероятность р0 того, что машина будет
исправна для стационарного режима. Решим соответствующую систему однородных алгебраических уравнений:
о = — ^0,1Ро+ *'2,0Р2“Ь ^3.0Pz ^4,0^1>
86