1055
.pdf2)области допустимых решений обеих задач пусты;
3)одна задача имеет неограниченную область допустимых решений, вторая – пустую.
Пример. Для данной задачи составить двойственную, решить ее графическим методом и, используя теоремы двойственности, найти решение данной задачи.
f x1 7x2 |
8x3 |
x4 4x5 |
max; |
|||||
x x |
|
x x |
|
|
x 1; |
|
||
1 |
|
2 |
3 |
|
4 |
|
5 |
|
x1 x2 2x3 x4 2x5 4; |
||||||||
|
|
|
|
|
|
|
|
|
0; |
j 1,5. |
|
|
|||||
xj |
|
|
Решение.
1.Составим матрицу из коэффициентов при неизвестных, свободных членах и коэффициентах целевой функции.
2.Транспонируем полученную матрицу (т. е. заменяем строки на столбцы).
3.По транспонированной матрице составляем двойственную систему ограничений и целевую функцию φ.
|
|
|
|
|
|
|
|
|
|
|
1 |
1 |
1 |
||
|
|
|
|
|
|
|
|
|
|
|
1 |
1 |
7 |
|
|
1 |
1 1 |
1 |
1 |
|
1 |
|
|
|
|||||||
|
|
||||||||||||||
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
1 |
2 1 |
2 |
|
|
|
|
1 |
2 |
8 . |
||||||
1 |
|
4 |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
1 |
1 |
1 |
|||
|
|
|
|
|
|
|
|
||||||||
1 |
7 |
8 |
1 |
4 |
|
f |
|
|
1 |
2 |
4 |
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
1 |
4 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Получаем двойственную задачу:
y1 4y2 |
|
min; |
||||||
y1 |
y2 |
1; |
|
(9) |
||||
y |
y |
2 |
7; |
|
(10) |
|||
1 |
|
|
|
|
|
|
|
|
|
|
2y2 |
8; |
(11) |
||||
y1 |
||||||||
y |
y |
2 |
1; |
|
(12) |
|||
|
1 |
|
|
|
4. |
(13) |
||
y |
2y |
2 |
||||||
|
1 |
|
|
|
|
|
|
Функция φ минимизируется, так как целевая функция исходной задачи максимизируется.
Поскольку на переменные исходной задачи наложены условия неотрицательности xj 0 ( j 1,5), то соотношения (9) – (13) в сис-
теме ограничений двойственной задачи являются неравенствами. Переменные y1, y2 не должны удовлетворять условию неотрицательности, т.к. они соответствуют ограничениям-неравенствам исходной задачи.
Решим полученную задачу графическим методом. На рис. 2 изображены область допустимых решений задачи и grad i 4j . Оптимальное решение задачи Y* = ( 2, 5) и φ(Y*) = 22.
y2
9
11
1 |
y1 |
8
7
10
Рис. 2. Область ограничений двойственной задачи
По первой теореме двойственности, |
f (X*) (Y*) 22. Под- |
|||
ставим оптимальное решение Y* |
( 2, 5) в систему ограничений. По- |
|||
лучим ограничения |
|
|
|
|
2 5 3, |
3 1, |
|
x 0; |
|
|
|
|
|
1 |
2 5 7; |
|
|
|
|
|
|
|
|
|
2 10 8; |
|
|
|
|
|
7 1 |
|
|
0; |
2 5 7, |
x4 |
|||
2 10 12, |
12 4 |
|
x 0. |
|
|
|
|
|
5 |
(9), (12), (13) выполняются как строгие неравенства. Согласно второй теореме двойственности соответствующие компоненты оптимального плана двойственной задачи, т.е. исходной задачи, равны ну-
лю: x1* x*2 x3* 0. Учитывая это, из системы ограничений исходной
задачи найдем ее оптимальное решение:
уравнения: 3x3* 3; x3* 1.
x1* 0;x2* 2;
*
Получаем решение x3 1;
x4* 0;x5* 0.
x* |
x* |
1; |
||
|
2 |
3 |
|
Вычитаем |
|
* |
|
* |
|
|
|
4. |
||
x2 |
2x3 |
X* = (0; 2; 1; 0; 0).
Ответ: fmax 22 при X* (0; 2; 1; 0; 0).
Вопросы для самопроверки
1.Какие задачи линейного программирования называются двойственными?
2.Укажите правила построения двойственной задачи.
3.Сформулируйте основные теоремы двойственности.
4.Приведите геометрическую интерпретацию двойственных за-
дач.
5.Укажите связь решений двойственных задач.
6.Приведите пример экономической задачи, решаемой методами линейного программирования, и экономическую интерпретацию ее двойственной задачи.
§ 2.5. Задачи для самостоятельного решения
Для задач 2.1─2.6 составить двойственную.
Задача 2.1.
f x1 2x2 max ;
5x1 |
2x2 |
3; |
|||||
x |
|
x |
2 |
1; |
|||
1 |
|
|
|
|
|
||
|
|
|
|
|
|
|
|
3x1 x2 3; |
|||||||
3x |
3x |
2 |
9; |
||||
|
|
1 |
|
|
|
j 1,2. |
|
x |
j |
0; |
|
|
|||
|
|
|
|
|
|
|
Задача 2.2.
f x1 2x2 x3 x4 x5 min;
x1 2x2 x3 3x4 2x5 8;
2x1 3x2 2x3 x4 x5 4;x1 3x2 4x5 8;
x1 0; x3 0; x5 0.
Задача 2.3. |
|
|
|
Задача 2.4. |
|
|
|
|
|||
f x1 2x2 3x3 x4 |
max; |
f |
2x1 |
x2 x3 3x4 |
x5 min; |
||||||
2x1 x2 2x3 3x4 5; |
2x1 x2 |
x3 3x4 |
x5 |
10; |
|||||||
|
2x2 |
x3 2x4 x5 8; |
|||||||||
|
2x2 x3 x4 3; |
|
x1 |
||||||||
x1 |
|
|
|
|
|
|
|
|
|||
x |
0, ...,x |
4 |
0. |
|
2x1 x2 |
3x3 x4 |
2x5 4; |
||||
1 |
|
|
|
x |
0; |
x 0; |
x |
4 |
0. |
|
|
|
|
|
|
|
1 |
|
3 |
|
|
|
Задача 2.5.
f 2x1 x2 5x4 min;
x1 3x2 x3 2x4 5;
2x1 x2 2x3 x4 2;x1 3x2 3x4 8;
x1 0, ... ,x4 0.
Задача 2.6.
f 3x1 x2 2x3 |
max; |
||||||
2x1 2x2 x3 2; |
|
||||||
|
|
|
|
2x3 |
6; |
||
3x1 3x2 |
|||||||
|
|
|
|
2x3 |
10; |
||
3x1 3x2 |
|||||||
x |
x |
2 |
2; |
|
|
||
1 |
|
x |
|
0; |
x |
0. |
|
x |
0; |
2 |
|||||
1 |
|
|
|
|
3 |
|
Для задач 2.7─2.12 составить двойственную, решить ее графическим методом и в случае разрешимости найти экстремальное значение целевой функции.
Задача 2.7. |
|
|
|
Задача 2.8. |
|
|||||
f 6x1 9x2 |
3x3 |
min; |
f |
2x1 2x2 |
x3 3x4 |
min; |
||||
x1 2x2 x3 2; |
2x1 2x2 x3 x4 1; |
|
||||||||
|
x2 |
x3 1; |
|
|
||||||
3x1 |
|
|
|
|
|
|
||||
|
|
|
|
|
|
x1 x2 3x3 x4 1; |
|
|||
3x1 |
x2 |
x3 1; |
|
x 0, ...,x |
4 |
0. |
|
|||
x 0; |
x |
2 |
0; x 0. |
|
1 |
|
|
|||
1 |
|
|
|
3 |
|
|
|
|
|
Задача 2.9. |
|
Задача 2.10. |
|
|
|||||||||
f 2x1 |
x2 3x3 |
min; |
f 2x1 4x2 |
23x3 4x4 |
min; |
||||||||
x x |
|
|
x 2; |
x 3x |
|
x |
|
|
2; |
|
|||
|
1 |
|
2 |
3 |
|
1 |
2 |
|
4 |
|
|
||
x1 |
3x2 2x3 1; |
x1 x2 2x3 3x4 2; |
|||||||||||
x |
0; x |
2 |
0; x |
0. |
x 0, |
...,x |
4 |
0. |
|
||||
1 |
|
|
|
3 |
|
1 |
|
|
|
|
|
Задача 2.11. |
|
|
Задача 2.12. |
|
|
|
|||
f 3x1 12x2 4x3 |
min; |
f x1 x2 |
min; |
||||||
x1 3x2 x3 2; |
|
x1 x2 x3 x4 8; |
|||||||
|
4x2 |
4x3 |
1; |
|
|
|
|
x4 0; |
|
x1 |
|
3x1 3x2 x3 |
|||||||
x |
0; x |
2 |
0; |
x 0. |
x 0, ...,x |
4 |
0. |
||
1 |
|
|
3 |
|
1 |
|
|
Для задач 2.13─2.20 составить двойственную, привести графическую интерпретацию решений пары двойственных задач.
Задача 2.13. Задача 2.14.
f x1 x2 |
max; |
x1 3x2 6;3x1 x2 6;x1 0; x2 0.
Задача 2.15.
f 2x1 x2 |
min; |
3x1 x2 1;x1 3x2 5;
x1 0; x2 0.
Задача 2.17.
f x1 2x2 |
max; |
x1 2x2 6;
3x1 5x2 15;
x2 0.
Задача 2.19.
f x2 min;
f 2x1 x2 |
|
|
max; |
|||
2x1 x2 |
2; |
|
||||
|
x2 |
2; |
|
|
|
|
x1 |
|
|
|
|||
x |
0; |
x |
2 |
0. |
|
|
1 |
|
|
|
|
|
|
Задача 2.16. |
|
|
||||
f 8x1 |
4x2 |
|
max; |
4x1 x2 9;2x1 x2 5;
x1 0; x2 0.
Задача 2.18.
f x1 |
2x2 |
|
max; |
||
3x |
2x |
|
6; |
|
|
|
1 |
|
2 |
|
|
x1 |
3x2 3; |
|
|||
|
0. |
|
|
|
|
x2 |
|
|
|
||
Задача 2.20. |
|
||||
f x1 3x2 |
|
|
min; |
2x x |
|
4; |
2x1 2x2 3; |
||||
|
1 |
|
2 |
|
|||
x1 |
x2 |
0; |
|
|
|
||
x 0. |
|
|
|
x1 |
x2 |
1. |
|
1 |
|
|
|
|
|
|
|
Для каждой из задач ЗЛП 2.21─2.30 выписать двойственную, решить ее графическим методом, используя вторую теорему двойственности, перейти от оптимального решения двойственной задачи к оптимальному решению исходной.
Задача 2.21.
f x1 |
6x2 |
2x3 |
|
min; |
|||
x1 4x2 x3 2; |
|
|
|||||
|
x2 |
2x3 |
1; |
|
|
||
x1 |
|
|
|||||
x |
0; x |
2 |
0; x 0. |
|
|||
1 |
|
|
|
3 |
|
|
|
Задача 2.23. |
|
|
|
||||
f x1 |
2x2 |
9x3 |
|
min; |
x1 x2 2x3 2;
4x1 2x2 3x3 3;x1 0; x2 0; x3 0.
Задача 2.25.
f 4x1 3x2 4x3 min;
x1 x2 x3 2;
4x1 x2 2x3 5;
x1 0; x2 0; x3 0.
Задача 2.27.
f x1 2x2 9x3 min;
x1 x2 2x3 1;
2x1 2x2 3x3 1,5;x1 0; x2 0; x3 0.
Задача 2.29
f 6x1 9x2 3x3 min;
Задача 2.22.
f x1 4x2 x3 |
|
min; |
|||||
x1 2x2 x3 2; |
|
||||||
|
x2 2x3 1; |
|
|
||||
x1 |
|
|
|||||
x |
0; x |
2 |
0; x |
0. |
|
||
1 |
|
3 |
|
|
|||
Задача 2.24. |
|
|
|||||
f x1 4x2 9x3 |
|
min; |
|||||
x1 x2 2x3 2; |
|
|
|||||
|
|
|
|
|
3x3 3; |
|
|
2x1 2x2 |
|
||||||
x |
0; x |
2 |
0; x |
0. |
|
||
1 |
|
3 |
|
|
|||
Задача 2.26. |
|
|
|||||
f 3x1 4x2 6x3 |
min; |
||||||
3x1 4x2 3x3 3; |
|
||||||
|
|
|
|
|
|
|
|
x1 x2 2x3 1; |
|
||||||
x |
0; x |
2 |
0; x |
0. |
|
||
1 |
|
|
3 |
|
|
||
Задача 2.28. |
|
|
|||||
f 3x1 2x2 6x3 |
|
min; |
|||||
x1 x2 3x3 1; |
|
|
|||||
|
|
|
|
|
|
|
|
3x1 2x2 2x3 4; |
|
||||||
x |
0; x |
2 |
0; x |
0. |
|
||
1 |
|
|
|
3 |
|
|
|
Задача 2.30. |
|
|
|||||
f 2x1 2x2 |
x3 3x4 min; |
x1 2x2 x3 2; |
2x1 2x2 x3 x4 1; |
||||
|
x3 1; |
|
|
1; |
|
3x1 x2 |
|
x1 x2 3x3 x4 |
|||
x 0; x |
2 |
0; x |
0. |
x 0; i 1,...,4. |
|
1 |
3 |
|
i |
|
Раздел 3. ТРАНСПОРТНАЯ ЗАДАЧА
§ 3.1. Математическая модель транспортной задачи
Пусть имеются m пунктов отправления A1, A2, …, Am, в которых находится однородный груз в количествах а1, а2, …, аm cоответственно, и n пунктов назначения B1, B2, …, Bn , потребности которых в данном грузе равны b1, b2, …, bn. Известны cij расходы на перевозку единицы груза из i-го пункта отправления в j-й пункт потребления. Требуется составить план перевозок так, чтобы запасы каждого поставщика были бы вывезены, спрос каждого потребителя удовлетворен и общая стоимость всех перевозок была минимальной.
Исходные данные транспортной задачи запишем в виде матрицы перевозок (табл. 1).
|
|
|
|
|
|
|
|
|
|
|
Таблица 1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Bj |
B1 |
|
B2 |
|
|
… |
Bn |
|
Запасы |
|||
Ai |
|
|
|
|
|
|
|
|
|
|
|
|
A1 |
С11 |
|
C12 |
|
|
… |
C1n |
|
a1 |
|||
A2 |
С21 |
|
C22 |
|
|
|
|
|
C2n |
|
a2 |
|
… |
… |
|
… |
|
|
… |
… |
|
… |
|||
An |
Сm1 |
|
Cm2 |
|
|
|
|
|
Cmn |
|
am |
|
Потребности |
b |
|
b |
|
|
|
… |
bn |
|
- |
||
|
1 |
|
2 |
|
|
|
|
|
|
|
|
|
Обозначим |
через xij 0 |
(i |
|
; |
j |
|
) |
количество |
единиц |
|||
1, m |
1, n |
груза, которое нужно перевезти из пункта Ai в пункт Bj .
Так как нужно перевезти весь груз из каждого пункта отправления Ai , то должны выполняться равенства
x11 x12 ... |
x1n a1; |
|
|
x22 |
x2n a2; |
x21 |
||
|
|
|
................................... |
||
|
|
|
|
xm2 |
xmn am. |
xm1 |
В каждый пункт назначения Bj должен быть завезен весь требуемый груз, потому
x11 x21 ... |
xm1 b1; |
|||||
|
x22 |
xm2 b2; |
||||
x12 |
||||||
|
|
|
|
|
|
|
................................... |
||||||
|
|
|
|
|
|
|
x |
x |
2n |
... |
x |
mn |
b . |
1n |
|
|
|
n |
Стоимость всех запланированных перевозок должна быть минимальной:
f c11x11 c12x12 ... cmnxmn |
min. |
Математическая модель транспортной задачи (ТЗ) в общем случае имеет вид
m |
|
n |
|
|
|
|
|
|
|
|
|
|
f |
cijxij |
|
min, |
(14) |
||||||||
i 1 |
j 1 |
|
|
|
|
|
|
|
|
|
|
|
n |
|
|
|
|
|
|
|
|
|
|
||
x |
|
a , |
i 1,m |
|
||||||||
ij |
|
i |
|
|
|
|
|
|
|
|||
j 1 |
|
|
|
|
|
|
|
|
|
(15) |
||
m |
|
|
|
|
|
|
|
|
; |
|||
|
|
|
|
|
|
|
|
|
|
|||
x |
b |
, |
j 1,n |
|
||||||||
ij |
|
j |
|
|
|
|
|
|
|
|
||
i 1 |
|
|
|
|
|
|
|
|
|
|
||
xij 0, |
i |
|
; |
j |
|
. |
(16) |
|||||
1,m |
1,n |
Таким образом, математически ТЗ формируется по следующей схеме. Заданы система ограничений (15) при условии (16) и целевая функция (14); требуется среди множества решений системы (15) найти такое неотрицательное решение, которое минимизирует функцию
(14).
В рассмотренной модели ТЗ предполагается, что суммарные запасы поставщиков равны суммарным запросам потребителей, т. е.
m |
n |
|
ai bj . |
(17) |
|
i 1 |
j 1 |
|
Такая задача называется задачей с правильным балансом, ее модель – закрытой. Для того, чтобы ТЗ линейного программирования имела решение, необходимо и достаточно выполнение равенства (17).
Всякое неотрицательное решение системы линейных уравнений
(13), определяемое матрицей |
X (xij ); |
i 1,m; |
j 1,n, называется |
||||
планом ТЗ. План X* (x* ); |
i |
|
; |
j |
|
, при котором целевая |
|
1,m |
1,n |
||||||
ij |
|
|
|
|
|
|
|
функция (14) принимает свое минимальное значение, называется оп-
тимальным планом ТЗ. Матрица (сij ); i 1,m; j 1,n называется матрицей тарифов (издержек или транспортных расходов).
§3.2. Свойства транспортной задачи
1.Ранг матрицы из коэффициентов при неизвестных системы ограничений ТЗ равен m + n – 1, где m и n – количество поставщиков и потребителей соответственно.
2.ТЗ всегда имеет оптимальный план.
3.В ТЗ всегда существуют допустимые планы, содержащие не более m + n – 1 положительных элементов.
4.Если в ТЗ все числа ai , bj целые, то она имеет оптимальный целочисленный план.
Решение (план перевозок) назовем допустимым, если оно удовлетворяет ограничениям (15), (16); опорным, если в нем отличны от нуля не более m + n – 1 базисных переменных, остальные равны нулю.
Решение ТЗ разобьем на три этапа:
1) определение первоначального допустимого решения;
2) проверка найденного решения на оптимальность (оценка плана по критерию стоимости). Если оно оптимальное, то ТЗ решена;
3) улучшение начального плана, т.е. последовательный переход от одного плана к другому, связанный с уменьшением суммарной стоимости перевозок.
§ 3.3. Методы решения транспортной задачи
Классическая транспортная задача решается симплекс-методом. Но для задач небольшой размерности часто проще и быстрее получить решение задачи иными способами. Укажем основные методы решения транспортной задачи.
1.Итерационное улучшение плана перевозок: требуется построить опорный план и последовательными итерациями получить оптимальное решение. Опорный план находят методами «северозападного угла», «наименьшего элемента», «двойного предпочтения», «аппроксимаций Фогеля». После нахождения опорного плана нужно применить один из алгоритмов его улучшения: «метод падающего камня», «метод потенциалов».
2.Решение с помощью теории графов: рассматривается двудольный граф, у которого пункты производства находятся в верхней доле, а пункты потребления – в нижней. Пункты производства и потребления попарно соединяются ребрами бесконечной пропускной способности с указанием стоимости перевозки единицы потока. К верхней доле искусственно присоединяется исток. Пропускная возможность ребер от истока до пунктов производства равна запасу продукта в этих пунктах. Цена за единицу потока этих ребер равна нулю. Аналогично к нижней доле присоединяется сток. Пропускная способность от потребителей к стоку равна потребности в продукте в этих пунктах. Цена за единицу потока равна нулю. Далее решается задача нахождения максимального потока минимальной стоимости. При решении этой задачи обычно используется алгоритм Беллмана-Форда.
§ 3.4. Методы нахождения начального плана перевозок
Клетки матрицы перевозок, где xij 0, называются базисными, а
остальные, где xij 0, – свободными.
В матрице есть m + n – 1 базисных клеток. Их число совпадает с числом независимых уравнений – ограничений.
Значение xij в матрице перевозок находится по формуле