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

Методы оптимизации эксперимента в химической технологии

..pdf
Скачиваний:
41
Добавлен:
15.11.2022
Размер:
13.27 Mб
Скачать

между средними результатами предыдущего цикла и новыми результа­ тами.

4. Встроке5 блока «Расчет средних» записывают новые суммы резуль­ татов наблюдений для каждой точки плана, т. е. суммы строк 1 и 3.

5. В строке 6 записывают новые средние для каждой точки плана как частные от деления новых сумм (строка 5) на число циклов к данному времени, т. е. в нашем случае частное от деления данных строки 5 на п =*=2.

6. По формулам (V. 185) —(V. 188) рассчитывают эффекты факторов, их взаимодействий и эффект изменения среднего. Полученные значения заносят в блок «Расчет эффектов».

7. В блоке «Расчет ошибки опыта» для второго цикла строки 1 и 2 не заполнены; в строку 3 записывают размах величин 5Ду, определен­ ный как разница между подчеркнутыми в строке 4 блока «Расчет сред­ них» наибольшим и наименьшим значениями разностей Ау. В нашем случае размах равен 5Ау —1,5 —(—2,0) =3,5.

8. В строку 4 записывают величину новой ошибки опыта, которую

считают методом размаха по формуле (V.183). Величину fNn

находим

из табл. 63. В нашем случае

5, л = 2, ^ 2=030. Величина

ошибки

s = 3,5 -0,30 = 1,5. К концу второго цикла еще нет оценок ошибки опыта из предыдущего цикла, поэтому в строки 5 и 6 блока «Расчет ошибки опыта» переписывают величину новой ошибки (строка 4).

9. По формулам (V.189) —(V.192) определяют доверительные интер­ валы для нового среднего, эффектов и эффекта изменения среднего, которые записывают, соответственно, в строки 1, 2 и 3 блока «Расчет доверительных интервалов».

10. Сравнивают величину эффектов с соответствующими доверитель­ ными интервалами и делают вывод о статистической значимости эффек­ тов, после чего принимают решение о следующем этапе.

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

Начиная с третьего цикла, заполняются все строки рабочей таблицы (табл. 68). Все вычисления проводятся так же, как для предыдущих циклов, только в блоке «Расчет ошибки опыта» в строки 1 и 2 пере­ писывают результаты из строк 5 и 6 того же блока предыдущего цикла. В строку 5 вносят результаты, полученные суммированием строк 2 и 4, строку 6 получают делением строки 5 на п - 1. Оценка результатов третьего цикла свидетельствует о значимом отличии от нуля эффекта кислотности. В данной ситуации принимается одно из двух возможных решений.

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

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

Врассматриваемом примере (табл. 68) было принято решение начать новую фазу с центром плана в точке 3 (рис. 57).

Вновой фазе планирования при необходимости можно изменить

 

Р а с ч е т с р е д н и х

 

 

 

 

 

Р а с ч е т о ш и б к и о п ы т а

 

 

Н о м е р т о ч е к э к с п е р и м е н т а

 

 

 

 

1

2

3

4

5

 

 

1. Сумма предыдущих

 

 

 

 

 

1. Сумма ошибок предыду­

циклов

12,0

12,0

9,5

10,5

12,0

щих циклов J —1,05

2.

Среднее предыду­

 

 

 

 

 

2.

Средняя ошибка предыду­

щих циклов

6,0

6,0

4,75

5,25

6,0

щих циклов Sy —1,05

3.

Новые наблюдения

5,5

6,5

4,2

4,0

6,0

3.

Размах —1,75

4.

Разности (2) —(3)

0,5

-0,5

0,55

1,25

0

4.

Новые значения ошибки

5.

Новые суммы (1)+ (3)

17,5

18,5

13,7

14,5

18,0

sy = размах fN n= 0,6

5.

Новая сумма ошибок s — 1,65

6.

Новые средние

 

 

 

 

 

