Применение процесса чистого размножения. Процесс чистого размножения

Марковский процесс с дискретными состояниями называется процессом гибели и размножения , если все состояния можно вытянуть в цепочку, в которой каждое из промежуточных состояний может переходить только в соседние состояния, а крайние состояния переходят лишь в состояния и соответственно. Граф состояний такой системы приведен на рис.4.

Название схемы взято из биологических задач, где состояние популяции означает наличие в ней особей.

На рис.4 переход вправо соответствует увеличению популяции, влево - ее уменьшению. Таким образом, можно определить как интенсивности размножения, а - как интенсивности гибели. Используется следующее соглашение: буквам и приписывается индекс того состояния, из которого выходит стрелка.

Марковским процессом гибели и размножения с непрерывным временем называется такой случайный процесс, исследуемый параметр которого может принимать только целые неотрицательные значения. Изменения рассматриваемого параметра могут происходить в любой момент времени, т.е. в любой момент времени он может либо увеличиться, либо уменьшиться на единицу.

Процессом чистого размножения называется такой процесс, у которого интенсивности всех потоков гибели равны нулю; аналогично процессом чистой «гибели» называется процесс, у которого равны нулю интенсивности всех потоков размножения.

Предельные (финальные) вероятности состояний для простейшего эргодического процесса гибели и размножения, находящегося в стационарном режиме, определяются по следующим формулам:

В качестве примера решения системы уравнений схемы гибели и размножения рассмотрим эксплуатацию автомобилей в крупной транспортной фирме.

Интенсивность поступления автомобилей на предприятие равна. Каждый поступивший на предприятие автомобиль списывается через случайное время. Срок службы автомобиля распределен по показательному закону с параметром. Процесс эксплуатации автомобилей является случайным процессом. - число автомобилей данной марки, находящихся в эксплуатации в момент времени.

Рассмотрим два случая: 1) нет ограничений на число эксплуатируемых автомобилей, 2) на предприятии может эксплуатироваться не более автомобилей.

Если в начальный момент на предприятии не было ни одного автомобиля, то решать систему уравнений нужно при начальных условиях:

Аналогично, если при эксплуатировалось автомобилей, то начальные условия имеют вид:

Решение системы дифференциальных уравнений Колмогорова при произвольном виде функции не может быть найдено в аналитическом виде. Однако при постоянных интенсивностях потоков гибели и размножения и конечном числе состояний будет существовать стационарный режим. Система в этом случае является простейшей эргодической системой.

Если интенсивности потока поступления и списания автомобилей постоянны, то оказываются справедливы формулы:

1. Максимальное число автомобилей не ограничено:

2. Математическое ожидание (среднее значение) числа эксплуатируемых автомобилей.

Простейшее обобщение пуассоновского процесса получается при предположении, что вероятности скачков могут зависеть от текущего состояния системы. Это приводит нас к следующим требованиям.

Постулаты. (i) Непосредственный переход из состояния возможен только в состояние .(ii) Если в момент времени система находится в состоянии , то (условная) вероятность одного скачка в последующем коротком интервале времени между и равна тогда как (условная) вероятность более чем одного скачка в этом интервале есть .

Отличительная черта этого предположения заключается в том, что время, которое система проводит в любом конкретном состоянии, не играет никакой роли; возможны внезапные изменения состояния, однако, пока система находится в одном состоянии, она не стареет.

Пусть снова будет вероятностью того, что в момент времени система находится в состоянии . Эти функции удовлетворяют системе дифференциальных уравнений, которую можно вывести при помощи рассуждений предыдущего параграфа с тем лишь изменением, что (5) в предыдущем параграфе заменяется на

Таким образом, мы получим основную систему дифференциальных уравнений

В пуассоновском процессе было естественно предполагать, что в момент времени 0 система выходит из начального состояния . Теперь мы можем допустить более общий случай, когда система выходит из произвольного начального состояния . Тогда получаем, что

