Исследование методов эволюционной оптимизации в задачах организационного управления производством Текст научной статьи по специальности «Компьютерные и информационные науки»

А.Н. Применение генетических алгоритмов для решения оптимизационных задач / А.Н. Г воздинский, Е.Г. Климно // Харьков: Радиоэлектроника и информатика. 2001. №2. С.109-113. 4. Гвоздинский А.Н. Применение методов эволюционной оптимизации для решения задач производственного планирования / А.Н. Гвоздинский, В. А. Ма-лышкин, С.В. Ткачев // АСУ и приборы автоматики. 2011. №154. С. 97-102.

Поступила в редколлегию 27.11.2013 Рецензент: д-р техн. наук, проф. Куземин А.Я.

Гвоздинский Анатолий Николаевич, канд. техн. наук, профессор кафедры искусственного интеллекта ХНУРЭ. Научные интересы: оптимизация процедур принятия решений в сложных системах управления. Адрес: Украина, 61166, Харьков, ул. Академика Ляпунова, 7, кв.9, тел. 702-38-23.

Руднев Вячеслав Константинович, студент группы КН-10-5 кафедры искусственного интеллекта ХНУРЭ. Научные интересы: методы принятия решения в системах искусственного интеллекта. Адрес: Украина, 61000, Харьков, пр. Ленина, 19, кв.80, тел. (093)-764-32-24.

УДК 519.7

ИССЛЕДОВАНИЕ МЕТОДОВ ЭВОЛЮЦИОННОЙ ОПТИМИЗАЦИИ В ЗАДАЧАХ ОРГАНИЗАЦИОННОГО УПРАВЛЕНИЯ ПРОИЗВОДСТВОМ

ГВОЗДИНСКИЙ А.Н., КОЗЛОВА А.Е.__________

Предлагается подход к решению проблем управления запасами, основанный на эволюционных методах, возникающих при поиске рационального плана раскроя. Демонстрируется возможность сравнительно легкой адаптации алгоритма к различным вариантам постановки задачи. В качестве объектов исследования рассматриваются системы интеллектуального анализа данных и объекты производственной деятельности предприятий и фирм различного назначения. Рассматривается задача управления запасами в системах производственного планирования на примере задачи оптимального раскроя ресурсов.

1. Введение

В данной работе приводится формальная постановка задачи, на решение которой изначально был направлен первый генетический алгоритм (ГА), и общая схема работы самого алгоритма. Описываются составляющие ГА и их назначение. Для решения конкретной задачи разработчик, исследователь может сконструировать свою конкретную схему работы генетического алгоритма из его составляющих. При этом разработчик руководствуется особенностями предметной области, формализации задачи, структурой используемых данных, возможно, даже результатами тестирования других схем ГА. Разумеется, в строении всех таких алгоритмов есть общие элементы и последовательность действий, которые являются основой генетических алгоритмов и отличают их от многих других. Данная статья посвящена опыту конструирования конкретных вариантов ГА.

2. Анализ методов эволюционной оптимизации

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

Нейронные сети относят к классу нелинейных адаптивных систем, строением они условно напоминают нервную ткань из нейронов. Это набор связанных друг с другом узлов, получающих входные данные,

которые осуществляют их обработку и вырабатывают на выходе определенный результат. На узлы нижнего слоя подаются значения входных параметров, на их основе производятся вычисления, необходимые для принятия решений, прогнозирования развития ситуации и т. д. Эти значения рассматривают как сигналы, которые передаются в вышележащий слой, усиливаясь или ослабляясь в зависимости от числовых значений (весов), приписываемых межнейронным связям. На выходе нейрона самого верхнего слоя вырабатывается значение, которое рассматривается как ответ, реакция всей сети на введенные начальные значения. Т ак как каждый элемент нейронной сети частично изолирован от своих соседей, у таких алгоритмов имеется возможность для распараллеливания вычислений. Нейроны сети обрабатывают входные данные, для которых известны и значения входных параметров, и правильные ответы на них. Обучение состоит в подборе весов межнейронных связей, которые обеспечивают наибольшую близость ответов сети к известным правильным ответам. После обучения на имеющихся данных сеть готова к работе и может быть использована для построения прогнозов поведения объекта в будущем. Нейронные сети используются для решения задач прогнозирования, классификации или управления.

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

Недостаток - необходимость иметь большой объем обучающей выборки. Окончательное решение зависит от начальных установок сети. Данные должны быть обязательно преобразованы к числовому виду. Полученная модель не объясняет обнаруженные знания.

