книги из ГПНТБ / Математическое программирование и производственные задачи
..pdfС другой стороны, необоснованное сокращение запасов может привести к срыву планов производства, так как неред ко предприятие-поставщик по различным причинам наруша ет заранее установленные сроки поставки материалов или же своевременно поставленные материалы поступают к потреби телю с недопустимой задержкой из-за неправильной организа ции работ в транспортных организациях. Хотя за задержку поставки соответствующая организация (предприятие-постав щик или транспортная организация) вынуждена платить штраф, однако эта сумма только частично компенсирует убы ток предприятия-потребителя. Кроме того, в результате ука занных нарушений наносится огромный ущерб всему народ ному хозяйству. Из изложенного видно, что существует опре деленный уровень запасов, при котором убытки предприя тия минимальны. Определением наивыгоднейшего уровня за пасов занимается теория управления запасами, являющаяся одной из наиболее молодых и одновременно бурно развиваю щихся дисциплин прикладной математики. Вопросам теории управления запасами посвящено большое количество работ отечественных и зарубежных авторов *.
Во многих работах задачи оптимального управления за пасами сводятся к задачам динамического программирова ния. Однако часто можно пользоваться более простыми ли нейными моделями. В настоящем параграфе разработана ли нейная экономико-математическая модель оптимального уп равления запасами и дан алгоритм решения поставленной за дачи.
Поскольку при планировании производства заранее уста навливаются нормы расходов средств производства, то при постановке математической задачи все нормативы предпола гаются известными. В задаче оптимального управления запа сами в качестве критерия принимается суммарный объем за пасов материалов (в денежном выражении), который рас сматривается не в динамике, а как постоянная величина в те чение определенного промежутка времени. Хотя такое пред
* Достаточно обширный перечень работ по задачам управления за пасами приведен в библиографии книги [49].
100
положение не полностью соответствует реальной действитель ности, однако результаты решения задачи пригодны для боль шинства предприятий крупносерийного и массового производ ства.
Для составления линейной модели |
оптимального управ |
||||||||
ления запасами |
введем следующие обозначения: |
|
|||||||
a u(i— 1, 2, . . |
т\ у'= 1, 2, |
. . |
л) —норма |
расхода мате |
|||||
риала у-го |
вида |
на изготовление |
продукции |
/-го вида; |
|||||
Ai ( i= l, |
2, . . |
т.) |
—объем |
выпускаемой |
продукции /-го |
||||
вида |
за год; |
|
|
|
|
|
|
|
|
A\{i= 1, 2, . . |
т\ |
/= 1, 2, |
. . |
L) —объем |
выпускаемой |
||||
продукции /-го вида за /-ый промежуток времени (нап |
|||||||||
ример, квартал или месяц) года; |
|
|
|
|
|||||
У/ (/= 1, 2, . . |
л )—допустимый минимальный объем запа |
||||||||
са у'-го вида материала (страховой запас); |
|
||||||||
Уу (у=1, 2, . . |
л)—допустимый |
максимальный |
объем за |
||||||
паса у-го вида материала; |
|
|
|
|
|
||||
ci (j—1, 2, . . |
л )—стоимость единицы материала у-го вида; |
||||||||
^у(/==1» 2, . . |
л; /=1, 2, . . .,£ )— лимитированный объем |
||||||||
материала у-го вида за /-ый промежуток времени; |
|||||||||
У;(У— 1, 2, . . |
т; 1=1, 2, . |
. |
Z.)—объем запаса материа |
||||||
ла у-го вида на складе |
в конце /-го |
промежутка вре |
|||||||
мени; |
|
|
|
|
|
|
|
|
|
■*{(/=1 , 2, . . |
т\ |
I —1 , 2, |
. . |
Z.) —объем |
выпускаемой |
||||
продукции /-го вида за /-ый промежуток времени; |
|||||||||
2, . . |
л; /= 1 , 2 , . - . , / - ) — объем материала у-го |
вида, ожидаемый от поставщика в течение /-го проме
жутка |
времени; |
|
|
|
Модель оптимального управления запасами включает сле |
||||
дующие требования: |
|
|
|
|
обеспечение выполнения |
плана |
выпуска |
продукции: |
|
v |
x ^ A t , /= 1, 2, |
. ; ., |
т- |
(1.1) |
I- 1
ограничение на расход материалов:
2 а ц х \ ^ Ь 1~ 1+ у ‘- 1, j = 1, 2 , |
п\ 1=\, 2, |
|
( 1.2) |
|
i-i |
|
|
|
|
сохранение |
допустимых |
пределов запасов |
материалов: |
|
У/ < У;'«=У; |
. У'=1, 2, . |
га; /=1, 2, . . |
Z.; |
(1.3) |
неотрицательность объема выпускаемой продукции*: |
||||
jc/:>0, /=1, 2, . . |
/га; /=1, 2, . . |
L. |
(1.4) |
|
Величина |
в ограничениях (1.2) выражает |
оста |
ток материала /-го вида на складе к началу года. Значения этих величин определяются во время составления сводных заявок.
Как указывалось, в качестве критерия принимается сум марный объем запасов материалов (в денежном выражении) :
/ i = 2 |
2 |
о уг |
(1.5) |
t~i |
j= I |
|
Таким образом, ставится следующая задача.
Задача 1.1. Минимизировать линейную функцию (1.5) при
условиях (1.1) — (1-4) |
или (1.Г), |
(1.2), |
(1.3). |
|
||
Задача |
1.1 |
является статической |
моделью оптимального |
|||
управления |
запасами |
(текущие |
изменения |
объемов запасов |
||
в модели не фигурируют). Если |
величина |
ЬУ детерминиро |
||||
вана для |
всех |
/ и I, |
то задача |
1.1 является обычной зада |
чей линейного программирования с двусторонними ограниче ниями на переменные. Однако на практике часто поставка
требуемых материалов задерживается по |
разным причинам |
|||||
(недовыполнение планов |
поставщиков, |
нарушение |
сроков |
|||
транспортировки и т. п.). |
Следовательно, |
Ь\ является не де |
||||
терминированной, а |
случайно)! |
величиной, |
причем |
|
||
* Более естественно |
вместо |
(1.1) |
и (1 4 ) ввести |
ограничения |
|
|
ф М , г, 1 |
= 1, 2 ........... ш; |
/ 1, 2 .................L , |
(1.1') |
т. е. требовать выполнение плана не в течение года, а по каждому про межутку времени. Однако такие ограничения гораздо увеличивали бы размеры задачи.
102
0K b j^ b j, j = |
1, 2, . . |
Л ; Z = l, 2, . . |
(1.6) |
Таким образом, |
задача |
1.1 становится задачей |
стохасти |
ческого линейного |
программирования. |
|
Как указывалось, в зависимости от преследуемых целей в стохастическом программировании существуют различные по
становки задач. |
В данном параграфе рассматривается |
сле |
|
дующая |
постановка. |
|
|
Если |
закон |
распределения величины b\ U— 1,2, . . |
., л; |
1 = 1, 2, . |
. .. L) известен, то, заменив ее математическим ожи |
||
данием |
|
решение стохастической задачи можно свести |
к решению детерминированной задачи. Но закон распределе ния случайных величин не всегда бывает известен. Поэтому для выявления закона распределения и определения матема тического ожидания заранее производится обработка стати стических данных.
|
Задача 1.2. Найти такие значения |
величин х\ и yj, ко |
||||||||
торые минимизируют функцию |
|
|
|
|||||||
|
|
|
Л = л * |
\ i-i |
2 Су У} ) |
(1 - 7 ) |
||||
|
|
|
|
|
|
j - 1 |
/ |
|
||
при |
условиях |
|
|
|
|
|
|
|
|
|
Р ф |
а и х ‘- |
у'--1 <. Ь '-у > -р ‘, |
/= 1, 2, |
., л; 1=1, 2, |
. . .,L , |
|||||
|
|
|
|
|
|
|
|
|
|
( 1.8) |
где р] (у= 1 , |
2, |
. . ., |
л; |
/=1, 2, . . ., |
L)—заданные вероят |
|||||
ности; |
|
|
|
|
|
|
|
|
|
|
|
х\>А\, |
i = 1, |
2, |
. . |
., |
т- |
1 = 1, |
2............ L-, |
(1.9) |
|
y/sSy'sSy., |
, |
7= 1, |
2, |
. . |
., |
л; |
/= 1, |
2, . . ., L. |
(1.10) |
В работе [47] для решения задачи (1.7) — (1.10) предла гается алгоритм, основанный на обработке статистических данных относительно случайных величин Ь}1 . Но для приме нения указанного алгоритма необходимо иметь достаточно большой объем статистического материала. С другой стороны,
103
при большой номенклатуре выпускаемой продукции число не известных в задаче (а также число ограничений) будет вели ко, что приводит к определенным затруднениям при реше нии задачи. Поэтому возникает необходимость в упрощении структуры задачи. Если принять
х\= А \ , 1= 1 , 2, . . ., т- / = } , 2, . . L, |
(1.11) |
то сразу можно избавиться от этих переменных. Тогда вместо задачи 1.1 получим следующую задачу.
Задача 1.3. Найти величины yl, минимизирующие |
ли |
|||||||
нейную функцию (1.5), при условиях (1.3) и |
|
|
|
|||||
у^1 > й ‘- Ь ‘-\ |
у = 1, 2, . |
. ., |
п- |
/= 1, 2 ........... L, |
|
(1.12) |
||
где |
|
|
|
|
|
|
|
|
гп |
|
. . |
., |
п- |
/=1, 2, . . ., |
L. |
|
|
е1‘ = ^ а и А\, /= 1, 2, |
|
|
||||||
<-1 |
|
|
|
|
|
|
|
|
Из постановки видно, что при фиксированных |
Ь1 |
задача |
||||||
имеет простую структуру. Нетрудно заметить, что для |
разре |
|||||||
шимости задачи |
необходимо |
выполнение условий |
|
|
|
|||
у/, /=1, 2............я; |
/=1, 2, . . ., L. |
|
(1.13) |
|||||
Эти условия являются и достаточными. Действительно, |
||||||||
при выполнении |
условий |
(1.13) |
можно непосредственно по |
|||||
строить решение задачи |
1.3: |
|
|
|
|
|
|
|
у}~_1 = |
max [dlj —blj - 1-, |
уу]. |
|
(1.14) |
||||
Ясно, что стохастический вариант этой задачи тоже будет |
||||||||
иметь простую структуру по сравнению с задачей |
1.2. |
|
|
|||||
Из (1.14) видно, |
что при £ j_1= 6 / = |
rfj имеем у;-_1= уу, |
т. е. |
|||||
если в течение (/—1 )-го промежутка |
времени от поставщика |
получается точно тот объем материала, который употребляет ся в производстве за /-ый промежуток, то достаточно в нача ле (/—1)-го промежутка времени иметь минимальный объем запаса. Однако в реальных условиях редко встречается такая ситуация. Поэтому величины ЬУ рассматриваются как слу
104
чайные, вследствие чего случайными будут и величины
Если закон распределения случайных величин Tjj изве стен, то стохастический вариант задачи 1.3 сводится к детер минированному. Но закон распределения случайной величины не всегда известен. В таких случаях целесообразно вместо стохастического варианта задачи 1.3 рассматривать параме трическую задачу, в качестве параметров принимая случай ные величины vjj. Нетрудно показать, что минимальное зна чение целевой функции (1.5) является неубывающей вогнутой функцией от параметров Учитывая вогнутость и кусочно линейную структуру этой функции, заключаем, что ее градиент принимает возможное наибольшее значение на подобласти,
близкой к вершине ^l.= d l. (у' = 1, 2, . . ., tv, /= 1, 2, . . ., L)
параллелепипеда 0^Сг/.£^й‘. Следовательно, на указанной подобласти функция имеет наибольшее приращение, т. е. вблизи указанной вершины целевая функция задачи 1.3 бо лее чувствительна относительно исходных данных. Из оп ределения vjj видно, что минимальное значение целевой функции (1.5) является невозрастающей функцией от Ь}1
(этот факт чувствуется и из экономической сущности задачи
1.3).
Таким образом, наилучшей подобластью (в смысле ми нимального значения целевой функции) является параллеле пипед 0<v;j<yy. Исходя из этого, отрезок [0, _уу ] можно принимать в качестве доверительного интервала случайной ве личины у*. Сужение доверительного интервала производится
на основе обработки статистических данных. |
Можно убедить |
||||
ся, что в качестве |
минимального объема запаса ( уу) нужно |
||||
брать |
верхнюю |
границу |
доверительного |
интервала ■/)' |
|
</= 1, |
2, . . ., L). |
Практически уу не может неограниченно |
|||
убывать. Через некоторое время |
(исходя из статистики реали |
||||
зации |
случайных |
величин |
f]j ) |
она станет неподвижной или |
с незначительной амплитудой будет колебаться вокруг данной стационарной точки, которую и можно принять в качестве уу.
Приведем еще одну постановку задачи управления запа сами. Предположим, что из-за нехватки единицы материала /-го вида в течение / -го промежутка времени предприятие не сет ущерб в сумме «у . Введем переменные
105
z J = * J - > + y r f j, |
|
7=1, |
2 |
|
/=1, |
2, . . |
L. |
|
||
При 2 j> 0 |
имеем, что zlj = |
y lj , |
а при |
zy< 0—потребнос |
||||||
ти производства в материале /-го |
вида не |
удовлетворяются |
||||||||
в объеме |zj|. |
|
|
|
|
|
|
|
|
|
|
Задача |
1.4. |
Минимизировать функцию |
|
|
||||||
|
|
|
|
Л = |
V v |
ау|2}1 |
|
|
(1.15) |
|
|
|
|
|
|
/-1 |
|
|
|
|
|
при условиях |
|
|
|
|
|
|
|
|
|
|
-1| J_ ?l- 1 |
|
|
|
|
|
|
|
|
|
|
I < |
+ |
W -1- |
-zr |
/— 1, 2............я; |
/=1, 2, |
. |
. ., Z.; |
|||
|
1 |
1 |
|
|
|
|
|
|
(1.16) |
|
|
|
|
|
|
|
|
|
|
|
|
- d ^ z ^ y j |
, 7= 1, |
2, . . ., я; /=1, |
2, . . |
L. |
(1.17) |
|||||
По существу |
задача |
(1.15)— (1.17) |
является |
нелиней |
||||||
ной (даже при фиксированных |
bJ ). Если величины |
Ь1. слу |
чайные, то она будет нелинейной стохастической задачей, ко торая более правильно отражает реальную действительность.
§ 2. ОПРЕДЕЛЕНИЕ ОПТИМАЛЬНЫХ СХЕМ ПЕРЕВОЗОК
Организационные трудности в строительстве, как и хозяй ственные неполадки, подчас возникают в результате непра вильного планирования материально-технического снабжения, в частности, планирования перевозок строительных материа лов.
Повышение качественного уровня планирования и управ ления экономикой является важным условием улучшения организации производства в строительстве.
Сложившиеся организационно-технические условия тре буют новых хозяйственных решений во всем комплексе вопро сов экономики строительства. Изо дня в день повышается актуальность систематического совершенствования организа ции производства на базе внедрения научной организации труда и управления в строительстве.
106
Старые формы управления и планирования строительствавытесняются, заменяются новыми. Для научно обоснованного планирования и управления строительством требуются мате матические методы экономического анализа строительства и его производственных подразделений (в строгой увязке между собой).
Современные экономико-математические методы позволя ют дать оптимальные решения по каждой конкретной органи зационной и планово-экономической задаче, встречающейся в строительном производстве.
Необходимо отметить, что транспортные расходы в се бестоимости промышленной продукции пока составляют зна чительную долю. В промышленности удельный вес транс портных затрат составляет, в среднем 13%, а в ее отдельных отраслях—еще выше. Так, например, в угольной промышлен ности удельный вес транспортных расходов составляет 27,8%, в промышленности строительных материалов—28,3%, а в целом в стране транспортные расходы составляют более 27 миллионов рублей в год. Отсюда и вытекает необходи мость научно обоснованного подхода к вопросам учета тран
спортных затрат, планирования |
и |
нормирования |
перевозок |
||||||
в материально-техническом снабжении. |
|
|
|
||||||
Постановка транспортной задачи. Классическая транс |
|||||||||
портная задача |
ставится следующим образом. |
Пусть имеются |
|||||||
т баз |
(карьеров) |
Л„, . . |
., |
А,„—с объемами строитель |
|||||
ного |
материала |
г-го |
вида а[, |
а 2,г |
. . |
а'т и п потребителей |
|||
(строительных |
объектов) Въ |
Я>, |
. . ., |
В п—с соответствую |
|||||
щими потоебностями: b[, Ьг2, . . |
., |
Ь'п. Стоимость |
перевозки |
||||||
единицы материала г -го вида из |
I-го пункта отправления |
||||||||
(базы) в у'-ый пункт назначения |
равна |
crt. и |
предполагается |
известной для всех комбинаций (г, у). Пусть х\}—количество материала r -го вида, перевозимого по маршруту (г, у). Экономический смысл задачи заключается в том, чтобы оп ределить такие величины x tjr (объемы перевозок) для всех маршрутов (г, у), при которых суммарная стоимость пере возок была бы минимальной.
Условия транспортной задачи для материала г -ro вида схематично можно представить в виде табл. 2.1.
107
|
|
|
|
Таблица 2.1 |
Пункты назначе |
|
|
|
Объемы мате |
ния (строительные |
|
|
• * - п |
|
объекты) |
1 2 |
• • • j |
риала r -го вида |
|
Пункты отправле |
|
|
|
на базе |
ния (базы) |
|
|
|
|
1 |
х[\ х 12 |
' • • 4 , • • • Х\п |
|
|
|
|
|||
2 |
х 2\ х 22 |
’ • • Xnj |
• • • х 2п |
а2 |
i
т
Объемы потреб ления материала г-го вида
хп хя • • • хЬ • * * Х1п |
|
а[ |
||
х т\ |
х т2 * * •xmj |
хтп |
|
ат |
ь' |
у2 . . • Ь) |
|
т |
п |
• • • К |
2 |
a't = 2 ь; |
Очевидно, что таких таблиц составляется |
1 штук ( г = |
|
= 1, 2, |
. . /). Элемент х\., стоящий на пересечении /-ой |
|
строки |
и у'-го столбца, показывает количество |
материала, |
перевозимого из /-го пункта отправления в у'-ый пункт пот ребления. Следовательно,
v x 'j= a rlt i = 1, 2, . . ., т, i~i
т. е. количество материала r -го вида, отправленного из /-ой
базы во все |
строительные |
объекты, должно быть равно |
запасу материала r -го вида на /-ой базе. |
||
С другой стороны, |
|
|
т |
х и = ьр У=1, |
2............я, |
2 |
т. е. количество строительного материала r -го вида, получен ного у'-ым строительным объектом из всех баз, должно быть равным его потребности.
108
Обычно предполагается, что сумма запасов материала г-го вида всех баз равна сумме потребностей всех строитель ных объектов в материале данного вида. Если это условие нарушается, то доба!Вляется либо фиктивный пункт потребле ния с объемом потребления
v а/ — v br= b r |
если |
|
либо фиктивная |
|
±ua i |
+ J uj—vn+1 |
|
|
|
<-i |
/-i |
т |
п |
т |
|
п |
|||
база с |
запасом У |
Ьг.—2 л \~аГт+\ |
( если 2 |
^-^>2 а \)- |
|
j - 1 |
г-i |
}-\ |
<- 1 |
Математическая формулировка транспортной задачи за ключается в следующем.
Задача 2.1. Требуется определить такие объемы пере
возок x ifr |
строительных материалов |
из баз на |
строительные |
объекты, |
которые минимизируют |
суммарные |
затраты пере |
возок |
|
|
|
|
2 |
|
|
|
<-1 /-1 сЬх и |
|
|
при условиях
п
V хг — :а1 * = i , 2, . . ., от;
т
II
У = 1 , 2, . . ., л ;
x tj>r |
0, *= 1 , |
2, . . ., т; j = 1, 2............п. |
Условия |
х'^-О |
означают, что объемы перевозок между |
базами и строительными объектами не могут быть отрица тельными, т. е. перевозки допустимы только в одном направ лении.
Существуют разные модификации постановки транспорт ной задачи, в которых, кроме вышеуказанных факторов, учи тываются такие, как пропускная способность дорог, увеличе ние (или уменьшение) потребности по отдельным видам строительных материалов в различных пунктах потребления и т. д.
109