Эти начальные условия единственным образом определяют решение системы (2). (В частности, ). Явные формулы для выводились независимо многими авторами, однако для нас они не представляют интереса.

Пример. Радиоактивный распад. В результате испускания частиц или -лучей радиоактивный атом, скажем урана, может превратиться в атом другого вида. Каждый вид представляет собой возможное состояние, и, когда процесс протекает, мы получаем последовательность переходов . Согласно принятым физическим теориям, вероятность перехода остается неизменной, пока атом находится в состоянии , и эта гипотеза находит выражение в нашем исходном предположении. Стало быть, этот процесс описывается дифференциальными уравнениями (2) (факт, хорошо известный физикам). Если – конечное состояние, из которого невозможны никакие другие переходы, то и система (2) обрывается при . (При мы автоматически получаем ).

Рассмотрим еще одну типичную схему непрерывных марковских цепей - так называемую схему гибели и размножения, часто встречающуюся в разнообразных практических задачах.

Марковский процесс с дискретными состояниями S 0 , S 1 , ..., S n называется процессомгибели и размножения , если все состояния можно вытянуть в одну цепочку, в которой каждое из средних состояний (S 1 , S 2 , ...,
S n -1 ) может переходить только в соседние состояния, которые, в свою очередь, переходят обратно, а крайние состояния (S 0 и S n ) переходят только в соседние состояния (рис. 3.7).

Название взято из биологических задач, где состояние популяции S k означает наличие в ней k единиц особей.

Переход вправо связан с размножением единиц, а влево - с их гибелью.

Рис. 3.7. Граф состояний для процесса гибели и размножения

l 0 (t), l 1 (t), l 2 (t), …, l n (t) - интенсивности размножения;

m 1 (t), m 2 (t), …, m n (t) - интенсивности гибели.

У l и μ индекс того состояния, из которою стрелка выходит.

С состоянием S k связана неслучайная величина Х k : если система S в момент времени t находится в состоянии S k , то дискретная случайная величина X(t) , связанная с функционированием системы, принимает значение k . Таким образом, получаем случайный процесс Х(t), который в случайные, заранее неизвестные моменты времени скачком изменяет свое состояние.

Марковским процессом гибели и размножения с непрерывным временем называется такой случайный процесс, который может принимать только целые неотрицательные значения. Изменения этого процесса могут происходить в любой момент времени, т. е. в любой момент времени он может либо увеличиться на единицу, либо уменьшиться на единицу, либо остаться неизменным.

В практике встречаются процессы чистого размножения и чистой гибели. Процессом чистого размножения называется такой процесс гибели и размножения, у которого интенсивности всех потоков гибели равны нулю; аналогично процессом чистой «гибели» называется такой процесс гибели и размножения, у которого интенсивности всех потоков размножения равны нулю.

Пример 1. Рассмотрим эксплуатацию моделей автомобилей одной марки в крупной транспортной фирме (на предприятии). Интенсивность поступления автомобилей на предприятие равна l(t) . Каждый поступивший на предприятие автомобиль списывается через случайное время T c . Срок службы автомобиля t распределен по показательному закону с параметром m . Процесс эксплуатации автомобилей является случайным процессом. A(t) - число автомобилей данной марки, находящихся в эксплуатации в момент t . Найдем одномерный закон распределения случайного процесса P i (t) = P{A(t) = i}, если: 1) нет ограничений на число эксплуатируемых машин, 2) на предприятии может эксплуатироваться не более n автомобилей.

Решение.

1. Случайный процесс эксплуатации автомобилей есть процесс гибели и размножения, размеченный граф которого представлен на рис. 3.8.

Рис. 3.8. Граф состояний

Система уравнений Колмогорова, соответствующая этому графу, имеет вид

где i = 1, 2, …

Если в начальный момент времени t = 0 на предприятии не было ни одного автомобиля, то решать эту систему уравнений нужно при начальных условиях P 0 (0) = 1, P i (0) = 0 (i = 1, 2, …). Если при t = 0 на предприятии было k автомобилей (k = 1, 2, ...), то начальные условия будут иметь вид

