Модификация метода последовательного симплексного планирования и ее применение к решению задач оптимизации Текст научной статьи по специальности «Математика»

39. Новый метод утилизации изношенных шин. - Экспресс - информация. Шинная промышленность. - М., ЦНИИТЭНефтехим, 1989, № 8, с. 22 - 24.

40. Gummibereifung, 1985, Bd. 61, № 5, s. 28.

41. Довженко И.В., Регишевская И.Д. Термическое обезвреживание кислородсодержащих отходов, пути совершенствования, экологическая оценка метода. - Тезисы докладов совещания «Состояние и перспективы использования изношенных шин и полученных из них продуктов» (Чехов 1989).- М., ЦНИИТЭНефтехим, 1989, с. 28 - 30.

42. Макаров В.М., Дроздовский В.Ф. Использование амортизованных шин и отходов производства резиновых изделий. - Л., «Химия», Ленинградское отделение, 1986 249 с.

Модификация метода последовательного симплексного планирования и ее применение к решению задач оптимизации

к.т.н. доц. Овсепян B.C., Дерцян A.C.

Ванадзорский государственный педагогический институт им. Ованеса Туманяна

(Армения) vardges193 7@mail. ru

Аннотация. Предложена модификация алгоритма Нелдера-Мида, при которой точка, определяющая направление отражения «наихудшей» вершины, выбирается с учетом значений минимизируемой функции в остальных вершинах симплекса. Исследована эффективность предложенной модификации на тестовых функциях.

Ключевые слова: симплекс, модификация, распределение, по времени, оптимизация.

Предлагаемая модификация

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

нелинейной оптимизации для поиска минимума целевой функции f (х), х е Rn, f : Rn ^ R.

Симплексом в n-мерном евклидовом пространстве называется многогранник с количеством вершин n+1, не лежащих в одной гиперплоскости. Для случая n=2 - это треугольник, для случая n=3 - тетраэдр. Суть метода заключается в последовательном перемещении и деформировании симплекса вокруг точки экстремума с помощью процедур отражения, растяжения, сжатия и редукции симплекса [2-4]. Каждая итерация начинается с

упорядочения вершин xi , х 2, ..., х n+i таким образом, чтобы было выполнено условие:

f (Xi) <f(X2) <■■.< f (X«+i). (1)

С точки зрения минимизации Xi считается наилучшей, a xn+i - наихудшей вершиной

симплекса. Затем вычисляется центр тяжести хс первых n точек симплекса, где

— 1 n — — _V

Xc — / Xi .

n i=i

На линии Xc +a (Xc - xn+i), б e R, проходящий через точки x„+i и Xc ,

nn T^n

рассматриваются точка отражения хг e R , точки внутренного и внешнего сжатия xc е R и Xoc е Rn, точка растижения хе е Rn. Исходя из соответствующих проверок, xn+i вершина симплекса заменяется одной из точек хг, хе, Xc, Xoc , для которой целевая функция уменьшает свое значение (рисунок 1).

Рисунок 1

Если данный шаг проходит неудачно (точки с «лучшим» значением функции на линии отображения не найдено), симплекс редукцируется: XI = Х1 + у(XI - Х1) сг = 0.5, г = 2, 3,...,п +1. Цикл поиска минимума целевой функции повторяется до тех пор, пока не выполняется критерий остановки, например | / (хя+1) - / (Х1)|<£, 0. Если требуемая точность достигнута, выводится значение / (Х1) как минимальное значение целевой функции.

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

— 1 п —

В методе Нелдера-Мида центр тяжести Хс = — V XI равномерно удален от первых п

п ы

вершин симплекса. Предлагаемая модификация основана на идее весовых коэффициентов. В этом случае по формуле (2) рассчитывается средневзвешенный центр грани, противостоящей «наихудшей» вершине.

X 'с = Л^Х 1 + Л2 X 2 + Хп , (2)

п

где = 1, Л, > 0, 1 = 1,п (3)

г=1

Предлагается выбрать весовые коэффициенты Лг, 1 = 1, п так, чтобы центр х'с был бы

ближе к той вершине симплекса, по направлению к которой целевая функция быстрее убывает.

