ЭММ. Вариант 5. Задача 1 Графический метод решения задач линейногопрограммирования. Составить математическую модель по условию задачи
Скачать 171.28 Kb.
|
Целевая функция F= 8750 . Задача 4 Теория игр Найти оптимальные стратегии и цену игры, заданной платежной матрицей графическим методом. Решение: Проверяем, имеет ли платежная матрица седловую точку. Если да, то выписываем решение игры в чистых стратегиях. Считаем, что игрок I выбирает свою стратегию так, чтобы получить максимальный свой выигрыш, а игрок II выбирает свою стратегию так, чтобы минимизировать выигрыш игрока I.
Находим гарантированный выигрыш, определяемый нижней ценой игры a = max(ai) = 4, которая указывает на максимальную чистую стратегию A2. Верхняя цена игры b = min(bj) = 5. Что свидетельствует об отсутствии седловой точки, так как a ≠ b, тогда цена игры находится в пределах 4 ≤ y ≤ 5. Находим решение игры в смешанных стратегиях. Объясняется это тем, что игроки не могут объявить противнику свои чистые стратегии: им следует скрывать свои действия. Игру можно решить, если позволить игрокам выбирать свои стратегии случайным образом (смешивать чистые стратегии). 2. Проверяем платежную матрицу на доминирующие строки и доминирующие столбцы. Иногда на основании простого рассмотрения матрицы игры можно сказать, что некоторые чистые стратегии могут войти в оптимальную смешанную стратегию лишь с нулевой вероятностью. Говорят, что i-я стратегия 1-го игрока доминирует его k-ю стратегию, если aij ≥ akj для всех j Э N и хотя бы для одного j aij >akj. В этом случае говорят также, что i-я стратегия (или строка) – доминирующая, k-я – доминируемая. Говорят, что j-я стратегия 2-го игрока доминирует его l-ю стратегию, если для всех j Э M aij ≤ ail и хотя бы для одного i aij il. В этом случае j-ю стратегию (столбец) называют доминирующей, l-ю – доминируемой. С позиции проигрышей игрока В стратегия B1 доминирует над стратегией B4 (все элементы столбца 1 меньше элементов столбца 4), следовательно, исключаем 4-й столбец матрицы. Вероятность q4 = 0.
В платежной матрице отсутствуют доминирующие строки. Мы свели игру 2 x 4 к игре 2 x 3. Так как игроки выбирают свои чистые стратегии случайным образом, то выигрыш игрока I будет случайной величиной. В этом случае игрок I должен выбрать свои смешанные стратегии так, чтобы получить максимальный средний выигрыш. Аналогично, игрок II должен выбрать свои смешанные стратегии так, чтобы минимизировать математическое ожидание игрока I. 3. Находим решение игры в смешанных стратегиях. Решим задачу геометрическим методом, который включает в себя следующие этапы: 1. В декартовой системе координат по оси абсцисс откладывается отрезок, длина которого равна 1. Левый конец отрезка (точка х = 0) соответствует стратегии A1, правый - стратегии A2 (x = 1). Промежуточные точки х соответствуют вероятностям некоторых смешанных стратегий S1 = (p1,p2). 2. На левой оси ординат откладываются выигрыши стратегии A1. На линии, параллельной оси ординат, из точки 1 откладываются выигрыши стратегии A2. Решение игры (2 x n) проводим с позиции игрока A, придерживающегося максиминной стратегии. Доминирующихся и дублирующих стратегий ни у одного из игроков нет. Выделяем нижнюю границу выигрыша B1NB3. Максиминной оптимальной стратегии игрока A соответствует точка N, лежащая на пересечении прямых B1B1 и B3B3, для которых можно записать следующую систему уравнений: y = 4 + (5 - 4)p2 y = 6 + (4 - 6)p2 Откуда p1 = 1/3 p2 = 2/3 Цена игры, y = 14/3 Теперь можно найти минимаксную стратегию игрока B, записав соответствующую систему уравнений, исключив стратегию B2, которая дает явно больший проигрыш игроку B, и, следовательно, q2 = 0. 4q1+6q3 = y 5q1+4q3 = y q1+q3 = 1 или 4q1+6q3 = 14/3 5q1+4q3 = 14/3 q1+q3 = 1 Решая эту систему, находим: q1 = 2/3. q3 = 1/3. Ответ: Цена игры: y = 14/3, векторы стратегии игроков: Q(2/3, 0, 1/3), P(1/3, 2/3) 4. Проверим правильность решения игры с помощью критерия оптимальности стратегии. ∑aijqj ≤ v ∑aijpi ≥ v M(P1;Q) = (4*2/3) + (3*0) + (6*1/3) = 4.667 = v M(P2;Q) = (5*2/3) + (6*0) + (4*1/3) = 4.667 = v M(P;Q1) = (4*1/3) + (5*2/3) = 4.667 = v M(P;Q2) = (3*1/3) + (6*2/3) = 5 ≥ v M(P;Q3) = (6*1/3) + (4*2/3) = 4.667 = v Все неравенства выполняются как равенства или строгие неравенства, следовательно, решение игры найдено верно. Поскольку из исходной матрицы были удалены и столбцы, то найденные векторы вероятности можно записать в виде: P(1/3,2/3) Q(2/3,0,1/3,0): Задача 5 Сетевые графики Пусть необходимо выполнить комплекс взаимосвязанных работ. Используя данные вариантов задания, постройте сетевой график, найдите критический путь, посчитайте критическое время выполнения комплекса работ и другие временные характеристики сетевого графика. Исходные данные представлены по вариантам. Постройте график с временными характеристиками. Сделайте выводы.
Решение: Резерв времени события показывает, на какой допустимый период времени можно задержать наступление этого события, не вызывая при этом увеличения срока выполнения комплекса работ. Для определения резервов времени по событиям сети рассчитывают наиболее ранние tp и наиболее поздние tп сроки свершения событий. Любое событие не может наступить прежде, чем свершаться все предшествующие ему события и не будут выполнены все предшествующие работы. Поэтому ранний (или ожидаемый) срок tp(i) свершения i-ого события определяется продолжительностью максимального пути, предшествующего этому событию: tp(i) = max(t(Lni)) где Lni – любой путь, предшествующий i-ому событию, то есть путь от исходного до i-ого события сети. Если событие j имеет несколько предшествующих путей, а следовательно, несколько предшествующих событий i, то ранний срок свершения события j удобно находить по формуле: tp(j) = max[tp(i) + t(i,j)] Задержка свершения события i по отношению к своему раннему сроку не отразится на сроке свершения завершающего события (а значит, и на сроке выполнения комплекса работ) до тех пор, пока сумма срока свершения этого события и продолжительности (длины) максимального из следующих за ним путей не превысит длины критического пути. Поэтому поздний (или предельный) срок tп(i) свершения i-ого события равен: tп(i) = tkp - max(t(Lci)) где Lci - любой путь, следующий за i-ым событием, т.е. путь от i-ого до завершающего события сети. Если событие i имеет несколько последующих путей, а следовательно, несколько последующих событий j, то поздний срок свершения события i удобно находить по формуле: tп(i) = min[tп(j) - t(i,j)] Резерв времени R(i) i-ого события определяется как разность между поздним и ранним сроками его свершения: R(i) = tп(i) - tp(i) Резерв времени события показывает, на какой допустимый период времени можно задержать наступление этого события, не вызывая при этом увеличения срока выполнения комплекса работ. Критические события резервов времени не имеют, так как любая задержка в свершении события, лежащего на критическом пути, вызовет такую же задержку в свершении завершающего события. Таким образом, определив ранний срок наступления завершающего события сети, мы тем самым определяем длину критического пути. При определении ранних сроков свершения событий tp(i) двигаемся по сетевому графику слева направо и используем формулы (1), (2). Расчет сроков свершения событий. Для i=0 (начального события), очевидно tp(0)=0. i=1: tp(1) = tp(0) + t(0,1) = 0 + 7 = 7. i=2: tp(2) = tp(0) + t(0,2) = 0 + 8 = 8. i=3: tp(3) = tp(1) + t(1,3) = 7 + 5 = 12. i=4: max(tp(1) + t(1,4);tp(2) + t(2,4)) = max(7 + 4;8 + 7) = 15. i=5: tp(5) = tp(2) + t(2,5) = 8 + 8 = 16. i=6: max(tp(3) + t(3,6);tp(4) + t(4,6);tp(5) + t(5,6)) = max(12 + 9;15 + 10;16 + 9) = 25. Длина критического пути равна раннему сроку свершения завершающего события 6: tkp=tp(6)=25 При определении поздних сроков свершения событий tп(i) двигаемся по сети в обратном направлении, то есть справа налево и используем формулы (3), (4). Для i=6 (завершающего события) поздний срок свершения события должен равняться его раннему сроку (иначе изменится длина критического пути): tп(6)= tр(6)=25 Далее просматриваются строки, оканчивающиеся на номер предпоследнего события, т.е. 5. Просматриваются все строчки, начинающиеся с номера 5. i=5: tп(5) = tп(6) - t(5,6) = 25 - 9 = 16. Далее просматриваются строки, оканчивающиеся на номер предпоследнего события, т.е. 4. Просматриваются все строчки, начинающиеся с номера 4. i=4: tп(4) = tп(6) - t(4,6) = 25 - 10 = 15. Далее просматриваются строки, оканчивающиеся на номер предпоследнего события, т.е. 3. Просматриваются все строчки, начинающиеся с номера 3. i=3: tп(3) = tп(6) - t(3,6) = 25 - 9 = 16. Далее просматриваются строки, оканчивающиеся на номер предпоследнего события, т.е. 2. Просматриваются все строчки, начинающиеся с номера 2. i=2: min(tп(4) - t(2,4);tп(5) - t(2,5)) = min(15 - 7;16 - 8) = 8. Далее просматриваются строки, оканчивающиеся на номер предпоследнего события, т.е. 1. Просматриваются все строчки, начинающиеся с номера 1. i=1: min(tп(3) - t(1,3);tп(4) - t(1,4)) = min(16 - 5;15 - 4) = 11. Далее просматриваются строки, оканчивающиеся на номер предпоследнего события, т.е. 0. Просматриваются все строчки, начинающиеся с номера 0. i=0: min(tп(1) - t(0,1);tп(2) - t(0,2)) = min(11 - 7;8 - 8) = 0. Таблица 1 - Расчет резерва событий
В данном случае имеются несколько критических путей: Критический путь №1:(0,2)(2,4)(4,6) Критический путь №2:(0,2)(2,5)(5,6) Продолжительность критического пути: 25 Рис. 1 - Сетевой график с нанесенными временными характеристиками. Список использованной литературы Вентцель, Е. С. Исследование операций. Задачи, принципы, методология / Е. С. Вентцель. - М. : Наука, 2016. – 345 с. Курицкий, Б. Я. Поиск оптимальных решений средствами EXCEL 7.0 / Б. Я.Курицкий. -СПб. : BHV-Санкт-Петербург, 2016. – 389 с. Математическое моделирование экономических процессов на железнодорожном транспорте / под ред. А. Б. Каплана. - М. : Транспорт, 2017. – 390 с. Математическое программирование / под ред. Н. Ш. Кремера. - М. : Финстатинформ, 2018. – 389 с. Мур, Дж. Экономическое моделирование в MicrosoftExcel / Дж. Мур, Л. Р. Уэдерфорд. – М: АСТ, 2019. – 239 с. |