6.

Новая средняя ошибка

 

s

 

 

 

 

 

II

|

 

15

5,8 6,2 4,6 4,8 6,0 * У - ) - 1 “ 0,81

Р а с ч е т э ф ф е к т о в

Р а с ч е т д о в е р и т е л ь н ы х

 

и н т е р в а л о в (Д И )

 

1.

Эффект кислотности

1. ДИ

для новых средних

В\

= - у

(уз + Ул - У2 - уъ) - 1 , 4

2

2 • 0,81

0,90

/2 ' у

Sb

 

 

Эффект нагрузки

 

 

2.

2. ДИ для эффектов

Я2 = - у

6* + Уб - У2 - ул) - -0,2

^ , - 0 , 9 4

 

 

3.

Эффект взаимодействия

3. ДИ для эффекта изменения

#12

=

(уг + уз - ул - уъ) —0,0

среднего—’l l sy

0,81

 

4 .

Эффект изменения среднего

 

Vn-

 

 

 

 

 

М

=

 

2 + уз +~УА + уъ - 4У0 = 0 ,3 2

 

 

 

 

Выводы: эффект кислотности значим

Априорная ошибка опыта

 

 

 

 

Sy “ 0,5

 

 

Решение: начать новую фазу

интервалы варьирования факторов, можно изменить и сами факторы: некоторые из прежних оставить на выбранном оптимальном уровне и включить в программу эксперимента новые факторы; можно добавить

 

<1

факторы к уже имеющимся. Наиболее часто

 

раза

при эволюционном планировании варьируют

13

 

два или три

фактора.

Расчеты

для трех

*

1 фаза

факторов во многом аналогичны двухфак­

<Г)

CL

— <

торной задаче. Для каждого фактора выби­

«о

 

рают по два уровня варьирования. Цикл

 

Кислотность (х,)

представляет

собой

полный

факторный

 

эксперимент 23 (рис. 58). Однако в данном

Рис. 57. Схема перехода ко

случае цикл делят на два блока, а в центре

 

второй фазе

плана проводят два опыта (табл. 69).

Каждый блок является дробной репли­ кой типа 23-1 и опытом в центре плана. Блок I представляет собой полуреплику 23- 1 с определяющим контрастом 1 = —х\Х2Х$.

Это позволяет получить следующую систему оценок:

“"*■ Э1 — Р2 3 »

^ ^ P a - P i a »

(V. 193)

Б

ad

 

1

I

1

1

1,6

/l 7-

П

 

O —

 

 

Нагрузка (х2)

 

^3

Рз —

Р12

 

 

Рис. 58. Схема цикла для /г—3

 

 

Т а б л и ц а

69. Матрица планирования

 

 

 

Н о м е р

X \

X I

X3

XiX2

X1X3

X2X3

X \X 2X 3

У

 

о п ы т а

 

 

 

 

 

 

 

 

 

1

0

0

0

0

0

0

0

yy

 

2

- 1

- 1

- 1

+ 1

+ 1

+ 1

- 1

У2

т

3

+ 1

+ 1

- 1

+ 1

- 1

- 1

- 1

 

1

У з

блок

4

+ 1

- 1

+ 1

- 1

+ 1

- 1

- 1

У4

 

5

- 1

+ 1

+ 1

- 1

- 1

+ 1

- 1

У5

 

6

0

0

0

0

0

0

0

У6

II

7

- 1

- 1

+ 1

+ 1

- 1

- 1

+ 1

У7

8

+ 1

+ 1

+ 1

+ 1

+ 1

+ 1

+ 1

У*

блок

9

+ 1

- 1

- 1

- 1

- I

+ 1

+ 1

Ув

 

10

- 1

+ 1

- 1

- 1

+ 1

- 1

+ 1

^ 1 0

 

Опыты блока II составляют другую полуреплику 23 1 с определяющим контрастом 1 Система оценок следующая:

В \ 1 -> Pi +

