Файл: Химмельблау Д. Анализ процессов статистическими методами.pdf

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

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

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

Добавлен: 09.04.2024

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

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

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

502

 

 

 

Глава

8

 

 

 

 

 

 

 

 

 

Таблица

П.8.1.16

 

Дисперсионный

анализ некодированных

данны х

 

 

 

 

 

 

Сумма

Число степе­

Средний

Отношение

Источник рассеяния

 

 

дисперсий

квадратов

ней свободы

квадрат

 

 

 

 

 

 

 

 

 

«2/31,3

Влияни е

исключения Ьі

 

504,1

1

504,1

 

16,1

Влияние

исключения Ьд

 

3125

1

3125

 

100

(после

исключения Ьі)

 

93,9

3

31,3

 

 

Отклонения от

эмпири­

 

 

 

ческой

линии

регрес­

 

 

 

 

 

 

сии

 

 

 

 

 

 

 

 

Общий

 

 

 

3723,0

5

 

 

 

Замечание: F0 ,95

(1, 3) =

10,13.

 

 

 

 

Д ля того чтобы решить, является ли некоторая модель подхо­

дящей,

можно

испытать

все виды моделей и

посмотреть, какая

из них дает наименьшие дисперсии для отклонений от эмпириче­

ской линии регрессии. Однако лучше, если это возмояшо,

выпол­

нить повторные

эксперименты

и таким образом получить

оценку

экспериментальной ошибки,

которая не зависит от математиче­

ской модели. Как объяснялось в гл. 4, можно использовать

F-кри-

терий, чтобы установить, является ли подгонка относительно оцен­ ки линии регрессии значимо отличной от экспериментальной ошибки. Всегда выгодно иметь достаточное число повторений опыта, чтобы обеспечить несколько степеней свободы для оцени­ вания ошибки, уменьшая тем самым средний квадрат ошибки до возможно малого значения, определяемого соображениями экономичности/

Пример 8.1.2. Одномерный эксперимент с повторением С целью расширения эксперимента из примера 8.1.1 были про­

ведены повторные

опыты

(ряд полностью завершенных опытов):

Температура,

Напряжение

Температура,

Напряжение

°С

сдвига, Н/м2

С

сдвига, Н/м2

60

