Размещено на http://www.allbest.ru/

Экономико-математические методы

1. Задача 1

перевозка продукция расход потребитель

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

Виды сырья	Расход сырья на единицу продукции	Общий запас сырья, ед.
	М1	М2	М3
П1	2	3	4	273
П2	1	4	3	200
П3	3	1	2	303
Уровень прибыли на ед. продукции	20	28	24

Решение.

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

Определим максимальное значение целевой функции

F(X) = 20x1 + 28x2 + 24x3 при следующих условиях-ограничениях:

2x1 + 3x2 + 4x3?273

x1 + 4x2 + 3x3?200

3x1 + x2 + 2x3?303

Для построения первого опорного плана систему неравенств приведем к системе уравнений путем введения дополнительных переменных (переход к канонической форме): В 1-м неравенстве смысла (?) вводим базисную переменную x4. В 2-м неравенстве смысла (?) вводим базисную переменную x5. В 3-м неравенстве смысла (?) вводим базисную переменную x6.

2x1 + 3x2 + 4x3 + 1x4 + 0x5 + 0x6 = 273

1x1 + 4x2 + 3x3 + 0x4 + 1x5 + 0x6 = 200

3x1 + 1x2 + 2x3 + 0x4 + 0x5 + 1x6 = 303

Матрица коэффициентов A = a(ij) этой системы уравнений имеет вид:

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

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

Решим систему уравнений относительно базисных переменных: x4, x5, x6

Полагая, что свободные переменные равны 0, получим первый опорный план:

X1 = (0,0,0,273,200,303)

Базисное решение называется допустимым, если оно неотрицательно.

Базис	B	x1	x2	x3	x4	x5	x6
x4	273	2	3	4	1	0	0
x5	200	1	4	3	0	1	0
x6	303	3	1	2	0	0	1
F(X0)	0	-20	-28	-24	0	0	0

Переходим к основному алгоритму симплекс-метода.

Итерация №0.

1. Проверка критерия оптимальности.

Текущий опорный план неоптимален, так как в индексной строке находятся отрицательные коэффициенты.

2. Определение новой базисной переменной.

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

3. Определение новой свободной переменной.

Вычислим значения Di по строкам как частное от деления: bi / ai2 и из них выберем наименьшее: min (273 : 3 , 200 : 4 , 303 : 1 ) = 50

Следовательно, 2-ая строка является ведущей. Разрешающий элемент равен (4) и находится на пересечении ведущего столбца и ведущей строки.

Базис	B	x1	x2	x3	x4	x5	x6	min
x4	273	2	3	4	1	0	0	91
x5	200	1	4	3	0	1	0	50
x6	303	3	1	2	0	0	1	303
F(X1)	0	-20	-28	-24	0	0	0	0

4. Пересчет симплекс-таблицы.

Формируем следующую часть симплексной таблицы.

Вместо переменной x5 в план 1 войдет переменная x2.

Строка, соответствующая переменной x2 в плане 1, получена в результате деления всех элементов строки x5 плана 0 на разрешающий элемент РЭ=4 На месте разрешающего элемента в плане 1 получаем 1.

В остальных клетках столбца x2 плана 1 записываем нули.

Таким образом, в новом плане 1 заполнены строка x2 и столбец x2.

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

НЭ = СЭ - (А*В)/РЭ

СТЭ - элемент старого плана, РЭ - разрешающий элемент (4), А и В - элементы старого плана, образующие прямоугольник с элементами СТЭ и РЭ.

Представим расчет каждого элемента в виде таблицы:

B	x 1	x 2	x 3	x 4	x 5	x 6
273-(200 * 3):4	2-(1 * 3):4	3-(4 * 3):4	4-(3 * 3):4	1-(0 * 3):4	0-(1 * 3):4	0-(0 * 3):4
200 : 4	1 : 4	4 : 4	3 : 4	0 : 4	1 : 4	0 : 4
303-(200 * 1):4	3-(1 * 1):4	1-(4 * 1):4	2-(3 * 1):4	0-(0 * 1):4	0-(1 * 1):4	1-(0 * 1):4
0-(200 * -28):4	-20-(1 * -28):4	-28-(4 * -28):4	-24-(3 * -28):4	0-(0 * -28):4	0-(1 * -28):4	0-(0 * -28):4

Получаем новую симплекс-таблицу:

Базис	B	x1	x2	x3	x4	x5	x6
x4	123	5/4	0	7/4	1	-3/4	0
x2	50	1/4	1	3/4	0	1/4	0
x6	253	11/4	0	5/4	0	-1/4	1
F(X1)	1400	-13	0	-3	0	7	0

Итерация №1.

1. Проверка критерия оптимальности.

Текущий опорный план неоптимален, так как в индексной строке находятся отрицательные коэффициенты.

2. Определение новой базисной переменной.

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

3. Определение новой свободной переменной.

Вычислим значения Di по строкам как частное от деления: bi / ai1

и из них выберем наименьшее: min (123 : 11/4 , 50 : 1/4 , 253 : 23/4 ) = 92

Следовательно, 3-ая строка является ведущей.

Разрешающий элемент равен (23/4) и находится на пересечении ведущего столбца и ведущей строки.

Базис	B	x1	x2	x3	x4	x5	x6	min
x4	123	11/4	0	13/4	1	-3/4	0	982/5
x2	50	1/4	1	3/4	0	1/4	0	200
x6	253	23/4	0	11/4	0	-1/4	1	92
F(X2)	1400	-13	0	-3	0	7	0	0

4. Пересчет симплекс-таблицы.

Формируем следующую часть симплексной таблицы. Вместо переменной x6 в план 2 войдет переменная x1. Строка, соответствующая переменной x1 в плане 2, получена в результате деления всех элементов строки x6 плана 1 на разрешающий элемент РЭ=23/4

На месте разрешающего элемента в плане 2 получаем 1.

В остальных клетках столбца x1 плана 2 записываем нули.

Таким образом, в новом плане 2 заполнены строка x1 и столбец x1.

Все остальные элементы нового плана 2, включая элементы индексной строки, определяются по правилу прямоугольника.

Представим расчет каждого элемента в виде таблицы:

B	x 1	x 2	x 3	x 4	x 5	x 6
123-(253 * 11/4):23/4	11/4-(23/4 * 11/4):23/4	0-(0 * 11/4):23/4	13/4-(11/4 * 11/4):23/4	1-(0 * 11/4):23/4	-3/4-(-1/4 * 11/4):23/4	0-(1 * 11/4):23/4
50-(253 * 1/4):23/4	1/4-(23/4 * 1/4):23/4	1-(0 * 1/4):23/4	3/4-(11/4 * 1/4):23/4	0-(0 * 1/4):23/4	1/4-(-1/4 * 1/4):23/4	0-(1 * 1/4):23/4
253 : 23/4	23/4 : 23/4	0 : 23/4	11/4 : 23/4	0 : 23/4	-1/4 : 23/4	1 : 23/4
1400-(253 * -13):23/4	-13-(23/4 * -13):23/4	0-(0 * -13):23/4	-3-(11/4 * -13):23/4	0-(0 * -13):23/4	7-(-1/4 * -13):23/4	0-(1 * -13):23/4

Получаем новую симплекс-таблицу:

Базис	B	x1	x2	x3	x4	x5
x4	8	0	0	13/11	1	-7/11
x2	27	0	1	7/11	0	3/11
x1	92	1	0	5/11	0	-1/11
F(X2)	2596	0	0	32/11	0	64/11

1. Проверка критерия оптимальности.

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

Окончательный вариант симплекс-таблицы:

Базис	B	x1	x2	x3	x4	x5
x4	8	0	0	13/11	1	-7/11
x2	27	0	1	7/11	0	3/11
x1	92	1	0	5/11	0	-1/11
F(X3)	2596	0	0	32/11	0	64/11

Оптимальный план можно записать так:

x2 = 27

x1 = 92

F(X) = 28*27 + 20*92 = 2596

Анализ оптимального плана.

В оптимальный план вошла дополнительная переменная x4. Следовательно, при реализации такого плана имеются недоиспользованные ресурсы 1-го вида в количестве 8

Значение 0 в столбце x1 означает, что использование x1 - выгодно.

Значение 0 в столбце x2 означает, что использование x2 - выгодно.

Значение 210/11> 0 в столбце x3 означает, что использование x3 - не выгодно.

Значение 59/11 в столбце x5 означает, что теневая цена (двойственная оценка) равна 59/11.

Значение 48/11 в столбце x6 означает, что теневая цена (двойственная оценка) равна 48/11.

2. Задача 2

Составить оптимальный план перевозок пищевых продуктов от 4-х поставщиков к 6-ти потребителям. Поставщики (П), потребители (М), объемы вывоза и завоза, кратчайшие расстояния между пунктами вывоза и завоза приведены в таблицах.

Поставщики	Потребители	Объемы вывоза, тонн
	М1	М2	М3	М4	М5	М6
П1	32	9	24	28	30	27	148
П2	40	24	22	23	20	21	76
П3	15	25	30	21	39	38	144
П4	39	11	36	8	27	30	132
Объемы завоза, тонн	112	92	84	36	80	96

Решение.

Математическая модель транспортной задачи:

F = ??cijxij, (1)

при условиях:

?xij = ai, i = 1,2,…, m, (2)

?xij = bj, j = 1,2,…, n, (3)

xij >= 0

Запишем экономико-математическую модель для нашей задачи.

Переменные:

x11 - количество груза из 1-го склада в 1-й магазин

x12 - количество груза из 1-го склада в 2-й магазин

x13 - количество груза из 1-го склада в 3-й магазин

x14 - количество груза из 1-го склада в 4-й магазин

x15 - количество груза из 1-го склада в 5-й магазин

x16 - количество груза из 1-го склада в 6-й магазин

x21 - количество груза из 2-го склада в 1-й магазин

x22 - количество груза из 2-го склада в 2-й магазин

x23 - количество груза из 2-го склада в 3-й магазин

x24 - количество груза из 2-го склада в 4-й магазин

x25 - количество груза из 2-го склада в 5-й магазин

x26 - количество груза из 2-го склада в 6-й магазин

x31 - количество груза из 3-го склада в 1-й магазин

x32 - количество груза из 3-го склада в 2-й магазин

x33 - количество груза из 3-го склада в 3-й магазин

x34 - количество груза из 3-го склада в 4-й магазин

x35 - количество груза из 3-го склада в 5-й магазин

x36 - количество груза из 3-го склада в 6-й магазин

x41 - количество груза из 4-го склада в 1-й магазин

x42 - количество груза из 4-го склада в 2-й магазин

x43 - количество груза из 4-го склада в 3-й магазин

x44 - количество груза из 4-го склада в 4-й магазин

x45 - количество груза из 4-го склада в 5-й магазин

x46 - количество груза из 4-го склада в 6-й магазин

Ограничения по запасам:

x11 + x12 + x13 + x14 + x15 + x16 <= 148

x21 + x22 + x23 + x24 + x25 + x26 <= 76

x31 + x32 + x33 + x34 + x35 + x36 <= 144

x41 + x42 + x43 + x44 + x45 + x46 <= 132

Ограничения по потребностям:

x11 + x21 + x31 + x41 >= 112

x12 + x22 + x32 + x42 >= 92

x13 + x23 + x33 + x43 >= 84

x14 + x24 + x34 + x44 >= 36

x15 + x25 + x35 + x45 >= 80

