При проведении опытных исследований различают пассивный и активный эксперимент Обработка результатов пассивного эксперимента проводится методами регрессионного и корреляционного анализа, и выбор вида эмпирической модели (уравнения регрессии), т.е. решение задачи структурной идентификации, является достаточно сложной задачей Для решения этой задачи для одной входной переменной х предложены эффективные методы, в которых предусматривается преобразование системы координат как для входной x, так и для выходной переменной у При большем числе входных переменных х 1..., х т надёжных методов определения вида уравнения регрессии в настоящее время не существует

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

Уравнения регрессии имеют следующий вид: вдали от экстремального значения выходной переменной у: вблизи экстремального значения выходной переменной у («в почти стационарной области»):

Они включают слагаемые с двойным взаимодействием входных переменных и не учитывают взаимодействия более высоких порядков (тройные, четверные и т.д.), вероятность которых существенно меньше Последнее уравнение включает слагаемые с квадратами входных переменных коэффициенты которого получаются при обработке результатов активных экспериментов II-го порядка Первое уравнение не включает слагаемые с квадратами входных переменных, и его коэффициенты получаются при обработке результатов активных экспериментов 1-го порядка например, ПФЭ полного факторного эксперимента Активный эксперимент планируется таким образом, чтобы упростить обработку его результатов методами регрессионного и корреляционного анализа

Метод полного факторного эксперимента (ПФЭ) служит для получения математического описания процесса в виде отрезка ряда Тейлора, содержащего линейные члены и парные взаимодействия переменных величин: Поверхность отклика

X 2 X 1

Поверхность отклика

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

Для упрощения записи условий эксперимента и обработки экспериментальных данных пользуются кодированными переменными Переход к кодированным переменным осуществляется по следующей формуле: где x j кодированное значение фактора; натуральное значение фактора; натуральное значение основного уровня; I j интервал варьирования; jномер фактора.

Пример

Общее число опытов ПФЭ определяется по формуле где k число факторов План проведения экспериментов, называется матрицей планирования Она может быть составлена в натуральных значениях соответствующих уровней факторов (этот вариант матрицы планирования необходим при постановке экспериментальных исследований) или в безразмерных (кодированных) значениях переменных (этот вариант матрицы используется главным образом для расчетов коэффициентов уравнения регрессии)

Матрица планирования эксперимента 2 2 x2x2 x1x Номер опытаx1x1 x2x2 y 1120ylyl 25 у 2 у узус 45 у 4 у 4 Матрица планирования эксперимента 2 2 Номер опытаx1x1 x2x2 y 1 ylyl 2+1 у 2 у узус у 4 у 4 ~~

Матрица планирования эксперимента 2 3 Номер опытаx1x1 x2x2 x3x3 y 1 ylyl 2+1 у 2 у узус у 4 у У5У у 6 у у 7 у у 8 у 8

Матрица планирования ПФЭ обладает следующими свойствами: где i номер опыта; j номер фактора Свойство, выраженное уравнением, называется ортогональностью матрицы Оно позволяет вычислять коэффициенты регрессии по простым формулам независимо друг от друга

Расчет коэффициентов регрессии для ПФЭ производится методом наименьших квадратов Благодаря свойствам матрицы полного факторного эксперимента получаются достаточно простые формулы для расчёта коэффициентов уравнения регрессии

Проверка гипотезы о значимости коэффициентов уравнения регрессии производится с помощью критерия Стьюдента Проверка адекватности уравнения регрессии экспериментальным данным проводится с помощью критерия Фишера

Пример Матрица планирования и результаты трёхфакторного эксперимента

Метод дробного факторного эксперимента (ДФЭ) предназначен для получения математического описания процесса в виде отрезка ряда Тейлора, содержащего линейные члены, а также в ряде случаев и парные взаимодействия переменных величин: Коэффициенты уравнения определяются по экспериментальным данным методом наименьших квадратов