Раз»

 

Я ^ - Р . + Ри.

<V-,94>

“*■Ре +

Эле*

 

Объединение опытов блоков I и II приводит к полному факторному эксперименту 23, который дает возможность получить несмешанные оценки для линейных эффектов и эффектов парного взаимодействия:

= y (* ' + * " ) - * Pl*

 

В*= ~2~(В2+ В2) Р*'

 

В3= J - ( B»+ S » ) ^ P 3.

(V.195)

Blt = Y (В? - Bl ) -+

* а з = “ j " ( * i I — * 1 )

Рад*

Рассмотрим порядок заполнения рабочих таблиц в трехфакторной задаче на примере промышленного эксперимента по минимизации потерь связанного азота. Вторую фазу эксперимента было решено начать с центром в точке 3 (см. рис. 58). Введем в рассмотрение третий фактор — температуру (хд). Основной уровень, интервалы варьирования и границы области исследования приведены в табл. 70.

Т а б л и ц а 70. Область исследования

Ф а к т о р ы

К и с л о т н о с т ь ,

Н а гр у зк а п о

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

 

Г / л

Ы Н з, м э /ч

 

Основной уровень

1,1

1930

65

Интервал варьиро­

 

 

 

вания

 

 

0 ,4

80

5

X :

=

+ 1

1,5

2010

70

х )

1

0,7

1850

60

 

 

 

В результате трех циклов второй фазы получены следующие значе­

ния отклика у ь г/л

(табл. 71).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Т а б л и ц а 71

Блок

 

 

I

 

 

 

 

11

 

 

Номера точек

1

2

3

4

5

6

7

8

9

10

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

 

 

 

 

 

 

 

 

 

 

Цикл 1

5,8

6,2

4,8

6,1

6,0

6,5

5,9

5,5

5,8

6,3

Цикл

2

6,2

5,6

5,9

5,5

5,6

4,9

5,7

5,1

5,0

6,3

Цикл 3

4,5

6,2

4,0

5,0

5,9

6,6

5,8

4,3

4,7

6,5

Расчеты для каждого цикла выполняются в двух рабочих таблицах, отдельно для каждого блока. В табл. 72 приведена рабочая таблица блока I третьего цикла второй фазы промышленного эксперимента по минимизации потерь связанного азота в производстве аммиачной селит­ ры, а в табл. 73 —блока II того же цикла той же фазы. В блоке «Расчет эффектов» (табл. 72, 73) вычисляются оценки эффектов в соответствии с (V. 193) и (V. 194). В блоке «Окончательный расчет эффектов» рабочей таблицы для блока II усредняются оценки для эффектов в соответствии с (V.195). Здесь же рассчитывается эффект изменения среднего по данным полного цикла усреднением соответствующих эффектов для блоков I и II.

Расчет ошибки опыта проводится по данным каждого из блоков, при этом размах определяется как сумма наибольших и наименьших

264

 

Р а с ч е т с р е д н и х

 

 

 

 

Р а с ч е т о ш и б к и о п ы т а

 

 

Н о м е р а т о ч е к э к с п е р и м е н т а

 

 

 

 

1

2

3

4

5

 

 

1. Сумма предыдущих

 

 

 

 

 

1. Сумма ошибок предыдущих

циклов блока I

12,0

11,8

10,7

11,6

11,6

циклов (все блоки) 5 —1,05

2.

Среднее предыду­

 

 

 

 

 

2.

Средняя ошибка предыду­

щих циклов блока I

6,0

5,9

5,4

5,8

5,8

щих циклов (все блоки) s y —0,52

3.

Новые наблюдения

4,5

6,2

4,0

5,0

5,9

3.

Размах - 1,8;/б^-0,35 .

4.

Разности (2) —(3)

1,5

-0,3

1,4

0,8

-0,1

4.

Новое значение ошибки

 

 

 

 

 

 

 

s у -

размах /ли- °,63

