
Задача 1
Для изготовления продукции двух видов А и B предприятие расходует ресурсы, а от реализации этой продукции получает доход. Информация о нормах затрат ресурсов на единицу выпускаемой продукции, запасах расходуемых ресурсов, имеющихся в распоряжении предприятия, и выручки от реализации готовой продукции приведены в таблице.
Наименование ресурсов |
Норма затрат на |
Объем ресурса |
|
Продукт А |
Продукт В |
||
Сырье (кг) |
4 |
1 |
391 |
Оборудование (ст.час.) |
1 |
3 |
232 |
Трудоресурсы (чел.час.) |
8 |
1 |
407 |
Цена реализации (руб.) |
407 |
232 |
|
Задача предприятия заключается в том, чтобы разработать программу выпуска, обеспечивающую получение максимальной выручки от реализации готовой продукции.
Требуется:
Построить математическую модель оптимизации выпуска продукции и записать ее в форме задачи линейного программирования.
Используя графический метод решения задачи линейного программирования, найти оптимальную программу выпуска продукции и максимум ожидаемой выручки.
Составив задачу, двойственную к задаче оптимизации выпуска продукции, найти её оптимальное решение, используя условия “дополняющей нежесткости”. Дать экономическую интерпретацию этого решения.
Решение
1. Пусть х1 – количество выпускаемого продукта А;
х2 – количество выпускаемого продукта В.
Искомая производственная программа X = (x1; x2) выпуска изделий А и В должна удовлетворять всем ресурсным ограничениям. Запишем их в математической форме.
4
х1
+ 1х2
<
391
1x1 + 3x2 < 232
8x1 + 1x2 < 407
x1 > 0, x2 > 0
Пусть z – выручка от продажи продуктов А и В. Задача состоит в таком выпуске продукции X = (x1; x2), который обеспечивает максимальную выручку, т.е. z = 407 x1 + 232 x2 max
2. Построим область допустимых решений (ОДР). Запишем уравнения граничных прямых для каждого из неравенств и по две точки на этих прямых.
4х1 + 1х2 = 391 (1) 1x1 + 3x2 = 232 (2) 8x1 + 1x2 = 407 (3)
x1 |
0 |
97,75 |
|
x1 |
0 |
232 |
|
x1 |
0 |
50,875 |
x2 |
391 |
0 |
|
x2 |
77,33 |
0 |
|
x2 |
407 |
0 |
При подстановке точки (0; 0) в левую часть неравенств они будут выполняться. Следовательно, искомые полуплоскости будут располагаться слева (ниже) граничных прямых. Получим ОДР в результате пересечения всех полуплоскостей в первом квадранте.
Находим градиент функции:
grad r = (dz/dx1; dz/dx2) = (407; 232)
Двигая линии уровня 407х1 + 232х2 = h вдоль вектора нормали, находим точку касания линии уровня и ОДР. Это и есть точка максимума функции z. В нашей задаче точка максимума X* лежит на пересечении граничных прямых (2) и (3). Находим ее координаты из системы:
1
х1
+ 3х2
= 232
8x1 + 1x2 = 407
x 1 = 43
x2 = 63
Оптимальная производственная программа X* = (43; 63) состоит в выпуске 43 изделий А и 63 изделий В. Ожидаемая выручка от их продажи:
z* = 407 * 43 + 232 * 63 = 32117 руб.
3. Исходная задача:
u1 4х1 + 1х2 < 391
u2 1x1 + 3x2 < 232
u3 8x1 + 1x2 < 407
x1 > 0, x2 > 0
z = 407 x1 + 232 x2 max
Двойственная задача:
4u1 + 1u2 + 8u3 407
1u1 + 3u2 + 1u3 232
u1 0, u2 0, u3 0
w = 391u1 + 232u2 + 407u3 min
Для того, чтобы допустимое решение X исходной задачи и допустимое решение U двойственной задачи были оптимальными, необходимо и достаточно, чтобы выполнялись следующие условия.
xj * vj = 0
uj * yj = 0
где vj = aijui – cj, yj = bi - aij*xj
Подставив найденные компоненты вектора X* = (43; 63) в условия получаем:
x1 * v1 = 0; x1 = 43 v1 = 0 4u1 + 1u2 + 8u3 = 407
x2 * v2 = 0; x2 = 63 v2 = 0 1u1 + 3u2 + 1u3 = 232
u1 * y1 = 0; y1 = 391 – 4*43 – 1*63 = 156 u1 = 0
u2 * y2 = 0; y2 = 232 – 1*43 – 3*63 = 0 u2 0
u3 * y3 = 0; y3 = 407 – 8*43 – 1*63 =0 u3 0
Получаем систему уравнений:
4
u1
+ 1u2
+ 8u3
= 407 u1*
= 0
1u1 + 3u2 + 1u3 = 232 u2* = 63
u1= 0 u3* = 43
Значение целевой функции двойственной задачи на этом решении:
w* = 391*0 + 232*63 + 407*43 =32117 руб.
Получены следующие результаты расчета модели:
x* = (43; 63); u* = (0; 63; 43)
z* = w* = 32117 руб.
Переменная ui характеризует абсолютный прирост оптимизируемого показателя z в случае увеличения объема i-го ресурса на одну единицу.
Переменная yj показывает, сколько остается i-го ресурса после выполнения производственной программы X.
Переменная vj интерпретируется как возможный убыток по полезности при выпуске j-го продукта.
Оценка u1* = 0 руб/кг. показывает, что объем сырья является избыточным. Уменьшение (в пределах интервала устойчивости) или увеличение объемов сырья не повлияет на величину ожидаемой выручки.
Оценка u2* = 63 руб/ ст.час показывает, что если оборудование увеличить на 1 ст.час, то при прочих равных условиях максимальная выручка увеличится на 63 руб., а если уменьшить на 1 ст.час, то снизится на 63 руб.
Оценка u3* = 43 руб/чел.час. показывает, что если объем используемого трудоресурса увеличить (уменьшить) на 1 чел.час, то максимальная выручка увеличится (уменьшится) на 43 руб.