Контрольная работа по "Экономико- математические методы и ПМ"
Контрольная работа, 21 Января 2013, автор: пользователь скрыл имя
Описание работы
Задача 1.
Решить графическим методом типовую задачу оптимизации.
1.4. На имеющихся у фермера 400 гектарах земли он планирует посеять кукурузу и сою. Сев и уборка кукурузы требует на каждый гектар 200 ден. ед. затрат, а сои – 100 ден. ед. На покрытие расходов, связанных с севом и уборкой, фермер получил ссуду в 60 тыс. ден. ед.. Каждый гектар, засеянный кукурузой, принесет 30 центнеров, а каждый гектар, засеянный соей – 60 центнеров. Фермер заключил договор на продажу, по которому каждый центнер кукурузы принесет ему 3 ден. ед., а каждый центнер сои – 6 ден. ед. Однако, согласно этому договору, фермер обязан хранить убранное зерно в течение нескольких месяцев на складе, максимальная вместимость которого равна 21 тыс. центнеров.
Фермеру хотелось бы знать, сколько гектар нужно засеять каждой из этих культур, чтобы получить максимальную прибыль.
Построить экономико-математическую модель задачи, дать необходимые комментарии к ее элементам и получить решение графическим методом. Что произойдет, если решать задачу на минимум и почему?
Задача 2
2.4 Использовать аппарат теории двойственности для экономико-математического анализа оптимального плана задачи линейного программирования.
Для изготовления трех видов продукции используют три вида сырья. Запасы сырья, нормы его расхода и цены реализации единицы каждого вида продукции приведены в таблице.
Файлы: 1 файл
Контр ЭММ.doc
— 793.00 Кб (Скачать файл)- Курсор в поле «Изменяя ячейки».
- Ввести адреса $C$3:$E$3.
7) Ввод ограничений.
- Курсор в поле «Добавить». Появится диалоговое окно Добавление ограничения (рис. 4).
Рис. 4. Ввод правых и левых частей ограничений.
- В поле «Ссылка на ячейку» ввести адрес $D$7.
- Ввести знак ограничения <=.
- Курсор в правое окно.
- Ввести адрес $F$7.
- Добавить. На экране опять диалоговое окно Добавление ограничения.
- Ввести второе ограничение.
- Добавить. На экране опять диалоговое окно Добавление ограничения.
- Ввести третье ограничение ОК.
На экране появится диалоговое окно Поиск решения с введенными условиями (рис. 5).
8) Ввод параметров для решения ЗЛП (рис. 6.).
- Открыть окно Параметры поиска решения.
- Установить флажок Линейная модель, что обеспечивает применение симплекс-метода.
- Установить флажок Неотрицательные значения.
- ОК. (На экране диалоговое окно Поиска решения).
- Выполнить. (На экране диалоговое окно Результаты поиска решения – рис. 7).
Рис. 5. Введены все условия для решения задачи.
Рис. 6. Ввод параметров.
Рис. 7. Решение.
Полученное решение означает, что максимум функции равен 1420, при Х1=0, Х2=74 и Х3=32.
Проверим, как удовлетворяется система функциональных ограничений оптимальным планом X* = (x1 = 0, x2 = 74, х3 = 32):
4*0+2*74+1*32 = 180 = 180
3*0+1*74+2*32 = 138 < 210 (*)
1*0+2*74+3*32 = 244 = 244
Значение целевой функции на этом плане равно
f(X*) = 10*0+14*74+12*32 = 1420
2. Двойственная задача имеет вид:
min g(Y) = 180y1+210y2+244y3
4y1+3y2+y3 ≥ 10
2y1+y2+2у3 ≥ 14
y1+2y2+3y3 ≥ 12
y1,2,3 ≥ 0.
Для нахождения оценок y1, у2, у3 используем вторую теорему двойственности. Поскольку второе ограничение в (*) выполняется как строгое неравенство, то у2* = 0. Так как :х2 > 0 и х3 > 0, то:
2y1*+y2*+2y3* = 14
y1*+2y2*+3y3* = 12.
Итак, для получения двойственных оценок имеем систему линейных уравнений:
у2* = 0
2y1*+y2*+2y3* = 14
y1*+2y2*+3y3* = 12.
т.е. у1* = 4,5; y2* = 0; y3* = 2,5.
Вычислим значение целевой функции двойственной задачи:
g(Y*) = 180*4,5+210*0+244*2,5 = 1420, т.е. f(X*) = g(Y*) = 1420
По первой теореме двойственности мы можем утверждать, что действительно найдены оптимальные значения двойственных переменных.
3. Нулевое значение переменной х1 в оптимальном плане означает, что изготовление этого вида продукции не выгодно, т.к. цена реализации этого вида продукции низкая, а нормы расхода сырья на изготовление одного изделия этого вида высокие.
4. Двойственные оценки
отражают сравнительную
В примере недефицитным ресурсом является II тип сырья, поскольку у2= 0.
Острее ощущается дефицитность ресурса I тип сырья (у1 = 4,5) - он более дефицитен, чем ресурс III тип сырья (у3 = 2,5).
Предположим, что запасы сырья I вида увеличили на 4 ед., т. е. теперь они составляют 180 + 4 = 184 единиц, и запасы сырья III вида увеличили на 4 ед., т. е. теперь они составляют 244 + 4 = 248 единиц. Из теоремы об оценках Δf(X) = yi • Δbi, известно, что колебание величины bi приводит к увеличению или уменьшению f(X). Оно определяется величиной yi в случае, когда при изменении величин bi, значения переменных yi в оптимальном плане соответствующей двойственной задачи остаются неизменными. Поэтому необходимо найти такие интервалы изменения каждого из свободных членов системы ограничений исходной ЗЛП, в которых оптимальный план двойственной задачи не менялся бы.
Пределы уменьшения (нижняя граница) определяются по тем xk (k = 1,..., т), для которых соответствующие dki >0:
Δbi(-) = min{xk/dki} для dki > 0 (1)
Пределы увеличения (верхняя граница) определяются по тем xk, для которых dki < 0:
Δbi(+) =| max{xk/dki} | для dki < 0 (2)
Определим интервалы устойчивости двойственных оценок в примере. Матрица А имеет вид:
4 2 1
A = 3 1 2
1 2 3
После приведения задачи к канонической форме матрица А примет следующий вид:
4 2 1 1 0 0
A = 3 1 2 0 1 0
1 2 3 0 0 1
С ненулевыми значениями в оптимальный план вошли x2* = 74, x3* = 32 и х5* = 72 (из (*) 210-138=72), следовательно, матрица А* будет составлена из второго, третьего и пятого столбцов матрицы А:
2 1 0
A* = 1 2 1
2 3 0
Для вычисления интервалов устойчивости необходимо найти матрицу D = А*-1 (обратную матрицу матрицы А*):
3/4 0 -1/4
D = -1/2 0 1/2
1/4 1 -3/4
При вычислении интервалов устойчивости по формулам (1) и (2) примем
х2* = 74 = xk=1, х3* = 32 = xk=2 и х5* = 72 = хk=3.
Интервалы устойчивости первого ресурса — «запасы сырья I типа»:
Δb1(-) = min{x1/d11, x3/d31} = 74:3/4 = 296/3
Δb1(+) = | max{x2/d21} | = | 32:(-1/2) | = | -64 | = 64
b1 = {b1 - Δb1(-) ; b1 + Δb1(+)} = {180 – 296/3; 180+64} = {244/3; 244}.
При изменении запасов сырья I типа в пределах от 244/3 до 244 единиц двойственная оценка его не изменится.
Интервалы устойчивости второго ресурса — «запасы сырья II типа»: этот ресурс в оптимальном плане используется не полностью и поэтому не имеет верхней границы интервалов устойчивости; нижняя граница определяется следующим образом:
Δb2(-) = min{x3/d32} = 72:1 = 72
b2 = {b2 – Δb2(-) ; b2 + Δb2(+) } = {210-72;210} = {138; 210}.
Интервалы устойчивости третьего ресурса — «запасы сырья III типа»:
Δb3(-) = min{x2/d23} = 74:1/2 = 148
Δb3(+) = |mах{x1/d13, x3/d33}| = | mах{-74:1/4, -72:3/4}| = |-96| = 96
b3 = {b3 – Δb3(-) ; b3 + Δb3(+)} = {244-148; 244+96} = {96; 340}.
В нашем примере определим величину изменения объема прибыли от реализации продукции при увеличении запасов I и III типа сырья на 4 ед. каждого. Эти изменения находятся в интервалах устойчивости двойственных оценок (244/3<184<244; 96<248<340), поэтому можно воспользоваться теоремой об оценках:
Δf(X) = 4·4,5+4·2,5 = 28
Объем стоимости выпускаемой продукции увеличится на 28 единиц.
Такой же ответ мы получили бы, если бы решили симплексным методом задачу с новыми ограничениями по запасам сырья I и III типа. Новый оптимальный план: Х*=(0, 76, 32, 0, 70, 0)
f(X) = 0*10+76*14+32*12 = 1448,
Δf(X) = 1448-1420 = 28
Структурных сдвигов в программе не произошло, но значения переменных в плане изменились. Значение целевой функции при новых ограничениях увеличится на 28 единиц.
Двойственные оценки служат инструментом определения эффективности отдельных хозяйственных решений (технологических способов), с их помощью можно определять выгодность производства новых изделий, эффективность новых технологических способов:
если Δ j = ∑ aij yi* - cj ≤ 0 - выгодно,
если Δ j > 0 – невыгодно.
Определим целесообразность включения в план изделия «Г» ценой 13 ед., на изготовление которого расходуется 1, 3 и 2 ед. каждого вида сырья соответственно.
Δ 4 = 1*4,5 + 3*0 + 2*2,5– 13 = -3,5 < 0, т.е. включение в план изделия «Г» ценой 13 ед. выгодно.
Определим целесообразность включения в план изделия «Д» ценой 12 ед., на изготовление которого расходуется по две единицы каждого вида сырья.
Δ 5 = 2*4,5 + 2*0 + 2*2,5 – 12 = 2 > 0, т.е. включение в план изделия «Д» ценой 12 ед. нецелесообразно.
Задача 3
3.4. Промышленная группа
предприятий (холдинг)
Требуется:
1) Проверить продуктивность
технологической матрицы A=(
2) Построить баланс (заполнить
таблицу) производства и распределения
продукции предприятий
В соответствии с номером Вашего варианта ниже в таблице 1 выберите числовые значения для таблицы 2.
Таблица 1
Вариант |
Для первой строки |
Для второй строки |
Для третьей строки | ||||||||||
№ |
1А |
2А |
3А |
4А |
1Б |
2Б |
3Б |
4Б |
1В |
2В |
3В |
4В | |
4 |
0,1 |
0,0 |
0,1 |
100 |
0,1 |
0,0 |
0,2 |
300 |
0,2 |
0,1 |
0,0 |
160 | |
Таблица 2
Предприятие |
Коэффициенты прямых затрат аi j |
Конечный продукт Y | ||
1 |
2 |
3 | ||
|
1 2 3 |
0,1 0,1 0,2 |
0,0 0,0 0,1 |
0,1 0,2 0,0 |
100 300 160 |
Решение:
Для определения общего (валового) выпуска продукции первого, второго и третьего видов воспользуемся моделью Леонтьева в виде
Определяем матрицу-разность :
С помощью функции МОБР Мастера функций Excel найдем обратную матрицу:
Делаем вывод о продуктивности матрицы А, поскольку матрица (Е–А) неотрицательно обратима. Значит, мы можем найти матрицу-столбец объёмов валовой продукции Х в соответствии с моделью Леонтьева.
С помощью функции МУМНОЖ Мастера функций Excel найдем матрицу Х как произведение матриц В и У:
Таким образом, общие объёмы производства продукции цехов: х1= 135,8885; х2 = 358,1882; x3=222,9965
Распределение продукции между предприятиями холдинга на внутреннее потребление определяем из соотношения:
,т.е.
Х11=0,1*135,8885=13,58885;
X12=0,0*358,1882=0;
X13=0,1*222,9965=22,2997;
X21=0,1*135,8885=13,5889;
X22=0,0*358,1882=0;
X23=0,2*222,9965=44,5993;
X31=0,2*135,8885=27,1777;
X32=0,1*358,1882=35,8188;
X33=0,0*222,9965=0
В итоге плановая модель – баланс производства и распределения продукции холдинга будет иметь следующий вид:
Задача 4
4.4 Исследовать динамику экономического показателя на основе анализа одномерного временного ряда.
В течение девяти последовательных недель фиксировался спрос Y(t) (млн. руб.) на кредитные ресурсы финансовой компании. Временной ряд Y(t) этого показателя приведен в таблице.