Исследование операций. Курс лекций. Учебное пособие
.pdf
личества n + m −1. Ноль можно поставить только в пустую клетку, которая не образует цикла вместе с уже заполненными.
Циклом называется такая последовательность клеток таблицы ТЗ (i1 , j1 ), (i1 , j2 ); ... (ik , j1 ), в которой две и только две соседние клетки располо-
жены в одной строке или столбце, причем первая и последние клетки также находятся в одной строке или столбце.
Цикл изображают в виде замкнутой ломаной линии. В любой клетке цикла происходит поворот звена ломаной линии на 900. Простейшие циклы изображены на рис. 6.1.
Рис. 6.1. Виды циклов
Построение циклов начинают с какой-либо занятой клетки и переходят по столбцу (строке) к другой занятой клетке, в которой делают поворот на 900 и движутся по строке (столбцу) к следующей занятой клетке и т.д., пытаясь возвратиться к первоначальной клетке.
Проверим вырожденность плана нашей задачи: мы имеем трех поставщиков ( т = 3) и пять потребителей ( п = 5), следовательно n + m −1 = 7. Заполненных клеток 6, т. е. план – вырожденный. Поставим нулевую поставку в клетку (3.3).
Вычислим значение целевой функции (стоимость перевозки) для этого
плана: Z = 50 9 + 40 1+ 20 4 +50 6 + 40 5 + 70 3 =1280 .
6.3. Проверка плана на оптимальность
При проверке плана на оптимальность воспользуемся методом потен-
циалов.
Каждому поставщику Ai и потребителю Bj поставим в соответствие некоторые числа Ui и Vj , называемые потенциалами таким образом, чтобы для
каждой занятой клетки сумма потенциалов равнялась бы стоимости перевозки, стоящей в этой клетке, т. е. Ui +V j = Cij .
Количество неизвестных потенциалов составит n + m , а уравнений в системе ограничений всего n +m −1. Поэтому один из потенциалов задается про-
извольно, после чего остальные потенциалы определяются однозначно. Обычно полагают U1 = 0 .
51
Критерий оптимальности: для того чтобы план ТЗ был оптимальным, необходимо и достаточно, чтобы для всех пустых (незанятых) клеток сумма потенциалов должна быть меньше или равна стоимости перевозки, стоящей в этой клетке, т. е. Ui +V j ≤ Cij .
Проверим на оптимальность построенный нами план. Для этого вначале построим систему потенциалов. Полагаем U1 = 0 , а далее Ui + Vj = сij для заня-
тых клеток таблицы.
(1.1): U1 |
+ V1 |
= 9 ; |
V1 |
= 9; |
(1.2): U1 |
+ V2 |
= 1; |
V2 |
= 1; |
(2.2): U2 |
+ V2 |
= 4; |
U2 |
= 3; |
(2.3): U2 |
+ V3 |
= 6; |
V3 |
= 3; |
(3.3): |
V3 +U3 =3; |
U3 = 0; |
||||
(3.4): |
U3 |
+ V4 |
= |
5 ; |
V4 |
= 5; |
(3.5): |
U3 |
+ V5 |
= 3 ; |
V5 |
= 3. |
|
Таблица 6.4
Пункт |
|
|
Пункты потребления |
|
Запасы |
||
поставки |
Ui Vj |
В1 |
В2 |
В3 |
В4 |
В5 |
|
|
V1 = 9 |
V2 = 1 |
V3 = 3 |
V4 = 5 |
V5 = 3 |
|
|
А1 |
U1 = 0 |
9 |
1 |
1 |
5 |
6 |
90 |
|
|
50 |
40 |
|
|
|
|
А2 |
U2 = 3 |
6 |
4 |
6 |
8 |
5 |
70 |
|
|
|
20 |
50 |
|
|
|
А3 |
U3 = 0 |
2 |
9 |
3 |
5 |
3 |
110 |
|
|
|
|
0 |
40 |
70 |
|
Потребности |
|
50 |
60 |
50 |
40 |
70 |
270 |
Проверим критерий оптимальности. Просматриваем строки и для каждой незанятой клетки проверяем выполнение условия Ui +V j ≤ Cij :
(1.3): U1 + V3 = 3 > 1 на 2
(1.4): U1 + V4 = 5 = 5
(1.5): U1 + V5 = 3 < 6
(2.1): U2 + V1 = 12 > 6 на 6
(2.4): U2 + V4 = 8 = 8 (2.5): U2 + V5 = 6 > 5 на 1
(3.1): U3 + V1 = 9 > 2 на 7
(3.2): U3 + V2 = 1 < 9
Поскольку условие оптимальности не выполнено для четырех клеток, то построенный план не является оптимальным и может быть улучшен.
6.4. Улучшение плана
Из условий, где критерий не выполняется, выбираем то условие, где разница максимальна. Это – клетка (3.1), и для нее строим цикл пересчета.
Перебросим в клетку (3.1) 50 единиц груза из клетки (1.1).
Чтобы компенсировать недостаток в первой строке, перебросим те же 50 единиц груза из клетки (2.3) в клетку (1.3).
52
Далее, чтобы компенсировать недостаток во второй строке, перебросим |
||||||||
из клетки (3,5) 50 единиц в клетку (2.5). |
|
|
|
|
||||
Таким образом, образовался цикл, показанный в табл. 6.5 пунктиром. |
||||||||
|
|
|
|
|
|
|
Таблица 6.5 |
|
Пункт |
|
|
|
Пункты потребления |
|
Запасы |
||
поставки |
|
|
В1 |
В2 |
В3 |
В4 |
В5 |
|
|
Ui |
Vj |
V1 = 9 V2 = 1 V3 = 3 V4 = 5 V5 = 3 |
|
||||
А1 |
U1 = 0 |
9 |
1 |
1 |
5 |
6 |
90 |
|
|
|
|
-50 50 |
40 |
+50 |
|
|
|
А2 |
U2 = 3 |
6 |
4 |
6 |
8 |
5 |
70 |
|
|
|
|
|
20 |
-50 50 |
|
+50 |
|
А3 |
U3 = 0 |
2 |
9 |
3 |
5 |
3 |
110 |
|
|
|
|
+50 |
|
0 |
40 |
-50 70 |
|
Потребности |
|
50 |
60 |
50 |
40 |
70 |
270 |
|
Клетка (3.1) – вершина цикла. Она находится в пустой клетке, остальные вершины – в заполненных. Клетка (3.1) помечена знаком «+», далее последовательно расставлены знаки «-» и «+». Среди клеток, помеченных знаком «-» обычно определяют наименьшее значение перевозки, и его перемещают по циклу. В нашем случае эта величина составила 50 единиц груза.
Получаем новую таблицу, для которой повторяем расчет потенциалов: Полагаем U1 = 0, а далее Ui + Vj = сij для занятых клеток таблицы.
U1 |
+ V2 |
= 1; V2 |
= 1 |
|
U3 |
+ V5 |
= 3; U3 |
= 1 |
|
||
U1 |
+ V3 |
= 1; V3 |
= 1 |
|
U3 |
+ V1 |
= 2; V1 |
= 5 |
|
||
U2 |
+ V2 |
= 4; U2 |
= 3 |
|
U3 |
+ V4 |
= 5 ; V4 |
= 4 |
|
||
U2 |
+ V5 |
= 5; V5 |
= 2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Таблица 6.6 |
|
Пункт |
|
|
Пункты потребления |
|
Запасы |
||||||
|
В1 |
В2 |
В3 |
|
|
В4 |
|
В5 |
|||
поставки |
Ui Vj |
V1 = 1 |
V2 = 1 |
V3 = 1 |
|
V4 = 4 |
V5 = 2 |
|
|||
А1 |
|
U1 = 0 |
9 |
1 |
1 |
|
|
5 |
|
6 |
90 |
|
|
|
|
40 |
50 |
|
|
|
|
|
|
А2 |
|
U2 = 3 |
6 |
4 |
6 |
|
|
8 |
|
5 |
70 |
|
|
|
|
20 |
|
|
|
|
|
50 |
|
А3 |
|
U3 = 1 |
2 |
9 |
3 |
|
|
5 |
|
3 |
110 |
|
|
|
50 |
|
|
|
40 |
|
20 |
|
|
Потребности |
|
50 |
60 |
50 |
|
40 |
|
70 |
270 |
||
Проверим критерий оптимальности: Ui + Vj ≤ сij для свободных клеток. |
|||||
U1 |
+ V1 |
= 1 < 9 |
U1 |
+ V5 |
= 2 < 6 |
U1 |
+ V4 |
= 4 < 5 |
U2 |
+ V1 |
= 4 < 6 |
53
U2 |
+ V3 |
=4 < 6 |
|
U3 |
+ V2 |
=2 < 9 |
|
|
||
U2 |
+ V4 |
=7 < 8 |
|
U3 |
+ V3 |
=2 < 3 |
|
|
||
Критерий выполнен, значит, полученное решение оптимально. |
|
|||||||||
|
|
|
0 |
40 |
50 |
0 |
0 |
|
||
|
|
|
|
|
|
|
|
|
|
|
Итак, оптимальный план перевозок: X * = |
0 |
20 |
0 |
0 |
50 |
|
||||
|
|
|
50 |
0 |
0 |
40 |
20 |
|
||
|
|
|
|
|||||||
Найдем минимальную стоимость перевозок:
Z = 40 +50 + 4 20 +5 50 + 2 50 +5 40 +3 20 = 780 .
6.5. Открытая модель транспортной задачи
Открытые модели ТЗ могут быть двух видов:
|
|
m |
n |
|
|
а) ∑ai |
>∑bj |
– запас превышает потребность; |
|
|
|
i=1 |
j=1 |
|
|
|
m |
n |
|
|
б) ∑ai <∑bj |
– потребность превышает запас. |
||
|
|
i=1 |
j=1 |
|
|
В случае а) вводят фиктивный (n +1)-й пункт назначения с потребно- |
|||
|
|
m |
n |
|
стью |
bn+1 |
= ∑ai −∑bj и полагают стоимость перевозок в этот пункт назначения |
||
|
|
i=1 |
j=1 |
|
ci,n+1 = 0 . |
Соответствующие значения в оптимальном плане будут означать |
|||
оставшийся неотправленным груз поставщиков. |
||||
|
В случае б) вводится фиктивный (m +1)-й пункт отправления с запасами |
|||
|
|
n |
m |
|
груза |
am+1= ∑b j −∑ai и стоимостью перевозок cm+1. j = 0 . Груз, соответствующий |
|||
|
|
j=1 |
i=1 |
|
поставкам от фиктивного поставщика, останется недополученным потребителями.
Таким образом, открытая модель ТЗ может быть сведена к закрытой мо-
дели.
Пример. Исходные данные транспортной задачи приведены в табл. 6.7. Таблица 6.7
Пункт |
Запасы |
|
Пункты потребления |
|
||||
В1 |
|
В2 |
В3 |
|
В4 |
|
||
поставки |
|
|
|
|
||||
А1 |
120 |
5 |
|
6 |
2 |
|
5 |
|
А2 |
100 |
6 |
|
9 |
3 |
|
4 |
|
Потребности |
|
40 |
|
180 |
60 |
|
170 |
|
|
m |
n |
|
|
|
|
|
|
Уравнение баланса не выполнено: ∑ai <∑bj . Таким образом, |
это задача |
|||||||
|
i=1 |
j=1 |
|
|
|
|
|
|
открытого типа. Поскольку спрос превышает предложение, добавляем фиктивного поставщика А3 с объем запасов 360 - 220 = 140 единиц груза.
54
Таблица 6.8
Пункт |
Запасы |
|
Пункты потребления |
|
|||
В1 |
|
В2 |
В3 |
|
В4 |
||
поставки |
|
|
|
||||
А1 |
120 |
5 |
|
6 |
2 |
|
5 |
А2 |
100 |
6 |
|
9 |
3 |
|
4 |
А3 |
140 |
0 |
|
0 |
0 |
|
0 |
потребности |
|
40 |
|
180 |
60 |
|
170 |
Получена закрытая транспортная задача: спрос равен предложению. Решение задачи такого типа рассмотрено выше.
ЗАКЛЮЧЕНИЕ
В предлагаемом учебном пособии изложен теоретический материал для изучения основ курса «Исследование операций», представлена методика практического решения оптимизационных задач.
Рассмотренные задачи и их экономико-математические модели имеют прикладное значение для выработки эффективных решений в планировании и управлении производством. Методы решения задач позволяют найти из множества альтернативных вариантов наилучший вариант производства, распределения или потребления. Ограниченные ресурсы при этом будут использованы наиболее эффективным образом для достижения поставленной цели.
Пособие охватывает лишь основные разделы исследования операций, в частности это линейное программирование. Дальнейшее изучение дисциплины касается моделирования нелинейности в экономических процессах, рассмотрения специальных задач исследования операций (целочисленные модели, дроб- но-линейные модели, сети и сетевые модели и др.).
55
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
а) основная литература
1.Исследование операций в экономике / под ред. Н. Ш. Кремера. – М.: ЮНИТИ, 2013 .
2.Грызина, Н. Ю. Математические методы исследования операций в экономике [Электронный ресурс]: учеб. пособие / Н. Ю. Грызина, И. Н. Мастяева, О. Н. Семенихина. – Электрон. текстовые данные. – М.: Евразийский открытый институт, 2009. – 196 c. – Режим доступа: http://www.iprbookshop.ru/10773. – ЭБС «IPRbooks».
3.Исследование операций в экономике [Электронный ресурс]: учеб. пособие / Г. Я. Горбовцов [и др.]. – Электрон. текстовые данные. – М.: Евразийский открытый институт, Московский государственный университет экономики, статистики и информатики, 2006. – 118 c. – Режим доступа: http://www.iprbookshop.ru/10690. – ЭБС «IPRbooks».
4.Брусенцев, А. Г. Исследование операций и теория игр [Электронный ресурс]: учеб. пособие / А. Г. Брусенцев, В. И. Петрашев, Ю. Д. Рязанов. – Электрон. текстовые данные. – Белгород: Белгородский государственный технологический университет им. В. Г. Шухова, ЭБС АСВ, 2012. – 258 c. – Режим до-
ступа: http://www.iprbookshop.ru/49709. – ЭБС «IPRbooks».
5.Математика в экономике. Ч. 1. Линейная алгебра, аналитическая геометрия и линейное программирование [Электронный ресурс]: учебник / А. С. Солодовников [и др.]. – Электрон. текстовые данные. – М.: Финансы и статистика,
2011. – 384 c. – Режим доступа: http://www.iprbookshop.ru/12434. – ЭБС
«IPRbooks».
б) дополнительная литература
1. Лунгу, К. Н. Линейное программирование. Руководство к решению задач [Электронный ресурс] / Лунгу К. Н. – Электрон. текстовые данные. – М.: ФИЗМАТЛИТ, 2009. – 132 c. – Режим доступа: http://www.iprbookshop.ru/12905. – ЭБС «IPRbooks».
2. Сборник задач по курсу «Математика в экономике». Ч. 1. Линейная алгебра, аналитическая геометрия и линейное программирование [Электронный ресурс]: учеб. пособие / С. В. Пчелинцев [и др.]. – Электрон. текстовые данные. – М.: Финансы и статистика, 2013. – 256 c. – Режим доступа: http://www.iprbookshop.ru/18835.— ЭБС «IPRbooks».
56
Учебное издание
Виктория Всеволодовна Потихонова Лариса Ивановна Король
Исследование операций
Учебное пособие
Издательский редактор Т. Н. Козлова
Учебное электронное издание сетевого распространения
Системные требования:
Электронное устройство с программным обеспечением для воспроизведения файлов формата PDF
Режим доступа: http://publish.sutd.ru/tp_get_file.php?id=2017152, по паролю
– Загл. с экрана.
Дата подписания к использованию 30.01.2017 г. Рег. № 152/17
ФГБОУВО «СПбГУПТД» Юридический и почтовый адрес:
191186, Санкт-Петербург, ул. Большая Морская, 18 http://sutd.ru
57
