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

pizhurin_a_a_pizhurin_a_a_modelirovanie_i_optimizaciya_proce

.pdf
Скачиваний:
271
Добавлен:
27.03.2016
Размер:
14.94 Mб
Скачать

дукции выполнен, ресурсы по всем видам оборудования не превышены, а полученная прибыль максимальна.

Составим математическую модель этой задачи. Обозначим через х\ количество изготавливаемых фабрикой шкафов; * 2 - столов; * 3 - стульев. Тогда требование выполнения плана запишется в виде неравенств:

*i> 150; х2 >200; *3>400.

(3.1)

Неравенства (3.1) представляют собой ограничения, обеспечиваю­ щие выполнение планового задания. Найдем выражение для длительности работы фрезерных станков. В процессе изготовления одного шкафа такой станок работает 0,25 ч (см. табл. 3.1). Длительность работы этих станков при изготовлении jci шкафов равна 0,25*1 ч. Аналогично, при изготовлении всех столов затрачивается 0 ,20*2 ч; всех стульев - 0,30*з ч машинного вре­ мени фрезерных станков. Всего затраты машинного времени по фрезерным станкам составляют 0,25*1 + 0 ,2 0 *2 + 0,30*3 ч. Поскольку суммарное время работы этого оборудования не должно превышать 250 ч, то ограничение по ресурсам машинного времени для фрезерных станков следует записать в виде неравенства

0,25*! + 0,20*2 + 0,30*з £ 250.

(3.2)

Подобным же образом легко составить аналогичные ограничения для сверлильных и шлифовальных станков. Они имеют вид соответственно

0,18*!+ 0,13*2+ 0,11*з ^

300;

(3.3)

0,24*! + 0,19*2 + 0,14*з ^

320.

(3.4)

Рассмотрим зависимость критерия оптимальности от переменных математической модели. Эта зависимость носит название целевой функ­ ции. В нашем примере критерием является прибыль.

Величина прибыли от изготовления х\ шкафов равна 5*i руб.; от из­ готовления * 2 столов - 3*2 руб.; от изготовления * 3 стульев - 2* 3 руб. Целе­ вая функция, таким образом, имеет вид W = 5*j + 3* 2 + 2*з. Требование достижения ее максимума записывают в виде

5*i + 3*2 + 2*з —>шах.

(3.5)

Таким образом, математическая модель данной задачи состоит из целевой функции (3.5) и ограничений (3.1) - (3.4) (ее решение приведено в п. 3.6). Существенно, что как целевая функция, так и левые части ограни­ чений являются линейными функциями переменных *ь *2, *з, то есть мно­ гочленами 1-го порядка. Поэтому составленная модель относится к классу задач линейного программирования.

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

требуется найти значения переменных х\, х2, х3, при которых функция (3.5) обращается в максимум и выполняются неравенства (3.1) - (3.4).

3.1.2. Задача оптимального раскроя древесностружечных плит (ДСтП)

Плиты размером 350 х 175 см подлежат раскрою на заготовки двух типоразмеров: 200x70 см и 160x90 см. Требуется получить не менее 300 заготовок первого и не менее 400 заготовок второго типоразмера. При этом суммарное по площади количество отходов должно быть минимально.

Рассмотрим все возможные варианты раскроя. На рис. 3.1, а пока­ зан вариант раскроя плиты на две заготовки 1-го и одну заготовку 2 -го ти­ поразмера, площадь отходов равна 18 850 см2. Часть плиты, уходящая в от­ ходы, заштрихована. Все другие варианты, содержащие эти же три заго­ товки, различаются только их расположением на плите и эквивалентны с точки зрения экономичности.

По варианту раскроя на рис. 3.1, б можно получить одну заготовку 1-го и две заготовки 2-го типоразмера, площадь отходов равна 18 450 см2. По третьему варианту раскроя получают три заготовки 2-го типоразмера (рис. 3.1, в) с площадью отходов 18050 см2.

