Исследование уравнений регрессии, полученных с помощью полного факторного эксперимента, проводится с использованием следующих формул:
- дисперсия погрешностей коэффициентов регрессии, вычисленных по формуле (3), равна
(3)
где
-
дисперсия ошибки опытов, оцененная по
параллельным опытам в лабораторной
работе №1;
- остаточная сумма квадратов вычисляется по формуле
(4)
где
-
значение зависимой переменной в u-ой
строке матрицы планирования;
-
значение зависимой переменной,
рассчитанное по модели (1) для u-ой
строки матрицы планирования;
- адекватность модели регрессии проверяется по критерию Фишера:
(5)
Если расчетное значение F меньше табличного значения критерия Фишера, выбранного для уровня значимости α=0,05 и числа степеней свободы fr и fy, то полученная модель считается адекватной;
- значимость коэффициентов регрессии проверяется по критерию Стьюдента:
(6)
Расчетное
значение сопоставляется с критическим
значением критерия Стьюдента, выбранным
для уровня значимости α=0,05 и числа
степеней свободы, с которым была
определена ошибка эксперимента
.
Если расчетное значение критерия больше
критического, коэффициент регрессии
считается значимым.
Незначимые коэффициенты отбрасываются без повторного пересчета значимых коэффициентов регрессии, так как матрица планирования ортогональная.
Если хотя бы один коэффициент регрессии при парном взаимодействии будет значимым, то в этом случае отвергается гипотеза линейности регрессионной модели.
Пример
Объект исследования задан математической моделью вида
y=200-4*(x1-5)2-4*(x2-5)2-4*(x3-5)2+e . (7)
В заданном интервале варьирования переменных (табл.2) изучаемый объект может описываться линейной моделью.
Таблица 2 – Интервалы варьирования переменных
Факторные переменные |
х1 |
х2 |
х3 |
Основной уровень "0" |
8 |
8 |
8 |
Шаг варьирования Δxi |
1 |
1 |
1 |
Верхний уровень +1 |
9 |
9 |
9 |
Нижний уровень -1 |
7 |
7 |
7 |
Имеем
дело с тремя независимыми переменными
и каждую из переменных варьируем на
двух уровнях, условно обозначаемых +1 и
-1. Все возможные комбинации для трех
факторов, варьируемых на двух уровнях,
будут исчерпаны, если поставить восемь
опытов. Матрица планирования эксперимента
с наложением ошибок измерений «е»,
описываемых нормальным законом
распределения с дисперсией, равной
=1,
приведена в таблице 3.
Таблица 3 – Матрица планирования эксперимента с наложением
ошибок измерений
Опыты |
х1 |
х2 |
х3 |
y |
е |
y+e |
|
1 |
-1 |
-1 |
-1 |
152 |
-0,30023 |
151,6998 |
151,9496 |
2 |
1 |
-1 |
-1 |
104 |
-1,27768 |
102,7223 |
104,5843 |
3 |
-1 |
1 |
-1 |
104 |
0,244257 |
104,2443 |
103,3869 |
4 |
1 |
1 |
-1 |
56 |
1,276474 |
57,27647 |
56,02168 |
5 |
-1 |
-1 |
1 |
104 |
1,19835 |
105,1984 |
104,0923 |
6 |
1 |
-1 |
1 |
56 |
1,733133 |
57,73313 |
56,72706 |
7 |
-1 |
1 |
1 |
56 |
-2,18359 |
53,81641 |
55,52967 |
8 |
1 |
1 |
1 |
8 |
-0,23418 |
7,765819 |
8,164406 |
Примечание: При расчете зависимой переменной «у» по формуле (7) вместо кодовых обозначений факторных переменных подставлялись соответствующие им фактические значения из таблицы 2.
По формуле (3) вычисляем коэффициенты уравнения регрессии:
b0 = 80,057; b1 = -23,68; b2 = -24,28; b3 = -23,93.
В результате модельного эксперимента получено уравнение регрессии, в котором факторные переменные представлены в кодированном виде:
y = 80,057 – 23,68x1 - 24,28 x2 – 23,93 x3 . (8)
Проведем исследование уравнения регрессии, полученного с помощью полного факторного эксперимента.
1) Дисперсия погрешностей коэффициентов регрессии, вычисленная по формуле (3), равна:
2) Проверяем значимость коэффициента регрессии b1 по критерию Стьюдента (6):
Критерий Стьюдента, выбранный для уровня значимости α=0,05 и числа степеней свободы 7, равен t=2,36.
Расчетное значение критерия больше критического, что подтверждает значимость коэффициент регрессии b1. Проверка остальных коэффициентов подтвердила их значимость.
По матрице планирования были вычислены коэффициенты парных взаимодействий, которые приняли значения b12= 0,428; b13= 0,303. Проверка по критерию Стьюдента показала их статистическую незначимость.
Незначимые коэффициенты отбрасываются без повторного пересчета значимых коэффициентов регрессии, так как матрица планирования ортогональная.
3) Коэффициенты при независимых переменных указывают на силу влияния фактора. Чем больше численная величина коэффициента, тем большее влияние оказывает фактор. Коэффициенты в уравнении регрессии (8) имеют знак минус, что указывают на обратную связь между факторами и зависимой переменной. Величина коэффициента соответствует вкладу данного фактора в величину зависимой переменной при переходе фактора с нулевого на верхний или нижний уровень.
4)
Для определения остаточной суммы
квадратов модели по уравнению регрессии
(8) вычисляются значения зависимой
переменной
для
всех строк матрицы планирования.
Результаты вычислений записаны в
последнем столбце таблицы 3. Остаточная
сумма квадратов модели, рассчитанная
по формуле
(4), составила величину
5) Адекватность модели регрессии проверяется по критерию Фишера:
Табличное значение критерия Фишера, выбранное для уровня значимости α=0,05 и числа степеней свободы fr=(8-3-1)=4 и fy=7 равно Fкр=4,12. Расчетное значение F меньше Fкр, поэтому полученная модель (8) считается адекватной.
1) Ознакомиться с методикой планирования полного факторного эксперимента и анализом полученного уравнения регрессии.
2) Получить у преподавателя вариант задания.
3) Составить матрицу планирования полного факторного эксперимента.
4) Реализовать модельный эксперимент, вычислить коэффициенты уравнения регрессии.
5) Исследовать уравнение регрессии, полученного с помощью полного факторного эксперимента:
- оценить дисперсию погрешностей коэффициентов регрессии;
- проверить значимость коэффициентов регрессии при линейных членах и парных взаимодействиях факторов;
- определить остаточную сумму квадратов модели регрессии;
- оценить адекватность модели регрессии.
6) Сделать выводы по результатам модельного эксперимента.
7) Оформить отчет по выполненной работе.
1) Название, цель работы и номер варианта задания.
2) Математическая модель объекта и интервалы варьирования переменных.
3) Матрица планирования эксперимента с наложением ошибок измерений.
4) Вычисления коэффициентов уравнения регрессии, оценка дисперсии погрешности коэффициентов уравнения регрессии.
5) Проверка значимости линейных и парных коэффициентов регрессии по критерию Стьюдента.
6) Расчетное определение остаточной суммы квадратов модели. Оценка адекватности модели.
7) Выводы по результатам выполненной лабораторной работы.
1) Что называется полным факторным экспериментом?
2) Какова процедура поиска коэффициентов линейной модели?
3) Как численно оцениваются эффекты взаимодействия?
4) В чем состоит исследование уравнений регрессии, полученных с помощью полного факторного эксперимента?
5) Как оценивается дисперсия погрешностей коэффициентов регрессии и их значимость?
6) Как рассчитывается остаточная сумма квадратов модели и проверяется адекватность линейности регрессионной модели?
1. В.В. Налимов, Н.А. Чернова. Статистические методы планирования экстремальных экспериментов. Изд. «Наука», Москва, 1965. – 340с.
2. Советов Б.Я., Яковлев С.А. Моделирование систем: Учеб. для вузов-3-е изд., перераб. и доп. –М.:Высш. Шк.., 2001.-343с. ISBN 5-06-003860-2.
Освоить методику проведения полного факторного эксперимента и исследование уравнений регрессии на модели объекта.
Вариант |
Модель объекта |
Диапазон изменения факторов |
Основной уровень |
Шаг изменения факторов |
Дисперсия ошибки измерений |
1 2 3 |
y=300-4*(x1-5)2-4*(x2-5)2-4*(x3-5)2 |
0-10 0-10 0-10 |
7 8 9 |
1,0 1,0 1,0 |
1,0 1,2 0,8 |
4 5 6 |
y=300-4*(x1-7)2-4*(x2-7)2-4*(x3-7)2 |
2-12 2-12 2-12 |
9 10 11 |
1,0 1,0 1,0 |
1,0 1,2 0,8 |
7 8 9 |
y=350-3*(x1-6)2-3*(x2-6)2-3*(x3-6)2 |
0-12 0-12 0-12 |
8 9 10 |
1,0 1,0 1,0 |
1,0 1,2 0,8 |
10 11 12 |
y=200-4*(x1-5)2-4*(x2-5)2-4*(x3-5)2 |
0-10 0-10 0-10 |
7 8 9 |
1,0 1,0 1,0 |
1,0 1,2 0,8 |