CC BY
69
ИНФОРМАЦИОННЫЕ ТЕХНОЛОГИИ И СИСТЕМЫ
УДК 519.816
И. Б. Бондаренко, Е. А. Каляева, Д. Н. Кокшаров
АДАПТАЦИЯ ПАРАМЕТРОВ ГЕНЕТИЧЕСКОГО АЛГОРИТМА ДЛЯ ОПТИМИЗАЦИИ СЛОЖНЫХ ФУНКЦИЙ
Исследовано влияние параметров разработанного генетического алгоритма на процесс оптимизации сложных функций. Показана взаимосвязь параметров при поиске глобального экстремума сложной функции. Даны рекомендации по использованию генетического алгоритма в подобных задачах.
Ключевые слова: генетический алгоритм, вещественное кодирование, пропорциональный отбор, глобальная оптимизация, тестовая функция, функция Растригина, функция Де Ионга, сложность алгоритма.
Генетические алгоритмы (ГА) относятся к методам эволюционного моделирования [1] и используются, наряду с методами нелинейного программирования и перебора, как для решения комбинаторных задач, так и для поиска глобального экстремума функций. В настоящей статье рассматривается концепция выбора основных параметров генетического алгоритма для оптимизации многомерных, мультимодальных функций и придания ему адаптирующих свойств.
Классический ГА базируется на следующих процедурах: отбор, скрещивание и мутация. Задачу оптимизации сложных функций произвольного вида выразим как
бопт = ех1г(б(Х)), б(Х) = 0(хьх2,...,хп), (1)
где п — количество показателей функции качества б(Х) готовой системы.
При вещественном кодировании хромосом и их пропорциональном отборе, т.е. когда вероятность отбора хромосомы вычисляется с использованием функции [2]
* = ^, (2)
Е / (Ч-)
г=1
где,/(Ьгг) — приспособленность г-й хромосомы; Лг11г — количество хромосом или размер популяции, возникает следующая дилемма: как подобрать параметры ГА таким образом, чтобы, во-первых, гарантированно найти решение задачи (1), а во-вторых, затратить на это меньше времени.
Параметрами ГА являются:
— вероятность мутации (*м);
— точность получения результата;
— количество итераций алгоритма или количество поколений (К);
— размер популяции (ЭДц-).
В работах, посвященных исследованиям генетического алгоритма [2, 3], авторы сходятся во мнении, что вероятность мутации должна выбираться из диапазона 0,5—1 %. Так как два последних параметра в известных публикациях „оставлены без внимания", то в настоящей статье именно их оценка и будет произведена. Необходимо отметить, что влияние обоих параметров на время вычислений прямо пропорционально.
В ходе исследований был разработан алгоритм, структурная схема которого приведена на рис. 1.
Рис. 1
Для исследования влияния параметров ГА на процесс поиска была выбрана тестовая функция Растригина [2]:
п
0(X) = 10п + £{*? -10 С08(2пху)}; (3)
1=1
X е (а, Ъ), а = - 5,12, Ъ = 5,12.
Наличие тригонометрической функции в выражении (3) приводит к образованию множества локальных экстремумов: см. рис. 2, где показана функция Растригина для п=2. Глобальный экстремум будет находиться в точке х{ = 0; г = 1, п .
Рис. 2
На основе соотношения (2) и генетического алгоритма (см. рис. 1) была создана программа и исследована тестовая функция вида (3) для ряда значений параметров п и Мл- при точности е = 0,1 и рм = 0,1. Результаты тестирования представлены на рис. 3. (При значениях п, равных 3 и 5, графики сливаются с осью абсцисс.)
1000'
-»-»--■- -е--»- ф-Ш-Ш-'Щ
800
600
400
200
I I
\ • я' '
1*1 ■ ■¡■¡»1И|И|И|И| М i"iHifriHiH|fr|ii»
5 30 60 90 120 150 180 250 400 550 700 850 1000 1300 1600 1900
Рис. 3
Значение K>1200 тыс. было признано неудовлетворительным ввиду значительных временных затрат.
Для сравнения с предыдущим результатом была исследована сферическая функция или тестовая функция Де Ионга 1 [2]
Q(X) = ^ х2, x £ (-5,12; 5,12): i=1
имеющая глобальный экстремум в той же точке, что и функция Растригина. Для тех же значений параметров, что и для функции (3), исследование сферической функции привело к результатам, представленным на рис. 4.
К-103
1000
800
600
400
200
0
5
30
60 90
(5)
(6)
120 150 180 250 400 550 700 Л^ Рис. 4
Исследования тестовых функций, схожих с функцией (4), например гиперэллипсоидной:
б(Х) = Е-х2, х- е (-5,12; 5,12),
г=1
или повернутой гиперэллипсоидной:
б(X) = ЕЕ х2, х- е (-65,536; 65,536),
г=1 ] =1
приводят к результатам, сходным с представленными на рис. 4.
При анализе рис. 3 и 4 возникает вопрос: что выгоднее — увеличивать число хромосом в популяции или количество итераций ГА, т.е. число поколений? Для ответа на этот вопрос была рассчитана сложность разработанного алгоритма:
2 (7)
(8)
при больших величинах Л^.
Следовательно, при увеличении числа хромосом, например, на величину н сложность алгоритма согласно выражению (8) определяется как
0{(ЛЫ + н)2}, (9)
а при увеличении числа поколений — как
0{н<}. (10)
Таким образом, так как сложность алгоритма, рассчитанная по выражению (10), с ростом н возрастает быстрее, чем при расчете по формуле (9), то выгоднее увеличивать число хромосом. Окончательно по результатам тестирования можно сделать следующие выводы.
Ясно, что
0{4 ЛЬг п + 8 ЛЬг + Л^}.
Иш(0{4 ЛЬг п + 8 ЛЬг + Л^}) ^ 0{Л2 }
1. Если исследуемая функция многомерная и унимодальная (вида (4)—(6)), то для гарантированного нахождения ее экстремума необходимо, например, 100 тыс. поколений, состоящих из 5—400 хромосом для 3—11 переменных.
2. Если функция многомерная и мультимодальная, например вида (3), то необходимо уже 5—1800 хромосом для 3—11 переменных при том же количестве поколений.
3. Уменьшение числа поколений приводит к необходимости увеличения числа хромосом. При поиске выгоднее увеличивать число хромосом, чем число поколений, несмотря на увеличение размерности массивов и необходимость дополнительной оперативной памяти.
В заключение следует отметить, что результатом данного исследования являются рекомендации по использованию классического генетического алгоритма при поиске экстремума сложных функций. С учетом вышеизложенных аспектов алгоритм можно адаптировать к степени сложности оптимизируемой модели, что позволит повысить скорость решения подобных задач.
список литературы
1. Бондаренко И. Б., Гатчин Ю. А., Соловьев Д. В. Интеллектуальная поддержка выработки оптимальных решений в САПР оптического производства // Тр. Междунар. науч.-техн. конф. „Интеллектуальные системы" (ЛШ'08) и „Интеллектуальные САПР" (СЛБ-2008). М.: Физматлит, 2008. Т. 1. С. 115—119.
2. Панченко Т. В. Генетические алгоритмы: Учеб. пособие / Под ред. Ю. Ю. Тарасевича. Астрахань: Изд. дом „Астраханский университет", 2007. 87 с.
3. Гладков Л. А., Курейчик В. В., Курейчик В. М. Генетические алгоритмы / Под ред. В. М. Курейчика. М.: ФИЗМАТЛИТ, 2006.
Игорь Борисович Бондаренко
Екатерина Анатольевна Каляева
Дмитрий Николаевич Кокшаров —
Сведения об авторах
канд. техн. наук, доцент; Санкт-Петербургский государственный университет информационных технологий, механики и оптики, кафедра проектирования компьютерных систем; E-mail: [email protected]
студентка; Санкт-Петербургский государственный университет информационных технологий, механики и оптики, кафедра проектирования компьютерных систем; E-mail: [email protected] канд. техн. наук; ЗАО „НАВИС", Санкт-Петербург; руководитель проекта; E-mail: [email protected]
Рекомендована кафедрой проектирования компьютерных систем СПбГУ ИТМО
Поступила в редакцию 08.02.10 г.
