- •Вінницький фінансово-економічний університет
- •Імітаційне моделювання конспект лекцій
- •Тема 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. Планування імітаційних експериментів у процесі оптимізації систем
- •Питання для самоконтролю
- •Список рекомендованої літератури
- •Навчальне видання:
2. Перевірка значущості коефіцієнтів регресії
Експериментальні дослідження проводять для того, щоб знайти оцінки коефіцієнтів. Теоретичні значення деяких коефіцієнтів можуть бути рівними нулю. Значустість коефіцієнтів лінійної регресії перевіряють окремо для кожного коефіцієнта за допомогою критерію Стьюдента. Послідовність перевірки складається з таких кроків:
знаходимо дисперсію похибки визначення коефіцієнтів (дисперсія є однаковою для всіх коефіцієнтів):
,
(9.5)
де
-
оцінка дисперсії функції відгуку; N, k
–
число спроб і повторень кожної спроби.
2) обчислюємо відношення абсолютного значення коефіцієнта регресії до середньоквадратичного відхилення похибки його визначення:
(9.6)
3) визначаємо число ступенів вільності:
(9.7)
4) обираємо рівень значущості (q =0,05);
5)
у
спеціальній таблиці для заданих
і q
знаходимо критичне значення tкр.;
6)
якщо обчислене за формулою (9.6) значення
відношення більше від критичного (
),
то коефіцієнт
можна вважати значущим. В протилежному
разі можна прийняти нуль-гіпотезу, тобто
коефіцієнт
можна вважати статистично не значущим.
В математичному описі функції відгуку лишаються змінні, коефіцієнти регресії при яких є статично значущими.
Значущість коефіцієнтів квадратичної регресії перевіряють за тими самими правилами, що й лінійної.
У загальному випадку дисперсії похибок визначення коефіцієнтів квадратичної регресії мають обчислюватися для кожного коефіцієнта окремо:
(9.8)
В даній формулі застосовується вже відомі позначення, котрі приводять рівняння квадратичної регресії до лінійного вигляду.
3. Перевірка адекватності моделі
Коефіцієнти, які визначені за методом найменших квадратів, можуть не відповідати (бути неадекватним) спостережуваним значенням ендогенної величини . Тому перед використанням математичної моделі для аналізу досліджуваної системи потрібно переконатися в її адекватності даним експерименту.
Гіпотеза адекватності моделі перевіряється оцінюванням відхилень передбачених значень функції відгуку від експериментально знайдених, усереднених за числом повторень в експериментальних точках факторного простору. Для оцінювання відхилень використовується критерій Фішера.
Найнадійніші результати перевірки гіпотези про адекватність математичної моделі спостережуваним даним дістають у рототабельних планах (повні та дробові факторні плани для лінійної регресії, рототабельні композиційні плани для квадратичних регресій), що забезпечують однакову точність передбачених значень функції відгуку в точках, що містяться на однаковій відстані від центра. Перевірку адекватності роблять у декілька етапів:
1)
Обчислюють статистичну оцінку дисперсії
адекватності
:
(9.9)
де
g – число членів апроксимуючого полінома
;
– значення функції відгуку, обчислене
за допомогою апроксимуючого полінома
в j-й точці плану ;
– експериментальне значення функції
відгуку в j-й
точці плану.
2) Знаходять значення F критерію Фішера:
(9.10)
3)
Визначають число ступенів вільності
і
:
= N – g; (9.11)
= N(k– 1). (9.12)
4) Обираєть рівень значущості q (як правило, q = 0,05).
5) У спеціальній таблиці за заданими q , f4 і f5 знаходять критичне значення параметра Fкр.
6) Якщо обчислене значення параметра F не перевищує табличного Fкр, тобто F<Fкр, то математичний опис функції відгуку рівнянням регресії вважається адекватним. У противному разі гіпотеза про адекватність відкидається і модель вважається не адекватною процесу, що вивчається.
Зауважимо, що перевірка гіпотези про адекватність можлива при f4>0, тобто коли число дослідних точок факторного простору більше від числа членів апроксимуючого полінома. Це необхідно враховувати як при визначенні структури апроксимуючого полінома, так і при виборі відповідного типу факторних планів.
Якщо гіпотеза про адекватність математичного опису досліджуваного процесу відхиляється, то необхідно або перейти до складнішої форми рівняння регресії, або зменшити інтервали варіювання факторів в експерименті.
Зменшення інтервалів варіювання з метою досягнення адекватності математичної моделі спричинюється до зменшення коефіцієнтів регресії, а через це зростає ризик прийняти помилкову гіпотезу про статистичну незначущість деяких коефіцієнтів.
У загальному випадку інтервал варіювання вибирається з умови забезпечення адекватності математичного опису досліджуваного процесу. Часто при виборі необхідних інтервалів варіювання проводяться попередні експрес-спроби, в яких крок варіювання становить 0,05 — 0,3 діапазону змінювання значень рівнів факторного простору.
