Для перевірки адекватності моделі розраховують число степенів вільності
, (2.8)
де число коефіцієнтів математичної моделі.
Далі знаходять залишкову дисперсію, або дисперсію адекватності
, (2.9)
де – кількість повторних дослідів параметра оптимізації, , – результати дослідів, середнє значення результатів дослідів.
Для перевірки гіпотези про адекватність використовується -критерій Фішера
, (2.10)
де дисперсія відтворювання:
, (2.11)
Перевірка значущості кожного коефіцієнта виконується незалежно. Спочатку розраховується дисперсія коефіцієнтів
. (2.12)
Із цієї формули видно, що дисперсії всіх коефіцієнтів однакові, оскільки вони залежать тільки від помилки та кількості дослідів. Знаходять надійний інтервал
(2.13)
де табличне значення критерію Стьюдента, при числі степенів вільності, з якими знаходили табл.(Д.2)
, (2.14)
Якщо абсолютна величина коефіцієнта більша ніж величина надійного інтервалу, то коефіцієнт суттєвий.
Прийняття рішень після побудови математичної моделі
Адекватна математична модель, яка побудована, має вигляд полінома першого степеня. Коефіцієнти моделі – це частинні похідні функції відклику або тангенси кутів нахилу гіперплощини до відповідної осі. Більший за абсолютною величиною коефіцієнт відповідає більшому куту нахилу та більшому впливу на параметр оптимізації. Величина коефіцієнта регресії – це кількісна міра цього впливу. Знак коефіцієнта свідчить про характер впливу. Якщо знаходять максимум функції впливу, то необхідне зростання величини коефіцієнтів, знак яких від’ємний та навпаки: коли необхідно знайти максимум функції впливу, необхідне зростання величин коефіцієнтів, що мають додатний знак. У деяких задачах важливо побудувати рівняння регресії для натуральних значень факторів. У цьому випадку, використовуючи формулу (2.4), можна отримати рівняння для натуральних значень факторів, але слід мати на увазі, що коефіцієнти регресії такого рівняння будуть залежати один від одного.
Практичне завдання № 3
Задача. Досліджується вплив температурних і технологічних параметрів на процес виробництва деякої речовини. Вивчення процесу дозволило виділити де-кілька факторів, які суттєво впливають на процес, що вивчається:
- частота обертання механізму від 300 до 1100 хв-1 збільшує вихід речовини (, хв.-1);
- зростання частки каталізатора збільшує вихід речовини (, л/год);
- зростання температури від 20 до 60 0С зменшує вихід речовини ( ,0С);
- зростання концентрації одного інгредієнта від 50 до 100 г/л збільшує вихід речовини (, г/л).
Як вихідна змінна прийнята частка виходу речовини. Необхідно знайти такі значення факторів, які забезпечать максимальний вихід речовини. Обрана така область зміни вказаних параметрів:
400 < < 1000 хв-1; 0.5 < < 3,5 л\ч);
10 < < 60 0С; 75 < < 150 г/л.
Усі необхідні дані зведені в табл. 2.3 та 2.4
Таблица 2.3
Вихідні дані для планування експерименту
Умови досліду | Позначення | Фактори | |||
, хв-1 | ,л\год | ,0С | ,г/л | ||
Нульовий рівень Інтервал варіювання Верхній рівень =+1 Нижній рівень = - 1 |
Таблиця 2.4
План та результати експерименту
Номер досліду | План експерименту | Вихідна змінна | ||||||
1 +1 +1 +1 +1 +1 90,2 93,5 91,85
2 +1 -1 +1 +1 +1 86,0 87,5 86,75
3 +1 +1 -1 +1 +1 85,0 86,0 85,50
4 +1 -1 -1 +1 +1 81,0 82,0 81,50
5 +1 +1 +1 -1 +1 99,0 98,5 98,75
6 +1 -1 +1 -1 +1 96,0 94,0 95,0
7 +1 +1 -1 -1 +1 97,0 97,5 97,25
8 +1 -1 -1 -1 +1 90,2 90,0 90,1
9 +1 +1 +1 +1 -1 71,0 73,0 72,0
10 +1 -1 +1 +1 -1 66,0 67,5 66,75
11 +1 +1 -1 +1 -1 64,5 63,5 64,0
12 +1 -1 -1 +1 -1 52,0 54,0 53,0
13 +1 +1 +1 -1 -1 89,0 88,0 88,5
14 +1 -1 +1 -1 -1 80,4 80,0 80,2
15 +1 +1 -1 -1 -1 85,0 82,0 83,5
16 +1 -1 -1 -1 -1 76,5 77,0 76,5
Завдання:
1. Розрахувати коефіцієнти регресійного рівняння (2.6) для повного факторного експерименту, виконати порівняльний аналіз розрахункових значень параметрів оптимізації з експериментом, зробити перевірку адекватності математичної моделі та значущості кожного коефіцієнта регресійного рівняння. Виконати, згідно з розд. 2.4, дослідження математичної моделі та знайти область з максимальним виходом продукту.
2. Побудувати план дробового факторного експерименту 24-1 для випадку лінійної моделі та моделі із взаємодією для таких пів реплік, які задані визначальними контрастами та наведені нижче у табл. 2.5 (вибір варіанта узгоджується з викладачем).
3. Побудувати математичну модель, обчислити коефіцієнти регресійного рівняння та виконати порівняльний аналіз обчислюваних значень параметрів оптимізації з експериментом.
Таблиця 2.5
Варіанти завдань
N варіанта | Визначальний контраст | N варіанта | Визначальний контраст |
1.0 | 2.0 | ||
1.1 | 2.1 | ||
1.2 | 2.2 | ||
1.3 | 2.3 | ||
1.4 | 2.4 | ||
1.5 | 2.5 | ||
1.6 | 2.6 | ||
1.7 | 2.7 | ||
1.8 | 2.8 |
ДОДАТОК
Таблиця Д.1
Таблица значений функций
0,00 | 0,45 | 0,1736 | 0,90 | 0,3159 | 1,35 | 0,4115 | |
0,01 | 0,0040 | 0,46 | 0,1772 | 0,91 | 0,3186 | 1,36 | 0,4131 ' |
0,02 | 0,0080 | 0,47 | 0,45 | 0,92 | 0,3212 | 1,37 | 0,4147. |
0,03 | 0,0120 | 0.48 | 0,46 | 0,93 | 0,3238 | 1,38 | 0,4162 |
0,04 | 0,0160 | 0,49 | 0,47 | 0,94 | 0,3264 | 1,39 | 0,4177 |
0,05 | 0,0199 | 0,50 | 0.48 | 0,95 | 0,3289 | 1,40 | 0,4192 |
0,06 | 0,0239 | 0,51 | 0,49 | 0,96 | 0,3315 | 1,41 | 0,4207 |
0,07 | 0,0279 | 0,52 | 0,50 | 0,97 | 0,3340 | 1,42 | 0,4222 |
0,08 | 0,00319 | 0,53 | 0,51 | 0,98 | 0,3365 | 1,43 | 0,4236 |
0,09 | 0,0359 | 0,54 | 0,52 | 0,99 | 0,3389 | 1,44 | 0,4251 |
0,10 | 0,0398 | 0,55 | 0,53 | 1,00 | 0,34Ш | 1,45 | 0,4265 |
0,11 | 0,0438 | 0,56 | 0,54 | 1,01 | 0,3438 | 1,46 | 0,4279 |
0,12 | 0,0478 | 0,57 | 0,55 | 1,02 | 0,3461 | 1,47 | 0,4292 |
0,13 | 0,0517 | 0,58 | 0,56 | 1,03 | 0,3485 | 1,48 | 0,4306 |
0,14 | 0,0557 | 0,59 | 0,57 | 1,04 | 0,3508 | 1,49 | 0,4319 |
0;15 | 0,0596 | 0,60 | 0,58 | 1,05 | 0,3531 | 1,50 | 0,4332 |
0,16 | 0,0636 | 0,61 | 0,59 | 1,06 | 0,3554 | 1,51 | 0,4345 |
0,17 | 0,0675 | 0,62 | 0,60 | 1,07 | 0,3577 | 1,52 | 0,4357 |
0,18 | 0,0714 | 0,63 | 0,61 | 1,08 | 0,3599 | 1,53 | 0, 4370 |
0,19 | 0,0753 | 0,64 | 0,62 | 1,09 | 0,3621 | 1,54 | 0,4382 |
0,20 | 0,0793 | 0,65 | 0,63 | 1,10 | 0,3643 | 1,55 | 0,4394 |
0,21 | 0,0832 | 0,66 | 0,64 | 1,11 | 0,3665 | 1,56 | 0,4406 |
0,22 | 0,0871 | 0,67 | 0,65 | 1,12 | 0,3686 | 1,57 | 0,4418 |
0,23 | 0,0910 | 0,68 | 0,66 | 1,13 | 0,3708 | 1,58 | 0,4429 |
0,24 | 0,0948 | 0,69 | 0,67 | 1,14 | 0,3729 | 1,59 | 0,4441 |
0,25 | 0,0987 | 0,70 | 0,68 | 1,15 | 0,3749 | 1,60 | 0,4452 |
0,26 | 0,1026 | 0,71 | 0,69 | 1,16 | 0,3770 | 1,61 | 0,4463 |
0,27 | 0,1064 | 0,72 | 0,70 | 1,17 | 0,3790 | 1,62 | 0,4474 |
0,28 | 0,1103 | 0,73 | 0,71 | 1,18 | 0,3810 | 1,63 | 0,4484 |
0,29 | 0,1141 | 0,74 | 0,72 | 1,19 | 0,3830 | 1,64 | 0,4495 |
0,30 | 0,1179 | 0,75 | 0,73 | 1,20 | 0,3849 | 1,65 | 0,4505 |
0,31 | 0,1217 | 0,76 | 0,74 | 1,21 | 0,3869 | 1,66 | 0,4515 |
0,32 | 0,1255 | 0,77 | 0,75 | 1,22 | 0,3883 | 1,67 | 0,4525 |
0,33 | 0,1293 | 0,78 | 0,76 | 1,23 | 0,3907 | 1,68 | 0,4535. |
0,34 | 0,1331 | 0,79 | 0,77 | 1,24 | 0,3925 | 1,69 | 0,4545. |
0,35 | 0,1368 | 0,80 | 0,78 | 1,25 | 0,3944 | 1,70 | 0,4554 |
0,36 | 0,1406 | 0,81 | 0,79 | 1,26 | 0,3962 | 1,71 | 0,4564 |
0,37 | 0,1443 | 0,82 | 0,80 | 1,27 | 0,3980 | 1,72 | 0,4573 |
0,38 | 0,1480 | 0,83 | 0,81 | 1,28 | 0,3997 | 1,73 | 0,4582 |
0,39 | 0,1517 | 0,84 | 0,82 | 1,29 | 0,4015 | 1,74 | 0,4591 |
0,40 | 0,1554 | 0,85 | 0,83 | 1,30 | 0,4032 | 1,75 | 0,4599 |
0,41 | 0,1591 | 0,86 | 0,84 | 1,31 | 0,4049 | 1,76 | 0,4608 |
0,42 | 0,1628 | 0,87 | 0,85 | 1,32 | 0,4066 | 1,77 | 0,4616 |
0,43 | 0,1664 | 0,88 | 0,86 | 1,33 | 0,4082 | 1,78 | 0,4625 |
0,44 | 0,1700 | 0,89 | 0,87 | 1,34 | 0,4099 | 1,79 | 0,4633 |
Продолжение таблицы
1,80 | 0,46410 | 2,00 | 0,4772 | 2,40 | 0,4918 | 2,80 | 0,4974 |
1,81 | 0,4649 | 2,02 | 0,4783 | 2,42 | 0,4922 | 2,82 | 0,4976 |
1,82 | 0,4656 | 2,04 | 0,4793 | 2,44 | 0,4927 | 2,84 | 0,4977 |
1,83 | 0,4664 | 2,06 | 0,4803 | 2,46 | 0,4931 | 2,86 | 0,4979 |
1,84 | 0,4671 | 2,08 | 0,4812 | 2,48 | 0,4934 | 2,88 | 0,4980 |
1,85 | 0,4678 | 2,10 | 0,4821 | 2,50 | 0,4938 | 2,90 | 0,4981 |
1,86 | 0,4686 | 2,12 | 0,4830 | 2,52 | 0,4941 | 2,92 | 0,4982 |
1,87 | 0,4693 | 2,14 | 0,4838 | 2,54 | 0,4945 | 2,94 | 0,4984 |
1,88 | 0,4699 | 2,16 | 0,4846 | 2,56 | 0,4948 | 2,96 | 0,4985 |
1,89 | 0,4706 | 2,18 | 0,4854 | 2,58 | 0,4951 | 2,98 | 0,4986 |
1,90 | 0,4713 | 2,20 | 0,4861 | 2,60 | 0,4953 | 3,00 | 0,49865 |
1,91 | 0,4719 | 2,22 | 0,4868 | 2,62 | 0,4956 | 3,20 | 0,49931 |
1,92 | 0,4726 | 2,24 | 0,4875 | 2,64 | 0,4959 | 3,40 | 0,49966 |
1,93 | 0,4732 | 2,26 | 0,4881 | 2,66 | 0,4961 | 3,60 | 0,499841 |
1,94 | 0,4738 | 2,28 | 0,4887 | 2,68 | 0,4963 | 3,80 | 0,499928 |
1,95 | 0,4744 | 2,30 | 0,4893 | 2,70 | 0,4965 | 4,00 | 0,499968 |
1,96 | 0,4750 | 2,32 | 0,4898 | 2,72 | 0,4967 | 4,50 | 0,499997 |
1,97 | 0,4756 | 2,34 | 0,4904 | 2,74 | 0,4969 | 5,00 | 0,499997 |
1,98 | 0,4761 | 2,36 | 0,4909 | 2,76 | 0,4971 | ||
1.99 | 0,4767 | 2,38 | 0,4913 | 2,78 | 0,4973 |
Таблиця Д.2