Математические модели в транспортных системах-1

Для оценки оптимальности текущего решения - поставок ресурса от потребителей поставщикам применяется ряд методов: распределительные методы (метод Хичкока, метод Креко, метод МОДИ), методы с разрешающими элементами (метод разрешающих слагаем^хх и метод разрешающих множителей) и др.

Наиболее широкое применение получил модифицированный метод (МОДИ). В распределительную таблицу вводят строку и столбец для расчета потенциалов. Метод, как и другие распределительные методы, основан на том, что если к стоимостям любой строки (столбца) распределительной таблицы прибавить (или отнять) одно и то же произвольное число, то оценка оптимальности относительно не изменится. Если, например, от стоимостей поставок от каждого 1-го поставщика отнять число uAi , а от стоимостей поставок каждому j- го потребителю - uBj, то относительной оценкой любой связи ij может служить оценочный

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

потребителей таким образом, чтобы соблюдалось принятое условие l.j - uAi- uBj = 0, т.е. стоимость для каждой загруженной связи должна быть равна сумме потенциалов по поставщику и потребителю.

Затем по вычисленным потенциалам uAi , uBj определяется, используя формулу (3.6), значение оценочного параметра s.j для каждой незагруженной связи (не входящей в базисный план). Величина параметра s.j характеризует общее изменение целевой функции, получаемое при включении в план единичной корреспонденции для связи ij по сравнению с рассматриваемым планом.

Если значение оценочного параметра свободной связи будет меньше нуля (sij <0), то это значит, что перераспределение корреспонденций с загрузкой такой связи, называемой потенциальной, уменьшит значение целевой функции на величину abs(sijXij), где X.j - размер корреспонденции, которой будет загружена связь ij. Отсутствие связей со значением параметра sij <0 означает, что проверяемый план закрепления потребителей за поставщиками является оптимальным. Если для какой-то свободной связи значение s.j равно нулю, то это означает, что можно получить другой план с тем же минимальным значением целевой функции.

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

При перераспределении загрузок по связям в распределительной таблице для наиболее потенциальной связи, как для загруженной, строится контур. Введение к m+n-1 загруженным связям еще одной образует ровно один определенный контур, присущий вводимой связи. После чего связи, соответствующие вершинам контура, нумеруются: номер 1 присваивается выбранной потенциальной связи, а дальнейшая нумерация ведется в порядке следования вершин по контуру (могут присваиваться связям контура положительные и отрицательные знаки).

Затем производится перераспределение по контуру корреспонденций следующим образом:

1) выявляется связь с четным номером, которой соответствует наименьшее значение корреспонденции Xм ;

91

2)значение объема ресурса Хм вычитается от значений объемов корреспонденций всех связей с четными номерами вершин. Если среди связей с четными номерами вершин окажутся две (или более) с одинаковой минимальной корреспонденцией, то из плана исключается только одна из них, для связи с большим расстоянием поставки, а вместо остальных оставляют условную корреспонденцию с нулевым значением, чтобы не допустить вырождения;

3)для всех связей с нечетными номерами вершин (включая и потенциальную) к значениям объемов корреспонденций прибавляется величина Хм.

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

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

При решении транспортной задачи линейного программирования на основе существующего базисного решения возможны варианты, когда число загруженных связей в таблице меньше m+n-1 (случай вырождения) или больше m+n-1 (случай наличия циклов). В первом случае нельзя оценить оптимальность решения, так как невозможно определить потенциалы для некоторых поставщиков и (или) потребителей, во втором - потенциалы определяются неоднозначно.

Вслучае вырождения необходимое число свободных связей загружают нулями с соблюдением правила: нулевые объемы поставок заносятся для связи поставщиков (потребителей), для которых нельзя рассчитать потенциал, с поставщиками (потребителями), для которых потенциалы уже определены. Целесообразно выбрать из данн^хх связей такие, которым соответствуют наименьшие стоимости поставок. Связи, загруженные нулями, не должны образовывать контуров (циклов) с другими загруженными связями.

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

1)строится один из возможн^1х контуров;

2)нумеруются углы контура;

3)определяется сумма стоимостей поставок для связей, в которых лежат углы контура с четными номерами и затем с нечетными номерами;

4)среди связей (с четными или нечетными номерами), для которых сумма стоимостей перевозок больше, выявляется связь с наименьшим значением размера корреспонденции Хм;

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

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

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

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

1) проверяется выполнение условия Xy <Xoij для связей ij. Если условие выполняется для всех связей, то решение закончено (переход на блок 10 общего алгоритма), иначе переход на пункт 2;

92

3)поставка транспортных средств 2-го предприятия на 2-й маршрут должна быть не менее 5 единиц;

4)поставка транспортных средств 1-го предприятия на 4-й маршрут запрещена. Требуется оптимальным образом распределить транспортные средства по маршрутам

перевозок. Расстояния между предприятиями и маршрутами (нулевые пробеги транспортных средств) приведены ниже в таблице 3.4.

Таблица 3.4 - Данные по нулевым пробегам транспортн^1х средств

Решение.

Проверяем задачу на сбалансированность, для чего рассчитываем суммарный объем наличия ресурса ХсА и суммарный объем потребности в ресурсе ХсВ

