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

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

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

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

Добавлен: 09.04.2024

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

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

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

526

 

 

 

 

Глава

8

 

 

 

 

 

 

 

 

 

 

Таблица

П.8.1.6а

 

 

 

 

 

Число

Сумма

Отношение

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

 

 

степеней

 

 

квадратов

дисперсий

 

 

 

 

 

свободы

 

 

 

Влияние

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

 

1

 

36 551

 

 

Влияние

исключения by

 

1

 

431,7

 

 

Влияние

исключения

Ь2

 

1

 

70,8

 

38,24

Отклонения от эмпири­

 

1

 

370,2

 

ческой

линии

регрес­

 

 

 

 

 

 

сии

 

 

 

 

1

 

9,68

 

 

Отклонения в точке по­

 

 

 

 

вторения опытов (ошиб­

 

 

 

 

 

 

ка)

 

 

 

 

 

 

 

 

 

 

Общий

2 У |

 

5

 

37 433

 

 

Замечание: F0

9 0 (1,1) =

39,86

 

 

 

 

наблюдений

таковы:

 

 

 

 

 

 

 

 

 

Вектор

 

 

Матрица независимых

 

 

 

наблюде­

 

 

 

 

 

переменных

 

 

 

ний

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Y

 

XQ

 

 

 

 

 

 

96,0'

 

1

1,000

0

 

 

 

78,7

 

1

0,500

0,866

 

 

 

76,7

 

1

- 0 , 5 0 0

0,866

 

 

 

54,6

 

1

- 1 , 0 0 0

0

 

 

 

64,8

 

1

—0,500

-0,866

 

 

 

78,9

 

1

0,500

-0,866

 

 

 

97,4

 

1

0

 

0

 

 

 

90,5

 

1

0

 

0

 

 

 

93,0

 

1

0

 

0

 

 

 

86,3

 

1

0

 

0

 

 

Оценка уравнения регрессии

имеет вид

 

 

 

Y = 91,80 +

16,48 ХІ

+

3,38 х2 -

16,50 х\ - 17,20 х\

-

— 6,98 хух^-

Предположим теперь, что во втором блоке экспериментов значе­ ния отклика уменьшены на 10 единиц, т. е. второе наблюдение имело бы значение 78,7 — 10 = 68,7 и т. д. Тогда

Вектор

Матрица

независимых

наблюде­

 

переменных

ний

 

Y

х0

Хі

Х2

Г96,0Т

1

1,0.

0

68,7

1

0,5

0,866

76,7

1

- 0 , 5

0,866

44,6

1

—1,0

0

64,8

1

—0,5

-0,866

68,9

1

0,5

-0,866

97,4

1

0

0

80,5

1

0

0

93,0

1

0

0

.76,3

Li

0

0


Стратегия

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

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

527

В действительности по блокам распределены две различные

груп­

пы данных. Оценка

уравнения

регрессии

была бы такова:

 

Y = 76,48 + 16,48 ХІ + 3,38 х2

-

16,50 х\

- 17,20 х\ -

 

 

 

 

 

— 6,98

хух2.

За исключением коэффициента Ь0, подогнанная модель тождест­ венно совпадает с первой. Таким образом, форма поверхности оди­ накового уровня остается неизменной, несмотря на то что во вто­ ром блоке все данные были уменьшены на 10.

Для того чтобы выявить различие, если оно имеет место, между данными из двух блоков, можно провести дисперсионный анализ. Сумма квадратов, обусловленная влиянием блоков, вычисляется по формуле

 

 

10

Z i

»І

î o

i=l

 

 

где B t — полная сумма квадратов для блока i, a.nt — число наблю­ дений в блоке і:

(427,9)2

, (339,0)2 (766,9)2

 

 

 

 

 

 

- — g — — ~ — • - —

= 790,3; число степеней свободы равно 1.

Результаты

расширенного

дисперсионного

анализа

приведены

в табл. П.8.1.6.б.

 

 

 

 

 

 

 

 

 

 

 

 

Таблица

П.8.1.66

 

 

 

Число

 

Сумма

Средний

Отношение

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

степеней

квадратов

квадрат

дисперсий

 

 

 

свободы

 

 

 

 

 

 

 

 

 

Обусловленный

в л и я ­

1

 

790,3

 

Значимо

нием блоков

 

 

 

 

 

 

»

Влияние

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

1

 

58 813

 

 

Влияние

исключения

2

 

