Елтаренко Исследование операцыи 2007
.pdfF |
F1 |
F2
у* |
r |
r' |
у |
Рис. 3.9. Целевая функция при y* |
r |
r |
|
Оптимальный размер заказа равен |
yопт |
r , а цена продукции рав- |
на q2.
Случай 3. Пороговый размер заказа – r r
Целевая функция в этом случае имеет вид (рис.3.10):
F |
F1 |
|
F2 |
у* |
r' r у |
Рис. 3.10. Целевая функция при r |
r |
Оптимальный размер заказа равен yопт |
r , а цена равна q2. |
140 |
|
Получим уравнение для определения величины r.
|
F ( y) |
|
T |
q |
|
K |
|
|
|
|
hy |
. |
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
1 |
|
|
|
|
1 |
|
|
y |
2 |
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Учитывая, что y* |
|
2K |
, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
h |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
F ( y* ) T q |
K h h 2K |
, |
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
1 |
|
|
|
1 |
|
2K |
2 |
|
|
h |
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
F (y*) T q |
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
2K h . |
|
|||||||||||||||||
|
1 |
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Условие для определения r : F ( y* ) |
F (r ) , |
|
||||||||||||||||||||
|
|
|
|
|
|
1 |
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
K |
hr . |
|
||||||||||||||
q |
|
|
2K h |
q |
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
1 |
|
|
|
|
|
|
2 |
|
|
|
|
r |
2 |
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
После упрощения получим уравнение для определения r:
(r )2 |
2r |
|
|
(q q ) |
|
|
2K |
|
|
2K |
0 . |
||||
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
h |
|
1 |
2 |
|
|
h |
|
|
h |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
||||||||||
С учетом y* |
2K |
|
|
уравнение перепишется в виде: |
|||||||||||
h |
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
(r )2 2r |
|
|
(q q ) y* |
|
( y* )2 0. |
||||||||||
|
|
|
|||||||||||||
|
|
|
|
|
|
h |
1 |
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Многопродуктовая статическая модель с ограничением на емкость склада
Число видов продукции – n, индексы продукции – i=1, 2, .., n. Цены на все виды продукции неизменны. Для каждого вида продукции заданы:
i – спрос в единицу времени,
Ki – затраты на размещение заказа,
hi – затраты на хранение 1 ед. продукции в единицу времени.
141
Кроме этого заданы объемные характеристики одной единицы продукции каждого вида – ai (i = 1, ... ,n), а также общая вместимость склада A (объем в м3 или площадь в м2).
Требуется определить оптимальные размеры заказов по каждому виду продукции yi* .
Целевая функция имеет вид:
n |
|
Ki i |
|
hi yi |
F |
|
|
||
|
yi |
2 |
||
i 1 |
|
|||
при ограничении |
|
|
|
|
n |
|
|
|
|
|
ai yi |
A . |
||
i 1 |
|
|
|
|
Для минимизации целевой функции воспользуемся методом Лагранжа. Составим функцию Лагранжа:
n |
Ki i |
|
hi yi |
n |
|
L( , y1...yn ) |
|
ai yi A , |
|||
yi |
2 |
||||
i 1 |
i 1 |
где – множитель Лагранжа (меньше нуля).
Для определения yi* составим систему уравнений:
dL |
|
Ki i |
hi |
ai 0 |
|
dy |
|
y2 |
2 |
|
|
i |
|
i |
|
(i 1,...,n). |
(3.2) |
dL |
n |
|
|||
|
|
|
|||
|
ai yi |
A |
0 |
|
|
d |
|
|
|||
i 1 |
|
|
|
Из (3.2) получим
y* |
|
2Ki i |
|
, |
0 . |
(3.3) |
|
||||||
i |
hi 2 ai |
|
|
|||
|
|
|
|
При = 0 получаем задачу без ограничения:
y*i |
|
2Ki i |
|
. |
|
||||
|
|
hi |
Если при этом выполняется ограничение
142
n
a y* A, |
(3.4) |
i i
i 1
то задача решена.
В случае невыполнения ограничения (3.4) систему уравнений решаем численными методами.
Для этого, на каждом шаге t уменьшая , рассчитываем по формуле (3.3) yi*(t) и проверяем ограничение. Итерационную процедуру
заканчиваем, когда будет выполнено ограничение (см. таблицу ниже).
|
|
|
|
|
n |
|
y1 |
y2 |
… |
yn |
ai yi A |
|
|
|
|
|
i 1 |
0 |
y*(1) |
y*(1) |
… |
y*(1) |
0 |
|
1 |
2 |
|
n |
|
-0.05 |
y*(2) |
y*(2) |
… |
y*(2) |
0 |
|
1 |
2 |
|
n |
|
Детерминированная динамическая однопродуктовая модель управления запасами
В данной модели спрос и другие стоимостные характеристики изменяются во времени. Разобьем весь период управления T на m интервалов (месяц, неделя), j=1, 2,…,m. Для каждого интервала должны быть заданы (или спрогнозированы):
j – спрос на продукцию в интервал j времени, Kj – затраты на размещение заказа,
hi – затраты на хранение одной единицы продукции в единицу времени,
qi – цена продукции в интервал j времени.
Кроме этого, необходимо задать уровень запасов на начало планового срока – P0 (j=1).
В модели не допускается дефицита, а время выполнения заказа (L) считается неизменным.
143
Ставится задача определения оптимального размера заказа в каждый из интервалов – y*j (j=1, 2,..., m), причем в некоторые интерва-
лы может быть |
y* |
0. |
|
j |
|
Целевая функция на отдельном интервале времени представляет собой затраты Fj ( j , Kj ,hj ,qj , Pj , yj ) , в зависимости от количества
запасов на начало периода Pj решается вопрос о размере заказа yj. Затраты в этот интервал времени определяются по формуле:
F( yj ) qj j |
K j j |
|
hj yj |
(см. однопродуктовую статическую |
|
yj |
2 |
||||
|
|
модель).
Интегральная целевая функция имеет вид:
|
m |
F |
Fj ( j , K j ,hj ,qj , Pj , yj ) , |
|
j 1 |
где Fj – затраты на функционирование системы управления запаса-
ми в интервал j.
Поставленная задача является сложной комбинаторной, но так как интегральная целевая функция представляет собой сумму, то для ее минимизации можно использовать динамическое программирование
[2].
3.3. Вероятностные модели управления запасами
В моделях этого класса спрос ( |
) – случайная величина с плот- |
ностью функции распределения ( |
) и математическим ожиданием |
M( ). Далее рассматриваются однопродуктовые вероятностные модели.
Упрощенная вероятностная модель
Будем считать, что время выполнения заказа (L) – постоянная величина. Обозначим через z суммарный спрос за время L, тогда f(z)
144
– плотность функции распределения вероятностей суммарного спроса, а F(z) – функция распределения вероятностей суммарного спроса.
На рис. 3.11 представлено изменение уровня запасов в вероятностных системах, РЗ – уровень резервного запаса.
P
Точки заказа
РЗ+M(z) |
|
z) |
y |
|
y |
РЗ |
|
|
t |
L |
L |
Рис. 3.11. График уровня запасов в вероятностных системах
Требуется определить РЗ при условии, что вероятность дефицита не больше, чем , т.е. P(z РЗ M(z)) .
Используя заданную функцию распределения вероятностей, определяем квантиль z , соответствующий вероятности 1- .
Искомый уровень запасов определяется в соответствии со следующим выражением: РЗ z M(z) .
Пример.
Дано: L = 4 недели – время выполнения заказа, плотность
функции распределения спроса f( |
)=N(100,10) – нормальное рас- |
пределение с M( )=100 в неделю и |
( )=10. Требуется определить |
уровень резервного запаса, при условии, что вероятность возможного дефицита не более 0,2.
Решение.
145
Математическое ожидание суммарного спроса за время выполне-
ния заказа: M(z) |
L M( ) |
400; дисперсия – |
||
D(z)=L 2( )=400, |
|
|
|
|
z |
400 |
20; f(z) – нормальное распределе- |
||
|
|
|
|
|
ние N(400, 20). |
|
|
|
|
Находим квантиль z , соответствующий вероятности 1- =0,8. По таблице нормального распределения F(Х )=0,8, Х =0,9.
Так как Хсоответствует нормированной случайной величине, то
|
z M (z) |
X , откуда |
z X |
|
M(z) , а уровень резервного |
|
|
z |
|||
|
|
||||
|
|
|
|
|
|
|
z |
|
|
|
|
запаса будет равен: |
|
|
|
||
|
|
РЗ = Х |
z = 0.9 20 = 18. |
Чтобы определить y* (уровень оптимального заказа), воспользуемся формулой Уилсона (3.1), подставив в нее M( ) (в этом заключа-
ется упрощенность модели):
y* 2KM( ) .
h
Отметим, что если L – случайная величина, и задана (L) – плотность функции распределения времени выполнения заказа, и, кроме того, если задана g(z/L) (условная функция плотности распределения спроса за время L), то безусловная функция плотности распределения f(z) определяется в соответствии с выражением:
f (z) g(z / L) (L)dL.
0
Имея f(z), можно определить необходимый резервный уровень запаса.
Вероятностная модель с контролем уровня запаса
В данной модели точка заказа определяется оптимальным уровнем резервного запаса R, который надо определить. Кроме этого, необходимо определить и оптимальный размер заказа y*.
Допущения, принятые в модели:
146
а) L – случайная величина;
б) плотность функции распределения вероятностей спроса f( ) за время выполнения заказа L не зависит от момента заказа;
в) неудовлетворенный спрос накапливается; г) в любой момент времени имеется не более одного невыполнен-
ного заказа.
В качестве исходных данных дано:
(L) |
f (z) g(z / L) (L)dL |
– функция плотности рас- |
|
g(z / L) |
|||
0 |
|
||
|
|
пределения спроса за время выполнения заказа; D – ожидаемый годовой спрос;
h – затраты на хранение в единицу времени на одну единицу продукции;
K – затраты на размещение заказа;
P – потери от дефицита одной единицы продукции.
Составим целевую функцию системы управления запасами, Dy –
количество заказов, которые будут сделаны в течение года. Целевая функция для оптимизации включает три составляющие:
|
|
|
|
|
|
F Cразм Cхран Сдефиц. |
|
|||||
Получим выражения для каждой из составляющих. |
||||||||||||
C |
D K – годовые затраты на размещение заказов; |
|||||||||||
разм |
y |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Cхран |
h H – затраты на хранение, где H – средний уровень за- |
|||||||||||
пасов. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
Будем определять H как среднее значение уровня запасов между |
||||||||||||
моментами получения заказа: |
|
|
|
|
|
|||||||
|
|
|
|
M(R z) y M(R z) |
|
y |
|
|||||
|
H |
M(R z) ; |
||||||||||
|
|
|
|
|
|
|
|
|
|
|||
|
2 |
|
2 |
|
||||||||
|
|
|
|
|
|
147
|
|
|
M(R z) |
(R z) f (z)dz R M(z) . |
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
Таким образом, получим: C |
h( |
y |
|
R |
M(z)) . |
|
||||||||||||||
2 |
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
хран |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
D |
|
|
|
|
|
|
|
|
|
|
|
||||||||
С |
dP – потери от возможного дефицита, где d |
– средняя |
||||||||||||||||||
дефиц |
y |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
величина дефицита; P – потери от дефицита. |
|
|
|
|
|
|
||||||||||||||
Дефицит (d) – случайная величина: |
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
d |
0 |
R |
|
z |
R |
, |
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
z |
|
z |
R |
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
d |
|
df (z)dz |
|
|
(z |
R) f (z)dz . |
(3.5) |
||||||||||
|
|
|
|
|
|
0 |
|
|
|
R |
|
|
|
|
|
|
|
|
|
|
В результате получаем целевую функцию: |
|
|
|
|
|
|
||||||||||||||
|
|
|
|
D |
|
y |
|
|
D |
|
|
|
|
|||||||
|
|
F( y, R) |
|
K h( |
|
R M(z)) |
Pd . |
|
||||||||||||
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
y |
|
2 |
|
|
|
|
|
y |
|
Чтобы найти оптимальные значения y и R, возьмем частные про-
изводные. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
F |
DK |
h |
D |
|
|
|
|
|
|
||||||
|
|
|
|
Pd |
0 ; |
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
y |
|
y2 |
2 |
y2 |
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
F |
|
h |
DP dd |
0 . |
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
R |
y |
dR |
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
||||||||||||
Остановимся подробнее на определении |
d(d ) . |
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dR |
|
|||
Заметим, что: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
d |
|
g(x)dx |
|
d |
|
G( ) G(R) |
g(R) , |
|||||||||||
|
dR R |
|
dR |
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
где g(x) – произвольная функция. Тогда, с учетом (3.5),
148
|
|
|
|
|
|
|
|
|
|
|
dd d |
d |
|
|
|
||||||
|
|
|
R zf (z)dz |
|
|
RR f (z)dz |
Rf (R) |
R f (z)dz Rf (R) |
||
dR |
dR |
dR |
f (z)dz .
R
Таким образом, для определения оптимальных y и R необходимо решить следующую систему уравнений:
|
|
|
|
DK |
|
h |
D |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
Pd |
0, |
|
|
|
|
|
(3.6) |
||||||||||||||||
|
|
|
|
|
y2 |
2 |
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
y2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
h |
DP |
|
f (z)dz |
0. |
|
|
|
|
|
|
(3.7) |
||||||||||||
|
|
|
|
y |
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
R |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
Из (3.6) получим: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
y* |
|
2D(K |
Pd ) |
, |
|
|
|
|
|
(3.8) |
|||||||||||||
|
|
|
|
|
|
|
|
|
h |
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
а из (3.7): |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
f (z)dz |
|
|
hy* |
, |
|
|
|
|
|
(3.9) |
|||||||||||
|
|
|
|
|
R* |
|
|
DP |
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
R* ) f (z)dz |
|
|
|
|
|
|
|
|
|
R*hy* . |
|
|||||||||||
|
d |
(z |
|
|
|
|
|
zf (z) |
(3.10) |
||||||||||||||||||
|
|
|
R* |
|
|
|
|
|
|
|
|
|
|
|
R* |
|
|
|
DP |
|
|||||||
Найдем y* |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
из (3.8), приняв d |
0, т.е. отсутствует дефицит, или |
|||||||||||||||||||||||||
min |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
уровень запасов R |
: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
y* |
|
|
|
2DK |
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
, d |
0. |
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
min |
|
|
|
|
|
h |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
||||||||||||||||||||
При уровне запасов R |
0, в соответствии с (3.5), d |
M(z) . |
|||||||||||||||||||||||||
Из (3.9) получим: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
yˆ* |
|
2D(K PM(z)) |
. |
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
h |
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
149 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|