Відмінності між версіями «Планування другого порядку»

 
(Не показані 48 проміжних версій 8 користувачів)
Рядок 1: Рядок 1:
{{Завдання|ihor_p|Назаревич О.Б.|05 березня 2010}}
+
{{Невідредаговано}}
 +
{{Студент | Name=Ігор | Surname=Пельц | FatherNAme= |Faculti=ФІС | Group=СНм-51 | Zalbook=}}
  
http://elartu.tstu.edu.ua/handle/123456789/404 Презентація доповіді (університетський репозиторій).
+
http://elartu.tstu.edu.ua/handle/123456789/415 Презентація доповіді (університетський репозиторій).
  
 
= Планування другого порядку =
 
= Планування другого порядку =
Рядок 7: Рядок 8:
 
Планування другого порядку застосовується для математичного опису об'єкта поблизу екстремальної точки статистичної характеристики або тоді, коли необхідний точніший опис в інших точках факторного простору. При цьому використовують  поліном  другого  порядку:
 
Планування другого порядку застосовується для математичного опису об'єкта поблизу екстремальної точки статистичної характеристики або тоді, коли необхідний точніший опис в інших точках факторного простору. При цьому використовують  поліном  другого  порядку:
  
<math>y={{a}_{0}}+\sum\limits_{i=1}^{n}{{{a}_{i}}{{x}_{i}}+}\sum\limits_{\begin{smallmatrix}
+
<center>
i=1 \\
+
<math>\operatorname{y}={{a}_{0}}+\sum\limits_{i=1}^{n}{{{a}_{i}}{{x}_{i}}+}\sum\limits_{i=1}^{n}{{{a}_{ij}}{{x}_{i}}{{x}_{j}}+}...+\sum\limits_{i=1}^{n}{{{a}_{i,i}}x_{i}^{2}}</math>
(j>i)
+
</center>
\end{smallmatrix}}^{n}{{{a}_{ij}}{{x}_{i}}{{x}_{j}}+}...+\sum\limits_{i=1}^{n}{{{a}_{i,i}}x_{i}^{2}}</math>
 
 
 
<center>[[Файл:Formula.jpg]]</center>
 
  
 
Задача, як і в ПФЕ, полягає у визначенні методом найменших квадратів за результатами спланованого експерименту коефіцієнтів цього рівня за умови, що виконуються передумови  регресійного аналізу.
 
Задача, як і в ПФЕ, полягає у визначенні методом найменших квадратів за результатами спланованого експерименту коефіцієнтів цього рівня за умови, що виконуються передумови  регресійного аналізу.
Рядок 18: Рядок 16:
  
 
<center>
 
<center>
 
+
<math>\sum\limits_{g=1}^{N}{{{z}_{gi}}{{z}_{gj}}=0;i\ne j.}</math>
 
</center>
 
</center>
  
Рядок 33: Рядок 31:
 
Ортогональність матриці композиційного планування забезпечується  виконанням  рівностей
 
Ортогональність матриці композиційного планування забезпечується  виконанням  рівностей
  
<center>[[Файл:formula3.gif]]</center>
+
<center>
  
де і-номер фактора;  j-номер  рядка; g-номер досліду. Для ортогоналізації першого з цих співвідношень застосовується перехід до нової змінної
+
<math>\sum\limits_{g=1}^{N}{z_{g,o}^{2}z_{y,i}^{2}=0,\sum\limits_{g=1}^{N}{z_{gi}^{2}z_{gj}^{2}=0,i\ne j}.}</math>
  
<center>[[Файл:formula4.gif]]</center>
+
</center>
  
 +
де і-номер фактора;  j-номер  рядка; g-номер досліду. Для ортогоналізації першого з цих співвідношень застосовується перехід до нової змінної
 +
<center>
 +
<math>\overset{-}{z}\,_{i}^{2}=z_{i}^{2}-\frac{\sum{z_{gi}^{2}}}{N}=z_{i}^{2}-\overset{-}{z}\,_{i}^{2}</math>
 +
</center>
 
де N-загальне число експериментів.  
 
де N-загальне число експериментів.  
 
Величина zi залежить тільки від числа факторів n і числа дослідів N, яке звичайно вибирається так, що
 
Величина zi залежить тільки від числа факторів n і числа дослідів N, яке звичайно вибирається так, що
  
<center>[[Файл:formula5.gif]]</center>
+
<center>
 +
 
 +
<math>\operatorname{N}={{N}_{n}}+{{N}_{a}}+{{N}_{0}}={{2}^{n}}+2n+1,</math>
 +
 
 +
</center>
  
 
де Nn = <math>2^n</math>-кількість вершин гіперкуба при ПФЕ; Nα = 2n-число зіркових точок; N0= 1- число дослідів у центрі плану. Якщо N вибрати так, то zi визначають за формулою
 
де Nn = <math>2^n</math>-кількість вершин гіперкуба при ПФЕ; Nα = 2n-число зіркових точок; N0= 1- число дослідів у центрі плану. Якщо N вибрати так, то zi визначають за формулою
  
<center>[[Файл:formula6.gif]]</center>
+
<center>
 +
 
 +
<math>\operatorname{z}_{i}^{2}=\frac{{{2}^{n}}+2{{\alpha }^{2}}}{N}.</math>
 +
 
 +
</center>
  
 
Тому ортогоналізація другої з вищенаведених умов досягається вибором бажаного α.
 
Тому ортогоналізація другої з вищенаведених умов досягається вибором бажаного α.
 
Для зручності підготовки і планування величини α, N, Nn, Nα обчислені і табульовані залежно від числа факторів (табл.  1).
 
Для зручності підготовки і планування величини α, N, Nn, Nα обчислені і табульовані залежно від числа факторів (табл.  1).
<center>Таблиця1 - Підготовка ОЦКП другого порядку</center>
+
 
<center>[[Файл:Табл1.gif]]</center>
+
 
 +
<table width="90%" border="1">
 +
  <caption>
 +
    Таблиця1 - Підготовка ОЦКП другого порядку.
 +
  </caption>
 +
  <tr>
 +
    <th scope="col">n&nbsp;</th>
 +
    <th scope="col"><math>{{\operatorname{N}}_{0}}</math>&nbsp;</th>
 +
    <th scope="col"><math>{{\operatorname{N}}_{n}}</math>&nbsp;</th>
 +
    <th scope="col"><math>{{\operatorname{N}}_{\alpha }}</math>&nbsp;</th>
 +
    <th scope="col">N&nbsp;</th>
 +
    <th scope="col"><math>\alpha </math>&nbsp;</th>
 +
  </tr>
 +
  <tr>
 +
    <td>2&nbsp;</td>
 +
    <td>1&nbsp;</td>
 +
    <td>4&nbsp;</td>
 +
    <td>2&nbsp;</td>
 +
    <td>9&nbsp;</td>
 +
    <td>1,000&nbsp;</td>
 +
  </tr>
 +
  <tr>
 +
    <td>3&nbsp;</td>
 +
    <td>1&nbsp;</td>
 +
    <td>8&nbsp;</td>
 +
    <td>6&nbsp;</td>
 +
    <td>15&nbsp;</td>
 +
    <td>1,215&nbsp;</td>
 +
</tr>
 +
<tr>
 +
    <td>5&nbsp;</td>
 +
    <td>1&nbsp;</td>
 +
    <td>16&nbsp;</td>
 +
    <td>8&nbsp;</td>
 +
    <td>25&nbsp;</td>
 +
    <td>1,414&nbsp;</td>
 +
</tr>
 +
<tr>
 +
    <td>5&nbsp;</td>
 +
    <td>1&nbsp;</td>
 +
    <td>32&nbsp;</td>
 +
    <td>10&nbsp;</td>
 +
    <td>43&nbsp;</td>
 +
    <td>1,596&nbsp;</td>
 +
</tr>
 +
