Графический метод решения задач линейного программирования
Наиболее простым и наглядным методом решения задачи линейного программирования (ЛП) является графический метод. Он применяется для решения задач ЛП с двумя переменными.
Рассмотрим задачу ЛП в стандартной форме записи:
Положим n=2, т.е. рассмотрим эту задачу на плоскости. Пусть система неравенств совместна (имеет хотя бы одно решение).
Каждое неравенство этой системы геометрически определяет полуплоскость с граничной прямой ai1 x1 + ai2 x2 = bi , i=1,2,…m. Условия неотрицательности определяют полуплоскости, соответственно, с граничными прямыми x1=0,x2 =0. Система совместна, поэтому полуплоскости, как выпуклые множества, пересекаясь, образуют общую часть, которая является выпуклым множеством и представляет собой совокупность точек, координаты каждой из которых являются решением данной системы. Совокупность этих точек называют многоугольником решений. Он может быть точкой, отрезком, лучом, многоугольником, неограниченной многоугольной областью.
Таким образом, геометрически задача линейного программирования (ЗЛП) представляет собой отыскание такой точки многоугольника решений, координаты которой доставляют линейной функции цели максимальное (минимальное) значение, причем допустимыми решениями являются все точки многоугольника решений.
Линейное уравнение описывает множество точек, лежащих на одной прямой. Линейное неравенство описывает некоторую область на плоскости. Определим, какую часть плоскости описывает неравенство 2х1+3х2£ 12. Во-первых, построим прямую 2х1+3х2=12. Эта прямая проходит через точки (6, 0) и (0, 4). Для того чтобы определить, какая полуплоскость удовлетворяет неравенству необходимо выбрать любую точку на графике, не принадлежащую прямой, и подставить ее координаты в неравенство. Если неравенство будет выполняться, то данная точка является допустимым решением и полуплоскость, содержащая точку, удовлетворяет
неравенству. Удобной для использования при подстановке в неравенство является начало координат. Подставим х1=х2=0 в неравенство 2х1+3х2£12. Получим 2´0+3´0£12. Данное утверждение является верным, следовательно, неравенству 2х1+3х2£12 соответствует нижняя полуплоскость, содержащая точку (0.0). Это отражено на графике, изображенном на рис.1.
Рис. 1. Неравенству 2х1+3х2£12 соответствует нижняя полуплоскость.
Аналогично можно изобразить графически каждое ограничение задачи линейного программирования.
Решением каждого неравенства системы ограничений ЗЛП является полуплоскость, содержащая граничную прямую и расположенная по одну сторону от нее. Пересечение полуплоскостей, каждая из которых определяется соответствующим неравенством системы, называется областью допустимых решений или областью определения. Необходимо помнить, что область допустимых решений удовлетворяет условиям неотрицательности (xj ³0, j=1,…,n). Координаты любой точки, принадлежащей области определения являются допустимым решением задачи.
Для нахождения экстремального значения целевой функции при графическом решении задач ЛП используют вектор–градиент, координаты которого являются частными производными целевой функции, т.е.
.
Этот вектор показывает направление наискорейшего изменения целевой функции. Прямая , перпендикулярная вектору–градиенту, является линией уровня целевой функции. В любой точке линии уровня целевая функция принимает одно и то же значение. Приравняем целевую функцию постоянной величине “a”. Меняя значение “a”, получим семейство параллельных прямых, каждая из которых является линией уровня.
Важное свойство линии уровня линейной функции состоит в том, что при параллельном смещении линии в одну сторону уровень только возрастает, а при смещении в другую сторону – убывает.
С геометрической точки зрения в задаче линейного программирования ищется такая угловая точка или набор точек из допустимого множества решений, на котором достигается самая верхняя (нижняя) линия уровня, расположенная дальше (ближе) остальных в направлении наискорейшего роста.
Графический метод решения ЗЛП состоит из следующих этапов.
1. Строится многоугольная область допустимых решений ЗЛП – ОДР,
2. Строится вектор-градиент ЦФ в какой-нибудь точке Х0 принадлежащей ОДР –
.
3. Линия уровня C1x1+C2x2 = а (а–постоянная величина) - прямая, перпендикулярная вектору –градиенту – передвигается в направлении этого вектора в случае максимизации f(x1,x2) до тех пор, пока не покинет пределов ОДР. Предельная точка (или точки) области при этом движении и является точкой максимума f(x1,x2).
4. Для нахождения ее координат достаточно решить два уравнения прямых, получаемых из соответствующих ограничений и дающих в пересечении точку максимума. Значение f(x1,x2), найденное в получаемой точке, является максимальным.
При минимизации f(x1,x2) линия уровня перемещается в направлении, противоположном вектору-градиенту. Если прямая при своем движении не покидает ОДР, то соответствующий максимум или минимум f(x1,x2) не существует.
Если линия уровня параллельна какому-либо функциональному ограничению задачи, то оптимальное значение ЦФ будет достигаться в любой точке этого ограничения, лежащей между двумя оптимальными угловыми точками, и, соответственно, любая из этих точек является оптимальным решением ЗЛП.
Рассмотрим графическое решение задач линейного программирования на следующем примере.
Задача. о планировании выпуска продукции пошивочному предприятию.
Намечается выпуск двух видов костюмов - мужских и женских. На женский костюм требуется 1 м шерсти, 2 м лавсана и 1 человеко-день трудозатрат. На мужской костюм - 3,5 м шерсти, 0,5 м лавсана и 1 человеко-день трудозатрат. Всего имеется 350 м шерсти, 240 м лавсана и 150 человеко-дней трудозатрат. Tребуется определить, сколько костюмов каждого вида необходимо сшить, чтобы обеспечить максимальную прибыль, если прибыль от реализации женского костюма составляет 10 денежных единиц, а от мужского - 20 денежных единиц. При этом следует иметь в виду, что необходимо сшить не менее 60 мужских костюмов.
Модель задачи.
Введем следующие обозначения: х1 - число женских костюмов; x2 - число мужских костюмов.
Прибыль от реализации женских костюмов составляет 10х1, а от реализации мужских 20х2, т.е. необходимо максимизировать целевую функцию
f(x) = 10´ х1 + 20´ х2 -> max.
Ограничения задачи имеют вид:
х1 + х2 £ 150
2 х1 + 0.5 х2 £ 240
х1 + 3.5 х2 £ 350
х2³ 60
х1 ³ 0
Первое ограничение по труду х1 + х2 £ 150. Прямая х1 + х2 = 150 проходит через точки (150, 0) и (0, 150).
Рис. 2. Решением первого неравенства является нижняя полуплоскость.
Второе ограничение по лавсану 2 х1 + 0.5 х2 £ 240. Прямая 2 х1 + 0.5 х2 = 240 проходит через точки (120, 0) и (0, 480). Третье ограничение по шерсти х1 + 3.5 х2 £ 350. Добавим четвертое ограничение по количеству мужских костюмов х2 ³ 60. Решением этого неравенства является полуплоскость, лежащая выше прямой х2 = 60. На рис.3. заштрихована область допустимых решений.
Рис. 3. Заштрихована область допустимых решений.
Для определения направления движения к оптимуму построим вектор-градиент Ñ, координаты которого являются частными производными целевой функции, т.е. = (10;20).
Что бы построить этот вектор, нужно соединить точку (10;20) с началом координат. При максимизации целевой функции необходимо двигаться в направлении вектора-градиента, а при минимизации — в противоположном направлении. Для удобства можно строить вектор, пропорциональный вектору Ñ. Так, на рис. 2.1.6. изображен вектор градиент (30;60).
В нашем случае движение линии уровня будем осуществлять до ее выхода из области допустимых решений. в крайней, угловой точке достигается максимум целевой функции. Для нахождения координат этой точки достаточно решить два уравнения прямых, получаемых из соответствующих ограничений и дающих в пересечении точку максимума:
х1 + 3.5 х2 = 350
х1 + х2 = 150 .
Отсюда легко записать решение исходной ЗЛП: max f(x) = 2300 и достигается при x1=70 и x2=80 (рис. 4.)
Рис.4. Максимум целевой функции достигается в точке (70, 80).
Задача.Для изготовления двух видов продукции А1 и А2 используют три вида ресурсов S1, S2, S3, запасы которых составляют 18, 16 и 5 усл.ед. Расход ресурсов на 1 ед. продукции приведен в таблице:
Виды ресурсов | Запасы ресурсов | Расходы ресурсов на 1 изд. | |
А1 | А2 | ||
S1 | 18 | 1 | 3 |
S2 | 16 | 2 | 1 |
S3 | 5 | - | 1 |
Прибыль | 2 руб. | 3 руб. |
Необходимо составить такой план производства продукции, который обеспечит наибольшую прибыль от ее реализации.
Составим экономико-математическую модель (ЭММ) задачи.
Пусть надо выпустить изделий A1 - x1 шт., а изделий А2 - x2 шт. Тогда прибыль F = 2x1 + 3x2 Þ max
x1 + 3x2 £ | 18 |
2x1 + x2 £ | 16 |
x2 £ | 5 |
x1 ³ 0, | x2 ³ 0 |
Решим задачу графически.
1) | x1 + 3x2 £ 18 | ||
x1 + 3x2 = 18 (0; 6) (18; 0) | |||
к.т. (0; 0), 0 + 3*0 < 18 (в) – входит | |||
2) | 2x1 + x2 £ 16 | ||
2x1 + x2 = 16 (0; 16) (8; 0) | |||
к.т. (0; 0), 2*0 + 0 < 16 (в) – входит | |||
3) | x2 £ 5 | ||
x2 = 5, x2 < 5 - ниже прямой | |||
4) | x1 ³ 0 - правее ОX2 | ||
5) | x2 ³ 0 - выше ОX1 | ||
Линия уровня | F = 2x1 + 3 x2 F = 0 | ||
2x1 + 3x2 = 0 (0; 0) (3; -2) | |||
q = {2; 3} - указывает направление возрастания F. | |||
max F достигается в т. С
т.С | x1 + 3 x2 = 18 | Þ - | 2 x1 + 6 x2 = 36 | Þ | 5 x2 = 20 | Þ |
2x1 + x2 = 16 | 2 x1 + x2 = 16 | x1 + 3 x2 = 18 | ||||
Þ | x2 = 4 | |||||
x1 = 6 |
Таким образом, необходимо выпустить x1 = 6 шт. изделий А1, x2 = 4 шт. изделий А2, чтобы получить max F = 2*6 + 3*4 = 24 ден.ед.
Симплексный метод решения задачи линейного программирования
Для решения ЗЛП существует универсальный метод – метод последовательного улучшения плана или симплекс-метод, который состоит из двух вычислительных процедур: симплекс-метода с естественным базисом и симплекс-метода с искусственным базисом (М-метод).
Выбор конкретной вычислительной процедуры осуществляется после приведения исходной ЗЛП к каноническому виду (КЗЛП):
В теории линейного программирования (ЛП) показано, что оптимальное решение ЗЛП связано с угловыми (крайними) точками многогранника решений, которым отвечают опорные планы (неотрицательные базисные решения системы уравнений КЗЛП). Каждый из опорных планов определяется системой m линейно независимых векторов, содержащихся в данной системе из n векторов А1, А2,…, Аn. Верхняя граница количества опорных планов, содержащихся в данной задаче, определяется числом сочетаний Сnm.
Задача. Получить решение по модели задачи об оптимальном использовании ограниченных ресурсов (решить ЗЛП):
Приведем задачу к каноническому виду путем введения дополнительных переменных x 3 и x4:
Найдем все опорные планы КЗЛП. Используя метод Жордана-Гаусса выписываем все базисные решения системы уравнений:
Последовательно будем иметь:
Х1 = (0,0,300,150); Х2= (150,0,150,0); Х3= (0,150,-150,0); Х4= (75,75,0,0); Х5= (300,0,0,-150); Х6= (0,100,0,50).
Среди этих базисных решений четыре опорных:
Х1 = (0,0,300,150); Х2= (150,0,150,0); Х4= (75,75,0,0); Х6= (0,100,0,50).
Указанным опорным планам на рис.5 отвечают соответственно следующие угловые точки и значения ЦФ в них:
А(0,0) и f(0,0)=0; Д(150,0) и f(150,0)=300; С(75,75) и f(75,75)=375; В(0,100) и f(0,100)=300.
Согласно теории ЛП оптимальное решение содержится среди опорных планов.
Таким образом, максимальное значение, равное 375, целевая функция f (x1,x2) достигает на опорном плане Х4= (75,75,0,0), т.е. оптимальное решение рассматриваемой ЗЛП: x1 = 75, x2 = 75.
Понятно, что при больших m и n найти оптимальный план, перебирая указанным выше способом (прямым перебором) все опорные ЗЛП весьма трудно. Поэтому необходимо иметь схему, позволяющую осуществлять упорядоченный переход от одного опорного плана к другому. Такой схемой является симплексный метод.
Симплекс-метод с естественным базисом
Для применения симплекс-метода с естественным базисом КЗЛП должна содержать единичную подматрицу размером mxm – в этом случае очевиден начальный опорный план (неотрицательное базисное решение системы ограничений КЗЛП).
Для определенности предположим, что первые m векторов матрицы системы уравнений составляют единичную матрицу. Тогда первоначальный опорный план очевиден – (b1, b2 ,…, bm ,0,…,0).
Проверка на оптимальность опорного плана проходит с помощью признака оптимальности, переход к другому опорному плану проводится с помощью преобразований Жордана-Гаусса при использовании математического признака оптимальности. Полученный опорный план снова проверяется на оптимальность и так далее. Процесс заканчивается за конечное число шагов, причем на последнем шаге либо выявляется неразрешимость задачи (конечного оптимума нет), либо получается оптимальный опорный план и соответствующее ему оптимальное значение ЦФ.
Математический признак оптимальности состоит из следующих двух теорем:
1. Если для всех векторов А1, А2,…, Аn выполняется условие
где ,
то полученный опорный план является оптимальным.
2. Если для некоторого вектора, не входящего в базис, выполняется условие , то можно получить новый опорный план, для которого значение ЦФ будет больше исходного, при этом могут быть два случая:
- если все компоненты вектора, подлежащего вводу в базис, неположительны, то ЗЛП не имеет решения (конечного оптимума нет);
- если имеется хотя бы одна положительная компонента у вектора, подлежащего вводу в базис, то можно получить новый опорный план.
На основании признака оптимальности в базис вводится вектор Ак , давший минимальную отрицательную величину симплекс разности:
Чтобы выполнялось условие неотрицательности значений опорного плана, выводится из базиса вектор Аr, который дает минимальное положительное оценочное отношение
Строка Аr называется направляющей, столбец Ак и элемент ar к – направляющими.
Элементы направляющей строки в новой симплекс-таблице вычисляются по формулам:
а элементы i-й строки – по формулам:
Значения нового опорного плана рассчитываются по формулам:
для i = r ;
Процесс решения продолжают либо до получения оптимального плана, либо до установления неограниченности ЦФ. Если среди симплекс-разностей (оценок) оптимального плана нулевые только оценки, соответствующие базисным векторам, то это говорит о единственности оптимального плана. Если же нулевая оценка соответствует вектору, не входящему, то в общем случае это означает, что оптимальный план не единственный.
Примечание. Для использования приведенной процедуры к минимизации линейной функции f (x1,x2,…, xn) следует искать максимум - f (x1,x2,…, xn), затем полученный максимум взять с противоположным знаком. Оптимальное решение то же.
Пример. Получить решение по модели:
Эта задача (модель) линейного программирования, приведем ее к каноническому виду путем введения дополнительных переменных x 3 и x4:
КЗЛП имеет необходимое число единичных столбцов, т.е. обладает очевидным начальным опорным планом (0,0,300,150). Решение осуществляется симплекс-методом с естественным базисом с оформлением расчетов в симплекс-таблицах:
Номер | В | |||||||
симплекс- | Базис | план | Q | |||||
таблицы | ||||||||
А3 | ||||||||
А4 | ||||||||
f(x) | -2 | -3 | ||||||
А2 | 1/3 | 1/3 | ||||||
I | А4 | 2/3 | -1/3 | |||||
f(x) | -1 | |||||||
А2 | 1/2 | -1/2 | ||||||
II | А1 | -1/2 | 3/2 | |||||
f(x) | 1/2 | 3/2 |
В симплекс-таблице II получен оптимальный опорный план, поскольку все симплекс-разности (оценки) j. Оптимальные значения переменных равны: x1*=75, x2* =75 (основные переменные), x3* =0, x4* =0 (дополнительные переменные). Максимальное значение целевой функции равно 375.
Таким образом, в рассмотренной выше задаче об оптимальном использовании ограниченных ресурсов, оптимальная производственная программа состоит в выпуске 75ед. изделий первого вида и 75ед. изделий второго вида. С этой программой связана максимальная выручка от реализации готовой продукции – 375 у.е.
Теория двойственности в задачах линейного программирования. Анализ полученных оптимальных решений.
С каждой задачей линейного программирования тесно связана другая линейная задача, называемая двойственной; первоначальная задача называется исходной или прямой.
Связь исходной и двойственной задач заключается, в частности, в том, что решение одной из них может быть получено непосредственно из решения другой.
Переменные двойственной задачи yi называют объективно обусловленными оценками, или двойственными оценками, или «ценами» ресурсов, или теневыми ценами.
Каждая из задач двойственной пары фактически является самостоятельной задачей линейного программирования и может быть решена независимо от другой.
Двойственная задача по отношению к исходной составляется согласно следующим правилам:
1) целевая функция исходной задачи формулируется на максимум, а целевая функция двойственной задачи— на минимум, при этом в задаче на максимум все неравенства в функциональных ограничениях имеют вид £, в задаче на минимум — вид ³;
2) матрица А, составленная из коэффициентов при неизвестных в системе ограничений исходной задачи и аналогичная матрица Ат в двойственной задаче получаются друг из друга транспонированием;
3) число переменных в двойственной задаче равно числу функциональных ограничений исходной задачи, а число ограничений в системе двойственной задачи — числу переменных в исходной задаче;
4) коэффициентами при неизвестных в целевой функции двойственной задачи являются свободные члены в системе ограничений исходной задачи, а правыми частями в ограничениях двойственной задачи — коэффициенты при неизвестных в целевой функции исходной задачи;
5) каждому ограничению одной задачи соответствует переменная другой задачи: номер переменной совпадает с номером ограничения; при этом ограничению, записанному в виде неравенства £, соответствует переменная, связанная условием неотрицательности. Если функциональное ограничение исходной задачи является равенством, то соответствующая переменная двойственнвой; задачи может принимать как положительные, так и отрицательные значения
Модель исходной (прямой) задачи в общем виде может быть записана следующим образом:
Модель двойственной задачи имеет вид:
Две приведенные задачи образуют пару симметричных двойственных задач. Основные утверждения о взаимно двойственных задачах содержатся в двух следующих теоремах.
Первая теорема двойственности.
Для взаимно двойственных задач имеет место один из взаимоисключающих случаев:
1. В прямой и двойственной задачах имеются оптимальные решения, при этом значения целевых функций на оптимальных решениях совпадают: f(x) = g(y).
2. В прямой задаче допустимое множество не пусто, а целевая функция на этом множестве не ограниченна сверху. При этом у двойственной задачи будет пустое допустимое множество.
3. В двойственной задаче допустимое множество не пусто, а целевая функция на этом множестве не ограничена снизу. При этом у прямой задачи допустимое множество оказывается пустым.
4.Обе из рассматриваемых задач имеют пустые допустимые множества.
Экономический смысл первой теоремы двойственности следующий. План производства Х и набор оценок ресурсов У оказываются оптимальными тогда и только тогда, когда прибыль от реализации продукции, определенная, при известных заранее ценах продукции равна затратам на ресурсы по «внутренним» (определяемым только из решения задачи) ценам ресурсов yi. Для всех же других планов Х и У обеих задач прибыль от продукции всегда меньше (или равна) стоимости затраченных ресурсов:
f(X) < g(Y}, т. е. ценность всей выпущенной продукции не превосходит суммарной оценки имеющихся ресурсов. Значит величина g(Y) - f(X) характеризует производственные потери в зависимости от рассматриваемой производственной программы и выбранных Оценок ресурсов.
Из первой теоремы двойственности следует, что при оптимальной производственной программе и векторе оценок ресурсов производственные потери равны нулю.
Вторая теорема двойственности
(теорема о дополняющей нежесткости)
Пусть X=(x1,x2,...xn) - допустимое решение прямой задачи, а Y= (y1,y2,...ym) - допустимое решение двойственной задачи. Для того чтобы они были оптимальными решениями соответственно прямой и двойственной задач необходимо и достаточно, чтобы выполнялись следующие соотношения:
*
**
Условия (*) и (**) позволяют, зная оптимальное решение одной из взаимно двойственных задач, найти оптимальное решение другой задачи.
Из второй теоремы двойственности в данном случае следуют такие требования на оптимальную производственную программу X=(X1,X2,...,Xn) и оптимальный вектор оценок Y=(Y1,Y2,...,Ym):
(4)
(5)
Условия (4) можно интерпретировать так: если оценка yi единицы ресурса i-го вида положительна, то при оптимальной производственной программе этот ресурс используется полностью, если же ресурс используется не полностью, то его оценка равна нулю. Из условия (5) следует, что если j-й вид продукции вошел в оптимальный план, то он в оптимальных оценках неубыточен, если же j-й вид продукции убыточен, то он не войдет в план, не будет выпускаться.
Рассмотрим еще одну теорему, выводы которой будут использованы в дальнейшем.
Теорема об оценках.
Значения переменных Yi в оптимальном решении двойственной задачи представляют собой оценки влияния свободных членов bi системы ограничений-неравенств прямой задачи на величину
Решая ЗЛП симплексным методом, мы одновременно решаем двойственную ЗЛП. Переменные двойственной задачи yi называют объективно обусловленными оценками.
Рассмотрим экономическую интерпретацию двойственной задачи на примере задачи оптимального использования ресурсов.
Пример.Сформулируем экономико-математическую модель двойственной задачи к задаче о коврах.
Прямая задача:
f(x) = 3Х1 +4Х2 +3Х3 +1Х4
Ограничения по ресурсам
7Х1 +2Х2 +2Х3 +6Х4 80
5Х1 +8Х2 +4Х3 +3Х4 480
2Х1 +4Х2 +Х3 +8Х4 130
Х1, Х2, Х3, Х4 0
Количество неизвестных в двойственной задаче равно числу функциональных ограничений в исходной задаче. В исходной задаче три ограничения – по труду, по сырью и по оборудованию. Следовательно, в двойственной задаче – три неизвестных:
Y1 – двойственная оценка ресурса труд, или «цена» труда;
Y2 – двойственная оценка ресурса сырье, или «цена» сырья;
Y3 – двойственная оценка ресурса оборудование, или «цена» оборудования.
Целевая функция двойственной задачи формулируется на минимум. коэффициентами при неизвестных в целевой функции двойственной задачи являются свободные члены в системе ограничений исходной задачи.
g = 80 ´Y1 + 480´Y2 + 130´Y3 ® min
Необходимо найти такие “цены” на ресурсы (Yi), чтобы общая стоимость используемых ресурсов была минимальной.
Ограничения. число ограничений в системе двойственной задачи равно числу переменных в исходной задаче. В исходной задаче четыре переменных, следовательно, в двойственной задаче четыре ограничения. правыми частями в ограничениях двойственной задачи являются коэффициенты при неизвестных в целевой функции исходной задачи. Левая часть ограничения определяет стоимость ресурсов, затраченных на производство единицы продукции. Каждое ограничение соответствует определенному виду продукции.
7 ´Y1 + 5´Y2 + 2´Y3 ³ 3
2 ´Y1 + 8´Y2 + 4´Y3 ³ 4
2 ´Y1 + 4´Y2 + 1´Y3 ³ 3
6 ´Y1 + 3´Y2 + 8´Y3 ³ 1
Y1 ,Y2 ,Y3 ³ 0
Задача о размещении производственных заказов
в планируемом периоде необходимо обеспечить производство 300 тыс. однородных новых изделий, которые могут выпускаться на четырех филиалах предприятия. Для освоения этого нового вида изделий выделены капитальные вложения в размере 18 млн. руб.. Разработанные для каждого филиала предприятия проекты освоения нового вида изделия характеризуются величинами удельных капитальных вложений и себестоимостью единицы продукции в соответствии с таблицей.
Необходимо найти такой вариант распределения объемов производства продукции и капитальных вложений по филиалам, при котором суммарная стоимость изделий будет минимальной.
Таблица
Показатель | Филиал предприятия | |||
Себестоимость производства изделия, руб. | ||||
Удельные капиталовложения, руб. |
В результате решения получен план распределения объемов производства по филиалам предприятия
Филиал предприятия | |||
100 тыс. штук | 200 тыс. штук |
Требуется:
1) Сформулировать экономико-математическую модель прямой и двойственной задачи.
2) Найти оптимальный план двойственной задачи, используя теоремы двойственности.
Решение
1) Экономико-математическая модель исходной задачи.
Xi - объем выпускаемой продукции на i-м филиале предприятия.
= 83X1+89X2+95X3+98X4 -> min,
Ограничения
X1+X2+X3+X4 ³ 300 (тыс. штук)
120X1+80X2+50X3+40X4 £ 18 (млн.руб.),
X1,2,3,4 ³0.
Y1 | |||||
Y2 |
Экономико-математическая модель двойственной задачи.
Y1 - двойственная оценка выпускаемой продукции, которая может быть ценой изделия;
Y2 - двойственная оценка капитальных вложений, которая может быть представлена как коэффициент эффективности капитальных вложений.
g =300000 Y1+18000000 Y2 -> mах
1 Y1+120Y2 £83
1 Y1+ 80Y2 £89
1 Y1+ 50Y2 £95
1 Y1+ 40Y2 £98
2) для определения оптимального плана двойственной задачи воспользуемся соотношениями второй теоремы двойственности. Если какое-либо ограничение исходной задачи выполняется как строгое неравенство, то соответствующая двойственная оценка равна нулю
( ).
0+100000+200000+0 = 300000
120´0+80´100000+50´200000+4´0 = 18000000
Если какая-либо переменная исходной задачи входит в оптимальный план, то соответствующее ограничение двойственной задачи выполняется как строгое равенство
).
В нашей задаче Х2=100000>0 и Х3=200000>0, поэтому второе и третье ограничения двойственной задачи обращаются в уравнения, решая которые найдем Y1и Y2 .
1 Y1+ 50Y2 = 95 Y1= 105 - средняя цена изделия
1 Y1+ 80Y2 = 89 Y2 = - 0.2 - двойственная оценка капитальных вложений.
105 =95 +50´0.2 = 105
105 =89+ 80´0.2 = 105
На втором и третьем филиалах выпускать новые изделия целесообразно так как затраты на его освоение и выпуск не превышают цену изделия.
Проверим выполнение первой теоремы двойственности.
g =300000 Y1+18000000 Y2 = 300000 ´105+18000000´(–0.2) = 279 000 000
= 83X1+89X2+95X3+98X4 =83´0+89´100000+95´200000+98´0 = 279 000 000.
Полученные оптимальные планы говорят о том, что в первом и четвертом филиалах размещать заказы по выпуску новых изделий невыгодно (Х1=0 и Х4=0), так как затраты на производство единицы изделия в этих филиалах больше цены изделия.
1 ´Y1+ 120´Y2 =83 Y1= 105 105+ 120´(-0.2) < 95 105< 95+24 = 119
1 ´Y1+ 40´Y2 = 98 Y2 = - 0.2105+ 40´(-0.2) <89 105<98+8 = 106.
Основным литературным источником являются [5,6,7], а дополнительным – [8-11, 13].
Дата добавления: 2016-11-22; просмотров: 1269;