x16 + x26 + x36 + x46 >= 96

Целевая функция:

32x11 + 9x12 + 24x13 + 28x14 + 30x15 + 27x16 + 40x21 + 24x22 + 22x23 + 23x24 + 20x25 + 21x26 + 15x31 + 25x32 + 30x33 + 21x34 + 39x35 + 38x36 + 39x41 + 11x42 + 36x43 + 8x44 + 27x45 + 30x46 > min

Стоимость доставки единицы груза из каждого пункта отправления в соответствующие пункты назначения задана матрицей тарифов

	1	2	3	4	5	6	Запасы
1	32	9	24	28	30	27	148
2	40	24	22	23	20	21	76
3	15	25	30	21	39	38	144
4	39	11	36	8	27	30	132
Потребности	112	92	84	36	80	96

Проверим необходимое и достаточное условие разрешимости задачи.

?a = 148 + 76 + 144 + 132 = 500

?b = 112 + 92 + 84 + 36 + 80 + 96 = 500

Условие баланса соблюдается. Запасы равны потребностям. Следовательно, модель транспортной задачи является закрытой.

Занесем исходные данные в распределительную таблицу.

	1	2	3	4	5	6	Запасы
1	32	9	24	28	30	27	148
2	40	24	22	23	20	21	76
3	15	25	30	21	39	38	144
4	39	11	36	8	27	30	132
Потребности	112	92	84	36	80	96

Этап I. Поиск первого опорного плана.

1. Используя метод наименьшей стоимости, построим первый опорный план транспортной задачи.

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

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

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

Искомый элемент равен 8

Для этого элемента запасы равны 132, потребности 36. Поскольку минимальным является 36, то вычитаем его.

x44 = min(132,36) = 36.

32	9	24	x	30	27	148
40	24	22	x	20	21	76
15	25	30	x	39	38	144
39	11	36	8	27	30	132 - 36 = 96
112	92	84	36 - 36 = 0	80	96	0

Искомый элемент равен 9

Для этого элемента запасы равны 148, потребности 92. Поскольку минимальным является 92, то вычитаем его.

x12 = min(148,92) = 92.

32	9	24	x	30	27	148 - 92 = 56
40	x	22	x	20	21	76
15	x	30	x	39	38	144
39	x	36	8	27	30	96
112	92 - 92 = 0	84	0	80	96	0

Искомый элемент равен 15

Для этого элемента запасы равны 144, потребности 112. Поскольку минимальным является 112, то вычитаем его. x31 = min(144,112) = 112.

x	9	24	x	30	27	56
x	x	22	x	20	21	76
15	x	30	x	39	38	144 - 112 = 32
x	x	36	8	27	30	96
112 - 112 = 0	0	84	0	80	96	0

x25 = min(76,80) = 76.

x	9	24	x	30	27	56
x	x	x	x	20	x	76 - 76 = 0
15	x	30	x	39	38	32
x	x	36	8	27	30	96
0	0	84	0	80 - 76 = 4	96	0

Искомый элемент равен 24

Для этого элемента запасы равны 56, потребности 84. Поскольку минимальным является 56, то вычитаем его.

x13 = min(56,84) = 56.

x	9	24	x	x	x	56 - 56 = 0
x	x	x	x	20	x	0
15	x	30	x	39	38	32
x	x	36	8	27	30	96
0	0	84 - 56 = 28	0	4	96	0

Искомый элемент равен 27

Для этого элемента запасы равны 96, потребности 4. Поскольку минимальным является 4, то вычитаем его. x45 = min(96,4) = 4.

x	9	24	x	x	x	0
x	x	x	x	20	x	0
15	x	30	x	x	38	32
x	x	36	8	27	30	96 - 4 = 92
0	0	28	0	4 - 4 = 0	96	0

Искомый элемент равен 30

Для этого элемента запасы равны 32, потребности 28. Поскольку минимальным является 28, то вычитаем его. x33 = min(32,28) = 28.

