Файл: Горюшко, В. Е. Планирование эксперимента в бытовой химии [обзор].pdf

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

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

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

Добавлен: 01.11.2024

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

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

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

Т а б л и ц а II.10

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

 

Концентрация отбеливателя, %

 

Итоги

Средние

Квадраты

 

 

 

 

 

значения

И Т О Г О В

 

0,1

0,2

0,3

0,4

0,5

Ai

 

А-

А]

 

 

 

 

 

 

 

20

96,20

105,40

104,60

118,23

92,80

517,24

103,50

267537,20

30

103,80

102,77

104,30

97,58

108,08

516,53

103,30

266803,24

40

97,35

83,60

122,30

125,00

119,40

547,65

109,50

299920,50

50

99,70

112,80

85,55

123,30

116,47

537,82

107,60

289250,35

60

84,72

115,87

114,22

110,20

123,90

548,90

109,80

301302,20

Итоги Bi

481,77

520,44

530,97

574,32

560,65

Средние значе­

96,35

104,09

106,20

114,86

112,13

ния B-t

 

 

 

 

 

 

 

 

Квадраты итогов

232102,3

270857,8

281929,1

329843,5

314328,4

В]

 

 

 

 

 

 

 

 


Т а б л и ц а II.11

 

Эффект латинских букв

 

 

Итоги

Средние

Квадраты

Фактор

 

 

 

 

 

значения

итогов

 

 

 

 

 

Си

Си

Г 2

 

 

 

 

 

 

 

 

А

96,20

108,08

122,30

123,30

115,87

565,75

113,15

320073,06

В

105,40

103,80

125,00

116,47

114,22

564,82

112,97

319100,71

С

104,60

102,77

119,40

99,70

110,20

536,67

107,30

288014,69

D

118,23

104,30

97,35

112,80

123,90

556,58

111,30

309781,30

Е

92,80

97,58

83,60

85,55

84,72

444,25

88,85

197358,06

 

Эффект греческих букв

 

Итоги

Средние

Квадраты

 

 

 

 

 

 

Фактор

 

 

 

 

 

значения

итогов

 

 

 

 

 

D m

D m

D 1

 

 

 

 

 

 

 

т

а

96,20

104,30

83,60

116,47

110,20

510,80

102,1

260885,99

р

105,40

97,58

122,30

99,70

123,90

548,90

109,8

301269,25

У

104,60

108,08

125,00

112,80

84,72

535,92

107,0

286439,04

8

118,23

103,80

119,40

85,55

115,87

542,85

108,6

287210,24

£

92,80

102,77

97,35

123,30

114,22

530,44

106,1

281366,59

Используя результаты, представленные в указанных таблицах, вычисляем вспомогательные суммы квадратов:

 

 

ЛГ

 

 

 

 

 

 

4

^ =

288515>35;

55.'>а__

2

Л?

 

 

 

 

*-1

1424813,59

284962,72;

 

 

п

5

 

 

 

 

I

*,2

1429061,14

 

01

55,

7=1

=

5

 

285812, 23;

 

 

п

 

~

 

5 5.

t

e l

 

 

 

 

U=1

1434327,82

= 286865,56;

 

 

 

5

 

 

 

й

*

т= 1

1424646,99

284929,40;

 

 

 

 

55e =

7119924,42 = 284760,98.

 

 

25

 

33


Суммы квадратов, характеризующие эффекты строк, столбцов, латинских и греческих букв, общая и остаточная суммы квадратов, соответственно равны:

SSa = SS2 - 556= 201,74;

SS„= SS3- SS6 = 1051,25;

SS£ = SSi — 55 6= 2104,58;

SSd = SSb- S S 6= 168,42;

5So6l4 = 5 5 1 - 5 5 e = 3754,37;

5 5 ^ = 228,38.

Результаты дисперсионного анализа сведены в табл. 11.12.

 

 

 

 

 

Т а б л и ц а

II. 12

 

Число

Сумма

Средний

 

Fтабл.

 

Источник дисперсии

степеней

квад­

Fэксп

 

 

квадрат

Я =95% Р=80%

 

свободы

ратов

 

Строки

4

201,74

50,40

1,76

3,8

1,9

Столбцы

4

1051,25

262,80

9,19

3,8

1,9

Латинские буквы

4

2104,58

526,14

18,39

