4439
.pdf11
Итак, процедура оптимизации затрат многошаговая с интервалами от t = 5 до t = 0. На каждом шаге оценивается каждое из возрастных состояний оборудования.
Формируем схему возрастных состояний и переходов из состояния в состояние на каждом к-ом году.
Рис. 1 Стрелки вверх показывают эксплуатацию оборудования, вниз и вправо - замену
старого и эксплуатацию нового оборудования.
Рассчитаем таблицу переходов (табл. 2.2.2) и нанесем полученные значения на схему (см. рис. 1).
Таблица 2.2.2
t |
0 |
1 |
2 |
3 |
4 |
5 |
(t) |
500 |
600 |
800 |
1 000 |
1 500 |
- |
|
|
|
|
|
|
|
g(t) |
- |
3 000 |
2 500 |
2 000 |
1 200 |
1 000 |
|
|
|
|
|
|
|
r(t) |
- |
1 500 |
2 000 |
2 500 |
3 300 |
- |
Выполним оптимизацию на каждом шаге. 5 шаг.
Каждому значению Ft,k 5 присваиваем значения ликвидной стоимости из таблицы
переходов. Значения записываем в кружках с отрицательным знаком.
Тогда F55 = -1 000, F45 = -1 200, F35 = -2 000, F25 = -2 500, F15 = -3 000.
4 шаг.
Запишем функциональное уравнение, используя (2.2.4), при к = 4:
(t) |
F |
|
, |
если |
Х к |
Х э |
Ft,4 min |
t,5 |
|
|
|
(2.2.4) |
|
|
|
|
|
Х к |
||
t r(t) |
F |
, |
если |
Х з . |
||
|
t,5 |
|
|
|
|
|
Подставим данные и получим следующие значения затрат:
F4,4 min |
1 |
500 |
( |
1 |
000) |
min |
500 |
300 ; |
4 |
3 |
300 |
( |
3 |
000 |
|
300 |
|
|
min |
1 |
000 ( |
1 |
2000) |
min |
|
200 |
500; |
|||
F3,4 |
2 |
500 |
( |
3 |
000) |
|
500 |
|||||
|
|
|
||||||||||
|
3 |
|
|
|
||||||||
F2,4 min |
800 |
( |
1 |
200) |
|
min |
1 |
200 |
1 200 ; |
|||
2 |
000 |
( |
3 |
000) |
1 |
000 |
||||||
2 |
|
|
|
|
|
12 |
|
|
|
|
|
F1,4 min |
600 ( |
2 |
500) |
min |
1 |
900 |
1 900. |
|
1 500 |
( |
3 000) |
1 |
500 |
||||
1 |
|
|
Нанесем на схему двойные стрелки, соответствующие минимальному планируемому значению затрат при каждом возрасте оборудования t.
3 шаг.
Запишем функциональное уравнение при к = 3:
(t) |
F |
, |
если |
Х к |
Х э , |
Ft,3 min |
t,4 |
|
|
|
(2.2.5) |
|
|
|
|
||
t r(t) |
Ft,4 , |
если |
Х к |
Х з . |
Подставим значения и получим следующие значения затрат:
F3,3 |
min |
1 |
000 |
300 |
900) |
min |
1 |
300 |
600; |
||
|
3 |
2 |
500 |
( |
1 |
|
600 |
|
|||
F2,3 |
min |
800 |
( 1 |
( |
500) |
min |
300 |
100; |
|||
|
2 |
|
2 |
000 |
1 900) |
|
|
100 |
|
||
F1,3 |
min |
600 |
( |
1 |
200) |
min |
|
600 |
600 . |
||
1 |
500 |
( |
1 |
900) |
|
400 |
|||||
|
1 |
|
|
|
Нанесем полученные результаты на схему и выделим двойными стрелками перспективное развитие на данном этапе.
2 шаг.
Запишем функциональное уравнение при к = 2:
(t) |
F |
, |
если |
Х к |
Х э |
Ft,2 min |
t,3 |
|
|
|
(2.2.6) |
|
|
|
|
||
t r(t) |
Ft,3 |
|
если |
Х к |
Х з . |
Подставим значения и получим следующие значения затрат:
F2,2 min |
800 |
600 |
600) |
min |
1 |
400 |
1 400 ; |
|||
2 |
|
2 |
000 |
( |
|
1 |
400 |
|
||
F1,2 |
min |
600 |
100 |
min |
700 |
700. |
||||
|
|
1 |
1 |
500 |
( 600) |
|
|
900 |
|
1 шаг.
Запишем функциональное уравнение при к = 1:
(t) |
F |
, |
если |
Х к |
Х э |
Ft,1 min |
t,2 |
|
|
|
(2.2.7) |
|
|
|
|
||
t r(t) |
Ft,2 , |
если |
Х к |
Х з . |
Подставим данные и получим следующие значения затрат:
|
F1,1 |
min |
600 1 |
400 |
min |
2 |
000 |
2 |
000. |
|
|
|
1 |
1 |
500 |
700 |
|
2 |
200 |
|
|
Общая сумма |
затрат, |
приходящаяся на эксплуатацию и замену устаревшего |
||||||||
оборудования за 5 |
лет составит: P0 |
(t 0) |
F1,1 |
4 |
000 |
500 |
2 000 6 500. |
Итак, предприятие рассматривает 2 варианта управления:
1.Приобрести новое оборудование и эксплуатировать его 3 года. Затем продать и приобрести новое, которое будет эксплуатироваться еще 2 года.
2.Приобрести новое оборудование и эксплуатировать его 2 года. Затем продать и купить новое и использовать его еще 3 года.
13
В результате эксплуатации оборудования по двум вариантам предприятие несет одинаковый объем затрат, равный 6 500 у.е.
3. Теория массового обслуживания
Теория массового обслуживания – область прикладной математики, занимающаяся анализом процессов в системах производства, обслуживания, в которых однородные события повторяются многократно. Такие системы носят название системы массового обслуживания (СМО), и в качестве примера могут выступать предприятия бытового обслуживания; системы приема, переработки и передачи информации и др.
Каждая СМО состоит из определенного числа обслуживающих единиц (приборов, устройств, пунктов), которые называются каналами обслуживания. По числу каналов СМО подразделяют на одноканальные и многоканальные. Заявки поступают в СМО обычно не регулярно, а случайно, образуя так называемый случайный поток заявок (требований).
Предметом теории массового обслуживания является установление зависимостей между характером потока заявок, числом каналов обслуживания, производительностью отдельного канала и эффективным обслуживанием с целью нахождения наилучших путей управления этими процессами.
СМО делят на два основных типа: СМО с отказами и СМО с ожиданием (очередью). В СМО с отказами заявка, поступившая в момент, когда все каналы заняты, получает отказ, покидает СМО и в дальнейшем процессе обслуживания не участвует. В СМО с ожиданием заявка, пришедшая в момент, когда все каналы заняты, не уходит, а становится в очередь на обслуживание.
СМО с ожиданием подразделяются на разные виды в зависимости от того, как организована очередь: с ограниченной или неограниченной длиной очереди, с ограниченным временем ожидания и др.
Важное значение имеет дисциплина обслуживания, т.е. порядок обслуживания поступивших заявок. Обслуживание заявки может быть организовано по принципу
«первая пришла – первая обслужена», «последняя пришла – первая обслужена» или обслуживание с приоритетом.
Последовательность однородных событий называется потоком событий. Чаще всего входящий поток требований описывается функцией распределение вероятностей Пуассона:
|
e * |
r |
(3.1) |
|
P(r) |
|
|
, |
|
r! |
|
|||
|
|
|
|
где r – число требований, поступивших в систему массового обслуживания за единицу времени;
P(r) – вероятность того, что в систему массового обслуживания за единицу времени поступило r требований;
– интенсивность потока требований. Рассмотрим некоторые модели задач ТМО.
3.1. Многоканальная СМО с отказами в обслуживании
На вход многоканальной СМО с отказами поступает пуассоновский поток заявок с интенсивностью . Интенсивность обслуживания заявок равна .
Примерами многоканальных систем массового обслуживания являются предприятия с несколькими телефонными каналами, автозаправочные станции и др.
14
Случайный процесс, протекающий в системе, описывается формулами Эрланга:
n |
k |
1 |
k |
|
P0 |
|
; Pk |
|
Р0 ; |
|
. |
(3.1.1) |
|
|
|
|||||
k 0 |
k! |
|
к! |
|
|
|
Вычислив все вероятности состояний n–канальной СМО с отказами p0, p1, p2,…, pn, можно найти характеристики системы обслуживания.
Вероятность отказа в обслуживании определяется вероятностью того, что поступившая заявка на обслуживание найдет все n каналов занятыми, система будет находиться в состоянии Sn:
|
n |
|
|
pотк pn p0 * |
|
; к = n. |
(3.1.2) |
|
|||
|
n! |
|
В системах с отказами события отказа и обслуживания составляют полную группу событий, поэтому
pотк pобс 1. |
(3.1.3) |
Относительная пропускная способность определяется как
Q pобс 1 ротк 1 рn . |
(3.1.4) |
Абсолютную пропускную способность можно определить по формуле
A |
* pобс . |
(3.1.5) |
Вероятность обслуживания (доля обслуженных заявок), определяет относительную пропускную способность СМО и может быть найдена по формуле
Q pобс |
nз |
. |
(3.1.6) |
|
Среднее число заявок находящихся под обслуживанием, или среднее число занятых обслуживанием каналов определяется по формуле
|
|
А |
. |
(3.1.7) |
|
nз |
* pобс |
||||
|
Коэффициент занятости каналов обслуживанием рассчитывается по формуле
К з |
nз |
|
|
|
* pобс . |
(3.1.8) |
n |
|
|
|
|||
|
|
|
n |
|
Вероятность занятости каналов обслуживанием определим как
pзан |
|
|
t з |
|
. |
(3.1.9) |
||
|
|
|
|
|
||||
t з |
tпр |
|||||||
|
|
|
|
Среднее время простоя каналов найдем по формуле
|
|
|
1 |
pзан |
|
|
||
tпр tобс * |
. |
(3.1.10) |
||||||
|
|
рзан |
||||||
|
|
|
|
|
|
|
Среднее время пребывания заявки в системе определим как
T |
nз |
|
. |
(3.1.11) |
|
|
|||
СМО |
|
|
|
|
Пример 1. Предприятие принимает заявки на выполнение работ по трем телефонным линиям. В среднем поступает 80 заявок в час. Среднее время переговоров составляет 3 мин. Дать характеристику СМО.
|
|
|
|
|
|
|
|
|
|
|
Решение. n = 3; |
80ед / ч; tобс 3 |
мин . |
|
|
||||||
|
|
|
|
|
|
|
|
|
3 |
|
Определим интенсивность загрузки: |
|
|
tобс 80 |
4. |
||||||
|
|
|
||||||||
|
|
60 |
||||||||
|
|
|
|
|
|
|
|
|
|
15
Находим вероятность того, что обслуживанием не занят ни один канал:
|
n |
k |
1 |
|
0 |
|
|
|
1 |
|
|
|
2 |
|
|
|
|
3 |
1 |
|
|
|
40 |
|
41 |
|
|
42 |
|
43 |
1 |
||||||||||
p0 |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
0,036. |
||||||||||||||||||||||||
k 0 |
к! |
0! |
|
|
|
1! |
|
|
|
|
2! |
|
|
|
|
3! |
|
|
0! |
|
1! |
|
2! |
3! |
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||
Найдем вероятность того, что занят один канал: |
p1 |
|
|
|
* p0 |
4 * 0,036 0,144; |
|||||||||||||||||||||||||||||||||||
1! |
|||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
- |
заняты два канала: |
|
p2 |
|
2 |
* p0 |
|
|
|
42 |
* 0,036 |
0,288; |
|
|
|
|
|
|
|||||||||||||||||||||||
|
2! |
|
|
|
2! |
|
|
|
|
|
|
||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
- |
заняты три канала: |
p3 |
|
3 |
* p0 |
43 |
|
* 0,036 |
0,38. |
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||
|
3! |
3! |
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
Относительная пропускная способность определяется как: |
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||
|
|
|
|
Q pобс |
1 |
|
|
ротк |
1 0,384 0,616. |
|
|
|
|
|
|||||||||||||||||||||||||||
Абсолютная пропускная способность: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
A |
|
pобс |
80 |
0,616 |
49,28 ед/ч . |
|
|
|
|
|
|||||||||||||||||||||||||||
Среднее число занятых обслуживанием каналов: |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
nз |
|
|
|
|
|
pобс |
4 |
0,616 |
2,464. |
|
|
|
|
|
|
|
|||||||||||||||||||
Коэффициент занятости каналов: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
К з |
|
|
|
nз |
|
|
|
2,464 |
0,82. |
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
n |
|
|
|
|
|
|
3 |
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Доля потерянных заявок составляет около 38%, а обслуженных – 62%. Абсолютная пропускная способность СМО – 49 заявок в час. Каждый из каналов занят обслуживанием всего на 82% рабочего времени.
3.2. Одноканальная СМО с ограниченной длиной очереди
В отличие от многоканальной СМО с отказами, число мест в очереди ограничено m. Следовательно, заявка, поступившая в тот момент, когда все места в очереди заняты, не принимается к обслуживанию, не встает в очередь и покидает систему.
Пример 2. На автомойку в среднем за час приезжают три автомобиля, если в очереди уже находятся два автомобиля, то вновь подъезжающие автомобили не желают терять время в ожидании обслуживания и покидают мойку, поскольку среднее время мойки одного автомобиля составляет 15 мин, а мест для мойки всего одно. Необходимо провести анализ работы системы обслуживания в течение рабочего дня (8 ч).
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Решение. |
3 авт / ч ; tобс |
|
15 мин ; m = 2; n = 1; tраб = 8 ч. |
||||||||||||||||||
Находим интенсивность обслуживания: |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
1 |
|
|
1 |
0,067 |
|
авт / мин |
4 авт / ч . |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
15 |
||||||||||||||
|
|
tобс |
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
Определяем интенсивность нагрузки: |
|
|
3 |
0,75. |
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
||||||||||||||
|
|
4 |
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
Предельные вероятности состояний равны: |
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
р0 р1 р2 рn m |
pз |
|
|
|
1 |
|
|
|
1 |
|
0,25. |
|||||||||
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
m 2 |
2 |
|
2 |
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
Вероятность отказа в обслуживании: pотк |
рз |
0,25. |
|
|
|
||||||||||||||||
Относительная пропускная способность: Q |
1 |
|
|
p з |
1 |
|
0,25 0,75. |
||||||||||||||
Абсолютная пропускная способность: А |
|
|
* Q |
3* 0,75 |
|
2,25авт / ч. |
|
|
|
16 |
|
|
|
|
|
|
|||
Средняя длина очереди: |
Lоч |
m * (m |
1) |
2 * 3 |
0,75авт. |
|||||||
|
|
|
|
|
|
|
||||||
2 * (m |
2) |
2 * 4 |
||||||||||
|
|
|
|
|
||||||||
Среднее время простоя в очереди: Tоч |
|
|
Lоч |
0,75 |
0,25ч. 15 мин . |
|||||||
|
|
|
|
|
|
3 |
|
|||||
|
|
|
|
|
|
|
|
|
|
|
Среднее число заявок, находящихся на обслуживании:
Lобс |
m |
1 |
3 |
0,75авт. |
||
|
|
|
|
|||
m |
2 |
4 |
||||
|
|
|||||
Среднее число заявок в системе: Lсмо |
|
Lобс Lоч 0,75 0,75 1,5 авт. |
Среднее время пребывания автомобиля в системе обслуживания мойки:
Тсмо |
Lсмо |
1,5 |
0,5ч 30 мин . |
|
|
|
3 |
||
|
|
|
Таким образом, средняя доля потерянных заявок составляет 25% из числа
поступивших, что при 8-часовом рабочем дне составит: 3 |
8 |
0,25 6 автомобилей. |
3.3. Многоканальная с ограниченной длиной очереди |
|
|
На вход СМО поступает пуассоновский поток заявок |
с интенсивностью , а |
|
интенсивность обслуживания каждого канала составляет |
. Максимально возможное |
число мест в очереди ограничено величиной m.
Вероятность отказа в обслуживании наступает, когда все n каналов и все m мест в очереди заняты:
|
|
n |
m |
|
pотк |
рn m |
|
|
* p0 |
nm |
|
|||
|
|
* n! |
Среднее число заявок, находящихся в очереди, будет равно:
|
|
|
|
|
|
m |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
n 1 |
1 |
|
* m 1 |
m * |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
n |
|
|
n |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
Lоч |
|
|
* |
|
|
|
|
* p0 , если |
|
|
|
1, |
иначе |
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
n * n! |
|
2 |
|
|
|
|
n |
||||||||||||||||||
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
n |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
n 1 |
m |
(m |
1) |
|
|
|
|
|
|
||||||
|
|
|
|
|
|
Lоч |
|
|
|
|
|
p0 . |
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
n * n! |
|
|
|
2 |
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
Среднее время ожидания в очереди: |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
T |
|
|
|
Lоч . |
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
оч |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Среднее время пребывания заявки в системе: |
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
ТСМО |
|
LСМО |
|
|
Lоч |
|
Lобс |
. |
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(3.3.1)
(3.3.2)
(3.3.3)
(3.3.4)
(3.3.5)
Пример 3. На автомойку в среднем за час приезжают 12 автомобилей, но если в очереди уже находятся 4 автомобиля, вновь подъезжающие клиенты не встают в очередь. Среднее время мойки автомобиля составляет 15 мин, а мест для мойки всего 2. Дать характеристику СМО.
|
|
|
|
|
|
|
|
|
|
|
|
Решение. Пусть |
12 авт / ч; tобс |
15 мин; m = 4; n = 2. |
|
||||||||
Определяем интенсивность нагрузки: |
|
|
|
|
|
15 |
|
|
|||
|
|
*tобс |
12* |
3. |
|||||||
|
|
60 |
|||||||||
|
|
|
|
|
|
|
|
|
|
17
Определим долю времени простоя автомойки:
|
|
|
2 |
3k |
|
32 |
1 |
|
|
|
|
3 4 |
1 |
|
|
||||||
|
p0 |
* 1 |
|
|
|
0,016. |
|||||||||||||||
|
k 0 |
к! |
|
2!(2 |
3) |
|
|
2 |
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
Вероятность отказа в обслуживании: pотк |
32 |
|
32 |
4 |
*0,016 |
0,36. |
|||||||||||||||
4 |
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
24 * 2! |
|
|
||||||
Относительная пропускная способность: Q |
1 |
|
pотк |
1 0,36 |
0,64. |
||||||||||||||||
Абсолютная пропускная способность: А |
* Q 12* 0,64 |
7,68авт / ч . |
|||||||||||||||||||
Среднее число автомобилей в очереди: |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
3 |
4 |
|
|
|
|
3 |
|
|
|
|
|
|
|
|
|
|
|
32 1 |
1 |
|
|
* 4 1 4 * |
|
|
|
|
|
|
|
|
|
|
||||||
|
2 |
2 |
|
|
|
|
|
|
|
|
|
||||||||||
L |
|
|
* |
|
|
|
|
|
|
|
|
|
|
* 0,016 |
2,62 авт. |
||||||
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
оч |
2 * 2! |
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1
n
Среднее число заявок, находящихся в обслуживании:
Lобс |
* Q |
3 * 0,64 |
1,92 авт. |
||
Среднее время ожидания в очереди: |
Tоч |
2,619 |
0,218ч. 13,1 мин. |
||
|
|||||
12 |
|||||
|
|
|
|
Среднее время пребывания автомобиля на мойке:
ТСМО |
2,619 |
1,92 |
0,378ч. 22,7 мин. |
|
|
|
|
||
12 |
|
|
||
|
|
|
|
4. Планирование и управление запасами
Фирмы часто делают различные запасы. Хранятся сырье, заготовки, готовая продукция, предназначенная для продажи.
Запасов не должно быть ни слишком много, ни слишком мало. В первом случае возникает необходимость неоправданных затрат на хранение, на амортизацию товара. Во втором случае может оказаться так, что на складе не будет нужного товара. Кроме того, малое количество запасов подразумевает их частое пополнение, что также требует затрат.
Задача управления запасами состоит в том, чтобы избежать обеих крайностей и сделать общие затраты по возможности меньше. Мы рассмотрим несколько моделей управления запасами.
4.1. Основная модель задачи управления запасами
Функция изменения запасов - это связь между количеством единиц товара на складе (Q) и временем t. Будем считать, что имеется один вид товара.
Если на товар имеется спрос, то функция изменения запаса Q = Q(t) убывает. Если товар, наоборот, завозят на склад, то эта функция возрастает.
Затраты, связанные с запасами, можно разделить на три вида:
1)стоимость товара;
2)организационные издержки, т.е. расходы связанные с оформлением товара, его доставкой, разгрузкой и т.д.;
3)издержки на хранение товара, т.е. затраты на аренду склада, амортизацию в процессе хранения и т.д.
18
Рассмотрим основные величины, которыми будем в дальнейшем пользоваться для решения задач.
1.Цена единицы товара – с.
2.Интенсивность спроса – d единиц товара в год.
3.Организационные издержки – s за одну партию товара. Организационные издержки не зависят от размера поставки, т.е. от количества единиц товара в одной партии.
4.Издержки на хранение запаса – h за единицу товара в год. Издержки хранения постоянны.
5.Размер одной партии товара постоянен – q единиц. Партия поступает мгновенно
втот момент, когда возникает дефицит, т.е. когда на складе запас становится равным нулю.
Параметры c, d, s, h считаются заданными. Задача управления запасами состоит в выборе параметра g, таким образом, чтобы минимизировать годовые затраты.
Для решения задачи определим:
1.Общая стоимость товара в год равна произведению годовой интенсивности
спроса (d) и цены единицы товара (с): c |
d . |
|
|
|
|
|
|
|
||||
2. Организационные издержки: |
d |
|
s , d - годовой спрос, q – объем товара в одной |
|||||||||
q |
||||||||||||
|
|
|
|
|
|
|
|
|
|
|||
партии. |
|
|
|
|
|
|
|
|
|
|
|
|
3. Общие издержки на хранение составят: |
q |
|
h , где |
q |
- средний уровень запаса. |
|||||||
2 |
2 |
|||||||||||
|
|
|
|
|
|
|
|
|
||||
Таким образом, общие издержки С рассчитаем по формуле |
||||||||||||
С с |
d |
|
s d |
|
|
|
q h |
. |
(4.1.1) |
|||
|
q |
|
|
2 |
||||||||
|
|
|
|
|
|
|
|
|
Для нахождения оптимального размера партии (q*) найдем производную от
функции С и приравняем к нулю, получим: |
|
|
|
|
||
|
s d |
|
h |
0 . |
(4.1.2) |
|
|
|
|
|
|
||
|
q |
2 |
2 |
|||
|
|
|
||||
|
|
|
|
|
|
Отсюда можно найти q*. Имеем:
|
|
|
|
q* |
2 s d / h . |
(4.1.3) |
Полученная формула называется формулой оптимального запаса или формулой Харриса.
Пример 1. Пусть интенсивность равномерного спроса составляет 2 000 единиц товара в год. Организационные издержки равны 20 у.е., издержки на хранение – 8 у.е на единицу товара в год, цена товара - 2 у.е.
Определить оптимальный размер партии, оптимальное число поставок за год, продолжительность цикла изменения запаса.
Решение. Имеем: d = 2 000, s = 20, h = 8, c = 2.
Общие затраты равны: С с d |
s d |
|
q h |
|
2 |
2 000 |
20 2 000 |
8 q |
. |
|||||||
q |
|
2 |
|
|
q |
|
2 |
|||||||||
|
|
|
|
|
|
|
|
|
|
|||||||
Тогда С (q) |
40 |
|
000 |
4 , |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
q |
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
а оптимальный размер поставки q* является решением уравнения |
|
|
||||||||||||||
|
|
|
|
|
40 |
000 |
4 |
0 , |
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
q2 |
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
19
|
|
|
|
|
|
|
|
|
|
|
|
|
т.е. q* 10 000 100 . |
|
|
|
|
|
|
|
|
|
|
||
Оптимальное число поставок найдем по формуле |
n* |
d |
2 000 |
20. |
||||||||
|
|
|
|
|
|
|
||||||
q * |
100 |
|
||||||||||
|
|
|
|
|
|
|
||||||
Продолжительность цикла изменения запаса равна: t* |
365 |
365 |
18 |
дней. |
||||||||
|
|
|
|
|
|
|
||||||
|
n * |
20 |
|
|||||||||
|
|
|
|
|
|
|
|
4.2.Модель производственных поставок
Вотличие от основной модели товары поставляются с работающей производственной линии. В качестве дополнительного параметра будем рассматривать
производительность производственной линии (p).
Издержки на хранение вычисляются следующим образом. Пусть – время поставки. В течение этого времени происходит как пополнение (с интенсивностью р), так и расходование (с интенсивностью p–d) запаса. Поэтому достигнутый к концу периода пополнения запаса максимальный его уровень М вычисляется по формуле
|
М |
|
( p |
d) |
|
|
|
|
|
|
|
|
|
(4.2.1) |
||
(заметим, что М < q). Однако p *t |
q |
(за время t |
при интенсивности производства р |
|||||||||||||
произведено q единиц товара). Из последних двух равенств следует, что |
|
|
||||||||||||||
|
М |
|
( p |
d ) |
|
q |
. |
|
|
|
|
|
|
|
(4.2.2) |
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
p |
|
|
|
|
|
|
|
||
Тогда издержки на хранение запаса будут равны |
( p |
d ) |
q h |
, где |
q |
- средний |
||||||||||
|
|
|
|
2 |
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
2 |
|
p |
|
||
уровень запаса. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Таким образом, общие издержки С рассчитаем по формуле |
|
|
|
|
||||||||||||
С с |
d |
|
s d |
|
( p |
|
d ) |
q |
h |
. |
|
(4.2.3) |
||||
|
q |
|
|
|
2 |
|
p |
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
Для нахождения оптимального размера партии (q*) найдем производную от
функции С и приравняем к нулю, получим: |
|
|
|
|
|
||
|
s d |
|
( p |
d ) h |
0 . |
(4.2.4) |
|
|
|
|
|
|
|
||
|
q |
2 |
2 |
p |
|||
|
|
|
|||||
|
|
|
|
|
|
|
Отсюда можно найти q*. Имеем:
|
|
|
|
q* |
2 p s d /(( p d ) h) . |
(4.2.5) |
Пример 2. Интенсивность равномерного спроса составляет 2 000 единиц товара в год. Товар поставляется с конвейера, производительность которого составляет 6 000 единиц в год. Организационные издержки равны 10 у.е., издержки на хранение – 4 у.е., цена единицы товара – 6 у.е.
Определить оптимальный объем партии, оптимальное число поставок за год, продолжительность поставки, продолжительность цикла изменения запаса.
Решение. Имеем: d =2 000, p = 6 000, s = 10, h = 4, c = 6.
Общие затраты равны:
С с d |
s |
d ( p |
d ) |
|
q |
h |
6 2 000 |
10 2 000 |
(6 000 |
2 |
000) 4 q |
. |
||||||||
|
q |
|
|
|
2 |
|
p |
|
|
q |
|
2 |
6 |
000 |
||||||
|
|
|
|
|
|
|
|
|||||||||||||
Тогда С (q) |
|
20 |
|
000 |
4 |
|
, |
|
|
|
|
|
|
|
|
|||||
|
|
q |
2 |
|
|
|
3 |
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
20 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
а оптимальный размер поставки q* является решением уравнения |
|
|
|
|
|
||||||||||||||||||||
|
|
|
20 |
|
000 |
4 |
|
0 |
, |
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
q |
2 |
|
3 |
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
т.е. q* |
20 000* |
3 |
122. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
4 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Оптимальное число поставок найдем по формуле |
n* |
|
|
d |
2 000 |
16. |
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
q * |
122 |
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
Продолжительность цикла изменения запаса равна t* |
365 |
365 |
23 |
дня. |
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
n * |
16 |
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
Продолжительность поставки определим по формуле |
|
q * |
365 |
122 |
365 7 |
||||||||||||||||||||
|
|
|
|
|
|
||||||||||||||||||||
|
p |
6 000 |
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
дней.
4.3. Модель поставок со скидкой
В этой задаче используется модель основной задачи, но рассматривается возможность поставки товара по льготной цене (со скидкой), если размер партии достаточно велик. Если размер партии q не менее заданного числа q0, товар поставляется по цене с0, где с0 < с.
Функция общих издержек С(q) будет иметь следующий вид:
|
|
|
c |
d |
s |
|
|
d |
|
|
q |
h |
, |
|
|
если q |
q0 , |
||||||||||||||
|
|
|
|
|
|
|
q |
|
2 |
|
|||||||||||||||||||||
|
С(q) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(4.3.1) |
|||||||||
|
|
|
|
s |
|
|
d |
|
|
q |
|
h |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
с0 |
d |
|
|
|
|
|
|
|
, |
|
|
если q |
q0 . |
|||||||||||||||
|
|
|
|
|
|
|
|
q |
|
|
|
2 |
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
Функция С(q) в точке q = q0 имеет разрыв. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
Обе функции |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
f (q) |
с |
|
d |
|
s |
|
d |
|
q |
|
h |
|
|
(4.3.2) |
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
q |
|
|
2 |
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
и |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
f (q0 ) |
с0 |
d |
|
s d |
|
|
|
q |
h |
(4.3.3) |
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
q |
|
|
2 |
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
имеют минимум в точке, где |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
f (q) |
|
f 0 (q) 0 , |
|
|
|
|
|
|
|
|
|
|
|||||||||||
т.е. в точке q |
2 s d / h . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
Для выяснения вопроса о том, какой размер партии оптимален, следует сравнить |
|||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||
значения функции С(q) в точках q и q , |
и та точка, |
|
где функция С(q) принимает |
||||||||||||||||||||||||||||
меньшее значение, будет оптимальным размером партии q* |
в модели поставок со |
||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||
скидкой. Может случиться так, что С(q) |
C(q) . Тогда, q* |
|
q |
|
q0 . |
||||||||||||||||||||||||||
Пример 3. Интенсивность равномерного спроса составляет |
2 000 единиц в год. |
Организационные издержки равны 15 у.е. , издержки на хранение – 5 у.е. Цена единицы товара равна 6 у.е., однако, если размер партии не менее 1 000 единиц, цена снижается до 4 у.е. Найти оптимальный размер партии.
Решение. Здесь d = 2 000, s = 15, h =5, c = 5, q0 = 1 000, c0 = 2.
Общие издержки определяются функцией С(q):