Файл: Кацев, П. Г. Статистические методы исследования режущего инструмента.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