Рис. 3.1. Варианты раскроя древесностружечной плиты

Для решения задачи следует выяснить, сколько плит надо раскро­ ить по каждому из рассмотренных вариантов при выполнении предъявляв-

мых требований. Обозначим через х\ количество плит, раскраиваемых по первому варианту; через х2 - по второму и через х3 - по третьему. Соста­ вим ограничение по выпуску заготовок 1-го типоразмера. Из одной плиты по первому варианту раскроя получаются две такие заготовки, из плит 2х) заготовок 1-го типоразмера. Кроме того, по одной такой заготовке по­ лучится при раскрое каждой плиты по второму варианту. Всего по этому варианту раскраивается х2 плит, из которых вырабатывается, следователь­ но, х2 заготовок 1-го типоразмера. По третьему варианту такие заготовки не получают. Значит, общее количество равно 2xi + х2. Поэтому имеем ог­ раничение

2*1+ *2 >300.

(3.6)

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

х\ + 2х 2 + Злг3 > 400.

(3.7)

Выражение для суммарного количества отходов при раскрое явля­ ется минимизируемой целевой функцией и имеет вид

18850*! + 18450*2 + 18050*3

min.

(3.8)

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

> 0; х2 > 0; хз > 0.

(3.9)

Совокупность соотношений (3.6) - (3.9) представляет собой мате­ матическую модель данной задачи, которая как и в предыдущем случае, является примером ЗЛП.

3.1.3. Задача о рациональном использовании сырья

Имеется сырье 1-го и 2-го сорта в количестве соответственно Ь\ и Ь2, м3. Три предприятия по его переработке, выпускающие однородную продукцию, располагают свободными мощностями, позволяющими пере­ работать соответственно С], с2 и с3, м3, сырья. Эти предприятия работают в разных условиях и имеют соответственно различные нормы выхода про­ дукции.

Пусть atj - объем продукции, вырабатываемой из 1 м3 сырья /-го сорта наj -м предприятии, i = 1, 2 ;/ = 1, 2, 3. Требуется определить, в каком количестве необходимо поставлять сырье каждого сорта на каждое из предприятий. При этом должен быть обеспечен максимальный объем вы­ пуска продукции, и учтены ограничения по запасам сырья и объемам его переработки на предприятиях. Элементы решения в данной задаче удобно обозначить переменными с двумя индексами: первый индекс - сорт сырья; второй - номер предприятия.

44

Таким образом, через Ху обозначен объем сырья /-го сорта, постав­ ляемого на у'-е предприятие. Задача поэтому содержит шесть элементов решения: х\\, xi2, *13, х2и *22, *23- Общий объем сырья 1-го сорта, поставляе­ мого на все три предприятия, равен сумме [х1х+ хп + х13) и не должен пре­

вышать имеющегося запаса Ь\9 м3. Следовательно, имеем группу из двух ограничений по запасам сырья каждого сорта:

*11 + *12 + *13 ^ * Ь 1

Г

(3.10)

*21 + *22 + *23 ^ ^ 2 -J

Объем поставок сырья 1-го и 2-го сорта на первое предприятие, равный (*ц+ х2\), не должен превышать мощность с\ этого предприятия по его переработке. Исходя из этого, составляются ограничения по мощности переработки сырья для первого и остальных предприятий:

*11 + *21 ^ с\\

*12 + *22 ^ С2;

(3 .11)

*13 + *23 ^ С3. ^

Общий объем выработки продукции на всех предприятиях равен

сумме

2 з

#11*11 + #12*12 + #13*13 + #21*21 + #22*22 + #23*23 =

М У=1

для которой надо найти максимум. Это выражение служит целевой функ­ цией задачи

^ = I Z * i7*i7 - ^ ах-

(3.12)

/=1 ;=1

Кроме того, надо учесть требование неотрицательности перемен­

ных

х,7 > 0 , i = l,2; j = 1,2,3.

(3.13)