Матрица дробного факторного эксперимента Одним из недостатков полного факторного эксперимента (ПФЭ) является то, что с увеличением количества факторов резко возрастает число опытов полного факторного эксперимента ( N=2 k ) Например: N=2 2 =4; N=2 5 =32; N=2 6 =64; N=2 15 =32768 Причем количество опытов ПФЭ значительно превосходит число определяемых коэффициентов уравнения регрессии Дробным факторным экспериментом называется система опытов, представляющая собой часть ПФЭ, позволяющая рассчитать коэффициенты уравнения регрессии и сократить объем экспериментальных работ Для нахождения математического описания процесса используется определенная часть ПФЭ 1/2, 1/4, 1/8 и так далее. Эти системы опытов называют дробными репликами

Возможные дробные реплики от ПФЭ типа 2 3 1/2 1/4 1/2 1/4 Предположим, что надо исследовать влияние на результаты химико-технологического процесса трех факторов и получить его математическое описание в виде линейного уравнения

Возьмем матрицу полного двухфакторного эксперимента и приравняем произведение х 1 х 2 к фактору х 3 Планирование типа Номер опыта х 0 х 0 х 1 х 1 х 2 х 2 х 1 х 2=х 1 х 2= х 3 х 3 х 1 х 3 х 2 х 3 х 2 х 3 х 1 х 2 х 3 х 1 х 2 х 3 Функция отклика ylyl у 2 у узус у 4 у 4 По данному плану мы можем определить коэффициенты регрессии b 0, b 1, b 2, b 3, однако, они будут смешаны с парными и тройными взаимодействиями:

Таким образом, сокращение числа опытов влечет за собой корреляцию между столбцами матрицы ДФЭ, что не позволяет раздельно оценивать эффекты факторов и эффекты взаимодействий В результате мы получаем так называемые совместные (смешанные) оценки Этот недостаток рассматриваемого плана является своеобразной «платой» за уменьшение общего числа опытов

Условия обозначения дробных реплик и число опытов Число факторов Дробная реплика Условное обозначение Число опытов для дробной- реплики для полного факторного эксперимента 3 1/2-реплика от /2-реплика от /4-реплика от /8-реплика от /16-реплика от /2-реплика от /4-реплика от /8-реплика от /16-реплика от /32-реплика от /64-реплика от /128-реплнна от /256-реплика от /512-реплика от /1024-реплика от /2048-реплика от

ПЛАНИРОВАНИЕ СО СМЕШИВАНИЕМ. ОПРЕДЕЛЕНИЕ СОВМЕСТНЫХ ОЦЕНОК КОЭФФИЦИЕНТОВ РЕГРЕССИИ Планирование эксперимента, при котором некоторые из факторов приравниваются к произведениям нескольких факторов, называется планированием со смешиванием Общее число опытов ДФЭ вычисляется по формуле N=2 k-p, где k общее число факторов р число факторов, приравненных к произведениям Существует правило, позволяющее определить, какие коэффициенты регрессии определяются совместно при планировании со смешиванием. Рассмотрим это правило на конкретном примере

Методом ДФЭ будем искать математическое описание процесса в виде уравнения регрессии Воспользуемся планированием типа и примем x 4 =-x 1 x 2 x 5 =x 1 x 2 x 3 Эти равенства называются генерирующими соотношениями Они показывают, какое из взаимодействий факторов принято незначимым, а поэтому заменено в матрице планирования новой независимой переменной Выбор генерирующих соотношений в общем случае произволен Умножив обе части генерирующих соотношений на х 4 и х 5, получим x 4 x 4 =1=-x 1 x 2 x 4 x 5 x 5 =1=x 1 x 2 x 3 x 5 Эти равенства называются определяющими контрастами

Перемножив почленное 1-й и 2-й определяющий контрасты, получим 1=(-x 1 x 2 x 4 )(x 1 x 2 x 3 x 5 ) 1=-x 3 x 4 x 5 Составим систему равенств из единицы и правых частей всех определяющих контрастов Получим выражение для обобщающего определяющего контраста Умножив фактор х 1 на обобщающий определяющий контраст, получим Отсюда следует, что коэффициент регрессии b 1 будет оценкой

Аналогично получим

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

