Математическое программирование-1
.pdf
Стратегия A4 строго доминируется стратегией A3.
|
B1 |
B2 |
A1 |
1 |
4 |
A2 |
3 |
1 |
A3 |
2 |
3 |
Получили матрицу выигрышей, где у игроков A и B нет доминирующих стратегий.
Сначала найдем оптимальную стратегию игрока B.
Заметим, что 0 и решим следующую задачу линейного программирования.
Максимизировать функцию Z y1 y2 max приограничениях:
y1 4 y2 13y1 y2 1 .2 y1 3y2 1y1, y2 0
Так как переменных всего две, то эту задачу линейного программирования проще всего решить графическим способом. Построим множество допустимых решений, т. е. область, описываемую этими неравенствами и ограниченную прямыми:
y1 4 y2 |
1 |
|||||
3y |
y |
2 |
1 |
|
||
|
1 |
|
|
|
1 . |
|
2 y |
3y |
2 |
|
|||
|
1 |
|
|
|
|
|
y |
0; y |
2 |
|
0 |
||
1 |
|
|
|
|
|
|
По оси абсцисс будем откладывать y1, а по оси ординат y2. Область допустимых значений закрашена серым цветом на рис. 7.5.
110
Рис. 7.5
Строим линии уровня Z y1 y2, которые имеют вид y1 y2 C,
где C – произвольная постоянная. Для увеличения C прямая должна занимать максимально «высокое» положение, но имея с областью допустимых решений хотя бы одну точку, как показано на рис. 7.6.
111
Рис. 7.6
Такое положение прямой – проходящее через точку . Найдем ее координаты как пересечение двух прямых.
y |
2 |
3y |
|
1 |
2 |
|
1 |
|
||
|
|
1 |
|
|
1 3y1 1 |
|
9 y1 3 2 y1 1 |
|||
y |
2 |
|
2 y |
|
3 y1 |
3 |
||||
|
|
3 |
1 |
3 |
|
|
|
|
||
|
|
|
|
|
|
|
|
|
||
7 y1 2 y1 72 ; y2 3 72 1 76 1 71 .
112
|
Получили y |
|
2 |
; y |
2 |
1 |
, тогда |
|
|
1 |
|
|
1 |
|
7 |
||||
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
1 |
7 |
|
7 |
|
|
|
|
y1 y2 |
2 |
|
1 |
3 |
||||
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
q y 2 |
7 |
2 ; |
|
|
|
1 |
7 1. |
7 |
7 |
|
||||||||
и |
q y |
2 |
|
|
|
|
|
||||||||||||
|
1 |
1 |
7 |
|
3 |
3 |
|
2 |
|
7 |
3 |
3 |
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
Таким образом, оптимальная смешанная стратегия игрока B имеет вид
q* 2 , 1 , 0, 0, 0 .
3 3
Для нахождения оптимальной стратегии игрока А решаем следующую задачу линейного программирования.
Минимизировать функцию F x1 x2 x3 min при ограничениях
x1 3x2 2x3 1.4x1 x2 3x3 1
Для решения этой задачи применим симплекс-метод.
Запишем систему ограничений в каноническом для симплексметода виде:
x1 3x2 2x3 1.4x1 x2 3x3 1
Введем такие дополнительные переменные, что неравенства преобразуются в равенства
x1 3x2 2x3 x4 1.4x1 x2 3x3 x5 1
Первый опорный план: X1 0, 0, 0, 1, 1 и целевая функция
F 0 0 0 0.
113
Занесем данные в симплекс-таблицу:
Базис |
x1 |
x2 |
x3 |
x4 |
x5 |
b |
x4 |
–1 |
–3 |
–2 |
1 |
0 |
–1 |
x5 |
–4 |
–1 |
–3 |
0 |
1 |
–1 |
F(x) |
–1 |
|
–1 |
–1 |
|
0 |
|
|
|
|
|
|
|
Значение базисных переменных отрицательны, нас такая ситуация не устраивает. Выберем среди значений базисного столбца наибольший по модулю. Оба таких значения равны 1, выберем первую строку (делим элементы целевой строки на элементы выбранной):
|
1 |
|
1 |
|
1 |
|
|
1 |
|
|
min |
|
; |
|
; |
|
|
|
. |
||
1 |
3 |
2 |
3 |
|||||||
|
|
|
|
|
|
Это значение соответствует второму столбцу, т. е. в базисе x4 следует вывести, а x2 – ввести. Разрешающий элемент (–3), стоя-
щий на пересечении первой строки и второго столбца. Пересчитаем элементы таблицы по методу Жордана-Гаусса.
Базис |
x1 |
|
x2 |
|
x3 |
|
x4 |
x5 |
b |
|
||
x4 |
13 |
|
1 |
2 3 |
|
|
13 |
0 |
13 |
|
||
x |
4 1 |
3 |
0 |
3 1 |
3 |
1 |
3 |
1 |
1 1 |
3 |
||
5 |
|
|
|
|
|
|
|
|
||||
F(x) |
1 13 |
0 |
1 13 |
0 13 |
0 |
0 13 |
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
Базис |
x1 |
|
x2 |
|
x3 |
|
x4 |
x5 |
b |
|
||
x2 |
13 |
|
1 |
2 3 |
|
|
13 |
0 |
13 |
|
||
x5 |
13 |
|
0 |
|
7 3 |
|
|
13 |
1 |
2 3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
F(x) |
2 3 |
|
0 |
|
13 |
|
13 |
0 |
13 |
|
||
114
Среди значений базисных переменных присутствуют отрицательные значения, план не является оптимальным. Максимальное значение столбца базисных переменных (по модулю) 2/3, это соответ-
|
|
|
|
|
|
|
2 |
|
|
1 |
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
3 |
|
3 |
|
3 |
|
|
|
|
2 |
|
|
|
1 |
|
|
1 |
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
ствует второй строке, |
min |
|
|
; |
|
; |
|
|
|
min |
|
|
|
; |
|
|
|
;1 |
7 |
, что |
||||||||||||
|
11 |
|
7 |
|
|
1 |
11 |
|
7 |
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
3 |
|
|
3 |
|
|
|
3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
соответствует третьему столбцу, т. е. выводим из базиса |
x5 |
и вво- |
||||||||||||||||||||||||||||||
дим туда x3. Пересчитываем |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
Базис |
x1 |
|
x2 |
|
|
|
|
x3 |
|
|
|
|
|
x4 |
|
|
|
|
x5 |
|
|
|
b |
|||||||||
x |
5 |
7 |
|
1 |
|
|
|
|
0 |
|
|
|
|
|
3 |
7 |
|
|
|
|
2 |
7 |
|
|
|
1 |
7 |
|||||
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
x |
11 |
7 |
|
0 |
|
|
|
|
1 |
|
|
|
|
|
|
1 |
7 |
|
|
|
3 |
7 |
|
|
2 |
7 |
||||||
3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
17 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
F(x) |
|
0 |
|
|
|
|
|
|
|
|
|
|
0 2 7 |
|
|
|
17 |
|
|
37 |
||||||||||||
Среди значений базисного столбца нет отрицательных переменных, т. е. получен оптимальный план:
x 0, |
x |
2 |
1 |
, |
x 2 |
, |
F 0 1 2 3 . |
|||
1 |
|
7 |
|
3 |
7 |
|
7 |
7 |
7 |
|
|
|
|
|
|
|
|||||
Тогда |
1 |
|
7 |
; |
p |
x 1 |
7 |
1; p |
x 2 |
7 |
2. |
||||||
|
|
||||||||||||||||
|
x1 x2 x3 7 |
|
2 |
|
2 |
|
7 3 3 |
3 |
3 |
7 |
3 |
3 |
|||||
|
|
|
|
|
|
|
|
|
|||||||||
Следовательно, |
оптимальная |
смешанная |
стратегия игрока A |
||||||||||||||
имеет вид: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
p |
* |
|
0, |
1 |
, |
2 |
, 0 |
|
|
|
|
|
|
|
|
|
|
|
|
3 |
3 |
. |
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
7.3. Контрольные вопросы
1. Какое основное назначение теории игр?
2. Что называется стратегией игры?
3. Вчем сутьчистыхи смешанныхстратегий икаковы ихсвойства?
115
4.Каковы понятия максиминной и минимаксной стратегий?
5.В чем заключается приведение матричной игры к задаче линейного программирования?
6.Каков способ нахождения оптимальной стратегии игрока с помощью пакета MS Excel?
7.Что означает оптимальная смешанная стратегия игрока?
8.Как определить оптимальные стратегии игроков геометрически и с помощью симплекс-метода?
7.4. Задания для самостоятельной работы
Решить матричные игры, имеющие платежные матрицы вида:
|
8 |
4 |
2 |
|
|
|
1 |
1 1 |
|
1 |
2 5 3 |
2 |
|||||||||||||||
1. |
|
2 |
8 |
4 |
|
|
|
2. |
|
|
2 |
2 |
|
|
2 |
|
3. |
|
1 |
4 7 2 |
4 |
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
1 |
2 |
8 |
|
|
|
|
|
|
3 |
3 |
|
|
|
|
|
|
|
|
5 1 1 |
1 3 |
|
||||
|
|
|
|
|
|
|
|
|
|
3 |
|
|
|
|
|||||||||||||
|
0 |
13 |
|
1 |
|
1 |
0 1 |
|
|
3 2 |
4 |
|
|
|
|
||||||||||||
4. |
|
|
|
|
|
|
|
5. |
|
0 |
2 1 |
|
|
6. |
|
4 3 |
2 |
|
|
|
|
|
|||||
13 0 |
|
13 |
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
1 13 |
|
0 |
|
|
|
1 |
1 |
|
3 |
|
|
|
|
2 4 |
3 |
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
3 |
6 |
0 |
|
|
|
|
|
3 |
0 |
7 |
|
|
|
|
|
203 |
403 |
103 |
|
|
|
||||
7. |
|
5 |
3 |
2 |
|
|
|
8. |
|
4 |
6 |
0 |
|
|
|
9. |
|
303 |
3 |
|
103 |
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
2 |
1 |
6 |
|
|
|
|
|
3 |
4 |
3 |
|
|
|
|
|
|
3 |
103 |
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
303 |
|
|
||||||||||||
|
2 |
|
11 |
1 |
|
|
7 |
5 4 |
|
|
|
16 |
0 14 |
|
|
|
|||||||||||
10. |
|
|
2 |
|
|
|
11. |
|
1 |
3 |
7 |
|
|
|
12. |
|
6 |
6 |
|
|
|
|
|
||||
15 |
|
|
11 |
|
|
|
|
|
16 |
|
|
|
|||||||||||||||
|
|
|
|
|
2 |
|
|
|
2 |
7 4 |
|
|
|
|
|
6 |
12 2 |
|
|
|
|
||||||
|
3 15 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
Матричную игру решить с помощью симплекс-метода и пакета
MS Excel.
116
ЗАНЯТИЕ 8. ТРАНСПОРТНАЯ ЗАДАЧА (ТЗ)
8.1. Краткие теоретические сведения
Одним из частных случаев задач линейного программирования является транспортная задача. Цель данной задачи состоит в определении оптимального плана перевозок, при котором минимальной была либо стоимость перевозок, либо минимальным было время доставки груза.
В общем виде транспортную задачу можно представить следующим образом: в m пунктах производства A1, A2 , , Am имеется
однородный груз в количествах, соответственно, a1, a2, , am. Этот груз необходимо доставить в n пунктов назначения B1, B2 , , Bn в количествах, соответственно, b1, b2 , , bn. Стоимость перевозки 1 ед. груза (тариф) из пункта Ai в пункт B j равна cij.
Требуется составить план перевозок так, чтобы:
1)мощности (запасы) всех поставщиков были реализованы;
2)заказы всех потребителей были удовлетворены;
3)суммарные затраты на перевозку были бы минимальны. Исходные данные транспортной задачи записываются в виде
таблицы (табл. 8.1).
Таблица 8.1
117
Виды транспортных задач
В зависимости от соотношения между суммарными запасами груза у поставщиков и суммарными потребностями в нем потребителей транспортные задачи делятся на два вида: закрытые и от-
крытые.
Определение 8.1. Если сумма запасов груза у поставщиков равна суммарной потребности в нем потребителей, т. е.
m |
n |
|
ai bj , |
(8.1) |
|
i 1 |
j 1 |
|
то транспортная задача называется закрытой. Определение 8.2. Если условие (8.1) не выполняется, т. е.
m |
n |
ai bj , |
|
i 1 |
j 1 |
то транспортная задача называется открытой.
Теорема 8.1 (необходимое и достаточное условие разрешимо-
сти ТЗ). Для разрешимости транспортной задачи необходимо и достаточно, чтобы запасы груза в пунктах отправления были равны потребностям в грузе в пунктах назначения, т. е. чтобы выполнялось равенство (8.1).
8.1.1. Математическая модель закрытой транспортной задачи
Рассмотрим закрытую транспортную задачу.
1. Неизвестными транспортной задачи являются xij – объемы
перевозок от каждого i-го поставщика каждому j-му потребителю. Эти переменные также можно записать в виде матрицы перевозок:
x11 |
x12 |
|
x1n |
|||
x |
21 |
x |
22 |
|
x |
|
X |
|
|
|
2n . |
||
|
|
|
|
|||
|
|
|
xm2 |
|
|
|
xm1 |
xmn |
|||||
118
2. Так как произведение cij xij определяет затраты на перевозку груза от i-го поставщика j-му потребителю, то суммарные затраты
m n
на перевозку всех грузов равны cij xij . Т. е. целевая функция
i 1 j 1
будет иметь вид
m n
F cij xij min. i 1 j 1
3. Система ограничений транспортной задачи состоит из двух групп уравнений.
Первая группа состоит из m уравнений и описывает тот факт, что запасы всех поставщиков A a1, a2, , am вывозятся полностью, т. е.:
n
xij ai ; i 1,m. (8.2)
j 1
Вторая группа состоит из n уравнений и выражает требование полностью удовлетворить запросы всех n потребителей:
B b1, b2 , , bn , т. е.
m |
|
|
|
|
|
|
||
xij bj ; |
j |
1,n |
. |
|
(8.3) |
|||
i 1 |
|
|
|
|
|
|
||
Кроме того, переменные xij должны быть неотрицательны: |
|
|||||||
xij 0; i |
|
; |
j |
|
. |
|
||
1,m |
1,n |
|
||||||
4. Оптимальным решением транспортной задачи является матрица размерности m n, удовлетворяющая системе ограничений
и обеспечивающая минимум целевой функции.
119
