Отриманий УВК має вичерпну інформацію про систему змішування дії всіх змінних, що входять в повну регресійну модель (табл. 8.9).
Таблиця 8.9
Система змішування дії змінних у ДФЕ
Список
ВК
УВК
Система
суттєвих
змішування
змінних
Дослідження систем методом планування експериментів
389
Аналогічно встановлюються системи змішування та розв’язувальні здатності реплік довільної дробовості. В табл. 8.10 наведена матриця планування одного із альтернативних планів ДФЕ для семи факторів.
Таблиця 8.10 Матриця планування ДФЕ для семи факторів
нів першого порядку зводиться до виконання таких процедур.
1. Вибір області планування. Вибір центру плану , тобто значення базових рівнів , здійснюється за формулою (8.4). Крок варіювання по кожному фактору визначається за формулою (8.5).
Оскільки вихідна інформація, за допомогою якої здійснюється вибір інтервалів варіювання, як правило, має обмежений характер і у більшості випадків є якісною, то по мірі проведення дослідів і опрацювання отриманих результатів може виявитися необхідність в коригуванні значень . У первісному варіанті доцільно вибирати припустимого діапазону зміни фактора . Якщо є переконливі підстави
вважати, що лінійна модель (8.7) або модель у вигляді (8.8) можуть статися недієздатними і необхідно перейти до планування другого порядку, доцільніше вибрати дещо меншим: цього діапазо-
ну. Визначені умови планування представляються у вигляді таблиці.
Незалежні
Рівні варіювання
Інтервали
390
Розділ 8
змінні
нижній –1
основний 0
верхній +1
варіювання
2. Складання плану експерименту у цілому і вибір елементів матриці дублювання U. План експерименту складається у відповідності з п. 8.3.1.
В плануванні першого порядку застосовується рівномірне дублювання дослідів, тобто, всі діагональні елементи дорівнюють один одному: . Це означає, що кожна комбінація значень факторів із певного рядка матриці плану повторюється в ході експерименту разів. Загальна кількість дослідів при цьому дорівнює . Зазвичай .
3. Проведення експерименту. Експеримент реалізується окре-
мими серіями. Кожна серія включає в себе проведення дослідів, що відповідають всім рядкам матриці плану. Кількість серій, очевидно, дорівнює . Усередині кожної серії порядок реалізації дослідів повинен бути випадковим, для цього застосовується процедура рандомізації: -кратне проведення даної процедури забезпечує різну (випадкову) послідовність реалізації рядків матриці плану в кожній серії дослідів. Безумовно, можливе і випадкове співпадіння цих послідовностей, особливо при малих .
Кінцевим результатом проведення експерименту є таблиця значень результативного показника .
Номер досліду
Серія дослідів
1
2
Дані таблиці являють собою вихідний матеріал для подальшого опрацювання.
Дослідження систем методом планування експериментів
391
8.4.3. Опрацювання експериментальних даних. Отримані дані оп-
рацьовуються у відповідності з загальною схемою регресійного аналізу. Однак специфічна форма використовуваних планів дає змогу знайти частинні варіанти загальних співвідношень регресійного аналізу, які є більш простими і зручними при розрахунках.
Алгоритм опрацювання даних включає в себе такі основні етапи. 1. Перевірка відтворюваності експерименту. Ця перевірка
здійснюється за результатами паралельних дослідів. Згідно з вимогами регресійного аналізу, правдиве опрацювання і використання результатів експериментальних досліджень можливе тільки у тому разі, коли дисперсії вимірювань функції відгуку в кожній точці досліду є однаковими. Така властивість називається однорідністю дисперсій. Вважається, що виконується, якщо справедлива гіпотеза : . Перевірка даної гіпотези за альтернативної : «хоча б одна дисперсія не дорівнює решті», здійснюється за допомогою критерію Кохрена (додаток Д9), оскільки кількість паралельних дослідів є однаковою (дорівнює ).
Дисперсія для -го досліду визначається за формулою
, (8.15)
де – -те вимірювання функції відгуку у -й точці плану;
– середнє арифметичне значення функції відгуку
в -му досліді.
Гіпотезу про однорідність вибіркових дисперсій за критерієм Кохрена перевіряють за такою схемою:
–серед обчислених за формулою (8.15) дисперсій знаходять найбільшу ;
–обчислюють значення критерію Кохрена
392
Розділ 8
;
(8.16)
–визначають кількість ступенів вільності: , ;
–вибирають рівень значимості (зазвичай );
–за значеннями , та знаходять величину табличного значення .
–порівнюють величини і . Якщо , тоді гіпотеза про однорідність дисперсій приймається. У протилежному випадку
() гіпотеза відкидається.
Якщо перевірка однорідності дисперсії дає негативний результат (гіпотеза про однорідність дисперсій відкидається), то здобутий емпіричний матеріал не рекомендується використовувати для апроксимації функції відгуку поліномами. Слід повторити експерименти, збільшуючи при цьому кількість паралельних дослідів.
сії, завдяки ортогональності матриці планування, обчислюються незалежно один від одного за формулами:
;
;
, (8.18)
де − середнє арифметичне функції відгуку з паралельних дослідів в -порядку матриці планування.
4. Перевірка значимості коефіцієнтів регресії. Значення деяких коефіцієнтів можуть дорівнювати нулю. Упевнитися в цьому можна за
Дослідження систем методом планування експериментів
393
допомогою оцінок коефіцієнтів регресії, перевіряючи їх значимість (перевірка нуль-гіпотези).
Значимість коефіцієнтів лінійної регресії перевіряють окремо для кожного коефіцієнта за допомогою критерію Стьюдента. Схема перевірки включає такі операції:
– знаходять похибку обчислення коефіцієнтів регресії (постійна для всіх коефіцієнтів) за формулою
;
(8.19)
– обчислюють статистику
;
(8.20)
–знаходять кількість ступенів вільності ;
–вибирають рівень значимості (зазвичай );
–у додатку Д6 для вибраних значень і знаходять табличне значення критерію Стьюдента ;
–порівнюють статистику з табличним значенням критерію . Якщо , відповідний коефіцієнт вважається значимим, у протилежному випадку коефіцієнт вважається статистично незначимим.
Статистична незначимість коефіцієнтів регресії може бути обумовленою декількома причинами, а саме: фактор, який відповідає незначимому фактору, не впливає на функцію відгуку; наявна велика похибка при визначенні функції відгуку; вибрано малий крок варіювання незалежної змінної; екстремум функції відгуку за відповідною змінною знаходиться близько до центру планування.
Якщо деякий коефіцієнт є незначимим, він виключається з рівняння регресії без перерахунку всіх інших коефіцієнтів, так як в плануванні лінійної моделі всі коефіцієнти є незалежними.
Кінцевим результатом після перевірки значимості всіх коефіцієнтів регресії, що входять в шукану модель, є рівняння регресії, що містить тільки змінні із значимими коефіцієнтами:
394
Розділ 8
,
де − загальна кількість значимих коефіцієнтів.
5. Перевірка адекватності моделі. Опис функції відгуку апрок-
симуючими поліномами, коефіцієнти яких визначені за методом найменших квадратів, може й не відповідати (бути неадекватним) спостережуваним значенням. Тому завжди, перш ніж використовувати математичну модель для дослідження виробничих систем, необхідно переконатися у її адекватності даним експерименту.
Гіпотеза адекватності моделі перевіряється оцінюванням відхилень значень функції , обчислених за рівнянням регресії, від експериментально встановлених . Для оцінювання відхилень використовується критерій Фішера. Перевірка адекватності виконується у такій послідовності:
–обчислюють статистичну оцінку дисперсії адекватності за формулою
;(8.21)
–знаходять розрахункове значення -критерію Фішера
;
(8.22)
– визначають кількість ступенів вільності і за формулами
, ;
–вибирають рівень значимості (зазвичай );
–у додатку Д8 для визначених , та знаходять табличне значення критерію Фішера ;
–якщо обчислене значення , то математичний опис функції відгуку прийнятим рівнянням регресії вважається адекватним.
Дослідження систем методом планування експериментів
395
Для перевірки адекватності необхідно, щоб виконувалась умова , тобто кількість дослідів повинна перевищувати кількість членів апроксимуючого полінома.
6. Перевірка дієздатності регресійної моделі
Як і в стандартній схемі регресійного аналізу, модель можна вважати дієздатною (придатною до практичного використання для цілей провіщання), якщо у цієї моделі коефіцієнт детермінації .
У цьому випадку похибка передбачення зменшується вдвічі, у порівнянні з передбаченням за середнім значенням функції відгуку без врахування впливу факторів на відгук .
Для отриманої регресійної моделі коефіцієнт детермінації розраховується за формулою:
(8.23)
Отримана за допомогою планів першого порядку дієздатна регресійна модель може бути використаною надалі для цілей передбачення значень відгуку в задачах інтерполяції і екстраполяції, а також для кращого з’ясування механізму явищ, що відбуваються в об’єкті дослідження.
Приклад 3. Досліджується процес централізованої доставки вантажів з центральної бази постачання клієнтам. За вихідний параметр (функція відгуку) прийнятий час , необхідний для обслуговування одного клієнта. Із незалежних параметрів, що впливають на процес доставки вантажів, апріорі вибрані такі:
–годинна інтенсивність надходження замовлень на транспортне обслуговування , од/год;
–кількість клієнтів, що потребують обслуговування ;
–відстань перевезення , км.
396
Розділ 8
Концептуальну модель досліджуваного процесу можна представити у вигляді функціональної залежності між структурними компонентами . Задача полягає у побудові емпіричної залежності між досліджуваними оцінками.
Розв’язок.
1. Встановлюємо область планування у вигляді основних характеристик незалежних змінних (табл. 8.11).
Таблиця 8.11
Характеристики незалежних змінних
Незалежна
Позна-
Рівні варіювання
Інтервал
нижній
основний
верхній
змінна
чення
варіювання
–1
0
+1
Інтенсивність над-
0,30
0,65
1,0
0,35
ходження, од/год.
Кількість клієнтів
3
4
5
1
Відстань перевезень, км.
10
55
100
45
2. Виконуємо кодування
змінних:
;
;
.
3. Спочатку для розв'язання задачі приймаємо лінійну модель, яка в кодованих змінних має вигляд
.
4. Для отримання у явному вигляді рівняння регресії застосовуємо ПФЕ типу з трьома () паралельними дослідами. Складаємо робочу таблицю планування (табл. 8.12) і у відповідністю з нею проводимо експеримент. Значення функції відгуку отримані в результаті імітаційного моделювання.
Таблиця 8.12
Матриця планування ПФЕ і результати експерименту
Номер досліду
1
2
3
4
5
6
7
Кодовані незалежні
Числові значення в точках плану
змінні
Незалежні змінні
Функція відгуку
+1
–1
–1
–1
0,3
3
10
5,95
6,35
6,15
6,15
0,0400
+1
+1
–1
–1
1,0
3
10
6,80
6,90
6,70
6,80
0,0100
+1
–1
+1
–1
0,3
5
10
8,88
8,83
8,75
8,81
0,0045
+1
+1
+1
–1
1,0
5
10
9,16
9,18
9,32
9,22
0,0152
+1
–1
–1
+1
0,3
3
100
3,84
3,68
3,82
3,78
0,0152
+1
+1
–1
+1
1,0
3
100
4,57
4,38
4,49
4,48
0,0182
+1
–1
+1
+1
0,3
5
100
7,82
7,90
7,83
7,85
0,0038
Дослідження систем методом планування експериментів
397
8
+1
+1
+1
+1
1,0
5
100
7,83
7,73
8,08
7,88
0,0330
5. Перевіряємо відтворюваність дослідів. За формулою (8.16) обчислюємо статистику
.
Табличне значення критерію Кохрена (додаток Д9) для рівня значимості , і ступенів вільності і дорівнює . Оскільки виконується умова , то досліди вважаються відтворюваними, тобто маємо однорідну дисперсію.
6. Обчислюємо коефіцієнти полінома за формулами (8.18):
7. Перевіряємо значимість коефіцієнтів регресії, для чого обчислюємо:
– оцінку дисперсії відтворюваності за формулою (8.17)
;
– похибку обчислення коефіцієнтів рівняння регресії за формулою (8.19)
;
– статистику для кожного коефіцієнта рівняння регресії за формулою (8.20):