-    

-    

-    

1.2       Метод северо-западного угла

1.    

2.      транспортной задачи.

3.     r=m+n-1.

4.     L=

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

Цикл в транспортной таблице – это ломаная с вершинами в клетках и звеньями, лежащих вдоль строк или столбцов  удовлетворяющая следующим требованиям:

-

-

Означенный цикл – цикл вершинам которого приписаны «+» и «-­­» поочередно. При переносе по означенному циклу k единиц перевозки в положительных вершинах добавляем  k единиц, а в отрицательных вершинах отнимаем k единиц. При таком переносе равновесие между запасами и заявками не нарушаются, следует план остается допустимым.

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

1.3       Распределительный метод или метод последовательного улучшения плана перевозок

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

Циклом пересчета данной свободной клетки называется цикл одна вершинка которого находится в этой свободной клетке, а остальные в базисных клетках. Для любой свободной клетки транспортной таблицы существует единственный цикл пересчета. Свободной клетке цикла пересчета присваивается знак «+». Если цена цикла пересчета некоторой свободной клетки отрицательна, то по этому циклу следует перенести   количество груза равного минимальному из перевозок в отрицательных вершинах.

1.4       Метод  потенциалов

Метод потенциалов позволяет, исходя из некоторого опорного плана, построить за конечное число итераций решение Т-задачи.

Метод потенциалов впервые предложили Л. В. Канторович и М. К. Гавурин в 1949 г. [18; 59]. Позже аналогичный метод разработал Г. Данциг, исходя из общих идей ЛП.

Общая схема метода такова. В данном начальном опорном плане перевозок каждому пункту ставят в соответствие некоторое число, называемое его предварительным потенциалом. Предварительные потенциалы выбирают так, чтобы их разность для любой пары пунктов A_i i B_j, связанных основной коммуникацией, была равна c_ij. Если окажется, что разность предварительных потенциалов для всех других коммуникаций не превосходит c_ij, то данный план перевозок - оптимальное решение задачи. В противном случае указывают способ улучшения текущего плана Т-задачи.

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

На предварительном этапе строят начальный опорный план и составляют матрицу

где  - предварительные потенциалы пунктов

 .

Предварительный этап. С помощью известного метода (например северо-западного угла или минимального элемента) определяют начальный опорный план Х₀ и вычисляют предварительные потенциалы .

Вычисление предварительных потенциалов производят так. По найденному опорному плану Х₀ строят схему перевозок Т-задачи из основных коммуникаций плана. Напомним, что основные коммуникации  плана Х₀ =  - это те, которым отвечают базисные компоненты плана, т.е. коммуникации для которых . Далее образуют следующие множества: J₁ - множество индексов всех пунктов B_j, которые связаны с пунктом А₁ основными коммуникациями; І₁ - множество индексов тех пунктов А_і, которые связаны основными коммуникациями с множеством J₁; J₂ - множество пунктов B_j, которые связаны основными коммуникациями с множеством І₁ и т.д. Образование таких множеств І_к продолжаем до тех пор, пока не получим  пустое множество.

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

Полагаем для определенности  и вычислим систему потенциалов относительно А₁. Тогда  где j J₁. Затем по значениям  определяем потенциалы пунктов  . Аналогично вычисляем потенциалы  (для и  .) и т.д. После того как потенциалы всех пунктов найдены, строим матрицу

Очевидно, позиции матрицы С₁, отвечающие базисным элементам плана Х₀, будут заняты нулями. Если матрица С₁ не содержит отрицательных элементов, то Х₀ - оптимальный план. В противном случае Х₀ - неоптимальный план, который может быть улучшен. Тогда переходим к выполнению однотипных итераций.

(k+1)-я итерация. Каждая итерация, кроме первой, где отсутствует первый этап, состоит из двух этапов. Предположим, что уже проведено k итераций (k=1,2,.),в результате которых получен план Х_k и вспомогательная матрицу С_k. Цель (k+1)-й итерации - построение матрицы С_k+1, а также либо установление оптимальности плана Х_k, либо нахождение более экономичного плана X_k+1.

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

Если все элементы матрицы С_k+1 окажутся неотрицательными, то X_k - оптимальный план, и на этом процесс заканчивается. В противном случае переходят ко второму этапу.

