Файл: Кацев, П. Г. Статистические методы исследования режущего инструмента.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 21.10.2024
Просмотров: 138
Скачиваний: 1
Эффект фактора будет значимым при tt > |
tKP с числом степеней |
свободы |
|
f = S ni — k, |
|
где k — число клеток в табл. 67. |
£0 05 = 2,365, ^01 = |
В нашем случае критическое значение |
=1,895.
Таким образом, оба выделенных эффекта значимы с вероят
ностью 95%.
После оценивания эффектов нескольких факторов произво дится корректировка результатов экспериментов для более чет кого выделения менее сильных факторов и парных взаимодействий. Корректировка производится путем прибавления с обратным знаком эффектов выделенных факторов к результатам отсеива ющих экспериментов: +6,0 ко всем результатам на уровне Х х + ;
— 13,7 ко всем результатам на уровне Х 4 + .
«Снимаем» действие выделенных эффектов Х 4 и Х 4. Корректи рованные результаты приведены в табл. 66 в графе у г. По скоррек
тированным результатам вновь строим диаграммы рассеяния (рис. 33), где помимо линейных эффектов приведены и эффекты парных взаимодействий. Оценивание эффектов и корректировка результатов производится до тех пор, пока эффекты факторов окажутся незначимыми для 10% уровня значимости. В нашем примере отсеивание факторов прекратили на третьем этапе (зна чения у 3 в табл. 66). После выделения факторов Х 7, Х 8 и оценки
их значимости по критерию Стьюдента убедились, что оба эффекта этих факторов незначимы (/8 = 0,87 и t7 = 0,79, критическое
значение ^0il = 1,865).
185
После серии корректировок результатов экспериментов строим точечную диаграмму распределения значений критерия оптими зации на различных стадиях выделения, значимых эффектов (рис. 34). После снятия значимых эффектов разброс точек значи
|
тельно снижается. |
|
факто- |
|||
У01 |
Результаты |
отсеивания |
||||
ров сводим в табл. 69j |
и |
строим |
||||
|
диаграмму |
эффектов |
(рис. 35). |
|||
|
Рис. 35 дает |
наглядное |
представ |
|||
|
ление о степени влияния |
отдель |
||||
|
ных факторов |
на стойкость сверл. |
||||
30 |
Проведение |
отсеивающего экспе |
||||
римента |
по |
методу случайного |
||||
|
баланса позволило выделить че |
|||||
|
тыре линейных эффекта, сильно |
|||||
|
влияющих на стойкость сверл диа |
|||||
g 2оL |
метром 5 |
мм. |
Кроме |
линейных |
||
эффектов выделено одно |
двойное |
|||||
|
взаимодействие |
заднего |
угла х 2 |
w
|
У, |
Уз |
|
Рис. 34. |
Точечная диаграмма |
Рис. 35. Диаграмма эффектов, |
|
распределения |
результатов на |
выделенных методом случайного |
|
блюдений |
на |
различных этапах |
баланса |
|
отсеяния |
|
с углом при вершине х3. Знак минус при эффекте этого взаимодей ствия показывает, что одновременное увеличение или уменьшение заднего угла и угла при вершине вызывает уменьшение стойкости и наоборот, — увеличение одного из факторов и уменьшение другого вызывает увеличение стойкости. Таким образом, для дальнейшего исследования надо включать пять факторов: длину рабочей части; задний угол; угол при вершине; длину поперечной кромки; толщину сердцевины.
Данная задача решалась и другим, более универсальным способом. Вместо построения таблиц с двумя входами пользова лись обычными приемами регрессионного анализа: строилась матрица X независимых переменных х1 по всем тем эффектам,
которые признаны значимыми на первом этапе по диаграммам
186
Таблица 69 |
|
|
|
|
|
Р езульт ат ы от сеивания |
ф акт оров |
|
|
||
Этап определения |
Выделенные |
Численное зна |
Расчетное зна |
||
эффекты |
чение эффектов |
чение t -крнте- |
|||
|
|
|
|
|
рия |
По исходным данным |
|
—6,0 |
2,37 |
||
После |
первой |
корректи- |
*4 |
13,74 |
5,18 |
*3 |
6,36 |
2,38 |
|||
ровки |
|
|
*5 |
3,20 |
1,91 |
|
|
|
*2* |
|
|
После |
второй |
корректи- |
*7 |
2,71 |
0,87 |
ровки |
|
|
*8 |
2,03 |
0,79 |
* хз выделен по двойному взаимодействию х2х, с эффектом - 7,645 (t23 =
= 1,97).
рассеяния. Затем с помощью ЭВМ находили коэффициенты нор мальных уравнений Х *Х и обратную матрицу {Х*Х)~1. При
этом на этапе «узнавания» по диаграммам рассеяния выделено 10 эффектов и по ним строилась матрица X. Эффекты, признанные
значимыми, снимались с использованием ветвящейся стратегии. Результаты обработки эксперимента по методу ветвящейся стра тегии оказались аналогичными с результатами ручного счета.
В более сложной задаче метод ветвящейся стратегии следует предпочесть как более эффективный. Это объясняется тем, что применение ЭВМ не только ускоряет расчеты, но и повышает разрешающую способность метода, так как позволяет одновре менно оценивать гораздо большее число эффектов. После расщеп ления модели можно использовать обычные методы факторного эксперимента, так как план из сверхнасыщенного по отношению ко всем эффектам становится ненасыщенным по отношению к зна чимым эффектам. В ряде случаев можно использовать уже вы полненный факторный эксперимент для движения к области оптимума по градиенту.
Планирование второго порядка*
[24—27, 43]
После достижения области оптимума исследователя обычно интересует более детальное описание поверхности отклика. По добное описание может быть использовано для локализации точки оптимума, а также для установления типа поверхности, что
* Раздел написан совместно с В. Г. Горским.
187
представляет большой интерес при интерпретации полученных результатов. Для выполнения этой задачи описания поверхности отклика полиномом первого порядка уже недостаточно. Во мно гих случаях вполне удовлетворительная аппроксимация может быть достигнута, если воспользоваться полиномом второго по рядка.
Математическая модель второго порядка имеет вид
k |
k |
к |
М Ы = т 1 = Р о + 2 Р л + |
2 РЦХ1Х, + |
2 Рих]. ■ (190) |
(=1 |
:</' |
<=1 |
Задача состоит в определении оценок коэффициентов этой модели по результатам спланированного эксперимента.
Основное требование, предъявляемое к плану второго порядка, состоит в том, что план должен допускать получение раздельных, не смешанных оценок коэффициентов регрессии. Для этого необ ходимо, чтобы число разных опытов было не меньше числа коэф фициентов в модели, равного числу сочетаний из /г + 2 по два, т. е.
N ^ |
(fe + y + |
]) — С |+2. |
|
Кроме того, требуется, чтобы |
каждый |
фактор варьировался |
|
не менее чем на трех |
уровнях. |
работы |
[581 появилось значи |
С момента публикации первой |
тельное число работ, посвященных теории планов второго порядка. Предложено много разных типов планов. Планы различаются между собой по числу точек; в факторном пространстве — числу опытов и т. д. Размещение точек плана влияет на статистичес кие свойства получаемой модели.
Симметричные планы второго порядка. Наиболее широкое рас пространение на практике находят симметричные планы второго порядка. Под симметричным планом второго порядка понимают план, отличающийся тем, что все нечетные моменты плана равны нулю:
N |
|
N |
xi£xjsxvg |
|
|
|
v = l , |
2, . . |
k |
i + i 4= Щ(191) |
|
2 |
xig— 0; |
2 |
|
9; t, |
/, |
||||||
g=i |
N |
g=i |
|
|
|
N |
|
|
|
|
|
|
xigxig ~ |
|
i =h /; |
xigxisxvg = |
0; |
i=fv- |
|||||
|
2 |
9; |
2 |
||||||||
|
g=i |
|
i, |
|
v = |
g=i |
|
|
|
||
а |
четные |
моменты |
/, |
1, |
2, . . . , |
К |
|
|
|||
|
|
|
|
|
|
|
|
||||
|
2 xh = |
(П = м *; |
|
2 xh xh = (?П = М з; i + /; (192) |
|||||||
|
g=i |
N |
|
|
|
g=i |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 *ig = |
(t4) = |
NU\ |
|
i, i = |
1, 2, .. ., |
k, |
|||
|
|
g=i |
|
|
|
|
|
|
|
|
188
N — общее число опытов, выполняемых по плану; g — порядко
вый |
номер |
опыта |
при |
сплошной |
нумерации |
опытов (g = 1, |
||||
2, . . |
N). |
|
|
|
|
|
|
|
|
|
Если план содержит повторяющиеся опыты, то четные моменты |
||||||||||
плана могут быть выражены иначе: |
|
|
|
|
||||||
NX2 = |
(t2) = |
£ |
NX3 = |
(t2/2) = |
2 j !'vXbx}v-, |
i =h /; |
||||
|
|
|
|
|
|
|
v = |
l |
|
|
|
|
|
|
|
|
|
|
|
|
(193) |
|
|
|
n k = |
(i4) = |
S |
h |
/ |
= |
1 .2 ........k, |
|
|
|
|
|
|
V = 1 |
|
|
|
|
|
где |
и — порядковый |
номер точки плана, |
отличающейся по усло |
виям проведения опытов от других точек плана. Если обозначить число одинаковых опытов, проводимых в каждой такой точке,
буквой |
(rv |
1), то общее число всех |
опытов |
|
|
S r „ = W . |
(194) |
|
|
»=i |
|
Величина |
п — число разных точек в |
плане. Иначе говоря, |
v — это порядковый номер серии, содержащей rv «параллельных»
опытов.
Матрица моментов, отвечающая симметричному плану второго порядка, в соответствии с уравнением (190) имеет вид
1 |
0 |
1 |
0 |
1 |
ЯгЦ, 1 |
|
0 |
XqE^ 1 |
0 |
1 |
0 |
(195) |
|
0 |
0 |
№ |
! |
|
0 |
|
1 |
0 |
! |
о |
! (Я4 — Я3) Еа-f- Я3Ektk |
|
где Ek, |
Ес2 — единичные „матрицы, |
соответственно |
порядка k и |
|||
|
|
|
|
t-A— * ( *2- 1 ) .’ |
|
|
где \kt i, |
\kk— матрицы |
соответственно размерами (&Xl) n (k x k), |
||||
составленные |
из |
( + 1); |
* — знак |
транспонирования. |
||
Матрица |
(195) |
в планировании |
экспериментов |
обычно назы |
вается информационной матрицей. Детерминант этой матрицы (точнее, обратная величина детерминанта) может служить мерой точности оценок коэффициентов регрессии:
detN~ l М = [Я4 |
+ |
k(k-i) |
(Я3 — А|)] (Я4— Я з ) * - ^ 2 . (196) |
189
Если план невырожденный, то
det N - m Ф О,
и при этом соблюдаются неравенства [46]
о ■‘С ^ 3 |
^ 4 ^== ^ 2 “С |
1 ; |
О <С[ Я4 — |
Я3 -|- /бЯз — |
( 1 9 7 ) |
kX^. |
||
Симметричные планы отличаются симметричным расположе |
||
нием точек в факторном пространстве. К симметричным планам |
||
относятся ортогональные, |
ротатабельные, D -оптимальные, сим- |
плексно-суммируемые планы и планы типа вращаемого эволю ционного планирования (РОВОП), используемые в эволюционном планировании.
Во многих случаях, особенно |
при проведении экспериментов |
в производственных условиях, |
трудно выдержать одинаковые |
условия выполнения всей совокупности намеченных опытов. В этой связи большой интерес представляют планы второго порядка с разбиением опытов на отдельные блоки, в которых однородность экспериментальных условий должна быть ограни чена лишь рамками блока. Группирование опытов в блоки позво ляет повысить точность получаемых оценок параметров (коэф фициентов) модели. При разбиении всех опытов на ортогональные блоки оценки параметров модели вычисляются независимо от различий между блоками в отношении однородности. К тому же и сами расчетные формулы для определения оценок не изменяются, небольшие коррективы приходится вносить лишь в формулы проверки гипотез.
Приведем (без вывода формул) общую последовательность регрессионного анализа результатов эксперимента, выполняемого по симметричному плану без разбиения на блоки и с разбиением
.на ортогональные блоки.
Регрессионный анализ без разбиения опытов на блоки. 1) Точеч ные оценки параметров модели (коэффициентов регрессии) при обычных предположениях о равноточности и некоррелирован ности отдельных наблюдений определяются по формулам:
k
Ь0 - аЫ~г {Оу) — ьы~х s т - , |
|
||
|
|
i=1 |
|
bu — — bN 1 (Оу) - f cN 1 (iiy) — dN 1 £ (iiy)\ |
( 1 9 8 ) |
||
bi = |
|
i—1 |
|
Я^-W -1 (((/); |
i Ф 0; |
|
|
btj = |
я - w 1 (ijy)\ |
i =j= j . |
|
190