<tr>
 +
    <td>6&nbsp;</td>
 +
    <td>1&nbsp;</td>
 +
    <td>64&nbsp;</td>
 +
    <td>12&nbsp;</td>
 +
    <td>77&nbsp;</td>
 +
    <td>1,706&nbsp;</td>
 +
</tr>
 +
<tr>
 +
    <td>7&nbsp;</td>
 +
    <td>1&nbsp;</td>
 +
    <td>128&nbsp;</td>
 +
    <td>14&nbsp;</td>
 +
    <td>143&nbsp;</td>
 +
    <td>1,909&nbsp;</td>
 +
</tr>
 +
  </table>
  
 
Сформуємо матриці ОЦКП для двох факторів. При n = 2 отримаємо α = 1,0 (див. табл. 1), для обчислення  скористаємося наведеними формулами. Оскільки  <math>z1^2=(4+2)/9=0,667</math>, то стовпець  добувається відніманням одного і того ж числа 0,67 від числа стовпця <math>zi^2</math>  того ж рядка (табл. 2).
 
Сформуємо матриці ОЦКП для двох факторів. При n = 2 отримаємо α = 1,0 (див. табл. 1), для обчислення  скористаємося наведеними формулами. Оскільки  <math>z1^2=(4+2)/9=0,667</math>, то стовпець  добувається відніманням одного і того ж числа 0,67 від числа стовпця <math>zi^2</math>  того ж рядка (табл. 2).
  
<center>Таблиця2 - Матриця ОЦКП другого порядку для двофакторного експерименту</center>
+
<table width="90%" border="1">
<center>[[Файл:Табл2.gif]]</center>
+
  <caption>
 +
    Таблиця2 - Матриця ОЦКП другого порядку для двофакторного експерименту
 +
  </caption>
 +
  <tr>
 +
    <th scope="col">Дослід&nbsp;</th>
 +
    <th scope="col">j&nbsp;</th>
 +
    <th scope="col"><math>z_0</math>&nbsp;</th>
 +
    <th scope="col"><math>z_1</math>&nbsp;</th>
 +
    <th scope="col"><math>z_2</math>&nbsp;</th>
 +
    <th scope="col"><math>z_1^2</math>&nbsp;</th>
 +
    <th scope="col"><math>z_2^2</math>&nbsp;</th>
 +
    <th scope="col"><math>z_{1c}^2</math>&nbsp;</th>
 +
    <th scope="col"><math>z_{2c}^2</math>&nbsp;</th>
 +
    <th scope="col"><math>z_1*z_2</math>&nbsp;</th>
 +
  </tr>
 +
  <tr>
 +
    <th scope="row">Вершини квадрата&nbsp;</th>
 +
    <td><center>1</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>0,33</center>&nbsp;</td>
 +
    <td><center>0,33</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
  </tr>
 +
  <tr>
 +
    <th scope="row">&nbsp;</th>
 +
    <td><center>2</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>-</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>0,33</center>&nbsp;</td>
 +
    <td><center>0,33</center>&nbsp;</td>
 +
    <td><center>-</center>&nbsp;</td>
 +
  </tr>
 +
  <tr>
 +
    <th scope="row">&nbsp;</th>
 +
    <td><center>3</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>-</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>0,33</center>&nbsp;</td>
 +
    <td><center>0,33</center>&nbsp;</td>
 +
    <td><center>-</center>&nbsp;</td>
 +
  </tr>
 +
  <tr>
 +
    <th scope="row">&nbsp;</th>
 +
    <td><center>4</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>-</center>&nbsp;</td>
 +
    <td><center>-</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>0,33</center>&nbsp;</td>
 +
    <td><center>0,33</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
  </tr>
 +
  <tr>
 +
    <th scope="row">Зіркові точки&nbsp;</th>
 +
    <td><center>5</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>-</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0,33</center>&nbsp;</td>
 +
    <td><center>-0,67</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
  </tr>
 +
  <tr>
 +
    <th scope="row">&nbsp;</th>
 +
    <td><center>6</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0,33</center>&nbsp;</td>
 +
    <td><center>-0,67</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
  </tr>
 +
  <tr>
 +
    <th scope="row">&nbsp;</th>
 +
    <td><center>7</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>-</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>-0,67</center>&nbsp;</td>
 +
    <td><center>0,33</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
  </tr>
 +
<tr>
 +
    <th scope="row">&nbsp;</th>
 +
    <td><center>8</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>-0,67</center>&nbsp;</td>
 +
    <td><center>0,33</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
  </tr>
 +
<tr>
 +
    <th scope="row">Центр&nbsp;</th>
 +
  <td><center>9</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>-0,67</center>&nbsp;</td>
 +
    <td><center>-0,67</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
  </tr>
 +
</table>
  
 
Реалізація експериментів за ОЦКП здійснюється за тією ж методикою, що і ПФЕ. Таким чином, через випадковий характер зміни вихідної величини у у кожній точці хg проводиться m паралельних дослідів і обчислюється середнє значення функції відклику
 
Реалізація експериментів за ОЦКП здійснюється за тією ж методикою, що і ПФЕ. Таким чином, через випадковий характер зміни вихідної величини у у кожній точці хg проводиться m паралельних дослідів і обчислюється середнє значення функції відклику
  
<center>[[Файл:formula7.gif]]</center>
+
<center><math>{{\overset{-}{y}}_{g}}=\frac{\sum\limits_{d=1}^{m}{{{y}_{g}}d}}{m}</math></center>
  
 
Перед реалізацією проводиться рандомізація рядків матриці планування.
 
Перед реалізацією проводиться рандомізація рядків матриці планування.
 
Перевірка відтворюваності, як і в ПФЕ, виконується за критерієм Кохрена, після чого обчислюється оцінка дисперсії відтворюваності
 
Перевірка відтворюваності, як і в ПФЕ, виконується за критерієм Кохрена, після чого обчислюється оцінка дисперсії відтворюваності
  
<center>[[Файл:formula8.gif]]</center>
+
<center>
 +
 
 +
 
 +
<math>S_{y}^{2}=S_{vidtv}^{2}=\frac{\sum\limits_{g=1}^{N}{S_{g}^{2}}}{N}.</math>
 +
 
 +
 
 +
</center>
  
 
Методика утворення математичної моделі незначно відрізняється від методики опису результатів ПФЕ.
 
Методика утворення математичної моделі незначно відрізняється від методики опису результатів ПФЕ.
 
Коефіцієнти регресії при ОЦКП обчислюються за формулою
 
Коефіцієнти регресії при ОЦКП обчислюються за формулою
  
<center>[[Файл:formula9.gif]]</center>
+
<center><math>{{b}_{i}}=\frac{\sum\limits_{g=1}^{N}{{{z}_{gi}}{{y}_{g}}}}{\sum\limits_{g=1}^{N}{z_{gi}^{2}}}.</math></center>
  
 
У цьому плануванні оцінки дисперсій коефіцієнтів bi (точність їхнього обчислення) не однакові, оскільки не однаковий знаменник у формулі дисперсії
 
У цьому плануванні оцінки дисперсій коефіцієнтів bi (точність їхнього обчислення) не однакові, оскільки не однаковий знаменник у формулі дисперсії
  
<center>[[Файл:formula10.gif]]</center>
+
<center><math>S_{bi}^{2}=\frac{S_{y}^{2}}{m\sum\limits_{g=1}^{N}{z_{gi}^{2}}}.</math></center>
  
 
В ПФЕ знаменник однаковий і дорівнює mN. Це істотний недолік ОЦКП, і ''тому часто надають перевагу складнішому за обчислювальними процедурами рототабельному плануванню''.
 
В ПФЕ знаменник однаковий і дорівнює mN. Це істотний недолік ОЦКП, і ''тому часто надають перевагу складнішому за обчислювальними процедурами рототабельному плануванню''.
Рядок 83: Рядок 283:
 