Построенная математическая модель (3.10) - (3.13) - это еще один пример ЗЛП.

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

Сформулируем задачу линейного программирования в общем виде. Имеется п переменных хи *2»--->*«. Требуется найти такие их значения, при которых целевая функция

W~C\X\ + с2х 2 + ... + спхп

(3.14)

обращается в максимум (минимум). При этом переменные хи х2 хп должны удовлетворять ряду ограничений, каждое из которых должно от­ носиться к одному из следующих типов:

dix1+ d 2x2 +... + dnxn= d 0;

(3.15)

ei*i + е 2д:2 + ... + е„х„ > е 0;

(3.16)

f lx[ + f 2x2 +... + f nx>l< f 0.

(3.17)

Кроме того, имеются ограничения вида

 

*!><); х2 > 0 ; хп > 0 ,

(3.18)

которые называют обычно простыми или тривиальными ограничения­ ми. Очевидно, что ограничения вида (3.17) легко сводятся к виду (3.16) умножением обеих частей исходного неравенства на ( - 1) и обратно - вид (3.16) к виду (3.17). От ограничений-равенств тоже всегда можно перейти к эквивалентному неравенству вида (3.16) посредством исключения какойлибо переменной.

Пусть, например, в ЗЛП имеется ограничение 2х\ + Зх2 - 3 = 4, где переменные х\9х2, х3неотрицательны. Выразим из этого равенства х3:

х з (2/5)дс, + (3 /5) х 2 - 4 /5 .

(3.19)

Тогда равенство (3.19) в сочетании с ограничением х3 >

0 эквива­

лентно неравенству (2/5)х} + (3/5)х2 - 4 / 5 > 0 . Отметим, наконец, что зада­

ча, в которой требуется найти максимум целевой функции W, эквивалентна задаче минимизации целевой функции W' = - W .

Запишем задачу линейного программирования общего вида в сле­ дующей форме:

W =с0 + с}хj + с2х2 +... + спхп min;

(3.20)

b l+ b ux ]+ b nx2 +... + buxn>0;

(3.21)

^!> 0 ; х 2>0; ... ; ^л > 0 ,

(3.22)

где т - число нетривиальных ограничений.

 

Совокупность значений переменных х\, х2

хп, удовлетворяющая

всем ограничениям ЗЛП, называется ее допустимым решением. Оптимальным решением ЗЛП называется такое ее допустимое решение, при

котором целевая функция достигает экстремума, т. е. максимума или ми­ нимума, в зависимости от условий задачи.

Пример.

W = 3xj+ 2 х 2 ->• max;

(3.23)

2 х 1+ 1 х 2 <14;

(3.24)

4jc,+ *2 < 4 ;

(3.25)

jc, > 0; х2 > 0 .

(3.26)

Возьмем хх= 0,5; x2 = 1. Подставив эти значения в неравенства (3.24) - (3.26), убеждаемся, что все они выполняются. Следовательно, совокупность значений х\ = 0,5; Х2 = 1 является допустимым решением данной ЗЛП. Значение целевой функции при

этом равно W = 3x0,5 + 2x1 = 3,5. Можно доказать, что оптимальным решением этой

задачи будет другое допустимое решение: х* =7/13; х* - 24/13, при котором значение

целевой функции оказывается равным 69/13.

3.2. Геометрическая интерпретация задачи линейного программирования

Для выяснения основных свойств решения задачи линейного про­ граммирования проведем ее геометрическую интерпретацию.

Требуется найти максимум функции

W = 3x j + 4 x 2

(3.27)

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

 

6 *!+ 5х 2 < 30;

(3.28)

4х,+ 1х2 < 28;

(3.29)

х,>0; х2 >0.

(3.30)

Рассмотрим координатную плоскость на рис». 3.2, по оси абсцисс которой будем откладывать значения переменной X], а по оси ординат - значения х2. Каждой точке плоскости соответствует некоторая пара значе­ ний *1 и х2. Множество точек, для которых выполняется ограничение х ,> 0 , представляет собой правую полуплоскость, т. е. часть плоскости