P k (0) = 1, P i (0) = 0 (i = 1, 2, …, i ¹ k ).

2. Если на предприятии может эксплуатироваться не более nавтомобилей моделей одной марки, то имеет место процесс гибели и размножения с ограниченным числом состояний, размеченный граф которого представлен на рис. 3.9.

Рис. 3.9. Граф состояний

Система уравнений Колмогорова для размеченного графа (рис. 3.9) имеет вид (3.4).

Эту систему надо решать при начальных условиях, рассмотренных выше. Решения систем уравнений (3.4) и (3.5) являются одномерными законами распределения Р i (t). Отыскание решений систем в общем виде при произвольном виде функции l(t) представляет значительные трудности и не имеет практических приложении.

При постоянных интенсивностях потоков гибели и размножения и конечном числе состояний будет существовать стационарный режим. Система S с конечным числом состояний (n + 1), в которой протекает процесс гибели и размножения с постоянными интенсивностями потоков гибели и размножения, является простейшей эргодической системой. Размеченный граф состояний для такой системы представлен на рис. 3.9.

Предельные (финальные) вероятности состояний для простейшего эргодического процесса гибели и размножения, находящегося в стационарном режиме, определяются по следующим формулам:

Правило. Вероятность k -гo состояния в схеме гибели и размножения равна дроби, в числителе которой стоит произведение всех интенсивностей размножения, стоящих левее S k , а в знаменателе - произведение всех интенсивностей гибели, стоящих левее S k , умноженной на вероятность кранного левого состояния системы P 0 .

В предыдущем примере для стационарного режима если интенсивность поступления автомобилей постоянная (l(t) = l = const ), то финальные вероятности состояний при условии, что нет ограничений на число автомобилей на предприятии, равны

При этом математическое ожидание числа эксплуатируемых автомобилей равно его дисперсии:

M = D = l /m. (3.10)

Если существует ограничение по числу автомобилей на предприятии (не более n ), то финальные вероятности можно записать в таком виде:

где ρ = l /m .

где k = 0, 1, 2, ..., n .

Математическое ожидание числа эксплуатируемых автомобилей в стационарном режиме

Пример 2. В состав поточной лини входит четыре станка. Бригада в составе четырех человек обслуживающего персонала проводит профилактический ремонт каждого из них. Суммарный поток моментов окончания ремонтов для всей бригады - пуассоновский с интенсивностью l(t). После окончания ремонта станок проверяется; с вероятностью Р он оказывается работоспособным (время проверки мало, и им можно пренебречь по сравнению со временем профилактики). Если станок оказывается неработоспособным, то вновь проводится его профилактика (время на которую не зависит от того, проводилась ли она ранее) и т. д. В начальный момент все станки нуждаются в профилактическом ремонте. Требуется:

1. Построить граф состояний для системы S (четыре станка).

2. Написать дифференциальные уравнения для вероятностей состояний.

3. Найти математическое ожидание числа станков M t , успению прошедших профилактику к моменту t .

Решение.

Граф состояний показан на рис. 3.10, в котором:

S 0 – все четыре станка нуждаются в профилактическом ремонте;

S 1 – один станок успешно прошел профилактику, а три нуждаются в профилактическом ремонте;

S 2 – два станка успешно прошли профилактику, а два нуждаются в профилактическом ремонте;

S 3 – три станка успешно прошли профилактику, один нуждается в профилактическом ремонте;

S 4 – все четыре станка успешно прошли профилактику.

Рис. 3.10. Граф состояний системы

Каждый профилактический ремонт успешно заканчивается с вероятностью P , что равносильно P -преобразованию потока окончаний ремонтов, после которого он останется пуассоновским, но с интенсивностью Pl(t) . В этом примере мы имеем дело с процессом чистого размножения с ограниченным числом состояний.

Уравнения Колмогорова имеют следующий вид:

Начальные условия P 0 (0) = 1, P 1 (0) = … = P 4 (0) = 0. При постоянной интенсивности l(t) = l и вероятности состоянии определяются по следующим формулам:

Математическое ожидание числа дисков, успешно прошедших профилактику к моменту t, равно

где n = 4.

Пример 3. Рассмотрим производство автомобилей на заводе. Поток производимых автомобилей - нестационарный пуассоновский с интенсивностью l(t). Найдем одномерный закон распределения случайною процесса X(t) - число выпушенных автомобилей к моменту времени t , если в момент t = 0 начат выпуск автомобилей.

Решение

Очевидно, что здесь процесс чистого размножения без ограничения на число состояний, при этом l i (t) = l(t) , так как интенсивность выпуска автомобилей не зависит от того, сколько их уже выпушено. Граф состояний такого процесса показан на рис. 3.11.

Рис. 3.11. Граф состояний

Одномерный закон распределения случайного процесса Х(t) для графа, изображенного на рис. 3.11, определяется следующей системой уравнений Колмогорова:

Так как число выпушенных автомобилей X(t) на любой фиксированный момент t распределено по закону Пуассона с параметром

M = D = a(t).

Рассмотренный в этом примере процесс X(t) называетсянеоднородным процессом Пуассона. Если интенсивность l(t) = l = const , то получим однородный процесс Пуассона . Для такого процесса при P 0 (0) = 1, P i (0) = 0 (i > 0)

Характеристиками процесса Пуассона будут

M = D = l×t.

Задача 1. Имеется прибор, который состоит из четырех узлов; поток отказов – простейший, среднее время безотказной работы каждого узла равно 11 час. Отказавший узел сразу начинает ремонтироваться; среднее время ремонта узла равно 2 час. (поток восстановления простейший). Найти среднюю производительность прибора, если при четырех работающих узлах она равна 100%, при трех 60%, при двух и менее прибор вообще не работает.

По мере развития количество клеток, из которых состоит зародыш, увеличивается. Деления клеток (дробление яйца) на самых ранних стадиях развития происходят равномерно (синхронно). Ho у одних видов раньше, у других позже эта синхронность нарушается и клетки, из которых образуются зачатки разных органов, начинают делиться с разной скоростью. Эти различия в скорости деления можно рассматривать как одно из первых проявлений их дифференцировки.

У зародышей млекопитающих уже после стадии 16–32 бластомеров большая часть клеток начинает делиться быстрее и образует трофобласт – зачаток будущей плаценты. Сам будущий зародыш состоит на этих ранних стадиях всего из нескольких клеток. Однако позже в ходе развития и роста зародыш и затем плод становятся во много раз больше плаценты.

У амфибий на стадии бластулы, состоящей из нескольких тысяч клеток, будущая мезодерма составляет менее одной трети всех клеток. Ho по мере развития мезодермальные производные – все мышцы, почти весь скелет, система кровообращения, почки и др. – занимают не менее 80 % всей массы головастика.

Особенно нагляден неодинаковый темп деления клеток в морфогенезе многих беспозвоночных. У видов с мозаичным развитием уже на стадии 30–60 клеток зачатки всех основных органов определены и представлены очень немногими клетками (иногда всего двумя). Далее деления клеток в каждом зачатке строго программируются. Так, например, ранний зародыш асцидий содержит 52 клетки эктодермы, 10 клеток энтодермы и всего 8 клеток мезодермы. В течение последующего развития число клеток эктодермы возрастает в 16 раз, энтодермы – в 20, а мезодермы – в 50. Благодаря программированности делений число клеток у некоторых взрослых беспозвоночных (например, у нематод) строго постоянно и каждый орган представлен определенным числом клеток. Далеко не всегда местоположение органа и место, где делятся составляющие его клетки, совпадают. Часто митозы происходят только в особой зоне размножения и оттуда клетки мигрируют к месту своей дифференцировки. Примеры такого рода мы уже видели при рассмотрении системы стволовых клеток. То же происходит, например, и при развитии головного мозга.