Щоб зробити план другого порядку рототабельним, вибирають для сфери, на якій розташовуються зіркові точки, радіус (зіркове плече) за формулою
 
Щоб зробити план другого порядку рототабельним, вибирають для сфери, на якій розташовуються зіркові точки, радіус (зіркове плече) за формулою
  
<center>[[Файл:formula4.jpg]]</center>
+
<center><math>\alpha ={{2}^{n/2}}.</math></center>
  
 
Інша умова рототабельності — збільшення числа дослідів на поверхні нульової сфери, тобто в центрі плану. У зв'язку з цим виникає повна назва методу: ''центральне композиційне  рототабельне  планування''  (ЦКРП).
 
Інша умова рототабельності — збільшення числа дослідів на поверхні нульової сфери, тобто в центрі плану. У зв'язку з цим виникає повна назва методу: ''центральне композиційне  рототабельне  планування''  (ЦКРП).
Рядок 89: Рядок 289:
 
Як і для ОЦКП, основні характеристики матриць рототабельного планування табульовані (табл. 3). Позначення тут ті самі, що і для ОЦКП (див. табл. 2). При ЦКРП,  починаючи  з  n = 5,  можна  застосувати  ДФЕ(дробовий факторний експеримент).
 
Як і для ОЦКП, основні характеристики матриць рототабельного планування табульовані (табл. 3). Позначення тут ті самі, що і для ОЦКП (див. табл. 2). При ЦКРП,  починаючи  з  n = 5,  можна  застосувати  ДФЕ(дробовий факторний експеримент).
  
<center>Таблиця 3 – Підготовка ЦКРП другого порядку</center>
+
<table width="90%" border="1">
<center>[[Файл:ТаблПід.jpg]]</center>
+
  <caption>
 +
    Таблиця 3 – Підготовка ЦКРП другого порядку
 +
  </caption>
 +
  <tr>
 +
    <th scope="col">n&nbsp;</th>
 +
    <th scope="col"><math>{{\operatorname{N}}_{n}}</math>&nbsp;</th>
 +
    <th scope="col"><math>{{\operatorname{N}}_{\alpha }}</math>&nbsp;</th>
 +
    <th scope="col"><math>{{\operatorname{N}}_{0}}</math>&nbsp;</th>
 +
    <th scope="col">N&nbsp;</th>
 +
    <th scope="col"><math>\alpha </math>&nbsp;</th>
 +
  </tr>
 +
  <tr>
 +
    <td>2&nbsp;</td>
 +
    <td>5&nbsp;</td>
 +
    <td>4&nbsp;</td>
 +
    <td>4&nbsp;</td>
 +
    <td>13&nbsp;</td>
 +
    <td>1,414&nbsp;</td>
 +
  </tr>
 +
  <tr>
 +
    <td>3&nbsp;</td>
 +
    <td>6&nbsp;</td>
 +
    <td>8&nbsp;</td>
 +
    <td>6&nbsp;</td>
 +
    <td>20&nbsp;</td>
 +
    <td>1,680&nbsp;</td>
 +
</tr>
 +
<tr>
 +
    <td>4&nbsp;</td>
 +
    <td>7&nbsp;</td>
 +
    <td>16&nbsp;</td>
 +
    <td>8&nbsp;</td>
 +
    <td>31&nbsp;</td>
 +
    <td>2,000&nbsp;</td>
 +
</tr>
 +
<tr>
 +
    <td>5&nbsp;</td>
 +
    <td>10&nbsp;</td>
 +
    <td>32&nbsp;</td>
 +
    <td>10&nbsp;</td>
 +
    <td>52&nbsp;</td>
 +
    <td>2,378&nbsp;</td>
 +
</tr>
 +
<tr>
 +
    <td>6&nbsp;</td>
 +
    <td>15&nbsp;</td>
 +
    <td>64&nbsp;</td>
 +
    <td>12&nbsp;</td>
 +
    <td>91&nbsp;</td>
 +
    <td>1,828&nbsp;</td>
 +
</tr>
 +
<tr>
 +
    <td>7&nbsp;</td>
 +
    <td>21&nbsp;</td>
 +
    <td>128&nbsp;</td>
 +
    <td>14&nbsp;</td>
 +
    <td>163&nbsp;</td>
 +
    <td>1,333&nbsp;</td>
 +
</tr>
 +
  </table>
  
 
При рототабельному плануванні для обчислення коефіцієнтів моделі і відповідних оцінок дисперсій знаходять спеціальні комплекси:
 
При рототабельному плануванні для обчислення коефіцієнтів моделі і відповідних оцінок дисперсій знаходять спеціальні комплекси:
  
<center>[[Файл:formula5.jpg]]</center>
+
<center><math>\begin{align}
 +
  & B=\frac{nN}{(n+2)(N-{{N}_{0}})}; \\
 +
& A=\frac{1}{2B[(n+2)B-n]}; \\
 +
& C=\frac{N}{N-{{N}_{0}}}, \\
 +
\end{align}</math></center>
  
 
де n-число факторів; N-загальне число дослідів у плануванні; N0-число дослідів у центрі плану.
 
де n-число факторів; N-загальне число дослідів у плануванні; N0-число дослідів у центрі плану.
 
За результатами експериментів обчислюють такі суми:
 
За результатами експериментів обчислюють такі суми:
  
<center>[[Файл:formula6.jpg]]</center>
+
<center><math>\begin{align}
<center>[[Файл:formula7.jpg]]</center>
+
  & {{S}_{0}}=\sum\limits_{g=1}^{N}{{{y}_{g}}}; \\
 +
& {{S}_{i}}=\sum\limits_{g=1}^{N}{{{y}_{g}}}{{z}_{gi}};i=1,2,...,n; \\
 +
& {{S}_{ik}}=\sum\limits_{g=1}^{N}{{{z}_{gi}}{{z}_{gk}}{{y}_{g}}};i\ne k; \\
 +
& {{S}_{ii}}=\sum\limits_{g=1}^{N}{z_{gi}^{2}{{y}_{g}}};i=1,2,...,n. \\
 +
\end{align}</math></center>
 +
 
  
 
Коефіцієнти  моделі  тут  розраховують  за  формулами
 
Коефіцієнти  моделі  тут  розраховують  за  формулами
  
<center>[[Файл:formula8.jpg]]</center>
+
<center><math>\begin{align}
 +
  & {{b}_{0}}=\frac{2AB}{N}[{{S}_{0}}B(n+2)-C\sum{{{S}_{ii}}}]; \\
 +
& {{b}_{i}}=\frac{C{{S}_{i}}}{N}; \\
 +
& {{b}_{ik}}=\frac{{{C}^{2}}{{S}_{ik}}}{BN},i\ne k; \\
 +
& {{b}_{ii}}=\frac{AC}{N}\{{{S}_{ii}}[B(n+2)-n]+C(1-B)\sum\limits_{i=1}^{n}{{{S}_{ii}}-2B{{S}_{0}}}\}. \\
 +
&  \\
 +
\end{align}</math></center>
  
 
Оцінки дисперсій для обчислених коефіцієнтів знаходять за такими формулами:
 
Оцінки дисперсій для обчислених коефіцієнтів знаходять за такими формулами:
  
<center>[[Файл:formula9.jpg]]</center>
+
<center><math>\begin{align}
 +
  & S_{b0}^{2}=\frac{2AB(n+2)}{N}S_{y}^{2}; \\
 +
& S_{bi}^{2}=\frac{S_{y}^{2}}{N-{{N}_{0}}};i=1,2,...,n; \\
 +
& S_{bik}^{2}=\frac{{{C}^{2}}S_{y}^{2}}{N},i\ne k; \\
 +
& {{S}_{bii}}=\frac{A{{C}^{2}}S_{y}^{2}}{N}[B(n+1)-(n-1)]. \\
 +
\end{align}</math></center>
  
 +
У цих формулах дисперсія відтворюваності <math>S_{y}^{2}</math> визначається за результатами дослідів у нульовій точці
 +
 +
