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

51

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 , где