Транспортная задача с ограничениями на пропускную способность

Некоторый однородный продукт, сосредоточенный у m поставщиков Ai в количестве ai (i = 1, 2, …, m) у каждого, необходимо доставить к потребителю Bj в количестве bj (j = 1, 2, …, n). Известна стоимость Cij перевозок единицы груза от i - го поставщика j - му потребителю и ограничения по пропускной способности dij.

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

Математическая модель задачи такова

при ограничениях:

,

где dij - предельное число единиц продукции, которое может быть перевезено от i- го поставщика j – му потребителю за время, указанное в условии задачи. Такая задача называется ТЗ с ограничениями на пропускную способность (Td). Для решения таких задач может быть приспособлен любой из вычислительных алгоритмов решения обычной ТЗ, например, венгерский метод, метод потенциалов.

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

Венгерский метод наиболее эффективен при решении ТЗ с целочисленными объемами производства и потребления, в этом случае число итераций не превышает величины D0/2, где D0 - суммарная невязка подготовительного этапа.

На подготовительном этапе строится матрица x0 =(xij[0])mn, элементы которой неотрицательны и удовлетворяют неравенствам:

.

Если эти условия являются равенствами, то матрица Х0- решение транспортной задачи. Если среди условий имеются неравенства, то осуществляется переход к первой итерации. На К- ой итерации строится матрица Хк = (Хij[К])mn, суммарная невязка которой DК

В результате первой итерации строится матрица x1, состоящая из неотрицательных элементов. При этом D1<D0. Если D1 = 0, то x1 - оптимальное решение задачи. Если D1 > 0, то переходят к следующей итерации. Они проводятся до тех пор, пока Dкпри некотором К не станет равным 0 . Соответствующая матрица xk является решением ТЗ.

Достоинством венгерского метода является возможность оценивать близость результата каждой итерации к оптимальному плану перевозок, что позволяет контролировать процесс вычислений и прекратить его при достижении определенных точностных показателей. Особенно это свойство существенно для задач большой размерности. Подробное изложение алгоритма см. в книге: Гольштейн Е.Г., Юдин Д.Б. Задачи линейного программирования транспортного типа.- М.: Наука, 1969. - 382 с..