CC BY
109
электронное научно-техническое издание
НАУКА и ОБРАЗОВАНИЕ
_Эл № ФС 77 - 30569. Государственная регистрация №0420900025. ISSN 1994-0408_
aüi Приближенное построение множества парето в задаче многокритериальной оптимизации методом роя частиц
Архив # 04, апрель 2010
авторы: Антух А. Э., Семенихин А. С., Хасанова Р. В.
ФОТОРЕПОРТАЖИ
Первые публикации Программы и программные системы Учебные программы Студенческая весна
Приложение к журналу
Общие проблемы инженерного образования Инженер в современной России Экобионика Зарубежное образование История технического прогресса Будущий инженер Вне рубрик
Расширеный поиск Подписаться на новости
МГТУ им. Н.Э. Баумана
Введение
В настоящее время при решении задач оптимизации все более широкое распространение получают стохастические поведенческие методы [1]. Одним из таких методов является метод роя частиц (Particle Swarm Optimization, PSO), основанный на закономерностях социального поведения [2-3]. Первоначальной целью исследований в области роя частиц было обнаружение базовых принципов, благодаря которым стая рыб или птиц синхронно меняет направление движения. К настоящему времени концепция роя частиц развилась в простой и перспективный оптимизационный многоагентный метод.
В методе PSO агентами являются частицы в пространстве параметров задачи оптимизации. В каждый момент времени (на каждой итерации) частицы имеют в этом пространстве некоторое положение и вектор скорости. Для каждого положения частицы вычисляется соответствующее значение целевой функции, и на этой основе по определенным правилам частица меняет свое положение и скорость в пространстве поиска [3].
Достаточно новым является применение метода PSO в задаче многокритериальной оптимизации (Multi-Objective Swarm Optimization, MOPSO) [1]. В работе рассматривается применение этого метода для приближенного построения множества Парето в указанной задаче.
1. Постановка задачи
Совокупность частных критериев оптимальности kt [1: s] образует векторный критерием оптимальности ^РЧ, где
JC - вектор варьируемых параметров. Положим, что ставится задача минимизации каждого из частных критериев в одной и той
же области допустимых значений ^ . Запишем задачу многокритериальной оптимизации в виде
Ф(Х)=Ф(Х")
Положим, что частные критерии оптимальности нормализованы с использованием, например, относительных отклонений частных критериев от их минимальных значений:
ПОИСК
Ред.совет
Специальности ф _ ^
Рецензентам = ЯП фк :
Авторам где Архив
Ф *
Сохраним за нормализованными критериями прежние обозначения.
Введем понятие пространства критериев . Пространство критериев имеет размерность 8 (по числу частных критериев) и образуется Б ортогональными осями координат, вдоль которых откладываются значения частных критериев
(1)
Логин
оптимальности
СОБЫТИЯ IV Международная конференция по вопросам обучения с применением технологий e-learning «MOSCOW Education Online 2010»
Вторая Международная научно-практическая конференция "Инфо-Стратегия 2010: Общество. Государство. Образование"
НОВОСТНАЯ ЛЕНТА
8.06.2010 XIII Всероссийская объединенная конференция "Интернет и современное общество " (IMS-2010)
2.06.2010
РБК > А.Фурсенко: В России в ближайшие три года количество студентов сократится на 2 млн человек »
1.06.2010
В список 500 самых мощных суперкомпьютеров мира вошли 11 российских систе
31.05.2010
Международная научно-техническая конференция "Актуальные проблемы анализа и построения информационных систем и процессов
23.04.2010 Журнал "Наука и образование: электронное научно-техническое издание" МГТУ им.Баумана был включен в Перечень ведущих периодических издани
Пароль Векторный критерий оптимальности выполняет отображение множества допустимых значений х С {X} в некоторую
. область -О® С {Ф}, где {X} — пространство варьируемых параметров.
ВХОД
Введем на множестве регистрация 2
1 А с Ll-ш
отношение предпочтения >- . Будем говорить, что вектор X С предпочтительнее
забыли пароль? вектора А и писать X.1 >-Х2/если среди равенств и неравенств имеется хотя бы
одно строгое неравенство [4].
Аналогично, на множестве введем отношение доминирования: будем говорить, что векторный критерий оптимальности
Ф^Х ) Е 1)Л доминирует векторный критерий оптимальности ) Е
если X1 >-ХЭ [4].
Ф(Х1)1><КХ2)
Пресс-релизы
Библиотека
Конференции
Выставки
Доска объявлений
Архив
Ассоциация технических Университетов Информация о проекте Авторы
Координационный совет
Выделим из множества подмножество точек, для которых нет точек, их доминирующих. Это множество
называется фронтом Парето (рис. 1). Множество ^ ^ ^ж соответствующее множеству А®, называется множеством Парето (переговорным множеством, областью компромисса). Поскольку множество на рисунке 1 является выпуклым, то множество
в данном случае представляет собой дугу АВ, в которой точка А соответствует , а точка В — (^)м . Среди
точек ) € ^ С Вф нет доминируемых, поскольку ФжО^ ) ФжС^ ), но ) >ф2(Х )
Рис. 1. К определению фронта Парето (Я = 2)
Ставится задача приближенного построения множества Парето (а, тем самым, и фронта Парето) в задаче многокритериальной оптимизации (1).
2. Канонический метод роя частиц (PSO)
Множество частиц обозначим ^ — i^i»' ^Р - ^П}, где Н - количество частиц в рое (размер популяции). Б дискретный момент времени f —[0. координаты частицы -^г определяются вектором ~ »—, а ее
скорость - вектором — СЧил^Уш»—»^¿tO. Множество частиц в дискретный момент времени обозначим ~ ;
m р У — У^ V —Vе
± - количество итерации. Начальные координаты и скорости частицы -"I равны ¿0 ri соответственно, где
X* - V*
ля - матрица случайных чисел, ri - нулевая матрица.
Итерации в каноническом методе PSO выполняются по следующей схеме:
" + КJti (3)
Здесь — XI представляет собой я -мерный вектор псевдослучайных чисел, равномерно распределенных в
интервале ® ■ символ покомпонентного умножения векторов; - вектор координат частицы ^й с наилучшим
значением целевой функции ФРЦ за все время поиска fl ; ^cf - вектор координат соседней с данной частицы с
наилучшим за время поиска [Q- f] значением целевой функции - свободные параметры алгоритма. Важнейшее
в методе PSO понятие соседства частиц зависит от используемой топологии соседства и определено, например, в работе [2].
Пересчет координат частиц по формулам (2), (3) может происходить по синхронной схеме (обновление координат частиц выполняется только после определения текущих скоростей всех Н частиц) или по асинхронной схеме (расчет новых координат части производится до завершения указанных вычислений).
у* р jг
В процессе итераций вектор & образует так называемый собственный путь (private guide) частицы а вектор fit
локальный путь (local guide) этой частицы.
Свободный параметр а определяет «инерционные» свойства частиц (при 4Z <1 движение частиц замедляется). Рекомендуемое значение параметра (L равно 0,7298 [2]. В процессе оптимизации может быть эффективным постепенное уменьшение коэффициента Q. от 0.9 до 0.4. При этом большие значения параметра обеспечивают широкий обзор пространства
поиска, а малые - точную локализацию минимума целевой функции. Рекомендуемые значения свободных параметров Д? равны 1,49618 [2].
Второй компонент в формуле (2) называется «когнитивным» компонентом (по социальной аналогии) и формализует тенденцию частицы вернуться в положение с минимальным значением целевой функции. Третий компонент в формуле (2) называется «социальным» компонентом. Компонент отражает влияние на данную частицу ее соседей.
3. Метод многокритериальной оптимизации роем частиц (MOPSO)
Важной частью метода MOPSO является определение глобально лучшей частицы (в смысле формулы (1)) для каждой частицы в
популяции. В многокритериальной задаче глобально лучшую частицу следует искать на множестве Парето. С этой целью в методе
у-
МОРБО используется архив частиц Д В нем хранятся координаты недоминируемых частиц . Схема метода МОРБО представлена на рисунке 2.
Шаг 1: t=0
Шаг 2 (инициализация):
Инициализация популяции Q: For i=l to Л?
^iO — ^t; По =
End;
Инициализация архива частиц:
Шаг 3: Обновление архива частиц: -= ;
Шаг 4:
Рог ¡=1А? Вычисление вектора глобально лучших частиц: Х*и ;= РЪЖЫкйВея^Л^
Foг to п
у** = аги+и№яф»
End; End;
Шаг 6 (проверка критерия останова итераций): t=t+1; IF t<T go to шаг 3; End
Рис. 2. Схема алгоритма MOPSO
В результате работы алгоритма MOPSO на итерации £ происходит обновление архива Функция Update выполняет сравнение
координат частиц из архива с координатами частиц, полученными на текущей итерации t. Если некоторая частица текущего
р р р Р
поколения доминирует частицу j из архива, то координаты частицы Ж1 заменяя.т в архиве координаты частицы j. Если р р р
частица j доминирует частицу то частицу не добавляется в архив.
Заметим, что из приведенной схемы функции Update следует, что на первой итерации, когда архив пуст, функция Update добавляет в архив все частицы текущего поколения, которые не доминируют друг друга.
Выбор глобально лучшей частицы осуществляет функция FindGlobalBest. Существует несколько способов реализации этой функции. В данной работе используется метод «меняющихся» соседей» Хью и Эберхарта [5]. Рассмотрим суть этого метода на примере задачи двухкритериальной оптимизации. Поиск глобально лучшей частицы для каждой частицы популяции
осуществляется в этом случае следующим образом: сначала вычисляем расстояние от частицы -^i до других частиц архива -4», используя значения первого («фиксированного») критерия оптимальности ФжР^- ). Таким образом, для частицы -^f находим к ее ближайших локальных соседей. Затем, используя второй критерий ФзР^- находим наиболее оптимальную частицу из этих числа этих соседей частицы -Ч. Это частица и принимается в качестве глобально лучшей частицы для частицы -Ч.
Итерации могут продолжаться до тех пор, пока множество недоминируемых решений не перестанет меняться, либо до достижения заданного числа итераций.
5. Исследование эффективности метода MOPSO
Исследование выполнено для следующей тестовой задачи:
множество допустимых значении
В1 = {Х\о<Ъ<\ие[1Л]}
(4)
частные критерии оптимальности
(5)
Известно, что множество -
> для задачи (4), (5) имеет вид, представлены на рис. 36, а множество Парето -
г представляет
собоИ отрезок прямоИ, соединяющий точки (0, 0), (1, 1) рис. 3а.
Рис. 3. Точные множество Парето (а) и фронт Парето (б) для тестовой задачи
На основе большого количества экспериментов было установлено, что оптимальными для задачи (4), (5) являются следующие
параметры метода МОРЭО: размер популяции К =100; количество итераций Т = 500.
На рисунке 4 изображена аппроксимация множества Парето и соответствующая аппроксимация фронта Парето, полученные для задачи (4), (5) с помощью алгоритма MOPSO. Получение указанных результатов потребовало примерно 10 минут процессорного времени (расчеты производились на персональном компьютере с процессором 2,16 Гц и 2 ГБ опертивной памяти).
Рис. 4. Аппроксимация множества Парето (а) и фронта Парето (б) для тестовой задачи
Необходимо отметить, что для вычислительно простых критериев оптимальности, решение задачи (1) полным перебором на достаточно точной сетке, покрывающей множество ^ж, требует затрат процессорного времени, сравнимых с указанным выше временем. Однако, при увеличении сложности критериев, метод MOPSO позволяет находить решение приемлемой точности за существенно меньшее время.
Поскольку известно точное решение тестовой задачи (4), (5), имеется возможность количественно оценить качество полученной аппроксимации множества Парето. Для этого было вычислено среднее отклонение частицы от идеального множества (отрезка прямой с концами в точках (0,0), (1,1), см. Рис.За). На рисунке 5 приведены полученные результаты.
Сплошной линией на рисунке 5 показано среднее отклонение от точного решения (М), как функция числа итераций (0. Рисунок показывает, что метод сходится уже на первых 100-200 итерациях, после чего имеет место стагнация итерационного процесса.
Пунктиром на рисунке 5 показана величина среднего отклонения (ст) частиц от точного решения. Величина (ст) в некоторой степени
демонстрирует то факт, что в архивном множестве 41 нет частиц, которые расположены суперблизко и супердалеко от точного решения.
Рис. 5. Оценка качества аппроксимации
Дополнительно было выполнено исследование производительности алгоритма. Для этого было выявлено, каким образом с течением итераций растет размер архива и как этот размер влияет на общее время решения задачи. Некоторые результаты экспериментов приведены на рисунке 6. На рисунке сплошная линия показывает размер архива | в функции количества итераций (0, а пунктирная линия - время вычислений Г. Заметим, что в данном варианте алгоритма был использован архив А, не имеющий ограничения на его размер. Это позволило накапливать в архиве ^ все Парето оптимальные точки. Для реальных задач многокритериальной оптимизации необходимо ограничивать размер архива, чтобы его неограниченный рост не привел к нехватке памяти компьютера.
Заключение
Работа демонстрирует подход к приближенному построению множества Парето в задаче многокритериальной оптимизации с помощью метода MOPSO. Проведенное исследование показало, что данный метод, будучи относительно простым (как математически, так и в реализации), обеспечивает решение задачи с приемлемой точностью.
1С. б. Время вычислений г и размер архива | в функции количества итераций £
К недостаткам метода можно отнести то, что выбор глобально лучшей частицы сильно зависит от выбора «фиксированного» критерия [1]. Для преодоления этого недостатка планируется использовать модифицированные методы, приведенные, например, в статье [1]. В развитие работы планируется также реализация критерия останова метода, основанного на отсутствии роста размера архива в течение заданного количества итераций.
Авторы выражают благодарность Карпенко А.П. и Селиверстову Е.Ю. за плодотворные обсуждения постановки задачи и методов ее решения.
Литература
1. Mostaghim S., Teich, J. Strategies for Finding Good Local Guides in Multi-Objective Particle Swarm Optimization (MOPSO) // Swarm Intelligence Symposium: Proceeding, 2003. - pp. 26-33.
2. Карпенко А.П., Селиверстов Е.Ю. Глобальная оптимизация методом роя частиц. Обзор // Информационные технологии, 2010, № 2, c. 25-34.
3. Субботин С.А., Олейник Ан.А., Олейник Ал.А. PSO-метод, «Интеллектуальные мультиагентные методы (Swarm Intelligence)», №3, 2006, с. 55-70.
4. Карпенко А.П. Методы оптимизации [Электронный ресурс].- (http://bigor.bmstu.ru).
5. Hu X., Eberhart R. Multiobjective optimization using dynamic neighborhood particle swarm optimization // World Congress on Computational Inelligence: Proceeding, 2002.- pp. 1677-1681.
Публикации с ключевыми словами: многокритериальная оптимизация, метод PSO, метод роя частиц, множество Парето Публикации со словами: многокритериальная оптимизация, метод PSO, метод роя частиц, множество Парето Смотри так же:
• Исследование канонического метода роя частиц (PSO) для топологий типа «кольцо» и «двумерный тор»
• Обзор методов роя частиц для задачи глобальной оптимизации (Particle Swarm Optimization)
• Исследование канонического метода роя частиц (PSO) для топологий типа «клика» и «кластер»
Тематические рубрики:
• Наука в образовании: Электронное научное издание
Ассоциация технических Университетов Вузы
Информационное агентство
Координационный совет Новости УМО Вузов
телефон (8499) 263-68-67 0 RSS
Q STACK GROUP
© 2003-2010 «Наука и образование: электронное научно-техническое издание»