<center><math>\begin{align}
 +
  & S_{y}^{2}=\frac{1}{{{N}_{0}}-1}\sum\limits_{g=1}^{{{N}_{0}}}{{{({{y}_{ge}}-\overset{-}{y}\,)}^{2}}}; \\
 +
& \overset{-}{y}\,=\frac{1}{{{N}_{0}}}\sum\limits_{g=1}^{{{N}_{0}}}{{{y}_{ge.}}} \\
 +
\end{align}</math></center>
 +
 +
Дисперсія адекватності оцінюється за формулою
 +
 +
<center><math>S_{adekv}^{2}=\frac{\sum\limits_{g=1}^{N}{{{({{y}_{ge}}-{{y}_{grozr}})}^{2}}-S_{y}^{2}({{N}_{0}}-1)}}{N-\frac{(n+2)(n+1)}{2}(N-1)},</math></center>
 +
 +
якшо число ступенів вільності
 +
 +
<center><math>{{f}_{adekv}}={{N}_{0}}-\frac{(n+2)(n+1)}{2}-({{N}_{0}}-1).</math></center>
 
==Приклад==
 
==Приклад==
  
 
Скласти матрицю ЦКРП на прикладі побудови математичної моделі технологічного процесу крупоутворення (див.: Пищевая технология.— 1976.— № 4.— С.  121—124).
 
Скласти матрицю ЦКРП на прикладі побудови математичної моделі технологічного процесу крупоутворення (див.: Пищевая технология.— 1976.— № 4.— С.  121—124).
'''Розв'язання'''. Як функції відклику прийнято у1, % — середня зольність крупи пшениці після перших трьох систем для дертя (швидкість обертання рифлених вальців усіх систем 6 м/с); у2, % — сумарний вихід всіх крупок, які добуваються в процесі крупоутворення; у3, кДж/(кг • %) — витрата енергії на одержання 1 % продукту з 1 кг зерна.
 
Незалежними змінними є, %: х1 — вихід крупи на першій системі для дертя; х2 — те ж на другій системі; х3 — те ж, для трьох систем для дертя. Інтервал варіювання для всіх xi, вибрано з умови охоплення області їхньої реальної зміни. Рівні змінних становили, %:
 
  
<center>[[Файл:ТаблП1.jpg]]</center>
+
'''Розв'язання'''. Як функції відклику прийнято <math>y_1</math>, % — середня зольність крупи пшениці після перших трьох систем для дертя (швидкість обертання рифлених вальців усіх систем 6 м/с); <math>y_2</math>, % — сумарний вихід всіх крупок, які добуваються в процесі крупоутворення; <math>y_3</math>, кДж/(кг • %) — витрата енергії на одержання 1 % продукту з 1 кг зерна.
 +
Незалежними змінними є, %: <math>x_1</math> — вихід крупи на першій системі для дертя; <math>x_2</math> — те ж на другій системі; <math>x_3</math> — те ж, для трьох систем для дертя. Інтервал варіювання для всіх <math>x_i</math>, вибрано з умови охоплення області їхньої реальної зміни. Рівні змінних становили, %:
 +
 
 +
<table width="90%" border="1">
 +
  <tr>
 +
    <th scope="col">Незалежні змінні&nbsp;</th>
 +
    <th scope="col">Нижній&nbsp;</th>
 +
    <th scope="col">Основний&nbsp;</th>
 +
    <th scope="col">Верхній&nbsp;</th>
 +
</tr>
 +
  <tr>
 +
    <th scope="row"><math>X_1</math>&nbsp;</th>
 +
    <td><center>5</center>&nbsp;</td>
 +
    <td><center>10</center>&nbsp;</td>
 +
    <td><center>15</center>&nbsp;</td>
 +
  </tr>
 +
  <tr>
 +
<th scope="row"><math>X_2</math>&nbsp;</th>
 +
    <td><center>30</center>&nbsp;</td>
 +
    <td><center>40</center>&nbsp;</td>
 +
    <td><center>50</center>&nbsp;</td>
 +
</tr>
 +
<tr>
 +
  <th scope="row"><math>X_3</math>&nbsp;</th>
 +
    <td><center>65</center>&nbsp;</td>
 +
    <td><center>70</center>&nbsp;</td>
 +
    <td><center>75</center>&nbsp;</td>
 +
</tr>
 +
  </table>
 +
 
  
 
У зв'язку з тим, що режими крупоутворення вивчалися досить детально, стало можливим ставити експерименти в області факторного простору, для якої значення всіх у близькі до оптимальних, а для опису цієї області застосувати відразу планування другого порядку. Було реалізовано центральний композиційний рототабельний план, який включає ПФЕ <math>2^3</math>, шість зіркових та шість центральних точок. Послідовність проведення дослідів була рандомізована, кожен дослід проводився тричі. У табл. 4 наведено матрицю планування та середні значення функцій відклику для кожного її рядка.
 
У зв'язку з тим, що режими крупоутворення вивчалися досить детально, стало можливим ставити експерименти в області факторного простору, для якої значення всіх у близькі до оптимальних, а для опису цієї області застосувати відразу планування другого порядку. Було реалізовано центральний композиційний рототабельний план, який включає ПФЕ <math>2^3</math>, шість зіркових та шість центральних точок. Послідовність проведення дослідів була рандомізована, кожен дослід проводився тричі. У табл. 4 наведено матрицю планування та середні значення функцій відклику для кожного її рядка.
 
За вищенаведеними формулами розраховані такі коефіцієнти в рівняннях регресії для всіх функцій відклику:
 
За вищенаведеними формулами розраховані такі коефіцієнти в рівняннях регресії для всіх функцій відклику:
  
<center>[[Файл:formula11.jpg]]</center>
+
<center>
 +
<math>\begin{align}
 +
  & {{y}_{1}}=0,65+0,0084{{z}_{1}}+0,0048{{z}_{2}}+0,0630{{z}_{3}}+0,0150{{z}_{1}}{{z}_{2}}-0,0050{{z}_{1}}{{z}_{3}}- \\
 +
& -0,0400{{z}_{2}}{{z}_{3}}+0,0038z_{1}^{2}+0,0076z_{2}^{2}+0,0314z_{3}^{2}; \\
 +
& {{y}_{2}}=43,5+1,37{{z}_{1}}+0,34{{z}_{2}}+0,89{{z}_{3}}-1,41{{z}_{1}}{{z}_{2}}-0,61{{z}_{1}}{{z}_{3}}+ \\
 +
& +0,74{{z}_{2}}{{z}_{3}}-0,83z_{1}^{2}-1,71z_{2}^{2}-1,52z_{3}^{2}; \\
 +
& {{y}_{3}}=6,4-0,28{{z}_{1}}-0,11{{z}_{2}}+0,61{{z}_{3}}+0,03{{z}_{1}}{{z}_{2}}-0,03{{z}_{1}}{{z}_{3}}- \\
 +
& -0,05{{z}_{2}}{{z}_{3}}+0,33z_{1}^{2}+0,68z_{2}^{2}+0,69z_{3}^{2}. \\
 +
\end{align}</math>
 +
</center>
  
 
Оцінки дисперсій для коефіцієнтів у цих рівняннях наведено в табл. 5.
 
Оцінки дисперсій для коефіцієнтів у цих рівняннях наведено в табл. 5.
 
Коефіцієнти при <math>z^2</math>  на порядок перевищують помилку в їхньому визначенні для всіх функцій відклику, отже, лінійними рівняннями описати їх не можна. Адекватність утворених нелінійних рівнянь було перевірено за F-критерієм.
 
Коефіцієнти при <math>z^2</math>  на порядок перевищують помилку в їхньому визначенні для всіх функцій відклику, отже, лінійними рівняннями описати їх не можна. Адекватність утворених нелінійних рівнянь було перевірено за F-критерієм.
  
