4.1. АНАЛИЗ РЕШЕНИЯ (МОДЕЛИ) НА ЧУВСТВИТЕЛЬНОСТЬ
Модель линейного программирования является как бы «моментаьным снимком» реальной ситуации, когда параметры модели (коэффициенты целевой функции и неравенств ограничений) предполагаются неизменными. Естесствено изучить влияние изменения параметров модели на полученное оптимальное решение задачи ЛП. Такое исследование называется анализом на чувствительность. В этом разделе анализ чувствительности основывается на графическом решении задачи ЛП.
Пример 4.1.
Компания производит краску для внутренних и наружных работ из сырья двух типов: М1 и М2.
Необходимая информация представлена в следующей таблице:
|
Расход на 1 тонну |
сырья краски |
Максимально возможный |
|
Для наружных работ |
Для внутренних работ |
ежедневный расход сырья |
Сырье М1 |
6 |
4 |
24 |
Сырье М2 |
1 |
2 |
6 |
Доход на тонну краски (тыс.дол) |
5 |
4 |
|
Отдел маркетинга компании ограничил ежедневное производство краски для внутренних работ до 2 т, а кроме того этот показатель не должен превышать более чем на тонну показатель выпуска краски для внешних работ.
Цель компании:
Определить оптимальное соотношение между видами выпускаемой продукции для максимизации общего ежедневного дохода.
Составленная математическая модель задачи выглядит следующим образом:
Максимизировать Z(x)=5X1+4X2
при выполнении ограничений
6Х1+4Х2 ≤24
Х1+2Х2≤6
Х2≤2
X2-X1≤1
Х1≥0
Х2≥0
В результате применения графического метода решения ЗЛП, рассмотренного в параграфе 3.2, получен график (рис.4.1.1)
Рис.4.1.1. Решение задачи
Решением задачи является точка с координатами: Х1=3;Х2=1,5. Целевая функция при таком решении принимает значение Z=21 тыс. дол.
Проведем для данной задачи анализ чувствительности. Рассмотрим два случая:
1) изменение коэффициентов целевой функции;
2) изменение значений констант в правой части неравенств ограничений.
1. Изменение коэффициентов целевой функции. В общем виде целевую функцию задачи ЛП можно записать следующим образом:
Максимизировать или минимизировать Z(x)=С1 X1+C2X2
Изменение значений коэффициентов С1 и С2 приводит к изменению угла наклона прямой Z. Графический способ решения показывает, что это может привести к изменению оптимального решения: оно будет достигаться в другой угловой точке пространства решений. Вместе с тем, очевидно, существуют интервалы изменения коэффициентов С1 и С2, когда текущее оптимальное решение сохраняется. Задача анализа чувствительности и состоит в получении такой информации. В частности, представляет интерес определение интервала оптимальности для отношения С1 /С2 (или, что то же самое, для С2 /С1); если значение отношения С1 /С2 не выходит за пределы этого интервала, то оптимальное решение в данной модели сохраняется неизменным.
На рис.4.1.1. видно, что функция Z(x)=5X1+4X2 достигает максимального значения в угловой точке С. При изменении коэффициентов целевой функции Z(x)=С1 X1+C2X2 точка С останется точкой оптимального решения до тех пор, пока угол наклона линии Z будет лежать между углами наклона двух прямых, пересечением которых является точка С. Этими прямыми являются 6Х1+4Х2 ≤24 (ограничение на сырье М1) и Х1+2Х2≤6 (ограничение на сырье М2). Алгебраически это можно записать следующим образом:
или
.
В первой системе неравенств условие означает, что прямая, соответствующая целевой
функции, не может быть горизонтальной. Аналогичное условие в следующей системе
неравенств означает, что эта же прямая не может быть вертикальной. Из рис.4.1.2
видно, что интервал оптимальности данной задачи (он определяется двумя
прякающимися в точке С) не разрешает целевой функции быть ни горизонтальной. Ни
вертикальной. Таким образом, получено две системы неравенств, определяющие
интервал оптимальности в данной задаче.
Рис.4.1.2. Интервал оптимальности
Итак, если коэффициенты С1 и С2 удовлетворяют приведенным выше неравенствам, оптимальное решение по прежнему будет достигаться в точке С. Отметим, если прямая Z(x)=С1 X1+C2X2 совпадет с прямой Х1+2Х2≤6, то оптимальным решением будет любая точка отрезка CD. Аналогично, если прямая, соответствующая с целевой функцией, совпадет с прямой 6Х1+4Х2 ≤24, тогда любая точка отрезка ВС будет оптимальным решением. Однако очевидно, что в обоих случаях точка С остается точкой оптимального решения.
Приведенные выше неравенства можно использовать при
определении интервала оптимальности для какого-либо одного коэффициента целевой
функции, если предположить, что другой коэффициент остается неизменным.
Например, зафиксируем значение коэффициента С2 (пусть С2=4),
тогда интервал оптимальности для коэффициента С1 получаем из
неравенств путем подстановки туда значения С2=4.
После выполнения элементарных арифметических опреаций получаем неравенства для
коэффициента С1: 2≤С1≤6.
Это означает, что при фиксированной цене на краску для внутренних работ, цена на краску для наружных работ может меняться в интервале от 2 тыс. дол.за тонну до 6 тыс.дол. за тонну, при том, что оптимальное соотношение (решение) останется неизменным.
Аналогично, если зафиксировать значение коэффициента С1
(пусть С1=5), тогда из неравенства получаем интервал оптимальности для
коэффициента С2:
.
2. Изменение значений констант в правой части неравенств ограничений. Стоимость ресурсов. Во многих моделях линейного программирования ограничения трактуются как условия ограниченности ресурсов. В таких ограничениях правая часть неравенств является верхней границей количества доступных ресурсов. Рассмотрим на примере чувствительность оптимального решения к изменению ограничений, накладываемых на ресурсы. Такой анализ задачи ЛП предлагает простую меру чувствительности решения, называемую стоимостью единицы ресурса; при изменении количества доступных ресурсов (на единицу) значение целевой функции в оптимальном решении изменится на стоимость единицы ресурса.
В примере 4.1. первые два неравенства представляют собой ограничения на использование сырья М1 и М2 соответственно. Определим стоимость единиц этих ресурсов.
В данной задаче оптимальное решение достигается в точке С, являющейся точкой пересечения прямых, соответствующих ограничениям на сырье М1 и М2. При изменении уровня доступности материала М1 (увеличение или уменьшение текущего уровня, равного 24 т) точка С оптимального решение «плывет» вдоль отрезка DG (рис.4.1.3).
Рис.4.1.3.
Любое изменение уровня доступности материала М1, приводящее к выходу точки пересечения С из этого отрезка, ведет к неосуществимости оптимального решения в точке С. Поэтому можно сказать, что концевые точки D=(2,2) и G=(6,0) отрезка DG определяют интервал осуществимости для ресурса М1. Количество сырья М1, соответсятвующего точке D=(2,2), равно 6Х1+4Х2=20 т. Аналогично, количество сырья, соответствующего точке G=(6,0), равно 36 т. Таким образом, интервал осуществимости для ресурса М1 составляет 20≤М1≤36. Если определить М1 как М1=24+D1, где D1– отклонение количества материала М1 от текущего уровня в 24 т, тогда последние неравенства можно переписатькак 20≤24+D1≤36 или -4≤D1≤12. Это означает, что текущий уровень ресурса М1 может быть уменьшен не более, чем на 4 т и увеличен не более, чем на 12 т. В этом случае структура оптимального решения не изменится.
Вычислим стоимость единицы материала М1. При изменении количества сырья М1 от 20 до 36 тонн, значения целевой функции Z будут соответствовать положению точки С на отрезке DG. Обозначив через y1 стоимость единицы ресурса М1, получим следующую формулу:
Если точка С совпадает с точкой D=(2,2), то Z=5*2+4*2=18 (тыс. дол.), если же точка С совпадает с точкой G=(6,0), тогда Z=5*6+4*0=30 (тыс.дол.). Отсюда следует, что
(тыс.дол.на тонну
материала М1).
Этот результат показывает, что изменение количества ресурса М1 на одну тонну приводит к изменению в оптимальном решении значсения целевой функции на 750 дол.
Рассмотрим ресурс М2. На рис. 4.1.4. видно, что интервал осуществимости для ресурса М2 определяется концевыми точками В и Н отрезка ВН, где В=(4,0) и Н=(8/3,2).
Рис.4.1.4.
Точка Н находится на пересечении прямых ЕD и ВС. Находим, что количество сырья М2,
соответствующего точке В, равно Х1+2Х2=4+2*0=4т, а в точке Н- 20/3т. Значение целевой функции в
точке В равно Z=5*4+4*0=20 тыс.дол., а в точке Н- Z=5*8/3+4*2=64/3 тыс.дол. Отсюда следует, что
количество сырья М2 может изменятся от 4 до 20/3 тонн, а стоимость единицы
ресурса М2, обозначенная как y2, равна (тысяч долларов на
тонну материала М2).
ДОМАШНЕЕ ЗАДАНИЕ 16
Решить задачу графическим или симплекс-методом и ответить на вопросы 1-10.
Вариант 1-5. Для приготовления двух видов продукции (A, B) используют три вида сырья. Ресурсы сырья, норма его расхода на единицу продукции и цена продукции заданы в соответствующей таблице.
1. Определить план выпуска продукции из условия максимизации его стоимости.
2. Определите статус, ценность каждого ресурса и его приоритет при решении задачи увеличения запаса ресурсов.
3. Определите максимальный интервал изменения запасов каждого из ресурсов, в пределах которого структура оптимального решения, то есть номенклатура выпускаемой продукции, остается без изменения.
4. Определите суммарную стоимостную оценку ресурсов, используемых при производстве единицы каждого изделия. Производство какой продукции нерентабельно?
5. На сколько уменьшится стоимость выпускаемой продукции при принудительном выпуске единицы нерентабельной продукции?
6. На сколько можно снизить запас каждого из ресурсов, чтобы это не привело к уменьшению прибыли.
7. Определите изменение стоимости продукции и количество выпускаемых .изделий при увеличении второго вида сырья на Z единиц.
8. Определите оптимальное решение
задачи для случая, когда вектор ресурсов задан в виде -строки.
9. Определите интервалы изменения цен на каждую продукцию, при которых сохраняется структура оптимального плана.
10. На сколько нужно снизить затраты каждого вида сырья на единицу продукции, чтобы сделать производство нерентабельного изделия рентабельным?
1.
|
Сырье |
Норма расходов |
Ресурсы |
|||||
|
A |
B |
|
|||||
I |
2 |
1 |
2400 |
|
||||
II |
1 |
5 |
1800 |
|
||||
III |
3 |
- |
2000 |
|
||||
Цена ( |
7,5 |
3 |
|
|
||||
Z=500, =(2000,1500,2000)
2.
Сырье |
Норма расходов |
Ресурсы |
||
A |
B |
|
||
I |
1 |
1 |
4500 |
|
II |
2 |
3 |
1200 |
|
III |
3 |
- |
2300 |
|
Цена ( |
7,5 |
3 |
|
|
Z=300, =(1500,2000, 2000)
3.
Сырье |
Норма расходов |
Ресурсы |
||
A |
B |
|
||
I |
4,5 |
1 |
2400 |
|
II |
1 |
5 |
820 |
|
III |
- |
10 |
2000 |
|
Цена ( |
10,5 |
3 |
|
|
Z=700, =(2000,2880,1500)
4.
Сырье |
Норма расходов |
Ресурсы |
|
A |
B |
|
|
I |
2 |
1 |
2600 |
II |
1,5 |
5 |
2200 |
III |
3 |
2 |
1000 |
Цена ( |
9 |
3 |
|
Z=450, =(2000,1500,700)
5.
Сырье |
Норма расходов |
Ресурсы |
|
A |
B |
|
|
I |
2 |
1 |
2700 |
II |
1 |
5 |
3200 |
III |
3 |
- |
1500 |
Цена ( |
13 |
3 |
|
Z=500, =(1000,2500,500)
Вариант 6-10. Из 2 видов кормов необходимо составить рацион, в состав которого должно входить не менее в1 ед. вещества А, в2 ед. вещества В и в3 ед. вещества С. Количество единиц вещества, содержащегося в 1 кг корма каждого вида, указано в соответствующей таблице. В ней же приведена цена 1 кг корма каждого вида.
1. Составить рацион, содержащий не менее нужного количества указанных питательных веществ и имеющий минимальную стоимость.
2. Определите, все ли виды кормов входят в рацион, ценность дополнительной единицы каждого питательного вещества и его приоритет при решении задач уменьшения стоимости рациона.
3. Определите суммарную стоимостную оценку питательных веществ в единице каждого корма, использование какого вида корма нерентабельно.
4. Содержание какого из питательных веществ превышает заданный минимальный уровень и на сколько?
5. Определите максимально возможное уменьшение содержания каждого из питательных веществ в рационе, при котором структура рациона остается без изменений.
6. На сколько уменьшится стоимость рациона и используемое количество кормов при снижении минимального уровня потребления питательного вещества В до Z ед.
7. Определите интервал изменения цен на каждый вид корма, при котором сохраняется структура рациона.
8. Возможно ли сделать выгодным использование корма, не вошедшего в рацион.
9. На сколько увеличится стоимость рациона при принудительном включении в рацион 1 кг нерентабельного вида корма.
10.На сколько нужно снизить минимальный уровень потребления каждого из питательных веществ, чтобы уменьшить стоимость рациона на 10%?.
6.
Вещество |
Количество единиц вещества, содержащегося в 1 кг корма каждого вида |
|
|
1 |
2 |
||
A |
10 |
5 |
|
B |
- |
10 |
|
C |
20 |
7 |
|
Цена 1 кг корма (руб) |
9 |
11 |
|
=(400,180,200); Z=70
7.
Вещество |
Количество единиц вещества, содержащегося в 1 кг корма каждого вида |
|
1 |
2 |
|
A |
12 |
5 |
B |
- |
4 |
C |
22 |
7 |
Цена 1 кг корма (руб) |
11 |
9 |
=(400,180,200); Z=30
8.
Вещество |
Количество единиц вещества, содержащегося в 1 кг корма каждого вида |
|
1 |
2 |
|
A |
10 |
- |
B |
20 |
14 |
C |
- |
7 |
Цена 1 кг корма (руб) |
9 |
11 |
=(400,180,200); Z=110
9.
Вещество |
Количество единиц вещества, содержащегося в 1 кг корма каждого вида |
|
1 |
2 |
|
A |
10.5 |
5 |
B |
- |
10 |
C |
20 |
- |
Цена 1 кг корма (руб) |
16 |
15 |
=(400,180,200); Z=60
10.
Вещество |
Количество единиц вещества, содержащегося в 1 кг корма каждого вида |
|
1 |
2 |
|
A |
10 |
5 |
B |
- |
7 |
C |
20 |
7 |
Цена 1 кг корма (руб) |
9 |
11 |
=(400,180,200); Z=30
Назад к разделу "3.7. РЕШЕНИЕ ЗЛП С ПОМОЩЬЮ MS EXCEL"
Вперед к разделу "4.2. АНАЛИЗ РЕШЕНИЯ ЗЛП НА ОСНОВЕ ОТЧЁТОВ MS EXCEL"