Планирование и техника эксперимента. Учебное пособие
.pdf
План и результаты полного двухфакторного эксперимента ( N 22 )
|
|
|
y |
№ точки |
x1 |
x2 |
|
|
|
|
y1 |
1 |
+1 |
+1 |
|
2 |
-1 |
+1 |
y2 |
3 |
+1 |
-1 |
y3 |
4 |
-1 |
-1 |
y4 |
|
|
|
|
План и результаты полного трехфакторного эксперимента ( N 23 ) |
|
||||
№ точки |
x1 |
x2 |
x3 |
|
y |
1 |
+1 |
+1 |
+1 |
|
y1 |
2 |
-1 |
+1 |
+1 |
|
y2 |
3 |
+1 |
-1 |
+1 |
|
y3 |
4 |
-1 |
-1 |
+1 |
|
y4 |
5 |
+1 |
+1 |
-1 |
|
y5 |
6 |
-1 |
+1 |
-1 |
|
y6 |
7 |
+1 |
-1 |
-1 |
|
y7 |
8 |
-1 |
-1 |
-1 |
|
y8 |
|
|
|
|
|
|
Отметим следующие два важных свойства рассмотренных планов: 1. Оси симметричны относительно центра эксперимента:
N |
|
|
хij 0; |
j=1, 2, 3,…N. |
|
j 1 |
i=1, 2, 3,…n; |
|
|
|
|
Сумма значений любого столбца равна нулю. |
||
2. Планы нормированы: |
|
|
N |
|
|
хij |
2 N; |
j= 0, 1, 2, 3,…N. |
j 1 |
i= 0, 1, 2, 3,…n; |
|
|
|
|
Сумма квадратов элементов любого столбца равна числу опытов.
Если известно, что изучаемую область факторного пространства можно описать линейной функцией, то в этом случае для построения модели нет необходимости ставить опыты по плану ПФЭ, так как экспериментальные затраты значительно превосходят затраты, необходимые для получения линейной модели. В подобных ситуациях можно воспользоваться планом, представляющим собой часть ПФЭ, который называется дробною факторным экспериментом (ДФЭ) или дробной репликой.
Эффективность применения дробных реплик зависит от удачного выбора системы смешивания линейных эффектов с эффектами взаимодейст-
51
вий. Реплики, у которых линейные эффекты смешаны с взаимодействиями наивысшего порядка, являются наиболее эффективными, они обладают наибольшей разрешающей способностью и называются главными.
Число точек в плане ДФЭ вычисляют по формуле
N 2n I ,
где n общее число факторов; I – величина, характеризующая дробность реплики (для полуреплики I 1, для четверть реплики I = 2 и т.д.).
Перед экспериментом по планам ДФЭ необходимо установить каков риск, т.е. определить разрешающую способность ДФЭ.
Например, по таблице найдем полуреплику от полного трехфакторного эксперимента следующим образом: оставим столбцы с факторами x1 и
x2 , а x3 рассчитаем по формуле
x3 x1 : x2 .
Это соотношение называется генерирующим.
|
Полуреплика от полного трехфазного эксперимента |
|
|||
|
|
|
|
|
x3 = x1 * x2 |
№ точки |
|
x1 |
x2 |
||
1 |
|
+1 |
+1 |
|
+1 |
2 |
|
-1 |
+1 |
|
-1 |
3 |
|
+1 |
-1 |
|
-1 |
4 |
|
-1 |
-1 |
|
+1 |
|
|
|
|
|
|
Композиционные планы второго порядка, представляют собой достройку планов первого порядка выполненную таким образом, чтобы по ним можно было получать модели второго порядка.
Алгоритм для полного факторного эксперимента 22 с условиями:
Уровень фактора |
А |
В |
|
|
|
Основной xi = 0 |
|
|
|
|
|
Интервал варьирования xi |
|
|
|
|
|
Нижний xi = -1 |
|
|
Верхний xi = +1 |
|
|
Кодовые обозначения |
X1 |
X2 |
и m=3 параллельными опытами.
52
1. Строят план эксперимента в кодированных элементах:
|
№ |
|
X1 |
|
X2 |
|
X1 |
|
Параллельные опы- |
|
уj |
|
j |
2 |
|
y j |
|
|
2 |
|
||
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
(y j y j ) |
|
|
||||||||||
|
опыта |
|
|
|
|
|
X2 |
|
|
ты |
|
|
|
|
|
|
|
|
|
|
|
|
|
j |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Y1 |
Y2 |
Y3 |
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
+ |
|
+ |
|
+ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
- |
|
+ |
|
- |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
3 |
|
+ |
|
- |
|
- |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
4 |
|
- |
|
- |
|
+ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2. |
Определяют y j и выборочные дисперсии j2 |
для каждого опыта по |
|||||||
|
|
1 |
|
m |
|
|
|
|
|
формуле 2j |
|
yij y j 2 . |
|
|
|
|
|||
m 1 |
|
|
|
|
|||||
|
|
i 1 |
|
|
|
|
|||
|
|
|
|
|
|
|
|
||
3. |
Проводят проверку воспроизводимости опытов. |
|
|||||||
|
|
|
|
|
|
1 |
|
m |
|
4. |
Находят дисперсию воспроизводимости 2j |
|
|
yij y j 2 . |
|||||
m 1 |
|||||||||
|
|
|
|
|
i 1 |
||||
|
|
|
|
|
|
|
|
||
Степень свободы дисперсии воспроизводимости равна f = N(m – 1). |
|||||||||
5. |
Определяем коэффициенты уравнения регрессии, которое в общем |
||||||||
случае имеет вид |
|
|
|
|
|
||||
yb0 b1x1 b2 x2 b12 x1x2 .
6.Для нахождения коэффициентов b0, b1, b2 и b12 используем соответственно зависимости:
|
|
|
|
|
|
1 |
N |
|
|
|
|
|
|
|
|
|
|
|
|
|
b0 |
y j |
|
|
|
|
y1 y2 ... yN |
|
|
||
|
|
|
|
|
N |
или |
b |
|
|||||||
|
|
|
|
|
|
N |
|||||||||
|
|
|
|
|
|
|
j 1 |
0 |
|
||||||
|
|
|
1 |
|
N |
|
|
|
xi1y1 xi2 y2 ... xiN yN |
|
|||||
bi |
|
xij y j |
|
bi |
|
||||||||||
N |
или |
||||||||||||||
|
|
||||||||||||||
|
|
|
|
|
j 1 |
|
|
|
|
N |
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
1 |
|
N |
|
|
|
|
|
|
|
|
|
|
||
bik |
xij xkj y j |
|
b |
xi1xk1 y1 xi2 xk2 y2 ... xiN xkN yN |
|||||||||||
N |
|
||||||||||||||
|
|
|
|||||||||||||
|
|
j 1 |
|
|
или |
ik |
|
|
N |
||||||
|
|
|
|
|
|
|
|
|
|
|
|||||
7. Проверяют гипотезу о значимости коэффициентов с помощью кри-
терия Стьюдента. Если выполняется следующее условие: если bi t ,f
в2 (f = N(m -1)), то данный коэффициент является статистически значимым,
2 |
|
2у |
2 |
|
2j |
где в |
|
|
– дисперсия коэффициентов, а y |
|
N – дисперсия воспро- |
N m |
изводимости. Все коэффициенты, для которых это условие не выполняется,
53
являются незначимыми, и они из уравнения регрессии исключаются вместе с соответствующим фактором, при этом оставшиеся коэффициенты не пересчитываются.
8. Проверяют гипотезу об адекватности полученной математической модели результатом эксперимента с использованием критерия Фишера (в случае, когда гипотеза об адекватности модели отвергается, необходимо переходить к более сложной форме математического описания, либо, если это возможно, провести эксперимент с меньшим интервалом варьирования).
9. Выполняют переход от кодированных значений факторов к натуральным:
|
|
|
0 |
|
|
|
|
0 |
|
|
|
0 |
|
|
0 |
|
|
y b |
b |
|
x1 |
x1 |
|
b |
|
x2 x2 |
|
b |
|
x1 x1 |
|
x2 x2 |
. |
||
|
x1 |
|
x |
|
|||||||||||||
0 |
1 |
|
|
|
2 |
x |
2 |
|
12 |
|
|
x |
2 |
|
|||
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|||
54
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1.Ахназарова, С.Л. Методы оптимизации эксперимента в химии и химической технологии / С.Л. Ахназарова, В.В. Кафаров. – М.: Высш. шк., 1985. – 418 с.
2.Адлер, Ю.П. Планирование эксперимента при поиске оптимальных условий / Ю.П. Адлер, Е.В. Маркова, Ю.В. Грановский. – М.: Наука, 1976. – 280 с.
3.Налимов, В.В. Теория эксперимента: учебник / В.В. Налимов. – М.:
Наука, 1971. – 208 с.
4.Джонсон, Н. Статистика и планирование эксперимента в технике и науке: Методы обработки данных / Н. Джонсон, Ф. Лион. – М.: Мир, 1980. – 610 с.
5. |
Хикс, Ч. Основные принципы планирования |
эксперимента / |
Ч. Хикс. – М.: Книга по требованию, 2013. – 203 с. |
|
|
6. |
Дрейпер, Н. Прикладной регрессионный анализ |
/ Н. Дрейпер, |
Г. Смит. – М.: Диалектика-Вильямс, 2007. – 912 с. |
|
|
7. |
Планирование и техника эксперимента: электронный обучающий |
|
курс / сост. О.П. Стебелева [и др.] // Система электронного обучения СФУ e.sfu-kras.ru. – Красноярск, 2014. – Режим доступа: http://e.sfu-kras.ru/ course/view.php?id=583.
55
Учебное издание
Кулагина Татьяна Анатольевна Стебелева Олеся Павловна
ПЛАНИРОВАНИЕ И ТЕХНИКА ЭКСПЕРИМЕНТА
Учебное пособие
Подписано в печать 2.06.2017. Печать плоская. Формат 60×84/16 Бумага офсетная. Усл. печ. л. 3,3. Тираж 100 экз. Заказ № 1346
Библиотечно-издательский комплекс Сибирского федерального университета 660041, Красноярск, пр. Свободный, 82а
Тел. (391) 206-26-67; http://bik.sfu-kras.ru, E-mail: publishing_house@sfu-kras.ru
56
