Планирование и организация эксперимента. Лабораторный практикум. Учебное пособие
.pdf
Окончание табл. 44
Вариант № 12
N |
x1 |
x2 |
x3 |
yэкс |
yэкс1 |
yэкс2 |
1 |
1 |
0 |
0 |
y1 |
68,72 |
69,88 |
2 |
0 |
1 |
0 |
y2 |
54,00 |
53,21 |
3 |
0 |
0 |
1 |
y3 |
30,26 |
29,35 |
4 |
1/2 |
1/2 |
0 |
y12 |
43,39 |
44,12 |
5 |
1/2 |
0 |
1/2 |
y13 |
33,34 |
32,36 |
6 |
0 |
1/2 |
1/2 |
y23 |
59,45 |
58,25 |
Таблица 45
Вид графика проекций линий равных значений свойств
1 |
Ternary Graph (Spreadsheet4 in Workbook001 4v*6c) |
|
|
Var4 = 65,85*x+42,94*y+20,905*z-154,08*x*y-10,01*x*z+56,39*y*z |
|
|
Var3 |
|
|
0,00 |
1,00 |
|
0,25 |
0,75 |
0,50 |
0,50 |
0,75
0,25
1,00 |
|
|
|
0,00 |
0,00 |
0,25 |
0,50 |
0,75 |
1,00 |
Var1 |
|
|
|
Var2 |
60
50
40
30
20
81
Продолжение табл. 45
Вид графика проекций линий равных значений свойств
2
3
Ternary Graph (Spreadsheet4 in Workbook001 4v*6c) |
|
Var4 = 75,97*x+71,135*y+75,02*z-131,27*x*y-202,24*x*z-125,53*y*z |
|
Var3 |
|
0,00 |
1,00 |
0,25 |
0,75 |
0,50
0,50
|
0,75 |
|
|
0,25 |
|
|
|
|
|
|
|
70 |
|
|
|
|
|
|
60 |
|
1,00 |
|
|
|
0,00 |
50 |
|
|
|
|
40 |
|||
0,00 |
0,25 |
0,50 |
0,75 |
1,00 |
||
30 |
||||||
Var1 |
|
|
|
Var2 |
||
|
|
|
|
Ternary Graph (Spreadsheet4 in Workbook001 4v*6c) |
|
Var4 = 27,09*x+73,625*y+30,175*z-74,05*x*y+49,85*x*z-38,94*y*z |
|
Var3 |
|
0,00 |
1,00 |
0,25 |
0,75 |
0,50
0,50
|
0,75 |
|
|
0,25 |
|
|
|
|
|
|
|
70 |
|
|
|
|
|
|
60 |
|
1,00 |
|
|
|
0,00 |
50 |
|
|
|
|
40 |
|||
0,00 |
0,25 |
0,50 |
0,75 |
1,00 |
||
30 |
||||||
Var1 |
|
|
|
Var2 |
||
|
|
|
|
82
Продолжение табл. 45
Вид графика проекций линий равных значений свойств
4 |
Ternary Graph (Spreadsheet4 in Workbook001 4v*6c) |
Var4 = 21,21*x+26,895*y+61,235*z+130,89*x*y+27,77*x*z+122,74*y*z
Var3
0,00 |
1,00 |
0,25 |
0,75 |
0,50
0,50
|
0,75 |
|
|
|
0,25 |
|
|
|
|
|
|
|
|
70 |
|
|
|
|
|
|
|
60 |
|
1,00 |
|
|
|
|
0,00 |
50 |
|
|
|
|
|
40 |
|||
0,00 |
0,25 |
0,50 |
|
0,75 |
1,00 |
||
|
30 |
||||||
Var1 |
|
|
|
|
Var2 |
||
|
|
|
|
|
|||
5 |
Ternary Graph (Spreadsheet4 in Workbook001 4v*6c) |
|
|||||
|
Var4 = 69,3*x+53,605*y+29,805*z-70,79*x*y-66,81*x*z+68,58*y*z |
|
|||||
|
|
Var3 |
|
|
|
||
|
|
0,00 |
1,00 |
|
|
|
|
|
0,25 |
|
|
0,75 |
|
|
|
|
0,50 |
|
|
0,50 |
|
|
|
|
0,75 |
|
|
|
0,25 |
65 |
|
|
|
|
|
60 |
|||
|
|
|
|
|
|
||
|
|
|
|
|
|
55 |
|
|
|
|
|
|
|
50 |
|
|
|
|
|
|
|
45 |
|
1,00 |
|
|
|
|
0,00 |
40 |
|
|
|
|
|
35 |
|||
0,00 |
0,25 |
0,50 |
0,75 |
1,00 |
|||
30 |
|||||||
Var1 |
|
|
|
|
Var2 |
||
|
|
|
|
|
|||
83
Продолжение табл. 45
Вид графика проекций линий равных значений свойств
6
7
Ternary Graph (Spreadsheet4 in Workbook001 4v*6c) |
|
Var4 = 19,015*x+57,865*y+26,025*z-25,44*x*y+105,12*x*z+23,42*y*z |
|
Var3 |
|
0,00 |
1,00 |
|
0,25 |
|
0,75 |
|
|
|
|
0,50 |
|
0,50 |
|
|
|
|
0,75 |
|
|
0,25 |
|
55 |
|
|
|
|
50 |
||
|
|
|
|
|
|
|
|
|
|
|
|
|
45 |
|
|
|
|
|
|
40 |
|
|
|
|
|
|
35 |
1,00 |
|
|
|
|
0,00 |
30 |
|
|
|
|
25 |
||
0,00 |
0,25 |
0,50 |
0,75 |
|
1,00 |
|
|
20 |
|||||
Var1 |
|
|
|
|
Var2 |
|
|
|
|
|
|
Ternary Graph (Spreadsheet4 in Workbook001 4v*6c) |
|
Var4 = 39,1*x+36,86*y+60,165*z+78,8*x*y-62,65*x*z-1,23*y*z |
|
Var3 |
|
0,00 |
1,00 |
0,25 |
0,75 |
0,50
0,50
|
0,75 |
|
|
0,25 |
|
|
|
|
|
|
|
60 |
|
|
|
|
|
|
55 |
|
|
|
|
|
|
50 |
|
1,00 |
|
|
|
0,00 |
45 |
|
|
|
|
40 |
|||
0,00 |
0,25 |
0,50 |
0,75 |
1,00 |
||
35 |
||||||
Var1 |
|
|
|
Var2 |
||
|
|
|
|
84
Продолжение табл. 45
Вид графика проекций линий равных значений свойств
8
9
Ternary Graph (Spreadsheet4 in Workbook001 4v*6c) |
|
Var4 = 35,29*x+37,725*y+55,415*z+58,17*x*y+38,05*x*z+32,58*y*z |
|
Var3 |
|
0,00 |
1,00 |
0,25 |
0,75 |
0,50
0,50
|
0,75 |
|
|
0,25 |
|
|
|
|
|
|
|
58 |
|
|
|
|
|
|
54 |
|
|
|
|
|
|
50 |
|
1,00 |
|
|
|
0,00 |
46 |
|
|
|
|
42 |
|||
0,00 |
0,25 |
0,50 |
0,75 |
1,00 |
||
38 |
||||||
Var1 |
|
|
|
Var2 |
||
|
|
|
|
Ternary Graph (Spreadsheet4 in Workbook001 4v*6c) |
|
Var4 = 60,66*x+23,235*y+45,14*z+66,91*x*y-47,96*x*z+180,21*y*z |
|
Var3 |
|
0,00 |
1,00 |
0,25 |
0,75 |
0,50
0,50
|
0,75 |
|
|
0,25 |
|
|
|
|
|
|
|
70 |
|
|
|
|
|
|
60 |
|
1,00 |
|
|
|
0,00 |
50 |
|
|
|
|
40 |
|||
0,00 |
0,25 |
0,50 |
0,75 |
1,00 |
||
30 |
||||||
Var1 |
|
|
|
Var2 |
||
|
|
|
|
85
Продолжение табл. 45
Вид графика проекций линий равных значений свойств
10
11
Ternary Graph (Spreadsheet4 in Workbook001 4v*6c) |
|
Var4 = 44,6*x+81,89*y+54,855*z-97,66*x*y-99,87*x*z-85,75*y*z |
|
Var3 |
|
0,00 |
1,00 |
0,25 |
0,75 |
0,50
0,50
|
0,75 |
|
|
0,25 |
|
|
|
|
|
|
|
80 |
|
|
|
|
|
|
70 |
|
|
|
|
|
|
60 |
|
1,00 |
|
|
|
0,00 |
50 |
|
|
|
|
40 |
|||
0,00 |
0,25 |
0,50 |
0,75 |
1,00 |
||
30 |
||||||
Var1 |
|
|
|
Var2 |
||
|
|
|
|
Ternary Graph (Spreadsheet4 in Workbook001 4v*6c) |
|
Var4 = 41,94*x+64,28*y+62,89*z-135,92*x*y-58,2*x*z-86,18*y*z |
|
Var3 |
|
0,00 |
1,00 |
0,25 |
0,75 |
0,50
0,50
|
0,75 |
|
|
0,25 |
|
|
|
|
|
|
|
60 |
|
|
|
|
|
|
50 |
|
1,00 |
|
|
|
0,00 |
40 |
|
|
|
|
30 |
|||
0,00 |
0,25 |
0,50 |
0,75 |
1,00 |
||
20 |
||||||
Var1 |
|
|
|
Var2 |
||
|
|
|
|
86
Окончание табл. 45
Вид графика проекций линий равных значений свойств
12 |
Ternary Graph (Spreadsheet4 in Workbook001 4v*6c) |
Var4 = 50,94*x+13,805*y+36,765*z+10,69*x*y+19,73*x*z+169,64*y*z Var3
0,00 |
1,00 |
0,25 |
0,75 |
0,50
0,50
|
0,75 |
|
|
0,25 |
1,00 |
|
|
|
0,00 |
0,00 |
0,25 |
0,50 |
0,75 |
1,00 |
Var1 |
|
|
|
Var2 |
60
50
40
30
20
Содержание отчета
Отчет оформляется так же, как в лабораторной работе № 1.
Лабораторная работа № 9. Планирование экспериментов при поиске оптимальных условий
Цель работы: изучение метода поиска оптимальных условий проведения эксперимента.
Краткие теоретические сведения
Во многих случаях перед исследователем возникает задача не только выявления характера связи между двумя или несколькими рядами наблюдений, но и нахождения таких численных значений факторов, при которых отклик достигает своего экстре-
87
мального значения (максимума или минимума). Эксперимент, решающий эту задачу, называется экстремальным.
В этом случае задача сводится к оптимизационной и формулируется следующим образом: требуется определить такие координаты экстремальной точки (х1*, х2*, …, xk*) поверхности отклика y = f (х1, х2, …, xk), в которой она максимальна (или минимальна).
Известно, что наиболее короткий путь – это движение по градиенту, т. е. перпендикулярно линиям равного уровня, на которых функция отклика принимает постоянные значения f (х1, х2, …, xk) = const. В связи с этим, при оптимизации процесса рабочее движение целесообразно совершать в направлении наиболее быстрого возрастания функции отклика, т. е. в направлении градиента функции у.
Пусть в окрестности точки А (центр плана) поставлен ПФЭ 22 и получено уравнение y = b0 + b1x1 + b2x2. (рис. 20).
Кратчайшее расстояние к оптимуму – по направлению градиента функции отклика
|
¶y |
r |
¶y |
r |
¶y |
r |
, |
|
grad y = |
|
i + |
|
j + ...+ |
|
k |
||
¶x1 |
¶x2 |
¶xk |
||||||
|
|
|
|
|
где i, j, …, k – единичные векторы в направлении факторных осей;
¶y – частная производная функции отклика по i-му фактору.
¶xi
у
х2
х1
А
Рис. 20
88
Движение в направлении градиента осуществляется путем изменения факторов пропорционально величинам коэффициентов регрессии. Коэффициент уравнения регрессии bi равен tg угла между линией регрессии и осью xi. Коэффициент b0 не оказывает влияние на расчет градиента. Шаг движения может быть выбран пропорционально рассчитанным составляющим градиента. Расчет новых точек ведут в натуральном масштабе.
Рекомендуется первоначально проделать мысленные опыты (расчет по уравнению регрессии), которые позволяют:
1)сократить объем реальных опытов;
2)оценить правильность выбора шага движения. Последовательность движения к оптимуму по направлению
градиента функции отклика (метод крутого восхождения) (рис. 20):
1)выбор нулевой точки А0;
2)расчет уравнения регрессии по результатам ПФЭ или ДФЭ;
3)расчет составляющих градиента и рабочих шагов;
4)постановка мысленных и реальных опытов в точках, полученных прибавлением к А0 рассчитанных шагов;
5)определение точки частного экстремумаВ0, которую следует принять в качестве новой нулевой точки;
6)организация второго цикла крутого восхождения(аналогично);
7)при необходимости выполнить третий цикл.
Признаком достижения точки экстремума является статистическая незначимость оценок коэффициентов в уравнении первого порядка.
Пример выполнения задания.
В результате ПФЭ 22 получено линейное уравнение регрессии в кодированных значениях
ˆy = 88 – 2x1 – 4,5x2.
Натуральные уровни факторов принимают следующие значения:
основной уровень |
х1, усл. ед. х2, % |
|
4,5 |
10,0 |
|
интервал варьирования (I) |
0,5 |
1,0 |
верхний уровень |
5,0 |
11,0 |
нижний уровень |
4,0 |
9,0 |
89
х2
|
|
|
|
В0 |
|
|
|
|
А3 |
А4 |
|
|
|
|
|
|
|
|
|
|
А0 |
|
|
|
|
|
А1 |
А2 |
х1 |
|
|
|
Рис. 21 |
|
|
Результаты опытов следующие (средние значения функции |
|||||
отклика по двум параллельным опытам): |
|
|
|||
N |
х0 |
х1 |
х2 |
у, °С |
|
1 |
+1 |
-1 |
-1 |
95,0 |
|
2 |
+1 |
+1 |
-1 |
90,0 |
|
3 |
+1 |
-1 |
+1 |
85,0 |
|
4 |
+1 |
+1 |
+1 |
82,0 |
|
Решение.
Рассчитаем составляющие градиента
b1 I1 = - 1,0 усл. ед., b2 I2 = - 4,5 %.
Преобразуем уравнение регрессии из кодированных в натуральные значения факторов с использованием следующего выражения
|
~ |
~ |
|
|
|
|
|
|
|
|
x |
= |
xi |
- xi0 |
, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
i |
|
|
Ii |
|
|
|
|
|
|
|
|
|
|
~ |
- 4,5 |
|
~ |
-10,0 |
|
||
ˆy |
= 88 – 2x1 – 4,5x2 = 88 – 2 |
x1 |
– 4,5 |
x2 |
= |
|||||
|
0,5 |
|
1,0 |
|||||||
|
|
|
|
|
|
|
|
|
||
~~
=151 – 4 x1 – 4,5 x2 .
90
