Скачать .docx |
Реферат: Экономико математические методы и модели 3
МИНИСТЕРСТВО ОБРАЗОВАНИЯ РЕСПУБЛИКИ БЕЛАРУСЬ
УЧРЕЖДЕНИЕ ОБРАЗОВАНИЯ
«ПОЛОЦКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ»
Кафедра: Высшей математики
Контрольная работа
по дисциплине «Экономико-математические методы и модели»
Вариант - 12
Студентки Финансово-экономического факультета
Заочного отделения
Группы У06-ЭПз-1
Галай Натальи Михайловны
Преподаватель:
Сороговец И.Б.
Новополоцк, 2008 г.
Содержание
По условиям контракта торгово-посредническая фирма должна поставить каждому из двух покупателей Bj (j = 1, 2) два вида товаров Tk (k = 1, 2) в количестве bj k по цене рj k за единицу товара. Эти товары можно закупить у трех производителей Ai (i = 1, 2, 3) по цене si k за единицу товара. Известны: количества ai k товара Tk , имеющегося у производителя Ai , а также стоимости ci j k перевозки единицы товара Tk от производителя Ai к покупателю Bj .
ТРЕБУЕТСЯ:
1. Построить математическую модель поставленной задачи, максимизирующую прибыль фирмы от реализации всех сделок в виде задачи линейного программирования.
2. Методом потенциалов найти оптимальный план закупок, перевозок и поставок по каждому товару от каждого производителя к каждому покупателю, а также сумму прибыли от реализации этого плана.
Исходные данные:
a11 | a12 | a21 | a22 | a31 | a32 |
400 | 410 | 480 | 550 | 420 | 480 |
s11 | s12 | s21 | s22 | s31 | s32 |
2 | 3 | 5 | 5 | 2 | 3 |
b11 | b12 | b21 | b22 |
480 | 130 | 270 | 320 |
c111 | c112 | c121 | c122 | c211 | c212 | c221 | c222 | c311 | c312 | c321 | c322 |
2 | 2 | 2 | 2 | 3 | 2 | 1 | 2 | 3 | 2 | 2 | 1 |
p11 | p12 | p21 | p22 |
16 | 14 | 15 | 15 |
РЕШЕНИЕ:
1) Для составления математической модели введем неизвестные – количество товара , покупаемое у производителя для перевозки потребителю . Индексы: i = 1, 2, 3 – номер производителя продукции;
j = 1, 2 – номер потребителя продукции; k = 1, 2 – номер товара. Найдем тарифы , т.е. прибыли на одну единицу товара , покупаемое у производителя для продажи потребителю . Эти прибыли состоят из цены продажи 1 единицы товара за вычетом цены покупки и стоимости перевозки, т.е. .
f111 = p11 – s11 – c111 = 16 – 2 – 2 = 12;
f112 = p12 – s12 – c112 = 14 – 3 – 2 = 9;
f121 = p21 – s11 – c121 = 15 – 2 – 2 = 11;
f122 = p22 – s12 – c122 = 15 – 3 – 2 = 10;
f211 = p11 – s21 – c211 = 16 – 5 – 3 = 8;
f212 = p12 – s22 – c212 = 14 – 5 – 2 = 7;
f221 = p21 – s21 – c221 = 15 – 5 – 1 = 9;
f222 = p22 – s22 – c222 = 15 – 5 – 2 = 8;
f311 = p11 – s31 – c311 = 16 – 2 – 3 = 11;
f312 = p12 – s32 – c312 = 14 – 3 – 2 = 9;
f321 = p21 – s31 – c321 = 15 – 2 – 2 = 11;
f322 = p22 – s32 – c322 = 15 – 3 – 1 = 11.
Значение полученных коэффициентов приведены в таблице 1.1:
fijk | f111 | f112 | f121 | f122 | f211 | f212 | f221 | f222 | f311 | f312 | f321 | f322 |
тариф | 12 | 9 | 11 | 10 | 8 | 7 | 9 | 8 | 11 | 9 | 11 | 11 |
Прибыль фирмы представляется выражением , где сумма берется по всем возможным значениям индексов i , j , k . По условию, выражение F следует максимизировать, т.е. F является целевой функцией поставленной задачи. Так как операции над товарами и можно производить по отдельности и выражение F представляется в виде суммы двух слагаемых , сгруппированных по товарам , , то поставленная задача сводится к решению двух оптимизационных задач. Ограничения для неизвестных диктуются наличием соответствующих товаров у производителей и потребностью в них покупателей. В результате приходим к двум задачам линейного программирования, которые относятся к задачам транспортного типа:
Задача 1 (по товару )
maxF1 = 12 * X111 + 11 * X121 + 8 * X211 + 9 * X221 + 11 * X311 + 11 * X321
X111 + X211 + X311 ≤ b11 ≤ 480 X121 + X221 + X321 ≤ b21 ≤ 270
X111 + X121 ≤ a11 ≤ 400 X211 + X221 ≤ a21 ≤ 480 X311 + X321 ≤ a31 ≤ 420
Xij 1 ≥ 0
Задача 2 (по товару )
maxF2 = 9 * X112 + 10 * X122 + 7 * X212 + 8 * X222 + 9 * X312 + 11 * X322
X112 + X212 + X312 ≤ b12 ≤ 130 X122 + X222 + X322 ≤ b22 ≤ 320
X112 + X122 ≤ a12 ≤ 410 X212 + X222 ≤ a22 ≤ 550 X312 + X322 ≤ a32 ≤ 480
Xij 2 ≥ 0
Как видно, решение поставленной задачи сводится к решению двух задач транспортного типа.
2) Для решения задач 1, 2 методом потенциалов, сопоставим суммарное наличие каждого товара у производителей и суммарные потребности покупателей.
= 400 + 480 + 420 = 1300, = 480 + 270 = 750;
1300 – 750 = 550
Наличие товара Т1 превышает потребности покупателей. Вводим фиктивного покупателя В3 с потребностью b31 = 550.
= 410 + 550 + 480 = 1440, = 130 + 320 = 990;
1440 – 450 = 990
Наличие товара Т2 превышает потребности покупателей. Вводим фиктивного покупателя В3 с потребностью b32 = 990.
Получаем закрытые модели двух транспортных задач. Для их решения составляем две таблицы. В верхних правых углах клеток выписаны тарифы и . Для фиктивных производителей и покупателей тарифы равны нулю. Последние строки и столбцы таблиц служат для записи потенциалов.
Таблица 1.2 (к задаче 1)
Производители | Покупатели | |||||
B1 | B2 | B3 | ai1 | ui | ||
A1 | 12 400 |
11 12 |
0 1 |
400 |
1 |
|
A2 | 8 11 |
9 11 |
0 480 |
480 |
0 |
|
A3 | 11 80 |
11 270 |
0 70 |
420 |
0 |
|
bj1 | 480 | 270 | 550 | 1300 | – | |
vj | 11 | 11 | 0 | – | – |
Таблица 1.3 (к задаче 2)
Производители | Покупатели | |||||
B1 | B2 | B3 | ai1 | ui | ||
A1 | 9 130 |
10 11 |
0 280 |
410 |
0 |
|
A2 | 7 9 |
8 11 |
0 550 |
550 |
0 |
|
A3 | 9 9 |
11 320 |
0 160 |
480 |
0 |
|
bj1 | 130 | 320 | 990 | 1440 | – | |
vj | 9 | 11 | 0 | – | – |
Начальные планы распределения товаров определены по методу максимальной прибыли, т.е. в первую очередь заполнялись по максимуму клетки с наибольшими тарифами. Более конкретно, просматривая таблицу 1.2, замечаем, что максимальный тариф 12 стоит в клетке (1,1). В эту клетку ставим число 400. При этом запасы производителя А1 исчерпан. Далее, в клетку (3,1) ставим 80, а в клетку (3,2) ставим 270. Из запасов производителя А3 осталось 70, так как 420-80-270=70, ставим их в клетку (3,3). Потребность покупателей В1 и В2 в товарах исчерпаны, следовательно, оставшиеся 480 товаров производителя А2 ставим в клетку (2,3). При этом товар производителей полностью распределён.
Полученный начальный план проверим на оптимальность. План невырожденный, так как число занятых клеток (3+3-1=5) равно m + n – 1 (m и n – число строк и столбцов распределительной матрицы). Обозначим через и потенциалы строк и столбцов. Для их нахождения отметим, что в занятых клетках сумма потенциалов строки и столбца должна равняться тарифу клетки. Получаем в данном случае 5 уравнений с 6-ю неизвестными:
v1 + u1 = 12; v2 + u3 = 11; v3 + u3 = 0.
v1 + u3 = 11; v3 + u2 = 0;
Полагая, что v3 = 0, последовательно получаем: u2 = 0, u3 = 0, v2 = 11, v1 = 11, u1 = 1.
Так как задача решается на максимум, то для оптимальности плана распределения, сумма потенциалов в незанятых клетках должна быть не меньше тарифов этих клеток. В нижних левых углах незанятых клеток выписаны суммы потенциалов. Все они превосходят соответствующие тарифы, т.е. начальный план закрепления покупателей за производителями по товару оптимален.
Аналогично, таблица 1.3 заполнена в следующей последовательности:
(3,2) – 320, (1,1) – 130, (1,3) – 280, (3,3) – 160, (2,3) – 550. Полученный план невырожденный, так как содержит 3 + 3 – 1 = 5 занятых клеток. Проверим его на оптимальность. Выпишем систему уравнений для нахождения потенциалов:
v1 + u1 = 9; v3 + u1 = 0; v3 + u3 = 0.
V2 + u3 = 11; v3 + u2 = 0;
Полагая, что u3 = 0, последовательно получаем: v3 = 0, u2 = 0, u1 = 0, v2 = 11, v1 = 9.
План распределения товара T 2 , заданный таблицей 2, оптимален.
Сумма прибыли = (12*400 + 11*80 + 11*270) + (9*130 + 11*320) = = 8650 + 4690 = 13340.
ОТВЕТ:
X111 = 400, X311 = 80, X321 = 270, X112 = 130, X322 = 320. Остальные = 0. Максимальная прибыль равна 13340.
С помощью алгоритма венгерского метода найти план закрепления работ за исполнителями, максимизирующий прибыль, связанную с выпуском всех пяти видов продукции. Матрица эффективности AN =, где – прибыль, получаемая при выполнении j -й работы i -м исполнителем, N - номер варианта, имеет вид:
40 | 28 | 44 | 38 | 46 | ||
36 | 52 | 51 | 43 | 30 | ||
A12= | 40 | 29 | 48 | 45 | 34 | , |
56 | 54 | 53 | 46 | 49 | ||
51 | 41 | 50 | 55 | 41 |
РЕШЕНИЕ:
I этап: приведение матрицы А12.
Алгоритм венгерского метода предназначен для решения задачи о назначениях по критерию минимизации суммарных затрат (задача на минимум). При решении задачи на максимум (так как – прибыль), ее следует свести к задаче на минимум. Для этого в каждом столбце матрицы определяем максимальный элемент и из него вычитаем все элементы столбца.
40 | 28 | 44 | 38 | 46 | 16 | 26 | 9 | 17 | 3 | |||
36 | 52 | 51 | 43 | 30 | 20 | 2 | 2 | 12 | 19 | |||
A12= | 40 | 29 | 48 | 45 | 34 | → A121= | 16 | 25 | 5 | 10 | 15 | → |
56 | 54 | 53 | 46 | 49 | 0 | 0 | 0 | 9 | 0 | |||
51 | 41 | 50 | 55 | 41 | 5 | 13 | 3 | 0 | 8 | |||
56 | 54 | 53 | 55 | 49 |
Так как в строках 1, 2, 3 нулей не оказалось, то вычитаем из элементов этих строк минимального из них, то есть вычитаем из строки 1 число 3, из строки 2 число 2, из строки 3 число 5. Получаем нули в этих строках.
13 | 23 | 6 | 14 | 0 | ||
18 | 0 | 0 | 10 | 17 | ||
→ A122= | 11 | 20 | 0 | 5 | 10 | → |
0 | 0 | 0 | 9 | 0 | ||
5 | 13 | 3 | 0 | 8 |
II этап: поиск назначения.
Выбираем один из нулей, помечаем его, например, точкой или звездочкой или обводим его другим цветом (в дальнейшем, звездочкой), а остальные нули строки и столбца, в которых стоит выбранный помеченный нуль, перечеркиваем. Далее переходим к следующему нулю. И так до тех пор, пока каждый нуль будет либо помечен, либо перечеркнут.
13 | 23 | 6 | 14 | 0* | |
18 | 0* | Ø | 10 | 17 | |
→ A122= | 11 | 20 | 0* | 5 | 10 |
0* | Ø | Ø | 9 | Ø | |
5 | 13 | 3 | 0* | 8 |
Помеченные нули составили полное назначение (количество помеченных нулей равно 5).
Следовательно, 1 исполнитель назначается на 5-ю работу, 2→2, 3→3, 4→1, 5→4.
Задача 3-1
Предприятие включает в себя три цеха по производству различной продукции и использует при этом четыре вида первичных ресурсов. Продукция, выпускаемая каждым цехом, частично отгружается за пределы предприятия (для удовлетворения конечного спроса), а частично распределяется внутри предприятия между цехами в качестве вторичных ресурсов. Баланс предприятия в натуральном выражении за прошедший год приведен в следующих двух таблицах 3.1.1 и 3.1.2:
Таблица 3.1.1
Производство товаров |
Внутреннее потребление | Конечный спрос |
||
цех 1 | цех 2 | цех 3 | ||
цех 1 | 240 | 72 | 140 | 348 |
цех 2 | 80 | 264 | 180 | 76 |
цех 3 | 0 | 120 | 400 | 480 |
Таблица 3.1.2
Первичные ресурсы |
Расходы ресурса за год | ||
цех 1 | цех 2 | цех 3 | |
А | 180 | 30 | 50 |
Б | 1200 | 1500 | 0 |
В | 400 | 1200 | 300 |
Г | 160 | 600 | 1000 |
ТРЕБУЕТСЯ:
1) найти матрицы коэффициентов прямых товаро-затрат и ресурсо-затрат на основании данных за предыдущий год;
2) найти план полных выпусков продукции каждого цеха на следующий год, обеспечивающих выполнение госзаказа по отгрузке продукции в объемах c1 =360, c2 =90, c3 =450 соответственно;
3)определить необходимый запас первичных ресурсов каждого вида.
РЕШЕНИЕ:
Если обозначить через полные выпуски продукции каждым цехом, то можно составить следующие соотношения
где – непосредственный натуральный расход продукции i -го цеха для обеспечения выпуска всей продукции j -го цеха. Числа называются коэффициентами прямых товаро-затрат. Их можно определить по статистическим данным за предыдущий год, т.е. Смысл коэффициента – количество продукции i -го цеха, используемое для производства 1 единицы продукции j -го цеха.
Аналогично, расход k -го ресурса j -м цехом представим в виде
Тогда коэффициенты называются коэффициен-тами прямых ресурсо-затрат. Они определяют количество k -го ресурса, необходимое для производства единицы продукции j -го цеха и находятся по результатам статистических данных за предыдущий год.
1. Для определения коэффициентов найдем полные выпуски продукции каждым цехом за предыдущий год хj :
x1 = 240 + 72 + 140 + 348 = 800;
x2 = 80 + 264 + 180 + 76 = 600;
x3 = 0 + 120 + 400 + 480 = 1000;
Тогда матрицы A и B коэффициентов и принимают вид:
240 | 72 | 140 | 0.30 | 0.12 | 0.14 |
800 | 600 | 1000 | |||
0.10 | 0.44 | 0.18 | , | ||
A= | 80 | 264 | 180 | = | |
800 | 600 | 1000 | 0.00 | 0.20 | 0.40 |
0 | 120 | 400 | |||
800 | 600 | 1000 |
180 | 30 | 50 | 0.23 | 0.05 | 0.05 |
800 | 600 | 1000 | |||
1.50 | 2.50 | 0.00 | , | ||
B= | 1200 | 1500 | 0 | = | |
800 | 600 | 1000 | 0.50 | 2.00 | 0.30 |
400 | 1200 | 300 | 0.20 | 1.00 | 1.00 |
800 | 600 | 1000 | |||
160 | 600 | 1000 | |||
800 | 600 | 1000 |
2. Заменяя выражения найденными коэффициентами получаем систему уравнений для определения искомых полных выпусков продукции:
x1 = 0.30*x1 + 0.12*x2 + 0.14*x3 + 360, x2 = 0.10*x1 + 0.44*x2 + 0.18*x3 + 90, x3 = 0.00*x1 + 0.20*x2 + 0.40*x3 + 450. |
Эту систему можно записать в матричной форме где E – единичная матрица, X – матрица-столбец из неизвестных, C – матрица-столбец из чисел c1 =360, c2 =90, c3 =450. Решая полученное матричное уравнение, найдем полные выпуски продукции. Его решение имеет вид: Строим обратную матрицу Для этого найдем алгебраические дополнения и определитель для матрицы Имеем:
0,70 | -0,12 | -0,14 | ||
Е–A= | -0,10 | 0,56 | -0,18 | , |
0,00 | -0,20 | 0,60 |
0,56 | -0,18 | ||||
A11= | = | 0,56*0,60 – (-0,18)*(-0,20) | = | 0,30 | , |
-0,20 | 0,60 |
Аналогично:
-0,10 | -0,18 | -0,10 | 0,56 | ||||
A12= – | = | 0,06 | , | A13= | = | 0,02 | , |
0,00 | 0,60 | 0,00 | -0,20 |
-0,12 | -0,14 | 0,70 | -0,14 | ||||
A21= – | = | 0,10 | , | A22= | = | 0,42 | , |
-0,20 | 0,60 | 0,00 | 0,60 |
0,70 | -0,14 | -0,12 | -0,14 | ||||
A23= – | = | 0,14 | , | A31= | = | 0,10 | , |
-0,10 | -0,18 | 0,56 | -0,18 |
0,70 | -0,14 | 0,70 | -0,12 | ||||
A32= – | = | 0,14 | , | A33= | = | 0,38 | , |
-0,10 | -0,18 | -0,10 | 0,56 |
При этом ∆ = 0,70*0,30 – 0,12*0,06 – 0,14*0,02 = 0,20,
1,5 | 0,5 | 0,5 | ||
= | 0,3 | 2,1 | 0,7 | , |
0,1 | 0,7 | 1,9 |
Умножая матрицу на C , найдем искомые полные выпуски продукции:
х1 | 1,5 | 0,5 | 0,5 | 360 | 1,5*360 | + | 0,5*90 | + | 0,5*450 | 810 | |||||
х2 | = | 0,3 | 2,1 | 0,7 | * | 90 | = | 0,3*360 | + | 2,1*90 | + | 0,7*450 | = | 612 | , |
х3 | 0,1 | 0,7 | 1,9 | 450 | 0,1*360 | + | 0,7*90 | + | 1,9*450 | 954 |
То есть, х1 = 810, х2 = 612, х3 = 954.
3. При определении запаса k -го вида ресурсов, необходимого для производства найденных полных выпусков продукции, достаточно умножить матрицу ресурсо-затрат B на матрицу-столбец из полных выпусков продукции:
b1 | 0,23 | 0,05 | 0,05 | 0,23*810 | + | 0,05*612 | + | 0,05*954 | 264,6 |
810 | |||||||||
b2 | 1,50 | 2,50 | 0,00 | 1,50*810 | + | 2,50*612 | + | 0,00*954 | 2745,0 |
= | * | 612 | = | = | , | ||||
b3 | 0,50 | 2,00 | 0,30 | 0,50*810 | + | 2,00*612 | + | 0,30*954 | 1915,2 |
954 | |||||||||
b4 | 0,20 | 1,00 | 1,00 | 0,20*810 | + | 1,00*612 | + | 1,00*954 | 1728,0 |
То есть запас ресурса следует иметь в количестве 264,6 ед., ресурса – в количестве 2745 ед., ресурса – в количестве 1915,2 ед., ресурса – в количестве 1728 ед.
Урожайность пшеницы зависит от количества внесенных удобрений и погодных условий. Фермер может вносить на 1 гектар , или центнеров удобрений. Погодные условия характеризуются тремя состояниями: , и . Урожайность пшеницы с одного гектара составляет центнеров при внесении центнеров удобрений и состоянии погоды . Рыночная цена на зерно составляет ден. ед., если было внесено ц/га удобрений. Стоимость одного центнера удобрений составляет S ден. ед.
Требуется определить, какое количество удобрений следует вносить в почву, чтобы получить как можно большую прибыль, если: а) известны вероятности состояний природы ; б) о вероятностях состояний природы ничего определенного сказать нельзя.
Указание. Составить платежную матрицу, рассчитав значении прибыли по формуле: , .
Исходные данные:
а1 | а2 | а3 | с1 | с2 | с3 | b11 | b12 | b13 | b21 | b22 | b23 | b31 | b32 | b33 | S | p1 | p2 | p3 | λ |
2 | 4 | 6 | 9 | 5 | 3 | 5 | 9 | 6 | 10 | 12 | 9 | 13 | 15 | 11 | 4 | 0,3 | 0,4 | 0,3 | 0,8 |
РЕШЕНИЕ:
Одним из участников рассматриваемой ситуации является фермер, который должен вносить удобрения в почву для получения хорошего урожая пшеницы. Если описанной ситуации придать игровую схему, то фермер выступит в ней в качестве сознательного игрока А, заинтересованного в максимизации прибыли с 1 гектара земли. Вторым участником является в буквальном смысле природа (игрок П), то есть внешние природные условия.
Так как фермер на 1 гектар земли может вносить разное количество центнеров удобрений, то чистыми стратегиями игрока А будут следующие стратегии:
– А1 : вносить 2 ц. удобрений на 1 гектар земли;
– А2 : вносить 4 ц. удобрений на 1 гектар земли;
– А3 : вносить 6 ц. удобрений на 1 гектар земли.
Природа может реализовать одно из трех состояний: П1 , П2 , П3 .
Таким образом, платежная матрица игры будет иметь размер 3х3.
Вычисляем значении прибыли по формуле:, .
h11 = 9*5 – 4*2 = 37; h23 = 5*9 – 4*4 = 29;
h12 = 9*9 – 4*2 = 73; h31 = 3*13 – 4*6 = 15;
h13 = 9*6 – 4*2 = 46; h32 = 3*15 – 4*6 = 21;
h21 = 5*10 – 4*4 = 34; h33 = 3*11 – 4*6 = 9;
h22 = 5*12 – 4*4 = 44;
Итак, платежная матрица принимает вид (таблица 4.1)
37 | 73 | 46 | |
34 | 44 | 29 | |
15 | 21 | 9 |
В платежной матрице нет доминируемых стратегий игрока А , поэтому матрица не требует упрощений.
а) для определения оптимальной стратегии игрока А по критерию Байеса вычислим среднее значение (математическое ожидание) выигрыша при использовании каждой из возможных стратегий по формуле: . Получаем:
= 37*0,3 + 73*0,4 + 46*0,3 = 54,1;
= 34*0,3 + 44*0,4 + 29*0,3 = 36,5;
= 15*0,3 + 21*0,4 + 9*0,3 = 15,6.
Оптимальной по критерию Байеса является стратегия , так как именно ей соответствует наибольшее из чисел :
max | { | 54.1 | ; | 73 | ; | 46 | } | = | 73 ; |
Таким образом, располагая информацией о возможных состояниях природы, наиболее выгодным для фермера будет использование стратегии А1 – вносить 2 ц. удобрений на 1 гектар земли. Среднее значение ожидаемой прибыли в этом случае составит 54,1 ден. ед.
б) для определения оптимальной стратегии игрока А с использованием максимаксного критерия , применим формулу: .
Получаем:
m1 = {37; 73; 46} = 73;
m2 = {34; 44; 29} = 44;
m3 = {15; 21; 9} = 21;
Оптимальной по максимаксному критерию является стратегия , так как именно ей соответствует наибольшее из чисел :
max | { | 73 | ; | 44 | ; | 21 | } | = | 73 ; |
Таким образом, в расчете на самое благоприятное стечение обстоятельств, наиболее выгодным для домовладельца будет использование стратегии – вносить 2 ц. удобрений на 1 гектар земли. Прибыль, потраченная при этом от продажи зерна, составит 73 ден. ед.
Определим оптимальную стратегию игрока А по критерию Вальда :
w1 = min {37; 73; 46} = 37;
w2 = min {34; 44; 29} = 29;
w3 = min {15; 21; 9} = 9.
max | { | 37 | ; | 29 | ; | 9 | } | = | 37 ; |
Следовательно, оптимальной по критерию Вальда является стратегия – вносить 2 ц. удобрений на 1 гектар земли. При этом минимальная прибыль составит 37 ден. ед.
Для определения оптимальной стратегии игрока А с использованием критерия Сэвиджа составим матрицу рисков. В каждом столбце платежной матрицы определим максимальный элемент и вычтем из него все элементы данного столбца. В первом столбце максимальным является элемент h11 = 37, во втором – h12 = 73, в третьем – h13 = 46.
Матрица рисков представлена в таблице 4.2.
Таблица 4.2
0 | 0 | 0 | |
3 | 29 | 17 | |
22 | 52 | 37 |
Определим максимальный риск при использовании каждой стратегии.
Получаем:
r1 = max {0; 0; 0} = 0,
r2 = max {3; 29; 17} = 29,
r3 = max {22; 52; 37} = 52.
min | { | 0 | ; | 29 | ; | 52 | } | = | 0 ; |
Таким образом, оптимальной по Сэвиджу является стратегия – вносить 2 ц. удобрений на 1 гектар земли.
Для определения оптимальной стратегии по критерию Гурвица найдем показатель критерия по формуле , .
Получаем:
γ1 = 0,8*37 + (1 – 0,8)*73 = 44,2;
γ2 = 0,8*29 + (1 – 0,8)*44 = 32,0;
γ3 = 0,8*9 + (1 – 0,8)*21 = 11,4.
max | { | 44,2 | ; | 32,0 | ; | 11,4 | } | = | 44,2 ; |
Следовательно, оптимальной является стратегия – вносить 2 ц. удобрений на 1 гектар земли.
Для определения оптимальной стратегии игрока А с использованием критерия Лапласа определим средние арифметические значения «выигрыша» домовладельца по формуле , .
Получаем:
= (37 + 73 + 46)/3 = 156/3 = 52;
= (34 + 44 + 29)/3 = 107/3 = ;
= (15 + 21 + 9)/3 = 45/3 = 15.
Оптимальной по критерию Лапласа является стратегия – вносить 2 ц. удобрений на 1 гектар земли, так как ей соответствует наибольшее из чисел :
max | { | 52 | ; | ; | 15 | } | = | 52 ; |
Таким образом, если все состояния природы представляются равновозможными, то для обеспечения средней прибыли в размере 52 ден. ед. фермеру следует придерживаться стратегии – вносить 2 ц. удобрений на 1 гектар земли.
Для наглядности все результаты вычислений приведем в сводной таблице 4.3.
Таблица 4.3
37 | 73 | 46 | 54,1 | 73 | 37 | 0 | 44,2 | 52 | |
34 | 44 | 29 | 36,5 | 44 | 29 | 29 | 32,0 | ||
15 | 21 | 9 | 15,6 | 21 | 9 | 52 | 11,4 | 15 |
Вывод: Проведенное по совокупности статистических критериев исследование возможных вариантов внесения удобрений на 1 гектар земли позволяет сделать следующее заключение:
а) при наличии достоверной информации о состоянии природы фермеру следует вносить 2 ц. удобрений на 1 гектар земли (стратегия ). При этом ожидаемая прибыль составит 54,1 ден. ед.
б) при отсутствии информации о состоянии природы фермеру также следует вносить 2 ц. удобрений на 1 гектар земли (стратегия ). Такое решение продиктовано на основании всех критериев, не использующих вероятностный подход. Заметим, что максимальная прибыль при выборе данной стратегии составит 37 ден. ед. при самых неблагоприятных погодных условиях.
Задача 5-2
Для реконструкции и модернизации производства выделены денежные средства в объеме 100 тыс. ден. ед., которые следует распределить между четырьмя цехами. По каждому из цехов известен возможный прирост выпуска продукции в зависимости от выделенной ему суммы x ≤ 100.000 (возможные значения x и приведены в таблице). Необходимо так распределить средства, чтобы максимально увеличить выпуск продукции производства в целом.
ЦЕХ № 1 | ЦЕХ № 2 | ||||||||||
x | 20 | 40 | 60 | 80 | 100 | x | 20 | 40 | 60 | 80 | 100 |
9 | 17 | 29 | 38 | 47 | 11 | 34 | 46 | 53 | 75 |
ЦЕХ № 3 | ЦЕХ № 4 | ||||||||||
x | 20 | 40 | 60 | 80 | 100 | x | 20 | 40 | 60 | 80 | 100 |
13 | 28 | 37 | 49 | 61 | 12 | 35 | 40 | 54 | 73 |
ТРЕБУЕТСЯ:
1. Основываясь на принципах динамического программирования, построить математическую модель поставленной задачи в виде функциональных уравнений Беллмана (числовые данные взять из таблиц).
2. Найти оптимальное распределение средств, обеспечивающее максимальный прирост выпуска продукции.
PЕШЕHИЕ.
1. Системой S в данном случае является предприятие из 4-х цехов, в которое вложена сумма 100.000 ед. Состояния и управления системы S однозначно взаимосвязаны – это способы распределения суммы между цехами. Для осуществления инвариантного погружения задачи будем считать, что вместо суммы 100.000 ед. вкладывается сумма y : 0≤у ≤100.000. Состояния системы искусственно разобьем на этапы: начальный (нулевой), первый, второй и третий этапы соответственно означают, что сумма y распределяется между четырьмя цехами, тремя цехами, двумя цехами и вся сумма y выделяется одному цеху. Нумерацию этапов удобнее проводить в обратном порядке: третий этап – m = 1,второй этап – m = 2, первый этап – m = 3, нулевой этап – m = 4. Тогда функция Беллмана, имеющая смысл максимальной прибыли при распределении суммы y между m цехами, запишется в виде:
Если при m = 1 … k –1 функция B (y , m ) уже построена, то функциональное уравнение Беллмана для данной задачи принимает вид:
Пpи m = 1 дополнительно имеем:
2. При m = 1 функция Беллмана уже построена, т.е.
y | 20 | 40 | 60 | 80 | 100 |
B(y, 1) | 9 | 17 | 29 | 38 | 47 |
При m = 2 уравнение из функционального уравнения Беллмана имеет вид:
Так как функции и заданы таблично, то для определения максимума функции при каждом y составляем таблицу значений этой функции:
x y |
0 | 20 | 40 | 60 | 80 | 100 | B(y, 2) | |
20 | 0 + 9 | 11 + 0 | 11 | 20 | ||||
40 | 0 + 17 | 11 + 9 | 34 + 0 | 34 | 40 | |||
60 | 0 + 29 | 11 + 17 | 34 + 9 | 46 + 0 | 46 | 60 | ||
80 | 0 + 38 | 11 + 29 | 34 + 17 | 46 + 9 | 53 + 0 | 55 | 60 | |
100 | 0 + 47 | 11 + 38 | 34 + 29 | 46 + 17 | 53 + 9 | 75 + 0 | 75 | 100 |
Подчеркнутые значения являются максимальными в строке, т.е. являются значениями функции Беллмана B (y , 2). Они выписаны в предпоследнем столбце. В последний столбец выписаны значения x , при которых достигается максимум функции Эти значения обозначены и их можно считать управлениями. Смысл – средства, выделяемые второму цеху, при оптимальном распределении суммы y между двумя цехами.
Аналогично, при m = 3
Составляем таблицу значений функции
x y |
0 | 20 | 40 | 60 | 80 | 100 | B(y, 3) | |
20 | 0 + 11 | 13 + 0 | 13 | 20 | ||||
40 | 0 + 34 | 13 + 11 | 28 + 0 | 34 | 0 | |||
60 | 0 + 46 | 13 + 34 | 28 + 11 | 37 + 0 | 47 | 20 | ||
80 | 0 + 55 | 13 + 46 | 28 + 34 | 37 + 11 | 49 + 0 | 62 | 40 | |
100 | 0 + 75 | 13 + 55 | 28 + 46 | 37 + 34 | 49 + 11 | 61 + 0 | 75 | 0 |
Как и выше – сумма средств, выделяемых третьему цеху, при оптимальном распределении суммы y между тремя цехами.
При m = 4
Составляем таблицу значений функции
x y |
0 | 20 | 40 | 60 | 80 | 100 | B(y, 4) | |
20 | 0 + 13 | 12 + 0 | 13 | 0 | ||||
40 | 0 + 34 | 12 + 13 | 35 + 0 | 35 | 40 | |||
60 | 0 + 47 | 12 + 34 | 35 + 13 | 40 + 0 | 48 | 40 | ||
80 | 0 + 62 | 12 + 47 | 35 + 34 | 40 + 13 | 54 + 0 | 69 | 40 | |
100 | 0 + 75 | 12 + 62 | 35 + 47 | 40 + 34 | 54 + 13 | 73 + 0 | 82 | 40 |
В двух последних столбцах этой таблицы получены значения функции Беллмана B (y , 4) и соответствующие им управления – т.е. количества средств, выделяемых четвертому цеху при распределении суммы y между четырьмя цехами. После этого составляем сводную таблицу значений функции Беллмана и соответствующих ей управлений:
y |
||||||||
0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
20 | 20 | 9 | 20 | 11 | 20 | 13 | 0 | 13 |
40 | 40 | 17 | 34 | 0 | 34 | 40 | 35 | |
60 | 60 | 29 | 60 | 46 | 47 | 40 | 48 | |
80 | 80 | 38 | 60 | 55 | 40 | 62 | 40 | 69 |
100 | 100 | 47 | 100 | 75 | 0 | 75 | 82 |
С помощью таблицы функции Беллмана для данной задачи можно произвести распределение любой суммы у от 0 до 100 между k цехами 1 ≤ k ≤ 4. В клетке стоит максимальная прибыль от этого распределения, а в клетке стоит сумма, выделяемая k -му цеху. Распределим сумму 100 между 4-мя цехами. По клетке максимально возможная прибыль равна 82. 4-му цеху следует выделить 40 тыс. $. На первые три цеха остается 60 тыс. $. По клетке 3-му цеху выделяется 20 тыс. $. На первые два цеха остается 40 тыс. $. По клетке 2-му цеху выделяется 40 тыс. $. Тогда 1-му цеху средства не выделяются.