Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

ЭМиКМ (пособие)

.pdf
Скачиваний:
42
Добавлен:
11.04.2015
Размер:
899.68 Кб
Скачать

46

ваний. Обозначим через x1 количество плит, раскраиваемых по первому вари- анту, x2 по второму варианту, x3 по третьему варианту. Составим ограни- чение по выпуску заготовок 1-го типоразмера. Из одной плиты по первому ва- рианту раскроя получаются две таких заготовки, из x1 плит 2x1 заготовок 1-го типоразмера. Кроме того, по одной такой заготовке получится при раскрое ка- ждой плиты по второму варианту. Всего по этому варианту раскраивается x2 плит, из которых вырабатывается x2 заготовок 1-го типоразмера. Таким обра- зом, общее количество заготовок 1-го типоразмера равно 2x1 + x2 и по усло- вию задачи оно не должно быть менее 300, т.е. имеем следующее ограничение:

2x1 + x2 ³ 300.

Аналогичным образом составляется ограничение по выработке заготовок 2-го типоразмера

x1 + 2x2 + 3x3 ³ 400.

Выражение для суммарного количества отходов при раскрое является це- левой функцией, имеющей вид:

Z(X ) =18850x1 +18450x2 +18050x3 .

Наконец, следует учесть естественные ограничения на неотрицательность

переменных

x1 ³ 0, x2 ³ 0, x3 ³ 0.

Таким образом, экономико-математическая модель задачи оптимального раскроя ДСтП на заготовки имеет вид: найти такой план раскроя

X = (x1, x2 , x3 ), при котором целевая функция, определяющая суммарную площадь отходов, минимальна

Z(X ) =18850x1 +18450x2 +18050x3 ® min

и выполняются следующие ограничения:

ì2x1 + x2 ³ 300,

ïíx1 + 2x2 + 3x3 ³ 400, ïîx1 ³ 0, x2 ³ 0, x3 ³ 0.

Задача №10. Определить оптимальный план раскроя 100 листов фанеры размером 2,5×1,5 м на прямоугольные заготовки размерами 2×1 м, 1×0,75 м и 0,5×0,5 м в ассортименте 1:4:12.

47

Общая постановка задачи о раскрое одного материала

На раскрой (распил, обработку) поступает материал одного образца в ко- личестве 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

 

 

 

åxiaik

= 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 принимает максимальное значение.

Общая постановка задачи о раскрое нескольких материалов

Задачу о раскрое можно легко обобщить на случай m раскраиваемых ма-

териалов.

Пусть каждая единица j -го материала ( j = 1, 2,...,m)

может быть

раскроена

n различными способами, причем использование

i -го способа

48

(i =1, 2,...,n) дает aijk единиц k -го изделия (k = 1, 2,...,l), а запас j -го мате- риала равен a j единиц.

Обозначим xij число единиц j -го материала, раскраиваемого i -м спо-

собом. Экономико-математическая модель задачи о раскрое в общей постанов-

ке примет вид: найти такое решение X = (x11, x12,..., xnm ), удовлетворяю-

щее системе

ì n

( j =1,2,...,m)

 

ïåxij £ a j

 

ïi=1

 

(2.21)

í n m

 

ïååxijaijk = bk x

(k =1,2,...,l)

 

ï

 

 

 

 

 

 

îi=1 j =1

 

 

и условию

 

 

 

 

xij ³ 0,

 

(2.22)

при котором целевая функция Z(

 

) принимает максимальное значение:

X

Z(

 

) = x ® max.

(2.23)

X

<Здесь должна быть задача о раскрое нескольких материалов. Например, про- изводственный участок мебельного цеха изготавливает столярные ящики, на которые расходуются доски (для каркаса), ДСтП (для передней стенки) и ДВП (для дна). Придумать!!! >

47

Глава 3. Способы решения задачи линейного программирования

В данной главе на примере решения простой задачи производственного планирования рассмотрены следующие два способа решения ЗЛП: графический и симплекс-метод.

Экономико-математическая модель задачи производственного планирования

Пусть

предприятие обладает

i видами производственных ресурсов

(i =1, 2,...,m), объем каждого из которых обозначим через

bi . Имеющиеся на

предприятии

ресурсы используются

для производства j

видов продукции

( j = 1, 2,...,n), причем известно количество i -го вида ресурсов, затрачиваемое на выпуск единицы j -го вида продукции, которое обозначим через aij . Кроме того, известна прибыль от реализации единицы каждого из j видов выпущен- ной продукции, которую обозначим через c j .

С учетом введенных обозначений экономико-математическая модель за- дачи формирования производственного плана, обеспечивающего получение максимальной прибыли, может быть представлена следующим образом: найти

план выпуска продукции

 

 

 

 

