Математические методы и модели / лекции12
.doc
2.1.5. Теория двойственности в анализе оптимальных решений экономических задач
Рассмотрим основные понятия и выводы специального раздела линейного программирования — теорию двойственности. На предыдущих лекциях мы выяснили, что любую задачу линейного программирования можно записать следующим образом:
-
f()=max ,
(2.18)
,
(2.19)
, j =.
(2.20)
С каждой задачей линейного программирования тесно связана другая линейная задача, называемая двойственной, первоначальная задача называется исходной или прямой.
Переменные двойственной задачи yi, называют объективно обусловленными оценками или двойственными оценками. Модель двойственной задачи имеет вид:
-
g()=min ,
(2.21)
,
(2.22)
, i =.
(2.23)
Каждая из задач двойственной пары фактически является самостоятельной задачей линейного программирования и может быть решена независимо от другой. Однако при определении симплексным методом оптимального плана одной из задач находится решение и другой задачи.
Двойственная задача по отношению к исходной составляется согласно следующим правилам:
1) целевая функция исходной задачи (2.18)-(2.20) формулируется на максимум, а целевая функция двойственной задачи (2.21)-(2.23 ) - на минимум, при этом в задаче на максимум все неравенства в функциональных ограничениях имеют вид , а в задаче на минимум – вид ;
2) матрица А=,
составленная из коэффициентов при неизвестных в системе ограничений (2.19) исходной задачи, и аналогичная матрица
=
в двойственной задаче получаются друг из друга транспонированием;
3) число переменных в двойственной задаче равно числу функциональных ограничений (2.19) исходной задачи, а число ограничений в системе (2.22) двойственной задачи - числу переменных в исходной задаче;
4) коэффициентами при неизвестных в целевой функции (2.21) двойственной задачи являются свободные члены в системе (2.19) ограничений исходной задачи, а правыми частями в ограничениях (2.22) двойственной задачи - коэффициенты при неизвестных в целевой функции (2.18) исходной задачи;
5) каждому ограничению одной задачи соответствует переменная другой задачи: номер переменной совпадает с номером ограничения; при этом ограничению, записанному в виде неравенства , соответствует переменная, связанная условием неотрицательности. Если функциональное ограничение исходной задачи является равенством, то соответствующая переменная двойственной задачи может принимать как положительные, так и отрицательные значения.
Математические модели пары двойственных задач могут быть симметричными и несимметричными. В несимметричных двойственных задачах система ограничений исходной задачи задается в виде равенств, а двойственной — в виде неравенств, причем в последней переменные могут быть и отрицательными. В симметричных задачах система ограничений как исходной, так и двойственной задачи задается неравенствами, причем на двойственные переменные налагается условие неотрицательности. В дальнейшем мы будем рассматривать только симметричные взаимодвойственные задачи линейного программирования.
Основные утверждения о взаимодвойственных задачах содержатся в двух следующих теоремах.
Первая теорема двойственности. Для взаимодвойственных ЗЛП имеет место один из взаимоисключающих случаев:
1. В прямой и двойственной задачах имеются оптимальные решения, при этом значения целевых функций на оптимальных решениях совпадают: max f()= min g().
2. В прямой задаче допустимое множество не пусто, а целевая функция на этом множестве не ограничена сверху. При этом у двойственной задачи будет пустое допустимое множество.
3. В двойственной задаче допустимое множество не пусто, а целевая функция на этом множестве не ограничена снизу. При этом у прямой задачи допустимое множество оказывается пустым.
4. Обе из рассматриваемых задач имеют пустые допустимые множества.
Вторая теорема двойственности (теорема о дополняющей нежесткости). Пусть = (х1, х2, ..., xп) - допустимое решение прямой задачи (2.18)-(2.20), a = (y1, y2, ..., ym) - допустимое решение двойственной задачи (2.21)-(2.23). Для того чтобы они были оптимальными решениями соответствующих взаимодвойственных задач, необходимо и достаточно, чтобы выполнялись следующие соотношения:
-
(2.24)
j =.
(2.25)
Условия (2.24) и (2.25) позволяют, зная оптимальное решение одной из взаимодвойственных задач, найти оптимальное решение другой задачи.
Рассмотрим еще одну теорему, выводы которой будут использованы в дальнейшем.
Теорема об оценках. Значения переменных yi в оптимальном решении двойственной задачи представляют собой оценки влияния свободных членов bi системы ограничений-неравенств прямой задачи на величину:
-
.
(2.26)
Решая ЗЛП симплексным методом, мы одновременно решаем двойственную ЗЛП.
Рассмотрим экономическую интерпретацию двойственной задачи на следующем примере.
▼
Пример 1. (Задача оптимального использования ресурсов). Пусть для выпуска четырех видов продукции P1, Р2, P3, Р4 на предприятии используют три вида сырья S1, S2 и S3. Объемы выделенного сырья, нормы расхода сырья и прибыль на единицу продукции при изготовлении каждого вида продукции приведены в табл. 6. Требуется определить план выпуска продукции, обеспечивающий наибольшую прибыль.
Составим экономико-математическую модель задачи оптимального использования ресурсов на максимум прибыли. В качестве неизвестных примем объем выпуска продукции j-го вида хj (j = 1,2,3,4).
Таблица 6
Вид сырья |
Запасы сырья |
Вид npодукции |
|||
P1 |
Р2 |
P3 |
Р4 |
||
S1 |
35 |
4 |
2 |
2 |
3 |
S2 |
30 |
1 |
1 |
2 |
3 |
S3 |
40 |
3 |
1 |
2 |
1 |
Прибыль |
14 |
10 |
14 |
11 |
Модель задачи: f()= 14х1 + 10х2 + 14х3 + 11х4 max
4х1+2х2+2х3+3х435
x1+x2+2х3+3х4 30
3х1+x2+2х3+x4 40
xj0, j= 1,2,3,4.
Теперь сформулируем двойственную задачу. Пусть некая организация решила закупить все ресурсы рассматриваемого предприятия. При этом необходимо установить оптимальную цену на приобретаемые ресурсы y1, y2, y3 исходя из следующих объективных условий:
1) покупающая организация старается минимизировать общую стоимость ресурсов;
2) за каждый вид ресурсов надо уплатить не менее той суммы, которую хозяйство может выручить при переработке сырья в готовую продукцию. Согласно первому условию общая стоимость сырья выразится величиной g()= 35y1 + 30y2 + 40y3 min. Согласно второму требованию вводятся ограничения: на единицу первого вида продукции расходуются четыре единицы первого ресурса ценой y1, одна единица второго ресурса ценой y2 и три единицы третьего ресурса ценой y3. Стоимость всех ресурсов, расходуемых на производство единицы первого вида продукции, равна 4y1 + y2 + 3y3 и должна составлять не менее 14, т. е. 4y1 + y2 + 3y314.
В результате аналогичных рассуждений относительно производства второго, третьего и четвертого видов продукции получаем систему неравенств:
4y1 + y2 + 3y314,
2y1 + y2 + y310,
2y1 + 2y2 + 2y314,
3y1 + 3y2 + y311.
По экономическому смыслу цены неотрицательные: ,,.
Получили симметричную пару взаимодвойственных задач. В результате решения данной задачи симплексным методом получен оптимальный план =(0;5;12,5;0); =(3;4;0).
▲
Экономический смысл первой теоремы двойственности следующий. План производства Х и набор оценок ресурсов Y оказываются оптимальными тогда и только тогда, когда прибыль от реализации продукции, определенная при известных заранее ценах продукции c1, с2, ..., сn, равна затратам на ресурсы по «внутренним» (определяемым только из решения задачи) ценам ресурсов y1, y2,…, yт. Для всех же других планов и обеих задач прибыль от продукции всегда меньше (или равна) стоимости затраченных ресурсов: f()g(), т. е. ценность всей выпущенной продукции не превосходит суммарной оценки имеющихся ресурсов. Значит величина g(Y) - f(X) характеризует производственные потери в зависимости от рассматриваемой производственной программы и выбранных оценок ресурсов.
Из первой теоремы двойственности следует, что при оптимальных производственной программе и векторе оценок ресурсов производственные потери равны нулю.
Экономический смысл первой теоремы двойственности можно интерпретировать и так: предприятию безразлично, производить ли продукцию по оптимальному плану и получить максимальную прибыль либо продать ресурсы по оптимальным ценам и возместить от продажи равные ей минимальные затраты на ресурсы.
Из второй теоремы двойственности в данном случае следуют такие требования на оптимальную производственную программу = (х1, х2, ..., xп) и оптимальный вектор оценок = (y1, y2, ..., ym):
-
если yi > 0, то;
если то yi=0, i =.
(2.27)
если xj > 0, то j =;
если то xj = 0, j =.
(2.28)
Условия (2.27) можно интерпретировать так: если оценка yi, единицы ресурса i-го вида положительна, то при оптимальной производственной программе этот ресурс используется полностью; если же ресурс используется не полностью, то его оценка равна нулю.
Из условия (2.28) следует, что если j-й вид продукции вошел в оптимальный план, то он в оптимальных оценках не убыточен; если же j-й вид продукции убыточен, то он не войдет в план, не будет выпускаться.
Кроме нахождения оптимального решения должно быть обеспечено получение дополнительной информации о возможных изменениях решения при изменении параметров системы. Эту часть исследования обычно называют анализом модели на чувствительность. Он необходим, например, в тех случаях, когда некоторые характеристики исследуемой системы не поддаются точной оценке.
Экономико-математический анализ решений осуществляется в двух основных направлениях: в виде вариантных расчетов по моделям с сопоставлением различных вариантов плана и в виде анализа каждого из полученных решений с помощью двойственных оценок.
Вариантные расчеты могут осуществляться при неизменной структуре самой модели (постоянном составе неизвестных, способов производства, ограничений задачи и одинаковом критерии оптимизации), но с изменением численной величины конкретных показателей модели. Вариантные расчеты могут проводиться также при варьировании элементов самой модели: изменении критерия оптимизации, добавлении новых ограничений на ресурсы или на способы производства их использования, расширения множества вариантов и т.д.
Одно из эффективных средств экономико-математического анализа — использование объективно обусловленных оценок оптимального плана. Такого рода анализ базируется на свойствах двойственных оценок. Выше мы установили общие математические свойства двойственных оценок для задач на оптимум любой экономической природы. Однако экономическая интерпретация этих оценок может быть совершенно различной для разных задач.
Перейдем к рассмотрению конкретных экономических свойств оценок уi оптимального плана. Сначала перечислим эти свойства, а затем проиллюстрируем их конкретными примерами.
Свойство 1. Оценки как мера дефицитности ресурсов и продукции.
Свойство 2. Оценки как мера влияния ограничений на функционал.
Свойство 3. Оценки как инструмент определения эффективности отдельных вариантов.
Свойство 4. Оценки как инструмент балансирования суммарных затрат и результатов.
▼
Пример 2. (Задача о планировании выпуска тканей). Пусть некоторая фабрика выпускает три вида тканей, причем суточное плановое задание составляет не менее 90 м тканей первого вида, 70 м — второго и 60 м— третьего. Суточные ресурсы следующие: 780 единиц производственного оборудования, 850 единиц сырья и 790 единиц электроэнергии, расход которых на один метр ткани представлен в табл. 7.
Цена за один метр ткани вида I равна 80 денежным единицам, II — 70 денежным единицам, III — 60 денежным единицам.
Таблица 7
Ресурсы |
Ткани |
||
I |
II |
III |
|
Оборудование |
2 |
3 |
4 |
Сырье |
1 |
4 |
5 |
Электроэнергия |
3 |
4 |
2 |
Необходимо определить, сколько метров ткани каждого вида следует выпустить, чтобы общая стоимость выпускаемой продукции была максимальной.
Составим модель задачи. Введем следующие обозначения. Неизвестными в задаче являются объемы выпуска ткани каждого вида:
хi — количество метров ткани вида I;
х2 — количество метров ткани вида II;
хз — количество метров ткани вида III.
-
f()=max ,
,
, j =,
.
С учетом имеющихся данных модель примет вид:
f()=80х1 + 70х2 + 60х3 max
2х1+3х2+4х3780
x1+4x2+5х3 850 Ограничения по ресурсам
3х1+4x2+2х3790
х190
x270 Плановое задание
х360
В результате решения задачи симплексным методом получен следующий оптимальный план: максимум общей стоимости продукции f()=19075 при
х1 = 112,5м — оптимальный план выпуска ткани вида I;
x2 = 70 м — оптимальный план выпуска ткани вида II;
х3 = 86,25 м — оптимальный план выпуска ткани вида Ш.
Решение двойственной задачи получим с использованием теорем двойственности. Введем обозначения:
y1 — двойственная оценка ресурса «оборудование»;
y2 — двойственная оценка ресурса «сырье»;
y3— двойственная оценка ресурса «электроэнергия»;
y4— двойственная оценка ткани вида I;
y5— двойственная оценка ткани вида II;
y6— двойственная оценка ткани вида III.
Модель двойственной задачи имеет вид:
g()= 780y1 + 850y2 + 790y3+90 y4+70 y5+60 y6 min
2y1 + y2 +3 y3+ y480,
3y1 + 4y2 + 4y3+ y570,
4y1 + 5y2 +2 y3+ y660.
,.
Из соотношений второй теоремы двойственности вытекают следующие условия:
для каждого ресурса:
если то yi=0,
если yi > 0, то
для задания по выпуску продукции:
-
если xj >Tj , то ym+j =0;
если ym+j < 0, то xj = Tj .
(2.29)
Для нашего примера в этих соотношениях m=3 (количество типов ресурсов).
Подставим значения х1 = 112,5, х2 = 70 и x3 = 86,25 в ограничения прямой задачи:
2*112,5+3*70 +4*86,25 = 780,
112,5 +4*70 + 5*86,25 = 823,75 < 850,
3*112,5+4*70 +2*86,25 = 790,
112,5 > 90,
70 = 70,
86,25 > 60.
Суточные ресурсы по оборудованию и электроэнергии использованы полностью. Сырье используется не полностью, имеется остаток в размере 850 - 823,75 = 26,25 (кг). План выпуска по тканям вида I и III перевыполнен.