5.

Новые суммы блока I

 

 

 

 

 

5.

Новая сумма ошибок (все

(1)4- (3)

16,5

18,0

14,7

16,6

17,5

блоки) 5 —1,68

6.

Новые средние

 

 

 

 

 

6.

Новая средняя ошибка

блока I У,- = ~ -

5,5

6,0

4,9

5,5

5,8

* -

2п - } = 0’56

 

' ' п

 

 

 

 

 

 

Р а с ч е т э ф ф е к т о в

 

 

 

 

 

1.

Эффект кислотности

 

 

 

 

 

 

 

в} =

у (7з + У4 - У2 - уь)---- 0,7

 

 

 

 

 

 

2.Эффект нагрузки

в\ = у (Уз + У5 - У2 - Ул) ——0,4

3.Эффект температуры

=у 2 + уз - ул - уь) ~~—0,2

4.Эффект изменения среднего блока I

М= у (уг + уз + У4 + уь - 4yij - 0,04

Выводы:

разностей данного блока. При определении сумм ошибок используют данные по всем уже реализованным блокам данной фазы.

Расчет доверительных интервалов выполняется только в рабочей таблице блока II (табл. 72). Доверительные интервалы в данном случае рассчитываются по следующим формулам:

для новых средних

Sy ,

(V. 196)

для эффектов

 

ДИ = 0,71 - 2 — sy.

(V. 197)

Уп

 

для эффекта изменения среднего

 

ДИ = 0,63 — — sv.

(V.198)

Уп

 

10-529

 

 

Расчет средних

 

 

 

 

 

 

Расчет ошибки опыта

 

 

 

 

 

Номера точек эксперимента

 

 

 

 

 

 

 

 

 

 

 

1

2

3

4

5

 

 

 

 

 

 

1. Сумма предыдущих

 

 

 

 

 

 

1. Сумма ошибок

предыдущих

циклов

 

 

 

11,4

11,6

10,6

10,8

12,6

циклов (все блоки) $ “

1,68

2.

Среднее предыду­

 

 

 

 

 

 

 

2.

Средняя

ошибка

предыду­

щих циклов

 

 

 

5,7

5,8

5,3

5,4

6,3

щих циклов (все блоки) s y “ 0,56

3.

Новые наблюдения

6,6

5,8

4,3

4,7

6,5

 

3.

Размах —1,9

 

 

4.

Разности (2) —(3)

 

 

0

0,7

-0,2

 

4.

Новое значение ошибки

5.

Новые суммы

 

 

 

 

 

 

 

s y

-

размах fNn—0,67

 

 

 

 

18,0

17,4

14,9

 

 

 

5.

Новая сумма ошибок (все

(1) +

(3)

 

 

 

15,5

19,1

блоки) s —2,35

 

 

 

6.

Новые средние

 

 

6,0

5,8

4,9

5,1

6,7

 

6.

Новая средняя ошибка

 

 

 

 

 

 

 

 

 

 

 

 

2л - 2 =

0 '59

 

 

 

 

Расчет эффектов

 

 

 

 

 

Расчет доверительных

 

 

 

 

 

 

 

 

 

 

 

 

интервалов (ДИ)

1. Эффект кислотности

 

 

 

 

 

 

1. ДИ для новых средних

В]^= - у - (уъ + УЪ -

У7 -

ую ) —1,25

 

 

 

- 7 ^ - 0 , 6 9

 

 

 

2.

Эффект нагрузки

 

 

 

 

 

 

 

2. ДИ для

эффектов

£:Р =

- у - (уъ + ую - у? -

уъ) “ + 0 ,3 5

 

 

 

0,71

—0,49

 

 

3.

Эффект температуры

 

 

 

 

 

 

3. ДИ для эффекта изменения

Вз^= - у (У7 + ув -

уъ -

ую ) — —0,55

 

 

 

среднего 0,63 —p= sy = 0,43

4.

