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

Метод разработан в 1940 г. академиком Л. В. Конторовичем. В 1951 г. американским ученым Дж. Б. Данцигом предложен распределительный метод (МОДИ), аналогичный методу потенциалов. В обоих методах при проверке допустимого плана на оптимальность определяются потенциалы (числа), с помощью которых вычисляются характеристики клеток без кружков (в них нет поставок).

Обозначив потенциалы строк через U_i, потенциалы столбцов через V_j, показатели C_ij в клетках с поставками и кружками через , характеристики клеток без кружков (без поставок) через E_ij, получим следующие соотношения::

метод потенциалов метод МОДИ

u_i = v_j – ; (8.10) u_i = – v_j;

v_j = + u_i; (8.11) v_j = – u_i;

= v_j – u_i; (8.12) = v_j + u_i;

E_ij = c_ij – (v_j – u_i); (8.13) E_ij = c_ij – (v_j + u_i).

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

Рассмотрим пример решения транспортной задачи, предложенный В. С. Михеевой (1981). Базисный допустимый план составлен способом наименьшего элемента в столбце, его первоначальный функционал 555 (табл. 8.4).

Вначале рассчитаем потенциалы строк и столбцов по методу потенциалов с использованием формул. Произвольно выбранную величину потенциала выбирает в том столбце или строке, где наибольшее количество клеток с кружками. В нашем примере – это третья строка. В качестве потенциала для нее возьмем число 0. По формуле (8.11) определяем потенциалы (v_j) первого (v_j = + u_i = 0 + 2 = 2), третьего (0 + 6 = 6), четвертого (0 + 4 = 4), пятого (0 + 2 = 2) столбцов. Зная потенциалы по четырем столбцам, можно вычислить потенциалы строк (u_i) по формуле (8.10): первой (u_i = v_j – = 2 – 1 = 1), четвертой (6 – 2 = 4). По полученным потенциалам новых строк вычисляем потенциалы новых столбов, а по ним – новых строк (см. табл. 8.4).

Определив потенциалы строк и столбцов, вычисляют характеристики клеток (E_ij) без кружков (в них нет поставок) по формуле (8.13). Приведем расчет характеристики клетки 1.2: E_ij = c_ij – (v_j – u_i) = 2 – (7 – 1) = – 4. Аналогично рассчитываем E_ij (показаны курсивом) для других клеток без кружков.

Таблица 8.4

Базисный допустимый план (матрица 1)

bj ai	50	85	35	25	20	ui
30	 30	2 –4	3 –2	4 1	5 4	1
40	4 6	 40	2 0	3 3	4 6	4
70	 20	5 –2	 5	 25	 20	0
75	4 6	 45	 30	1 1	1 3	4
vj	2	7	6	4	2

Среди вычисленных характеристик клеток (курсив в клетках) отрицательные величины получены в клетках матрицы 1.2, 1.3, 3.2 (их величины соответственно: – 4, –2, –2), поэтому составленный первичный базовый план не оптимален. Проводят следующее перераспределение поставок по правилам цепи.

Выбираем клетку с наибольшей отрицательной абсолютной величиной характеристики (E₁₂) равную – 4. К клетке 1.2 строится цепь по перемещению наименьшей поставки 5 из клетки 3.3, так как функционал стремится к минимуму. Путь перемещения следующий: 3.3 (–5)4.3 (+5)4.2 (–5)1.2 (+5)1.1 (–5)3.1 (+5)3.3. К имеющейся поставке в клетке прибавляется или отнимается 5 с целью сохранения баланса между поставщиками и потребителями.

Получаем новую матрицу с измененными поставками (8.5). В ней повторяем алгоритм расчетов, как в табл. 8.4: рассчитываем потенциалы строк и столбцов, характеристику клеток без поставок, производим перераспределение поставок с использованием другой минимальной поставки 25 в клетке 3.4. Другая минимальная поставка 25 в клетке 3.1 перемещаться не может по цепи, так как в свободной клетке 4.4 с максимальной отрицательной абсолютной характеристикой (–3) ее следует вычитать из несуществующей поставки.

Таблица 8.5