Материал: 4571

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

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].

Метод градиентного поиска подразумевает движение в направлении наиболее быстрого изменения целевой функции у (отклика, рассчитанного по регрессионной

F(x1,x2)

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 – случайная величина, а

– константа.

Источник: https://studfile.net/preview/16414606/