Программа клеточных делений не всегда очень строга и предопределяет точное их число. Чаще, вероятно, деления происходят до тех пор, пока количество клеток или размер органа не достигнет определенной величины. Речь идет, таким образом, о двух принципиально различных механизмах регуляции клеточных делений.

В одном случае (как в яйцах с мозаичным развитием) он, по‑видимому, заключен в самой делящейся клетке, которая должна «уметь отсчитывать» свои деления. В другом же случае должна существовать некоторая «петля обратной связи», когда масса органа или число клеток, достигая некоторой величины, начинает тормозить дальнейшие деления.

Оказалось, что число делений в нормальных клетках, не трансформированных в злокачественные, вообще не беспредельно и обычно не превышает 50–60 (большинство клеток делится меньше, так как если бы яйцо равномерно разделилось 60 раз, то число клеток в организме (260) оказалось бы в тысячи раз выше, чем в действительности). Однако ни механизм такого предела числа клеточных делений (называемого по имени открывшего его ученого предел Хайфлика), ни его биологический смысл пока непонятен.

Что же является «датчиком» в системе регуляции – размер органа или число клеток? Однозначный ответ на этот вопрос дают опыты с получением животных с измененной плоидностью – гаплоидные, триплоидные или тетраплоидные. Их клетки соответственно в 2 раза меньше или в 1,5 или 2 раза больше нормальных диплоидных. Тем не менее и размер самих животных, и размер их органов, как правило, нормальные, т. е. они содержат больше или меньше клеток, чем в норме. Регулируемой величиной, следовательно, является не количество клеток, а масса органа или всего организма.

Иначе обстоит дело у растений. Клетки тетраплоидных растений, как и у животных, соответственно больше диплоидных. Но и размеры частей тетраплоидных растений – листьев, цветков, семян – часто оказываются больше обычных почти в 2 раза. Похоже, что у растений «датчиком» при определении числа клеточных делений является не размер органа, а само число клеток.

Механизмы, регулирующие клеточные деления – пролиферацию клеток, изучаются очень интенсивно и с разных сторон. Одним из стимулов такой активности ученых является то, что отличия раковых клеток от нормальных во многом и состоят в нарушении регуляции клеточных делений, в выходе клеток из‑под такой регуляции.

Примером одного из механизмов регуляции клеточных делений может служить поведение клеток, посеянных на дно флакона с питательной средой, – клеточной культуры. Их деления в хороших условиях происходят до тех пор, пока они не покроют все дно и клетки не коснутся друг друга. Далее наступает так называемое контактное торможение, или торможение, зависимое от плотности клеток. Его можно нарушить, как это делал Ю. М. Васильев, расчистив от клеток небольшое окошко на поверхности стекла. В это окошко со всех сторон устремляются клетки, вокруг него проходит волна клеточных делений. Можно думать, что и в организме контакты с соседними клетками являются механизмом, сдерживающим клеточные деления.

У опухолевых клеток эта регуляция нарушается – они не подчиняются контактному торможению, а продолжают делиться, громоздясь друг на друга. Аналогично, увы, они ведут себя и в организме.

Ho контактное торможение не является единственным механизмом регуляции: ее барьер может быть преодолен и у вполне нормальных клеток. Так, например, плотно прижатые друг к другу клетки печени у молодого животного тем не менее делятся и печень растет вместе с ростом всего животного. У взрослых животных эти деления практически прекращаются. Однако если две доли печени удалить, то в оставшейся доле очень быстро начнутся массовые деления клеток – регенерация печени. Если удалить одну почку, то в течение немногих дней вторая почка за счет клеточных делений увеличится вдвое. Очевидно, что в организме существуют механизмы, способные стимулировать клеточные деления в органе, активировать его рост и приводить размеры органа тем самым в некоторое количественное соответствие с размерами всего организма.