<center>Таблиця 4 – Реалізація матриці ЦКРП другого порядку</center>
 
<center>[[Файл:Табл.jpg]]</center>
 
  
<center>Таблиця 5 – Оцінка дисперсій коефіцієнтів рівняння регресії за ЦКРП</center>
+
<table width="90%" border="1">
<center>[[Файл:Табл2.jpg]]</center>
+
  <caption>
 +
    Таблиця 4 – Реалізація матриці ЦКРП другого порядку
 +
  </caption>
 +
  <tr>
 +
    <th scope="col">&nbsp;</th>
 +
    <th scope="col"><math>z_0</math>&nbsp;</th>
 +
    <th scope="col"><math>z_1</math>&nbsp;</th>
 +
    <th scope="col"><math>z_2</math>&nbsp;</th>
 +
    <th scope="col"><math>z_3</math>&nbsp;</th>
 +
    <th scope="col"><math>z_1^2</math>&nbsp;</th>
 +
    <th scope="col"><math>z_2^2</math>&nbsp;</th>
 +
    <th scope="col"><math>z_3^2</math>&nbsp;</th>
 +
    <th scope="col"><math>z_1*z_2</math>&nbsp;</th>
 +
    <th scope="col"><math>z_1*z_3</math>&nbsp;</th>
 +
    <th scope="col"><math>z_2*z_3</math>&nbsp;</th>
 +
    <th scope="col"><math>y_1c</math>&nbsp;</th>
 +
    <th scope="col"><math>y_2c</math>&nbsp;</th>
 +
    <th scope="col"><math>y_3c</math>&nbsp;</th>
 +
  </tr>
 +
  <tr>
 +
    <td><center>1</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>0,75</center>&nbsp;</td>
 +
    <td><center>40,5</center>&nbsp;</td>
 +
    <td><center>8,4</center>&nbsp;</td>
 +
  </tr>
 +
<tr>
 +
    <td><center>2</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>-</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>-</center>&nbsp;</td>
 +
    <td><center>-</center>&nbsp;</td>
 +
    <td><center>0,68</center>&nbsp;</td>
 +
    <td><center>36,7</center>&nbsp;</td>
 +
    <td><center>7,3</center>&nbsp;</td>
 +
  </tr>
 +
<tr>
 +
    <td><center>3</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>-</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>-</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>-</center>&nbsp;</td>
 +
    <td><center>0,78</center>&nbsp;</td>
 +
    <td><center>41,3</center>&nbsp;</td>
 +
    <td><center>8,7</center>&nbsp;</td>
 +
  </tr>
 +
<tr>
 +
    <td><center>4</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>-</center>&nbsp;</td>
 +
    <td><center>-</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>-</center>&nbsp;</td>
 +
    <td><center>-</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>0,61</center>&nbsp;</td>
 +
    <td><center>42,7</center>&nbsp;</td>
 +
    <td><center>7,3</center>&nbsp;</td>
 +
  </tr>
 +
<tr>
 +
    <td><center>5</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>-</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>-</center>&nbsp;</td>
 +
    <td><center>-</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>0,72</center>&nbsp;</td>
 +
    <td><center>41,0</center>&nbsp;</td>
 +
    <td><center>8,7</center>&nbsp;</td>
 +
  </tr>
 +
<tr>
 +
    <td><center>6</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>-</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>-</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>-</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>-</center>&nbsp;</td>
 +
    <td><center>0,61</center>&nbsp;</td>
 +
    <td><center>37,0</center>&nbsp;</td>
 +
    <td><center>7,9</center>&nbsp;</td>
 +
  </tr>
 +
<tr>
 +
    <td><center>7</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>-</center>&nbsp;</td>
 +
    <td><center>-</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>-</center>&nbsp;</td>
 +
    <td><center>-</center>&nbsp;</td>
 +
    <td><center>0,78</center>&nbsp;</td>
 +
    <td><center>38,2</center>&nbsp;</td>
 +
    <td><center>9,2</center>&nbsp;</td>
 +
  </tr>
 +
<tr>
 +
    <td><center>8</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>-</center>&nbsp;</td>
 +
    <td><center>-</center>&nbsp;</td>
 +
    <td><center>-</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>0,62</center>&nbsp;</td>
 +
    <td><center>34,9</center>&nbsp;</td>
 +
    <td><center>8,0</center>&nbsp;</td>
 +
  </tr>
 +
<tr>
 +
    <td><center>9</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>+1,68</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>+2,83</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0,67</center>&nbsp;</td>
 +
    <td><center>44,7</center>&nbsp;</td>
 +
    <td><center>6,8</center>&nbsp;</td>
 +
  </tr>
 +
<tr>
 +
    <td><center>10</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>-1,68</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>+2,83</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0,66</center>&nbsp;</td>
 +
    <td><center>39,4</center>&nbsp;</td>
 +
    <td><center>7,7</center>&nbsp;</td>
 +
  </tr>
 +
<tr>
 +
    <td><center>11</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>+1,68</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>+2,83</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0,86</center>&nbsp;</td>
 +
    <td><center>40,7</center>&nbsp;</td>
 +
    <td><center>9,3</center>&nbsp;</td>
 +
  </tr>
 +
 
 +
<tr>
 +
    <td><center>12</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>-1,68</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>+2,83</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0,65</center>&nbsp;</td>
 +
    <td><center>37,8</center>&nbsp;</td>
 +
    <td><center>8,5</center>&nbsp;</td>
 +
  </tr>
 +
<tr>
 +
    <td><center>13</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>+1,68</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>+2,83</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0,86</center>&nbsp;</td>
 +
    <td><center>40,7</center>&nbsp;</td>
 +
    <td><center>9,3</center>&nbsp;</td>
 +
  </tr>
 +
<tr>
 +
    <td><center>14</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>-1,68</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>+2,83</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0,63</center>&nbsp;</td>
 +
    <td><center>39,3</center>&nbsp;</td>
 +
    <td><center>7,0</center>&nbsp;</td>
 +
  </tr>
 +
<tr>
 +
    <td><center>15</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0,65</center>&nbsp;</td>
 +
    <td><center>41,6</center>&nbsp;</td>
 +
    <td><center>6,4</center>&nbsp;</td>
 +
  </tr>
 +
<tr>
 +
    <td><center>16</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0,63</center>&nbsp;</td>
 +
    <td><center>42,7</center>&nbsp;</td>
 +
    <td><center>6,6</center>&nbsp;</td>
 +
  </tr>
 +
<tr>
 +
    <td><center>17</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0,66</center>&nbsp;</td>
 +
    <td><center>44,5</center>&nbsp;</td>
 +
    <td><center>6,2</center>&nbsp;</td>
 +
  </tr>
 +
<tr>
 +
    <td><center>18</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0,66</center>&nbsp;</td>
 +
    <td><center>42,9</center>&nbsp;</td>
 +
    <td><center>6,1</center>&nbsp;</td>
 +
  </tr>
 +
<tr>
 +
    <td><center>19</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0,65</center>&nbsp;</td>
 +
    <td><center>44,5</center>&nbsp;</td>
 +
    <td><center>6,8</center>&nbsp;</td>
 +
  </tr>
 +
<tr>
 +
    <td><center>20</center>&nbsp;</td>
 +
    <td><center>+</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0</center>&nbsp;</td>
 +
    <td><center>0,65</center>&nbsp;</td>
 +
    <td><center>44,0</center>&nbsp;</td>
 +
    <td><center>6,5</center>&nbsp;</td>
 +
  </tr>
 +
</table>
 +
 
 +
 
 +