x	9	24	x	x	x	0
x	x	x	x	20	x	0
15	x	30	x	x	38	32 - 28 = 4
x	x	x	8	27	30	92
0	0	28 - 28 = 0	0	0	96	0

Искомый элемент равен 30

Для этого элемента запасы равны 92, потребности 96. Поскольку минимальным является 92, то вычитаем его.

x46 = min(92,96) = 92.

x	9	24	x	x	x	0
x	x	x	x	20	x	0
15	x	30	x	x	38	4
x	x	x	8	27	30	92 - 92 = 0
0	0	0	0	0	96 - 92 = 4	0

Искомый элемент равен 38

Для этого элемента запасы равны 4, потребности 4. Поскольку минимальным является 4, то вычитаем его.

x36 = min(4,4) = 4.

x	9	24	x	x	x	0
x	x	x	x	20	x	0
15	x	30	x	x	38	4 - 4 = 0
x	x	x	8	27	30	0
0	0	0	0	0	4 - 4 = 0	0

	1	2	3	4	5	6	Запасы
1	32	9[92]	24[56]	28	30	27	148
2	40	24	22	23	20[76]	21	76
3	15[112]	25	30[28]	21	39	38[4]	144
4	39	11	36	8[36]	27[4]	30[92]	132
Потребности	112	92	84	36	80	96

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

2. Подсчитаем число занятых клеток таблицы, их 9, а должно быть m + n - 1 = 9.

Следовательно, опорный план является невырожденным.

Значение целевой функции для этого опорного плана равно:

F(x) = 9*92 + 24*56 + 20*76 + 15*112 + 30*28 + 38*4 + 8*36 + 27*4 + 30*92 = 9520

Этап II. Улучшение опорного плана.

Проверим оптимальность опорного плана. Найдем предварительные потенциалы ui, vj. по занятым клеткам таблицы, в которых ui + vj = cij, полагая, что u1 = 0.

u1 + v2 = 9; 0 + v2 = 9; v2 = 9

u1 + v3 = 24; 0 + v3 = 24; v3 = 24

u3 + v3 = 30; 24 + u3 = 30; u3 = 6

u3 + v1 = 15; 6 + v1 = 15; v1 = 9

u3 + v6 = 38; 6 + v6 = 38; v6 = 32

u4 + v6 = 30; 32 + u4 = 30; u4 = -2

u4 + v4 = 8; -2 + v4 = 8; v4 = 10

u4 + v5 = 27; -2 + v5 = 27; v5 = 29

u2 + v5 = 20; 29 + u2 = 20; u2 = -9

	v1=9	v2=9	v3=24	v4=10	v5=29	v6=32
u1=0	32	9[92]	24[56]	28	30	27
u2=-9	40	24	22	23	20[76]	21
u3=6	15[112]	25	30[28]	21	39	38[4]
u4=-2	39	11	36	8[36]	27[4]	30[92]

Опорный план не является оптимальным, так как существуют оценки свободных клеток, для которых ui + vj > cij

(1;6): 0 + 32 > 27; ?16 = 0 + 32 - 27 = 5

(2;6): -9 + 32 > 21; ?26 = -9 + 32 - 21 = 2

max(5,2) = 5

Выбираем максимальную оценку свободной клетки (1;6): 27

Для этого в перспективную клетку (1;6) поставим знак «+», а в остальных вершинах многоугольника чередующиеся знаки «-», «+», «-».

	1	2	3	4	5	6	Запасы
1	32	9[92]	24[56][-]	28	30	27[+]	148
2	40	24	22	23	20[76]	21	76
3	15[112]	25	30[28][+]	21	39	38[4][-]	144
4	39	11	36	8[36]	27[4]	30[92]	132
Потребности	112	92	84	36	80	96

Цикл приведен в таблице (1,6 > 1,3 > 3,3 > 3,6).