По формуле (4) сначала вычисляют весовые коэффициенты ¡л{,1 = 1,п.

/(Хп+1 ) - /(X г) , - ,

^ = ^--^-/ , где р(Хп+1, х.) =

р( Хп+1, X г) \

К^ -х^)2 (4)

J=1

Здесь р(Хп+1, хг) - расстояние между вершинами Хп+1 и XI = 1,п, а числа (хг1 ,хг2,...,хгп) - координаты XI= 1,п +1. Весовой коэффициент д. есть относительное

изменение целевой функции /(х), на вершинах симплекса Хп+1 и х.. Чем больше это относительное изменение д, тем больше вероятность того, что значения целевой функции

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

больше, г = 1, п. Отметим, что весовые коэффициенты д не удовлетворяют условию

^Мг = 1, поэтому по формуле (5) рассчитывают из коэффициентов ¡л{ новые нормированные весовые коэффициенты Л. -1, п.

п

Л1 /¡л, г = 1,п; М = • (5)

1=1

Для всяких 1, j, если д. < /и^, то Л. < и соотношение весовых коэффициентов не изменяется. Полученные коэффициенты Л. полностью удовлетворяют условию (3). После вычисления весовых коэффициентов, определяется х'с средневзвешенный центр грани симплекса:

х'с = Л1Х1 +Л2 х 2 + ... + Лп х п

Весовой центр х'с ближе к той вершине симплекса, для которой весовой коэффициент

больше. Затем вершина хп+1 заменяется точкой линии Хс +а (х'с-хя+1), б е Я, уменьшающей значение целевой функции. Данная модификация меняет направление поиска новой вершины симплекса, что повышает эффективность алгоритма. В качестве примера

рассмотрим тестовую функцию "Тпё": /(х) = /(х1 ,х2 )= (х1 -1)2 +(х2 -1)2 - х1 х2.

Данная функция принимает свое минимальное значение /(х ) = —2 на точке х =(2;2). На рисунке 2 представлен симплекс с вершинами х 1, х2, хз е Я2, где Х1 =(2.5;0.3); Х2 =(-1,1.2); Хз =(0.6;-2.3):/(Х) = 1.99; /(Х2) = 5.24; /(Хз) = 12.43.

Поскольку /(Х1) < /(х2) < /(х3), то вершина Х1 считается наилучшей вершиной симплекса, ах 3 — наихудшей вершиной. Коэффициенты ¡и1

будут рассчитываться

следующим образом:

А = /(Х3)-1(Х1) «3.24; = /(Х3)-1(Х2) «1.87

р(Х 3 , Х1 ) р( Х 3 , Х 2 )

Рисунок 2

С помощью формулы (5) определяем коэффициенты Л1,Л2, для которых Л1 +Л2 = 1: Л1 /(д) = 3.24 / 5.11« 0.634; Л2 = /л2 + /л2 ) = 1.87 / 5.11« 0.366

Тогда х'с весовой центр симплекса будет определяться следующим образом:

х'с = 0.634X1 + 0.366X2.

Из рисунка 2 видно, что полученный х'с весовой центр симплекса находится ближе к — _* —

наилучшей вершине Х1 и к точке минимума целевой функции х , чем центр тяжести Хс, определяющийся методом Нелдера-Мида. Тем самым, с большой вероятностью можно сказать, что на точках линии 12 целевая функция примет меньшие значения, чем на точках линии 11, (линия 11 проходит через точки Хс и х3, а линия 12 - через точки х'с и х3). Для достижения глобального минимума целевой функции нужно повторять процедуру поиска несколько раз.

Описание алгоритма

Шаг 1.1. В области определения целевой функции /(х), / : Яп ^ Я строят начальный регулярный симплекс £ = {х1, х 2..Хп+1}. Для этого сначала случайным образом генерируют вершину Х1 е Яп, а остальные вершины симплекса определяют следующим образом:

XI = х 1 + Не I -1, 1 = 2, п +1, И > 0,

где е е Яп, 1 = 1,п - единичные координатные векторы, И - ребро симплекса. Вычисляют значения целевой функции /(Х{ ),г = 1, п +1. Выбирают £> 0 число, £ = 108.

Шаг 1.2. Вершины симплекса сортируют так, что /(Х1) </(х2) <...< /(Хп+1).

Основной цикл

Пока /(хп+1)-/(х11 > £, повторяется последовательность 2-6 шагов, в противном случае работа алгоритма завершается и выводится вершина симплекса Х1 и значение / (х1). Шаг 2. Определяют весовые коэффициенты д - 1,п ,

^ = .1(х"+0—1(хг) , где п1, х г) =

р(хп+1, хг) \

X (*п+1,: хг ,з) ,

J=1

р{хп+1, хг) - расстояние Хп+1 и XI вершин, (хг1 ,хг2,... ,х^п) - координаты вершин

хг-, г = 1, п + 1.

Шаг 3. Нормированием коэффициентов д вычисляют новые весовые коэффициенты

А

Лг = д / ¡л, г = 1, п; ¡л = ^ д

г=1

Коэффициенты Лг удовлетворяют условию = 1, Лг > 0, г = 1, п .

г=1

Шаг 4. Определяют х'с е Яп средневзвешенный центр первых п вершин симплекса.

х'с = Л1Х1 +Л2 X 2 + ... + Ап Хп

Шаг 5. На линии 12 = {х'с +а (х'с - Хп+1)}, а<Е Я , определяют точку отражения хг е Яп.

хг - х'с +а{х'с - Хп+1) при а = 1.

Находят значение целевой функции / (хг) на этой точке. Дальше, как в методе Нелде-

ра-Мида, в зависимости от значения /(хг), наихудшая вершина симплекса хп+\ шагами 5.15.4 заменяется новой точкой, находящейся на линии 12.

Шаг 5.1. Если /(хг) е[/(Х1), /(хп)), то наихудшая вершина симплекса хи+1 заменяется точкой отражения хг, делается переход на шаг 2.

Шаг 5.2. Если /(хг) < /(Х1), проводят операцию растяжения. хе - Хс +а(х'с -Хп+1), при а - 2

Если /(хе) < /(хг), то наихудшую вершину симплекса Хп+1 заменяют точкой хе и

переходят на шаг 2. В противном случае вершину Хп+1 заменяют точкой хг и переходят на шаг 2.

Шаг 5.3. Если /(хг) е [/(хп), /(Хп+1)), делают операцию внешнего сжатия

Хос - х'с +а (х'с - Хп+1), при а = 0.5 Если точка Хос е Я" лучше, чем точка отражения, то есть /(Хос) < /(хг), вершина Хп+1 заменяется точкой Хос, переходят на шаг 2. В противном случае переходят на шаг 6. Шаг 5.4. Если /(хг) > /(Хп+1), делается операция внутреннего сжатия

хс - х'с +а( х'с - Хп+1), при а = -0.5 Если точка хс лучше, чем вершина симплекса хя+1, то есть /(хс) < /(Хп+1), вершина

Хп+1 заменяется точкой хс, переходят на шаг 2. В противном случае переходят на шаг 6.

Шаг 6. Проводим операцию редукции симплекса.

Х1 = Х1 + у (Х1 - Х1), у = 0.5, [ = 2,п +1.

Тестирование предложенной модификации

Приведенные ниже тестовые функции оптимизации были тестированы методом Нелде-ра-Мида и предложенной модификацией весовых коэффициентов (ВК), каждый по 100 раз. Каждый раз первая вершина начального симплекса была выбрана из области поиска тестовой функции произвольным образом. В таблице 1 приведено сравнение методов Нелдера-Мида и предложенной модификации по среднему процессорному времени (СРи-сек.) и среднему числу расчета значений тестовой функции (вызовов функции) для одного теста. Тестирование проводилось на компьютере с 2,66 ГГц процессорной мощностью и 3 вЬ оперативной памятью.

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

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

л

средние

минимального значения тестовой функции "Тпё" /(х-1)2xixi_

V г=1 1=2 )

затраты процессорного времени для п=2 сокращаются на 2.08%-ов, для п=4 - на 8.77%-ов, а для п=6 - на 18.34%-ов.

Таблица 1

Сравнение методов Нелдера-Мида и предложенной модификации (п - число переменных)

N Название тестовой функции n Нелдер-Мид Модификация ВК

CPU-сек. Число выз. функции CPU-сек. Число выз. функции

1 Trid Function 2 0.048 101.58 0.047 102.25

2 Trid Function 4 0.171 285.69 0.156 253.26

3 Trid Function 6 0.529 567.08 0.432 471.5

4 Zakharov Function 2 0.048 109.69 0.046 108.75

5 Zakharov Function 4 0.181 291.62 0.163 263.14

6 Zakharov Function 6 0.583 619 0.479 515.71

7 Helical valley Function 3 0.131 255.49 0.141 281.52

8 Gaussian Function 3 0.123 183.54 0.113 177.2

9 Box 3 dim. Function 3 0.127 267.05 0.099 266.1

10 Colville Function 4 0.372 603.63 0.368 601.57

11 Branin Function 2 0.047 108.67 0.043 107.39

12 Sphere Function 3 0.081 165.87 0.076 159.01

13 Sphere Function 5 0.241 338.97 0.209 291.45

14 Sphere Function 10 2.347 1144.98 1.132 682.38

15 Sum Squares Function 3 0.087 177.16 0.079 168.29

16 Sum Squares Function 5 0.259 360.28 0.224 313.69

17 Sum Squares Function 10 2.939 1276.25 1.386 781.67

18 Rotated hyper-ellipsoid 3 0.106 198.44 0.105 193.21

19 Rotated hyper-ellipsoid 5 0.305 420.46 0.276 380.31

Список тестовых функций

1. Trid Function: f (х)=^(xt -1)2 -

i=l i=2 .2 ^ „ / „2

Область поиска: — n < xt < n , i = 1, 2,...,n.

Минимум функции ( x* - точки минимума приведенных тестовых функций): при n = 2, /(x*)= -2 ; при n = 4, /(х*)= -16 ; при n = 6, f (x*)= -50

n f n Л2 / n Л4

2. Zakharov Function: /(X) = ^ хг2 + ^ 0.5ixi I + ^ 0.5гхг I ,

i=i V г=1 J V г=1 J

Область поиска: -15 < хг < 15, i = 1,n. Минимум функции: /(x*) = 0

,12 , ! r, Si 2

3. Helical valley Function: f (X) = 100[x3 - 10d{x1 ,x2 )]2 +100 A/x12 + x22 -1

+ x.

^(x1 ,x2 ) =

— arctg(x2/x1 ), когда x1 > 0 2n

— arctg(x2/x1) + 0.5, когда x1 < 0 „ 2n

Область поиска: -10 < xi < 10, i=1,2,3. Минимум функции: X* =(1,0,0), f (х*)= 0

метод Нелдера-Мида ■ модифицированный метод весовых коэффициент»

W JM4

Ihdttl ijlturmOl en^ifViJl ttrr^jlvJIt-'f (Jj СипкиЦ) hxlAin 01 Www! И SunWnOf hu hfprf rliipk |J|

A MS OW

47% Itw 0U1

0U1

011»

ООП

о \n

я 00Й1

0 076

■ OCt 7 ООП

I 0 104

f> W

T(J\

Lift

1№4

ni4

11=2 3

itJfmc ilu

10*",

u=4. S ъ

Я H-IC

Itm

flTM*

n = 10

Mi»

Рисунок 3. Сравнение средних затрат процессорного времени (сек.) для поиска

минимума тестовых функций

=■

4. Gaussian Function: f (х) = ^

i=i

8 - i

x1exp

x 2 (7i x3)

2 Л

V

■yt

2

i = 1,..., m;m = 15

y=(0.0009, 0.0044, 0.0175, 0.054, 0.1295, 0.242, 0.3521, 0.3989, 0.3521, 0.242, 0.1295, 0.054, 0.0175, 0.0044, 0.0009);

2

2

Серия 4. Химическое машиностроение и инженерная экология Область поиска: -1.5 < xi < 1.5, i=1,2,3. Минимум функции: f (jé ) = 1.12793 *10 8

5 Box 3 dimensional Function: f (X)= ]Г[e^'1"1 -e^ -x3(e_a1i -e~iJ]2,

i =1

Область поиска: - 50 < xt < 50, i=1,2,3 . Минимум функции: f (x * )= 0 6. Colville Function:

f (x ) = 100(x12-X2)2+(x1-1 )2+(x3-1 )2+90(x32-X4)2+10.1 ((X2-1 )2+(X4-1 )2)+19.8(X2-1 )(X4-1) Область поиска: -10 < xt < 10, i = 1,2,3,4. Минимум функции: f (x*) = 0

( 5 5 Y ( 1 ^

7. Branin Function: f (X) = [ x2---Xj2 + — x1 -6 +10 1--Icos{x1)+10,

4p P ; ^ 8pj

Область поиска: - — < xj < 10, 0< x2 < 15. Минимум функции: f (x*)= 0.397887

x2.

8. Sphere Function: / (X) = ^;

г =1

Область поиска: -2.56 < хг < 5.12, i = 1,2,____ n. Минимум функции: /(х*)= 0

n

9. Sum Squares Function: / (X) = ^ ixt 2,

i=1

Область поиска: - 5 < xt < 10, i = 1,2,..., n. Минимум функции: /(х *) = 0

10. Rotated hyper-ellipsoid function: f(X^

/ л2 x,

J

i=1 VJ=1

Область поиска: - 65 < xi < 65, i = 1, n. Минимум функции: f (х* )= 0

Применение метода симплексного планирования в задачах оптимизации

Предложенная модификация последовательного симплексного планирования может быть применена для решения задачи равномерного распределения ресурсов по времени. Пусть организация осуществляет комплексную программу, состоящую из n проектов. Известны функции f (t), t e[0,di], i=l,...,n, характеризующие потребность ресурсов в момент времени t для реализации проекта i. Задана продолжительность каждого проекта -di, самый ранний момент начала - ci и самый поздний момент окончания - vi, i -1, n. Если реализация i-ro проекта начнется с задержкой ti е \ct, vi - di ], распределение ресурсов будет соответствовать функции fi (t - tt), где t e[ti, tt + dt ], i - 1,n . Суммарная потребность в

n

ресурсах в момент времени t равна ^ f (t - tt). Требуется выбрать такие моменты начала

i=1

выполнения проектов t1 ,t2,...,tn, чтобы распределение суммарных ресурсов по времени было

наиболее равномерным. Под равномерным распределением ресурсов понимается такое распределение, когда значения суммарной функции как можно меньше отличаются от средней потребности M, где

1 b n

M = T— iE f (t - ti )dt,

b - a J 7~r

a i-1

a = min tt = 0, b = max tt + dt, (b-a) - продолжительность комплексной программы.

i i

В [1] в качестве оценки близости потребности ресурсов к среднему значению выбран квадратичный критерий. Тогда задача равномерного распределения ресурсов сходится к задаче

min

ti Aj

b - a

IE / (t - и)-M

dt.

Литература

1. Цирлин A.M. Вариационные методы расчета химических аппаратов.-М.: «Машиностроение», 1978

2. Nelder J.A, Mead R. A simplex method for function minimization. Computer J. 7:308-313, 1965

3. Hedar A.R, Fukushima M. Simplex coding genetic algorithm for the global optimization of nonlinear functions, Graduate School o/In/ormatics, Kyoto University, Kyoto 606-8501, Japan, 2002

4. Овсепян B.C., Дерцян A.C. Об одной модификации последовательного симплексного метода. Сборник трудов международной II конференции Горисского университета, Горис, 2011

2

b

1

a

Полимеризация порошковых красок с использованием ИК-излучателей с наноструктурированным керамическим покрытием

Шкарин Н.Ю.1, д.т.н. Рахимов Р.Х.2, к.ф.-м.н. в.н.с. Казенин Д.А.1

1 Университет машиностроения 2Институт материаловедения НПО «Физика-Солнца» АНРУз

8(488)267-16-69, n_shkarin@mail. ги

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

Ключевые слова: порошковые краски, ИК-излучатели, ИК-полимеризация, энергосбережение.

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

Такие загрязнения уменьшают устойчивость защитно-декоративных покрытий, повышают скорость коррозионных процессов и в конечном итоге служат одной из причин, приводящей к снижению надежности машин и агрегатов.

По сравнению с традиционными лакокрасочными материалами порошковые краски обеспечивают практически безотходную технологию производства покрытий, также физико-механические свойства покрытий из порошковых красок по многим факторам превосходят покрытия из жидких лакокрасочных материалов [1].

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

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