Файл: Шепелев, И. Г. Математические методы планирования и управления в строительстве конспект лекций.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 29.10.2024
Просмотров: 52
Скачиваний: 0
В общем случае, если число поставок за месяц равно т , то средний запас равен сумме интегралов
|
Y = |
Мк- |
(Ук- 1 |
+ |
Мк_! — Хк t) dt—j— |
||
|
|
"k-i |
|
|
|
|
|
|
га—1 |
Ti+j |
|
|
hi- |
|
|
|
|
|
|
|
|||
|
H-----— |
j |
(yk-i H ~ |
M k+1 |
+ |
qi — |
^k t)dt-j-■ |
|
j=1 |
Ti+j+ l |
|
|
i= 1 |
|
|
|
|
Mk |
|
|
|
ra |
|
+ |
—-----— -------г |
i" |
(Ук-i |
Mk_i + |
Q — Akt)dt. (8.3.11) |
||
|
Mk - Ti + m - 1 |
. J |
|
|
|
|
|
|
Месячный объем поставок qk |
определяем как сумму поста |
|||||
вок за k-й месяц |
|
П1 |
|
|
|
|
|
|
|
|
|
|
|
(8.3.12) |
|
|
|
|
qk = 2 Q i > |
|
|
||
|
|
|
i=i |
|
|
|
|
|
Минимизация |
функции L, |
определяемой |
соотношением |
(8.3.1), проводится методом динамического программирования. Поскольку начальный уровень запаса Уо и конечный у12 == у?а
заданы, то минимизацию можно проводить как в «прямом» вре мени, так и в обратном. Мы рассмотрим алгоритм решения за дачи в прямом времени. Представим задачу минимизации в ви
де 1 2 -шагового процесса. |
|
1,..., 12), определяе |
Разобьем диапазон изменения qk, (k = |
||
мый соотношением (8.3.3) |
на ряд значений с некоторым шагом |
|
Aq. На первом этапе всех |
допустимых значений q находим Уг, |
|
г/2 и 1Х. |
|
|
Li (Yj) = min (у0; qt) |
(8.3.13) |
|
|
д |
|
q
л
Lj — дает минимальные издержки за первый месяц, q, соответ ствующее этому значению, условно называется оптимальным значением. На втором этапе для каждого у\ и всех допустимых
значений q2 определяем У2; г/2; к (У). Находим
L2 (Y2) = |
min [l2 (yu q2) + L, (Yt)]. |
(8.3.14) |
||
Для нахождения |
оптимума за год применяем |
принцип |
||
Р. Веллмана (см. гл. V), |
решив к — 1 |
рекуррентных соотноше |
||
ний. |
|
|
|
|
Lk (Yk) = |
min [/k (yk_i, qk) + |
Lk_! (Yk_i)]. |
(8.3.15) |
112
Здесь Lk(yic) — минимально возможные издержки за первые к месяцев при условии, что в конце k-го месяца уровень запасов
равен |
На последнем шаге мы найдем L* = L12 (У12) |
и соот |
||
ветствующие этому значению |
у*р q*2, а |
следовательно |
найдем |
|
все оптимальные значения . |
q* . .. q* 2 |
и у*... y:j2, а также |
||
—* |
—* |
|
|
|
У1 • • • У12* |
значения |
поставок строительных |
||
Полученные оптимальные |
||||
материалов могут быть использованы при планировании |
мате |
риально-технического снабжения строек, а оптимальные значе ния запасов Ук для планирования загрузки оборотных фондов.
§ 8.4. Оптимальное планирование поставок строительных материалов
в стохастической постановке
Выше был рассмотрен метод решения задач по управлению запасами в детерминированной постановке. На практике усло вия, соответствующие такой постановке, встречаются редко, ча ще всего при управлении запасами наблюдаются неопределен ности, при учете которых задача становится стохастической. Природа этих неопределенностей обычно вызывается:
а) ошибками при определении потребности в строительных материалах;
б) неопределенностью поставок строительных материалов; в) источником неопределенности является и целевая функ
ция, как правило, корреляционная, например, формула (8.2.16), отражающая только средние значения, связанных формулой ве личин и имеющая ошибку аппроксимации, а, следовательно, и возможность ошибки расчета по этой формуле. Для решения задач с учетом этих неопределенностей необходимо применять методы стохастического программирования (см. гл. V).
В качестве метода решения задачи используется метод ди намического программирования, причем при подсчете среднего
запаса и в соотношении (8.3.2) |
величины qk и Sk |
рассматрива |
||
ются как случайные. Тогда |
|
|
|
|
f (Ук) |
= f (У к -0 |
+ f (qk) — f(S k), |
(8.4.1) |
|
где f (ук) — случайная |
величина остатка за k-ый период; |
|||
f (qk) — случайная |
поставка за k-ый период; |
|
||
f (Sk) — случайный спрос за k-ый период; |
|
|||
f (г/k-i) — случайный остаток за |
предыдущий период. |
|||
Меняется метод определения |
среднего запаса |
У, который |
5 И. Шепелев |
113 |
становится функцией |
случайных переменных |
f(t/k-i) и f(qk) |
при случайном спросе f (Sk): |
|
|
y k = |
F[f(yk_1); f (qk); f(Sk)].- |
(8.4.2) |
Минимизируемая функция L за k периодов выражается:
L* = minL{y1[f(y0); f(qi)5 f(SO] ...
•••yk[f(yk-i); f(qk); f(Sk)]}. |
(8.4.3) |
Если учесть, что функция Li, будучи корреляционной, имеет в своем составе случайную компоненту е, то затраты зависят не только от уровня запасов, но и сами являются случайными, то есть
Тогда |
L = |
I (У; |
е); |
( |
(8.4.4) |
|
|
|
|
|
|
L* = min L {... Ук [f (yk-0 ; |
f (qk); f (Sk)] ... *}• |
(8.4.5.) |
|||
Выражение (8.4.5) является общим видом целевой функции |
|||||
календарного |
планирования поставок |
материалов при полной |
|||
стохастической |
постановке |
задачи. Конкретные методы реше |
|||
ния этой задачи зависят от |
вида случайных функций |
f (yk-i); |
f (qk); f (Sk) и e.
Оптимизация должна производиться с применением принци па Р. Веллмана и реализацией задачи на каждом шаге метода ми стохастического программирования.
В работах [12, 13] приведен анализ поставок строительным трестам Южного Урала цемента, арматурной стали, кирпича. Статистический анализ . фактических поставок строительных материалов показал, что поставки могут быть аппроксимирова ны факторными моделями, цепями Маркова и другими случай ными законами.
В § 8.5 приведены основные положения и примеры аппрок симации поставок факторными моделями и цепями Маркова.
§ 8.5. Планирование поставок с аппроксимацией цепями Маркова
Рассмотрим постановку и метод решения задачи оператив ного оптимального календарного планировния поставок строи тельных материалов, при этом будем считать, что спрос (по требность) в строительных материалах является детерминиро ванным и нестационарным (определяется при календарном пла нировании строительно-монтажных работ), поставки имеют ве роятностный характер, причем, в качестве модели процесса принята простая неоднородная цепь Маркова.
114
Для построения модели, как и в § 8.3, введем следующие обозначения:
qn — поставка материала в течение n-ой недели; Ул — остаток материала на конец n-ой недели.
Требуется определить значения поставок и остатков для каждой недели планового периода, удовлетворяющие ограниче ниям
|
|
УN = |
Уопт» |
|
|
(8.5.1) |
|
Sn ' |
|
|
N |
|
|
(8.5.2) |
|
Уп-1 |
Q n |
2 |
^ |
Уп—1 |
Уопт> |
||
|
|
|
1= п |
|
|
|
|
|
|
|
п |
|
|
|
(8.5.3) |
|
О ^ У п " ^ 2 |
^ "1” Уопт» |
|
||||
|
|
|
1=п-Н |
|
|
|
|
при достижении минимума функции |
|
|
|
||||
|
|
N |
(Чп> Уп)> |
|
(8.5.4) |
||
|
|
2 |
|
||||
причем, |
|
П = |
1 |
|
|
|
|
Уп = Уп-1 “f* Qn |
Sn, |
|
(8.5.5) |
||||
|
|
||||||
где Sn — потребность за n-ю неделю; |
на |
конец месяца, |
|||||
Уолт— оптимальное |
значение |
остатка |
|||||
определенное в оптимальном годовом плане; |
|||||||
N — число недель в плановом периоде; |
|
||||||
/п (qn, Уп) — функция потерь для n-ой недели. |
|
|
|||||
Ограничение |
(8.5.1) |
обуславливается необходимостью при |
вести в соответствие оптимальный оперативный и оптимальный годовой планы.
В качестве критерия оптимальности принимается минимум функции потерь (см. § 8 .2 ).
Оптимизация плана осуществляется методом динамического программирования. Вычислительный процесс строится от конеч ного состояния к исходному, где под состоянием понимается ве личина остатка на начало недели. Тогда общее рекуррентное соотношение можно записать в следующем виде
Ln (Уп-i) = min [In (qn; yn)-f Ln_i (У„)], |
(8.5.6) |
Л |
|
q |
|
где Ln(«/n_ i) — минимальные потери для n оставшихся |
недель |
планового периода при условии, что остаток материала на нача ло п-ой недели равен уп-и
В процессе решения задачи для вычисления значений функ
5* |
115 |
ции потерь требуется определить величину среднего недельного
запаса, который |
является аргументом целевой функции. |
||
Из |
графика |
изменения |
уровня запаса в течение недели |
(рис. |
2 2 ) видно, |
что величина |
среднего запаса численно равня |
ется площади заштрихованной фигуры, деленной на число дней
Рис. 22. Схема формирования среднего запаса (стохастический подход)
в неделе. Тогда средний запас в n-ю неделю может быть опре делен по формуле:
Уп |
|
Sn (2) |
|
— -4-Чп f |
|
|
|
(8.5.7) |
( 1 ) ( 2 ) |
( 5 ) |
' |
qn; qn ■• • qn |
— поставки в течение первого, второго, ..., пя |
того дня недели;
m — число дней в неделе. .
Для определения объемов поставок используется простая не
однородная |
цепь Маркова, переходные вероятности которой, |
в |
||||
|
|
отличие от однородной цепи [13], меня |
||||
Т а б л и ц а 18 |
ются с течением времени. Под состоянием |
|||||
Состояние |
%выполнения |
системы |
здесь понимается процент вы |
|||
полнения |
недельного заказа на какой- |
|||||
системы |
недельного |
|||||
|
заказа |
либо день недели (табл. |
18). |
|
||
1 |
0 |
Математическая модель процесса пред |
||||
ставляет собой матрицы начальных и пе |
||||||
2 |
20 |
реходных |
вероятностей, |
полученные |
в |
|
3 |
40 |
результате обработки статистических дан |
||||
4 |
60 |
|||||
5 |
80 |
ных по фактическим поставкам материа |
||||
6 |
100 |
лов за два года. |
|
|
116