- •Міністерство освіти і науки україни національний транспортний університет Дослідження операцій в моделюванні управлінських рішень
- •Isbn 978-966-632-185-8
- •I. Предмет і завдання дослідження операцій
- •1.1. Що таке дослідження операцій
- •1.2. Основні поняття та принципи дослідження операцій
- •1. 3. Математичні моделі операцій
- •1.4. Багатокритеріальні задачі дослідження операцій
- •1.5. Системний підхід до задач дослідження операцій
- •1.6. Типові класи задач дослідження операцій
- •1. Задачі управління запасами
- •2. Задачі розподілу ресурсів
- •3. Задачі ремонту і заміни обладнання
- •4. Задачі масового обслуговування
- •5. Задачі впорядкування
- •6. Задачі сітьового планування і управління (спу)
- •7. Задачі вибору маршруту (сітьові задачі на транспорті)
- •8. Комбіновані задачі
- •II. Класифікація методів оптимізації
- •2.1. Математичне формулювання загальної задачі оптимізації
- •2.2. Геометрична інтерпретація задач оптимізації
- •2.3. Класифікація методів оптимізації по виду цільової функції
- •2.4. Характеристика градієнтних методів пошуку екстремуму
- •III. Лінійне програмування
- •3.1. Задачі лінійного програмування
- •3.2. Основна задача лінійного програмування
- •3.3. Транспортна задача лінійного програмування
- •IV. Динамічне програмування
- •4.1. Метод динамічного програмування
- •4.2. Приклади розв’язання задач динамічного програмування
- •4.3. Завдання динамічного програмування в загальному вигляді. Принцип оптимальності
- •V. Елементи теорії масового обслуговування
- •5.1. Основні поняття теорії масового обслуговування
- •5.2. Класифікація систем масового обслуговування
- •5.3. Основні елементи систем масового обслуговування
- •5.4. Параметри смо і загальна методика дослідження
- •5.5. Системи масового обслуговування з відмовами
- •5.6. Кількісні показники смо з відмовами
- •5.7. Короткий опис основних типів пуассонівських смо
- •5.7.1. Смо з відмовами:
- •5.7.2. Смо з очікуванням і обмеженим потоком заявок
- •5.7.3. Смо змішаного типу з обмеженням по довжині черги
- •5.7.4. Смо змішаного типу з обмеженням за часом перебування заявки в черзі і на обслуговуванні
- •5.7.5. Приклади розв’язку задач тмо
- •VI. Теорія ігор
- •6.1. Ігрові моделі прийняття рішень
- •6.2. Прямокутні матричні ігри
- •6.3. Аналіз матричних ігр
- •6.4. Елементарні методи розв’язку ігор
- •6.4.1. Загальна схема розв’язання
- •6.4.2. Методи розв’язку гри 2 X 2
- •6.4.3. Методи розв’язання ігор 2хn
- •Ордината точки n дорівнює ціні гри ν, а абсциса дорівнює частоті застосування стратегії а1.
- •Аналітичний метод розв’язання гри 2хn
- •6.5. Загальний розв’язок гри m X n методом лінійного програмування
- •6.6. Наближені методи розв’язання матричних ігр
- •6.7. Приклади розв’язання задач теорії ігор
- •VII. Моделювання на пк
- •7.1. Поняття моделі і моделювання
- •7.2. Метод статистичних випробувань
- •7.3. Імітація випадкових впливів на пк
- •7. 4. Методи формування в пк базових впливів
- •7.5. Оцінка точності характеристик, отриманих методом статистичних випробувань. Необхідна кількість реалізацій
- •7.6. Приклади розв’язку задач методом статистичних випробувань
- •VIII. Метод cітьового планування
- •8.1. Поняття про сітьове планування та управління
- •8.2. Основні визначення
- •8.3. Основні елементи сітьового графіка
- •8.4. Правила побудови сітьового графіка
- •8.5. Часові параметри сітьових графіків
- •8.6. Способи розрахунку сітьових графіків
- •8.7. Оптимізація сітьових графіків
- •8.8. Імовірнісні тимчасові оцінки сітьових моделей
- •8.9. Укрупнення сітьових моделей
- •8.10. Зшивання сітьових моделей
- •Література
- •Для нотаток
8.8. Імовірнісні тимчасові оцінки сітьових моделей
Тимчасові оцінки тривалості виконання окремих робіт комплексу можуть бути детермінованими і імовірнісними. Детермінована оцінка задається одним числом — tij.
На практиці оцінки тривалості безлічі робіт не можуть бути визначені однозначно і є невизначеними. У цьому випадку використовуються імовірнісні оцінки тривалості робіт, розрахунок цих оцінок можна здійснювати двома способами.
Перший спосіб полягає в тому, що час виконання роботи розглядається як неперервна випадкова величина, що заповнює інтервал, обмежений межами, апріорно задається відповідальними виконавцями у вигляді мінімальної (оптимістичній) оцінки тривалості tmin=аi максимальної (неоптимістичні) tmax = b.
Одночасно вони дають і найбільш ймовірну оцінку тривалості робіт tнв=m
Мінімальна, або оптимістична, оцінка тривалості визначається виходячи з припущення, що роботи будуть виконуватися за найсприятливішого збігу обставин. Це означає, що швидше завершити роботу практично неможливо.
Максимальна, або неоптимістична, оцінка визначається для самого несприятливого збігу обставин виконання роботи. Ця оцінка передбачає більшу, ніж зазвичай, кількість невдач, зривів, перешкод і т.д. Однак стихійні лиха та нещасні випадки при цьому до числа несприятливих умов не відносяться. Найбільш ймовірна оцінка тривалості визначається для умов виконання робіт, які найчастіше зустрічаються.
Як правило, характеристики а, b і m визначаються на основі експертних оцінок, згідно з якими розраховується математичне сподівання тривалості роботи
,
а також дисперсія для тривалості операції
.
Як показує досвід, визначення оцінки т викликає труднощі у відповідальних виконавців, у той же час визначення оцінок а і b таких ускладнень не викликає, тому в багатьох випадках кращі формули, запропоновані видомим вченим Д. І. Голенко:
и
З цих формул видно, що чим більший розмах (b-а), тим більша невизначеність, і навпаки.
Другий спосіб охоплює більш загальний випадок, коли характеристики а і b невідомі. Відхилення випадкової величини ti (часу виконання роботи аi) від її заздалегідь заданого значення ti) може бути в обидві сторони, як у більшу (запізнення), так і в меншу (випередження).
При цьому виникають наступні запитання:
— яка ймовірність того, що фактичний час виконання комплексу робіт Т не перевершить заданої величини Т0?
— як слід організувати комплекс робіт для того, щоб величина Т не перевершила заданого Т0 з досить високою ймовірністю?
Розглянемо перше питання (тим більше, що для відповіді на друге, перш за все, треба відповісти на перше). Припустимо, що часи виконання роботі t1, t2,..., tn являють собою випадкові величини з відомими законами розподілу і для простоти будемо вважати, що ці випадкові величини незалежні і щільності їх рівні
Розглянута функція цих випадкових величин — загальний час виконання всього комплексу робіт
Поставлена задача буде вирішена, якщо вдасться знайти функцію розподілу випадкової величини Т:
(8.4)
Тоді підставляючи в неї замість t величину Т0, ми знайдемо шукану ймовірність.
Функція F(t) у загальному випадку є досить складною, так як сам критичний шлях випадковий і залежить від тих значень, які приймають випадкові величини ti. При одних значеннях ti може бути один критичний шлях, при інших-другий. Однак, якщо обмежитися тільки порівняно малими відхиленнями випадкових величин ti від своїх номінальних значень ti(0)(так, що критичний шлях залишається тим же), то завдання сильно спрощується. Тоді у формулі (8. 4) фігурують тільки декілька цілком певних випадкових величин ti. Закон розподілу випадкової величини Т являє собою в цьому випадку не що інше, як композицію законів розподілу випадкових величин ti, що відносяться до критичних робіт.
Надалі нам приходить на допомогу сама складність плану і наявність на критичному шляху багатьох робіт. Ми знаємо, що при додаванні досить великого числа незалежних випадкових величин, розподілених за будь-яких законів і порівняльних по порядку дисперсій, закон розподілу суми виявляється близьким до нормального (центральна гранична теорема). Тому, якщо на критичному шляху стоїть досить велика кількість робіт (порядку 5-6 і більше), то на практиці можна вважати наближено величину Т нормально розподіленою з математичним очікуванням
,
де mti — математичне сподівання часу виконання i-ї роботи, і середнім квадратичним відхиленням
,
де
— середнє квадратичне відхилення часу
виконання i-ї
роботи. Щільність ймовірності випадкової
величини Т буде мати вигляд
.
У даному випадку для знаходження закону розподілу часу виконання комплексу робіт немає потреби знати закони розподілу fi(t) окремих часів ti, досить знати їх математичні очікування і середні квадратичні відхилення. Якщо ці величини відомі, ймовірність виконання комплексу в термін Т знайдеться за відомою формулою
,
де
—
функція Лапласа.
Якщо при випадкових змінах часів ti, може змінюватися і сам критичний шлях, задача обчислення ймовірності Р (Т <Т0) неможлива. При порівняно малому числі робіт у комплексі це завдання може бути вирішена аналітичним способом, але при великому їх числі розрахунки стають громіздкими, тому на практиці виявляється зручніше визначати ці ймовірності методом Монте-Карло на ЕОМ. При цьому розігруються значення випадкових часів ti. Для кожної сукупності отриманих значень визначається час Т виконання комплексу робіт тим способом, який застосовується для невипадкових часів. Отримавши достатньо велику кількість N таких реалізацій, ми можемо безпосередньо знайти математичне сподівання, дисперсію і середнє квадратичне відхилення випадкової величини T. Як наближеного значення ймовірності Р (Т <To) можна взяти частоту цієї події в серії реалізації, оцінивши при цьому довірчі ймовірності та довірчі інтервали відомими методами математичної статистики.
Слід звернути увагу ще на один момент, пов’язаний із застосуванням ЕОМ при сітьовому плануванні. Зазвичай при виконанні складних комплексів робіт спочатку намічені плани не виконуються; їх доводиться по ходу роботи переглядати. При цьому надзвичайно зручно тримати всі дані про комплекс (як початковий план, так і поступає інформацію про його порушення) в пам’яті ЕОМ, яка час від часу заново переглядає план робіт, знаходить для кожного моменту часу новий критичний шлях («загрозливі» за термінами роботи) і оптимізує план, зазначаючи, які саме роботи і якою мірою слід форсувати.