После этого сравниваем объемы ХСА и ХСВ. Поскольку они не равны, то задача несбалансирована. Так как ХсА >ХсВ, то вводится фиктивный потребитель В5 (n=n+1) и по не нему назначается объем ресурса в размере ХВ5 = abs(XcА - ХсВ)=80-70=10. Для фиктивного потребителя назначаются стоимости (расстояния) равными постоянной величине, например, 1ф=10 max 1ij = 8*10=80(1i5 = 80, 1=1,2,3). Однако в связи с тем, что транспортные средства 2-го

1, J J

предприятия необходимо использовать в полном объеме, то расстояние 12-5 следует принять значительно больше других, например, 12-5=800. Затем по фиктивному пункту назначаются ограничения равными постоянной величине, например, Х0ф > abs(XcА - ХсВ) > 10 (для фиктивного потребителя Xoi5 = 10, 1 = 1, 2, 3).

Для учета запрещения поставки ресурса от поставщика А1 потребителю В4 заменяем реальную стоимость (расстояние) 11.4 на большое число, например 114 > 100max 1ij > 800.

i,j

Принимаем 11-4 = 800.

Обеспечение поставки от поставщика А2 потребителю В2 в гарантированном объеме ХГ2 - 2=5 производится следующим образом:

в окончательное решение вводится корреспонденция X2-2 = 5; корректируются соответствующие размеры наличия и потребности в ресурсе

Х А 2 = Х А 2 - X 2 . 2 =30-5=25; Х В 2 = Х В 2 - X 2 . 2 =10-5=5.

С учетом исходных условий и последующих изменений составляем распределительную таблицу (таблица 3.5) и производим первоначальное базисное распределение по методу минимального элемента. Для этого назначаем корреспонденции для связей в порядке возрастания по ним стоимостей (расстояний). Связи загружаются в следующей последовательности: А1-В3 объемом 10, А3-В2 объемом 5, А3-В3 объемом 10; А3-В4 объемом 25, А2-В1 объемом 10, А2-В4 объемом 5, А2-В5 объемом 10.

Полученное первоначальное распределение проверяем на оптимальность, для чего рассчитываем значения потенциалов uAi и uBj. При этом значение одного из потенциалов, например для Aj, задаем равным нулю (UAJ =0). После того, как определены значения потенциалов, для каждой свободной связи рассчитываем оценочный параметр sij по формуле (3.6).

94

Таблица 3.8 - Распределительная таблица (3-е улучшенное базисное распределение)

Согласно укрупненному алгоритму решения задачи следует с новыми данными вернуться на составление первоначального базисного решения (таблица 3.9).

Производим проверку полученного распределения на оптимальность. Так как в данном распределении для отдельн^хх свободных связей оценочные параметры s.j < 0, то данное решение является неоптимальным.

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

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

97

Полное окончательное решение содержит:

ранее включенные корреспонденции в окончательное решение

X2-2 = 5;

X1.3 = 5;

и последнее решение для реальных пунктов

X j - j = 5; X2-1 = 5; X2-2 = 5; X2-3 = 15; X3.4 = 30.

Ответ: исходя из минимума нулевых пробегов и других дополнительных условий оптимально закрепить за 1-м предприятием 1-й и 3-й маршрут с выделением по 5 единиц, за 2-м - 1-й, 2-й и 3-й с выделением соответственно 5, 10 и 15 единиц, за 3-м - 4-й маршрут с выделением 30 единиц транспортных средств.

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

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

Одной из наиболее простых является следующая нелинейная однопродуктовая распределительная задача:

необходимо оптимальным образом распределить ресурс в количестве Хсо по n вариантам (объектам, схемам, этапам) при целевой функции

98

Z Xi < X со ,

1=1

где X. - количество ресурса, назначаемое по 1-му варианту (0 <xi < xco);

f(x) - нелинейная функция, определяющая эффект (затраты) в зависимости от значений x; n - общее число вариантов;

xco - общее количество распределяемого ресурса.

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

Обозначим через Xc. - количество ресурса, назначаемого суммарно по 1-му варианту и всем предыдущим этапам (вариантам). При этом количество ресурса, входящее в xci, распределено по предыдущим вариантам оптимально исходя из минимума затрат или максимума эффекта.

Тогда целевая функция Z. при рассмотрении 1-го этапа решения определяется по выражению

Zi(Xci, Xi)= fi(Xi) + f*i-i (Xci - Xi),

где 0 <xci < xco при i < n-1 и xci= xco при i = n;

f*i - i - функция, определяющая эффект (затраты) оптимальным образом по предыдущим вариантам в зависимости от количества ресурса

Значение f*. определяется по Z. :

f l * ( x c i ) = f 1 ( x c i ) 1 = 1 ( x c i = x i ) ;

fi (xci) = max (min)Zi(xc1, x . ) при х. , i = 2, n

Xl

Определение fi*{xci) и соответствующих им х* при данном xci производятся поэтапно от i=2 до i=n.

По результатам завершенных поэтапных расчетов формируется оптимальное решение

{ X o п т i } :

Вариант алгоритма решения однопродуктовой задачи динамического программирования приведен на рисунке 3.25. Решение предусмотрено для целевой функции, подлежащей оптимизации на максимум. Для решения задачи на минимум необходимо ввести исходные значения s^ij. с отрицательным знаком или внести изменения в блоках 12 (заменить 1E+12 на -1E+12) и 15 (условие < на >).

Выводимые j, kj.dx , snm представляют соответственно номер варианта, количество ресурса по варианту и значение целевой функции для оптимального решения.

99