Учитывая эти формулы, из равенства (294) получаем выра жения для поправок:
ЛЬ0 |
(</0 — b0) b N 0a . |
|
|
= |
N -у ДN0a |
’ |
|
|
ь |
Aft0; |
» = |
1, 2, |
----- к\ |
\ |
(299) |
Аft,, = ----- и- |
Дft, = Aft,у = 0; |
t |
/, t, |
/ ' = 1 , 2 , . . . . |
/г. |
|
Из выражений (296), (297) и (298) получаем формулы для поправок к дисперсиям и ковариациям оценок параметров:
|
As2 {bn} N |
__ |
( |
ДN0a |
|
|
(ЗСО) |
|
s 2 ( г / ) |
|
\ /V + ДЛ'0я ) |
’ |
|
|
|
Д-^2 {6») N |
А соу [Ьц, bjj] N |
_ |
/ |
ДЛУ> |
(301) |
s3 ({/1 |
|
s 3 { / / } |
|
” |
V |
N 4 - ДЛ/0а ) |
и ' |
|
Д coy |
N = |
( — Ал/°а— 'j ь |
(302) |
|
s2 {У} |
|
\ N + bN0a ) и- |
|
Дисперсии и ковариации остальных оценок параметров не изменяются при добавлении опытов в центр симметричного плана второго порядка или в центр плана Хартли.
Анализ поверхности отклика на экстремум. Следующим этапом после построения регрессионной модели является ее анализ. Регрессионную модель второго порядка обычно приводят к кано ническому виду, что существенно облегчает ее исследование. Приведение модели к каноническому виду достигается путем линейного преобразования факторов. Модель в канонической форме не содержит смешанных произведений (взаимодействий)
ивключает минимально возможное число линейных членов. Процесс приведения модели второго порядка к канонической
форме разбивается на два этапа: 1) поворот координатных осей и совмещение их с направлениями собственных векторов (новые координатные оси называют каноническими); 2) перенос начала координат в особую точку.
Первый этап позволяет исключить из уравнения взаимодей ствия, второй — свести до минимума число линейных членов. Такие преобразования обстоятельно описаны в руководствах по аналитической геометрии.
Поверхности второго порядка, приведенные к каноническому виду, можно классифицировать в зависимости од математической структуры уравнения.
На рис. 37 представлены некоторые типы поверхностей второго порядка. На рис. 38 изображены сечения эллиптического и ги перболического параболоидов.
После того, как проведено каноническое преобразование по верхности отклика и установлен его тип, производится отыскание
оптимальной точки. Это сочетается с логически профессиональным анализом характера поверхности, при проведении которого не редко используются сечения. Процесс отыскания точки оптимума во многом предопределяется тем, насколько близко расположен
Рис. 37. Поверхности второго порядка:
а — эллипсоид', б — однополостный гиперболоид; о — двухполостный гиперболоид', г — эллиптический параболоид', д — гиперболический параболоид
центр поверхности к центру плана. Если особая точка поверх ности (центр поверхности) расположена вблизи центра плана, то процесс отыскания оптимальной точки зависит от типа поверх ности. Допустим, что наша задача состоит в отыскании точки,
X |
|
2| |
|
У h h |
h.k |
|
\ _____
|
|
*! |
|
в) |
Рис. 38. Линии равных значений отклика: |
а — сечения эллиптического параболоида; б — сечения |
гиперболического параболоида; |
в — контурные линии для стационарного возвышения; г — контурные линии для возрастаю щего возвышения
отвечающей наибольшему значению у . Если поверхность унимо-
дельна (имеет единственную точку максимума), тогда эту точку можно найти исходя из необходимых условий экстремума (при
равнивая частные производные ■^ нулю). В альтернативном
случае, когда поверхность отклика является поверхностью типа минимакс, т. е. относится к классу гиперболических параболоидов, точка максимума может быть найдена только при наложении
определенных ограничений на переменные. На практике такие ограничения чаще всего задают в виде куба или сферы некоторых размеров, центры которых совмещены с особой точкой. Поиск условного экстремума производится с применением метода не определенных множителей Лагранжа или поисковых процедур нелинейного программирования. Если центр поверхности (особая точка) существенно удален от центра плана, то для отыскания точки максимума применяются следующие приемы.
Первый ■— поиск условного максимума при ограничениях в виде куба или сферы. Однако в отличие от предыдущей ситуации такие ограничения надо формировать, принимая за центр ограни
|
|
|
|
|
|
|
|
|
|
|
чивающей |
области |
не |
особую |
Таблица 85 |
|
|
точку |
поверхности, |
а |
центр |
|
|
плана. Для поиска условного |
Уровни факторов для |
|
экстремума |
при |
ограничениях |
центровочных сверл |
|
в виде сферы разработан спе |
|
|
Факторы |
|
циальный |
вычислительный |
ме |
Уровни
|
2Ф |
а |
к |
тод «риджанализ» |
[65]. |
|
|
|
|
|
Второй |
прием |
состоит |
в |
|
X2 |
л'з |
|
|
движении |
вдоль |
канонической |
|
|
|
|
оси, которая определяет напра |
+ i |
151,5 |
36,0 |
1.4 |
вление |
наиболее |
сильного |
воз |
0 |
147,5 |
32,0 |
1,2 |
растания функции отклика. |
—1 |
143,5 |
28,0 |
1 |
Иногда |
исследуемое явление |
|
|
|
|
характеризуется |
несколькими |
|
|
|
|
показателями. В подобных ситуациях прибегают к одновременному построению нескольких моделей второго порядка (по каждому из показателей). Затем возникает задача отыскания наилучшего режима. Один из путей ее решения состоит в том, что выбирается
наиболее важный показатель, скажем, у г и затем формулируется
задача отыскания точки х , отвечающей максимуму у ъ |
при усло |
вии, что наложены ограничения на факторы х г, х 2, |
■ |
хк и |
отклики у о, г/3, . . . Решение таких задач проводится с исполь
зованием метода неопределенных множителей Лагранжа и раз личных поисковых методов нелинейного программирования.
Количественному исследованию поверхности отклика второго порядка посвящено много работ. Детальное и строгое описание методов нелинейного программирования содержится в работе [51 ].
Исследование режущего инструмента с помощью планов вто рого порядка. Решалась задача установления зависимости пока зателя стойкости центровочных сверл диаметром 4 мм от трех факторов: угол в плане 2ср, задний угол а, толщина сердцевины k.
Факторы и их уровни представлены в табл. 85.
Для получения математической модели зависимости показа теля стойкости сверл от параметров в виде равенства (190) был реализован ортогональный план второго порядка для k — 3
согласно табл. 85. Испытания повторяли по 3 раза. Условия испы-
танин — см. с. 169. План и результат испытаний представлены
в табл. 86. Там же-приведены расчетные стойкости yv по полу
ченной модели (стойкость дана в количестве просверленных отверстий).
Таблица 86
Ортогональный план второго порядкадля k = 3 с результатами опытов а оценок при испытаниях центровочных сверл
S |
|
vi |
|
,v0 |
|
л’з |
|
«V |
|
1 |
—1 |
—1 |
|
—1 |
|
1627 |
1512 |
2 |
- И |
—1 |
|
—1 |
|
535 |
655 |
3 |
—1 |
-м |
|
—1 |
|
942 |
1046 |
4 |
+ 1 |
-1-1 |
|
—1 |
- |
521 |
440 |
5 |
—1 |
—1 |
|
+1 |
1421 |
1466 |
6 |
+ |
1 |
—1 |
|
+ 1 |
|
968 |
917 |
7 |
—1 |
+ 1 |
|
И |
|
1712 |
1471 |
8 |
+ |
1 |
+ |
1 |
|
+ 1 |
|
919 |
998 |
9 |
—1,215 |
0 |
|
0 |
|
1415 |
1575 |
10 |
+ |
1,215 |
0 |
|
0 |
|
880 |
819 |
11 |
|
0 |
—1,215 |
|
0 |
|
1145 |
1139 |
12 |
|
0 |
т-1,215 |
|
0 |
|
855 |
960 |
13 |
|
0 |
0 |
|
—1,215 |
1620 |
1516 |
14 |
|
0 |
0 |
|
+ 1,215 |
1626 |
1827 |
15 |
|
0 |
0 |
|
0 |
|
1694 |
1542 |
Оценки коэффициентов регрессии находили по |
формулам |
(241) с учетом численных |
значений моментов и вспомогательных |
коэффициентов |
из табл. |
73. |
|
|
|
|
|
|
|
|
15 |
_ |
|
15 _ |
|
|
|
|
|
Ьо |
о = 1 |
|
о =1 |
|
1192; |
|
|
|
|
45 ~ |
|
15 |
|
|
|
|
|
|
|
|
|
|
|
|
Ь„ = |
|
|
15 |
|
|
|
|
|
3,45- |
45 |
Е З * ?о"/Лм |
|
|
|
|
|
|
|
0 = 1 |
|
|
|
|
Ьа = 233,9; |
Ь22 = — 333,7; |
Ью = 8,0; |
|
~10,954
|
Ьг = —311,3; |
Ь2 = — 73,7; |
^ = 1 2 8 ,0 ; |
|
Ь ц = = о |
3 X i vX j gy 0t |
|
|
|
|
|
0=1 |
|
|
|
|
*„ = 41,4; |
b13 = 33,4; |
b23= |
117,6. |
Оценку |
коэффициента |
b 0, |
входящего |
в |
исходную модель, |
находим по формуле (243): Ь0 = 1192— 0,7303 |
(233,9 — 333,7 + |
+ 88,0) = |
1201. |
|
|
|
|
|