Допомога у написанні освітніх робіт...
Допоможемо швидко та з гарантією якості!

Метод імовірнісних твірних функцій (аналітичний метод) (реферат)

РефератДопомога в написанніДізнатися вартістьмоєї роботи

Із ненадійним приладом обслуговування В одноканальну систему обслуговування надходить пуассонівський потік вимог із інтенсивністю час обслуговування яких є випадковою величиною, що має експоненціальний закон розподілу з параметром Канал (прилад) обслуговування з деяких причин виходить із ладу. При цьому поломка каналу як подія відбувається у випадкові моменти часу, утворюючи пуассонівський потік… Читати ще >

Метод імовірнісних твірних функцій (аналітичний метод) (реферат) (реферат, курсова, диплом, контрольна)

РЕФЕРАТ На тему:

Метод імовірнісних твірних функцій (аналітичний метод) У темі досліджувалися системи обслуговування, до яких над­ходив один пуассонівський потік вимог, а тому ймовірнісні моделі цих систем мали просту структуру й основні числові характеристики для них визначалися шляхом нескладних алгебраїчних перетворень із використанням умови нормування.

Отже, особливих труднощів під час розв’язування задач такого класу не виникало.

Задача значно ускладнюється, якщо до системи обслуговування надходять два і більше пуассонівські потоки, для котрих, здебільшого, вводиться певна черговість обслуговування вимог, що можуть здійснювати не один, а кілька приладів. Тоді у стаціонарному режимі ймовірнісну модель буде подано системою лінійних алгебраїчних рівнянь високого порядку відносно ймовірностей станів системи. Розв’язати таку систему рівнянь методами, які застосовувалися в попередніх задачах, а також класичними методами лінійної алгебри досить складно, а іноді й неможливо.

Існує аналітичний метод розв’язування таких систем — метод твірних функцій, які називають імовірнісними твірними функціями.

Метод імовірнісних твірних функцій дає змогу звести систему лінійних алгебраїчних рівнянь високого порядку до системи функціональних рівнянь відносно ймовірнісних твірних функцій значно нижчого порядку. Розв’язуючи цю спрощену систему, знаходимо аналітичний вираз зазначених функцій, за допомогою якого визначаються основні числові характеристики.

Спинимося докладніше на понятті ймовірнісних твірних функцій та їх основних властивостях. Поняття ймовірнісних твірних функцій уже розглядалось у темі 5. Тепер приділимо йому значно більше уваги.

6.1. Імовірнісні твірні функції та їх властивості.

Як уже наголошувалося в темі 5, збіжний степеневий ряд виду.

A ( x ) = p 0 + xp 1 + x 2 p 2 + x 3 p 3 + . . . + x k p k + . . . = k = 0 x k p k (245).

називають імовірнісною твірною функцією, де | x | <= 1 , p k  — імовірність того, що система містить k вимог.

Основні властивості A ( x ) :

1. Оскільки.

A ( x ) = p 0 + xp 1 + x 2 p 2 + x 3 p 3 + . . . + x k p k + . . . , .

то p 0 = A ( 0 ) - p 1 = A ' ( 0 ) - p 2 = A ' ' ( 0 ) 2 ! - p 3 = A ' ' ' ( 0 ) 3 ! - . . . .

p k = A ( k ) ( 0 ) k ! . (246).

2. Оскільки A ( x ) = k = 0 x k p k , то.

A ' ( x ) = kx k - 1 p k і за х = 1 дістаємо, що A ' ( 1 ) = kp k . З того, що k = x k , випливає:

M ( x ) = k = 0 kp k = A ' ( 1 ) . (247).

3. Оскільки A ' ' ( 1 ) = k ( k - 1 ) p k = k = 0 k 2 p k - k = 0 kp k ,.

то.

M ( x 2 ) = k = 0 k 2 p k = A ' ' ( 1 ) + A ' ( 1 ) . (248).

Тоді.

D ( x ) = A ' ' ( 1 ) + A ' ( 1 ) - ( A ' ( 1 ) ) 2 . (249).

Якщо в систему надходять два пуассонівські потоки, то ймовірнісна твірна функція в цьому разі має такий вигляд:

A ( x , y ) = i = 1 j = 0 x i y j p ij , (250).

