- •"Томский политехнический университет"
- •Предисловие
- •Тема 1 Предмет экономико-математического моделирования
- •Моделирование как метод научного познания
- •Классификация экономико-математических моделей
- •Этапы экономико-математического моделирования
- •Взаимосвязи этапов
- •Моделирования
- •Тема 2 Системный подход к изучению экономических явлений Системный анализ как научная дисциплина
- •Вычислительная техника в системном анализе
- •Системный подход Основные определения: элементы, связи, система
- •Принципы системного подхода
- •Об использовании принципов системного подхода
- •Тема 3 Математические методы
- •И основные классы задач оптимизации
- •Общая постановка математической модели задач
- •Оптимизации
- •Тема 4 Линейное программирование
- •Пример решения станковой задачи
- •Симплекс-метод решения задач линейного программирования
- •Свойства опорных решений
- •Решение задач линейного программирования симплекс-методом
- •Конечность симплекс-метода
- •Метод искусственного базиса для отыскания начального опорного решения
- •Двойственность в линейном программировании
- •Виды математических моделей двойственных задач
- •Тема 5 Целочисленное программирование
- •Постановка задачи и метод решения
- •Метод Гомори
- •Составление дополнительного ограничения (сечения Гомори)
- •Тема 6 Транспортная задача
- •Построение первоначального опорного плана
- •Метод минимальной стоимости
- •Определение оптимального плана транспортных задач, имеющих некоторые усложнения в их постановке
- •Тема 7 Нелинейное программирование
- •Теорема Куна – Таккера
- •Тема 8 Регрессионный анализ
- •Тема 9 Игровые методы обоснования решений
- •Основные термины
- •Постановка задачи и выбор критерия оптимизации
- •Построение математической модели
- •Исследование математической модели
- •0Ропт.1; 0qопт.1.
- •Упрощение платёжной матрицы
- •Тема 10 Основы сетевого планирования и управления
- •Параллельности работ
- •Временные параметры сетевого графика
- •Алгоритм расчёта ранних сроков начал и окончаний работ
- •Критическое время и критический путь
- •Алгоритм построения критического пути
- •Исследование сетевой модели
- •Оптимизация сетевых моделей
- •Тема 11 Задачи упорядочения. Задачи управления запасами. Задачи замены оборудования
- •Классификация задач упорядочения
- •Детерминированная задача упорядочения Постановка задачи и выбор критерия оптимизации
- •Выявление основных особенностей, взаимосвязей и количественных закономерностей
- •Построение математической модели
- •Исследование математической модели
- •Задачи управления запасами
- •Классификация задач управления запасами
- •Однопродуктовая детерминированная задача управления запасами Постановка задачи и выбор критерия оптимизации
- •Выявление основных особенностей, взаимосвязей и количественных закономерностей
- •Построение математической модели
- •Исследование математической модели
- •Задача управления запасами с учётом убытков
- •Постановка задачи
- •Выявление основных особенностей, взаимосвязей и количественных закономерностей
- •Выявление основных особенностей, взаимосвязей и количественных закономерностей
- •Построение математической модели
- •Исследование математической модели
- •Задачи замены оборудования
- •Классификация задач замены оборудования
- •Задача замены оборудования длительного пользования Постановка задачи. Выбор критерия оптимизации
- •Выявление основных особенностей, взаимосвязей и количественных закономерностей
- •Построение математической модели
- •Исследование математической модели
- •Задача замены оборудования с целью предупреждения отказа Постановка задачи и выбор критерия оптимизации
- •Выявление основных особенностей, взаимосвязей и количественных закономерностей
- •Построение математической модели
- •Исследование и решение математической модели
- •Тема 12 Задачи массового обслуживания
- •Классификация смо
- •Задачи анализа одноканальных систем массового обслуживания
- •Задача анализа детерминированной системы Постановка задачи
- •Выявление основных особенностей, взаимосвязей и количественных закономерностей
- •Построение математической модели
- •Выявление основных особенностей, взаимосвязей и количественных закономерностей
- •Построение математической модели
- •Исследование математической модели
- •Задача анализа замкнутой системы с ожиданием (потоки требований пуассоновские) Постановка задачи
- •Выявление основных особенностей, взаимосвязей и количественных закономерностей
- •Построение математической модели
- •Исследование и решение математической модели
- •Тема 13 Балансовые методы согласования
- •Ресурсов и потребностей
- •Анализ хозяйственных связей с помощью моделей
- •Межотраслевого баланса
- •Принципиальная схема межотраслевого баланса
- •Экономико-математическая модель межотраслевого баланса
- •2. Определить объёмы валовой продукции отраслей x1, x2,…, Xn по заданным объёмам конечного продукта y1, y2,…,Yn по формуле
- •Пример построения экономико-математической модели межотраслевого баланса и его расчёта для случая трёх отраслей
- •Экономическая природа коэффициентов прямых и полных затрат и их расчёт
- •Тема 14 Многокритериальные задачи
- •Классификация методов многокритериальной оценки альтернатив
- •Пример определения конкурентоспособности наукоемкой продукции на основе показателя “значимость технического решения” порогами несравнимости
- •Тема 15 Моделирование в условиях нечеткой информации
- •Нечеткие высказывания Нечеткими высказываниями называют высказывания следующего вида:
- •Тема 16 Моделирование процесса принятия решений
- •Интегральная модель определения конкурентоспособности продукции
- •Определение нечетких коэффициентов весомости критериев оценки конкурентоспособности продукции
- •Математическая модель рейтинговой оценки конкурентоспособности продукции
- •Отбор кандидатов в эксперты методом многокритериального выбора альтернатив с использованием правила нечеткого логического вывода
- •Заключение
- •Список литературы
- •Оглавление
Постановка задачи и выбор критерия оптимизации
Пусть в игре принимают участие две стороны: А и В. Условия игры заданы платёжной матрицей (табл. 9.2).
Таблица 9.2
Общий вид платёжной матрицы игры
-
A
B
B1
…
Bj
…
Bn
A1
c11
…
c1j
…
c1n
…
…
…
…
…
…
Ai
ci1
…
cij
…
cin
…
…
…
…
…
…
Am
cm1
…
cmj
…
cmn
Требуется определить оптимальные вероятности использования стратегий, максимизирующих средний выигрыш первой стороны.
Выбранную стороной А стратегию будем соответственно обозначать А1, А2 ,…, Аm, аналогично стратегии стороны В – символами В1, В2, …, Bn; pi -- вероятность использования стратегии i первой стороной A, qj – вероятность использования стратегии j второй стороной В, при этом .
Построение математической модели
Средний выигрыш первой стороны А:
.
Ограничения:
.
Исследование математической модели
Для получения оптимальных вероятностей использования стратегий А1,…,Аi, …,Аm (B1,,…,Bj,…,Bn) необходимо взять соответствующие частные производные от целевой функции, приравнять их к нулю и решить систему уравнений:
Для примера исследуем матрицу размерности 2х2:
Дана матрица выигрышей:
-
Аi
Bj
B1
B2
A1
с11
c12
A2
с21
c22
Требуется определить оптимальные смешанные стратегии:
.
Математическая модель этой задачи:
,
при этом p1+p2=1, q1+q2=1.
Обозначим р1=р, а р2=1-р, аналогично q1=q; q2=1-q. Подставляя полученные выражения в математическую модель и беря соответствующие частные производные, находим:
При этом должны выполняться необходимые условия:
0Ропт.1; 0qопт.1.
Значения ропт., qопт., меньшие нуля, необходимо считать равными нулю, а большие единицы – равными единице.
Подставляя оптимальные значения ропт., qопт., 1-ропт., 1-qопт. в целевую функцию
,
получаем максимальный средний выигрыш первой стороны:
Анализируя эту формулу, видим, что при с11с22-с12с21=0, т.е. с11/с12=с21/с22, игра становится безобидной. Существует так называемая основная теорема теории игр, состоящая в следующем: каждая конечная игра имеет, по крайней мере, одно решение, возможно, в области смешанных стратегий.
Применение смешанной стратегии позволяет получить выигрыш, равный цене игры: аvb.
Для оптимальных стратегий игроков имеет место соотношение
.
Применение игроком А оптимальной стратегии должно обеспечивать ему при любых действиях игрокаВ выигрыш не меньше цены игры v.
Поэтому должны выполняться следующие соотношения:
. |
(9.1) |
Аналогично, для игрока В оптимальная стратегия игрока должна обеспечивать при любых стратегиях игрокаА проигрыш, не превышающей величину v, т.е. справедливо соотношение:
. |
(9.2) |
В дальнейшем соотношения (9.1) и (9.2) используются для решения игры. Вообще, задача решения игры, если её матрица не содержит седловой точки, тем сложнее, чем больше значения m и n.
Стратегии игроков А и В, для которых вероятности рi и qj отличны от нуля, называются активными.
Методы решения задач теории игр во многом зависят от условий задачи и от вида матрицы выигрышей первого игрока.
Как уже упоминалось, если матрица С имеет седловую точку, то решение игры сводится к нахождению седловой точки матрицы С. Оптимальные стратегии игроков определяются при этом координатами (i0, j0) седловой точки матрицы С, а цена игры – элементом .
Если матрица С имеет размер mx2 или 2xn, то решение задачи может быть получено графически.
Поиллюстрируем на примере.
Пример 9.2.
C=
В плоскости переменных (p v) построим vj(p) – ожидаемый средний выигрыш первого игрока, применяющего первую стратегию с вероятностью р при условии, что второй игрок отвечает чистой стратегией
vj(j=1,2,3): v1(p)=3p+1(1-p)=2p+1,
v2(p)=2p+9(1-p)=-7p+9, 0p1
v3(p)=p+20(1-p)=-19p+20
v
20
с23 v2(p)
с11 р1+р2+р3=1
10 A 3 2р+1=7р+9,
9р=8, р=8/9.
с22 v1(p) 2 с12
1 P2 P3 р*=(1/9, 1/9);
1 q*=(7/9, 2/9, 0)
с21 P1 1-P1 с13 v*=25/9.
P
0 1
В точке А пересечения прямых v1(p) и v2(p) гарантированный выигрыш первого игрока, изображённый жирной линией, достигает наибольшего значения.
Следовательно, игрок А применяет стратегию А1 с вероятностью 8/9, а стратегию А2 - с вероятностью 1/9. Игрок В применяет стратегию В1 с вероятностью 7/9 или В2 с вероятностью 2/9. При этом его выигрыш (игрока А) в среднем составляет 25/9 единиц. Столько же составляет средний проигрыш игрока В.
Пример 9.3. Найти решение игры, заданной матрицей:
С=.
Матрица имеет размерность 4х2, поэтому решение задачи находим для игрока В. Рассматриваем задачу на минимизацию верхней границы проигрыша для игрока В. Следует отметить, что геометрические построения имеет смысл использовать только для определения активных стратегий игроков, а не для определения численных значений. Затем решение игры можно получить с помощью формул (*), (**) и (***). Формулы эти можно использовать, так как из соответствующей матрицы исключаются все стратегии, кроме активных, тогда она содержит 2 строки и 2 столбца. Отрезок NK определяет цену игры. Активными стратегиями для игрока А являются первая и четвёртая.
v
q1=
q2=Q*=;p1=
p2=
P*=
v*=
6 А4
5 А3
А1 4 K C42 L 4 А2
3 А1
А2 2 v*
C12
А3 N
0 q2 q1 1 q
-1
Pешение игр с матрицей mxn, если min {m; n}>2, рассмотрим позднее.
Рассмотрим метод Брауна приближённого решения задач теории игр, он подходит и для матриц, размерностью больше 2.
Метод Брауна представляет собой модель практического взаимного обучения игроков, при котором каждый игрок, анализируя способ поведения противника, старается ответить наилучшим способом. Этот метод относится к итерационным, т.е. получающим решение за какое-то количество шагов, причём точность вычисления зависит от числа ходов и от выбора начальной строки.
Первый игрок выбирает одну из своих стратегий. Второй игрок отвечает стратегией, которая минимизирует выигрыш первого игрока. Каждый игрок отвечает на очередной ход противника той своей стратегией, которая является оптимальной относительно всех предыдущих ходов противника.
Пример 9.4.
Рассмотрим решение задачи с матрицей
С=.
Первый игрок пусть выбирает i0=3. Второй игрок, просматривая 3-ю строку, выбирает с33=1 – наименьший элемент строки и отмечает его «*». Справа от матрицы выписываем элементы 3-го столбца. Первый игрок, зная выбор второго и просматривая выделенный столбец, выбирает в нём максимальный элемент (отмеченный звёздочкой) и назначает i1=1. Второй строкой под матрицей выписываем суммарный выигрыш первого игрока за 2 хода.
Минимизируя выигрыш первого игрока за два хода, второй игрок выбирает j2=2 и т.д.
Частота выбора строк и столбцов определяет компоненты приближённого решения, элементы со звёздочкой определяют суммарный выигрыш за n ходов.
Вычисления в итерационном методе Брауна удобно располагать, как в приведённой таблице:
Таблица 9.3
20 9 6 9* 9 9 9 15 21* 23* 23* 23* 23 23 23 25 28
1 3 6 0 6 9* 12* 15* 15* 15 16 19 22 25* 28* 31* 32* 31
4 2 1 3 1 3 5 7 10 13 17 19 21 23 25 27 31 35*
4 2 1* 3 В левом верхнем углу записана матрица игры – в
6 2* 10 9 данном случае матрица размера 3х4. Строки ниже
7 5* 16 9 матрицы содержат сложенный за n партий возможный
8 8* 22 9 выигрыш первого игрока, где n – количество партий
911 28 9* и номер строки. Столбцы правее матрицы представ-
10 14 34 9* ляют собой сложенный за n партий возможный проиг-
12* 14 43 15 рыш второго игрока, причём n – количество партий и 14 14*52 21 номер столбца. Пусть после n партий первый игрок
16 14* 62 27 обнаружил, что его противник выбрал j-ю чистую -
18 14* 70 33 стратегию (по количеству «звёздочек» в j-м столбце
19 17* 76 33 под матрицей) Sj раз (j=). На этом основании пер- 20 20* 82 33 вый игрок допускает, что его противник придержива-
21* 23 88 33 ется смешанной стратегии
22*26 94 33
и выбирает в последующей партии чистую стратегию, дающую максимальный средний выигрыш qn. Номер чистой стратегии – это номер максимальной компоненты n-го столбца.
Аналогично, второй игрок после n партий предполагает, что смешанная стратегия соперника , гдеti – частота появления i-й стратегии в предыдущих партиях. В (n+1)-й партии второй игрок выбирает чистую стратегию, обеспечивающую ему минимальный средний проигрыш. Номер этой чистой стратегии – это номер столбца минимального элемента n-й строки.
Степень приближения к решению зависит от выбора начальной строки и от числа ходов, поэтому:
при n=10 при n=14
Можно решить задачу теории игр, сведя математическую игру к задаче линейного программирования.
В самом деле, рассмотрим игру, матрица С которой имеет размерность mxn:
C=.
Пусть матрица не содержит седловой точки, поэтому решение игры представлено в смешанных стратегиях: Р=(p1, p2,…, pm); Q=(q1, q2,…, qn). При оптимальной стратегии игрока А выполняется условие, как мы уже выяснили (9.1, 9.2):
,
а оптимальной стратегии игрока В удовлетворяет условие:
.
Таким образом, можно рассмотреть задачу отыскания оптимальной стратегии игрока А, для которой имеют место следующие ограничения:
|
(9.3) |
Величина v (цена игры) неизвестна, однако можно предположить, что v0, если все элементы платежной матрицы неотрицательны, а этого всегда можно достигнуть, как уже упоминалось, прибавляя ко всем элементам матрицы некоторое положительное число.
Преобразуем систему ограничений, разделив все члены неравенств на v. В результате получим:
|
(9.4) |
где ti = pi/v, i=1,2, …, m.
Из условия p1+p2+ … +pm=1 следует, что t1+t2+ … +tm=1/v.
Решение игры должно максимизировать значение v, значит, функция должна принимать минимальное значение.
Таким образом, получена задача линейного программирования: при ограничениях (9.4) и дополнительных условиях неотрицательности переменныхti(i=). Решая её, находим значение ti и величину 1/v, затем отыскиваем значение pi=vti.
Для определения стратегии игрока В запишем следующие условия:
|
(9.5) |
Разделив все члены неравенств на v, получим:
|
(9.6) |
где uj=qj/v, j=1, ... , n. Переменные u1, u2, …, un должны быть выбраны так, чтобы выполнялись условия (9.6) и достигался максимум функции
W=u1+u2+…+un=1/v, uj0 (i=1,…,.m). |
(9.7) |
Таким образом, для решения игры имеем пару двойственных симметричных задач линейного программирования. Используя свойство симметричности, можно решить одну из них, требующую меньших вычислений, а решение второй задачи найти на основании оптимального плана двойственной.