Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

5066

.pdf
Скачиваний:
2
Добавлен:
13.11.2022
Размер:
932.47 Кб
Скачать

21

опорные оптимальные планы X 1 , X 2 , ..., X k и записать оптимальное решение в виде выпуклой линейной комбинации этих планов:

 

 

 

 

 

 

 

k

 

X опт. t1 X 1 t2 X 2 ... tk X k , где t j 0,

t j 1.

 

 

 

 

 

 

 

j

1

Решить симплексным методом:

1.

Z

x1

max

 

 

4x1

3x2

12,

 

 

x1

x2

2,

 

 

x1

x2

2,

 

 

x1,2

0 .

 

3.

Z

x1

2x2

max

x1 x2 5, 2x1 6, x2 5,

x1,2 0.

5.

Z

x2

x3

max

 

 

 

x1

x2

 

x3

 

1,

 

 

 

x2

 

2x3

x4

2,

 

 

x j

0, j

1,...,4.

 

 

7.

Z

2x1

 

3x2

5x3

min

 

 

2x1

x2

x3

 

4,

 

 

x1

2x2

 

x4

12,

 

 

x j

0 , j 1,...,4.

 

 

9.

Z

x1

2x2

2x3

x4

max

 

x1

x3

 

0.5x4

1,

 

 

x2 x3

 

x4

1,

 

 

x j

0,

j

1,...,4.

 

 

2. Z

2x1

x2

3x3 2x4

x5

max

 

x1

x2

x3

 

1,

 

 

x1

x2

 

x4

1,

 

 

x1

x2

 

 

x5

1,

 

 

x j

0 , j

1,...,5.

 

 

 

4. Z

4x1

2x2

max

 

 

 

2x1

x2

14,

 

 

 

 

 

x1

x2

10,

 

 

 

 

 

x1

 

5,

 

 

 

 

 

x1,2

0 .

 

 

 

 

 

6. Z

2x1

3x2

5x3

max

 

x1

x3

1,

 

 

 

 

 

x2

x3

 

6,

 

 

 

 

x j

0,

j

1, 2, 3.

 

 

 

8. Z

x1

2x2

x3

2x4

 

x5

min

x1

2x2

 

x3

 

 

2,

 

2x1

x2

 

x4

 

0,

 

x1

3x2

 

 

x5

 

6,

 

x j

0 , j

1,...,5 .

 

 

 

10. Z

x1

2x2

3x3

x4

2x5

min

 

x1

3x2

 

4x3

 

 

6,

 

 

 

2x2

 

5x3

x4

 

4,

 

 

 

x2

 

2x3

x5

1,

 

 

x j

0,

j

1,...,5.

 

 

 

22

11. Z

x1

x2

max

 

 

12. Z

2x1

2x2

max

 

 

2x1

x2

2,

 

 

 

x1

x2

5,

 

 

 

x1

2x2

2,

 

 

2x1

x2

2,

 

 

 

x1

 

x2

5,

 

 

 

 

x1

2x2

2,

 

 

 

x1,2

0.

 

 

 

 

x1,2

0.

 

 

 

 

 

 

13. Z

2x1

3x2

5x3

max

 

14. Z

2x1

3x2

6x3

3x4

max

x1

 

x2

 

x3

3,

 

 

2x1

x2

 

x3

x4

1,

2x1

 

3x2

x3

5,

 

 

2x1

x2

 

x3

x4

2,

2x1

 

2x2

3x3

6,

 

 

 

x1

4x2

2x3

2x4

3,

x j

 

0,

j

1, 2, 3.

 

 

 

x j

0,

j

1,...,4.

 

 

15. Z

x2

2x3 2x5

min

 

16. Z

x4

 

x5

 

max

 

 

x1

 

3x2

 

x3

2x5

7,

 

x1

 

 

x4

2x5

1,

 

 

 

2x2

 

4x3

x4

 

2,

 

x2

 

 

2x4

x5

2,

 

 

 

4x2

 

3x3

 

8x5 x6

6,

 

 

x3

 

3x4

x5

3,

 

x j

 

0,

j

1,...,6.

 

 

 

x j

0,

j

1,...,5.

 

 

5. МЕТОД ИСКУССТВЕННОГО БАЗИСА (М-задача)

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

Рассмотрим задачу в общем виде:

n

 

 

 

 

 

 

 

 

 

aij x j

ai0

i 1, m ,

