Кодирование варьируемых параметров
| Кодовые обозначения факторов |
|
|
|
| Варьируемые параметры |
| В |
|
| Единица измерения |
|
|
|
Основной уровень ( =0)
Единица варьирования
Верхний уровень ( =+1)
Нижний уровень ( =-1)
| 0,065 0,015 0,05 0,08 |
Результаты двух параллельных опытов приведены в табл. 2.
Таблица2
Результаты эксперимента по линейному плану типа 23
|
|
|
|
|
|
|
|
| 11,0 15,3 31,7 23,5 | 14,8 15,4 35,7 22,5 | 12,8 17,3 39,3 25,1 | 11,0 16,2 36,3 28,6 | 12,5 17,2 35,6 26,1 | 13,4 18,9 38,2 27,0 |
Вычисляем средние арифметические
и построчные выборочные дисперсии
параллельных опытов (табл. 3).









Таблица 3
Построчные средние арифметические и выборочные дисперсии
|
|
|
|
|
|
| 12,86 35,56 23,7 | 3,61 1,27 14,45 1,72 | 12,3 17,43 36,7 27,23 | 1,47 1,86 1,81 1,47 |
Вычисляем оценки коэффициентов регрессии (табл. 5).
Составляем расширенную матрицу планирования эксперимента (табл.4)
Таблица 4.
Расширенная матрица планирования ПФЭ 23.
| Номер опыта | х0 | x1 | x2 | x3 | x1x2 | x1x3 | x2x3 | x1x2x3 |
| 1. | +1 | -1 | -1 | -1 | +1 | +1 | +1 | -1 |
| 2. | +1 | +1 | -1 | -1 | -1 | -1 | +1 | +1 |
| 3. | +1 | -1 | +1 | -1 | -1 | +1 | -1 | +1 |
| 4. | +1 | +1 | +1 | -1 | +1 | -1 | -1 | -1 |
| 5. | +1 | -1 | -1 | +1 | +1 | -1 | -1 | +1 |
6.
| +1 | +1 | -1 | +1 | -1 | +1 | -1 | -1 |
| 7. | +1 | -1 | +1 | +1 | -1 | -1 | +1 | -1 |
| 8. | +1 | +1 | +1 | +1 | +1 | +1 | +1 | +1 |
Рассчитаем линейные коэффициенты регрессии. Любой коэффициент уравнения регрессии bj определяется скалярным произведением столбца y на соответствующий столбец xj , отнесенным к числу опытов в матрице планирования N:
Таблица 5
Оценки коэффициентов регрессии
|
|
|
|
|
|
|
|
| 22,72 | -1,63 | 8,075 | 0,69 | -3,7 | 4,38 | 0,47 |
Проверяем гипотезу об однородности построчных выборочных дисперсий:
G=
=
,
где
- табличное значение критерия Кохрена при уровне значимости
и числах степеней свободы
и N=8. (см. Приложение)
Условие =
выполняется, т.е. гипотеза о том, что расхождения между построчными выборочными дисперсиями незначимые, не противоречит экспериментальным данным, поэтому можно, усреднив
, вычислить дисперсию воспроизводимости опытов
=
=28,66/8=3,58
при числе степеней свободы
.
Проверяем значимость коэффициентов регрессии, вычисляя их доверительный интервал
,
где
=2,12 - критическое значение t- распределения при двустороннем ограничении, доверительном уровне
и числе степеней свободы
.
Коэффициент регрессии статистически значимый, если

Следовательно, незначимы и не должны включаться в уравнение регрессии коэффициенты
.
Таким образом, уравнение регрессии имеет вид:
Проверяем гипотезу о наличии зависимости между функцией отклика и факторами.
Вычисляем среднее арифметическое всех результатов эксперимента (уравнение нулевого порядка):

Вычисляем остаточную дисперсию для уравнения нулевого порядка