Из грузов хij стоящих в минусовых клетках, выбираем наименьшее, т.е. у = min (3, 6) = 4. Прибавляем 4 к объемам грузов, стоящих в плюсовых клетках и вычитаем 4 из Хij, стоящих в минусовых клетках. В результате получим новый опорный план.

	1	2	3	4	5	6	Запасы
1	32	9[92]	24[52]	28	30	27[4]	148
2	40	24	22	23	20[76]	21	76
3	15[112]	25	30[32]	21	39	38	144
4	39	11	36	8[36]	27[4]	30[92]	132
Потребности	112	92	84	36	80	96

Проверим оптимальность опорного плана. Найдем предварительные потенциалы ui, vj. по занятым клеткам таблицы, в которых полагая, что u1=0.

ui + vj = cij,

	v1=9	v2=9	v3=24	v4=5	v5=24	v6=27
u1=0	32	9[92]	24[52]	28	30	27[4]
u2=-4	40	24	22	23	20[76]	21
u3=6	15[112]	25	30[32]	21	39	38
u4=3	39	11	36	8[36]	27[4]	30[92]

u1 + v2 = 9; 0 + v2 = 9; v2 = 9

u1 + v3 = 24; 0 + v3 = 24; v3 = 24

u3 + v3 = 30; 24 + u3 = 30; u3 = 6

u3 + v1 = 15; 6 + v1 = 15; v1 = 9

u1 + v6 = 27; 0 + v6 = 27; v6 = 27

u4 + v6 = 30; 27 + u4 = 30; u4 = 3

u4 + v4 = 8; 3 + v4 = 8; v4 = 5

u4 + v5 = 27; 3 + v5 = 27; v5 = 24

u2 + v5 = 20; 24 + u2 = 20; u2 = -4

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

ui + vj > cij

(2;6): -4 + 27 > 21; ?26 = -4 + 27 - 21 = 2

(4;2): 3 + 9 > 11; ?42 = 3 + 9 - 11 = 1

max(2,1) = 2

Выбираем максимальную оценку свободной клетки (2;6): 21

Для этого в перспективную клетку (2;6) поставим знак «+», а в остальных вершинах многоугольника чередующиеся знаки «-», «+», «-».

	1	2	3	4	5	6	Запасы
1	32	9[92]	24[52]	28	30	27[4]	148
2	40	24	22	23	20[76][-]	21[+]	76
3	15[112]	25	30[32]	21	39	38	144
4	39	11	36	8[36]	27[4][+]	30[92][-]	132
Потребности	112	92	84	36	80	96

Цикл приведен в таблице (2,6 > 2,5 > 4,5 > 4,6).

Из грузов хij стоящих в минусовых клетках, выбираем наименьшее, т.е. у = min (2, 5) = 76. Прибавляем 76 к объемам грузов, стоящих в плюсовых клетках и вычитаем 76 из Хij, стоящих в минусовых клетках. В результате получим новый опорный план.

	1	2	3	4	5	6	Запасы
1	32	9[92]	24[52]	28	30	27[4]	148
2	40	24	22	23	20	21[76]	76
3	15[112]	25	30[32]	21	39	38	144
4	39	11	36	8[36]	27[80]	30[16]	132
Потребности	112	92	84	36	80	96

Проверим оптимальность опорного плана. Найдем предварительные потенциалы ui, vj. по занятым клеткам таблицы, в которых полагая, что u1=0.

ui + vj = cij,

u1 + v2 = 9; 0 + v2 = 9; v2 = 9

u1 + v3 = 24; 0 + v3 = 24; v3 = 24

u3 + v3 = 30; 24 + u3 = 30; u3 = 6

u3 + v1 = 15; 6 + v1 = 15; v1 = 9

u1 + v6 = 27; 0 + v6 = 27; v6 = 27

u2 + v6 = 21; 27 + u2 = 21; u2 = -6