(1)

j 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x j

0( j

1,n),

 

ai0

0.

 

n

 

 

 

 

 

 

 

Z

c j x j

max.

(2)

 

j 1

 

 

 

 

 

 

 

Пусть система (1) не является системой с базисом. Прибавим к левой части каждого уравнения системы (1) переменную уi ≥ 0, которую назовем искусственной. Система примет вид:

n

 

 

 

(3)

aij x j yi ai 0 i 1, m .

j1

(3)– система с базисом.

23

Составим новую целевую функцию:

n

 

m

 

 

T

c j x j M

yi

max.

(4)

j 1

 

i 1

 

 

Задача нахождения максимума функции (4) при ограничениях (3) называется М-задачей.

Замечание 1. Если исходная задача решается на минимум, то целевая функция М-задачи составляется так:

n

m

T

c j x j M yi min.

j 1

i 1

В обоих случаях М может принимать сколь угодно большое положительное значение.

Замечание 2. Искусственные неизвестные следует вводить только в те ограничения, которые не содержат базисных неизвестных.

Связь между решениями исходной и М-задачей устанавливается следующими теоремами.

Теорема 1. Если в оптимальном плане Y 1 , 2 ,..., n , 0,...,0 М-задачи все искусственные переменные равны нулю, то соответствующее решение X 1 , 2 ,..., n исходной задачи также является оптимальным.

Теорема 2. Если в оптимальном плане М-задачи хотя бы одна из искусственных переменных отлична от нуля, то исходная задача решения не имеет.

Алгоритм метода искусственного базиса имеет свои особенности:

1)симплексная таблица имеет две оценочные строки: М-строку и Z- строку. Оценка в М-задаче имеет вид: а + bМ, где М > 0 сколь угодно большое число. Следовательно, знак оценки определяется знаком коэффициента b. Число а записываем в Z-строку (первую строку оценки), а коэффициент b – в М-строку (вторую строку);

2)разрешающий столбец выбирается по оценкам М-строки;

3)если все искусственные переменные вышли из базиса, задача решается дальше обычным симплекс-методом;

4)если М-задача решена, но искусственные переменные не вышли из базиса, то исходная задача решения не имеет.

24

Пример 1.

Z 5x1

2x2 x3

max

2x1

x2

x3

5,

3x1

2x2

x3

6,

5x1

3x2

4x3

1,

x j

0,

j

1, 2, 3.

Преобразуем систему ограничений к системе уравнений:

2x1

x2

x3 x4

5,

3x1

2x2

x3

6,

5x1

3x2

4x3

x5 1.

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

2x1

x2

x3 x4

 

5,

3x1

2x2

x3

y1

6,

5x1

3x2

4x3

x5

y2 1.

Целевая функция М-задачи:

 

 

T 5x1

x2

x3 M y1 y2

max .

 

 

 

Составляем симплексную таблицу:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Сj

Б

0

 

5

 

2

 

-1

0

0

θ

ai 0

 

x1

 

x2

 

x3

x4

x5

 

 

 

 

 

 

0

x4

5

 

2

 

1

 

1

1

0

5/2

M

y1

6

 

3

 

2

 

1

0

0

2

M

y2

1

 

5

 

3

 

4

0

-1

1/5

 

Z

0

 

-5

 

-2

 

1

0

0

 

 

M

-7

 

-8

 

-5

 

-5

0

1

 

0

x4

23/5

 

0

 

-1/5

 

-3/5

1

2/5

23/2

M

y1

27/5

 

0

 

1/5

 

-7/5

0

3/5

27/3

5

x1

1/5

 

1

 

3/5

 

4/5

0

-1/5

-

 

Z

1

 

0

 

1

 

5

0

-1

 

 

M

-27/5

 

0

 

-1/5

 

7/5

0

-3/5

 

0

x4

1

 

0

 

-1/3

 

1/5

1

0

 

0

x5

9

 

0

 

1/3

 

-7/5

0

1

 

5

x1

2

 

1

 

2/3

 

1/3

0

0

 

 

Z

10

 

0

 

4/3

 

8/3

0

0

 

25

Оптимальный план: X îïò . 2, 0, 0 ,

Zmax 10.

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

Пример 2.

Z 2x1 3x2 max x1 x2 1,

3x1 2x2 6,

x1,2 0.

Вводим балансовые переменные:

x1

x2 x3

1,

3x1

2x2