Деревья решений используют разбиение данных на группы на основе значений переменных. В результате получается иерархическая структура операторов «Если..., то...», которая имеет вид дерева. Для классификации объекта или ситуации нужно ответить на вопросы, стоящие в узлах этого дерева, начиная от его корня. Если ответ положительный, переходят к правому узлу следующего уровня, если отрицательный - к левому узлу и т.д. Заканчивая ответы, доходят до одного из конечных узлов, где указывается, к какому

РИ, 2013, № 4

91

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

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

Недостаток - проблема значимости. Данные могут разбиваться на множество частных случаев, возникает “ кустистость” дерева, которое не может давать статистически обоснованных ответов. Полезные результаты получают только в случае независимых признаков.

Генетические алгоритмы на данный момент являются наиболее известным представителем эволюционных методов оптимизации. Они содержат все основные генетические операции. ГА получены в процессе обобщения и имитации в искусственных системах свойств живой природы:

- приспосабливаемость к изменениям среды;

- естественный отбор;

- наследование потомками наиболее «ценных» свойств родителей.

С помощью Г А можно улучшать работу поисковых систем, которые требуют обработки больших массивов информации; оптимизировать работу нефтяных трубопроводов; распределять инструменты в металлообрабатывающих цехах; осуществлять оптимизацию профилей балок в строительстве. Решение задач комбинаторной оптимизации является одной из основных областей применения ГА. Г енетические алгоритмы имитируют процесс естественного отбора в природе. Для нахождения более оптимального критерия все решения описываются набором чисел или величин нечисловой природы. Поиск оптимального решения похож на эволюцию популяции индивидов, которые представлены наборами их хромосом. В этой эволюции действуют три механизма:

- отбор сильнейших наборов хромосом, которым соответствуют наиболее оптимальные решения;

- скрещивание - получение новых индивидов при помощи смешивания хромосомных наборов отобранных индивидов;

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

Наиболее популярное приложение генетического алгоритма - оптимизация многопараметрических функций. Реальные задачи формируются как поиск оптимального значения сложной функции, зависящей от некоторых n - входных параметров. Сила ГА - в способности манипулировать одновременно многими параметрами. В одних случаях получено точное

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

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

Существует множество модификаций генетического алгоритма, которые отличаются методами отбора, скрещивания, хромосомной кодировки. Традиционный (основной) ГА работает с двоичной строкой определенной длины, использует следующие свойства:

- на каждом поколении реализуется отбор пропорционально функции приспособленности;

- одноточечный кроссовер (скрещивание);

- мутация.

Как метод оптимизации Г А обладает внутренним параллелизмом: различные существенные комбинации генов отыскиваются параллельным образом, одновременно для всех комбинаций. Чем меньше комбинация, тем легче ее определить. В связи с тем, что алгоритм в процессе поиска использует некоторую кодировку значений, а не сами значения, он эффективно может применяться для решения задач дискретной оптимизации, которые могут быть определены как на числовых, так и на конечных множествах произвольной природы. Так как в качестве информации об оптимизируемой функции для работы алгоритма используются лишь значения в рассматриваемых точках пространства поиска, то данный алгоритм применим к широкому классу функций. Генетический алгоритм успешно справляется с поставленной задачей даже там, где не существует общеизвестных алгоритмов решения или высока степень априорной неопределенности.

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

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

92

РИ, 2013, № 4

3. Формулировка задачи и построение - 19*110 мм, необходимость в которых Пь П2, ПЗ, П4

математической модели единиц соответственно.

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

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

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

Шаг 0. Создать случайную начальную популяцию {Sk(0)}.

Шаг 1. Вычислить приспособленность f(S ) каждой

k

особи (S ) популяции {S (t)} . k k

Шаг 2. Производя отбор особей из {S^t)} в соответствии с их приспособленностями f(Sk ) и применяя генетические операторы (скрещивания и мутации) к отобранным особям для получения потомства, сформировать популяцию следующего поколения {Sk(t +1)}.

Шаг3. Повторять шаги 1,2 для t = 0,1,2,... до тех пор, пока не выполнится некоторое условие окончания эволюционного поиска (прекращается рост максимальной приспособленности в популяции, число поколений t достигает заданного предела).

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

- цепочки символов в хромосомах бинарные (символы S^ принимают значения 0 либо 1), длина цепочек

постоянна (N = const);

- метод отбора - пропорционально-вероятностный;

- скрещивания производятся по схеме одноточечного.

Рассмотрим распределение ресурсов на примере решения задачи оптимального раскроя [2].

На заготовочный участок радиоцеха поступили листы фольгированного гетинакса размером 11хІ2 мм. Его необходимо разрезать на заготовки для печатных плат: П1 размером 1зхІ4 мм, П2 - 1зх1§ мм, П3 - 1ух1§ мм, П4