правее оси 0Х2, а все точки, для которых х2 > 0, располагаются в верхней

полуплоскости - выше оси 0Хх. Точки, для которых удовлетворяются оба эти неравенства, jtj>0 и х2 >0, лежат, очевидно, в первом квадранте. До­

бавим к этим ограничениям неравенство (3.28). Рассмотрим сначала соот­ ветствующее равенство:

6 х х+ 5 х 2 =ЪЪ.

(3.31)

Это уравнение прямой /, отсекающей отрезок в 5 единиц на оси 0Х\ и в 6 единиц на оси 0Х2. Учитывая знак неравенства (3.28), можно сделать

вывод о том, что точки,

 

 

удовлетворяющие

ему,

 

 

лежат ниже прямой / или

 

 

на ней самой. Штрихов­

 

 

кой возле прямой отме­

 

 

чена область, в которой

 

 

соответствующее

нера­

 

 

венство не

выполняется.

 

 

Одновременно трем

ог­

 

 

раничениям, х х> 0 ; х2> О

 

 

и (3.28), удовлетворяют

 

 

точки, лежащие внутри и

 

 

на границе треугольника

 

 

ОСВ. Аналогично учи­

 

 

тывается

ограничение

 

 

(3.29). Все точки, для ко­

 

 

торых

оно

выполняется,

Рис. 3.2. Геометрическая интерпретация

 

лежат

ниже прямой

т

 

или на ней. Уравнение

ЗЛП, имеющей единственное решение

 

 

 

этой прямой

 

 

 

 

 

 

 

 

 

 

4х,+ 1х2 =28.

(3.32)

Поэтому точки, удовлетворяющие всем четырем ограничениям задачи, ле­ жат внутри и на границах выпуклого четырехугольника ОМКВ (на рис. 3.2 заштрихован). Здесь К - точка пересечения прямых / и т.

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

Среди множества допустимых решений следует найти оптимальное. Для этого введем в рассмотрение целевую функцию (3.27). Будем прида­ вать ей различные значения и исследовать, где находятся соответствующие им точки в ОДР. Положим сначала значение целевой функции равным ну­ лю: W = Зх,+ 4х2 = 0. Это уравнение прямой п на рис. 3.2. На ней имеется

единственная точка {*]= 0; х2 = О}, принадлежащая ОДР. Пусть теперь

W = Ъхх+ 4лг2 = 6 . Это уравнение прямой /?, параллельной предыдущей и

расположенной выше нее. Стрелкой указано направление смещения пря­ мой, соответствующее возрастанию целевой функции. Среди всех точек этой прямой области допустимых решений принадлежат только точки от­ резка AD, для каждой из которых W = 6. Еще предпочтительнее с точки зрения приближения к максимуму целевой функции положение отрезка EF прямой q, уравнение которой 3*j+4* 2 =12. Очевидно теперь, что макси­

мальное значение целевой функции на множестве допустимых решений достигается при самом верхнем положении прямой W, когда она все еще проходит хотя бы через одну точку, принадлежащую ОДР. Это прямая г, проходящая через точку К. Ее координаты определяются из решения сис-

темы уравнений (3.31) и (3.32)

2

2

х<= 3—;

х2 = 2—. Эта пара чисел и есть оп-

 

п

и

тимальное решение данной ЗЛП. Значение целевой функции для него рав­ но

W = 3-3—+ 4 - 2 —= 18—.

и и и

В приведенном примере ЗЛП имела единственное решение. Пока­ жем, что возможны и другие случаи.

Рассмотрим ЗЛП

W = 5хх+Ъх2 ->max;

(3.33)

4х,+ 8*2 <32;

(3.34)

5x,+ 3x2 <15;

(3.35)

Xj>0; х 2>0 .

(3.36)