x4 6.

Система не каноническая. Составляем М-задачу:

T 2x1 3x2 My max

x1

x2

x3

1,

3x1

2x2

x4

y 6,

x j

0, j

1,...,4,

 

y 0.

Решаем М-задачу симплексным методом:

Сj

Б

ai 0

2

3

0

0

θ

x1

x2

x3

x4

 

 

 

 

0

x3

1

1

1

1

0

1

M

y

6

3

2

0

-1

2

 

Z

0

-2

-3

0

0

 

 

M

-6

-3

-2

0

1

 

2

x1

1

1

1

1

0

 

M

y

3

0

-1

-3

-1

 

 

Z

2

0

-1

2

0

 

 

M

-3

0

1

3

1

 

М-задача решена (нет отрицательных оценок в М-строке), но в этом решении искусственная неизвестная y осталась в базисе, следовательно,

26

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

Решить следующие задачи:

1.

Z

x1

x2

max

 

2x1

x2

2,

 

 

x1

2x2

2,

 

 

x1

x2

5,

 

 

x1,2

0.

 

3.

Z

x1

2x2

min

x1 x2 8, x1 x2 2,

x1,2 0.

5.

Z

3x1

x2

x3

max

 

2x1

 

x2

x3

2,

 

 

x1

2x2

 

8,

 

 

x1

 

x2

 

5,

 

 

x j

0,

j

1, 2, 3.

 

7.

Z

4x1

2x2

x3

min

 

x1

 

3x2

x3

4,

 

 

 

 

2x2

x3

1,

 

 

x j

0,

j 1, 2, 3.

 

9.

Z

x1

2x2

min

 

x1 x2 8, x1 x2 2, x1 3,

x1,2 0.

11. Z 2x1 3x2

min

3x1 2x2 6, x1 4x2 4,

x1,2 0.

2. Z

3x1

x2 3x3

34x4 min

x1

2x2

x3

x4

0,

2x1

2x2

3x3 3x4

9,

x1

x2

2x3

x4

6,

x j

0, j

1,...,4.

 

 

4. Z

x1

max

 

 

x1 2x2 0, x1 x2 1, x1 x2 1,

x1,2 0.

6. Z x1 2x2 max

 

x1

x2

1,

 

 

 

x1

x2

5,

 

 

 

x1

 

 

1,

 

 

 

x1,2

0.

 

 

 

8. Z

2x1

2x2

3x3

max

x1

x2

2x3

1,

 

 

2x2

3x3

8,

 

x j

0,

j

1, 2, 3.

 

10. Z

x1

2x2

 

x3

max

 

x1

x2

2x3

4,

 

x1

2x2

 

 

3,

 

x1

4x2

 

 

3,

 

x j

0,

j

1, 2, 3.

12. Z

x1

4x2

x3

max

 

x1

x2

x3

14,

2x1

5x2

x3

0,

 

x j

0,

j

1, 2, 3.

27

13. Z 2x1 x2 max

14.

2x1 x2 1, x1 x2 4,

x1,2 0.

15. Z

2x1

x2

5x3

min

16.

x1

x2

x3

4,

 

 

x1

5x2

x3

5,

 

 

x j

0,

j 1, 2, 3.

 

 

Z 6x1 4x2 4x3 min x1 x2 2x3 1,

x1 2x2 2x3 6,

x j

0, j 1, 2, 3.

Z 3x1 2x2 4x3 min x1 x2 2x3 4,

3x1 x2 4x3 7,

x j

0, j 1, 2, 3.

6. ДВОЙСТВЕННОСТЬ

Пример 1. Рассмотрим задачу об оптимальном плане выпуска продукции: для изготовления 4 видов продукции используется 2 вида сырья. Запасы сырья и его расход на изготовление единицы каждого вида продукции даны в таблице:

Виды сырья

 

Запасы

 

Виды продукции

 

 

I

II

III

 

IV

 

 

 

 

А

 

160

4

3

1

 

1

Б

 

900

-

4

9

 

12

Доход

 

12

5

4

 

1

 

 

 

 

 

 

 

 

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

Математическая формулировка (модель) задачи: Максимизировать функцию

Z

12x1

5x2

4x3

x4 max,

при ограничениях

 

 

 

 

 

 

4x1

3x2

x3

x4

160,

 

y1

 

 

4x2

9x3

12x4

900,

 

y2

x j

0