В этом случае действуют не контактные механизмы, а какие‑то химические факторы, может быть связанные с функцией печени или почек. Можно представить, что недостаточность функции этих органов, при удалении части их или при отставании их роста от роста всего организма, так нарушает весь метаболизм в организме, что это вызывает компенсаторную стимуляцию клеточных делений именно в данных органах. Есть и другие гипотезы, объясняющие, например, подобные явления действием особых ингибиторов клеточных делений – кейлонов, выделяемых самим органом; если орган меньше, то меньше и кейлонов и больше клеточных делений в этом органе. Если такой механизм и существует, то действует он не везде. Например, потеря одной ноги не приводит сама по себе к увеличению размеров другой ноги.

Деления стволовых и дифференцирующихся клеток крови стимулируются, как мы уже говорили, гормонами, такими, как, например, эритропоэтин. Гормоны стимулируют клеточные деления и во многих других случаях. Например, стимуляция роста числа клеток яйцевода у кур активируется женским половым гормоном. Существуют химические факторы – обычно это небольшие белки, которые действуют не как гормоны, т. е. не разносятся с кровью по всему организму, а влияют более ограниченно, на соседние ткани. Это известные сейчас факторы роста – эпидермальный и др. Однако в большинстве случаев конкретные химические факторы регуляции клеточных делений и механизмы их действия нам неизвестны.

Еще меньше мы знаем о регуляции клеточных делений во время основных процессов морфогенеза – в эмбриональном развитии. Мы уже говорили, что здесь способность одних клеток делиться быстрее, чем другие, является проявлением их дифференцировки. В то же время нельзя не заметить, что дифференцировка и клеточные деления в определенном смысле противостоят друг другу и иногда даже исключают друг друга. В некоторых случаях это связано с невозможностью деления при далеко зашедшей, терминальной дифференцировке клеток. Может ли, например, разделиться эритроцит с его очень специализированной структурой, жесткой оболочкой и почти полной утратой большинства клеточных функций, а у млекопитающих еще и с потерей ядра? Нервные клетки хотя и сохраняют очень высокий темп метаболизма, но их длинный аксон и дендриты, связанные с другими клетками, служат очевидными препятствиями к делению. Если бы такое деление у нервной клетки все же произошло, это привело бы к потере связи этой клетки с другими и, следовательно, к потере ее функции.

Поэтому обычной последовательностью событий является сначала период пролиферации клеток, а уже затем дифференцировка, носящая терминальный характер. Более того, ряд ученых предполагают, что как раз во время клеточных делений хромосомы как бы «освобождаются» для следующего этапа дифференцировки, – последнему митозу перед дифференцировкой придается особое значение. Эти представления носят пока во многом умозрительный характер п не имеют на молекулярном уровне хороших экспериментальных оснований.

Ho и не зная конкретных механизмов регуляции клеточных делений, мы вправе рассматривать их программированный характер как такое же проявление программы развития, каким являются и все остальные его процессы.

В заключение мы кратко остановимся и на явлении, как бы обратном размножению клеток, – их гибели, которая в определенных случаях формообразования является необходимым этапом развития. Так, например, при образовании пальцев в зачатках кисти передних и задних конечностей клетки мезенхимы собираются в плотные тяжи, из которых потом формируются хрящи фаланг. Среди клеток, оставшихся между ними, происходит массовая гибель, за счет которой отчасти пальцы отделяются друг от друга. Нечто похожее происходит и при дифференцировке зачатка крыла у птиц. Механизмы гибели клеток в этих случаях – факторы, внешние по отношению к клеткам, и события внутри клеток – остаются малоизвестными. А. С. Уманский предполагает, например, что гибель клетки начинается с деградации ее ДНК.

Размножение клеток, несмотря на всю его важность, нельзя считать основным механизмом морфогенеза: в создании формы оно участвует все же косвенно, хотя такие важные параметры, как общая форма органа и его относительные размеры, могут регулироваться именно на уровне клеточных делений. Еще меньшую роль играет в морфогенезе программированная гибель клеток. Ho тем не менее они являются в нормальном развитии совершенно необходимыми компонентами. В регуляции этих явлений участвуют практически все компоненты клетки и ее генетический аппарат. Это показывает нам, что в развитии не бывает простых процессов. Попытка до конца разобраться в любом из них заставляет нас обращаться к основным молекулярным механизмам работы клетки. А здесь еще много нерешенного.