Второй этап. Цель этого этапа - построить более экономичный план Х_k+1. Выбирают наибольший по модулю отрицательный элемент матрицы С_k+1. Пусть это элемент . Строят цепочку из положительных элементов плана, которая замыкается на . После того, как цепочка построена, в ней находят минимальный нечетный по порядку следования элемент:

Прибавляют  ко всем четным элементам (по порядку следования) цепочки и к элементу и вычитают  из всех нечетных элементов. Остальные элементы Х_k оставляют без изменения.

Новый план Х_k+1 построен. Он является базисным, так как число его ненулевых элементов не изменилось.

Пусть L_k - транспортные издержки, отвечающие плану Х_k. Тогда новое значение целевой функции, отвечающее плану X_k+1, находят по соотношению

 .                             (3.2.1)

Так как  и , то . Поэтому Х_k+1 - улучшенный опорный план.

Затем производят аналогично (k+2)-ю итерацию.

Поставим в соответствие каждому пункту A_i  некоторое число  и каждому пункту назначения B_j  некоторое число

 и    называются потенциалами,  - это псевдостоимость.

В базисных  клетках c_ij=  перевозок является оптимальным если

-     c_ij=

-     c_ij, для всех свободных клеток.

Алгоритм метода потенциалов

1.     Строим  исходный оптимальный план, в котором  r=m+n-1 базисных клеток.

2.     Одной из неизвестных   присваиваем  произвольное численное значение (к примеру, 0) и по формуле  для базисных клеток находим потенциалы   и 

3.     Вычисляем псевдостоимости   для всех свободных клеток, если псевдостоимость ≤ стоимости ( c_ij), то план перевозок оптимальный.

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

5.     Подсчитываем новые  потенциалы.

2         ПОСТАНОВКА ЗАДАЧИ

2.1       Предметная область и общая постановка задачи

Объектом данной работы будет  отдел крупной торговой фирмы ООО «ТОНВИДЕО»( пер.Доломановский 183) в Ростове-на-Дону, который занимается распределением и сбытом  в Ростове-на-Дону инсекцицидное средство «КРА ДЕО СУПЕР» для уничтожение летающих насекомых, которое поставляется в Ростов-на-Дону из Казани (ул 3-я Кленовая 9)  железнодорожными путями.

Товар принимается в Ростове на трех складах: Можайская 167(в р-не авто рынка «Алмаз»), Врубова 32 и Доватора 44/3,  и уже оттуда распределяется на рынки: рынок «Лидер»(р-н александровка), «Нахичеванский», Ц.Рынок, «Привоз», «Военвед», «Темерник», в которых арендуются небольшие складские помещения специально под донный товар.

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

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

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

количеству упаковок  которые можно поместить на складе.

              «Лидер» - 30 

«Нахичеванский»-40

           «Ц.Рынок»-50

             «Привоз»-40

            «Военвед»-20

          «Темерник»-60

2.2       Математическая постановка задачи

Имеются  3 пункта отправления товара Можайская 167 (А₁), Врубова 32(А₂) и Доватора 44/3 (А₃), в которых сосредоточено 90, 80 и 80 упаковок соответственно, предназначенных для доставки, и 6 пунктов назначения: «Лидер» (В₁), «Нахичеванский», (В₂), Ц.Рынок (В₃), «Привоз»(В₄), «Военвед»(В₅), «Темерник» (В₆), которые подали заявки на некоторое количество товара, которое описано выше. Известны  время перевозки из каждого склада на каждый рынок.

Требуется составить план перевозок, при котором все заявки  были бы  удовлетворены и суммарное время перевозок  была бы минимальна.

Обозначим x_ij-количество товара, которое надо отправить из склада  на рынок. Тогда наша задача выглядит следующим образом L= min, где   j=(1,6), i=(1,3) (n=6, m=3). План перевозок x_ij, будет опорным, если в нем не равны нулю не более чем r=m+n-1 перевозок x_ij.Так как 90+80+80=30+40+50+40+20+60, следует транспортная задача закрытая.

Транспорт перевозит товар из  А₁ в В₁ за 20 минут