<table width="90%" border="1">
 +
  <caption>
 +
    Таблиця 5 – Оцінка дисперсій коефіцієнтів рівняння регресії за ЦКРП
 +
  </caption>
 +
  <tr>
 +
    <th scope="col"><math>{{\operatorname{y}}_{i}}</math>&nbsp;</th>
 +
    <th scope="col"><math>{{\operatorname{S}}_{b0}}</math>&nbsp;</th>
 +
    <th scope="col"><math>{{\operatorname{S}}_{b1}}</math>&nbsp;</th>
 +
    <th scope="col"><math>{{\operatorname{S}}_{b2}}</math>&nbsp;</th>
 +
    <th scope="col"><math>{{\operatorname{S}}_{b3}}</math>&nbsp;</th>
 +
  </tr>
 +
  <tr>
 +
    <td><math>y_1</math>&nbsp;</td>
 +
    <td>0,0053&nbsp;</td>
 +
    <td>0,0035&nbsp;</td>
 +
    <td>0,0034&nbsp;</td>
 +
    <td>0,0046&nbsp;</td>
 +
  </tr>
 +
<tr>
 +
    <td><math>y_2</math>&nbsp;</td>
 +
    <td>0,48&nbsp;</td>
 +
    <td>0,31&nbsp;</td>
 +
    <td>0,30&nbsp;</td>
 +
    <td>0,41&nbsp;</td>
 +
  </tr>
 +
<tr>
 +
    <td><math>y_3</math>&nbsp;</td>
 +
    <td>0,13&nbsp;</td>
 +
    <td>0,8&nbsp;</td>
 +
    <td>0,8&nbsp;</td>
 +
    <td>0,11&nbsp;</td>
 +
  </tr>
 +
 
 +
  </table>
  
 
= Перелік використаних джерел =
 
= Перелік використаних джерел =

Поточна версія на 10:03, 20 березня 2012

Невідредагована стаття
Цю статтю потрібно відредагувати.
Щоб вона відповідала ВИМОГАМ.


{{{img}}}
Імя Ігор
Прізвище Пельц
По-батькові
Факультет ФІС
Група СНм-51
Залікова книжка


http://elartu.tstu.edu.ua/handle/123456789/415 Презентація доповіді (університетський репозиторій).

Планування другого порядку

Планування другого порядку застосовується для математичного опису об'єкта поблизу екстремальної точки статистичної характеристики або тоді, коли необхідний точніший опис в інших точках факторного простору. При цьому використовують поліном другого порядку:

[math]\operatorname{y}={{a}_{0}}+\sum\limits_{i=1}^{n}{{{a}_{i}}{{x}_{i}}+}\sum\limits_{i=1}^{n}{{{a}_{ij}}{{x}_{i}}{{x}_{j}}+}...+\sum\limits_{i=1}^{n}{{{a}_{i,i}}x_{i}^{2}}[/math]

Задача, як і в ПФЕ, полягає у визначенні методом найменших квадратів за результатами спланованого експерименту коефіцієнтів цього рівня за умови, що виконуються передумови регресійного аналізу. ПФЕ типу [math]2^n[/math] дає змогу дістати роздільні оцінки як лінійних коефіцієнтів bi (після переходу до безрозмірних z), так і коефіцієнтів парних взаємодій bij. Точки ПФЕ лежать у вершині n-вимірного куба. Вектор-стовпці лінійних факторів матриці планування ортогональні між собою, тобто виконується умова

[math]\sum\limits_{g=1}^{N}{{{z}_{gi}}{{z}_{gj}}=0;i\ne j.}[/math]

З теорії інтерполяції (апроксимації) відомо, що для розв'язання задачі знаходження роздільних оцінок параметрів апроксимуючого виразу число рівнів для кожної із змінних повинно бути на одиницю більше ступеня апроксимуючого полінома, тобто для полінома другого порядку число рівнів дорівнює трьом. Однак, як показали дослідження, ПФЕ типу [math]3^n[/math] (планування на трьох рівнях) не є раціональним через велике число дослідів. Задача розв'язується іншим способом. До ПФЕ типу [math]2^n[/math] додають центральну точку з координатами (0, 0, ..., 0) і зіркові точки з координатами (0, 0, ..., ±α), які лежать на сфері діаметра 2α (рис. 1). Зіркові точки будують на осях факторного простору. Вибір відстані від нульової точки до зіркової, яка визначається плечем α, залежить від критерію оптимальності плану.


ОЦКП.jpg
Рисунок 1 - Ортогональне центральне композиційне планування

Розглянемо планування, оптимальне з точки зору незалежності оцінок bi,i Його називають ортогональним центральним композиційним плануванням (ОЦКП), тобто планом, в якому критерієм оптимальності є ортогональність стовпців матриці планування. Композиційним таке планування, як й інші форми планування другого порядку, називається тому, що новий план дістають шляхом компонування первинного двофакторного плану з деякою кількістю додаткових точок. Оскільки в числі цих додаткових точок обов'язково фігурує центральна, в якій всі змінні хi мають середній рівень, а zi=0, плани називають центральними. Ортогональність матриці композиційного планування забезпечується виконанням рівностей

[math]\sum\limits_{g=1}^{N}{z_{g,o}^{2}z_{y,i}^{2}=0,\sum\limits_{g=1}^{N}{z_{gi}^{2}z_{gj}^{2}=0,i\ne j}.}[/math]

де і-номер фактора; j-номер рядка; g-номер досліду. Для ортогоналізації першого з цих співвідношень застосовується перехід до нової змінної

[math]\overset{-}{z}\,_{i}^{2}=z_{i}^{2}-\frac{\sum{z_{gi}^{2}}}{N}=z_{i}^{2}-\overset{-}{z}\,_{i}^{2}[/math]

де N-загальне число експериментів. Величина zi залежить тільки від числа факторів n і числа дослідів N, яке звичайно вибирається так, що

[math]\operatorname{N}={{N}_{n}}+{{N}_{a}}+{{N}_{0}}={{2}^{n}}+2n+1,[/math]

де Nn = [math]2^n[/math]-кількість вершин гіперкуба при ПФЕ; Nα = 2n-число зіркових точок; N0= 1- число дослідів у центрі плану. Якщо N вибрати так, то zi визначають за формулою

[math]\operatorname{z}_{i}^{2}=\frac{{{2}^{n}}+2{{\alpha }^{2}}}{N}.[/math]

Тому ортогоналізація другої з вищенаведених умов досягається вибором бажаного α. Для зручності підготовки і планування величини α, N, Nn, Nα обчислені і табульовані залежно від числа факторів (табл. 1).


Таблиця1 - Підготовка ОЦКП другого порядку.
[math]{{\operatorname{N}}_{0}}[/math]  [math]{{\operatorname{N}}_{n}}[/math]  [math]{{\operatorname{N}}_{\alpha }}[/math]  [math]\alpha[/math] 
1,000 
15  1,215 
16  25  1,414 
32  10  43  1,596 
64  12  77  1,706 
128  14  143  1,909 

Сформуємо матриці ОЦКП для двох факторів. При n = 2 отримаємо α = 1,0 (див. табл. 1), для обчислення скористаємося наведеними формулами. Оскільки [math]z1^2=(4+2)/9=0,667[/math], то стовпець добувається відніманням одного і того ж числа 0,67 від числа стовпця [math]zi^2[/math] того ж рядка (табл. 2).

Таблиця2 - Матриця ОЦКП другого порядку для двофакторного експерименту
Дослід  [math]z_0[/math]  [math]z_1[/math]  [math]z_2[/math]  [math]z_1^2[/math]  [math]z_2^2[/math]  [math]z_{1c}^2[/math]  [math]z_{2c}^2[/math]  [math]z_1*z_2[/math] 
Вершини квадрата 
1
 
+
 
+
 
+
 
+
 
+
 
0,33
 
0,33
 
+
 
 
2
 
+
 
-
 
+
 
+
 
+
 
0,33
 
0,33
 
-
 
 
3
 
+
 
+
 
-
 
+
 
+
 
0,33
 
0,33
 
-
 
 
4
 
+
 
-
 
-
 
+
 
+
 
0,33
 
0,33
 
+
 
Зіркові точки 
5
 
+
 
