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

Елтаренко Исследование операцыи 2007

.pdf
Скачиваний:
95
Добавлен:
16.08.2013
Размер:
2.72 Mб
Скачать

F

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