Для того чтобы оценить всю сложность развития многоклеточного организма, надо представить себе этот процесс происходящим как бы в многомерном пространстве. Одну ось составляет длинная цепь этапов реализации генетической информации – от гена до признака. Второй такой осью можно назвать всю совокупность генов в хромосомах. В ходе развития продукты различных генов взаимодействуют друг с другом. Развертывание событий но двум осям образует как бы сеть на плоскости. Однако существует п третья ось – разнообразие событий, происходящих в разных частях зародыша. События эти могут происходить относительно автономно, как у животных с мозаичным развитием. Ho частично и у них, а в полной мере у видов с регуляционным типом развития между частями организма осуществляются большие или меньшие взаимодействия и всегда сложные перемещения клеток. Рассматривать их все как одну ось можно, только идя на значительные упрощения. И наконец, все развитие (гаметогенез, эмбриогенез и постэмбриональное развитие) происходит во времени, масштаб которого совершенно иной, чем время, измеряемое на пути от гена до белка. По этой (условно четвертой) оси вся многомерная картина радикально изменяется – яйцо превращается в размножающийся организм. Эта многомерность иллюстрирует сложность всех процессов и их взаимоотношений и трудности их понимания.

У части вирусов роль наследственного вещества выполняет не ДНК, а сходная с ней по строению РНК.

В предыдущем параграфе мы убедились, что зная размеченный граф состояний системы, можно сразу написать алгебраические уравне­ния для предельных вероятностей состояний. Таким образом, если две непрерывные цепи Маркова имеют одинаковые графы состояний и раз­личаются только значениями интенсивностей λ ij , то нет надобности находить предельные вероятности состояний для каждого из графов в от­дельности: достаточно составить и решить в буквенном виде уравнения для одного из них, а затем подставить вместо λ ij , соответствующие зна­чения. Для многих часто встречающихся форм графов линейные урав­нения легко решаются в буквенном виде.

Рис. 29

В данном параграфе мы познакомимся с одной очень типичной схемой непрерывных марковских цепей – так называемой «схемой гибели и размножения».

Марковская непрерывная цепь называется «процессом гибели и размножения», если ее граф состояний имеет вид, представленный на рис.29, т. е. все состояния можно вытянуть в одну цепочку, в кото­рой каждое из средних состояний (S 2 , ..., S n –1) связано прямой и об­ратной связью с каждым из соседних состояний, а крайние состояния (S 1 , S n) – только с одним соседним состоянием.

Пример 1. Техническое устройство состоит из трех одинаковых узлов; каждый из них может выходить из строя (отказывать); отказавший узел немедленно начинает восстанавливаться. Состояния систе­мы нумеруем по числу неисправных узлов:

S 0 – все три узла исправны;

S 1 – один узел отказал (восстанавливается), два исправны;

S 2 – два узла восстанавливаются, один исправен;

S 3 – все три узла восстанавливаются.

Граф состояний показан на рис.30. Из графа видно, что про­цесс, протекающий в системе, представляет собой процесс «гибели и размножения».

Рис. 30

Схема гибели и размножения очень часто встречается в самых раз­нообразных практических задачах; поэтому имеет смысл заранее рас­смотреть эту схему в общем виде и решить соответствующую систему алгебраических уравнений с тем, чтобы в дальнейшем, встречаясь с кон­кретными процессами, протекающими по такой схеме, не решать зада­чу каждый раз заново, а пользоваться уже готовым решением.

Итак, рассмотрим случайный процесс гибели и размножения с гра­фом состояний, представленным на рис.31

Рис. 31

Напишем алгебраические уравнения для вероятностей состояний. Для первого состояния S 1 имеем:

λ 12 p 1 = λ 21 p 2 (7.1)

