1 ... 18 19 20 21 22 23 24 25 ... 36 Рисунок 3.8Знайдемо частинний розв’язок системи диференціальних рівнянь Колмогорова (3.17), яка для цього випадку має вигляд p1 t p1 t, p t k1p t kp t, k k1 k (3.25) k 1, 3, , n1, pn t n1 pn1 t. Згідно з (3.18) 1 pt et, t 0. Використовуючи рекурентну формулу (3.22), отримаємо t t t k 2 p t e2t p e2d e2t e e2d e2t ed 2 1 0 0 0 e2t et1 et1 et ; t t k 3 p3 t e3t 2p2 e3d e3t 2 e1 ee3d 0 0 t t 2e3t e2 1 ed 2e3t e2 ed 2e3t 0 0 1 e2t 1 1 et 1 et1 2et e2t et1 et2 . 2 2 Аналогічно pt et1 etk 1 , t 0, 1 k n1. k Таким чином, процесчистогорозмноженнязінтенсивністю kв момент часу tє випадкове число народжень Xt в інтер- валі 0; t, яке розподілене за геометричним законом з ймовірністю ус- піху p et. Отже, n1 pnt 1 pkt, k 1 n ; M Xt 1 p et; p2 D Xt 1p et1 et. Процесичистоговимирання Розглянемо процес експлуатації деяких виробів (банкоматів, авто- машин, комп’ютерів тощо) на підприємстві. Нехай Xt – число виробів, які знаходяться в експлуатації на момент часу t, при цьому X0 nта Xt набуває значень k n, n 1, , 0. Припустимо, що нові вироби на підприємство не надходять, а ті, що виходять з ладу, – списуються. Вважаємо, що потік списаних виробів – найпростіший з параметром k. Потрібно знайти одновимірний закон розподілу випадкового про- цесу Xt, а також математичне сподівання mxt та дисперсію Dxt випадкового процесу Xt. Зазначений процес моделюється процесом чистого вимирання, розмічений граф якого зображено на рис. 3.6. Система диференціаль- них рівнянь Колмогорова (3.16) матиме вигляд p t pt, n n n p t p t p t, k k 1 k 1 k k (3.26) k 1, 3, , n1, p0 t 1 p1 t, з початковою умовою pn 0 1. Знайдемо розв’язок n-го рівняння системи, що задовольняє почат- кову умову pn0 1: dpndt npn, dpnpn n dt, dpnpn n dt, ln pn nt ln C, ln pn C nt, n p t C ent, тоді, враховуючи початкову умову, маємо n p0 1 1 C e0 C 1. гляд Отже, частинний розв’язок n-го рівняння системи (3.26) має ви- n pt ent, t 0. (3.27) Розглянемо k-те рівняння системи (3.26), k 1, 3, , n 1 та знай- демо його розв’язок методом варіації довільної сталої. Розв’яжемо спо- чатку відповідне однорідне рівняння pk t kpkt 0, pk t kpkt, dpk dt, k pk ln pk ln pk Ck kt ln Ck, kt, k k pt C ekt. Нехай Ck Ckt, тоді k k pt Ct ekt, (3.28) звідки k k k k p t C t ekt C t ekt. Після підстановки значень pkt і pk t в k-те рівняння маємо k k k 1 C t ekt Ck ekt pk1 t kCk t ekt, k k 1 C t ekt pk1 t, Таким чином, Ck t k1 pk1 tekt. Ck t k1 pk1 t ekt. (3.29) Проінтегрувавши обидві частини рівності (3.29) від 0 до t, мемо мати- Ckt Ck 0 t k 1 0 t pk 1 ek d , (3.30) Ckt Ck 0 k 1 0 pk1 ek d . Ураховуючи значення Ckt (3.30), розв’язок (3.28) матиме вигляд t pkt Ck 0 k 1 0 pk 1 e kd e kt. Згідно з початковою умовою pk0 Ck0 0. Отже, для 1 k n 1 отримаємо рекурентну формулу t k pt ekt 0 k 1 pk 1 ek d . (3.31) яка разом з (3.27) дозволяє послідовно знайти ймовірності pnt, pn1 t, ..., p1 t, а також n1 p0 t 1 pkt. k 0 (3.32) Математичне сподівання та дисперсію випадкового процесу знаходять, використовуючи формули Xt n mxt M Xt kpkt; k1 t D Xt kp n D 2 x k k 1 t M2 Xt. Розглянемо окремі випадки процесу чистого вимирання. І. Процес чистого вимирання з інтенсивністю k , 1 k n. |