Файл: Кацев, П. Г. Статистические методы исследования режущего инструмента.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 21.10.2024

Просмотров: 146

Скачиваний: 1

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

шера и Бартлета. Использование F-критерия Фишера при числе дисперсий более двух неэффективно, так как при этом в оценке участвуют только наибольшая и наименьшая дисперсии. Крите­ рий Кохрена пригоден для случаев, когда число повторных опы­ тов во всех точках плана одинаково. Из всех дисперсий sj нахо­ дится наибольшая s?, таХ) которая делится на сумму всех диспер­ сий по точкам. Критерий Кохрена — это отношение максимальной дисперсии к сумме всех дисперсий. В нашем случае, пользуясь табл. 45, найдем

55,44

=

0, 22.

(148)

249,0

 

 

 

По соответствующей таблице в работе [11]

находим для /отах =

= 2, /З..ам = N степеней свободы и

 

уровня значимости 5%. Кри­

тическое значение GKP = 0,61. Гипотеза об однородности диспер­ сий принимается, если, как в нашем случае, экспериментальное значение критерия Кохрена не превышает табличного значения.

Регрессионный анализ. На основе метода наименьших квадра­ тов найдено уравнение связи или математическая модель. Метод наименьших квадратов использован как вычислительный прием. Теперь предстоит выполнить статистические оценки полученной модели.

Обычный регрессионный анализ основан на следующих пред­ посылках.

1. Результаты наблюдений у ъ у 2, . . ., уп параметра опти­ мизации в и точках факторного пространства представляют собой

независимые, нормально распределенные случайные величины. О том, что означает это условие и как оно проверяется, подробно сказано в гл. II.

2.Дисперсия величины у не зависит от абсолютной величины у

изначений факторов, т. е. дисперсии в разных точках плана одинаковы. Проверка выполнения этого условия была показана выше.

3.Значения факторов суть неслучайные величины. Практи­

чески это означает, что независимые переменные л^, а 2 , .. ., хк

измеряются с пренебрежимо малой ошибкой по сравнению с ошиб­ кой воспроизводимости для этих факторов. Нарушение этого условия приводит к трудностям при реализации матрицы плани­ рования и поэтому легко обнаруживается в процессе выполнения эксперимента.

Проверка значимости коэффициентов модели. Проверка значи­ мости каждого коэффициента проводится независимо. Для этого можно использовать проверку по ^-критерию Стыодента. При использовании полного факторного эксперимента или регуляр­ ных реплик (о них ниже) доверительные интервалы для всех коэф­ фициентов равны друг другу. Прежде всего находим дисперсию коэффициента регрессии s2 {!?,•}. При равномерном дублировании

142


опытов по точкам с числом повторных опытов г она определяется

по формуле

* \ ь , !

l i M

(149)

nr

сfE = п (/•— 1) степенями свободы.

Внашем случае

