Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

4571

.pdf
Скачиваний:
2
Добавлен:
08.01.2021
Размер:
1.17 Mб
Скачать

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

– константа.

16

Ротатабелный план - это план, обеспечивающий равноточность исследования по всем направлениям пространства факторов.

А - оптимальный план - это план, обеспечивающий минимальное значение средней дисперсии оценок параметров модели.

Д - оптимальный план – обеспечивает минимальный объем эллипсоида рассеяния оценок параметра.

G - оптимальный план минимизирует максимальные значения дисперсии предсказываемых моделью значений параметра оптимизации. Для G -оптимальных планов в области планирования эксперимента не окажется точек, в которых точность предсказания значений параметра оптимизации будет хуже, чем в точках плана эксперимента.

Всем этим критериям удовлетворяют полные факторные планы (ПФП) и дробные факторные планы (ДФП).

ПФП – план, содержащий все возможные комбинации всех факторов на определенном числе уровней равное число раз.

ДФП – план, содержащий часть комбинаций ПФП.

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

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

Двухуровневый ПФП содержит общее количество опытов, выражаемое формулой:

N 2K ,

(8)

где к – количество факторов.

Этот план пригоден для построения линейчатых моделей вида:

k

k k

k k k

 

y b0 x0 bi xi bij xi x j ... bijk xi x j xn ...

(9)

i 1

i 1 j i 1

i 1 j i 1 n j 1

 

Количество неизвестных коэффициентов этого уравнения регрессии также определяется по формуле (8). Следовательно, двухуровневый ПФП может быть использован для нахождения всех коэффициентов уравнения (9).

17

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

Для регрессионных моделей (9), содержащих смешанные произведения факторов хi, xj, xn, …, избыточность ПФП И определяется по формуле:

И 2k (k 1) m ,

(10)

где m – количество в регрессионной модели членов со смешанными произведениями факторов.

В случае линейной модели неизвестно к+1 коэффициентов. Следовательно, для линейной модели ПФП является избыточным на величину 2к – (к+1).

Избыточные опыты плана могут быть использованы либо для уменьшения случайных погрешностей линейной модели, либо для уменьшения погрешности адекватности модели путем включения в нее произведений факторов: (aij xi x j , aijn xi x j xn , ...), либо и для того и другого вместе.

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

Методика построения ПФП и ДФП

Для унификации планов факторы в них записываются в безразмерном

нормализованном виде. Для этого вводятся величины:

 

1)

Центр эксперимента для данного фактора:

 

 

 

X

 

 

xmax

xmin

.

 

(11)

 

0

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2)

Шаг варьирования данного фактора:

 

 

 

 

 

 

X

xmax

xmin

,

(12)

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

Тогда безразмерная нормализированная величина фактора х определяется по формуле:

X

x x0

.

(13)

 

 

X

 

18

Нетрудно увидеть, что при x=xmin, Xmin= –1; при x=x0, X0=0; при x=xmax, – Xmax=1.Следовательно, в таблицах активного эксперимента используются уровни –1 ,0 , +1, в которых цифру "1" опускают.

Рассмотрим пример, когда количество факторов к=3. В этом случае Non 23 8 , а план эксперимента имеет вид, представленный в табл. 3 [4].

Строки плана представляют собой условия проведения опытов при соответствующих фиксированных значениях уровней факторов Х1, Х2, Х3.

В таблицу вводится столбец Х 0, который всегда равен +1.Это так называемый постоянный фактор, который является сомножителем свободного члена уравнения регрессии (y=b0∙X0+…), где Х0=1.

Правило заполнения столбцов ПФП

Заполнение начинают с левого столбца (X1) сверху вниз чередованием знаков "+" и "–". Следующий столбец заполняют чередованием по 2 знака "+" и "–". Третий столбец – по 4 знака "–" и 4 знака “+” Затем (если бы у нас была n-факторная модель) по 8 знаков "–" и 8 знаков "+", затем по 16 и т.д. Столбец yi – отклик.

Данный ПФП имеет избыточность, равную 23–(3+1)=4. Поэтому для нахождения параметров достаточно использовать половину из опытов ПФП.

Правило заполнения столбцов ДФП

Таблица 4 – Полный 3-факторый план эксперимента

Non

Х 0

Х 1

Х 2

Х 3

yi

1

+

y1

2

+

+

y2

 

 

 

 

 

 

3

+

+

y3

4

+

+

+

y4

 

 

 

 

 

 

5

+

+

y5

6

+

+

+

y6

 

 

 

 

 

 

7

+

+

+

y7

8

+

+

+

+

y8

 

 

 

 

 

 

Заголовки столбцов здесь те же, а количество опытов – 4 (табл. 4).