Эффект изменения среднего

 

 

 

 

 

 

 

 

 

Цч = - у - ( у ? + у в + у ъ + у ю -

4ув)

—0,3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Окончательный расчет эффектов

 

 

 

 

 

В ,

=

- у

( Я ?

+

- “ 0 , 9 8 ;

 

=

—у

(Д з — Дз) “

—0,15

 

 

 

Вч

=

 

№ * + В ч )

- - 0 ,0 2 5 ;

 

Дю = - j -

(Дг**- Вч)

+ 0 ,3 7 5

 

 

 

 

=

—jjn (Въ

+ В ^ )

— —0,375;

 

Вчз = - j -

( В Г - В}) -

-0,275

 

 

 

 

Эффект

изменения

среднего

ц«=—j- (Mi + Мг) “ -0,13

 

 

Выводы: эффект

кислотности значим

 

 

 

 

Решение: начать новую фазу

 

В рассматриваемом примере после третьего цикла второй фазы

(табл. 72) статистически

значимым

оказался

эффект

Ву

(0,98 >0,49),

поэтому было принято решение перейти к новой фазе.

 

 

 

 

Таким образом, метод эволюционного

планирования дает возмож­

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

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

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

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

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

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

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

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

Упражнения

1.Какими свойствами обладает матрица планирования полного факторного экспери­ мента 2*? Построить матрицы планирования для к —4 и к —5.

2.Построить линейные ортогональные планы для к —6, используя дробные реплики

от полного факторного эксперимента и симплексные планы. Какими свойствами обла­

дают эти планы?

3. Полный факторный эксперимент 22, приведенный в таблице:

Номер

21

Z2

 

У

опыта

 

 

 

 

1

80

11,34

83,1

85,2

2

20

9,75

60,6.

62,5

3

80

9,75

71,8

73,9

4

20

11,34

83,7

81,9

использовался для изучения зависимости соотношения между водной и общей форма­ ми РгОб (у,%) от температуры процесса аммонизации (zi, °С) и содержания воды в спирто­ вой фазе (z2, %) при получении монокальцийфосфата кислотным разложением фосфатов

с применением жидкостной экстракции. Каждый опыт повторялся 2 раза.

„ пеГпессии в безразмерном и натуральном масштабах;

б) — — регрессйИ

экспериментх^^^^^^ С^^Р^о/С^о™’в^^ениМхран^нма™Т%ут):Л риВ^ л и ч н о ^ И|^Двр1-

 

Z\

Z2

^min

5

60

Z j max

20

72

Был использован ортогональный план 2-го порядка (табл. 43), получены следующие значения коэффициента разложения:

Н о м е р

у, %

Н о м е р

у, %

о п ы т а

 

о п ы т а

 

 

 

 

1

96

6

85

2

85

7

93

3

76

8

81

4

91

9

88

5

1

 

 

 

 

 

Дисперсия воспроизводимости s?BOCnp равна 4%, /воспр —3.

Определить: а) уравнение регрессии второго порядка; б) значимость коэффициентов и адекватность уравнения регрессии эксперименту; в) тип поверхности; г) оптимальные условия, обеспечивающие практически полное разложение (у =* 100%); Д) проанализировать

чувствительность коэффициента разложения к изменению времени хранения и содер-

жаниюНаБОд.

5. По условиям примера 12 для исследования кинетики сополимеризации; д) соста­ вить D-оптимальный симплекс; б) используя результаты эксперимента (табл. 61) и свой­

ства симплексного плана, начать движение в сторону увеличения скорости сополиме­ ризации.

ГЛ АВ А VI

ПЛАНИРОВАНИЕ ЭКСПЕРИМЕНТА ПРИ ИЗУЧЕНИИ ДИАГРАММ СОСТАВ - СВОЙСТВО

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

редставляет собой правильный (q - 1)-мерный симплекс. Для систем выполняется соотношение

 

 

Ч

 