Например, при выборе полу реплики типа возможны 8 вариантов решений Наибольшая разрешающая способность у реплик 7-й и 8-й При наличии априорной информации о значимости взаимодействий факторов можно разработать наилучшую систему смешивания оценок Если таких сведений нет, то выбирают реплику с наибольшей разрешающей способностью, так как тройные взаимодействия менее важны, чем двойные

РАСЧЕТ КОЭФФИЦИЕНТОВ РЕГРЕССИИ Вычисление коэффициентов регрессии при использовании ДФЭ производится методом наименьших квадратов по формулам, которые применяются и для обработки данных ПФЭ При расчете коэффициентов регрессии следует учесть, что подлежат вычислению только те коэффициенты при взаимодействиях факторов, столбцы уровней которых в матрице планирования не коррелированны со столбцами отдельных факторов Игнорирование этого правила ведет к нарушению ортогональности плана

РАСЧЕТ КОЭФФИЦИЕНТОВ РЕГРЕССИИ В общем случае формулы для расчета коэффициентов регрессии по результатам ДФЭ имеют вид

Пример Матрица планирования и результаты ДФЭ типа 2 5-2

Планы второго порядка используют в тех случаях, когда функция отклика существенно нелинейна и не может быть аппроксимирована линейным приближением в рассматриваемой области факторного пространства Это имеет место, например, в окрестности экстремума функции отклика Эту часть поверхности отклика принято называть «почти стационарной» областью Обычно исследователь обнаруживает «почти стационарную» область в результате оптимизации по методу крутого восхождения или симплекс- планирования Для адекватного математического описания этой области часто используют полином второго порядка

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

Исходя из этих соображений были разработаны так называемые композиционные планы второго порядка Структура данных планов представляет собой композицию из плана первого порядка и некоторого количества добавочных опытов При этом один или несколько опытов проводятся в центре плана Благодаря своей структуре такие планы экспериментов называются центральными композиционными (ЦКП)

Матрица двухфакторного композиционного плана второго порядка (0, - α) (0,α) (α,0) ( - α,0) ( 0,0) (-1, 1) (1, 1) (-1, -1) (1,-1) x2x2 x1x1 Номер опыта Фрагмент планаx1x1 x2x2 у 1Ядро плана+ 1 у 1 у у 2 у 2 3 у 3 у 3 4 у 4 у 4 5Опыты в «звездных» точках+α0 у 5 у 5 6- α0 у 6 у αу 7 у αу 8 у 8 9Опыт в центре плана 00 у 9 у 9

Матрица трехфакторного композиционного плана второго порядка

Если число факторов больше четырех, то в качестве ядра плана целесообразно использовать дробный факторный эксперимент Общее число опытов ЦКП рассчитывается по формуле N цкп =N я +N зв +N о где N я число опытов в ядре плана; N зв число опытов в «звездных» точках; N о число центральных опытов. Очевидно, N зв =2k, т. е. вдвое превышает количество факторов Из структуры ЦКП следует, что каждый фактор варьируется на пяти уровнях: -α, -1, 0, +1, +α. Из структуры ЦКП следует, что каждый фактор варьируется на пяти уровнях: -α, -1, 0, +1, +α.

Известны два вида центрального композиционного планирования - ортогональное - ротатабельное

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

Например, для ОЦКП с двумя факторами имеем Для обеспечения ортогональности всех столбцов матрицы планирования вместо квадратов значений факторов вводят новые переменные величины Из условия ортогональности

4α 4 +4α 2 N яN B (N 3B +N 0 ) = 0. получено уравнение для звездного плеча α: Решение этого уравнения имеет вид Основные характеристики ортогонального центрального композиционного планирования nNaNa N 3B NoNo Nα , , , ,547

Матрица ОЦКП для двух факторов Ном ер опыт а Фрагмен т плана х 1 х 1 х 2 х 2 х 1 х 2 х 1٭х 1٭ х 2٭х 2٭ y 1 Ядро +1 +0,33 y1y ,33 y2y ,33 y3y ,33 y4y ,33-0,67y5y5 6 « Звездные» 00+0,33-0,67y6y6 7 точки ,67+0,33y7y ,67+0,33y8y8 9Центр 000-0,67 y9y9

Матрица планирования трёхфакторного ОЦКП и результаты эксперимента