Построить модель, на ее основании составить экстремальную задачу, решение которой дает минимальное количество разрезанных выходных листов. Найти решение задачи в условиях минимизации отходов раскроя. Предусмотреть возможность максимизации количества комплектов печатных плат при условии выполнения комплектности П1: П2: Пз :Щ= =k^k2:k3 :k4 (таблица).

Параметры Значение

11 136

12 158

Із 92

14 36

15 51

1б 88

І7 46

18 40

19 50

І10 50

П1 270

П2 200

Пз 310

П4 90

Компетентность 3:2:1:5

Исходные данные: 11х12, П1, П2, ПЗ, П4, k1, k2, k3, k4.

Используемые переменные: П j = a^ - необходимое

количество заготовок i-го типа; C ■ - отход при

раскрое единиц исходного материала j-м способом;

j = 1m - количество вариантов раскроя; a.. - коли’ У

чество заготовок i-го типа (i = 1,4), полученных при раскрое единицы исходного материалаj -м вариантом; i - вид заготовки; q - количество изделий (П1, П2, П3, П4); a - количество разрезанных листов; K - количество типов листов, поступающих на раскрой (K = 1); 1i - компетентность-количество заготовок i-го вида в одном комплекте; x . -количество единиц исходного материала, который 'раскраивается j-м вариантом;

4. Формализация исследуемой задачи

Сформируем ограничения для данной задачи:

m

Е x ■ = a j = 1 j

(1)

РИ, 2013, № 4

93

x . > 0, j = 1,m J

m

z

J = 1

a..x. 1J J

= П.

і •

Функция цели имеет вид:

m

F(x) = Z x . ^ min,

J = 1 J

F(x) = a ^ min . Функция цели имеет вид:

(2) при условиях существования ограничении вида

AX < B , (10)

(3) A =ll aij llmn,X =11 xi,x2,...,xn ll,B =11 Ь15Ь2,...,ЬШ ||.

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

(4)

(5)

x , ,x 2 , n +1 n + 2 приобретает вид:

’ xn + m . Вследствие этого система

n

Z C x , ^ min . (6)

j = 1 j j

Необходимо предусмотреть возможность максимизации количества комплектов печатных плат при условии выполнения комплектности П1:П2:Пз:П4=

=к1:к2:кз:к4.

Функция цели имеет вид:

F(x) = X ^ max, (7)

1 m

где X = П1Г( Z a.ix.).

Пlkl.=1J1J

А ограничения (3) примут следующий вид:

p p

Z ax. = П1k1X Z ax. = П2k2X . j1 J 1 1 , . j2 j 22,

J = 1 J = 1

P P

Z a x. = П к X Z a x. = Пk X . . J3 J 33, . j4 j 44.

J = 1 J = 1

(8)

5. Алгоритмы решения исследуемых задач

Для начала рассмотрим решение данной задачи на основе методов линейного программирования [3].

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

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

Общая постановка задачи линейного программирования - определить величины xj > 0,J = 1,n, которые максимизируют линейную форму

n

Z =

Z

J=1

C x . ^ max JJ

(9)

a x + a x +. ..a x +x = b

11 1 12 2 1n n n + 1 1

a x + a x +. ..a x +x =b

21 1 22 2 2n n n+1 2

am1x1 + am2x2 +. ..a x mn n +x n +1 = b m

(11)

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

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

во-первых, задача раскроя любого вида является ОТ-полной;

во-вторых, время получения решения для заказчиков предполагается порядка десятки секунд для размерности 100, это достаточно высокие требования ко времени;

Для реализации ГА особью предлагается считать любое размещение заданных предметов на полосе, удовлетворяющее условию задачи. Пусть все предметы некоторым образом пронумерованы. Тогда пусть хромосома — это перечисление всех номеров предметов в некотором порядке, отражающем их размещение. Физический внешний вид особи, конкретные координаты размещения предметов на полосе можно получить, выложив предметы на полосу в порядке их перечисления в хромосоме по правилам переработки — декодирования. Эти правила называются процедурой декодирования, или декодером. Различные декодеры могут работать по разной схеме, но суть одна — выложить предметы на полосу согласно некоторому правилу в порядке, указанном в хромосоме. Правила являются эвристическими и содержательно строятся так, чтобы размещение было по возможности плотным.

В качестве критерия оптимизации, критерия оценки пригодности особи в поставленных задачах рассматривается длина полосы, занятая размещением предметов.

94

РИ, 2013, № 4

Пустоты учитыв аются и нумеруются, начиная с числа, на единицу большего количества предметов, но со знаком минус. Декодер, работающий с такой структурой данных, называется блочным. Размещение представлено в виде совокупности кортежей (блоков): (1., 2., ...)х, где j = 1,2,...

