- •Вінницький фінансово-економічний університет
- •Імітаційне моделювання конспект лекцій
- •Тема 1 Вступ до курсу "Імітаційне моделювання"
- •1. Історія виникнення та розвиток імітаційного моделювання
- •2. Види моделювання та особливості їх використання
- •3. Основні напрями використання машинної імітації
- •4. Сфери застосування імітаційних моделей
- •5. Програмна реалізація імітаційних моделей
- •6. Приклад імітаційної моделі
- •III. Gpss — програма імітаційної моделі еом
- •Питання для самоконтролю
- •Тема 2 Сутність, розвиток і застосування імітаційного моделювання
- •1. Сутність та поняття імітаційного моделювання
- •2. Загальна схема і цілі імітаційного моделювання
- •3. Переваги і недоліки методу машинної імітації
- •4. Поняття імітаційного моделювання у вузькому та широкому сенсі
- •5. Імітація еволюційних процесів
- •Питання для самоконтролю
- •Тема 3 Основні етапи побудови імітаційних моделей
- •1. Види робіт при розробці імітаційної моделі.
- •2. Основні етапи побудови імітаційних моделей. Визначення задачі та її аналіз як етап побудови імітаційних моделей
- •3. Визначення вимог до інформації як етап побудови імітаційних моделей
- •4. Збирання інформації як етап побудови імітаційних моделей
- •5. Висунення гіпотез і прийняття припущень як етап побудови імітаційних моделей. Встановлення основного змісту моделі
- •6. Визначення параметрів, змінних і критеріїв ефективності як етап побудови імітаційних моделей. Описання концептуальної моделі і перевірка її вірогідності
- •7. Створення логічної структурної схеми як заключний етап побудови імітаційної моделі
- •Питання для самоконтролю
- •Тема 4 Імітаційна модель керування запасами
- •1. Поняття та сутність оптимального керування запасами
- •2. Основні параметри моделювання задачі керування запасами
- •3. Сутність та характеристика детермінованих моделей керування запасами
- •4. Керування багатопродуктовими запасами
- •5. Концептуальна імітаційна модель керування запасами
- •6. Блок-схема імітаційної моделі керування запасами
- •7. Аналіз результатів машинної реалізації імітаційної моделі
- •Питання для самоконтролю
- •Тема 5 Поняття про Метод Монте-Карло
- •1. Метод Монте-Карло в імітаційному моделюванні
- •2. Приклад застосування методу Монте-Карло
- •3. Точність оцінки ймовірності за допомогою відносної частоти.
- •4. Рівномірна випадкова послідовність чисел рвп [0,1]
- •5. Властивості рівномірної випадковості послідовності чисел
- •Питання для самоконтролю
- •Тема 6.Генерування рвп [0,1]
- •1. Поняття про генератори (датчики) випадкових чисел.
- •2. Табличний спосіб одержання рвп [0,1]
- •3. Фізичний спосіб одержання рвп [0,1]
- •4. Програмні способи одержання рвп [0,1]
- •5. Загально статистичні методи перевірки якості псевдовипадкових чисел
- •Питання для самоконтролю
- •Тема 7 Генерування випадкових подій і випадкових величин під час машинної імітації
- •1. Імітація випадкових подій
- •2. Стандартний метод імітації дискретно розподілених випадкових величин
- •3. Спеціальні методи імітації деяких дискретних розподілів
- •Питання для самоконтролю
- •Тема 8. Генерування неперервних випадкових величин
- •1. Суть проблеми імітації неперервних розподілів
- •2. Стандартний метод імітації неперервних випадкових величин
- •3. Метод добору (відбракування) неперервних випадкових величин
- •4. Наближене формування розподілів неперервних випадкових величин
- •5. Генерування нормально розподілених випадкових величин
- •Питання для самоконтролю
- •Тема 9 Планування імітаційних експериментів
- •1. Загальна характеристика планування імітаційних експериментів
- •2. Апроксимуючий поліном функції відгуку
- •3. Дворівнева система вимірювання факторів
- •4. Повні факторні плани та їхні властивості
- •5. Дробові факторні плани і умови доцільності їх застосування
- •6. Засоби планування експерименту системи statgraphics Планування експерименту (Experimental Design)
- •Питання для самоконтролю
- •Тема 10 Одержання апроксимуючих поліномів
- •2. Апроксимуючий поліном другого плану
- •1. Одержання коефіцієнтів лінійної регресії
- •2. Апроксимуючий поліном другого плану
- •3. Побудова композиційних планів
- •4. Ортогональний центральний композиційний експеримент
- •5. Рототабельний композиційний експеримент
- •Питання до самоконтролю:
- •Тема 11 Узагальнення та статистична перевірка результатів імітаційних експериментів
- •1. Перевірка результатів імітаційних експериментів за допомогою властивості однорідності дисперсій
- •2. Перевірка значущості коефіцієнтів регресії
- •3. Перевірка адекватності моделі
- •Питання для самоконтролю
- •Тема 12 Планування імітаційних експериментів у процесі дослідження та оптимізації систем
- •1. Планування імітаційних експериментів у процесі дослідження систем
- •2. Планування імітаційних експериментів у процесі оптимізації систем
- •Питання для самоконтролю
- •Список рекомендованої літератури
- •Навчальне видання:
6. Приклад імітаційної моделі
Сутність ідеї машинної імітації з’ясуємо на простому прикладі. Нехай обчислювальний центр має одну ЕОМ, що обслуговує користувачів за допомогою системи розподілених терміналів. Проміжки часу між послідовними надходженнями заявок на обслуговування є рівномірно розподіленою випадковою величиною, яка набуває значення від 1, 2,..., 10 хв з однаковою ймовірністю 0,1. Вважатимемо, що час обслуговування однієї заявки (час виконання замовлення) — це випадкова величина, що з однаковою ймовірністю набуває значень 1, 2, ..., 6 хв. Потрібно визначити середній час перебування заявки в обчислювальній системі (очікування в черзі і обслуговування) і відносний час (у відсотках) простою ЕОМ. Обмежимося випадком обслуговування 10 заявок.
Для моделювання системи потрібний штучний експеримент, що відбиває основні умови розглядуваної ситуації. Насамперед фізично реалізуємо випадкові величини: t — час між появами послідовних заявок на обслуговування і — час обслуговування заявки, скориставшись «жеребкуванням».
Візьмемо 10 фішок із числами від 1 до 10, а також кубик, на гранях якого записано цифри від 1 до 6. Витягаючи навмання фішку з будь-якої посудини, реалізуємо величину t (реалізацією є число, записане на цій фіщці). Після цього фішку знову кладемо в посудину. Підкидаючи кубик, на верхній його грані читаємо число, що є реалізацію величини .
І. Головні передпосилки (концептуальна модель)
1. ЕОМ з допомогою системи розподілених терміналів обслуговує користувачів за правилом черги FIFO «Першим прийшов — першим обслужений».
2. Проміжки часу
між послідовними появами заявок —
випадкова величина
,
яка рівномірно розподілена на відрізку
[1, 10], набуваючи на ньому всіх цілих
значень.
3. Час обслуговування
і-ї
заявки
— випадкова величина, що з однаковою
ймовірністю набуває значень 1, 2, 3, 4, 5,
6.
4. Скориставшись імітаційною моделлю, визначити:
1)
— середній час перебування заявок в
обчислювальній системі (час очікування
і час обслуговування),
(1.1)
тут
— час перебування і-ї
заявки в системі;
n — кількість
заявок.
2) L — відносний час простою ЕОМ , %,
(1.2)
тут
— час простою ЕОМ безпосередньо перед
обслуговуванням і-ї
заявки;
— час закінчення обслуговування
останньої n-ї
заявки.
5. Основні співвідношення:
1)
— поточний
(системний) час — момент надходження
і-ї
заявки
(Т1=
0),
(1.3)
2) Toi — час початку обслуговування і-ї заявки,
(1.4)
3) Tki — час закінчення обслуговування і-ї заявки,
(1.5)
Звідси маємо:
(1.6)
(1.7)
Пояснювальну схему (графік Ганта) до модельованої ситуації зображено на рис. 1.2.
Рис. 1.2. Графік Ганта
Проміжні результати
моделювання, що відбивають перехід
імітованої системи з одного стану до
іншого, записують у вигляді табл. 1.1.
Результати одноразового прогону моделі
такі:
Таблиця 1.1
Результати імітації системи
i |
ti |
i |
Ti |
Toi |
Tki |
i |
i |
1 |
0 |
3 |
0 |
0 |
3 |
3 |
0 |
2 |
10 |
1 |
10 |
10 |
11 |
1 |
7 |
3 |
1 |
4 |
11 |
11 |
15 |
4 |
0 |
4 |
3 |
3 |
14 |
15 |
18 |
4 |
0 |
5 |
6 |
4 |
20 |
20 |
24 |
4 |
2 |
6 |
1 |
2 |
21 |
24 |
26 |
5 |
0 |
7 |
5 |
1 |
26 |
26 |
27 |
1 |
0 |
8 |
9 |
3 |
35 |
35 |
38 |
3 |
8 |
9 |
8 |
3 |
43 |
43 |
46 |
3 |
5 |
10 |
8 |
2 |
51 |
51 |
53 |
2 |
5 |
ІІ. Схема алгоритму (ручний варіант)
1. Заносимо в табл. 1.1 початкові значення розглядуваних величин: i=1, t1= 0, T1= 0, T01= 0, 1= 0.
2. Підкинувши кубик, читаємо на верхній його грані значення 1.
3. Знаходимо
4. Відшукуємо
5. Беремо значення і = і+1.
6. Витягаємо фішку з номером ti .
7. Підкидаємо кубик і знаходимо і .
8. Записуємо в табл. 1.1 значення елементів і-го рядка за формулами (1.3) — (1.7).
9. Перевіряємо умову i = n; якщо вона виконується, переходимо до п. 10, у противному разі — до п. 5.
10. Обчислюємо середній час перебування заявки в системі за формулою (1.1).
11. Знаходимо відносний час простою ЕОМ за (1.2).
Зауваження 1. Під час заповнення табл. 1.1 використовувалися рівномірно розподілені випадкові числа.
Зауваження 2. Наведена модель є найпростішою. Її можна розширити, наприклад, розглянутим далі способом.
1. Оскільки в моделі використовуються випадкові величини, то для визначення середнього арифметичного (оцінки математичного сподівання) часу перебування заявки в системі і відносного часу простою ЕОМ потрібно багато разів «прогнати» модель при ідентичних початкових умовах. Щоб дістати статистично достовірні результати, необхідно багато (не менш як 1000 разів) продублювати експеримент.
2. Можна включити вартісні показники, вводячи до системи можливі витрати через затримку виконання замовлення і простій ЕОМ.
3. В імітаційній моделі застосовується правило черги FIFO — «Першим прийшов — першим обслужений». Можна поставити задачу дослідити й інші пріоритети, зокрема правило LIFO — «Останнім прийшов — першим обслужений», SIO — правило найкоротшої операції.