Сначала процедура заполнения такая же, как и для случая ПФП. Она проводится до того столбца, где по правилу должны быть поставлены только знаки "–" и чередование знаков отсутствовало бы (в данном случае столбец X3 должен был бы состоять из знаков "–"). В этом и последующих столбцах знаки находятся по правилу перемножения уже имеющихся знаков

данной строки (в данном случае X3= X2∙∙X1).

ДФП называют также дробной репликой ПФП.

Если дробление ПФП производится многократным последовательным делением числа его опытов на 2, то реплика называется регулярной. Число раз "р"

19

такого последовательного деления называется дробностью реплики. При Р =1 ДФП называют полурепликой, при Р =2 –четвертьрепликой, при Р =3 – 1/8 реплики, при Р =4 – 1/16 реплики и т.д.

Таблица 5 Пример дробного факторного плана

Non

X0

X1

X2

X3

yi

 

 

 

 

31∙Х2)

 

1

+

+

y1

 

 

 

 

 

 

2

+

+

y2

 

 

 

 

 

 

3

+

+

y3

4

+

+

+

+

y4

В этом случае расчет количества опытов ДФП осуществляется по формуле

N 2k p ,

(14)

а избыточность плана определяется по формуле

И 2k p (k 1) m . (15)

Методика построения плана эксперимента второго порядка

Как указывалось ранее, уравнение регрессии 1-го порядка зачастую не является адекватным и для обеспечения адекватности модели к исходному уравнению 1-го порядка добавляют члены со сменными произведениями факторов, типа bi j xi x j , bij xi x j xn и т.д., где bi j , bijn – коэффициенты регрессии смешанных произведений факторов xi , x j , xk и т.д.

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

Однако практика показывает, что во многих случаях такое добавление членов

k

смешенных произведений факторов к правой части линейной модели y b0 bi xi

i 1

и не приводит к желаемому результату обеспечения адекватности модели. В этом случае либо уменьшают интервалы варьирования одного из факторов (или, может быть, всех факторов), либо переходят к моделям 2-го, затем 3-го и т.д. порядков, в которые вводят члены, содержащие факторы во 2-й, 3-й и т.д. степенях.

Рассмотрим построение оптимального плана эксперимента 2-го порядка на примере двухфакторной модели:

y f x , x

2

b

b

x

b

x

2

b

x2

b

x

x

2

b

x2 .

(16)

1

0

1

1

2

 

11

1

12

1

 

22

2

 

20

Функция отклика у этой модели представляет собой поверхность 2-го порядка: цилиндр, конус, эллипсоид, параболоид, гиперболоид, параболический цилиндр, параболический гиперболоид или сферу.

Сечения этой поверхности, параллельные каждой из осей X1 и X2, в общем случае представляют собой кривые также второго порядка (окружности, эллипсы, параболы, гиперболы). Известно, что для однозначного определения такой кривой необходимо располагать координатами не менее 3-х ее точек, т.е. факторы X1 и X2 должны варьироваться не менее чем на 3-х уровнях.

Для удовлетворения этого требования план эксперимента 2-го порядка помимо опытов плана эксперимента 1-го порядка, называемого ядром плана 2-го порядка, должен включать опыты в звездных точках и в центре плана [4]

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

Звездной точкой плана эксперимента 2-го порядка называется точка, лежащая на координатной оси факторного пространства.

Звездная точка удалена от центральной точки на расстояние звездного плеча. Для центральных композиционных ротатабельных планов звездное плечо

рассчитывается по формуле:

 

 

 

 

 

 

 

 

 

Таблица 13

 

 

 

 

 

 

 

 

 

 

n p

 

План эксперимента 2-го порядка

 

 

 

 

 

2 2 .

 

 

 

 

 

 

 

 

 

 

 

(17)

Non

X0

X1

X2

 

yi

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

+

+

+

 

y1

 

 

 

 

 

 

 

 

2

+

+

 

y2

Опыты в

Общее количество звездных

3

+

+

 

y3

ядре

точек,

расположенных

 

на

4

+

 

y4

плана

 

 

 

 

от центра плана

5

+

0

 

 

y5

 

расстояниях

6

+

0

-

 

y6

Опыты в

равно

2k.

Все координаты

7

+

 

0

 

y7

звездны

 

звездной

точки

за исключением

8

+

-

0

 

y8

х

 

 

 

 

 

 

точках

одной равны

нулю. Например,

7

+

0

0

 

y9

 

 

 

координата

звездной

точки,

8

+

0

0

 

y10

Опыты в

 

 

 

 

 

 

 

 

9

+

0

0

 

y11

центре

лежащей

 

на

положительной

10

+

0

0

 

y12

плана

полуоси Х1, записываются так: (α,

11

+

0

0

 

y13

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0, 0, 0, … 0).

 

 

 

 

 

 

Для 2-х факторного плана эксперимента с ядром ПФП

 

 

 

 

 

 

 

 

 

 

 

2 0

 

 

 

 

 

 

 

2 4 2 .

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]