u4 + v6 = 30; 27 + u4 = 30; u4 = 3

u4 + v4 = 8; 3 + v4 = 8; v4 = 5

u4 + v5 = 27; 3 + v5 = 27; v5 = 24

	v1=9	v2=9	v3=24	v4=5	v5=24	v6=27
u1=0	32	9[92]	24[52]	28	30	27[4]
u2=-6	40	24	22	23	20	21[76]
u3=6	15[112]	25	30[32]	21	39	38
u4=3	39	11	36	8[36]	27[80]	30[16]

Опорный план не является оптимальным, так как существуют оценки свободных клеток, для которых ui + vj > cij

(4;2): 3 + 9 > 11; ?42 = 3 + 9 - 11 = 1

Выбираем максимальную оценку свободной клетки (4;2): 11

Для этого в перспективную клетку (4;2) поставим знак «+», а в остальных вершинах многоугольника чередующиеся знаки «-», «+», «-».

	1	2	3	4	5	6	Запасы
1	32	9[92][-]	24[52]	28	30	27[4][+]	148
2	40	24	22	23	20	21[76]	76
3	15[112]	25	30[32]	21	39	38	144
4	39	11[+]	36	8[36]	27[80]	30[16][-]	132
Потребности	112	92	84	36	80	96

Цикл приведен в таблице (4,2 > 4,6 > 1,6 > 1,2). Из грузов хij стоящих в минусовых клетках, выбираем наименьшее, т.е. у = min (4, 6) = 16. Прибавляем 16 к объемам грузов, стоящих в плюсовых клетках и вычитаем 16 из Хij, стоящих в минусовых клетках. В результате получим план.

	1	2	3	4	5	6	Запасы
1	32	9[76]	24[52]	28	30	27[20]	148
2	40	24	22	23	20	21[76]	76
3	15[112]	25	30[32]	21	39	38	144
4	39	11[16]	36	8[36]	27[80]	30	132
Потребности	112	92	84	36	80	96

Проверим оптимальность опорного плана. Найдем предварительные потенциалы ui, vj. по занятым клеткам таблицы, в которых ui + vj = cij, полагая, что u1 = 0.

u1 + v2 = 9; 0 + v2 = 9; v2 = 9

u4 + v2 = 11; 9 + u4 = 11; u4 = 2

u4 + v4 = 8; 2 + v4 = 8; v4 = 6

u4 + v5 = 27; 2 + v5 = 27; v5 = 25

u1 + v3 = 24; 0 + v3 = 24; v3 = 24

u3 + v3 = 30; 24 + u3 = 30; u3 = 6

u3 + v1 = 15; 6 + v1 = 15; v1 = 9

u1 + v6 = 27; 0 + v6 = 27; v6 = 27

u2 + v6 = 21; 27 + u2 = 21; u2 = -6

	v1=9	v2=9	v3=24	v4=6	v5=25	v6=27
u1=0	32	9[76]	24[52]	28	30	27[20]
u2=-6	40	24	22	23	20	21[76]
u3=6	15[112]	25	30[32]	21	39	38
u4=2	39	11[16]	36	8[36]	27[80]	30

Опорный план является оптимальным, так все оценки свободных клеток удовлетворяют условию

ui + vj <= cij.

Минимальные затраты составят:

F(x) = 9*76 + 24*52 + 27*20 + 21*76 + 15*112 + 30*32 + 11*16 + 8*36 + 27*80 = 9332

Анализ оптимального плана.

Из 1-го склада необходимо груз направить в 2-й магазин (76), в 3-й магазин (52), в 6-й магазин (20)

Из 2-го склада необходимо весь груз направить в 6-й магазин

Из 3-го склада необходимо груз направить в 1-й магазин (112), в 3-й магазин (32)

Из 4-го склада необходимо груз направить в 2-й магазин (16), в 4-й магазин (36), в 5-й магазин (80).

Размещено на Allbest.ru

...