Ekonomiko-matematicheskoe_modelirovanie_v_servise / курс лекций по ЭММ / лекции ЭММ!!!
.pdf51
Z(X ) =18850x1 +18450x2 +18050x3 ® min
и выполняются следующие ограничения:
ì2x1 + x2 ³ 300
ïíx1 + 2x2 + 3x3 ³ 400
ïîx1 ³ 0, x2 ³ 0, x3 ³ 0.
Общая постановка задачи о раскрое одного материала
На раскрой (распил, обработку) поступает материал одного образца в ко- личестве a единиц. Требуется изготовить из него l разных комплектующих из- делий в количествах, пропорциональных b1,b2 ,...,bl (условие комплектности). Каждая единица материала может быть раскроена n различными способами, причем использование i -го способа (i =1,2,...,n) дает aik единиц k -го изделия
(k = 1, 2,...,l).
Необходимо найти план раскроя, обеспечивающий максимальное число комплектов.
Составим экономико-математическую модель задачи. Обозначим xi - число единиц материала, раскраиваемых i -м способом, и x - число изготавли- ваемых комплектов изделий.
Так как общее количество материала равно сумме его единиц, раскраи- ваемых различными способами, то
|
n |
|
|
|
åxi = a. |
|
(2.18) |
|
i=1 |
|
|
Требование комплектности выразится уравнениями |
|
||
n |
|
|
|
åxi aik = bk x |
(k = 1, 2,...,l). |
(2.19) |
|
i=1 |
|
|
|
Очевидно, что |
|
|
|
xi ³ 0 |
(i = 1, 2,...,n). |
(2.20) |
|
Экономико-математическая |
модель |
задачи: найти такое |
решение |
X = (x1, x2 ,..., xn ) , удовлетворяющее системе уравнений (2.18) и (2.19) и условию (2.20), при котором функция Z(X ) = x принимает максимальное значение.
52
Общая постановка задачи о раскрое нескольких материалов
Задачу о раскрое можно легко обобщить на случай m раскраиваемых ма- териалов. Пусть каждая единица j -го материала ( j = 1, 2,...,m) может быть рас-
кроена n |
различными способами, |
причем использование |
i -го способа |
|
(i = 1,2,...,n) |
дает aijk |
единиц k -го изделия (k = 1, 2,...,l), а запас |
j -го материала |
|
равен aj единиц. |
|
|
|
|
Обозначим xij |
- число единиц |
j -го материала, раскраиваемого i -м спо- |
собом. Экономико-математическая модель задачи о раскрое в общей постанов-
ке примет вид: найти такое решение |
|
= (x11, x12 ,..., xnm ) , |
удовлетворяющее |
|||||
X |
||||||||
системе ограничений |
|
|
|
|
|
|
|
|
n |
|
|
|
|
|
|
|
|
ïìåxij £ a j |
|
( j = 1,2,...,m) |
|
|||||
ï i=1 |
|
|
|
|
|
|
|
(2.21) |
í n |
m |
|
|
|||||
ïïååxij aijk = bk x |
|
(k =1,2,...,l) |
|
|||||
î i=1 |
j=1 |
|
|
|
||||
и условию |
|
|
|
|
|
|
|
|
|
|
|
|
|
xij ³ 0, |
|
|
(2.22) |
при которых целевая функция Z( |
|
|
) принимает максимальное значение: |
|||||
X |
||||||||
|
Z ( |
|
) = x ® max . |
(2.23) |
||||
|
X |
53
Глава 3. Способы решения задачи линейного программирования
В данной главе на примере решения простой задачи производственного планирования рассмотрены следующие два способа решения ЗЛП: графический и симплекс-метод.
Экономико-математическая модель задачи производственного планирования
Пусть |
предприятие обладает |
i видами производственных ресурсов |
||
(i = 1, 2,...,m) , |
объем каждого из которых обозначим через |
bi . Имеющиеся на |
||
предприятии |
ресурсы используются |
для |
производства j |
видов продукции |
( j = 1, 2,...,n) , |
причем известно количество |
i -го вида ресурсов, затрачиваемое |
на выпуск единицы j -го вида продукции, которое обозначим через aij . Кроме того, известна прибыль от реализации единицы каждого из j видов выпущен- ной продукции, которую обозначим через cj .
С учетом введенных обозначений экономико-математическая модель за- дачи формирования производственного плана, обеспечивающего получение максимальной прибыли, может быть представлена следующим образом: найти
план выпуска продукции |
|
|
= (x1, x2 ,..., xj ,..., xn ), удовлетворяющий системе ог- |
||||||||||||||||||||||||
X |
|||||||||||||||||||||||||||
раничений |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ìa11x1 + a12 x2 |
+ ... + a1 j x j |
+ ... + a1n xn £ b1 |
|
||||||||||||||||||||||||
ïa |
x |
+ a |
22 |
x |
2 |
|
+ ... + a |
2 j |
x |
j |
|
+ ... + a |
2n |
x |
n |
£ b |
|
||||||||||
ï |
21 1 |
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
||||||||||||
ï |
|
|
|
... |
|
|
|
|
|
... ... |
|
|
|
|
... ... |
|
|
|
|
|
... |
(3.1) |
|||||
í ... |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
ïa |
|
x |
+ a |
m2 |
x |
2 |
+ ... + a |
mj |
x |
j |
+ ... + a |
mn |
x |
n |
£ b |
||||||||||||
ï |
m1 1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
m |
||||||||||||
ï |
|
|
|
|
|
|
³ 0, ..., x j |
³ 0, ..., xn |
³ 0, |
|
|
||||||||||||||||
îx1 ³ 0, x2 |
|
|
при которых целевая функция задачи (прибыль от реализации произведенной продукции) принимает максимальное значение
Z( |
|
) = c1x1 + c2 x2 + ...+ cj xj + ...+ cn xn → max. |
(3.2) |
X |
Как отмечалось выше, данная задача относится к ЗЛП, поскольку и целе- вая функция, и ограничения являются линейными.
54
Выражения (3.1) и (3.2) можно записать в краткой форме
|
|
|
|
|
|
|
|
|
n |
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
åaij x j £ bi , |
i = |
|
|
; |
(3.3) |
|||
|
|
|
|
|
|
|
|
|
1, m |
||||||||
|
|
|
|
|
|
|
|
|
j=1 |
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
x j ³ 0, |
j = |
|
; |
(3.4) |
||||
|
|
|
|
|
|
|
|
|
1,n |
||||||||
|
|
|
|
|
|
|
|
|
|
|
n |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Z( |
|
) = åc j xj ® max . |
(3.5) |
|||||
|
|
|
|
|
|
|
|
|
X |
||||||||
|
|
|
|
|
|
|
|
|
|
|
j=1 |
|
|
|
|
|
|
Другим кратким способом записи выражений (3.1) и (3.2) является сле- |
|||||||||||||||||
дующая матричная форма [9]: |
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
AX ≤ B , |
|
|
|
|
|
(3.6) |
|
|
|
|
|
|
|
|
|
|
|
X ³ 0, |
|
|
|
|
|
(3.7) |
|
|
|
|
|
|
|
|
|
|
CX ® max , |
(3.8) |
||||||
æa |
|
a |
|
... a |
|
... a |
|
ö |
|
|
|
|
|
|
|
|
|
ç 11 |
12 |
1 j |
1n |
÷ |
|
|
|
|
|
|
|
|
|
||||
ça |
21 |
a |
22 |
... a |
2 j |
... a |
2n |
÷ |
- так называемая технологическая матрица; |
|
|||||||
где A = ç |
|
|
|
÷ |
|
||||||||||||
ç ... ... ... ... ... ... |
÷ |
|
|
|
|
|
|
|
|
|
|||||||
ç |
|
|
|
|
|
|
|
÷ |
|
|
|
|
|
|
|
|
|
èam1 am2 ...amj ...amn |
ø |
|
|
|
|
|
|
|
|
|
æb |
ö |
|
æ x |
|
ö |
|
ç 1 |
÷ |
|
ç 1 |
÷ |
|
|
çb2 |
÷ |
- вектор ресурсных ограничений; |
ç x2 |
÷ |
- вектор-план; |
|
B = ç |
÷ |
X = ç |
|
÷ |
||
ç ... |
÷ |
|
ç ... |
÷ |
|
|
çb |
÷ |
|
ç x |
n |
÷ |
|
è m |
ø |
|
è |
ø |
|
С = (с1 с2 ... cn ) - вектор стоимости.
Разновидностью матричной формы записи математической модели зада- чи является векторная форма [20]
P1x1 + P2 x2 + ...+ Pj xj + ...+ Pn xn ≤ P , |
(3.9) |
X ³ 0, |
(3.10) |
CX ® max , |
(3.11) |
где CX - скалярное произведение векторов С = (с1 с2 ... cn ) и |
X = (x1 x2 ... xn ) , |
т.е. число, равное сумме произведений соответствующих координат этих век- торов: CX = c1x1 + c2 x2 + ...+ cj xj + ...+ cn xn (по определению); векторы
|
æ a |
|
ö |
|
|
æa |
|
ö |
|
|
æa |
|
ö |
|
|
æb |
ö |
|
ç 11 |
÷ |
|
|
ç 12 |
÷ |
|
|
ç 1n |
÷ |
|
|
ç 1 |
÷ |
|||
P = |
ç a21 ÷ |
, |
P = |
ça22 ÷ |
, …, |
P = |
ça2n ÷ |
, |
P = |
çb2 ÷ |
|||||||
ç |
|
÷ |
ç |
|
÷ |
ç |
|
÷ |
ç |
÷ |
|||||||
1 |
|
|
2 |
|
|
n |
|
|
|
||||||||
|
ç... |
÷ |
|
|
ç... |
÷ |
|
|
ç... |
÷ |
|
|
ç... |
÷ |
|||
|
ç a |
m1 |
÷ |
|
|
ça |
m2 |
÷ |
|
|
ça |
|
÷ |
|
|
çb |
÷ |
|
è |
ø |
|
|
è |
ø |
|
|
è |
mn ø |
|
|
è m ø |
55
состоят соответственно из коэффициентов при переменных и свободных чле- нов.
Пример построения экономико-математической модели задачи производственного планирования
Задача № 10. Предприятие производит два вида продукции P1 и P2 , ис- пользуя для этого три вида сырья: S1 , S2 , S3 . Нормы расхода сырья каждого вида приведены в табл. 8. Найти такой план выпуска продукции, при котором прибыль от ее реализации является максимальной.
|
|
|
|
Таблица 8 |
|
Расход сырья на единицу продукции |
Прибыль от |
||
|
|
|
|
реализации |
Продукция |
|
|
|
ед. продукции |
S1 |
S2 |
S3 |
||
P1 |
1 |
2 |
1 |
3 |
P2 |
1 |
1 |
2 |
4 |
Кол-во сырья |
5 |
9 |
7 |
|
в наличии |
|
|
|
|
|
|
|
|
|
Экономико-математическую модель данной задачи можно записать сле-
дующим образом: определить объемы производства продукции P1 и P2 , при
которых достигается максимизация целевой функции |
|
|||||
Z( |
|
) = 3x1 + 4x2 ® max |
(3.12) |
|||
X |
||||||
при ограничениях |
|
|
|
|
||
|
|
ìx |
+ x |
|
£ 5 |
(3.13) |
ï 1 |
|
2 |
|
(3.14) |
||
|
|
ï2x1 + x2 £ 9 |
||||
ï |
+ 2x2 £ 7 |
(3.15) |
||||
|
|
ïíx1 |
||||
|
|
ïx1 ³ 0 |
|
|
(3.16) |
|
ï |
³ 0, |
|
(3.17) |
|||
|
|
îx2 |
|
где x1 и x2 - количество единиц продукции P1 и P2 .
Так как в математической модели задачи использованы только две пере- менные x1 и x2 , ее можно решить графически.
56
Графический способ решения ЗЛП
Графический (геометрический) способ решения ЗЛП обычно предполага- ет последовательное выполнение следующих действий:
1.Запись математических выражений, представляющих целевую функ- цию и ограничения, в виде равенств (уравнений).
2.Построение на графике прямых для уравнений, соответствующих огра- ничениям.
3.Определение области допустимых решений (ОДР) для задачи.
4.Построение на графике прямой, соответствующей целевой функции.
5.Параллельный перенос (перемещение) прямой, построенной для целе- вой функции, в одну из крайних точек ОДР для получения оптимально- го решения.
Пример решения задачи производственного планирования графическим методом
Решим задачу производственного планирования, рассмотренную выше, графическим способом. Вначале определим ОДР, в которой одновременно вы- полняются все ограничения, представленные в модели задачи. Для этого пред- ставим каждое из ограничений (3.13)−(3.17) в виде равенства (уравнения) и по- лучим для них по паре точек, через которые проведем на графике отрезки пря- мых.
1-е ограничение (3.13): x1 |
+ x2 = 5 ; при |
x1 = 0 , x2 = 5 ; |
|
|
|
при |
x2 = 0 , x1 = 5, т.е. получим пару |
|
|
точек с координатами (0; 5) и (5; 0). |
|
2-е ограничение (3.14): 2x1 + x2 |
= 9 ; при |
x1 = 0 , x2 = 9 ; |
|
|
|
при x2 = 0 , x1 = 4,5 , т.е. получим па- |
|
|
|
ру точек с координатами (0; 9) и |
|
|
|
(4,5; 0). |
|
3-е ограничение (3.15): x1 |
+ 2x2 |
= 7; при |
x1 = 0 , x2 = 3,5; |
|
|
при x2 = 0 , x1 = 7 , т.е. получим пару |
|
|
|
точек с координатами (0; 3,5) и (7; 0). |
|
4-е ограничение (3.16): x1 |
= 0 ; |
прямая, совпадающая с осью Ox2 . |
|
5-е ограничение (3.17): x2 |
= 0 ; |
прямая, совпадающая с осью Ox1 . |
57
Таким образом, ОДР задачи представляет собой выпуклый пятиугольник ABCDE, показанный на рис. 7.
x2 |
|
|
|
|
9 |
|
|
|
|
8 |
|
|
|
|
7 |
|
|
|
|
6 |
|
|
|
|
5 |
|
|
|
|
4 |
B |
Ñ Z |
|
|
3 |
|
C |
|
|
2 |
|
|
|
|
|
|
|
|
|
1 |
|
D |
|
|
|
E |
|
|
|
A |
|
|
|
|
O |
1 |
2 3 4 5 6 |
7 |
8 9 x1 |
Рис. 7. Графический способ решения ЗЛП
В качестве уравнения для целевой функции можно использовать выраже- ние 3x1 + 4x2 = 12, правая часть которого равняется произведению коэффициен- тов при неизвестных. Полагая x1 = 0 , получим x2 = 3. Аналогично, при x2 = 0 , получим x1 = 4 . Следовательно, прямая для целевой функции проходит через пару точек с координатами (0; 3) и (4; 0). На рисунке эта прямая обозначена пунктиром.
Направление параллельного переноса прямой, представляющей целевую функцию, определяется вычислением градиента функции Z :
|
|
æ |
¶F |
|
¶F |
ö |
|
||
|
|
|
|||||||
ç |
; |
÷ |
= (3;4) , |
||||||
¶x |
¶x |
|
|||||||
ÑZ(X ) = ç |
2 |
÷ |
|||||||
è |
1 |
|
|
ø |
|
показывающего направление наискорейшего возрастания функции. Вектор, на- чало которого совпадает с началом системы координат, а конец имеет коорди-
58
наты (3; 4), задает направление переноса. Имеется строгое математическое до- казательство того, что оптимальное решение задачи находится в одной из угло- вых (крайних) точек ОДР (см., например, теорему 3.3 в [9]). Перемещая прямую для целевой функции параллельно самой себе в направлении, заданном векто- ром-градиентом, получим, что такой точкой является вершина пятиугольника C, имеющая координаты (3; 2). Таким образом, оптимальный план для данной задачи составит 3 единицы продукции P1 и 2 единицы продукции P2 , от реали- зации которых предприятие получит максимум прибыли, равный 17 денежным единицам.
Анализ чувствительности модели задачи производственного планирования
После того, как найдено оптимальное решение задачи производственного планирования, можно выполнить анализ его чувствительности к изменениям исходных данных модели. В частности, целесообразно выяснить следующее
[16]:
1.На сколько можно увеличить запас некоторого вида сырья для улучше- ния полученного оптимального плана?
2.На сколько можно снизить запас некоторого вида сырья при сохране- нии полученного оптимального плана?
3.Увеличение запасов какого вида сырья наиболее выгодно?
4.Каков диапазон изменения того или иного коэффициента целевой функции, при котором не происходит изменение оптимального реше- ния?
Прежде всего выполним классификацию ограничений, разделив их на
активные (или связывающие) и неактивные (или несвязывающие). Активные ограничения соответствуют прямым, проходящим через точку, представляю- щую оптимальное решение (на рис. 8 это вершина C пятиугольника ABCDE). Неактивные ограничения соответствуют прямым, которые ограничивают ОДР, но не проходят через точку оптимального решения.
В данном случае активными являются 1-е и 3-е ограничения, поскольку соответствующие им прямые проходят через точку, координаты которой и дают оптимальный производственный план. Активные ограничения лимитируют за- пасы сырья S1 и S3 , которые логично отнести к разряду дефицитных, тогда как
59
неактивное 2-е ограничение ассоциируется с видом сырья S2 , запас которого имеется в некотором избытке, т.е. оно является недефицитным.
x2 |
|
|
|
|
|
|
9 |
S2 |
|
|
|
|
|
8 |
|
|
|
|
|
|
|
|
|
|
|
|
|
7 |
|
|
|
|
|
|
6 |
L |
|
|
|
|
|
5 |
|
A3 |
|
|
|
|
S1 |
|
ÑZ(X ) |
|
|
||
|
|
|
|
|
||
4 |
B |
|
|
A1 |
|
|
3 |
|
|
|
|
|
|
|
C |
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
K |
|
|
|
1 |
|
|
|
D |
|
S3 |
|
|
|
E |
|
||
A |
|
|
|
|
||
|
|
|
|
|
||
O |
1 |
2 3 |
4 5 6 |
7 |
8 9 x1 |
Рис. 8. Анализ чувствительности модели задачи
Вначале рассмотрим влияние увеличения запасов сырья S1 (1-е ограниче- ние) на улучшение полученного оптимального плана. Из рисунка видно, что при увеличении запаса этого вида сырья соответствующая прямая (в частности, отрезок СD) перемещается вверх параллельно самой себе, постепенно стягивая в точку треугольник CKD. В точке K, соответствующей оптимальному реше- нию, активными являются 2-е и 3-е ограничения, а 1-е ограничение становится неактивным. Поэтому запас сырья S1 не следует увеличивать сверх того преде- ла, когда ограничение становится неактивным. Это предельное значение рас- считывается путем подстановки координат точки K(11/3; 5/3) в левую часть 1-го ограничения:
1×113 +1× 53 = 163 ,
60
Z (K ) = 3×113 + 4 × 53 = 533 .
Аналогично рассматривается вопрос о целесообразности увеличения де- фицитного сырья S3 . Из рисунка видно, что при увеличении этого вида сырья соответствующая прямая (в частности отрезок BC), перемещается вверх парал- лельно самой себе. В точке L 3-е ограничение становится неактивным, поэтому объем сырья S3 не следует увеличивать сверх этого предела. Это предельное значение рассчитывается путем подстановки координат точки L(0; 5) в левую часть 3-го ограничения:
1×0 + 2 ×5 =10,
Z (L) = 3× 0 + 4 ×5 = 20 .
Таким образом, предельные запасы для сырья S1 и S3 составляют соот-
ветственно 513 и 10 единиц соответственно. Дальнейшее повышение этих запа-
сов нецелесообразно, поскольку не приведет к улучшению оптимального реше- ния задачи.
Рассмотрим теперь вопрос об уменьшении запаса недефицитного сырья S2 (2-е ограничение). Из рисунка видно, что, не изменяя оптимального реше- ния, соответствующую прямую (в частности, отрезок KE) можно перемещать параллельно самой себе до пересечения с точкой C. Уменьшение запасов сырья S2 до величины 2 ×3 +1× 2 = 8 никак не повлияет на оптимальность полученного ранее решения.
Результаты проведенного анализа представлены в табл. 9 [16].
|
|
|
|
|
Таблица 9 |
Сырье |
Тип сырья |
Макс. изменение запаса, |
Макс. изменение прибыли, |
||
|
|
Si |
DZ( |
|
) |
|
|
X |
|||
S1 |
Дефицитное |
16 / 3 - 5 =1/ 3 |
53/ 3 - 51/ 3 = 2 / 3 |
||
S2 |
Недефицитное |
9 -8 =1 |
17 -17 = 0 |
||
S3 |
Дефицитное |
10 - 7 = 3 |
20 -17 = 3 |
Для ответа на третий вопрос введем характеристики ценности каждой до- полнительной единицы сырья. Обозначим ценность дополнительной единицы i -го вида сырья через Yi , где