Здесь область допустимых решений - многоугольник ОАВС (рис. 3.3, а). Отрезок ВС является частью прямой /, уравнение которой получено из неравенства (3.35) и имеет вид 5хх+Ъх2 =15. Все прямые, на которых

целевая функция (3.33) принимает фиксированное значение W = 5xj+ + 3;t2 = const, будут параллельны этой прямой. Поэтому максимум ее дос­

тигается в любой точке отрезка ВС. Таким образом, эта ЗЛП имеет бесчис­ ленное множество решений.

Возможны случаи, когда задача линейного программирования не имеет решений. Это возможно, когда множество допустимых решений пусто. Например, в ЗЛП

9 х - 2 -» min;

(3.37)

х 1+2x2 <2;

(3.38)

2хх+ Зх2 > 6 ;

(3.39)

*!>(); * 2 - 0

(3.40)

система неравенств (3.38) - (3.40) не имеет решений, поэтому допустимых решений задачи не существует (см. рис. 3.3, б). Кроме того, такой случай возможен, если ОДР существует, но не ограничена, например, в ЗЛП:

W= х }+х 2 —> ш ах;

х}- 2 х 2 < 1;

49

-2х}+ х2 < 2 ;

*j> 0 ; *2 > 0 .

Вэтом случае целевая функция может принимать сколь угодно большие значения (см. рис. 3.3, в).

Рис 3.3. Графическое представление различных ЗЛП:

а - имеющей множество решений; б - имеющей пустое множество допустимых решений; в - с неограниченной ОДР; г - область допустимых решений

3.3. Основная задача линейного программирования и свойства ее решений

Рассмотрим так называемую основную задачу линейного програм­ мирования (ОЗЛП), отличающуюся от задачи линейного программирова­ ния общего вида (3.20) - (3,22). В ней требуется минимизировать целевую

функцию вида (3.14), а все нетривиальные ограничения имеют вид ра­ венств. Общий вид 03 ЛП следующий:

W = с0 +

CJXJ+ с

2 +... + спхп -> min;

(3.41)

^ 21*

1+ Ь22х2

+••• + £ 2пХп ~ ^2 >

 

 

 

>

(3.42)

* i> 0 ; х2 > 0 ; ... ; хп > 0 .

(3.43)

Некоторые теоремы линейной алгебры сразу позволяют решить во­ прос о существовании допустимого решения ОЗЛП. Предположим, что из т уравнений системы (3.42) только г линейно независимы. Это значит, что каждое из остальных (т - г) уравнений можно получить, преобразуя рас­ смотренные г уравнений. Число г называется рангом системы уравнений.

Например, в системе

2 х ]+ З х 2 + * 3 = 1 ;

< х }- 2 х 2 - 4 * 3 = 2 ;

4 * j - х2 - 7*3 = 5

только первые два уравнения линейно независимы. Третье уравнение можно получить, если сложить первое уравнение с удвоенным вторым. Поэтому ранг данной системы уравнений равен 2.

Величина г не превышает числа переменных п: г < п . Случай, когда г = п мало интересен. При этом система уравнений (3.42) имеет единствен­ ное решение {**, xj,..., **}. Если все х* при этом неотрицательны (г = 1,

2, ..., и), то найденное решение системы уравнение (3.42) является одно­ временно допустимым и оптимальным решением ОЗЛП (3.41) - (3.43). Ес­ ли же хотя бы одна из величин х* отрицательна, то ОЗЛП не имеет допус­

тимых решений. Поэтому в дальнейшем будем считать, что г строго мень­ ше п. Тогда, если система (3.42) совместна, то она имеет бесчисленное множество решений. При этом (гг - г) переменным можно придавать про­ извольные значения. Эти переменные называются свободными. Остальные г переменных выражаются через свободные. Они называются базисными. Для существования допустимых решений ОЗЛП необходимо, чтобы среди множества решений системы (3.42) нашлись неотрицательные, т. е. удов­ летворяющие ограничениям (3.43).

Пример. Дана ОЗЛП

х\+ 4*2

min;

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]