= (x1, x2 ,..., x j ,..., 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

 

x

 

£ b

,

ï

m1 1

 

 

 

 

 

 

 

 

 

 

 

mn n

m

 

ïx

 

³ 0,

x

2

³ 0,

..., x

j

³ 0, ..., x

n

³ 0,

 

 

î 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

при которых целевая функция задачи (прибыль от реализации произведенной продукции) принимает максимальное значение

Z(

 

) = c1x1 + c2 x2 + ...+ c j x j + ...+ cn xn → max .

(3.2)

X

Как отмечалось выше, данная задача относится к ЗЛП, поскольку и целе- вая функция, и ограничения являются линейными.

Выражения (3.1) и (3.2) можно записать в краткой форме:

n

 

 

 

 

 

 

åaij x j £ bi ,

i =

1, m

;

(3.3)

j =1

 

 

 

 

 

 

x j ³ 0,

j =

 

;

(3.4)

1, n

48

n

Z(X ) = åc j x j ® max . (3.5)

j=1

Другим кратким способом записи выражений (3.1) и (3.2) является мат- ричная форма [9]:

 

 

 

 

 

 

 

AX B,

 

 

 

(3.6)

 

 

 

 

 

 

 

X ³ 0 ,

 

 

 

 

(3.7)

 

æa

a

 

... a

... a

CX ® max ,

 

 

 

(3.8)

 

 

ö

 

 

 

 

 

 

ç

11

12

1 j

1n

÷

 

 

 

 

 

где

ça21

a22 ... a2 j ... a2n ÷

 

 

 

 

 

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 x j

+ ...+ 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 +... + c j x j + ... + 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 ø

состоят соответственно из коэффициентов при переменных и свободных чле- нов.

Пример построения экономико-математической модели задачи производственного планирования

49

Задача №10. Предприятие производит два вида продукции P1 и P2 , используя для этого три вида сырья: S1, S2 , S3. Нормы расхода сырья каждого

вида приведены в таблице 3.1. Найти такой план выпуска продукции, при кото- ром прибыль от ее реализации является максимальной.

 

 

 

 

Таблица 3.1

 

Расход сырья на единицу продукции

Прибыль от

 

 

 

 

реализации

Продукция

 

 

 

ед. продукции

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

при ограничениях

 

 

 

 

 

 

 

ìx1 + x2 £ 5,

(3.13)

 

 

ï2x + x

2

£ 9,

(3.14)

ï

 

1

 

 

ï

 

+ 2x2 £ 7,

(3.15)

 

 

íx1

 

 

ïx

 

³ 0,

 

 

(3.16)

ï 1

³ 0,

 

 

(3.17)

 

 

ïx

2

 

 

î

 

 

 

 

где x1 и x2 количество единиц продукции P1 и P2 .

Так как в математической модели задачи использованы только две пере- менные x1 и x2 , ее можно решить графически.

Графический способ решения ЗЛП

Графический (геометрический) способ решения ЗЛП обычно предполага- ет выполнение следующих действий:

1.Запись математических выражений, представляющих целевую функ- цию и ограничения, в виде равенств (уравнений).

2.Построение на графике прямых для уравнений, соответствующих огра- ничениям.

3.Определение области допустимых решений (ОДР) для задачи.

4.Построение на графике прямой, соответствующей целевой функции.

50

5.Параллельный перенос (перемещение) прямой, построенной для целе- вой функции, в одну из крайних точек ОДР для получения оптимально- го решения.

Пример решения задачи производственного планирования графическим методом

Решим задачу производственного планирования, рассмотренную выше, графическим способом. Вначале определим ОДР, в которой одновременно вы- полняются все ограничения, представленные в модели задачи. Для этого пред- ставим каждое из ограничений (3.13)−(3.17) в виде равенства (уравнения) и по- лучим для них по паре точек, через которые проведем на графике отрезки пря- мых.

1-е ограничение (13): x1

+ x2 = 5;

при x1 = 0, x2 = 5;

 

 

при x2 = 0,

x1 = 5, т.е. получим па-

 

 

ру точек с координатами (0; 5) и

 

 

(5; 0).

 

2-е ограничение (14): 2x1 + x2 = 9;

при x1 = 0,

x2 = 9;

 

 

при x2 = 0,

x1 = 4,5, т.е. получим

 

 

пару точек с координатами (0; 9) и

 

 

(4,5; 0).

 

3-е ограничение (15): x1

+ 2x2 = 7 ;

при x1 = 0,

x2 = 3,5;

 

 

при x2 = 0,

x1 = 7 , т.е. получим па-

 

 

ру точек с координатами (0; 3,5) и

 

= 0;

(7; 0).

 

4-е ограничение (16): x1

прямая, совпадающая с осью Ox2 .

5-е ограничение (16): x2 = 0;

прямая, совпадающая с осью Ox1.

Таким образом, ОДР задачи представляет собой выпуклый пятиугольник ABCDE, показанный на рисунке 6.

В качестве уравнения для целевой функции можно использовать выраже- ние 3x1 + 4x2 =12, правая часть которого равняется произведению коэффици- ентов при неизвестных. Полагая x1 = 0, получим x2 = 3. Аналогично, при x2 = 0, получим x1 = 4 . Следовательно, прямая для целевой функции прохо-

дит через пару точек с координатами (0; 3) и (4; 0). На рисунке эта прямая обо- значена пунктиром.

x2

9

8

7

6

5

4 B

3

2

1

A

51

Ñ Z

C

D

E

O

1

2 3

4 5 6

7

8 9 x1

 

Рисунок 6 − Графический способ решения ЗЛП

Направление параллельного переноса прямой, представляющей целевую функцию, определяется вычислением градиента функции Z

 

 

æ

F

 

F

ö

 

ç

;

÷

= (3;4),

x

x

 

ÑZ(X ) = ç

2

÷

è

1

 

 

ø

 

показывающего направление наискорейшего возрастания функции. Вектор, на- чало которого совпадает с началом системы координат, а конец имеет коорди- наты (3; 4), задает направление переноса. Имеется строгое математическое до- казательство того, что оптимальное решение задачи находится в одной из угло- вых (крайних) точек ОДР (см. например, теорему 3.3 в [9]). Перемещая прямую для целевой функции параллельно самой себе в направлении, заданном векто- ром-градиентом, получим, что такой точкой является вершина пятиугольника C, имеющая координаты (3; 2). Таким образом, оптимальный план для данной

задачи составит 3 единицы продукции P1 и 2 единицы продукции P2 , от реали-

зации которых предприятие получит максимум прибыли, равный 17 денежным единицам.

52

Анализ чувствительности модели задачи производственного планирования

После того, как найдено оптимальное решение задачи производственного планирования, можно выполнить анализ его чувствительности к изменениям исходных данных модели. В частности, целесообразно выяснить следующее

[16]:

1.На сколько можно увеличить запас некоторого вида сырья для улучше- ния полученного оптимального плана?

2.На сколько можно снизить запас некоторого вида сырья при сохране- нии полученного оптимального плана?

3.Увеличение запасов какого вида сырья наиболее выгодно?

4.Каков диапазон изменения того или иного коэффициента целевой

функции, при котором не происходит изменение оптимального реше- ния?

Прежде всего, выполним классификацию ограничений, разделив их на

активные (или связывающие), и неактивные (или несвязывающие). Активные ограничения соответствуют прямым, проходящим через точку, представляю- щую оптимальное решение (на рисунке 7 это вершина C пятиугольника ABCDE). Неактивные ограничения соответствуют прямым, которые ограничивают ОДР, но не проходят через точку оптимального решения.

В данном случае активными ограничениями являются (1) и (3), поскольку соответствующие им прямые проходят через точку, координаты которой и дают оптимальный производственный план. Активные ограничения лимитируют за-

пасы сырья S1 и S3, которые логично отнести к разряду дефицитных, тогда

как неактивное ограничение (2) ассоциируется с видом сырья S2 , запас которо- го имеется в некотором избытке, т.е. оно является недефицитным.

Вначале рассмотрим влияние увеличения запасов сырья S1 (ограничение

1) на улучшение полученного оптимального плана. Из рисунка видно, что при увеличении запаса этого вида сырья соответствующая прямая (в частности, от- резок СD) перемещается вверх параллельно самой себе, постепенно стягивая в точку треугольник CKD. В точке K, соответствующей оптимальному решению, активными являются ограничения (2) и (3), а ограничение (1) становится неак-

тивным. Поэтому запас сырья S1 не следует увеличивать сверх того предела,

когда ограничение становится неактивным. Это предельное значение рассчиты- вается путем подстановки координат точки K(11/3; 5/3) в левую часть ограни- чения (1):

1×113 +1× 53 = 163 ,

Z(K) = 3×113 + 4 × 53 = 533 .

 

 

 

 

 

 

 

53

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

9

 

 

S2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

8

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

7

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6

 

 

L

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S1

 

 

 

 

 

A3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ñ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

Рисунок 7 − Анализ чувствительности модели задачи

Аналогично рассматривается вопрос о целесообразности увеличения де- фицитного сырья S3. Из рисунка видно, что при увеличении этого вида сырья

соответствующая прямая (в частности отрезок BC), перемещается вверх парал- лельно самой себе. В точке L ограничение (3) становится неактивным, поэтому

объем сырья S3 не следует увеличивать сверх этого предела. Это предельное

значение рассчитывается путем подстановки координат точки L(0; 5) в левую часть ограничения (3):

1×0 + 2 ×5 =10,

Z (L) = 3×0 + 4 ×5 = 20 .

Таким образом, предельные запасы для сырья S1 и S3 составляют соот-

ветственно 5 13 и 10 единиц соответственно. Дальнейшее повышение этих за-

пасов нецелесообразно, поскольку не приведет к улучшению оптимального ре- шения задачи.

Рассмотрим теперь вопрос об уменьшении запаса недефицитного сырья S2 (ограничение 2). Из рисунка видно, что, не изменяя оптимального решения,

соответствующую прямую (в частности, отрезок KE), можно перемещать па- раллельно самой себе до пересечения с точкой C. Уменьшение запасов сырья