Тандартизация масштаба факторов
Для удобства расчетов масштаб факторов выбирают так, чтобы значение верхнего уровня было равно +1, а нижнего -1. С этой целью делают преобразование начала координат факторов и переходят к нормированному (стандартному) масштабу
, (2.2)
где x0i – базовый или начальный уровень i-й переменной, Δxi – интервал варьирования по i-й переменной.
2.2. Составление матрицы планирования ПФЭ
План ПФЭ изображают в виде таблицы, столбцы которой отражают уровни факторов, а строки — номера опытов. Эти таблицы называют матрицами планирования (МП) эксперимента. Поскольку значения уровней факторов по модулю всегда равны единице, то обычно в МП записывают только знак уровня (т. е. «+» вместо «1» и «-» вместо «-1»)., В таблицах 2.2 и 2.3 в качестве примера приведены МП для ПФЭ типа 22 и 23 соответственно.
Таблица 2.2 — МП ПФЭ типа 22 | Таблица 2.3 — МП ПФЭ типа 23 | ||||||||
N | x1 | x2 | y | N | x1 | x2 | x3 | y | |
- | - | y1 | - | - | - | y1 | |||
+ | - | y2 | + | - | - | y2 | |||
- | + | y3 | - | + | - | y3 | |||
+ | + | y4 | + | + | - | y4 | |||
- | - | + | y5 | ||||||
+ | - | + | y6 | ||||||
- | + | + | y7 | ||||||
+ | + | + | y8 |
Геометрической интерпретацией ПФЭ 22 является квадрат в факторной плоскости (рисунок 2.1), а ПФЭ 23 — куб.
Здесь нормированные координаты x] и х2 проходят через точку пересечения основных уровней факторов, и масштабы их осей выбраны так, чтобы интервал варьирования равнялся 1. Тогда условия проведения опытов в МП эксперимента будут соответствовать вершинами квадрата, центром которого является основной уровень. Если п>3, то фигуру, задающую в многомерном пространстве область эксперимента, называют гиперкубом.
Влияние факторов на выходной параметр может зависеть от уровня, на котором находится другой фактор, или от сочетания уровней нескольких факторов. Если априорно не известно, что такой зависимости между факторами нет, то строят развернутую МП, учитывающую не только факторы, но и их взаимодействия.
Рисунок 2.1 — Геометрическая интерпретация ПФЭ типа 22.
Пример развернутой МП типа для ПФЭ 23 представлен в табл. 2.4. Фиктивный фактор x0 вводят для удобства компьютерного расчета свободного члена b0 (для идентичности формул).
Таблица 2.4 — Развернутая МП ПФЭ типа 23. | |||||||||
N | x0 | x1 | x2 | x3 | x1 x2 | x1 x3 | x2 x3 | x1x2 x3 | y |
+ | - | - | - | + | + | + | - | y1 | |
+ | + | - | - | - | - | + | + | y2 | |
+ | - | + | - | - | + | - | + | y3 | |
+ | + | + | - | + | - | - | - | y4 | |
+ | - | - | + | + | - | - | + | y5 | |
+ | + | - | + | - | + | - | - | y6 | |
+ | - | + | + | - | - | + | - | y7 | |
+ | + | + | + | + | + | + | + | y8 |
Для МП эксперимента должны выполняться следующие условия:
а) симметричности относительно центра эксперимента , где i — номер фактора, j — номер опыта, N — число опытов;
б) нормировки ;
в) ортогональности ,если f≠i.
г) рототабельности — условия равноточного предсказания исследуемого параметра на равных расстояниях от центра эксперимента вне зависимости от направления.
Условие ортогональности позволяет упростить вычисления и получить независимые оценки коэффициентов регрессии. Это означает, в частности, что замена нулем любого коэффициента в уравнении ММ не изменит оценок остальных коэффициентов. Это свойство может быть полезным, когда точный вид модели не известен и требуется по экспериментальным данным отобрать факторы, существенно влияющие на исследуемый параметр. Если условие ортогональности не выполняется, после исключения каждого незначимого коэффициента необходимо пересчитывать оценки оставшихся коэффициентов и их дисперсии. При этом могут измениться как доверительные интервалы, так и выводы относительно коэффициентов значимости;
Матрица, удовлетворяющая условиям симметричности, нормировки и ортогональности, называется оптимальной.
МП ПФЭ является оптимальной для линейных ММ. Если же ММ содержит взаимодействия, то условие рототабельности не выполняется.
2.3. Порядок постановки ПФЭ
Для оценки точности эксперимента для каждой i- точки факторного пространства (для каждого сочетания уровней факторов МП) проводят n опытов. В результате получают значения уi1, yi2, ..., уin исследуемого параметра, для которых находят среднее значение
. (2.3)
При этом опыты в одной точке проводят не подряд, а обходят все точки в первой серии опытов, затем во второй, и так далее до n-й. Для уменьшения влияния внешней среды и неконтролируемых факторов внутри каждой серии точки факторного пространства обходят случайным образом – рандомизируют последовательность опытов. Рандомизацию опытов можно провести с помощью генератора случайных чисел.
2.4. Проверка воспроизводимости опытов (однородности дисперсий)
Опыт считается воспроизводимым, если дисперсиявыходного параметра у, однородна в каждой точке факторного пространства. Оценка s2дисперсии определяется для каждой точки факторного пространства по формуле:
. (2.4)
Если сравниваемое количество дисперсий больше двух, можно воспользоваться критерием Кохрена (Приложение В). Этот критерий пригоден для случаев, когда во всех точках имеется одинаковое число повторных опытов. Для этого формулируют нулевую гипотезу Н0: s2y1= s2y2=…= s2yn – дисперсии значений отклика в каждом опыте однородны.
Для этого подсчитывается дисперсия значений отклика для каждого из опытов эксперимента, а затем из всех дисперсий находится наибольшая , которая делится на сумму всех дисперсий:
, (2.5)
Если полученное значение дисперсионного отношения Gрасч>Gтабл(α;γ1;γ2) больше приведенного в таблице, это означает, что дисперсии значимо отличаются друг от друга, т. е. что они неоднородны. Здесь α — уровень значимости, γ1 и γ2 — степени свободы, при этом γ1=m-1, γ2=N.