849,3

 

 

»

линейных членов

 

 

 

 

 

 

Влияние

исключения

3

 

718,4

 

 

»

членов

второго

поряд­

 

 

 

 

 

 

ка

 

 

1

 

1,60

 

 

 

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

эмпири­

 

1,60

Незначимо

ческой

л и н и и

регрес­

 

 

 

 

 

 

сии

 

 

2

 

18,5

9,25

 

 

Отклонения в точке по­

 

 

 

вторения опытов (ошиб­

 

 

 

 

 

 

ка)

 

 

 

 

 

 

 

 

 

Общий

 

10

 

61 191

 

 

 


528

Глава 8

Число различных наборов данных минус число коэффициентов составляет 1 степень свободы (7 6 = 1) для суммы квадратов отклонений от эмпирической линии регрессии. Сумма квадратов ошибок подсчитывается сначала внутри каждого из блоков и затем полученные результаты объединяются. Так как в каждом блоке имеется по два повторных измерения, сумма квадратов ошибок дается выражением 2 d%!2 или

(4,4)' + (4,2)» = 1 8 ) 5

с двумя степенями свободы. Видно, что неадекватность не является значимой, члены второго порядка оказываются значимыми и что влияние блоков дает значимый вклад. Однако различие между блоками не влияет на форму поверхности отклика, так как все оценки коэффициентов, связанных с переменными хх2, х\, х\ и Хіх2, остались прежними.

 

8.1.5. Практические

замечания

 

Для

некоторого непрерывного

процесса,

который

является

рабочим

процессом и должен таким

остаться,

при точном

исполь­

зовании методов поверхности отклика, описанных выше, можно

столкнуться с рядом

практических

трудностей.

1. Использование

ортогональных

планов может оказаться

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

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

Две другие проблемы, связанные с методами поверхности откли­ ка, касаются предположений, лежащих в основе используемого метода линейной регрессии. Эти проблемы еще не разрешены и по­

этому

будут

лишь кратко упомянуты здесь.

3.

Влияние

ошибок в уровнях переменных. Основное предположе­

ние регрессионного анализа, обсуждавшееся в гл. 4, состоит в том,

что все ненаблюдаемые ошибки содержит зависимая переменная

Y,

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

что

характер ошибок в переменных x-t совершенно неизвестен, и экспе­ риментальные данные анализируются так, как будто ошибок нет.


Стратегия

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

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

529

Бокс [7] показал, что стандартный статистический анализ справед­ лив, даже если экспериментальные уровни содержат ошибки, до тех пор, пока отклик линеен по х. Каждый случай квадратичной (по х) функции и других нелинейных функций должен быть рас­ смотрен отдельно, ибо величина ошибки, переносимой из пере­ менных X; в отклик, зависит от наклона поверхности отклика в со­ ответствующей точке плана. Двухуровневый факторный план и некоторые дробные факторные планы, по-видимому, весьма устойчивы по отношению к ошибкам в переменных х-г.

4. Использование планов для коррелированных данных. Второе основное предположение регрессионного анализа, используемое в методах поверхности отклика, состоит в том, что эксперимен­ тальные отклики некоррелированы; однако хорошо известно, что данные, полученные для непрерывных процессов, весьма часто подвержены влиянию различных трендов или сдвигов среднего уровня (низкочастотные изменения), а также сериальной или вре­ менной корреляции, несмотря на постоянные попытки избежать этого. Некоторые трудности анализа коррелированных данных описаны в разд. 5.4.

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

8.2. К А Н О Н И Ч Е С К И Й А Н А Л И З

В этом разделе будет дана графическая и аналитическая интер­ претация поверхности отклика (т. е. поверхности, описываемой оценкой уравнения регрессии). Поверхности, уравнения которых выражаются функциями двух экспериментальных переменных, нетрудно изобразить на двумерной плоскости, на которой отклик представляется в виде контуров постоянного значения. Поверхно­ сти, которые описываются уравнениями, содержащими три экспе­ риментальные переменные, также можно изобразить графически, хотя это и не так просто, как в предыдущем случае. Линд и Янг [8] предложили трехмерное устройство для воспроизведения повер­ хности отклика; можно составить машинную программу для построения графиков поверхностей отклика и воспроизведения их с помощью осциллографов. Если уравнение поверхности откли­ ка содержит больше трех экспериментальных переменных, то гра­ фическое изображение этой поверхности связано с большими труд­ ностями.

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