Планирование эксперимента в химической технологии
..pdf№ п/п
1
2
3
4
5
Границы |
Число эле |
|
№ п/п |
Границы |
Число эле |
Частота |
||
ментов в |
Частота |
ментов в |
||||||
интервала |
интервала |
|||||||
интервале |
|
|
интервале |
|
||||
0 |
со |
16 |
0,Ю7 |
6 |
50—59 |
19 |
0,127 |
|
|
1 |
|
0,100 |
7 |
60—69 |
14 |
0,093 |
|
10— 19 |
15 |
|||||||
20—29 |
19 |
0,127 |
8 |
70—79 |
И |
0,073 |
||
30—39 |
13 |
0,087 |
9 |
80—89 |
13 |
0,087 |
||
40—49 |
14 |
0,093 |
10 |
90—99 |
16 |
0,107 |
||
|
|
|
|
|
|
|
Таблица |
15. |
Результаты наблюдений за среднесуточной температурой |
воздуха |
|||||||||||
№ |
|
|
Границы |
|
Число эле |
№ |
|
Границы |
|
Число эле |
||||
п/п |
|
|
интервалов |
|
ментов в |
|
п/п |
|
интервалов |
ментов в |
||||
|
|
|
|
|
|
интервале |
|
|
|
|
|
|
интервале |
|
1 |
|
от |
—40 до —30 |
5 |
|
6 |
ОТ +10 |
до +20 |
|
81 |
||||
2 |
|
от |
—30 до |
—20 |
11 |
|
7 |
от |
+20 |
до |
+30 |
|
36 |
|
3 |
|
от —20 до — 10 |
25 |
|
8 |
от +30 |
до +40 |
|
20 |
|||||
4 |
|
от |
—10 до |
0 |
42 |
|
9 |
от |
+ 4 0 |
до +50 |
|
8 |
||
5 |
|
от |
|
0 до |
+10 |
88 |
|
10 |
от +50 |
до +(Ю |
|
4 |
||
Таблица |
16. |
Результаты ранжирования |
факторов |
|
|
|
|
|
|
|||||
Эксперты |
|
|
|
|
|
|
Факторы |
|
|
|
|
|
|
|
|
|
|
x t |
|
Х4 |
ХЬ |
х* |
|
*7 |
X* |
|
xt |
*10 |
|
|
|
|
|
|
|
|
||||||||
1 |
|
1 |
|
4,5 |
2 |
4,5 |
3 |
7,5 |
|
6 |
9 |
|
7,5 |
10 |
2 |
|
2,5 |
|
1 |
2,5 |
4,5 |
4,5 |
8 |
|
9 |
6,5 |
10 |
6,5 |
|
3 |
|
2 |
|
1 |
4,5 |
4,5 |
4,5 |
4,5 |
|
8 |
8 |
|
8 |
10 |
Таблица |
17. |
Характеристика |
Таблица |
18. |
Данные |
о |
качестве |
|||||||
прочности на разрыв |
для |
четырех |
продукта |
в опытах |
|
|
|
|
||||||
сплавов |
|
|
|
|
|
|
Микро |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
2 |
|
|
3 |
|
Номер |
|
|
|
|
|
примесь |
|
|
|
|||||
|
|
Уш |
Ур |
У, |
на уровне |
|
|
|
|
|
|
|||
сплава |
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
1 |
0,15 |
0,16 |
0,15 |
|||
1 |
|
0,12 |
|
0,15 |
0,33 |
0,9 |
|
2 |
0,22 |
0,17 |
0,20 |
|||
2 |
|
1,2 |
|
5,3 |
2,1 |
2,5 |
|
3 |
0,21 |
0,23 |
0,20 |
|||
|
|
|
4 |
0,35 |
0,30 |
0,24 |
||||||||
3 |
|
6,1 |
|
3,5 |
1,0 |
4,0 |
|
|||||||
|
|
|
5 |
0,26 |
0,25 |
0,29 |
||||||||
4 |
|
3,1 |
|
2,0 |
2,5 |
0,9 |
|
6 |
0,30 |
0,31 |
0,89 |
12. Исследовалось четыре сплава, содержащие разное количе ство добавки А (с равным шагом). Прочность на разрыв оценива лась по четырем параллельным опытам (табл. 17).
№ |
У |
п/п |
1 |
|
+ |
_ |
|
+ |
+ |
|
+ |
27,64 |
2 |
|
|
+ |
|
+ |
+ |
|
+ |
8,85 |
3 |
|
— |
|
|
+ |
|
|
— |
20,00 |
4 |
|
+ |
+ |
|
+ |
+ |
|
+ |
12,75 |
5 |
|
+ |
— |
|
— |
+ |
|
+ |
13,18 |
6 |
|
|
+ |
|
+ |
|
— |
20,50 |
|
7 |
|
__ |
Л- |
__, |
|
|
+ |
16,50 |
|
8 |
|
+ |
— |
— |
|
17,45 |
|||
|
1 |
|
|
||||||
Таблица |
20. |
Результаты |
измерений плотности |
и состава |
продукта |
|
|||
N°. |
|
|
V |
|
схл |
№ |
V |
СКЛ |
|
п/п |
|
|
|
п/п |
|
||||
1 |
|
|
1,611 |
|
80,1 |
9 |
1,617 |
80,8 |
|
2 |
|
|
1,613 |
|
81,1 |
10 |
1,626 |
82,6 |
|
3 |
|
|
1,620 |
|
81,3 |
11 |
1,626 |
82,2 |
|
4 |
|
|
1,617 |
|
81,1 |
12 |
1,625 |
83,0 |
|
5 |
|
|
1,616 |
|
81,2 |
13 |
1,615 |
81,1 |
|
6 |
|
|
1,620 |
|
81,4 |
14 |
1,626 |
81,8 |
|
7 |
|
|
1,626 |
|
82,7 |
15 |
1,622 |
82,0 |
|
8 |
|
|
1,620 |
|
80,9 |
16 |
1,620 |
81,3 |
|
Таблица |
21. |
Результаты измерений расходов |
и плотности продукта |
|
|||||
№ |
|
|
|
|
У |
№ |
|
|
|
п/п |
|
|
|
|
п/п |
|
|
|
|
1 |
|
290 |
394 |
|
1,433 |
9 |
316 |
401 |
1,432 |
2 |
|
300 |
410 |
|
1,429 |
10 |
320 |
392 |
1,429 |
3 |
|
279 |
374 |
|
1,428 |
11 |
294 |
393 |
1,428 |
4 |
|
296 |
394 |
|
1,438 |
12 |
300 |
390 |
1,429 |
5 |
|
324 |
390 |
|
1,441 |
13 |
310 |
392 |
1,432 |
6 |
|
330 |
390 |
|
1,436 |
14 |
312 |
396 |
1,433 |
7 |
|
300 |
395 |
|
1,434 |
15 |
296 |
408 |
1,434 |
8 |
|
308 |
397 |
|
1,435 |
16 |
304 |
390 |
1,432 |
Определить значимость влияния добавки на прочность сплава» используя метод дисперсионного анализа,
13. Для оценки влияния микропримеси вещества А на продукт В, качество которого измеряется, был поставлен эксперимент по схеме дисперсионного анализа. Примесь изменялась с шагом 0,5% на 6 уровнях. Производилось три параллельных опыта (табл. 18).
Определить значимо ли влияние микропримеси на качество про дукта.
14. Отсеивающий эксперимент был поставлен на установке кри сталлизации фармакопейной трихлоруксусной кислоты (ТХУК) [60]. В отсеивающий эксперимент были включены следующие фак торы:
хх — скорость мешалки, об!мин, |
|
|
|
|
||
х.2 — скорость охлаждения раствора, град/ч, |
|
|
|
|||
х3 — время кристаллизации, |
ч, |
|
|
|
|
|
х4 — начальная концентрация ТХУК в растворе, %, |
|
|
|
|||
*s — температура |
кристаллизации, °С. |
|
|
|
||
Переменная |
состояния |
процесса — весовой |
выход |
продукта. |
||
Матрица эксперимента |
представлена в табл. |
19. |
Уровни факторов |
|||
распределялись |
случайным |
образом. |
|
|
|
|
Выделить наиболее |
значимые факторы |
методами |
случайного |
|||
баланса. |
|
|
|
|
|
|
15. Плотность продукта хлорирования этилового спирта заме ряется прибором. С целью выяснения возможности использования плотномера для измерения состава продукта были замерены его плотность и содержание хлораля в продукте. Результаты измере ний представлены в табл. 20.
Определить корреляцию между переменными и оценить ее зна чимость.
16.Плотность продукта хлорирования этилового спирта зави сит от расходов спирта хг и хлора х2. Данные измерений этих пере менных приведены в табл. 21.
Оценить связь между переменными, рассчитав коэффициенты корреляции.
17.Определялась зависимость растворимости тиосульфита нат рия у от температуры. Экспериментальные данные приведены ниже:
U °С |
Ю |
20 |
30 |
40 |
50 |
60 |
70 |
80 |
о/0 |
37,0 |
41,2 |
46,1 |
50,0 |
52,0 |
56,3 |
64,3 |
69,9 |
Определить коэффициент корреляции между переменными.
Гла ва III. ОСНОВНОЙ ЭКСПЕРИМЕНТ, ПЛАНЫ ПЕРВОГО ПОРЯДКА
Задача основного эксперимента — получение матема тической модели исследуемого объекта, которая используется для оптимизации объекта исследования (экстремальный эксперимент) или для целей аппроксимации. Чтобы получить математическую модель, используют факторный эксперимент, суть которого заклю чается в варьировании всех факторов объекта исследования по опре деленному плану. Построение планов факторного эксперимента может быть различным. Рассмотрим один из вариантов.
Таблица |
22. |
Матрица планирования |
Предположим, |
что объектом |
|||
эксперимента для двух факторов на |
исследования является |
реактор, |
|||||
двух уровнях |
|
|
|
||||
|
|
|
в котором выход продукта у за |
||||
|
|
Планирование |
Перемен |
||||
Опыты |
|
висит от температуры х1 и давле |
|||||
|
|
|
ная состоя |
||||
|
|
|
хг |
ния у |
ния х2 в реакторе. Дополнитель |
||
|
|
|
|
|
но известно, что изменение тем |
||
1 |
+ i |
+ i |
+ 1 |
У\ |
пературы от 60 до 80° С и давле |
||
2 |
+ i |
— 1 |
+ 1 |
Уч |
ния от 1 до 1,5 атм изменяет |
||
3 |
+ i |
+ i |
—1 |
Уз |
выход продукта. Обозначим мак |
||
4 |
+ i |
—1 |
—1 |
У* |
симальные и минимальные значе |
||
|
|
|
|
|
ния факторов хх и х2 символами |
||
комбинации факторов |
|
+ 1 и —1. Тогда все возможные |
|||||
при варьировании на двух уровнях (мини |
|||||||
мальном и |
максимальном) будут |
определены четырьмя |
опытами. |
||||
Такой план эксперимента принято записывать в виде |
матрицы пла |
||||||
нирования |
(табл. |
22). |
|
|
|
|
Во второй графе таблицы приведены значения фиктивной пере менной х0 (тождественно равной +1), которая понадобится при вы числении свободного члена полинома. В первой строке таблицы спла нирован первый опыт, когда факторам хг и х2 придают максимальные значения; во второй строке — когда фактору хх придают ми нимальное значение, а фактору х2 — максимальное, и т. д. Ока зывается подобное планирование имеет ряд достоинств и поэто му широко применяется для получения моделей. Например, поль
зуясь планом — табл. 22, |
можно |
после проведения |
эксперимента |
определить коэффициенты линейного уравнения .регрессии |
|||
У = |
К + |
+ b2x2. |
(S7) |
В главе III будут рассмотрены различные планы получения ли нейных моделей объектов исследования, которые могут служить основой для поиска области оптимума.
§ 1. Полный факторный эксперимент первого порядка
Сущность факторного эксперимента первого порядка состоит в одновременном варьировании всех факторов при его про ведении по определенному плану, представлении математической модели (функции отклика) в виде линейного полинома и исследова нии последнего методами математической статистики.
Теперь, когда принцип построения плана факторного экспери мента разъяснен, введем несколько определений.
Уровнем фактора называют определенное значение фактора, которое будет фиксироваться при проведении эксперимента. В пре дыдущем примере уровнями факторов будут 60 и 80° С для фактора «температура», а также 1 и 1,5 атм— для фактора «давление». Очевидно, уровнем факторов можно назвать и средние значения рассматриваемых интервалов, т. е. 70° С и 1,25 атм. Эти значения
факторов называются нулевыми уровнями, они определяют неко торую точку факторного про странства, которая в предвари тельном эксперименте была оце нена наилучшей по максимуму (или по минимуму) переменной состояния. Обозначим нулевой уровень i-го фактора, выражен ного в натуральных единицах (в данном примере в °С и атм), через Х/0. Введем еще одно по нятие — интервал варьирования. Это такое значение фактора в натуральных единицах, прибав ление которого к нулевому уров
ню дает верхний, а вычитание — нижний уровень фактора. Обозначим его АХ,.
Наконец, экстремальные значения, которые могут принимать факторы, не меняя своих физико-химических свойств и не искажая сути исследуемого процесса, назовем границами существования факторов, а интервал (Ximax — Ximjn^ — областью определения
факторов (область L на рис. 8). Очевидно, что интервал варьирова ния факторов должен составлять часть о'бласти определения фак торов, если решается задача оптимизации. Это необходимо для того, чтобы осуществить движение к оптимуму в области определения факторов. На рис. 8 область проведения эксперимента обозначена буквой М. В задачах же аппроксимации (или интерполяции) ин тервал варьирования охватывает всю описываемую область, т. е. для двухфакторной задачи верхними уровнями факторов Х х и Х2 являются Ximax, Х2тах, а нижними уровнями — Ximln и Х2т1п.
Тогда область L можно назвать интерполяционной, область М — областью постановки экстремального эксперимента.
Из определений следует, что областей М может быть несколько (в общем случае конечное множество). Можно также предположить несколько областей оптимума. Область определения факторов для данной задачи исследования одна. Обозначение верхних и нижних уровней факторов символами «+1», «—Ь> фактически соответствует
кодированию факторов по формуле
Xiщ—XiQ
*t = — s x r - |
(57а) |
Для рассмотренного примера (табл. 22) кодированные значения факторов (верхние и нижние уровни) следующие:
XlB |
80 — 70 |
_ , |
и |
v |
_ |
60 — 70 |
— 1; |
||
1Л |
— +I 1 ; |
х1н |
ю |
= |
|||||
|
10 |
|
|
|
|
|
|
|
|
П-- |
1,5— 1,25 |
+ |
1; |
|
1, 0 - 1,25 |
= — 1. |
|||
0,25 |
= |
л2н |
|
0,25 |
|
Рис. 9. Геометрическая интерпретация плана 22 |
Рис. 10. Геометрическая ин |
на плоскости: |
терпретация плана З2. |
а — в натуральных координатах; б — в кодирован |
|
ной форме. |
|
Кодирование факторов, по сути, означает переход от системы координат в натуральных единицах (рис. 9, а) к системе координат в кодированной форме (рис. 9, б). Каждая точка факторного про странства — (+1, +1), (—1, +1), (+1, —1), (—1, —1) — это опыт
висследованиях.
Вобщем случае эксперимент, в котором реализуются все воз можные сочетания уровней факторов, называетея-ншншм фактор
ным экспериментом |
(ПФЭ). Если каждый фактор варьируется на |
двух уровнях, то |
получается ПФЭ типа 2П Для двух факторов |
(п = 2) число опытов N = 22 = 4, что видно из табл. 23 и рис. 9. |
|
Можно осуществлять планирование эксперимента на трех уров |
нях (верхний, средний, нижний), тогда ПФЭ будет типа Зп и |
для |
п = 2 общее число опытов будет N = З2 = 9. Это положение |
ил |
люстрируется рис. 10.
Выбор уровней и интервалов варьирования факторов. В некото рых книгах по планированию эксперимента [1; 3] этот этап выде ляют как этап принятия решений перед составлением плана экспе римента. Действительно, построение плана эксперимента начинают с выбора определяющих его характеристик. Обычно первой рас сматривают область определения факторов. Ранее уже упоминалось (см. гл. II), что область определения факторов фиксируется в пред варительном эксперименте. Для этого используются результаты опытов и теоретические представления о процессе.
Далее из области определения факторов выбором нулевых уров ней и интервалов варьирования факторов выделяется часть области для планирования эксперимента (область М, рис. 8). Правильный выбор нулевых уровней (центра эксперимента) и интервалов варьи рования факторов имеет решающее значение для действенности математической модели.
Идеальным случаем при выборе нулевых уровней факторов яв ляется «попадание» центра эксперимента в область оптимальных
значений переменной состояния. Но такое выгодное обстоятельство возможно лишь при очень высоком уровне априорной информации, на который трудно рассчитывать при современном темпе исследова ния и внедрения технологических процессов. Поэтому если имеет ся некоторый опыт управления объектом исследования, можно при нять в качестве нулевых уровней те величины факторов, которые дали наилучшее значение переменной состояния. Такой подход таит в себе опасность получения лишь локального оптимума при нескольких экстремумах функции отклика, хотя вряд ли могут существовать априорные сведения об их наличии. Огромная роль при изучении неисследовавшихся ранее объектов с целью полу чения их математических моделей принадлежит аналогии и ин туиции.
Основное требование к интервалу варьирования состоит в том, чтобы он превышал удвоенную квадратическую ошибку фактора:
2s < |
Д Х с (Х - |
— X. . ), |
Г58) |
|
*1 |
1 |
‘ max |
‘ min' |
l 0 6 / |
где sx. — среднеквадратическое |
отклонение фактора |
X , ; А Х , — |
||
интервал варьирования; |
(Х,тах — Xtmin) — область |
определения |
фактора.
Это требование связано с тем, что интервал, между двумя сосед ними уровнями должен значимо (неслучайно) влиять на перемен ную состояния. Обычно интервал варьирования выбирается на ос новании априорной информации (или интуитивно) и затем уточ няется (если он выбран неудачно) после получения математической модели. Цена уточнения ощутима, так как повторение эксперимента резко увеличивает число опытов. Удачный выбор интервала варьи рования факторов гарантирует получение достоверной математиче ской модели объекта.
Как упоминалось ранее, определенные сведения о нулевых уров нях и интервалах варьирования удается получить на этапе предва рительного эксперимента.
Пример 1. Воспользуемся примером выбора области определе ния интервалов варьирования и нулевых уровней факторов, приве денных в работе [3, с. 82].
Изучалось ионообменное разделение смесей группы редкозе мельных элементов растворами иминодиуксусной кислоты. Пере менная состояния — содержание (в %) неодима в выходном рас творе. Предварительный эксперимент выделил всего два фактора — концентрацию (в вес. %) входного раствора Х г и pH раствора Х2. Область определения фактора строилась из следующих сообра
жений.
Известно, что при Хг > 3 работать нельзя, так как это предел растворимости данного вещества при нормальной температуре. Таким образом, X i ^ = 3. При выборе нижней границы области
определения фактора учитывалось то, чтсГ чем ниже концентрация, тем дольше идет процесс. При XimIn = 0,5 время протекания
процесса находится еще в допустимых пределах; дальнейшее сн.ц„ жение его уже нецелесообразно.
При выборе области определения Х г исходили из теоретического положения, что ионообменное разделение происходит благодаря одновременному присутствию в системе двух соединений: моно- ди-комплексов. Предварительный эксперимент показал, что прц
pH < |
3 |
кислота находится в недиссоциированном состоянии, а прц |
|
pH > |
8 |
оба |
соединения разрушаются. Следовательно, Х 2 |
= 8. |
|
“ |
3- |
В качестве нулевых уровней были приняты значения Xi0 = 1,5 |
|||
X ^ = |
7. |
В точке факторного пространства с такими координатам^ |
был получен наилучший результат предварительного эксперимент та. Важно также то, что она лежит внутри области определение факторов.
Изучая результаты предварительных опытов, экспериментатор пришел к выводу о том, что:
1) точность фиксирования факторов средняя (по результата^ ряда опытов);
2)поверхность отклика линейная (по однофакторным эксперт ментам);
3)диапазон изменения переменной состояния небольшой.
Тогда экспериментатор принимает решение — выбрать широ. кий (до 20% от области определения) интервал варьирования, чтобы его изменение было заметно по изменению переменной состояния: ДХ, = 0,5, ДХ2 = 1,0.
Построение матрицы планирования. План, содержащий запись всех комбинаций факторов или их части в кодированной форме, называется матрицей планирования. Табл. 22, например, является матрицей планирования для двух факторов на двух уровнях. Для построения матрицы планирования с большим числом факторов используют ряд приемов. Первый из них состоит в том, что элемен тарное сочетание первого фактора (+1, —1) повторяется для каж дого следующего фактора на верхнем и нижнем уровнях (табл. 23). Этот прием распространяется на построение матриц любой размернрсти. Столбец х0 — это столбец значений фиктивной переменной. Доказано, что его участие в матрице планирования делает расчеты коэффициентов математической модели более общими.
Иногда указанный прием построения матриц планирования трак туется как прием чередования знаков. Действительно, в первом
столбце знаки не |
меняются, |
во втором — меняются поочередно, |
в третьем — они |
чередуются |
через два, в четвертом — через 4 |
и т. д. (по показателям степеней двойки).
Свойства матрицы планирования. Рассмотренные матрицы пла нирования обладают такими свойствами, которые позволяют счи тать, что их построение выполнялось оптимально с точки зрения получаемой по результатам реализации матрицы планирования математической модели. Если мы ищем модель в виде уравнения
Таблица 23. Организация матриц планирования ПФЭ от 2а до 26
ЛЪ п/п |
|
1ИП эксперимента |
|
|
Ф ак торы |
|
|
|
|
|
|
|
|
|
|
||
|
|
|
*0 |
|
** |
|
*4 |
|
1 |
|
£ |
+ i |
+ i |
+ 1 |
+ 1 |
+ 1 |
+ 1 |
2 |
|
(Г) |
+ i |
|
+ 1 |
+ 1 |
- и |
+ 1 |
|
е |
|
||||||
3 |
|
с |
+ i |
+ i |
— 1 |
+ 1 |
+ 1 |
+ 1 |
4 |
|
ся |
+ 1 |
— 1 |
— 1 |
+ 1 |
+ 1 |
+ 1 |
|
|
(Г) |
|
|
||||
|
|
|
|
|
|
|
|
|
5 |
|
е |
+ i |
+ i |
+ 1 |
—1 |
+ 1 |
+ 1 |
|
с |
|||||||
6 |
|
|
+ i |
— 1 |
+ 1 |
—1 |
+ 1 |
+ 1 |
7 |
|
О) |
+ i |
+ i |
— 1 |
—1 |
+ 1 |
+ 1 |
8 |
|
+ i |
— 1 |
— 1 |
—1 |
+ 1 |
+ 1 |
|
|
CD |
|||||||
|
|
О |
|
|
|
-И |
|
|
9 |
|
а |
+ i |
+ i |
+ 1 |
— 1 |
+ 1 |
|
10 |
|
|
+ i |
— 1 |
+ 1 |
+ 1 |
— 1 |
+ 1 |
И |
|
|
+ 1 |
+ i |
— 1 |
+ 1 |
— 1 |
+ 1 |
12 |
ыО |
|
+ i |
— I |
— 1 |
+ 1 |
— 1 |
+ 1 |
13 |
|
+ i |
+ 1 |
+ 1 |
— 1 |
— 1 |
+ 1 |
|
14 |
<М |
|
+ i |
— 1 |
+ 1 |
— 1 |
— 1 |
+ 1 |
О) |
|
|||||||
15 |
е |
|
+ i |
-и |
— 1 |
— 1 |
—1 |
+1 |
16 |
с |
|
+ | |
—! |
—1 |
—1 |
—1 |
+1 |
|
|
|||||||
17 |
|
|
+ i |
+ 1 |
+1 |
+1 |
+1 |
—1 |
18 |
|
|
+1 |
—1 |
+1 |
+1 |
+1 |
— 1 |
19 |
|
|
+ i |
-Т-1 —1 -Т-1 |
+1 |
—1 |
||
20 |
|
|
+1 |
—1 |
—1 |
+1 |
+1 |
—1 |
21 |
|
|
+ i |
+ 1 |
+1 |
—1 |
+1 |
—1 |
22 |
|
|
+ i |
—1 |
+1 |
—1 |
+1 |
—1 |
23 |
|
|
-и |
+1 |
—1 |
—1 |
+ 1 |
—1 |
24 |
|
|
—1 |
—1 —1 |
+1 |
—1 |
||
25 |
|
|
+ i |
+1 |
+1 |
+1 |
—1 |
—1 |
26 |
|
|
+1 |
—1 |
+1 |
+1 |
—1 |
— 1 |
27 |
|
|
+ i |
+1 |
—1 |
+1 |
—1 |
—1 |
28 |
|
|
+ i |
—1 |
—I |
+1 |
—1 |
—1 |
29 |
|
|
+ | |
+1 |
+1 |
—1 |
—1 |
—1 |
30 |
|
|
+ i |
^1 |
+1 |
—1 |
— 1 |
—1 |
31 |
|
|
+ i |
+1 |
—1 |
—1 —1 —1 |
||
32 |
|
|
+ | |
—1 —1 —1 —1 —1 |
регрессии, то коэффициенты должны быть наилучшим^ и точность предсказания значений переменной состояния одинакова в любом направлении факторного пространства. Эти требования формули руются как условия ортогональности и рототабельности.
Рассмотрим эти условия подробнее. Вначале заметим, что из построения матрицы планирования вытекают два следствия, кото рые математически можно представить так:
N |
|
|
|
и2= 1 хш = 0; |
|
(59) |
|
N |
|
|
|
2 xlu = N (i = 1.2........... |
п; и = 1,2............. |
N), |
(60) |
U = I
где п — число факторов; N — число опытов (или строк матрицы планирования).
Эти условия легко проверить по табл. 23. Называются они со ответственно свойство симметричности и свойство нормировки.
Условие ортогональности предполагает равенство нулю суммы почленных произведений любых двух столбцов матрицы:
N |
|
|
2 XiuXju — 0 (£, / = 1, 2, |
/2, i "=£=■/). |
(61) |
и=1 |
|
|
Это условие также легко проверить по табл. 23. Действительно, полный факторный эксперимент типа 2п является ортогональным.
же |
Ортогональные планы ПФЭ (для линейных моделей) имеют так |
еще одно свойство— рототабельность. Последнее предполага |
|
ет |
равенство и минимальность дисперсий предсказанных значений |
переменной состояния для всех точек факторного пространства. По закону накопления ошибок можно записать для дисперсии пред сказанных уравнением регрессии значений переменной состояния:
s~ = |
4> + |
4,* 1+ |
+ |
, |
<б2> |
где si. — дисперсия |
коэффициентов |
модели |
Ь(. |
|
|
Из условия (61) вытекает, что дисперсии коэффициентов регрес |
|||||
сии равны между собой. Тогда можно записать: |
|
||||
|
п |
х2 — р2, где р радиус сферы: |
|
||
или с учетом того, что 2 |
|
||||
|
1=1 |
|
|
|
|
|
s~ = Sty (1 + |
р2)- |
|
(64) |
Отсюда ясно, что дисперсия предсказанного значения выход ной переменной зависит только от радиуса сферы. Это свойство рототабельности эквивалентно независимости дисперсии выходной переменной от вращения координат в центре плана и оправдано при поиске оптимума градиентными методами. Ясно, что исследо вателю удобно иметь дело с такой информацией, содержащейся в уравнении регрессии, которая равномерно распределена по сфере с радиусом р2. Действительно, такое положение можно признать разумным, ибо с помощью уравнения регрессии будут предприни маться попытки предсказать положение еще неизвестных участков факторного пространства. Равноценность этих участков в смысле ошибки предсказания, по-видимому, является необходимой.
Основное преимущество факторного эксперимента заключается в том, что в эксперименте варьируются одновременно все факторы. Это приводит к тому, что дисперсия в оценке коэффициентов регрес сии оказывается в N раз меньше ошибки опыта.
При классическом подходе эксперименты ставятся в определен ной последовательности: все факторы фиксируются на некотором