метод планирования эксперимента в исследовательских работах

advertisement
ОСНОВА И ПОТЕНЦИАЛ КАФЕДРЫ
МЕТОД ПЛАНИРОВАНИЯ ЭКСПЕРИМЕНТА
В ИССЛЕДОВАТЕЛЬСКИХ РАБОТАХ
М.И. СТАЛЬНАЯ, Т.А. МОШКИНА, К.А. БЕЖАНОВ
Метод планирования эксперимента был
разработан в 1951 году Боксом и Уилсоном в
Англии.
Экспериментатору
предлагается
ставить последовательные небольшие серии
опытов, в каждой из которой одновременно
варьируются по определённым правилам все
факторы. Серии организуются таким образом,
чтобы после математической обработки
предыдущей можно было выбрать условие
проведения
следующей
серии.
Так
последовательно
достигается
область
оптимума.
Судя по этим опытам можно судить об
оценке эффективности функционирования
системы, но важно помнить, что речь идёт о
системе в целом. Часто система состоит из
ряда подсистем, каждая из которых может
оцениваться своим локальным параметром
оптимизации. Путь к единому параметру
оптимизации часто лежит через обобщение. Из
многих откликов, определяющих объект, очень
часто трудно выбрать один, самый важный.
Одним из наиболее удобных способов
построения обобщенного отклика является
обобщенная
функция
желательности
Харрингтона. В основе построения этой
обобщенной
функции
лежит
идея
преобразования
натуральных
значений
частных откликов в безразмерную шкалу
желательности или предпочтительности.
Чтобы получить шкалу желательности,
удобно
пользоваться
готовыми
разработанными
таблицами
соответствий
предпочтения
в
между
отношениями
эмпирической и числовой (психологической)
системах (табл. 1).
После
того,
как
выбрана
шкала
желательности,
и
частные
отклики
преобразованы
в
частные
функции
желательности, можно приступить к основной
задаче – построению обобщенного показателя
D, названного Харрингтоном обобщенной
функцией желательности. Обобщать, то есть
переходить от di к D, предлагается по формуле
(1).
D=n
n
∏d
u
.
(1)
u =1
Способ задания обобщенной функций,
желательности таков, что если хотя бы одна
ПОЛЗУНОВСКИЙ АЛЬМАНАХ №3 2004
частная желательность du = 0, то обобщенная
функция тоже будет равна
нулю, с другой
стороны D = 1 тогда и только тогда, когда все du =
l (u = 1, 2, …, n). Обобщенная функция
желательности весьма чувствительна к малым
значениям частных желательностей.
После того как выбран объект исследования
и параметр исследования, нужно включить в
рассмотрение все существенные факторы,
которые могут влиять на процесс.
Наиболее подходящим в этом случае
является планирование эксперимента [1], в
котором при выборе области эксперимента,
прежде всего надо оценить границы областей
определения факторов. После чего выбираются
интервалы варьирования факторов.
Интервалом
варьирования
факторов
называется некоторое число, прибавление
которого к основному уровню даёт верхний, а
вычитание – нижний уровни фактора.
Для
упрощения
записи
условий
эксперимента и обработки экспериментальных
данных масштабы по осям выбираются так,
чтобы верхний уровень соответствовал +1,
нижний –1, а основной – нулю. Для факторов с
непрерывной областью определения это всегда
можно сделать с помощью преобразования:
xj =
~
xj −~
x j0
Ij
,
(2)
где: x j – координатное значение фактора;
~
x j – натуральное значение фактора;
~
x j0 – натуральное значение основного
уровня;
I j – интервал варьирования;
j – номер фактора.
Запись матрицы планирования, особенно
для многих факторов, громоздка. Для её
сокращения удобно ввести условные буквенные
сокращения строк. Это делается следующим
образом. Порядковый номер фактора ставится в
соответствие
строчной
букве
латинского
алфавита: x1 – a, x2 – b и т. д. Если теперь для
строки
матрицы
планирования
выписать
латинские
буквы
только
для
факторов,
находящихся на верхних уровнях, то условия
опыта будут заданы однозначно. Опыт со всеми
43
Т.А. МОШКИНА, К.А. БЕЖАНОВ
факторами на нижних уровнях условимся
обозначать (1). Пример матрицы планирования
вместе
с
принятыми
эксперимента
22
буквенными обозначениями приведена в табл.
2.
После построения матрицы планирования
и проведения эксперимента необходимо найти
по его результатам значения неизвестных
коэффициентов модели. Можно утверждать,
что эксперимент проводится для проверки
гипотезы о том, что линейная модель η = β0 +
β1x1 + β2x2 адекватна. Греческие буквы
используются для обозначения «истинных»
генеральных
значений
соответствующих
неизвестных.
Эксперимент,
содержащий
конечное число опытов, позволяет только
получить
выборочные
оценки
для
коэффициентов уравнения y = b0 + b1x1 + … +
bkxk. Их точность и надёжность зависят от
свойств выборки и нуждаются в статистической
проверке. Коэффициенты можно вычислять по
формуле:
N
bj =
∑x
i =1
ji
N
yi
,
(3)
где j = 0, 1, …, k.
Количество опытов в полном факторном
эксперименте значительно превосходит число
определяемых
коэффициентов
линейной
модели, то есть он обладает большой
избыточностью опытов. Сократить их число
возможно за счёт той информации, которая не
очень существенна при построении линейной
модели. При этом нужно стремиться, чтобы
матрица планирования не лишилась своих
оптимальных свойств. Это лежит в основе
дробного факторного эксперимента.
К примеру, пользуясь таким планированием
для эксперимента 22, можно вычислить четыре
коэффициента
и
представить
результаты
эксперимента в виде неполного квадратичного
уравнения, показанного в формуле:
y=b0x0+b1x1+b2x2+b12x1x2.
(4)
Если имеется основание считать, что в
выбранных интервалах варьирования процесс
может быть описан линейной моделью, то
достаточно определить три коэффициента b0, b1
и b2. Остаётся одна степень свободы. Употребим
её для минимизации числа опытов. При
линейном приближении b12 → 0 и вектор-столбец
x1x2 можно использовать для нового фактора x3.
Здесь уже не будет раздельных оценок, которые
мы имели в полном факторном эксперименте 2k.
Оценки смешиваются следующим образом: b1 →
β1 + β23, b2 → β2 + β13, b3 → β3 + β12. Так как
постулируется линейная модель, все парные
взаимодействия незначимы. То есть, чтобы
сократить число опытов, необходимо новому
фактору присвоить вектор-столбец матрицы,
принадлежащий
взаимодействию,
которым
можно пренебречь. Тогда значение нового
фактора в условиях опытов определяется
знаками этого столбца. Поставив четыре опыта
для оценки влияния трёх факторов, мы
воспользовались половиной полного факторного
эксперимента 23 или «полурепликой», но с её
помощью смогли оценить эксперимент 23 (в
полном факторном эксперименте понадобилось
бы ставить восемь опытов).
Таким образом, при применении метода
планирования эксперимента целью является
нахождение
аналитической
зависимости,
определяющей взаимосвязь между фактором и
исследуемым параметром. На основе найденного
уравнения можно производить не только
оптимизацию исследуемого процесса, но и
прогнозировать поведение этого процесса при
различных изменяющихся условиях.
Таблица 1
Стандартные отметки на шкале желательности
Желательность
Очень хорошо
Хорошо
Удовлетворительно
Плохо
Очень плохо
44
Отметки на шкале желательности
1 – 0,8
0,8 – 0,63
0,63 – 0,37
0,37 – 0,2
0,2 – 0
Таблица 2
ПОЛЗУНОВСКИЙ АЛЬМАНАХ №3 2004
МЕТОД ПЛАНИРОВАНИЯ ЭКСПЕРИМЕНТА В ИССЛЕДОВАТЕЛЬСКИХ РАБОТАХ
Матрица планирования эксперимента 22 с принятыми буквенными обозначениями
Номер
опыта
1
2
3
4
x1
x2
Буквенные обозначения строк
y
–1
+1
–1
+1
–1
–1
+1
+1
(1)
a
b
ab
y1
y2
y3
y4
СПИСОК ЛИТЕРАТУРЫ
1.
Адлер Ю.П., Маркова Е.В., Грановский Ю.В.
Планирование эксперимента при поиске
ПОЛЗУНОВСКИЙ АЛЬМАНАХ №3 2004
оптимальных условий. – М.: Наука, 1976. –
280 с.
45
Download