Математические основы теории систем. Методы оптимизации
.pdf2x1 − x2 + x3 = 4,− + =
x1 2x2 x4 2,+ + =
x1 x2 x5 5,
G= − x1 + x2 → min.
2.Второе ограничение было домножено на (–1). Составляется симплекс-таблица:
|
1 |
2 |
|
|
|
|
3 |
2 |
–1 |
4 |
|
|
|
|
|
|
|
|
|
разрешающая строка |
4 |
1 |
–2 |
2 |
|
|
|
|
|
|||||
5 |
1 |
1 |
5 |
|
|
|
|
|
|
|
|
|
|
|
–1 |
1 |
0 |
|
|
|
|
|
|
|
|
|
|
разрешающий столбец
Поскольку имеются отрицательные коэффициенты целевой функции, то вершина, которой соответствует симплекс-таблица, неоптимальна. Максимальный по модулю отрицательный коэффициент ЦФ – (–1), следовательно, первый столбец является разрешающим.
Для определения разрешающей стройки находятся отношения правых частей ограничений к положительным коэффициентам разрешающего столбца и выбирается минимальное отношение:
4 |
|
|
|
|
|
= 2 |
|
||
2 |
|
|||
|
|
|
||
2 |
|
|
|
|
|
= 2 |
можно выбрать любое отношение (допустим, второе). |
||
1 |
||||
|
|
|
||
5 |
= 5 |
|
|
|
|
|
|
||
1 |
|
|||
|
|
|
3. Составляется новая симплекс-таблица: разрешающий элемент = 1; новый разрешающий элемент: 1:1 = 1;
новая разрешающая строка: 1 (–2) = –2; 1 2 = 2;
новый разрешающий столбец: 2 (–1) = –2; 1 (–1) = –1; –1 (–1) = 1.
31
Элементы других столбцов: |
|
|
|
|
|
|
|||||
|
–1 |
|
2 |
|
3 |
|
4 |
|
2 |
|
0 |
|
|
|
|
|
|
||||||
|
1 |
+2 |
1 |
= |
3 |
|
5 |
–2 |
1 |
= |
3 |
|
1 |
|
–1 |
|
–1 |
|
0 |
|
–1 |
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
Новая симплекс-таблица: |
|
|
|
|
|||
|
|
4 |
2 |
|
|
|
|
|
|
|
|
|
|
||
3 |
–2 |
3 |
0 |
|
|
разрешающая строка |
|
|
|
||||||
1 |
1 |
–2 |
2 |
|
|
|
|
5 |
–1 |
3 |
3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
–1 |
2 |
|
|
|
|
|
|
|
разрешающий столбец |
Проверка ограничений: 2·2 – 0 – 0 = 4, 2 – (–2)·0 + 0 = 2,
2 + 0 + 3 = 5.
Проверка целевой функции:
–2 + 0 = –2.
Как видим, полученная вершина неоптимальна, так как среди коэффициентов целевой функции есть отрицательный коэффициент (–1). Требуется перейти к следующей вершине.
4. Составляется новая симплекс-таблица:
|
4 |
3 |
|
2 |
–2/3 |
1/3 |
0 |
1 |
–1/3 |
2/3 |
2 |
5 |
1 |
–1 |
3 |
|
|
|
|
|
1/3 |
1/3 |
2 |
|
|
|
|
Проверка ограничений: 2·2 – 0 + 0 = 4, 2 – 2·0 + 0 = 2, 2 + 0 + 3 = 5.
Проверка целевой функции: –2 + 0 = –2.
32
Все коэффициенты целевой функции положительны, следовательно, найдено оптимальное и единственное решение задачи.
Ответ: координаты вершины: x1 = 2, x2 = 0, x3 = 0, x4 = 0, x5 = 3,
Gmin = –2.
Задание для самостоятельной работы
Дано: целевая функция и ограничения в виде неравенств. Требуется: решить ЗЛП графическим и симплекс-методами.
1.В системах координат Х1, Х2 на основании заданных ограничений построить область допустимых решений.
2.Найти решение ЗЛП с помощью построения графика целевой функции.
3.Найти решение ЗЛП посредством вычисления целевой функции в граничных точках области допустимых решений. Сравнить результаты двух методов.
4.Привести ЗЛП к каноническому виду. Если ограничение имеет
вид неравенства типа «≥», а в правой части ограничения стоит отрицательное число, то следует левую и правую часть ограничения умножить на –1, при этом неравенство приобретет вид «≤».
5.Построить симплекс-таблицу и решить ЗЛП симплексметодом.
6.Если в правой части некоторой строки симплекс-таблицы стоит «0», то данную строку следует принять в качестве разрешающей.
7.После каждого пересчета коэффициентов симплекс-таблицы следует проверять выполнение ограничений-равенств. Если ограничения не выполняются, то пересчет произведен неверно.
8.Если в столбце свободных членов симплекс-таблицы после пересчета коэффициентов появились отрицательные числа, то либо неверно произведен пересчет коэффициентов, либо неправильно выбран разрешающий элемент.
1. Q = –3x1 – 2x2 → min |
2. Q = x1 – 2x2 → min |
x1 + 2x2 ≤ 7 |
–x1 + x2 ≤ 0 |
2x1 + x2 ≤ 8 |
2x1 + x2 ≤ 3 |
x2 ≤ 3 |
– x1 + x2 ≥ –1 |
x1, x2 ≥ 0 |
x1, x2 ≥ 0 |
33
3. Q = –x1 – 3x2 → min 2x1 + x2 ≤ 2
x1 – x2 ≥ 0 x1 – x2 ≤ 1 x1, x2 ≥ 0
5. Q = x2 – x1 → min 2x1 – x2 ≤ 4
–x1 + 2x2 ≥ –2 x1 + x2 ≤ 5 x1, x2 ≥ 0
7. Q = x2 – x1 → min –2x1 + x2 ≤ 2
x1 – 2x2 ≤ 2 x1 + x2 ≤ 5 x1, x2 ≥ 0
9. Q = 2x1 + x2 → max 2x1 + x2 ≤ 1
3x1 – x2 ≥ –1 x1 – 4x2 ≤ 2 x1, x2 ≥ 0
11. Q = –9x1 – 11x2 → min 4x1 + 3x2 ≤ 10
x1 ≤ 5
x1 + 2x2 ≤ 8 x1, x2 ≥ 0
13. Q = x1 – 10x2 → min 3x1 + x2 ≤ 12
–8x1 + 3x2 ≤ 24 x1, x2 ≥ 0
34
4. Q = –x1 – x2 → min x1 ≤ 3
x2 ≤ 2
x1 +x2 ≤ 1 x1, x2 ≥ 0
6. Q = 2x1 + 3x2 → max x1 + x2 ≤ 5
x1 + 3x2 ≤ 9 x1 ≤ 4
x1 + 2x2 ≤ 8 x1, x2 ≥ 0
8. Q = 4x1 + 6x2 → max x1 + x2 ≤ 18
0,5x1 + x2 ≤ 12 x1 ≤ 12
x2 ≤ 9 x1, x2 ≥ 0
10.Q = x1 + x2 → max x1 – x2 ≥ –1
x1 – x2 ≤ 1 x1 ≤ 2
x2 ≤ 2 x1, x2 ≥ 0
12.Q = –4x1 – 3x2 → min 4x1 + x2 ≤ 10
2x1 + 3x2 ≤ 8 x1, x2 ≥ 0
14.Q = x1 – 20x2 → min
–x1 + 10x2 ≤ 40 4x1 + 2x2 ≤ 29 x1, x2 ≥ 0
15. Q = –x1 – 2x2 → min |
16. Q = 9x1 + 4x2 → mах |
x1 + 2x2 ≤ 7 |
2х1 + 4х2 ≥ –14 |
2x1 + x2 ≤ 8 |
2х1 + х2 ≤ 8 |
x2 ≤ 3 |
2x2 ≤ 6 |
x1, x2 ≥ 0 |
x1, x2 ≥ 0 |
17. Q = – 2x1 – 4x2 → min |
18. Q = 2x1 + 6x2 → mах |
x1 – x2 ≥ – 7 |
–2х1 – х2 ≥ –2 |
4x1 + 2x2 ≤ 3 |
2х1 – 2х2 ≥ 0 |
x1 – x2 ≤ 1 |
x1 – x2 ≤ 1 |
x1, x2 ≥ 0 |
x1, x2 ≥ 0 |
19. Q = x1 + x2 → mах |
20. Q = –x1 + x2 → mах |
2x1 ≥ –6 |
–2х1 – х2 ≥– 4 |
x2 ≤ 2 |
–2х1 + 4х2 ≤ 4 |
x1 + x2 ≤ 1 |
x1 + x2 ≤ 5 |
x1, x2 ≥ 0 |
x1, x2 ≥ 0 |
2.5. Приведение ЗЛП к каноническому виду
Процесс приведения задачи к каноническому виду называется нахождением начальной вершины.
1-й случай. Ограничения имеют вид неравенств типа (≤), ЦФ стремится к минимуму:
a11x1 + a12 x2 + ... + a1n xn |
≤ b1, |
a21x1 + a22 x2 + ... + a2n xn |
≤ b2 , |
|
(2.10) |
... |
|
|
|
am1x1 + am2 x2 + ... + amn xn ≤ bm , |
|
|
|
b1,b2 ,...,bm ≥ 0, |
|
Q = p1x1 + p2 x2 + ... + pn xn |
→ min. |
Задача приводится к каноническому виду путем введения искусственныхпеременных, которыеявляются базисными. ЦФнеменяется.
35
a x + a x |
+ ... + a x |
+ x |
= b , |
|
|
11 1 12 2 |
1n n |
n+1 |
|
1 |
|
a21x1 + a22 x2 + ... + a2n xn + xn+2 |
= b2 |
, |
|||
|
|
|
|
|
(2.11) |
... |
|
|
|
|
|
am1x1 + am2 x2 + ... + amn xn + [x] |
= bm . |
||||
|
|
|
|
n+m |
|
2-й случай. Ограничения имеют вид неравенств (≥), ЦФ минимизируется. Кроме этого, все коэффициенты ЦФ неотрицательны, а среди bi имеется хотя бы один положительный.
Вводятся искусственные переменные, которые вычитаются из правых частей ограничений. Все коэффициенты ограничений домножаются на (–1):
−a11x1 − a12 x2 − ... − a1n xn + xn+1 = −b1, |
|
||||||||
|
|
− a22 x2 |
− ... − a2n xn + xn+ |
2 = −b2 |
, |
||||
−a21x1 |
|||||||||
|
|
|
|
|
|
|
|
|
(2.12) |
... |
|
|
|
|
|
|
|
|
|
−a |
x |
− a |
m |
x |
− ... − a x |
+ x |
+m |
= −b . |
|
|
m1 1 |
|
2 2 |
mn n |
n |
|
m |
Далее задача решается двойственным симплекс-методом. Двойственный симплекс-метод применяется, когда имеются базисные переменные в ограничениях-равенствах, все коэффициенты ЦФ положительны, ЦФ минимизируется, а среди правых частей ограничений имеется хотя бы одно отрицательное значение. Алгоритм двойственного симплекс-метода рассмотрен в п. 2.6.
3-й случай. Ограничения имеют вид равенств, ЦФ минимизируется.
a11x1 + a12 x2 + ... + a1n xn |
= b1, |
|
a21x1 + a22 x2 + ... + a2n xn |
= b2 , |
|
|
|
(2.13) |
... |
|
|
|
+ ... + amn xn = bm , |
|
am1x1 + am2 x2 |
||
|
|
|
xi ≥ 0, |
|
|
Q = p1x1 + p2 x2 + ... + pn xn |
→ min. |
36
К такому виду сводятся все задачи линейного программирования, т.е. этот случай является общим.
Существует несколько методов приведения таких задач к каноническому виду. Один из них – метод искусственного базиса.
2.5.1. Метод искусственного базиса
Суть метода: вводятся искусственные переменные xn+1, xn+2,…, xn+m. Эти искусственные переменные добавляются к левым частям ограничений:
a |
x + a |
x + ... + a |
x + x |
= b , |
|
||||
11 1 |
12 2 |
1n n |
n+1 |
|
1 |
|
|||
a21x1 + a22 x2 + ... + a2n xn + xn+ |
2 = b2 |
, |
|||||||
|
|
|
|
|
|
|
|
|
(2.14) |
... |
|
|
|
|
|
|
|
|
|
a |
x |
+ a |
x |
+ ... + a |
x |
+ x |
+m |
= b . |
|
|
m1 1 |
|
m2 2 |
|
mn n |
n |
|
m |
Вводится искусственная целевая функция:
G = xn+1 + xn+2 + ... + xn+m → min, xi ≥ 0,
i = 1,(n + m).
Искомые переменные исключаются из выражения искомой целевой функции, при этом искомая целевая функция оказывается выраженной через искомые переменные x1,…,xn.
Полученная задача приведена к каноническому виду. В ней искусственные переменные являются базисными. Далее эта задача решается симплекс-методом.
Чтобы ограничения остались прежними, в результате решения должны быть получены нулевые значения искусственных переменных и искусственной целевой функции.
Особенности решения: если в процессе решения искусственная переменная переходит из базисной в свободную, то столбец, соответствующий этой переменной, исключается из рассмотрения, так как эта переменная становится равной нулю.
37
Анализ полученного решения.
1.Оптимальное решение полученной целевой функции положительно.
В этом случае задача не имеет решения.
2.Оптимальное решение полученной целевой функции равно нулю и среди базисных переменных нет ни одной искусственной.
В этом случае исходные переменные оказались поделенными на две группы: свободные и базисные. После этого осуществляется переход к искомой целевой функции. Для этого в исходном выражении целевойфункциибазисные переменные выражаются через свободные.
Пример 1
3x1 |
+ 2x3 − x5 |
= 12, |
x1 − x2 + x3 = 5, |
||
|
|
|
x1 + x3 + x4 = 6, |
||
Q = 2x1 + x2 |
− x3 + 3x4 |
− x5 → min. |
1. Вводим искусственные переменные x6 , x7 , x8 . Искусственная целевая функция: G = x6 + x7 + x8 → 0 ,
G = x6 + x7 + x8 = 12 − 3x1 − 2x3 + x5 − x1 + x2 − x3 + 6 − x1 − x3 − x4 =
=23 − 5x1 + x2 − 4x3 − x4 + x5 → min .
2.Составляем симплекс-таблицу:
|
1 |
2 |
3 |
4 |
5 |
|
6 |
3 |
0 |
2 |
0 |
–1 |
12 |
|
|
|
|
|
|
|
7 |
1 |
–1 |
1 |
0 |
0 |
5 |
8 |
1 |
0 |
1 |
1 |
0 |
6 |
|
|
|
|
|
|
|
|
5 |
1 |
–4 |
–1 |
1 |
–23 |
|
|
|
|
|
|
|
Далее по симплекс-методу находим разрешающий элемент 3. Меняем местами переменные x1 и x6, причем столбец, соответствующий x6, вычеркиваем из таблицы.
38
Новая симплекс-таблица:
|
2 |
3 |
4 |
5 |
|
1 |
0 |
2/3 |
0 |
–1/3 |
4 |
7 |
–1 |
1/3 |
0 |
1/3 |
1 |
|
|
|
|
|
|
8 |
0 |
1/3 |
1 |
1/3 |
2 |
|
|
|
|
|
|
|
1 |
–2/3 |
–1 |
–2/3 |
–3 |
|
|
|
|
|
|
Теперь свободной становится переменная x8, следовательно, вычеркивается соответствующий столбец.
|
2 |
3 |
5 |
|
1 |
0 |
2/3 |
–13 |
4 |
7 |
–1 |
1/3 |
1/3 |
1 |
|
|
|
|
|
4 |
0 |
1/3 |
1/3 |
2 |
|
|
|
|
|
|
1 |
–1/3 |
–1/3 |
–1 |
|
|
|
|
|
Из равноценных столбцов желательно выбрать такой, чтобы искусственная переменная стала свободной.
В результате всех преобразований получена таблица:
|
2 |
3 |
|
1 |
–1 |
1 |
5 |
5 |
–3 |
1 |
3 |
4 |
1 |
0 |
1 |
|
|
|
|
|
0 |
0 |
0 |
|
|
|
|
Как видно из таблицы, в результате тождественного преобразования исходные переменные разделены на базисные и свободные.
3. Переход к исходной целевой функции:
Q = 2x1 + x2 − x3 + 3x4 − x5 .
Запишемограничения, полученныеизитоговойсимплекс-таблицы:
x1 − x2 + x3 = 5,− + =
x5 3x2 x3 3,+ =
x4 x2 1.
39
Выразим базисные переменные через свободные и подставим эти выражения в целевую функцию:
Q = 2x1 + x2 − x3 + 3x4 − x5 =
= 2(5 + x2 − x3 ) + x2 − x3 + 3(1− x2 ) − (3 + 3x2 − x3 ) = 10 − 3x2 − 2x3.
В итоге получена начальная симплекс-таблица, т.е. задача приведена к каноническому виду:
|
2 |
3 |
|
1 |
–1 |
1 |
5 |
5 |
–3 |
1 |
3 |
4 |
1 |
0 |
1 |
|
|
|
|
|
–3 |
–2 |
–10 |
|
|
|
|
Пример 2
x1 + x2 + x3 + x4 = 7,− + + + =
3x1 x2 2x3 x4 6,+ + − =
2x1 x2 x3 x4 2,
Q = − x1 + x2 + 2x3 − x4 → min.
1. Вводятся искусственные переменные x5 , x6 , x7 :
x1 + x2 + x3 + x4 + x5 = 7,− + + + + =
3x1 x2 2x3 x4 x6 6,+ + − + =
2x1 x2 x3 x4 x7 2.
Искусственная целевая функция: G = x5 + x6 + x7 → min.
G = x5 + x6 + x7 = 7 − x1 − x2 − x3 − x4 + 6 + 3x1 − x2 − 2x3 − x4 − − x6 + 2 − 2x1 − x2 − x3 − x4 = 15 − 3x2 − 4x3 − x4 → min.
40