s2{M = - % r = 1-296; s{M = V T 2 9 6 = l,1 3 .

Из формулы видно, что дисперсии всех коэффициентов равны друг другу, так как они зависят только от ошибкиопыта и числа опытов. Теперь рассчитаем значения ^-критерия по формуле

 

 

,

_

Ibi |

 

 

 

(150)

 

 

 

'■

s{M

 

 

 

 

 

 

 

 

Для нашей задачи получаем:

 

 

 

 

h :

34,47

tx-. "

5,59

4,95;

t2-

2,01

=1,78;

1,13 = 30,49;

1,13

1 1,13

^ з = т т г ^ 6’76;

 

0,65

0,57;

/1Я=

0,65

: 0,57;

112 -

1,13

1,13

 

t23 --- 0,39

= 0,34;

‘123

0,5

= 0,44.

 

 

1,13

 

 

 

1,13

 

 

 

Критическое значение tKP находится по таблице работы [11]- при п (г — 1) = 16 степенях свободы и заданном уровне значи­ мости а — 5%. В нашем случае 4 р — 1,74. Если tL> fKP, то гипотеза отвергается и коэффициент bt признается значимым.

Впротивном случае bt считается статистически незначимым, т. е.

=0. В нашем случае такими незначимыми коэффициентами являются b12, b13, bад, Ь123. Теперь можно построить довери­

тельный интервал длиной 2А bh где

Д6, = V {М = 1,74-1,296 = 2,26.

(151)

Коэффициент значим, если его абсолютная величина больше половины длины доверительного интервала. Ортогональное пла­ нирование позволяет определять доверительные границы для каждого из коэффициентов регрессии в отдельности, если какойлибо из коэффициентов окажется незначимым, он может быть отброшен без пересчета всех остальных. После этого математи- . ческая модель объекта составляется в виде уравнения связи вы­

ходного параметра у и переменных х(, включающего только зна­

чимые коэффициенты. Для нашего примера получаем

у = 34,47 + 5,59*! + 2,01х3 — 7,64х3.

(152)

143


Чтобы получить модель в натуральных переменных, надо в уравнение (152) подставить выражения х( из формулы преобра­

зования (141).

Получаем

 

 

 

у = 22,6

l,4otj

0,34у2 — 25,5ср3.

(153)

Статистическая незначимость коэффициента /?£ может быть

обусловлена следующими причинами:

 

 

 

1)

уровень

базового режима х[0 близок

к точке

частного эк-

 

 

 

 

 

 

 

А <%■

 

стремума

по

переменной

xit т.

е.

=

- s O

;

 

 

 

 

 

 

 

dxi

 

2)

шаг

варьирования

Дх£

выбран

малым;

 

3)

данная

переменная

(произведение переменных) не имеет

функциональной связи с выходным параметром у;

4) велика ошибка эксперимента вследствие наличия неуправ­ ляемых и неконтролируемых переменных.

Если имеет место первая или третья причина, значение фак­ тора стабилизируется на определенном уровне, во втором случае увеличивают интервал варьирования. Если имеет место четвертая

причина, следует

принять меры

к уменьшению ошибки экспери­

мента.

 

 

В нашем случае незначимость коэффициентов, оценивающих

взаимодействие

(Ь12\ й13; й23;

вызвана, вероятно, отсут­

ствием этих эффектов. В качестве величины, характеризующей вклад коэффициентов регрессии в уравнение, иногда используют

множественный

коэффициент корреляции

R:

 

R =

(yv — tjv) 2

1 Л ___

 

1 — °

----_

0,95,

 

оУ)2

V

790,44

 

 

 

 

 

где у = 34,47; £

(у0У) 2

= 39

(см. табл.

45);

£ (Уи ~ У) 2 =

= 790,44 (табл.

46).

 

 

 

 

Таким образом, уравнение регрессии практически полностью описывает результаты эксперимента. Чем больше по абсолютному значению величина коэффициента регрессии bt, тем сильнее

влияние его на критерий оптимизации в заданном интервале варь­ ирования факторов. Если Ь{ > 0, то увеличение х£ вызывает

увеличение критерия оптимизации. В противном случае, т. е. при Ь£ < 0, увеличение х{ приводит к уменьшению критерия опти­

мизации. Подобные рассуждения справедливы для модели линейной по факторам. Если же некоторый коэффициент при сме­ шанном произведении факторов, например, Ьц является значимым,

то это свидетельствует о том, что Действие одного из этих фак­ торов, скажем xt, зависит от уровня, на котором находится другой

144


фактор, Xj. Вклад слагаемого bijxixj в величину критерия опти­ мизации при bti > 0 будет положительным, если оба фактора

находятся на верхних или нижних уровнях. И наоборот, вклад этого слагаемого будет отрицательным, если факторы находятся

на разных уровнях. При

b(j <

0 картина будет обратной.

 

Проверка

адекватности

мо­

 

 

 

 

 

 

дели,

После

вычисления коэф-

Таблица 46

 

 

 

 

фициентов

модели

следует,

расчет V

(yv—y) 2

 

 

прежде всего,

проверить ее при-

 

^

 

 

 

 

годность или

адекватность

мо­

 

 

 

 

 

дели.

Чтобы

проверить

гипо­

точки

 

yv

1 у 0 - у )

(у 0

и ) г

тезу об адекватности представ­

V

 

 

 

 

 

 

 

 

 

 

 

ления

результатов

эксперимен­

1

36,90

2,44

 

5,95

та найденным уравнениям связи,

 

достаточно оценить

отклонение,

2

44,07

9,61

92,35

3

34,83

0,37

 

0,14

предсказанное

уравнением

ре­

4

52,62

18,16

329,78

грессии

выходной

величины у

5

18,11

—16,35

267,32

6

30,72

—3,74

13,99

от результатов

эксперимента у

7

25,62

—8,84

78,15

вПразличных точках факторного

8

32,80

—1,66

 

2,76

пространства.

 

 

 

 

 

 

 

 

 

 

 

 

Рассеяние результатов экспе­

У =

34,47

2 =

790,44

римента

относительно

уравне­

ния связи,

аппроксимирующего

 

 

 

 

 

 

искомую

функциональную

за­

 

 

 

 

 

дис­

висимость,

можно

охарактеризовать с помощью остаточной

персии

или дисперсии

адекватности а |д,

оценка

которой

нахо-

дитаГ^по формуле,

которая справедлива лишь при равном числе

дублирующих

 

опытов

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ад ■

 

~

1~

Ъ ( У » ~ Ъ \

 

 

 

(154)

 

 

 

 

 

 

 

 

п

т

0=1

 

 

 

 

 

т — число

 

членов

аппроксимирующего

полинома

(включая

свободный

член).

 

 

 

 

 

 

 

 

 

 

Дисперсия адекватности определяется с числом степеней сво­ боды *

/ад = П — т-

Проверка адекватности состоит в выяснении соотношения между дисперсией адекватности s|fl и дисперсией воспроизводи­ мости s2 {г/}. Проверка гипотезы об адекватности модели про­ водится с использованием F-критерия Фишера.

* Числом степеней свободы в статистике называется разность между числом опытов и числом коэффициентов (констант), которые уже вычислены по резуль­ татам этих опытов независимо друг от друга.

Ю П. Г. Кацев

145


Критерий

Фишера позволяет проверить нуль-гипотезу о ра­

венстве двух

генеральных дисперсий ст!д и а 2

\у\ в

том

случае,

если F-критерий формируется как отношение

F =

s3^ }

' ^ слн

вычисленное значение критерия меньше критического FKP, опре­ деляемого по таблице работы [11] для соответствующих степе­ ней свободы

f a n =

> 1 — и171 f E

= n ( r

1)

при заданном уровне

значимости — а

%, то

нуль-гипотеза при­

нимается. В противном случае гипотеза отвергается, и описание признается неадекватным объекту. Если выборочная дисперсия неадекватности s;^ не превосходит оценки дисперсии воспроизво­

димости

s2 {г/}, тогда F-отношение будет меньше

(или

равно)

единице

и неравенство F <

FKP выполняется

для

любого

числа

степеней свободы /ад и / £, т.

е. гипотеза а1д с

а2 {*/) не противо­

речит выборочным данным и математическая модель адекватно представляет объект.

Проверка адекватности возможна лишь при /ад > 0. Если число точек в плане равно числу оцениваемых коэффициентов регрессии уравнения связи (я = т), то, следовательно, не остается

степеней свободы (/ад = 0) для проверки нуль-гипотезы об аде­ кватности представления экспериментальных данных выбранной формой аппроксимирующего полинома. В этом случае можно

воспользоваться тем обстоятельством,

что

свободный член есть

k

 

 

совместная оценка, т. е. Ь0 —>р 0 + ^

Рп .

Если имеются повтор-

г=1

 

 

ные опыты в центре планирования, т. е. в

точке с координатами

(0, 0, 0,

. . ., 0), то средний отклик дает

несмещенную оценку

Уо—> Ро-

Поэтому, если разность Ь0 у 0 окажется статисти­

чески значимой, то это будет указывать на неадекватность линей­ ной модели, а также модели, содержащей взаимодействия, т. е. на то, что хотя бы часть коэффициентов |3l7 не равна нулю. Если некоторые коэффициенты регрессии оказались незначимыми или ими можно пренебречь ввиду их малости или требуется доказать правомерность линейной аппроксимации в заданном интервале варьирования, то число членов проверяемого уравнения в этом случае будет, как правило, меньше числа точек в плане и одна или несколько степеней свободы останутся для проверки гипотезы адекватности.

Проведем оценку адекватности в нашем примере, пользуясь

данными расчетной табл.

45:

 

 

s ap,

п

т

1 >(У»~У«)2 = # ^

= 29,2.

 

и = 1

8 — 4

 

 

 

 

 

 

Так как в нашем случае я|д < s2 \у\ (29,2 < 31,12), то аде­

кватность модели очевидна без расчета критерия Фишера. Если бы

146