где
– расчетное значение функции отклика для u-го варианта.
Предварительно заполняем таблицу 6.
Таблица 6
| 1. | 2. | 3. | 4. | 5. | 6. | 7. | 8. | |
| 20,405 | 7,505 | 20,395 | 31,055 | 13,025 | 17,645 | 21,775 | 41,195 |
| 2,315 | 15,215 | 2,325 | -8,335 | 9,695 | 5,075 | 0,945 | -18,475 |
| 5,36 | 231,496 | 5,405 | 69,47 | 93,99 | 25,755 | 0,89 | 341,325 |
|
|

с числом степеней свободы
Вычисляем дисперсию адекватности для полученного уравнения
,
предварительно заполнив таблицу 7.
Таблица 7
| 1. | 2. | 3. | 4. | 5. | 6. | 7. | 8. | |
| 20,405 | 7,505 | 20,395 | 31,055 | 13,025 | 17,645 | 21,775 | 41,195 |
| -7,545 | 8,495 | 15,165 | -7,355 | -0,725 | -0,215 | 14,925 | -13,965 |
| 56,93 | 72,16 | 229,97 | 54,09 | 0,525 | 0,046 | 222,75 | 195,02 |
|
|

с числом степеней свободы 
n1 – число коэффициентов уравнения регрессии (в нашем случае n1 =6, т.к.
)
Находи отношение большей из найденных дисперсий к меньшей.

Табличное значение критерия Фишера
=19,35 при m=
=7 и n= 
Если
то полученная модель описывает поверхность отклика не лучше, чем среднее арифметическое
, т.е. не имеет информационной ценности.
Проверяем приемлемость линейного уравнения,
Линейное уравнение приемлемо, если разность
статистически незначима, т.е. выполняется неравенство
,
|
где
- средневзвешенное двух дисперсий с числом степеней свободы
=8+2-2=8
- дисперсия коэффициентов регрессии;
- дисперсия среднего значения
;
- критическое значение t- распределения при двустороннем ограничении, доверительном уровне
и числе степеней свободы
.
=2. при v=8 и g=0,95.
Ставим опыты в центре плана:
и


= 


Отсюда,

Т.к. условие 
Не выполняется, то гипотезу о приемлемости линейной модели функции отклика надо отвергнуть и строить уравнение регрессии второго порядка.
Таблица 8
Результаты дополнительных опытов по плану второго порядка
|
|
|
|
|
|
|
|
| 13,2 | 11,2 | 13,7 | 32,5 | 33,3 | 34,1 | ||
| 12,6 | 12,5 | 12,1 | 12,3 | 13,3 | 14,9 | ||
| 13,8 | 13,0 | 13,9 | 14,7 | 15,7 | 16,1 | ||
| 12,0 | 15,7 | 13,1 |
Вычисляем средние арифметические 







по формуле







построчные выборочные дисперсии параллельных опытов (табл. 9).
Проверяем гипотезу об однородности всех построчных выборочных дисперсий в табл. 3 и 8, вычисляя статистику

Таблица 9
Построчные средние арифметические и выборочные дисперсии
дополнительных опытов по композиционному плану
| u |
|
| u |
|
|
| 12,7 | 1,75 | 33,3 | 0,64 | ||
| 12,4 | 0,7 | 13,5 | 1,72 | ||
| 13,56 | 0,242 | 15,5 | 0,52 | ||
| 13,6 | 3,61 |
Критерий Кохрена при уровне значимости
=0,05, числах степеней свободы
и
=15 определяем путём интерполирования:

Поскольку условие
G=
,
выполняется, можно, усредняя дисперсии
, вычислить дисперсию воспроизводимости опытов

при числе степеней свободы 
По формуле

вычисляем оценки коэффициентов регрессии (табл. 10).
Таблица 10
Оценки коэффициентов уравнения регрессии)
|
|
|
|
|
|
| 19,65 | -0,11 | 0,53 | 0,045 | 0,018 |
|
|
|
|
|
|
| 0,013 | 0,010 | -3,68 | 0,55 | 0,47 |
=0)
Единица варьирования
Верхний уровень (
,