-
 
0
 
+
 
0
 
0,33
 
-0,67
 
0
 
 
6
 
+
 
+
 
0
 
+
 
0
 
0,33
 
-0,67
 
0
 
 
7
 
+
 
0
 
-
 
0
 
+
 
-0,67
 
0,33
 
0
 
 
8
 
+
 
0
 
+
 
0
 
+
 
-0,67
 
0,33
 
0
 
Центр 
9
 
+
 
0
 
0
 
0
 
0
 
-0,67
 
-0,67
 
0
 

Реалізація експериментів за ОЦКП здійснюється за тією ж методикою, що і ПФЕ. Таким чином, через випадковий характер зміни вихідної величини у у кожній точці хg проводиться m паралельних дослідів і обчислюється середнє значення функції відклику

[math]{{\overset{-}{y}}_{g}}=\frac{\sum\limits_{d=1}^{m}{{{y}_{g}}d}}{m}[/math]

Перед реалізацією проводиться рандомізація рядків матриці планування. Перевірка відтворюваності, як і в ПФЕ, виконується за критерієм Кохрена, після чого обчислюється оцінка дисперсії відтворюваності


[math]S_{y}^{2}=S_{vidtv}^{2}=\frac{\sum\limits_{g=1}^{N}{S_{g}^{2}}}{N}.[/math]


Методика утворення математичної моделі незначно відрізняється від методики опису результатів ПФЕ. Коефіцієнти регресії при ОЦКП обчислюються за формулою

[math]{{b}_{i}}=\frac{\sum\limits_{g=1}^{N}{{{z}_{gi}}{{y}_{g}}}}{\sum\limits_{g=1}^{N}{z_{gi}^{2}}}.[/math]

У цьому плануванні оцінки дисперсій коефіцієнтів bi (точність їхнього обчислення) не однакові, оскільки не однаковий знаменник у формулі дисперсії

[math]S_{bi}^{2}=\frac{S_{y}^{2}}{m\sum\limits_{g=1}^{N}{z_{gi}^{2}}}.[/math]

В ПФЕ знаменник однаковий і дорівнює mN. Це істотний недолік ОЦКП, і тому часто надають перевагу складнішому за обчислювальними процедурами рототабельному плануванню.

Рототабельне планування

У зв'язку з тим, що дисперсії коефіцієнтів рівняння регресії при ОЦКП нерівномірні, ортогональність матриці часто не є досить сильним критерієм оптимальності планування другого порядку. Його заміняють критерієм ротоптабельності, тобто однаковості дисперсій коефіцієнтів при повороті координатних осей на будь-який кут. Зазначимо, що при плануванні першого порядку ортогональність матриці просто збігається з її рототабельністю, тому ПФЕ доцільно називати рототабельним. Щоб зробити план другого порядку рототабельним, вибирають для сфери, на якій розташовуються зіркові точки, радіус (зіркове плече) за формулою

[math]\alpha ={{2}^{n/2}}.[/math]

Інша умова рототабельності — збільшення числа дослідів на поверхні нульової сфери, тобто в центрі плану. У зв'язку з цим виникає повна назва методу: центральне композиційне рототабельне планування (ЦКРП). Таким чином ЦКРП багато в чому нагадує ортогональне планування, проте метод рототабельного планування експерименту дає змогу дістати точніший математичний опис поверхні відклику порівняно з ОЦКП, завдяки збільшенню числа дослідів у центрі плану і спеціальному вибору величини зіркового плеча α. Як і для ОЦКП, основні характеристики матриць рототабельного планування табульовані (табл. 3). Позначення тут ті самі, що і для ОЦКП (див. табл. 2). При ЦКРП, починаючи з n = 5, можна застосувати ДФЕ(дробовий факторний експеримент).

Таблиця 3 – Підготовка ЦКРП другого порядку
[math]{{\operatorname{N}}_{n}}[/math]  [math]{{\operatorname{N}}_{\alpha }}[/math]  [math]{{\operatorname{N}}_{0}}[/math]  [math]\alpha[/math] 
13  1,414 
20  1,680 
16  31  2,000 
10  32  10  52  2,378 
15  64  12  91  1,828 
21  128  14  163  1,333 

При рототабельному плануванні для обчислення коефіцієнтів моделі і відповідних оцінок дисперсій знаходять спеціальні комплекси:

[math]\begin{align} & B=\frac{nN}{(n+2)(N-{{N}_{0}})}; \\ & A=\frac{1}{2B[(n+2)B-n]}; \\ & C=\frac{N}{N-{{N}_{0}}}, \\ \end{align}[/math]

де n-число факторів; N-загальне число дослідів у плануванні; N0-число дослідів у центрі плану. За результатами експериментів обчислюють такі суми:

[math]\begin{align} & {{S}_{0}}=\sum\limits_{g=1}^{N}{{{y}_{g}}}; \\ & {{S}_{i}}=\sum\limits_{g=1}^{N}{{{y}_{g}}}{{z}_{gi}};i=1,2,...,n; \\ & {{S}_{ik}}=\sum\limits_{g=1}^{N}{{{z}_{gi}}{{z}_{gk}}{{y}_{g}}};i\ne k; \\ & {{S}_{ii}}=\sum\limits_{g=1}^{N}{z_{gi}^{2}{{y}_{g}}};i=1,2,...,n. \\ \end{align}[/math]


Коефіцієнти моделі тут розраховують за формулами

[math]\begin{align} & {{b}_{0}}=\frac{2AB}{N}[{{S}_{0}}B(n+2)-C\sum{{{S}_{ii}}}]; \\ & {{b}_{i}}=\frac{C{{S}_{i}}}{N}; \\ & {{b}_{ik}}=\frac{{{C}^{2}}{{S}_{ik}}}{BN},i\ne k; \\ & {{b}_{ii}}=\frac{AC}{N}\{{{S}_{ii}}[B(n+2)-n]+C(1-B)\sum\limits_{i=1}^{n}{{{S}_{ii}}-2B{{S}_{0}}}\}. \\ & \\ \end{align}[/math]

Оцінки дисперсій для обчислених коефіцієнтів знаходять за такими формулами:

[math]\begin{align} & S_{b0}^{2}=\frac{2AB(n+2)}{N}S_{y}^{2}; \\ & S_{bi}^{2}=\frac{S_{y}^{2}}{N-{{N}_{0}}};i=1,2,...,n; \\ & S_{bik}^{2}=\frac{{{C}^{2}}S_{y}^{2}}{N},i\ne k; \\ & {{S}_{bii}}=\frac{A{{C}^{2}}S_{y}^{2}}{N}[B(n+1)-(n-1)]. \\ \end{align}[/math]

У цих формулах дисперсія відтворюваності [math]S_{y}^{2}[/math] визначається за результатами дослідів у нульовій точці

[math]\begin{align} & S_{y}^{2}=\frac{1}{{{N}_{0}}-1}\sum\limits_{g=1}^{{{N}_{0}}}{{{({{y}_{ge}}-\overset{-}{y}\,)}^{2}}}; \\ & \overset{-}{y}\,=\frac{1}{{{N}_{0}}}\sum\limits_{g=1}^{{{N}_{0}}}{{{y}_{ge.}}} \\ \end{align}[/math]

Дисперсія адекватності оцінюється за формулою

[math]S_{adekv}^{2}=\frac{\sum\limits_{g=1}^{N}{{{({{y}_{ge}}-{{y}_{grozr}})}^{2}}-S_{y}^{2}({{N}_{0}}-1)}}{N-\frac{(n+2)(n+1)}{2}(N-1)},[/math]

якшо число ступенів вільності

[math]{{f}_{adekv}}={{N}_{0}}-\frac{(n+2)(n+1)}{2}-({{N}_{0}}-1).[/math]

Приклад

Скласти матрицю ЦКРП на прикладі побудови математичної моделі технологічного процесу крупоутворення (див.: Пищевая технология.— 1976.— № 4.— С. 121—124).

