11
Общий вывод: не имеет смысла уменьшать С.И . более, чем до С.И . 0 / 3, а
увеличивать объем выборки n имеет смысл до тех пор, пока 
2С.И . 20 / n М или систематической составляющей погрешности СИ.
Некоторые практические рекомендации [3]:
начинающим исследователям не следует рассматривать разработанные в ТПЭ положения как догму, а видеть эту теорию как руководство к действию, как аналог, пример возможных решений, приглашающий к размышлению и принятию собственных эвристических или волевых решений;
иногда необходимо бывает построить план с более сложными, не рассматриваемыми здесь зависимостями. В этом случае необходимо помнить, что если несколько худшее решение можно получить сегодня, а оптимальное – неизвестно когда, то целесообразно исходить из положения «лучшее – враг хорошего» и использовать положения ТПЭ как примеры, достойные подражания для принятия эвристических решений в ситуациях, оптимальное решение которых еще неизвестно.
Определение оптимальных условий процессов. Стратегия поиска оптимума
Применение методов планирования эксперимента позволяет эффективно находить оптимальные условия и режимы протекания процессов. Чаще всего для этого используется градиентный метод крутого восхождения.
Постановка задачи нахождения оптимума
Многомерная оптимизация – это поиск экстремума функции нескольких аргументов (функции векторного аргумента).
Количество переменных (координат), определяющих функцию, будем обозначать k. Задача многомерной оптимизации имеет место при количестве факторов k ≥ 2.
Задача: найти экстремальное (минимальное или максимальное) значение отклика y и соответствующий ему векторный аргумент (факторы) x ={x1, x2…xk}:
y F(x1, x2 , ...,xk ) extr .
Функция k переменных задает гиперповерхность в (k+1)-мерном пространстве. Для функций одной или двух переменных легко получить графическую интерпретацию задачи (рис. 3) [3].
Метод градиентного поиска подразумевает движение в направлении наиболее быстрого изменения целевой функции у (отклика, рассчитанного по регрессионной
12
модели). Метод основан на последовательном нахождении вектора-градиента целевой функции, который указывает направление к максимуму (минимуму), и сводится к определению градиента отклика у в очередной точке факторного пространства (x1,x2,…,xk), перемещении в направлении градиента (или против градиента при минимизации) в новую точку и повторении процедуры.
Перемещения прекращаются, когда в очередной точке модуль градиента становится близким к нулю.
Fmin |
|
x2min |
X2 |
|
x1min |
|
Х1 |
а |
|
у4
x1min
у2 у3
x2min
Рис.1. Пример функции (отклика) y двух аргументов (факторов) Х1, X2 [3]:
а – поверхность в трехмерном пространстве; б – график в изолиниях; вид сверху на трехмерную поверхность (а); Х1min, X2min – искомые оптимальные значения факторов
градиента (при поиске минимума).
5.Перейти в пункт 2.
Алгоритм поиска экстремума градиентным методом (рис. 4) [7]:
1. Переместиться в произвольную начальную точку х1.
2.Найти градиент отклика у в текущей точке.
3.Если модуль градиента близок к нулю (не превышает допустимую погрешность), то текущая точка – это точка экстремума, тогда осуществляется выход из процедуры.
4.Перейти из текущей точки в новую х2, сделав шаг по направлению градиента (при поиске максимума) или против
13
Метод крутого восхождения (наискорейшего спуска) требует еще меньших затрат времени, т.к. в отличие от обычного градиентного метода, здесь в направлении выбранного градиента делается не один, а несколько шагов. Движение в выбранном направлении продолжается до тех пор, пока целевая функция (отклик)
X1 |
|
|
|
|
не |
начнет |
меняться |
в |
другую |
|
|
|
|
|
сторону, |
и |
только |
тогда |
|||
|
|
|
|
|
||||||
|
|
|
|
|
вычисляется новый |
градиент и |
||||
|
|
х2 |
|
|
направление |
поиска |
меняется |
|||
|
|
х1 |
х3 |
|
[3]. Так продолжается до выхода |
|||||
x10 |
х |
0 |
хmin |
|
отклика на экстремум, при этом |
|||||
|
|
|
в |
найденной |
зоне |
оптимума |
||||
|
|
|
|
|
||||||
0 |
x20 |
|
|
факторов |
|
коэффициенты |
||||
|
X2 |
регрессии при линейных членах |
||||||||
|
|
|
|
|
||||||
|
Рис. 2 Пример градиентного метода |
|
уравнения |
регрессии |
должны |
|||||
поиска экстремума [3] |
|
|
быть незначительными. В этой |
|||||||
|
|
|
|
|
зоне |
уже |
необходимо |
|||
использовать план эксперимента 2-го порядка для нахождения уравнения регрессии 2-го порядка, с помощью которых уточнить оптимальные факторы.
Метод многофакторного планирования (ПФП) и метод крутого восхождения по поверхности отклика дают возможность достаточно достоверно и экономично определять оптимум по двум-четырем исследуемым факторам. При большем количестве факторов более эффективным становится метод разбиения всех факторов на группы по 3-4 фактора в каждой группе. Разбивать надо таким образом, чтобы факторы, которые являются взаимодействующими, или могут быть взаимодействующими, не попали в разные группы. Если это удается сделать, то тогда оптимальная область всех факторов представляет собой совокупность их оптимальных значений в каждой группе [4].
В случае, если все факторы не возможно разбить на такие группы,
используется метод «последовательного преодоления лимитирующих факторов» с
использованием полного факторного эксперимента, представляющего собой циклическое воспроизведение метода разбиения с уточнением факторов [5].
14
Использование метода дробных реплик и отсеивающего эксперимента, основанного на ранжировании факторов, не целесообразно.
5. Построение оптимальных планов эксперимента первого и второго порядка
Эффективность и основные понятия
План эксперимента – совокупность данных, определяющих число, условия и порядок реализации опытов.
Теория планирования эксперимента позволяют существенно повысить эффективность экспериментальных исследований, повысить точность регрессионной модели или уменьшить необходимое число опытов.
Проиллюстрируем это преимущество на простом примере. С помощью пружинных весов требуется определить вес 3-х предметов: А, В, С. Обозначим через x неизвестную массу чашки весов. Обычный способ определения весов предметов состоит во взвешивании сначала чашки, а затем последовательно каждого предмета
(табл. 1) [4].
|
|
|
|
|
|
Таблица 1 |
|
|
Обычный способ взвешивания |
|
|
||||
|
|
|
|
|
|
|
|
№ взвешивания |
x |
А |
В |
С |
|
Результат взвешивания |
|
1 |
+ |
– |
– |
– |
|
y0 |
|
2 |
+ |
+ |
– |
– |
|
y1 |
|
3 |
+ |
– |
+ |
– |
|
y2 |
|
4 |
+ |
– |
– |
+ |
|
y3 |
|
Знак “+” или “–” означает наличие или отсутствие соответствующего груза на
весах.
Масса предмета, например, предмета А определится по формуле:
Ра = y1 – y0.
Погрешность измерения оценивается дисперсией D или среднеквадратическим отклонением σ.1 Так как дисперсия суммы двух случайных величин равна сумме их дисперсий , то дисперсия ошибки определения массы груза А равна: D(ра) = D(у1) + D(у0) = 2D, где D – дисперсия ошибки одного взвешивания.
Теперь рассмотрим эффективный способ определения весов грузов (табл. 10).
1 Как правило, результаты измерения подчиняются нормальному закону распределения ошибок. В этом случае разброс замеров одного и того же предмета (например, предмета А) будет находиться в интервале (М(Ра )–3σ; М(Ра )+3σ), где М(Ра ) – математическое ожидание веса предмета А. При этом вероятность того, что вес предмета А окажется вне этого интервала, будет всего лишь 0,27%.
15
Таблица 2 – Эффективный способ определения весов грузов
№ взвешивания |
х |
А |
В |
С |
Результат взвешивания |
|
|
|
|
|
|
1 |
+ |
+ |
– |
– |
y 1 |
|
|
|
|
|
|
2 |
+ |
– |
+ |
– |
y2 |
|
|
|
|
|
|
3 |
+ |
– |
– |
+ |
y3 |
|
|
|
|
|
|
4 |
+ |
+ |
+ |
+ |
y4 |
|
|
|
|
|
|
В этом случае масса предмета А определяется по формуле:
Р y1 y4 y2 y3 . |
|
а |
2 |
|
|
Дисперсия ошибки определения веса предмета А будет равна2:
y1 y4 |
y2 |
y3 |
|
|
4 D |
|
|
D(Ра ) D |
|
|
|
|
|
|
D . |
|
2 |
|
4 |
||||
|
|
|
|
|
|||
Таким образом, дисперсия ошибки определения массы груза в этом случае в 2 раза меньше, т.е. мы получаем интервал разброса (рассеяния) замеров на 12,28% меньше. Этот пример показывает, что при одних и тех же трудозатратах (проводится 4 замера) за счет использования оптимального плана эксперимента получаем более точный результат! В общем случае ТПЭ позволяет повысить точность регрессионной модели или уменьшить необходимое число опытов.
Наиболее эффективный план эксперимента позволяет находить теория планирования эксперимента (ТПЭ) [4]. В этой теории указывается, что эффективный план эксперимента должен удовлетворять определенным критериям. В соответствии с различными критериями имеется более 20 различных типов оптимальных планов. Наиболее распространенными критериями оптимальности планов, обеспечивающими удобство статистической обработки эксперимента и его эффективности (повышение точности и /или уменьшения количества опытов), являются ортогональность, ротатабельность, А, Д и G –оптимальность.
Следствием ортогональности плана эксперимента является то, что в случае незначительности каких-либо параметров модели изучаемого объекта, нет необходимости повторения процедуры построения и статистического анализа модели после исключения из нее незначительных членов. Также ортогональные планы эксперимента обеспечивают минимальную дисперсию предсказываемого по модели значения параметра оптимизации.
2Из теории вероятностей известны свойства дисперсии: D(Σyi)=ΣD(yi), D(a∙yi)=a2∙D(yi), где yi – случайная величина, а
– константа.