Для изучения влияния некоторых технологических факторов на исследуемый процесс были поставлены эксперименты по схеме полного факторного эксперимента для двух входных факторов на двух уровнях.
Требуется построить уравнение регрессии, учитывая все взаимодействия факторов, проверить значимость коэффициентов для уравнения регрессии (со статистической надежностью 95%), проверить полученную модель на адекватность и произвести ее интерпретацию.
Таблица 1. Значения выходных показателей y.
Наименование опытов
Выходной показатель y
Значения выходного показателя y
Полный факторный эксперимент для двух входных факторов на двух уровнях типа 22
y1 опыт
y2 опыт
y3 опыт
y4 опыт
Наименование опытов
Выходной показатель y
Значения выходного показателя y
Дополнительные параллельные опыты из центра плана
y1
y2
y3
Решение
Составление исходной матрицы планирования по схеме полного факторного эксперимента для двух входных факторов на двух уровнях.
Чтобы рассчитать математическую модель, составим план проведения эксперимента.
По условию задания рассмотрим полный факторный эксперимент для двух входных факторов на двух уровнях.
Число опытов в полном факторном эксперименте определяется по формуле:
,
где N – число опытов;
k – число изучаемых входных факторов.
Значит, число опытов для двух входных факторов будет равно:
.
Таким образом, необходимо провести 4 опыта.
Составим матрицу планирования полного факторного эксперимента для двух входных факторов (на двух уровнях).
Таблица 2 - План эксперимента
№
опыта
Условное обозначение уровней изменения входных факторов
-1
-1
+1
-1
-1
+1
+1
+1
Уравнение регрессии при проведении полного факторного эксперимента для двух входных факторов имеет вид:
,
где - среднее арифметическое значение выходного показателя во всех опытах; - характеризует влияние входного фактора на выходной показатель ; - характеризует влияние входного фактора на выходной показатель ; - характеризует влияние парного взаимодействия входных факторов и на выходной показатель .
Найдем значение коэффициентов по формулам:
,
,
где - значение выходного показателя в i - ом опыте (i = 1, 2,…N); N - число опытов в плане эксперимента; - значение фактора в i - ом опыте.
Расширение и преобразование исходной (заданной) матрицы планирования
Для удобства подсчета всех коэффициентов вj составим расширенную матрицу планирования (табл. 2).
Таблица 2 - Расширенная матрица планирования.
№
опыта
+1
-1
-1
+1
+1
+1
-1
-1
+1
-1
+1
-1
+1
+1
+1
+1
Найдем коэффициенты «вj» с помощью формул:
Для удобства воспользуемсяся матричной записью
Для :
∙ = .
Для :
∙ = .
Для :
∙ = .
Для :
∙ = .
Следовательно, полная математическая модель имеет вид:
Проверим существенности влияния (значимости)коэффициентов для уравнения регрессии (со статистической надежностью 95%).
Проверка существенности влияния (значимости) коэффициентов «bj» для уравнения регрессии проводит с помощью t - критерия Стьюдента по формуле:
где tj- критерий Стьюдента для «bj» коэффициента; «bj» - коэффициент уравнения;
Sв- показатель, характеризующий ошибку в определении коэффициентов «bj» и определяемый по формуле:
где Sу – стандартное отклонение, характеризующее погрешность результатов измерений выходного показателя у (величина была определена ранее Sу =2); N – число опытов в плане эксперимента (N=4).
Найдем : .
Рассчитаем опытные значения критерия Стьюдента tjопытное по формуле :
Находим критерий Стьюдента каждого коэффициента.
Если величина tопыт. меньше величины tтабл., то соответствующий коэффициент «вj» можно считать незначимым и им можно пренебречь в уравнении регрессии.
Если величина tопыт. больше величины tтабл., то соответствующий коэффициент «вj» можно считать значимым и его нужно учитывать в уравнении.
Для уравнения значимости р = 0,05; P = 0,95; f = ; t табл. (f) = .
Сравним величины t2 опытное с величиной t табл. (приложение 1) (t табл. = ).
Таким образом, путем обработки результатов эксперимента методом регрессионного анализа и после проверки значимости коэффициентов уравнения регрессии (с помощью t – критерия Стьюдента) была получена математическая модель.
Уравнение регрессии в окончательном виде:
Доказательство соответствия (адекватности) полученной математической модели (уравнения регрессии) исследуемому объекту
Докажем соответствие (адекватность) полученной математической модели (уравнения регрессии) исследуемому объекту путем расчета величины F - критерия Фишера
Дисперсия адекватности (или остаточная дисперсия)
где l - число коэффициентов в уравнении регрессии (при линейной регреcсии l = 2),
Сведем полученные расчетные значения выходного показателя в табл. 3.
Таблица 3 - Расчетные значения выходного показателя.
№
опыта
Входные факторы
Выходной показатель
X1
Х2
уэкспер.
урасчётн.
-1
-1
+1
-1
-1
+1
+1
+1
Дисперсия воспроизводимости найдена ранее:
Sy =
=
Полученное значение сравниваем с Fтабл, которое находим по таблице №2 приложения, задаваясь статистической надежностью р = 95% и учитывая величины f1 и f2.
Определим табличное значение Fтабл. по таблице 2 приложения, задаваясь статической надежностью ρ = 95% и учитывая величину
f1 = N - l = 4 -……=……
и величину
f2 = m – 1 = 3 – 1 = 2
fтабл. (f1=1, f2=2, ρ = 95% ) =
Если
то уравнение адекватно исследуемому объекту.
Поскольку величина Fonыт меньше величины Fтабл, то полученное уравнение адекватно исследуемому объекту.