- •Моделирование химико-технологических процессов
- •Оглавление
- •Введение, основные понятия и определения
- •1. Методы моделирования и область их применения
- •1.1. Физическое моделирование (фм)
- •1.2. Математическое моделирование (мм)
- •2. Основные принципы и направления при построении и решении математических моделей
- •2.1. Составление математической модели
- •2.2. Нахождение решения математической модели
- •2.3. Проверка моделей на адекватность
- •3. Математическое описание структуры потока в аппарате – основа построения моделей
- •3.1. Методы исследования структуры потоков
- •3.2. Основные характеристики функции распределения потока по времени пребывания в аппарате
- •4. Типовые модели структуры потока
- •4.1. Модель идеального перемешивания
- •4.2. Модель идеального вытеснения
- •4.3. Ячеечная модель
- •4.4. Ячеечная модель с рециркуляцией
- •4.5. Диффузионная модель
- •4.6. Сравнение аппаратов соответствующих моделям ип и ив
- •5. Методы статистического анализа эксперимента
- •5.1. Основные характеристики случайных величин
- •5.2. Равномерное распределение
- •5.3. Нормальное распределение
- •5.4. Доверительные интервалы и доверительная вероятность, распределение Стьюдента
- •5.5. Определение общей дисперсии для серии параллельных опытов
- •5.6. Оценка дисперсии нормально распределенной случайной величины
- •5.7. Проверка однородности результатов измерений
- •6. Планирование эксперимента
- •6.1. Основные понятия и определения
- •6.2. Выбор области проведения эксперимента
- •6.3. Полный факторный эксперимент (пфэ)
- •6.4. Дробный факторный эксперимент
- •7. Оптимизация эксперимента
- •8. Планы второго порядка
- •8.1. Центральное композиционное планирование
- •8.2. Ортогональный план второго порядка
- •8.3. Ротатабельный план второго порядка
- •Список литературы
- •Приложение 1. Квантили нормального распределения
- •Приложение 2. Квантили распределения Пирсона
- •Приложение 3. Значения параметра для различных уровней значимости и степеней свободы
- •Приложение 4. Квантили распределения Стьюдента
- •Приложение 6. Пример использования модели ип для описания процесса непрерывной массовой кристаллизации
- •Приложение 7. Примеры использования типовых моделей для описания процесса теплообмена
- •Лицензия лр № 020370
6.4. Дробный факторный эксперимент
С увеличением числа факторов в соответствии с формулой N = 2k количество опытов полного факторного эксперимента резко возрастает. Однако если при получении уравнения можно ограничиться линейным приближением, то число опытов можно значительно уменьшить путем использования дробных реплик от полного факторного эксперимента. Полученный план в этом случае будет называться дробным факторным экспериментом (ДФЭ). Чтобы дробная реплика представляла собой ортогональный план, в качестве реплики следует брать полный факторный эксперимент для меньшего количества факторов. Число опытов при этом должно быть больше или равно числу неизвестных коэффициентов в уравнении регрессии. Допустим, что требуется получить линейное приближение небольшого участка поверхности отклика при трех независимых факторах:
= b0 + b1 x1 + b2 x2 + b3 x3 . (6.13) |
Для решения этой задачи можно ограничиться четырьмя опытами, если в планировании для ПФЭ 22 (табл.6.4) использовать столбец x1x2 в качестве плана для x3 (табл. 6.5).
Сокращенный план (табл. 6.5) называется полурепликой от ПФЭ 23. Пользуясь таким планированием, можно оценить свободный член и коэффициенты уравнения регрессии при линейных членах.
Таблица 6.4 |
|
Таблица 6.5 |
||||||||
ПФЭ 22 |
|
Полуреплика от ПФЭ 23 |
||||||||
№ оп. |
x0 |
x1 |
x2 |
x1x2 |
|
№ оп. |
x0 |
x1 |
x2 |
x3 |
1 |
+1 |
+1 |
+1 |
+1 |
|
1 |
+1 |
+1 |
+1 |
+1 |
2 |
+1 |
–1 |
–1 |
+1 |
|
2 |
+1 |
–1 |
–1 |
+1 |
3 |
+1 |
–1 |
+1 |
–1 |
|
3 |
+1 |
–1 |
+1 |
–1 |
4 |
+1 |
+1 |
–1 |
–1 |
|
4 |
+1 |
+1 |
–1 |
–1 |
Однако надо иметь в виду, что использовать столбец x3 вместо столбца x1x2 можно только при равенстве нулю эффектов взаимодействия (в рассматриваемом случае коэффициент b12). На практике не всегда удается установить, равны ли нулю эффекты взаимодействия, но часто имеются основания полагать, что некоторые из них малы по сравнению с линейными эффектами. Если коэффициенты регрессии при парных взаимодействиях не равны нулю, то полученные коэффициенты будут смешанными оценками для генеральных коэффициентов:
b1 1 + 23; b2 2 + 13; b3 3 + 12, (6.14) |
где – математические ожидания для соответствующих коэффициентов.
Эти генеральные коэффициенты не могут быть раздельно оценены по плану, включающему всего четыре опыта, так как при этом столбцы для линейных членов и парных коэффициентов одинаковы. Если, например, в дополнение к столбцам табл. 6.5 вычислить столбец для произведения х1х3, то окажется, что элементы этого столбца в точности равны элементам столбца х2. Таким образом, сокращение числа опытов приводит к получению смешанных оценок для коэффициентов. Чтобы определить, какие генеральные коэффициенты смешаны, удобно пользоваться следующим приемом: поставив х3 на место х1х2 (табл.6.5), получаем соотношение
x3 = x1 x2, (6.15) |
называемое генерирующим соотношением. Умножив обе части генерирующего соотношения на x3, получим единичный столбец
= x1 x2 x3 = 1. (6.16) |
Произведение x1x2x3 называется определяющим контрастом, при его помощи удобно определять, в каких столбцах одинаковые элементы. Умножив поочередно определяющий контраст на x1, x2 и x3, получим
х1 = = х2х3; х2 = х1х3; х3 = х1х2. (6.17) |
При использовании ДФЭ необходимо иметь четкое представление о так называемой разрешающей способности дробной реплики, т.е. определить заранее, какие коэффициенты являются несмешанными оценками для соответствующих генеральных коэффициентов. Тогда в зависимости от поставленной задачи подбирается дробная реплика, при помощи которой можно извлечь максимальную информацию из эксперимента. Необходимо выбирать такие реплики, у которых смешанные коэффициенты минимальны. При этом следует иметь в виду, что в реальных задачах взаимодействия большего порядка бывают равны нулю значительно чаще, чем меньшего, т.е. более вероятно, что нулю равно тройное, а не двойное взаимодействие. Таким образом, заменяя столбцы с минимальными смешанными коэффициентами новыми факторами, мы можем существенно уменьшить количество проводимых опытов при минимальной потере точности модели.