Теория эксперимента. Учебное пособие
.pdf
SS |
Q1 Q2 . |
(4.32) |
8 Рассчитывают оценку дисперсии S2 x :
S2x SSx |
(4.33) |
m 1. |
9 Рассчитывают оценку дисперсии S2 :
S2 |
|
|
SS |
|
|
|
m n 1 . |
(4.34) |
|||||
|
|
|
||||
10 Оценивается влияние фактора Х по зависимости (4.18) - (4.25).
Результаты расчета, как правило, представляются в виде следующей таблицы дисперсионного анализа.
Таблица 4.2 – Однофакторный дисперсионный анализ с равным числом повторений опытов
Компоненты |
Число |
Сумма |
|
|
Средний |
Математическ |
|||||||||
дисперсии |
степеней |
квадратов |
|
|
квадрат |
ое ожидание |
|||||||||
|
свободы |
|
|
|
(оценка |
среднего |
|||||||||
|
|
|
дисперсии) |
квадрата |
|||||||||||
Междууровн |
|
SSx Q2 Q3 |
|
|
2 |
|
|
|
SSx |
2 |
2 |
||||
евая |
m-1 |
S |
|
x |
|
|
|
|
n x |
|
|||||
|
m 1 |
||||||||||||||
|
|
|
|
|
|
|
|
|
|
||||||
Внутриуров |
|
|
S2 |
|
|
|
|
SS |
|
|
2 |
||||
m(n-1) |
SS Q1 Q2 |
m n 1 |
|
||||||||||||
невая |
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|||||
Общая |
mn-1 |
SSобщ. Q1 Q3 |
S |
2 |
|
|
|
SSобщ |
– |
||||||
(полная) |
|
y |
|
|
|
|
|
|
|||||||
|
|
mn 1 |
|||||||||||||
|
|
|
|
|
|
|
|
|
|||||||
В случае, когда на различных уровнях фактора проводится разное число параллельных опытов, можно, ориентируясь на уровень с меньшим числом,
отбросить лишние наблюдения в остальных уровнях. Однако такое
171
отбрасывание нежелательно, так как существенно снизит точность проводимого анализа. Тем более что новая схема вычислений лишь не многим будет отличаться от случая равных столбцов. Пусть на уровне хi проведено ni
параллельных наблюдений. Общее число всех наблюдений будет равно:
|
m |
|
|
|
N ni . |
|
|
(4.35) |
|
|
i 1 |
|
|
|
Определим: |
|
|
|
|
1 Итоги по столбцам – суммы наблюдений |
по каждому |
уровню |
||
фактора хi: |
|
|
|
|
|
n |
|
|
|
Yi yji . |
|
|
(4.36) |
|
|
j 1 |
|
|
|
2 Суммы квадратов всех наблюдений: |
|
|
||
n |
m |
|
|
|
Q1 |
yji |
2 . |
|
(4.37) |
j 1 |
i 1 |
|
|
|
3 Сумму квадратов итогов по |
столбцам, |
деленных на |
число |
|
наблюдений в соответствующем столбце: |
|
|
|
|
m |
Y2 |
|
Q2 |
i |
. |
|
||
|
n |
|
i 1 |
i |
|
4 Квадрат общего итога, деленный на число всех наблюдений:
Q3 1 m Yi 2.
i 1
N
(4.38)
(4.39)
172
Дальнейшие расчеты проводятся по формулам (4.30) - (4.34). Если оценки дисперсий S2 x :и S2 :значимо отличаются друг от друга, то дисперсию фактора Х можно вычислить по формуле:
x2 |
m 1 N |
Sx2 S 2 . |
(4.40) |
m |
|||
|
N 2 ni |
|
|
i 1
При этом соотношение для чисел степеней свободы следующее:
fSобщ2 . N 1 N m m 1 fS2. |
fSx2. . |
(4.41) |
Задача 4.1.
Пороховой завод изготовил заряды к артиллерийским боеприпасам на четырех однотипных поточных линиях. С каждой поточной линии случайным образом было выбрано и отстрелено в одинаковых условиях по пять боеприпасов. Результаты замеров начальных скоростей полета снарядов V0,
м/с, для каждой группы приведены в таблице 4.3.
Таблица 4.3 – Результаты испытаний
Номер опыта |
Уровни фактора Х (поточные линии), i=1,m |
|||
j=1,n |
x1 |
x2 |
x3 |
x4 |
1 |
600 |
570 |
690 |
450 |
2 |
420 |
450 |
570 |
510 |
3 |
510 |
630 |
600 |
450 |
4 |
435 |
450 |
570 |
510 |
5 |
495 |
450 |
600 |
540 |
|
|
|
|
|
Стоит задача: выяснить, существенно ли влияние данных поточных
линий на величину начальной скорости полета снарядов V0.
Решение.
1 Проверяем предпосылки ДА.
173
Пусть из априорных данных известно, что результаты наблюдений являются независимыми случайными величинами, имеющие нормальный закон распределения; опыты, проведенные нами, равноточны, а случайные ошибки наблюдений подчиняются нормальному закону. Однородность выборочных дисперсий проверим по критерию Кохрена.
– определяем средние по столбцам:
y1 |
|
1 yj1 |
600 420 510 435 495 492 ; |
||||
|
|
|
|
n |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
n |
j 1 |
5 |
|
|
y2 510; |
y3 606; |
y4 492; |
– определяем оценки дисперсий для каждого столбца:
|
|
|
|
|
|
n |
|
|
2 |
|
|
|
|
|
|
|
|
yj |
y |
|
|
||
|
|
|
|
Sy2 |
|
j 1 |
|
|
|
||
|
|
|
|
n 1 |
|
|
|||||
|
|
|
|
1 |
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
600 492 2 420 492 2 510 492 2 |
|||||||||
|
|
|
|
|
|
|
|
|
|
|
5107; |
5 1 |
|
|
|
|
|
|
|
||||
|
435 492 2 495 492 2 |
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
7200; |
2 |
|
|
2 |
1620; |
||
|
|
|
Sy2 |
Sy3 2430; |
|
Sy4 |
|||||
– определяем расчетное значение критерия Кохрена:
|
2 |
max |
7200 |
|
|
Syi |
|
||||
Gрасч |
|
|
|
|
0.4402 ; |
m |
|
|
|||
Sy2i |
5107 7200 2430 1620 |
||||
i1
–определяем табличное значение G-критерия по таблице приложения 4:
174
Gтабл. 0,05; |
n 4; |
f 5 1 4 0,6287. |
Так как Gрасч<Gтабл., то оценки дисперсий однородны и можно приступить к проведению ДА.
2 Подсчитываем итоги по столбцам:
n
Yi yji , j 1
Y1 600 420 510 435 495 2460;
Y2 2550;
Y3 3030;
Y4 2460.
3 Вычисляем сумму квадратов всех наблюдений:
n |
m |
|
Q1 |
yji |
2 ; |
j 1 |
i 1 |
|
Q1 6002 4202 ... 5402 5622750.
4 Вычисляем сумму квадратов итогов по столбцам, деленную на число наблюдений в столбце:
Q2 1 m Yi 2 ; n i 1
175
Q |
1 |
24602 25502 |
30302 |
24602 5557320. |
|||||
5 |
|||||||||
2 |
|
|
|
|
|
|
|||
5 Рассчитаем квадрат общего итога, деленный на число всех |
|||||||||
наблюдений: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 m |
|
2 |
||
|
|
|
|
Q3 |
|
Yi |
; |
||
|
|
|
|
||||||
|
|
|
|
|
mn i 1 |
|
|
||
Q |
|
|
1 |
2460 2550 3030 2460 2 5512500. |
|||||
|
|
|
|||||||
3 |
|
4 5 |
|
|
|
|
|||
6Определяем сумму квадратов для столбца:
SSx Q2 Q3 ;
SSx 5557320 5512500 44820.
7Определяем общую сумму квадратов:
SSобщ. Q1 Q3 ;
SSобщ. 5622750 5512500 110250.
8 Определяем остаточную сумму квадратов для оценки ошибки эксперимента:
SS Q1 Q2 ;
SS 5622750 5557320 65430.
9 Рассчитывают оценку дисперсии S2 x :
176
S2x SSx ; m 1
S2x 44820 14940. 4 1
10 Рассчитывают оценку дисперсии S2 :
S |
2 |
|
SS |
; |
|
m n 1 |
S |
2 |
|
|
65430 |
|
4089. |
|
4 5 1 |
|
11 Вычисляем расчетное значение F - критерия:
Fрасч |
|
S2x |
; |
|
S2 |
||||
|
|
|
14940
Fрасч 4089 3,6537.
12 Определяем |
критическое |
табличное значение F - критерия |
(приложение Б): |
|
|
Fтабл. 0,05; |
f1 m 1 3; |
f2 m n 1 16 3.2389. |
13 Сравним Fрасч и Fтабл:
177
Fрасч 3,6537 Fтабл 3,2389.
Таким образом, влияние поточных линий на величину начальной скорости полета снаряда V0 следует признать значимым.
14 Определяем величину влияния номера поточной линии:
x2 |
|
Sx2 S2 |
|
14940 4089 |
2170. |
n |
|
||||
|
|
5 |
|
||
15 Производим сравнение влияния номеров поточных линий по критерию Дункана. С этой целью:
- вычисляем нормированную ошибку среднего:
S |
|
|
S2 |
|
4089 |
28.60; |
|
n |
5 |
||||
|
||||||
|
y |
|
|
|||
|
|
|
|
|
||
- располагаем значения средних по столбцам (смотреть таблицу 2.3) в
порядке возрастания их величин:
y1 y4 y2 y3 ;
492 492 510 606
- выписываем из таблицы Дункана (приложение В)
nD fS2 m n 1 16; |
0.05; |
|
p 2.3...,m |
все m |
1 значения значимых рангов |
178
р |
2 |
3 |
4 |
ранги |
3,00 |
3,15 |
; |
3,23 |
- вычисляем (m-1) наименьшие значимые ранги (НЗР), умножая ранги на нормированную ошибку среднего:
р |
2 |
3 |
4 |
|||
НЗР S |
|
|
|
90.09 |
; |
|
|
|
85.80 |
92.38 |
|||
|
||||||
y |
||||||
|
||||||
- проверяем значимость различия между средними по столбцам для всех m(m-1)/2 парами:
y3 |
y1 606 492 114 92,38 |
различия значимо; |
|||||||
|
3 |
|
|
|
4 |
606 |
492 114 92,09 |
различия значимо; |
|
y |
|
y |
|||||||
|
3 |
|
|
|
2 |
606 |
510 96 85.80 |
различия значимо; |
|
y |
|
y |
|||||||
|
2 |
|
|
1 |
510 492 18 85.80 |
различия незначимо; |
|||
y |
y |
||||||||
|
2 |
|
|
4 510 |
492 18 85.80 |
различия незначимо; |
|||
y |
y |
||||||||
|
4 |
|
|
1 |
492 492 0 85.80 различия незначимо. |
||||
y |
y |
||||||||
Таким образом, для выпуска однородной продукции усовершенствование целесообразно начать с третьей поточной линии.
4.4Двухфакторный дисперсионный анализ
Дисперсионный анализ особенно эффективен при одновременном изучении нескольких факторов. В этом случае при постановке активного эксперимента одновременно варьируют все факторы, и каждое наблюдение выходного параметра Y используют для их одновременной оценки. При этом немаловажен и тот факт, что в ряде случаев можно не делать параллельных
179
наблюдений, ограничиваясь лишь одним наблюдением для каждого сочетания уровней изучаемых факторов.
Пусть требуется изучить одновременное влияние на процесс двух качественных факторов ХА и ХВ. Фактор ХА исследуется на уровнях ХА1, ХА2, …,
ХАi, …, ХАК, фактор ХВ – на уровнях ХВ1, ХВ2, …, ХВj, …, ХВm. Допустим, что при каждом сочетании уровней исследуемых факторов было проведено по n
параллельных наблюдений (равное количество взято для упрощения вкладок).
Матрицу наблюдений в рассматриваемом случае можно представить в виде таблице 4.4.
Общее число наблюдений равно:
N n k m.
Таблица 4.4 – Исходные данные для двухфакторного ДА с равным числом повторений опытов
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
XA |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
ХВ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
y j |
|||||||||||||
|
|
|
xA1 |
|
|
xA2 |
… |
|
|
|
|
xAi |
… |
|
|
|
xAK |
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
y111, y112, |
y211, y212, |
… |
|
yi11, yi12, |
… |
yK11, yK12, |
|
|
|
|
|
|
|
|||||||||||||||||||||||
|
xB1 |
…, y11n |
…, y21n |
|
…, yi1n |
…, yK1n |
|
|
y 1 |
||||||||||||||||||||||||||||||
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
… |
|
|
|
|
|
|
|
|
… |
|
|
|
|
|
|
|
|
K1 |
|
|
|
|
|
|
|
|
|
|
|
|
y |
|
11 |
|
|
y |
21 |
|
|
|
|
y |
|
i1 |
|
|
|
y |
|
|
|
|
|
|
|
||||||||||
|
|
|
y121, y122, |
y221, y222, |
… |
|
yi21, yi22, |
… |
yK21, yK22, |
|
|
|
|
|
|
|
|||||||||||||||||||||||
|
xB2 |
…, y12n |
…, y22n |
|
…, yi2n |
…, yK2n |
|
|
y 2 |
||||||||||||||||||||||||||||||
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
… |
|
|
|
|
|
|
|
|
… |
|
|
|
|
|
K2 |
|
|
||||||||||||
|
|
|
|
|
y |
|
12 |
|
|
y |
22 |
|
|
|
y |
|
i2 |
|
|
y |
|
|
|
|
|
|
|
|
|||||||||||
|
… |
|
|
… |
|
|
… |
… |
|
|
|
|
… |
… |
|
|
|
|
… |
|
|
… |
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
y1j1, y1j2, |
y2j1, y2j2, |
… |
|
yij1, yij2, |
… |
yKj1, yKj2, |
|
|
|
|
|
|
|
|||||||||||||||||||||||
|
xBj |
…, y1jn |
…, y2jn |
|
…, yijn |
…, yKjn |
|
|
y j |
||||||||||||||||||||||||||||||
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
y1j |
|
|
y2 j |
… |
|
|
|
|
y |
ij |
… |
|
|
|
y |
Kj |
|
|
|
|
|
|
|
|||||||||||
|
… |
|
|
… |
|
|
… |
… |
|
|
|
|
… |
… |
|
|
|
|
… |
|
|
… |
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
y1m1, y1m2, |
y2m1, y2m2, |
… |
|
yim1, yim2, |
… |
yKm1,yKm2, |
|
|
|
|
|
|
|
|||||||||||||||||||||||
|
xBm |
…, y1mn |
…, y2mn |
|
…, yimn |
…, yKmn |
|
|
y m |
||||||||||||||||||||||||||||||
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
|
y1m |
|
|
y2m |
… |
|
|
yim |
… |
|
|
|
|
|
|
Km |
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
y |
|
|
|
|
|
|
|
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
… |
|
|
|
|
|
|
… |
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
y |
|
|
|
y |
|
|
|
y |
|
|
|
|
|
y |
|
|
|
|
y |
|
|
|
|
|
|
|
||||||||||||
|
i |
|
|
1 |
|
|
2 |
|
|
|
|
i |
|
|
|
K |
|
|
|
|
|
y |
|
||||||||||||||||
180
