Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Теория эксперимента. Учебное пособие

.pdf
Скачиваний:
0
Добавлен:
15.08.2026
Размер:
1 Мб
Скачать

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

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]