Для второго состояния S 2 суммы членов, соответствующих входя­щим и выходящим стрелкам, равны:

λ 23 p 2 + λ 21 p 2 = λ 32 p 3 + λ 12 p 1

Но, в силу (7.1), можно сократить справа и слева равные друг дру­гу члены и; получим:

λ 23 p 2 = λ 32 p 3

λ 34 p 3 = λ 43 p 4

. . . . . . . . .

Одним словом, для схемы гибели и размножения члены, соответ­ствующие стоящим друг над другом стрелкам, равны между собой:

λ k -1,k p k -1 = λ k,k -1 p k (7.2)

где k принимает все значения от 2 до n.

Итак, предельные вероятности состояний р 1 , р 2 , ..., р n в любой схеме гибели и размножения удовлетворяют уравнениям:

λ 12 p 1 = λ 21 p 2

λ 23 p 2 = λ 32 p 3

λ 34 p 3 = λ 43 p 4

. . . . . . . . . . . (7.3)

λ k -1,k p k -1 = λ k,k -1 p k

. . . . . . . . . . .

λ n -1, n p n -1 = λ n , n -1 p n

и нормировочному условию:

р 1 + р 2 + ... + р n = l (7.4)

Будем решать эту систему следующим образом: из первого урав­нения (7.3) выразим р 2:

Из второго, с учетом (7.5), получим:

(7.6)

Из третьего, с учетом (7.6):

(7.7)

Эта формула справедлива для любого k от 2 до n .

Обратим внимание на ее структуру. В числителе стоит произведе­ние всех плотностей вероятности перехода (интенсивностей) λ ij , стоя­щих у стрелок, направленных слева направо, с начала и вплоть до той, которая идет в состояние Sk ; в знаменателе – произве­дение всех интенсивностей λ ij , стоящих у стрелок, идущих справа налево, опять-таки, с начала и вплоть до стрелки, исходящей из состояния S k . При k =n в числителе будет стоять произведение интен­сивностей λ ij , стоящих у всех стрелок, идущих слева направо, а в знаменателе – у всех стрелок, идущих справа налево

Итак, все вероятности р 1 , р 2 , ..., р n выражены через одну из их: p 1. Подставим эти выражения в нормировочное условие: р 1 + р 2 + ... + р n = l Получим:

Остальные вероятности выражаются через p 1:

(7.9)

Таким образом, задача «гибели и размножения» решена в общем виде: найдены предельные вероятности состояний.

Пример 2. Найти предельные вероятности состояний для процесса гибели и размножения, граф которого показан на рис. 32.

Рис. 32

Решение. По формулам (7.8) и (7.9) имеем:

Пример 3. Прибор состоит из трех узлов; поток отказов – простейший, среднее время безотказной работы каждого узла равно . Отказавший узел сра­зу же начинает ремонтироваться; среднее время ремонта (восстановления) узла равно ; закон распределения этого времени показательный (поток восста­новлений простейший). Найти среднюю производительность прибора, если при трех работающих узлах она равна 100%, при двух – 50%, а при одном и менее – прибор вообще не работает.

Решение. Перечень состояний системы и граф состояний уже приводились в примере 1 данного параграфа. Разметим этот граф, т. е. проставим у каждой стрелки соответствующую интенсивность λ ij (см. рис. 33.).

Рис. 33.

Так как поток отказов каждого узла – простейший, то промежуток вре­мени между отказами в этом потоке распределен по показательному закону с па­раметром , где –среднее время безотказной работы узла.

По стрелкам вправо систему переводят отказы. Если система на­ходится в состоянии S 0 , то работают три узла; каждый из них подвергается по­току отказов с интенсивностью ; значит, поток отказов, действующий на всю систему, в три раза более интенсивен: = 0,5 и

Средняя производительность прибора в установившемся режиме:

Номинала.

Происхождение термина «схема гибели и размножения» ведет начало от биоло­гических задач, где подобной схемой описывается процесс изменения числен­ности популяции