3,8

1,9

Греческие буквы

4

168,42

42,10

1,47

3,8

1,9

Остаток

8

228,38

28,60

Таким образом, на отбеливающую способность оказывают су­ щественное влияние тип (латинские буквы) и концентрация (столб­ цы) отбеливателя. Наиболее сильно влияние типа отбеливателя. Из рассмотрения эффектов латинских букв следует, что наиболее велики итоги по эффекту действия отбеливателей типа А и В (см.

табл. 11.11, итоги 565,75 и 564,82).

Из табл. 11.12 видно, что эффект греческих букв значимо не проявился. Очевидно, что отбеливающая способность не зависит от типа ПАВ.

Итак, в главе II рассмотрен метод дисперсионного анализа. В исследовательской работе влияние «неуправляемых» факторов на эксперимент пытаются уменьшить, пользуясь методом рандомиза­ ции. Однако при известных источниках неоднородности более эф­ фективны планы дисперсионного анализа. Применение этих планов целесообразно, когда необходимо построить наиболее экономные схемы планирования эксперимента с качественными и количествен­ ными факторами.

34


Г Л А В А III

ПЛАНИРОВАНИЕ МНОГОФАКТОРНОГО ЭКСПЕРИМЕНТА

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

На рис. 2 представлена поверхность гипотетической функции y = f(x ь х2), описываемая горизонталями. Сравним два метода дви­ жения к оптимуму.

Рис. 2. Поверхность функции отклика и пути движения к оптимуму.

При классическом подходе к планированию будем поочередно изменять каждый параметр до достижения частного максимума. Так, двигаясь по прямой АВ, фиксируем точку максимума С. При движении параллельно оси х2 достигаем максимума в точке D и т. д.

Очевидными недостатками этого метода являются длинный путь движения к оптимуму и необходимость фиксировать п— 1 па-

35

раметров. При сложной геометрии функции отклика метод вовсе не приводит к цели.

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

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

1. Требования к параметру оптимизации, факторам и модели

Математически задачу оптимизации можно сформулировать как задачу отыскания наибольшего (наименьшего) значения функции нескольких переменных при определенных ограничениях по дру­ гим параметрам:

У = / { х 1, Xj,..., .*,•)= max (min),

где Xi — независимые переменные

объекта исследова­

ния, или факторы;

 

f(xь х2, ..., Xi) — функция отклика;

параметра оптими­

у — количественная оценка

зации.

 

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

Объект оптимизации и, естественно, все факторы должны быть

управляемыми, а опыты — воспроизводимыми.

факторов являются

Основными требованиями к совокупности

совместимость и отсутствие линейной зависимости между ними.

 

Например, при оптимизации состава герметика факторами мо­

гут служить количества парахинондиоксима,

двуокиси

марганца,

цемента, растворителя БР-1, гипериза и нафтената

кобальта

в

граммах на 100 г бутилкаучука.

 

 

 

Под математической моделью процесса понимают функцию от­

клика вида

 

 

 

У= ft)+ ?1Х1 + &5-*2+ • •' + Pi2-*l-*2 + ?13Х1Х3~Ь • •' +

§jk...mXjXk----Хт- (

1)

Описываемая этим уравнением нелинейная поверхность назы­ вается поверхностью отклика.

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

(III.1).

36


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

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

2. Полный факторный эксперимент

Процедура планирования эксперимента — это последователь­ ность этапов, результаты которых детально анализируются с целью принятия решений.

Перед тем как запланировать эксперимент, необходимо оценить области допустимого изменения факторов и увязать их с техноло­ гией процесса и техникой эксперимента.

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

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

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

уровень фактора.

Масштабы факторов выбирают так, чтобы верхний уровень со­ ответствовал + 1, нижний —1, основной — нулю.

Кодированное значение фактора

 

 

 

(III. 2)

где Xj и Xjo — натуральные значения соответственно

фактора и

основного уровня;

 

Xj — шаг варьирования.

нижних Xj =

Для верхних

уровней Xj= + \ и Xj = Xj0 + Xj, для

-----1 и Xj Xjo

Xj>

 

Например, при основном уровне содержания цемента в компо­

зиции герметика Xjo= 54 г и интервале варьирования

12 г верхний

уровень Xj= 5 4 + 12= 66 г и нижний х3-=54—12= 42 г.

 

37