(VI-. 1)

 

 

2 *i =

i.

 

 

(=1

 

 

Дел";я^ 1 Г

Г

еНТРаЦИЯ К0МП0Нента; Я-количество компонентов.

нентов До Г

л Г Г Г ЫХ °ИСТеМ симплекс - прямая содержание компо-

c^rS T eK r-nf

соотношением

отрезков. При q - Ъ

правильный

вечает одномуНг т ^ Р°ННИИ треугольникКаждая точка треугольника от­ вечает одному определенному составу тройной системы и, наоборот,

268

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

В

В

л л

г. / у v л ' л ' л я

 

 

ш

 

.

 

 

 

/ \ 7 \ Ж Ж

 

 

 

 

 

/

Nr

^ у \ Y v &

J \

 

у у \

А

\

 

 

 

П4/

\ Л

V

 

 

V V i V

>

 

 

 

 

A

 

b

 

d

С

АЛ

 

 

 

 

 

 

 

 

 

 

 

 

С ,%

Рис.

59. Концентрационный треуголь­

Рис. 60. Концентрационный треуголь­

 

ник Гиббса

 

 

ник Розебума

10%. На рис. 59 на соответствующих высотах тпеугольника указано содержание каждого из компонентов (в точке d 30% А и 70% С). Практи­ чески большей частью не прибегают к построению высот, можно откладывать содержание компонентов непосредственно на сторонах треугольника. Такой способ отсчета принят в треугольнике Гиббса. В треугольнике Розебума состав тройной системы отсчитывается по трем отрезкам одной стороны треугольника (рис. 60). В концентрацион­ ном треугольнике точки, лежащие на прямой, выходящей из вершины треугольника, соответствуют смесям с постоянным отношением содер­ жаний компонентов, изображаемых двумя другими вершинами. Свойство

(у) обычно представляют проекциями линий равного значения на плос­ кость концентрационного треугольника.

При 0е 4 правильный симплекс —тетраэдр, каждая вершина которого соответствует чистым компонентам. Ребро представляет собой двух­ компонентную систему, грань —трехкомпонентную. Точки внутри тет­ раэдра соответствуют четырехкомпонентным системам. Так, компонент х } отсутствует на грани х2, х& х4, а по сечениям тетраэдра, приближаю­ щимся к верщине х,, содержание компонента х} увеличивается.

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

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

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

Поверхности отклика в многокомпонентных системах имеют, как правило, очень сложный характер. Для адекватного описания таких поверхностей необходимы полиномы высоких степеней и, следователь­ но, большое количество опытов. Обычный полином степени п от q переменных имеет Cqn+n коэффициентов:

у = 6 0 + 2 <W +

2

bUxi*1+

2 ' bUh*lx)xh + " ’ (У1-2)

\< i< q

K i < J < q

 

\< t < f < k < q

Я

Соотношение 2 xt = 1позволяет исключить q-й компонент и снизить

i= 1

число коэффициентов до Cq+n_v Однако по существу задачи желательно ввести в модель все lq компонентов.

Шеффе предложил описывать свойства смесей приведенными поли­ номами, получаемыми из (IV.2) с учетом условия нормированное™ суммы независимых переменных (VI. 1). Покажем, например, как полу­ чить такой приведенный полином второй степени для тройной системы. Общий вид полинома

У = Ьц+ Ьххх+ Ь2х2+ М з + Ъ12ххх2+ Ьхзххх8+ Ьаз*а*з +

+ &ii*i +

^22*2 "Ь Ья з 4 -

(VI.3)

Так как

xt + *3 = 1.

(VI. 4)

*1 +

7 0

 

 

V i + V * + &o*s = bo’

(VI. 5)

Умножив (VL4) последовательно на х,, х2 и х3, получим

 

*1 = * 1 — •*!** — *1 *8 .

 

*2= **— *i** — *2*в•

(VI.6)

•^ = *3 — *1*8 — *2 *3 •