(

14

 

Г 26

I

15*

120

{

26

 

I

16

 

I

28*

 

 

1

6

 

(

44

80

f

140

і

45*

17'

 

I

18

 

I

46

20

10020*

21

Здесь символом* обозначены старые значения.


Стратегия

эффективного

экспериментирования

5 03

Модель первого порядка Y =

ß 0 - f ßi x + e можно

подогнать

к кодированным данным (как в предыдущем примере),

вычисляя

следующие матрицы:

 

 

 

 

 

 

матрица а

 

 

Г1.500

-101

7,500-Ю1 !,

 

 

t-7,500-101

4,950.102 J

 

обратная матрица а - 1 ( м а т р и ц а с)

 

Г

2,270

-Ю-1

—4,166-Ю-*"!,

 

L—

4,166

-Ю-2

8,333-Ю-3 J

 

матрица (вектор) g 11,470-10* 1,153-103],

коэффициенты регрессии

Ъ0 = -7,616, by = 3,483.

Результаты соответствующего дисперсионного анализа приве­ дены в табл. П.8.1.2а.

 

 

 

 

 

 

Таблица

П.8.1.2а

 

 

 

 

Сумма

Число степе­

Средний

Отношение

Источник

рассеяния

 

дисперсий

квадратов

ней свободы

квадрат

 

 

 

 

 

 

 

 

 

«V0.933

Влияние

исключения Ь\

1456

1

1456

 

1560*

Влияние

исключения bo

1440

1

1440

 

1500*

(после

исключения by)

309,0

3

103,0

 

110*

Отклонения

от

эмпири­

 

ческой

линии

р е г р е с ­

 

 

 

 

 

сии

 

внутри на­

9,33

10

0,933

 

 

Отклонения

 

 

боров

 

 

 

 

 

 

 

 

Общий

 

 

 

3215

15

 

 

 

Замечания: F 0 ,95

(1, 10) = 4,96; FQ>96

(3, 10) = 3 , 7 1 .

 

 

 

Символом * отмечены значимые величины.

Простой расчет дает результат, аналогичный результату в при­ мере 8.1.1, а именно остаточную сумму квадратов

,309,0 + 9,33

318,33


504

Глава 8

которая приводит к отношению

318,33 = 24,4

13

вотличие от 31,3 в примере 8.1.1.

^-критерий показывает, что отклонение от линии регрессии первого порядка значительно больше экспериментальной ошибки и что модель можно существенно улучшить. Так как линейная модель первого порядка оказалась неудовлетворительной, можно испытать другую модель первого порядка, полученную, возможно, с помощью некоторого преобразования. Аддитивная ненаблюдае­

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

ане к первоначальному отклику. Если такие модели окажутся неэффективными, можно испытать модель второго порядка. При этом лучше руководствоваться соображениями простоты модели,

ане выбирать непродуманно сложную модель.

Для

модели

 

 

 

 

 

 

 

 

Y =

ß 0

+

ßi* +

ß2 tf2

+ г

матрицы

имеют вид

 

 

 

 

 

 

 

 

 

 

матрица

а

 

 

1,500-101

 

7,500 - 10 1

 

4,950 - Ю 2 -

 

7,500 - 10 1

 

4,950 - 10 2

 

3,675 - 10 3

 

|_4,950-102

 

3,675 - 10 3

 

2,900 - 10 4 J

 

обратная матрица а - 1

(матрица с)

 

7 , 0 5 0 - Ю - 1

— 2,946 . Ю - 1

 

2 , 5 2 9 - Ю - 3 ^

 

-2,946-10"1

 

1 , 5 7 1 - К Г 1

— 1 , 4 8 8 - « Г 2

 

2 , 5 2 9 - Ю - 2

— 1 , 4 8 8 - Ю - 2

 

1 , 4 8 8 - « Г 8 J

 

матрица

(вектор)

g

 

[1,470-102

1,153.10» 9,458-10»],

 

коэффициенты регрессии

 

 

 

Ь0

=

3,210,

 

 

 

 

h

=

- 2 , 8 8 5 ,

 

 

 

 

Ьг

=

0,636.

 

Результаты дисперсионного анализа приведены в табл. П.8.1.26. Модель второго порядка также оказалась неадекватной. Если из модели исключать Ь2 и затем bit то первые две суммы квадратов

в дисперсионном анализе равнялись бы соответственно 1676 и 53. Такое различие имеет место вследствие того, что независимые переменные я и ж2 не являются ортогональными, однако вывод о неадекватности модели остается в силе. В качестве модели можно испытать полином третьего порядка

Y = ß 0 + ß,si + ß 2 . z 2

ß 3 z 3

e.


 

 

Стратегия

эффективного

экспериментирования

 

505

 

 

 

 

 

 

 

Таблица

П.8.1.26

 

 

 

 

 

Сумма

Число степе­

Средний

Отношение

Источник рассеяния

 

дисперсий

квадратов

ней свободы

квадрат

 

 

 

 

 

 

 

 

 

 

S2/0.933

Влияние

исключения by

1456

1

1456

 

1560*

Влияни е исключения Ъ2

273

1

273

 

292*

(после

исключения

by)

 

1

 

 

1500*

Влияние

исключения

Ь0

1440

1440

 

(после

исключени я by

 

 

 

 

 

и Ъ2)

 

от

эмпири­

36

2

18,2

 

19,5*

Отклонения

 

ческой

линии

регрес­

 

 

 

 

 

сии

 

внутри

на­

9,33

10

0,933

 

 

Отклонения

 

 

боров

 

 

 

 

 

 

 

 

 

Общий

 

 

 

 

3215

15

 

 

 

Замечания:

F0,95

(1, 10) =

4,96; F 0 M (2,

10) = 4,10.

 

 

 

Символом * отмечены значимые величины.

Возможно, что линеаризованная форма некоторой нелинейной (по коэффициентам) модели типа

In Y = а + ßz + e

лучше описывала бы эти данные и привела бы к незначимым, соглас­ но F-критерию, отклонениям от эмпирической линии регрессии. Ради экономии места соответствующие расчеты здесь не приводят­ ся, но общий метод теперь должен быть ясен.

8-1-2. Двумерные модели и планы экспериментов

Теперь рассмотрим модели, содержащие две геометрические координаты и одну зависимую переменную. Метод поверхности отклика можно сделать наиболее эффективным, если использовать ортогональные планы. Когда оценки коэффициентов являются орто­ гональными 1 ) , при данном объеме экспериментов можно получить максимальный объем информации.

1 ) Термин «ортогональные оценки» требует уточнения. Оценки пара ­ метров являютс я линейными функциями наблюдаемых значений зависимой переменной У г . Каждой линейной функции можно поставить в соответствие вектор коэффициентов, входящих в данную функцию. Под ортогональными оценками автор понимает оценки, которым соответствуют ортогональные векторы. Легко показать, что ортогональные оценки представляют собой случайные величины, которые взаимно некоррелирова'ны.— Примред.


506

Глава 8

Ортогональные планы экспериментов представляют собой такие схемы независимых переменных, при которых для всех пар разных индексов J ж к суммы по наборам данных і — 1, 2, . . ., п обра­ щаются в нуль 1 ) , т. е. 2 xijxik — 0 П Р И 1 Ф к. В последующем

і

изложении под термином размерность будет пониматься число геометрических координат, необходимое для описания отклика, исключая координату самого отклика. Термин порядок модели будет означать степень полинома, составляющего модель (первая степень — первый порядок, вторая степень — второй порядок и т. д.). Порядок плана связан с порядком модели, т. е. план пер­ вого порядка можно использовать для оценивания линейных моде­ лей первого порядка с размерностью, равной 1, 2 и больше, но не для моделей второго порядка. Хотя планы второго порядка и мож­ но использовать для оценивания моделей первого порядка, они потребовали бы большего числа опытов, чем действительно необ­ ходимо.

Наиболее простой двухуровневый факторный 2 ) план экспери­ мента (план 22 ), пригодный для двумерной модели первого порядка

Y = ß 0 z 0

+ ß A + ß 2 x 2 + s,

(8.1.1)

имеет вид (в кодированной

форме)

 

Уровни й == 1 для всех уровней)

—1

—1

1

—1

—1

1

1

1

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

г )

Это определение справедливо только для планов первого

порядка

при условии одинакового числа повторяющихся опытов в точках

п л а н а . —

Прим.

ред.

 

2 ) Фактором называется любая экспериментальная переменная, которую

изменяют от эксперимента к эксперименту (например, температура,

давление,

время, концентрация). Если этой экспериментальной переменной нельзя приписать определенное количественное значение, фактор может носить

качественный

характер, быть

«большим» или

«малым», «присутствовать»

или

«отсутствовать».

Факторным

экспериментом

называется эксперимент,

в котором используются все возможные комбинации факторов.

Например,

если

к а ж д а я

из двух

переменных фиксируется

на

двух уровнях,

то имеют

место четыре возможные комбинации условий эксперимента; такой экспе­ римент называется двухуровневым факторным экспериментом, или, сокра­ щенно, планом 22 .