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

Blue check.png Дана стаття являється неперевіреним навчальним завданням.
Студент: ihor_p
Викладач: Назаревич О.Б.
Термін до: 05 березня 2010

До вказаного терміну стаття не повинна редагуватися іншими учасниками проекту. Після завершення терміну виконання будь-який учасник може вільно редагувати дану статтю і витерти дане попередження, що вводиться за допомогою шаблону.



 .................... Презентація доповіді (університетський репозиторій).

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

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

Formula.jpg

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

Formula2.gif

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

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

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

Formula3.gif

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

Formula4.gif

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

Formula5.gif

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

Formula6.gif

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

Таблиця1 - Підготовка ОЦКП другого порядку
Табл1.gif

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

Таблиця2 - Матриця ОЦКП другого порядку для двофакторного експерименту
Табл2.gif

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

Formula7.gif

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

Formula8.gif

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

Formula9.gif

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

Formula10.gif

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

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

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

Formula4.jpg

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

Таблиця 3 – Підготовка ЦКРП другого порядку
ТаблПід.jpg

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

Formula5.jpg

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

Formula6.jpg
Formula7.jpg

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

Formula8.jpg

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

Formula9.jpg