А именно: (1; 2) хь (3; 4; -9; 2) Х2, (3; 5; 6) хз, (7; -10; 5; 6) х4, (7; 8; 6) х5, (-11) х6. Длина занятой части полосы M = Ах. Одна из особенностей блочного декодера в том, что размещая очередной предмет, он сначала проверяет все образовавшиеся ранее пустоты, подойдут ли они по размеру, для достижения большей плотности размещения (рис.1).

Рис. 2. Формальная запись края размещения

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

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

Для решения задачи предлагается модификация генетического алгоритма, использующая процедуру декодирования как параметр. Сам генетический алгоритм (т.е. схема его работы, смена поколений), структура хромосомы, операторы мутации, скрещивания, селекции не чувствительны к геометрии размещаемых предметов. Например, оператор мутации случайным образом переставляет два номера предметов в хромосоме. Оператор скрещивания двух хромосом

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

Например, оператор скрещивания (для случайных параметров р = 3, q = 4) для родителей Sp1 = (1, 2, 3,

4, 5, 6, 7, 8, 9) и Spi = (4, 7, 9,1, 2, 5, 3, 8, 6) порождает потомка Sc1 = (3, 4, 5, 6, 7, 9, 2, 8) и для тех же родителей (для удобства понимания по-другому записанных) Sp2 = (4, 7, 9,1, 5, 3, 8, 6) и Sp2 = (1, 2, 3, 4,

5, 6, 7, 8, 9) - потомка Sc2 = (9, 1, 2, 5, 3, 4, 6, 7, 8). Оператор мутации (для случайных параметров а = 2, b = 7) исходную особь S = (9, 8, 7, 6, 5, 4, 3, 2, 1) преобразует в особь вида Sm = (9,3, 7, 6, 5,4, 8, 2, 1).

Г енетический алгоритм показал, что оценочная полиномиальная сложность алгоритма и относительно невысокие времена работы программы не зависят от типа задачи. В то же время оценочная сложность блочного декодера пропорциональна квадрату размерности задачи (числу предметов), а краевого — линейно зависит от размерности.

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

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

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

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

6, Выводы

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

РИ, 2013, № 4

95

системе управления запасами. Исследованы методики поддержки принятия решений в процессах управления поступлением, хранением ресурсов, основанных на предложенной модели расчета затрат. Разработана методика принятия решений в целях оптимизации системы управления запасами материально-технических ресурсов, используемых на предприятиях. Проведен анализ эффективности разработанной методики управления запасами;

Научная новизна исследования заключается в обосновании экономического механизма управления запасами материально-технических ресурсов и разработке методического обеспечения процесса принятия решений по оптимизации объемов запасов.

Основными результатами исследования являются: обоснование места, роли и рациональной структуры запасов в системе планирования сбытовой деятельности предприятия в современных условиях; разработка концептуальных подходов к решению оптимизационных задач в области управления запасами материально-технических ресурсов, использующих инструментарий анализа и процесса принятия решений; определение факторов стратегии управления запасами, основанных на различных параметрах системы закупок, таких как периодичность заказа, постоянная партия объема пополнения заказа, уровень восполнения запаса. Предложена методика анализа состояния запасов, основанная на анализе изменения структуры и объема затрат предприятия, позволяющая контролировать состояние запасов в реальном режиме времени.

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

Литература: 1. Гвоздинский А.Н. Методы аналитической обработки информации / А.Н. Г воздинский, Е.Г. Климко // Радиоэлектроника и информатика. 2000. №4. С.111-112. 2. Гвоздинский А.Н. Принципы и методы оптимизации в современных организационных системах управления // Х: Бионика интеллекта, 2013. №2(81). С.66-70. 3. Гвоздинсь-кий А.М. Методи оптимізації в системах прийняття рішень // Х: ХНУРЕ, 2006. С.327.

Поступила в редколлегию 23.11.2013

Рецензент: д-р техн. наук, проф. Куземин А.Я.

Гвоздинский Анатолий Николаевич, канд. техн. наук, профессор кафедры искусственного интеллекта ХНУРЭ. Научные интересы: оптимизация процедур принятия решений в сложных системах управления. Адрес: Украина, 61166, Харьков, ул. Академика Ляпунова, 7, кв.9, тел. 702-38-23.

Козлова Алина Евгеньевна, студентка группы КН-10-4 кафедры искусственного интеллекта ХНУРЭ. Научные интересы: методы принятия решения в системах искусственного интеллекта. Адрес: Украина, 61000, Харьков, пр. Героев Сталинграда, кв.47, тел. (050)-61 -30-071.

96

РИ, 2013, № 4