де | x | <= 1, | y | <= 1 , p ij є ймовірність того, що система містить і вимог першого пуассонівського потоку і j вимог другого пуассонівського потоку.

Для системи обслуговування, до якої надходять три пуассонівські потоки, імовірнісна твірна функція подається таким збіжним степеневим рядом.

A ( x , y , z ) = i = 1 j = 0 k = 0 x i y j z k p ijk , (251).

де | x | <= 1, | y | <= 1, | z | <= 1 , p ijk  — імовірність того, що система містить i вимог першого пуассонівського потоку, j вимог другого і k вимог третього такого потоку.

2. Імовірнісна модель (М/М/1).

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

У стаціонарному режимі роботи система подається ймовірнісною моделлю такого виду:

1 . 0 = 1 , 2 . ( + ) p k = k + 1 + k - 1 , k = 1, 2, 3, . . . . { (252).

Безперечно, для визначення основних числових характеристик системи можна скористатися методами, що застосовувалися в попередніх задачах. Проте з метою наочності скористаємося для цієї задачі методом імовірнісних твірних функцій.

Нехай A ( x ) = k = 1 x k p k .

Тоді A ( 1 ) + p 0 = 1 , оскільки це умова нормування.

Помноживши друге рівняння системи на x k , дістанемо:

1 . 0 = 1 , 2 . ( + ) x k p k = k p k + 1 + k p k - 1 , k = 1, 2, 3, . . . . { (253).

Додамо перше та друге рівняння системи і виконаємо підсумування за k:

0 + ( + ) k = 1 x k p k = x k = 1 x k p k + xp 0 + k = 1 x k p k -> .

-> 0 + ( + ) A ( x ) = x A ( x ) + xp 0 + xA ( x ) -> .

-> [ ( 1 - x ) + ( 1 - 1 x ) ] A ( x ) = ( x - 1 ) p 0 -> .

-> A ( x ) = ( x - 1 ) p 0 ( 1 - x ) + ( 1 - 1 x ) -> .

A ( x ) = 0 1 x - = 0 1 x - = xp 0 1 - , = .

За х = 1 маємо.

A ( 1 ) = 1 - p 0 .

Згідно з умовою нормування маємо:

A ( 1 ) + p 0 = 1 -> 1 - p 0 + p 0 = 1 -> p 0 = 1 - ,.

A ( 1 ) = 1 - ( 1 - ) = .

Отже, дістали:

імовірність того, що система простоюватиме (відсутні вимоги), така:

p 0 = 1 -  — (254).

імовірність того, що система зайнята обслуговуванням вимог, становить.

A ( 1 ) = . (255).

Щоб визначити М, скористаємося (247):

M = A ' ( x ) | x = 1 = ( xp 0 1 - ) x = 1 = ( ( 1 - ) + 2 x ( 1 - ) 2 p 0 ) x = 1 = .

= - 2 + 2 ( 1 - ) 2 p 0 = ( 1 - ) 2 p 0 = ( 1 - ) 2 ( 1 - ) = 1 - .

Отже, математичне сподівання кількості вимог, що перебувають у системі:

M = 1 - .

Для систем обслуговування розглядуваного класу визначаємо такі числові характеристики:

  • ередню кількість вимог потоку, що перебуває в черзі:

L = M - A ( 1 ) -> L = 1 - - - (256).

  • ередній час перебування вимог у системі:

T = M = - (257).

  • ередній час перебування вимоги в черзі:

T = L . (258).

3. Імовірнісна модель системи (М/М/1).

із ненадійним приладом обслуговування В одноканальну систему обслуговування надходить пуассонівський потік вимог із інтенсивністю час обслуговування яких є випадковою величиною, що має експоненціальний закон розподілу з параметром Канал (прилад) обслуговування з деяких причин виходить із ладу. При цьому поломка каналу як подія відбувається у випадкові моменти часу, утворюючи пуассонівський потік із параметром Час, що витрачається на відновлення каналу до роботоздатного стану, розглядається як випадкова величина, що має експоненціальний закон розподілу із параметром /p>

Позначимо p k ( t ) імовірність того, що в момент часу t в системі перебуває k вимог і канал перебуває в роботоздатному стані- R k ( t )  — імовірність того, що в системі перебуває k вимог і канал перебуває у стані ремонту. При цьому, у стані ремонту каналу, до системи продовжують надходити нові вимоги. Імовірнісна модель має такий вигляд:

p 0 ' ( t ) = - 0 ( t ) + 1 ( t ) + 0 R 0 ( t ) - p k ' ( t ) = - ( + ) p k ( t ) + k + 1 ( t ) + k - 1 ( t ) + 0 R k ( t ) , k = 1, 2, 3 . . . - R 0 ' ( t ) = - ( + 0 ) R 0 ( t ) + 0 p 0 ( t ) - R k ' ( t ) = - ( + 0 ) R k ( t ) + 0 p k ( t ) + k - 1 ( t ) , k = 1, 2, 3 . . . . { { { (259).

У стаціонарному режимі система (259) набирає такого вигляду:

1 . 0 = 1 + 0 R 0 , 2 . ( + ) p k = k + 1 + k - 1 + 0 R k , 3 . ( + 0 ) R 0 = 0 p 0 , 4 . ( + 0 ) R k = 0 p k + k - 1 . { { { (260).

Нехай A 1 ( x ) = k = 1 x k p k , A 2 ( x ) = k = 1 x k R k .

Помножимо друге і четверте рівняння на x k і просумуємо рівняння за k. Дістанемо систему функціональних рівнянь відносно ймовірнісних твірних функцій A 1 ( x ) , A 2 ( x ) такого вигляду:

[ 0 + ( 1 - x ) + ( 1 - 1 x ) ] A 1 ( x ) - 0 A 2 ( x ) = [ ( x - 1 ) - 0 ] p 0 , [ 0 + ( 1 - x ) ] A 2 ( x ) - 0 A 1 ( x ) = 0 p 0 . { (261).

Розв’язуючи систему (261), маємо:

A 1 ( x ) = [ 0 - 2 ( x - 1 ) + 0 1 ] xp 0 [ 0 + ( 1 - x ) + 2 ( x - 1 ) x - 0 - 0 ] , (262).

A 2 ( x ) = 0 p 0 [ 0 + ( 1 - x ) + 2 ( x - 1 ) x - 0 ] - 0 . (263).

Із (262) і (263) знаходимо:

A 1 ( 1 ) = ( 1 + 0 ) 1 - - 0 p 0  — A 2 ( 1 ) = 0 p 0 1 - - 0 . .

За умовою нормування маємо:

A 1 ( 1 ) + A 2 ( 1 ) + p 0 = 1 ,.

звідки.

( 1 + 0 ) 1 - - 0 p 0 + 0 1 - - 0 p 0 + 1 = 1 -> .

-> p 0 = 1 - - 0 1 + 0 = K Г ( 1 - - 0 ) . (264).

Тут K = 1 1 + 0 називають коефіцієнтом готовності системи. Скориставшись виразом для p 0 , дістанемо:

A 1 ( 1 ) = ( 1 + 0 ) 1 - - 0 p 0 1 - - 0 1 + 0 = - .

A 2 ( 1 ) = 0 1 - - 0 1 - - 0 1 + 0 = 0 1 + 0 = K 0 . .

Отже, маємо:

  • мовірність того, що система перебуває у роботоздатному стані.

A 1 ( 1 ) = - (265).

  • мовірність того, що система перебуває у стані ремонту.

A 2 ( 1 ) = K 0 . (266).

Для визначення решти числових характеристик застосуємо формули (262) — (265).

Так, для визначення математичного сподівання кількості вимог, що перебувають у системі, необхідно знайти A 1 ' ( 1 ) , A 2 ' ( 1 ) :

A 1 ' ( 1 ) = ( 1 + 0 ) - ( 2 + 0 ) ( 1 - - 0 ) ( 1 + 0 ) , (267).

A 2 ' ( 1 ) = 0 - 0 2 + 0 2 + 0 ( 1 - - 0 ) ( 1 + 0 ) . (268).

Отже, одержали:

M = A 1 ' ( 1 ) + A 2 ' ( 1 ) . (269).

Тут 0 = 0  — =  — = 0 .

4. Імовірнісна модель (М/М/1).

із двома пуассонівськими потоками вимог:

1) із надійним каналом (приладом) обслуговування.

В одноканальну систему надходять два пуассонівські потоки вимог із параметрами 1 , 2 , час обслуговування яких є випадковою величиною, що має експоненціальний закон розподілу з відповідними параметрами 1 , 2 . Параметри 1 , 1 стосуються вимог першого потоку, 2 , 2  — другого потоку.

Вимоги першого потоку користуються так званим відносним пріоритетом щодо вимог другого потоку, який називають простим, а саме: коли в системі перебувають вимоги обох потоків і канал обслуговує вимогу простого потоку, то після його завершення перевагу в наступному обслуговуванні надають вимогам першого, пріоритетного потоку. Перший канал обслуговує їх до повної ліквідації черги цього потоку й лише потім починає обслуговувати вимоги простого потоку, коли вони перебувають у черзі. Якщо в системі немає вимог першого і другого потоків, то канал налаштовується на обслуговування вимог простого потоку.

Позначимо p i , k ( t ) імовірність того, що в момент часу t у системі перебуває i вимог простого потоку та k вимог пріоритетного і канал обслуговує вимогу простого потоку;

Q i , k ( t )  — імовірність того, що в момент часу t у системі перебуватиме і вимог простого потоку, k вимог пріоритетного і канал обслуговує вимогу пріоритетного потоку.

Імовірнісна модель динамічного стану системи набере такого вигляду:

1. p 00 ' ( t ) = - ( 1 + 2 ) p 00 ( t ) + 1 Q 01 ( t ) + 2 p 1,0 ( t ) ,.

2. p 1,0 ' ( t ) = - ( 1 + 2 + 2 ) p 1,0 ( t ) + 1 Q 1,1 ( t ) + 2 p 2,0 ( t ) + 2 p 0,0 ( t ) ,.

3. p 1, k ' ( t ) = - ( 1 + 2 + 2 ) p 1, k ( t ) + 1 p 1, k - 1 ( t ) , k = 1, 2, 3, . . . , (270).

4. p i + 1, k ' ( t ) = - ( 1 + 2 + 2 ) p i + 1, k ( t ) + 1 p i + 1, k - 1 ( t ) + 2 p i , k ( t ) , .

i = 1, 2, . . . k = 1, 2, . . . ,.

5. Q 0, k ' ( t ) = - ( 1 + 2 + 1 ) Q 0, k ( t ) + 1 Q 0, k - 1 ( t ) + 1 Q 0, k + 1 ( t ) + 2 p 1, k ( t ) , (260).

6. Q i , 1 ' ( t ) = - ( 1 + 2 + 1 ) Q i , 1 ( t ) + 1 Q i , 2 ( t ) + 2 Q i - 1,1 ( t ) + 2 p i + 1,1 ( t ) ,.

7. Q i , k + 1 ' ( t ) = - ( 1 + 2 + 1 ) Q i , k + 1 ( t ) + 1 Q i , k + 2 ( t ) + 2 Q i - 1, k + 1 ( t ) + .

+ 1 Q i , k ( t ) + 2 p i + 1, k + 1 ( t ) .

У стаціонарному режимі диференціально-різницева система рівнянь (270) переходить у лінійну однорідну систему алгебраїчних рівнянь, яка набирає такого вигляду:

1. ( 1 + 2 ) p 00 = 1 Q 0,1 + 2 p 1,0 ,.

2. ( 1 + 2 + 2 ) p 1,0 = 1 Q i , 1 + 2 p i + 1,0 + 2 p i - 1,0 ,.

3. ( 1 + 2 + 2 ) p 1, k = 1 p 1, k - 1 , .

4. ( 1 + 2 + 2 ) p i + 1, k = 2 p i , k + 1 p i + 1, k - 1 ,.

5. ( 1 + 2 + 1 ) Q 0, k = 1 Q 0, k - 1 + 1 Q 0, k + 1 + 2 p 1, k , (271).

6. ( 1 + 2 + 1 ) Q i , 1 = 1 Q i , 2 + 2 Q i - 1,1 + 2 p i + 1,1 ,.

7. ( 1 + 2 + 1 ) Q i , k + 1 = 1 Q i , k + 2 + 2 Q i - 1, k + 1 + 1 Q i , k + 2 p i + 1, k + 1 .

Систему (271) розв’язуємо методом імовірнісних твірних функцій.

Нехай.

A ( x , y ) = A 1 ( x , y ) + A 2 ( x , y ) + p 0,0 , (272).

де.

A 1 ( x , y ) = i = 1 k = 0 x i y k p ik , (273).

A 2 ( x , y ) = i = 0 k = 1 x i y k Q i , k . (274).

Тут A 1 ( x , y ) , A 2 ( x , y )  — так звані часткові ймовірні твірні функції.

A 1 ( 1, 1 ) = i = 1 k = 0 p ik  — імовірність того, що в системі перебуває і вимог простого потоку та k вимог пріоритетного потоку і система обслуговує вимогу простого потоку;

A 2 ( 1, 1 ) = i = 0 k = 1 Q i , k  — імовірність того, що в системі перебуває і вимог простого потоку та k вимог пріоритетного і система обслуговує вимоги цього потоку.

Помножимо друге рівняння системи (271) на x i , третє - на xy k , четверте — на x i + 1 y k , п’яте — на y k , шосте та сьоме — відповідно на x i y та x i y k + 1 і підсумуємо за індексами і та k. У результаті дістанемо систему функціональних рівнянь відносно часткових імовірнісних твірних функцій виду:

[ 1 ( 1 - y ) + 2 ( 1 - x ) + 2 ] A 2 ( x , y ) = [ 2 ( x - 1 ) - 1 ] p 0,0 + 1 A ( x ) , [ 1 ( 1 - y ) + 2 ( 1 - x ) + 1 ( 1 - 1 y ) ] A 1 ( x , y ) - 2 x A 2 ( x , y ) = 1 yp 0,0 - 1 A ( x ) , { (275).

де A ( x ) = i = 0 x i Q i , 1 + 21 i = 1 x i p i 0 . Тут 21 = 2 1 , що настає внаслідок формування часткових імовірнісних твірних функцій A 1 ( x , y ) , A 2 ( x , y ) . .

Зі здобутої системи функціональних рівнянь (275) дістаємо:

A 1 ( x , y ) = [ 1 2 ( 1 - x ) - 1 1 x ( 1 - y ) - 1 2 x ( 1 - x ) ] A ( x ) [ 1 ( 1 - y ) + 2 ( 1 - x ) + 1 ( 1 - 1 y ) ] [ 1 ( 1 - y ) + 2 ( 1 - x ) + 2 ] x + .

+ [ 2 2 ( x - 1 ) - 1 2 xy ( 1 - x ) + 1 2 xy ( 1 - y ) + 2 1 ( xy - 1 ) ] p 00 [ 1 ( 1 - y ) + 2 ( 1 - x ) + 1 ( 1 - 1 y ) ] [ 1 ( 1 - y ) + 2 ( 1 - x ) + 2 ] x  — (276).

A 2 ( x , y ) = 1 A ( x ) + [ 2 ( x - 1 ) - 1 ] p 00 1 ( 1 - y ) + 2 ( 1 - x ) + 2 . (277).

За x -> 1, y -> 1 знаменник виразу для A 1 ( x , y ) набирає нульового значення. Але оскільки A 1 ( x , y ) є аналітичною функцією, поданою у вигляді збіжного степеневого ряду за | x | <= 1, | y | <= 1 , то корені знаменника (268) будуть також і коренями його чисельника. Справді, знаменник A 1 ( x , y ) набирає нульового значення за x -> 1, y -> 1 , коли.

1 ( 1 - y ) + 2 ( 1 - x ) + 1 ( 1 - 1 y ) | x -> 1 y -> 1 = 0 . (278).

Із (270) випливає, що за x -> 1 значення y також прямує до одиниці. Отже, у буде функцією від х, тобто y = у (х).

Далі нам знадобиться значення y ' ( 1 ) , а тому продиференціюємо (270) за х і спрямуємо x -> 1 . Тоді.

( - 1 y ' ( x ) - 2 + 1 y ' ( x ) y 2 ) | x -> 1 = - 1 y ' ( 1 ) - 2 + 1 y ' ( 1 ) = 0 -> .

-> y ' ( 1 ) = 2 1 - 1 = 2 1 ( 1 - 1 ) . (279).

Отже, доходимо висновку, що.

[ 1 2 ( 1 - x ) - 1 1 x ( 1 - y ( x ) ) - 1 1 ( 1 - x ) ] A ( x ) + .

+ [ 2 2 ( x - 1 ) + 1 2 xy ( x ) ( 1 - x ) + 1 2 xy ( x ) ( 1 - y ( x ) ) + 2 1 ( xy ( x ) - 1 ) ] p 0,0 = 0 .

Звідси.

A ( x ) = [ 2 2 ( 1 - x ) - 1 2 x y ( x ) ( 1 - x ) - 1 2 x y ( x ) ( 1 - y ( x ) ) + 1 2 ( 1 - x ) - 1 1 x ( 1 - y ( x ) ) - + 2 1 ( 1 - x y ( x ) ) ] p 0,0 - 1 2 ( 1 - x ) -> .

-> A ( x ) = [ 2 2 - 1 2 x y ( x ) - 1 2 x y ( x ) 1 - y ( x ) 1 - x + 2 1 1 - x y ( x ) 1 - x ] p 0,0 1 2 - 1 1 - 1 2 1 - y ( x ) 1 - x . .

Далі.

lim x -> 1 A ( x ) = [ 2 2 - 1 2 - 1 2 lim x -> 1 1 - y ( x ) 1 - x + 2 1 lim x -> 1 1 - x y ( x ) 1 - x ] p 0,0 1 2 - 1 1 - 1 2 lim x -> 1 1 - y ( x ) 1 - x = .

оскільки lim x -> 1 1 - y ( x ) 1 - x = ( 0 0 ) , lim x -> 1 1 - xy ( x ) 1 - x = ( 0 0 ) , то згідно з правилом Лопіталя дістаємо: lim x -> 1 1 - y ( x ) 1 - x = lim x -> 1 - y ' ( x ) - 1 = y ' ( 1 ) - lim x -> 1 1 - xy ( x ) 1 - x = lim x -> 1 - y ( x ) - x y ' ( x ) - 1 = 1 + y ' ( 1 ) rli | | | | = .

= 2 2 - 1 2 - 1 2 y ' ( 1 ) + 2 1 ( 1 + y ' ( 1 ) ) 1 2 - 1 2 - 1 1 y ' ( 1 ) p 00 = .

= 2 2 - 1 2 - 1 2 2 1 ( 1 - 1 ) + 2 1 + 2 1 2 1 ( 1 - 1 ) 1 2 - 1 2 - 1 1 2 1 ( 1 - 1 ) p 00 = .

= 2 2 1 2 - 1 2 1 2 + 2 1 1 2 + 2 1 2 1 2 2 1 1 - 1 1 - 1 2 1 2 - 1 1 2 1 2 ( 1 - 1 ) p 00 = .

= 21 2 - 1 2 + 1 + 21 1 2 - 1 2 2 1 - 1 1 - 2 - 1 2 1 - 1 p 0,0 = .

= 21 2 - 21 1 2 - 1 2 + 1 2 2 + 1 - 1 2 + 21 1 2 - 1 2 2 1 - 1 - 2 + 1 2 - 1 2 p 0,0 = .

= 21 2 - 1 2 + 1 - 1 2 1 - 1 - 2 p 0,0 .

Таким чином, маємо:

A ( 1 ) = 21 2 - 1 2 + 1 - 1 2 1 - 1 - 2 p 0,0 . (280).

Із (276) і (277) визначаємо:

A 1 ( 1, 1 ) = 1 1 - 1 - 2 p 0,0 , (281).

A 2 ( 1, 1 ) = 2 1 - 1 - 2 p 0,0 . (282).

Оскільки A 1 ( 1, 1 ) + A 2 ( 1, 1 ) + p 0,0 = 1 , то.

1 1 - 1 - 2 p 0,0 + 2 1 - 1 - 2 p 0,0 + p 0,0 = 1 -> .

-> p 0,0 = 1 - 1 - 2 . (283).

Тоді A 1 ( 1, 1 ) , A 2 ( 1, 1 ) набирають такого вигляду:

A 1 ( 1, 1 ) = 1  — (284).

імовірність того, що система зайнята обслуговуванням вимог пріоритетного потоку;

Показати весь текст
Заповнити форму поточною роботою