Розв'язання. Як функції відклику прийнято [math]y_1[/math], % — середня зольність крупи пшениці після перших трьох систем для дертя (швидкість обертання рифлених вальців усіх систем 6 м/с); [math]y_2[/math], % — сумарний вихід всіх крупок, які добуваються в процесі крупоутворення; [math]y_3[/math], кДж/(кг • %) — витрата енергії на одержання 1 % продукту з 1 кг зерна. Незалежними змінними є, %: [math]x_1[/math] — вихід крупи на першій системі для дертя; [math]x_2[/math] — те ж на другій системі; [math]x_3[/math] — те ж, для трьох систем для дертя. Інтервал варіювання для всіх [math]x_i[/math], вибрано з умови охоплення області їхньої реальної зміни. Рівні змінних становили, %:

Незалежні змінні  Нижній  Основний  Верхній 
[math]X_1[/math] 
5
 
10
 
15
 
[math]X_2[/math] 
30
 
40
 
50
 
[math]X_3[/math] 
65
 
70
 
75
 


У зв'язку з тим, що режими крупоутворення вивчалися досить детально, стало можливим ставити експерименти в області факторного простору, для якої значення всіх у близькі до оптимальних, а для опису цієї області застосувати відразу планування другого порядку. Було реалізовано центральний композиційний рототабельний план, який включає ПФЕ [math]2^3[/math], шість зіркових та шість центральних точок. Послідовність проведення дослідів була рандомізована, кожен дослід проводився тричі. У табл. 4 наведено матрицю планування та середні значення функцій відклику для кожного її рядка. За вищенаведеними формулами розраховані такі коефіцієнти в рівняннях регресії для всіх функцій відклику:

[math]\begin{align} & {{y}_{1}}=0,65+0,0084{{z}_{1}}+0,0048{{z}_{2}}+0,0630{{z}_{3}}+0,0150{{z}_{1}}{{z}_{2}}-0,0050{{z}_{1}}{{z}_{3}}- \\ & -0,0400{{z}_{2}}{{z}_{3}}+0,0038z_{1}^{2}+0,0076z_{2}^{2}+0,0314z_{3}^{2}; \\ & {{y}_{2}}=43,5+1,37{{z}_{1}}+0,34{{z}_{2}}+0,89{{z}_{3}}-1,41{{z}_{1}}{{z}_{2}}-0,61{{z}_{1}}{{z}_{3}}+ \\ & +0,74{{z}_{2}}{{z}_{3}}-0,83z_{1}^{2}-1,71z_{2}^{2}-1,52z_{3}^{2}; \\ & {{y}_{3}}=6,4-0,28{{z}_{1}}-0,11{{z}_{2}}+0,61{{z}_{3}}+0,03{{z}_{1}}{{z}_{2}}-0,03{{z}_{1}}{{z}_{3}}- \\ & -0,05{{z}_{2}}{{z}_{3}}+0,33z_{1}^{2}+0,68z_{2}^{2}+0,69z_{3}^{2}. \\ \end{align}[/math]

Оцінки дисперсій для коефіцієнтів у цих рівняннях наведено в табл. 5. Коефіцієнти при [math]z^2[/math] на порядок перевищують помилку в їхньому визначенні для всіх функцій відклику, отже, лінійними рівняннями описати їх не можна. Адекватність утворених нелінійних рівнянь було перевірено за F-критерієм.


Таблиця 4 – Реалізація матриці ЦКРП другого порядку
  [math]z_0[/math]  [math]z_1[/math]  [math]z_2[/math]  [math]z_3[/math]  [math]z_1^2[/math]  [math]z_2^2[/math]  [math]z_3^2[/math]  [math]z_1*z_2[/math]  [math]z_1*z_3[/math]  [math]z_2*z_3[/math]  [math]y_1c[/math]  [math]y_2c[/math]  [math]y_3c[/math] 
1
 
+
 
+
 
+
 
+
 
+
 
+
 
+
 
+
 
+
 
+
 
0,75
 
40,5
 
8,4
 
2
 
+
 
+
 
+
 
-
 
+
 
+
 
+
 
+
 
-
 
-
 
0,68
 
36,7
 
7,3
 
3
 
+
 
+
 
-
 
+
 
+
 
+
 
+
 
-
 
+
 
-
 
0,78
 
41,3
 
8,7
 
4
 
+
 
+
 
-
 
-
 
+
 
+
 
+
 
-
 
-
 
+
 
0,61
 
42,7
 
7,3
 
5
 
+
 
-
 
+
 
+
 
+
 
+
 
+
 
-
 
-
 
+
 
0,72
 
41,0
 
8,7
 
6
 
+
 
-
 
+
 
-
 
+
 
+
 
+
 
-
 
+
 
-
 
0,61
 
37,0
 
7,9
 
7
 
+
 
-
 
-
 
+
 
+
 
+
 
+
 
+
 
-
 
-
 
0,78
 
38,2
 
9,2
 
8
 
+
 
-
 
-
 
-
 
+
 
+
 
+
 
+
 
+
 
+
 
0,62
 
34,9
 
8,0
 
9
 
+
 
+1,68
 
0
 
0
 
+2,83
 
0
 
0
 
0
 
0
 
0
 
0,67
 
44,7
 
6,8
 
10
 
+
 
-1,68
 
0
 
0
 
+2,83
 
0
 
0
 
0
 
0
 
0
 
0,66
 
39,4
 
7,7
 
11
 
+
 
0
 
+1,68
 
0
 
0
 
+2,83
 
0
 
0
 
0
 
0
 
0,86
 
40,7
 
9,3
 
12
 
+
 
0
 
-1,68
 
0
 
0
 
+2,83
 
0
 
0
 
0
 
0
 
0,65
 
37,8
 
8,5
 
13
 
+
 
0
 
0
 
+1,68
 
0
 
0
 
+2,83
 
0
 
0
 
0
 
0,86
 
40,7
 
9,3
 
14
 
+
 
0
 
0
 
-1,68
 
0
 
0
 
+2,83
 
0
 
0
 
0
 
0,63
 
39,3
 
7,0
 
15
 
+
 
0
 
0
 
0
 
0
 
0
 
0
 
0
 
0
 
0
 
0,65
 
41,6
 
6,4
 
16
 
+
 
0
 
0
 
0
 
0
 
0
 
0
 
0
 
0
 
0
 
0,63
 
42,7
 
6,6
 
17
 
+
 
0
 
0
 
0
 
0
 
0
 
0
 
0
 
0
 
0
 
0,66
 
44,5
 
6,2
 
18
 
+
 
0
 
0
 
0
 
0
 
0
 
0
 
0
 
0
 
0
 
0,66
 
42,9
 
6,1
 
19
 
+
 
0
 
0
 
0
 
0
 
0
 
0
 
0
 
0
 
0
 
0,65
 
44,5
 
6,8
 
20
 
+
 
0
 
0
 
0
 
0
 
0
 
0
 
0
 
0
 
0
 
0,65
 
44,0
 
6,5
 


Таблиця 5 – Оцінка дисперсій коефіцієнтів рівняння регресії за ЦКРП
[math]{{\operatorname{y}}_{i}}[/math]  [math]{{\operatorname{S}}_{b0}}[/math]  [math]{{\operatorname{S}}_{b1}}[/math]  [math]{{\operatorname{S}}_{b2}}[/math]  [math]{{\operatorname{S}}_{b3}}[/math] 
[math]y_1[/math]  0,0053  0,0035  0,0034  0,0046 
[math]y_2[/math]  0,48  0,31  0,30  0,41 
[math]y_3[/math]  0,13  0,8  0,8  0,11 

Перелік використаних джерел

  1. http://tstu.edu.ua/(березень2010)
  2. Аністратенко В.О., Федоров В.Г. Математичне планування експериментів в АПК http://tstu.edu.ua/(березень2010)