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