j 1,...,4 .

 

 

 

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

28

Обозначим соответственно, через y1 и y2 цену единицы сырья А и Б. Производство продукции вида I приносит предприятию доход 12

денежных единиц. При этом расходуется 4 единицы сырья А и 0 единиц сырья Б. Выручка от продажи сырья, расходуемого на единицу продукции I по ценам y1 и y2 , составит

4 y1 0 y2 .

Эта величина должна быть не меньше тех доходов, которые предприятие получит от реализации продукции вида I, следовательно,

4 y1 0 y2 12.

Аналогичные рассуждения в отношении единицы продукции вида II, III, IV приводят к следующим неравенствам:

3y1

4

y2

5,

y1

9

y

2

4,

y1

12y

2

1.

Общая стоимость всех запасов сырья, приобретаемого организацией, составит W 160y1 900y2 .

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

Получим задачу:

W

160y1

900y2

 

min

4 y1

0 y2

12,

 

x1

 

3y1

4 y2

5,

 

x2

y1

9 y2

4,

 

x3

y1

12y2

1,

 

x4

y1

0, y2

0.

 

 

Получили задачу, двойственную данной. Следовательно, для стандартной задачи нужно выполнить следующие действия, для того чтобы получить ей двойственную:

1)число неизвестных в двойственной задаче равно числу ограничений в исходной;

2)неравенства в системе ограничений двойственной задачи будут противоположного смысла, чем неравенства в системе ограничений исходной задачи; сохраняется неотрицательность переменных;

3)свободные члены ограничений исходной задачи становятся коэффициентами целевой функции двойственной задачи, а коэффициенты целевой функции исходной задачи превращаются в свободные члены двойственной задачи;

4)в исходной задаче целевая функция минимизируется, а в двойственной – максимизируется.

29

По решению одной из задач можно сразу определить решение другой.

Решим исходную задачу симплекс-методом:

Сj

Б

0

12

5

4

1

0

0

θ

ai 0

x1

x2

x3

x4

x5

x6

 

 

 

0

x5

160

4

3

1

1

1

0

40

0

x6

900

0

4

9

12

0

1

-

 

Z =

0

-12

-5

-4

-1

0

0

 

 

 

 

 

 

 

 

 

 

 

12

x1

40

1

3/4

1/4

1/4

1/4

0

160

0

x6

900

0

4

9

12

0

1

100

 

Z =

480

0

4

-1

2

3

0

 

12

x1

15

1

23/36

0

-1/12

1/4

-1/36

 

4

x3

100

0

4/9

1

4/3

0

1/9

 

 

Z =

580

0

40/9

0

10/3

3

1/9

 

 

 

 

 

 

 

 

y1

y2

 

 

 

 

 

 

15, 0,100,0 .

Zmax 580 ден. ед.при X

Следовательно, для двойственной задачи

 

 

 

 

3,1/9 .

Wmin 580 ед.при Y

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

переменной (Сj), т.е. y1

3 0

3; y2

1/ 9

0

1/ 9 .

Проверим:

 

 

 

 

 

 

W 160y1

900y2

min .

При y1 3, y2 1/ 9 , Wmin

160 3

900 1/ 9

480

100 580 .

Пример 2. В двойственной задаче к основной переменные могут иметь любой знак. Составим двойственную задачу к основной:

Z 3x1 x2 x3 max

2x1

x2

x3

6,

y1

x1

2x2

x3

4,

y2

x j

0, j

1, 2, 3.

 

30

Двойственная задача имеет следующий вид:

W 6 y1 4 y2 min

2 y1

y

2

3,

x1

y1

2 y

2

1,

x2

y1

y2

1.

x3

Решим двойственную задачу графически

y2

3

 

А

2

n

y1

1

Координаты точки А дают значения неизвестных y1 и y2 , при которых функция W принимает минимальное значение.

Найдем координаты этой точки:

y1

2 y2

3,

y1

2 y2 3,

y1

y2

1,

y1

1 y2 ,

y1

3,

 

 

 

y2

2,

A

3, 2 ,

 

Wmin

6 3

4 2

18 8

26.

По решению двойственной задачи найдем решение исходной по второй

теореме двойственности (теореме равновесия).

 

Подставим координаты точки

 

A 3, 2

в систему ограничений

двойственной задачи:

 

 

 

2 y1

y2

3,

 

y1

2 y2

1,

 

y1

y2

1.

 

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]