Использование экономико-математических методов и моделей в логистике
.pdf
|
q s |
|
C2 |
|
|
|
||||
|
P(0) |
|
|
|
; |
|
||||
|
C |
C |
2 |
|
||||||
|
q 0 |
1 |
|
|
|
|
||||
Ss |
P 0 d 2 |
1 |
0,714; |
|||||||
|
d 1 |
|
||||||||
|
|
|
|
|
|
|
|
|
||
0,91s 1 |
0,09·0,714 0,91s |
0,286; |
||||||||
|
0,09 |
|
|
|
|
|
|
|
||
S |
ln 0,286 |
13,27 т. |
|
|||||||
ln 0,91 |
|
|||||||||
|
|
|
|
|
|
|
|
|
||
Принимаем, что s = 14 т. Значит, для хранения товара необходим склад вместимостью 14 т.
Найдем суммарныеиздержки при данномобъемезапаса на складе:
P(q) P(0)·d q 1; d 1 P(0);
C S |
|
q S |
S |
q |
·P 0 ·d q 1 C |
· |
q |
|
q S ·P 0 ·d q 1. |
||||||
C · |
|
|
|
||||||||||||
|
|
1 |
q 0 |
|
|
|
|
|
2 |
q S 1 |
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
q |
|
|
|
q S |
|
|
|
|
|
|
|
||
Так как |
|
d q 1 |
1 d q 1, то |
|
|
|
|
|
|
||||||
|
|
q S 1 |
|
|
|
q 0 |
|
|
|
|
|
|
|
||
|
|
|
|
|
|
C S P 0 C1 C2 |
|
|
|||||||
|
q S |
|
|
|
|
|
q S |
|
|
|
|
|
2 |
|
|
|
|
|
P 0 |
q 1 |
|
|
|
q 1 |
|
·P 0 1 S , |
|||||
S· |
1 |
|
|
q 1 P 0 |
|
|
C |
||||||||
|
q 0 |
|
|
|
|
|
q 0 |
|
|
|
|
|
|
|
|
50
C S 0,09 100 250
|
q 14 |
0,91 q 1 |
q 14 |
|
250·0,09 13 2252,39 у. е. |
|
14· |
|
q 0,91 q 1 |
|
|||
|
|
q 0 |
|
q 0 |
|
|
|
|
|
|
|
||
6.2. Определение параметров складской системы с учетом затрат на хранение товара
Спрос на хранение на складе товара имеет случайный характер. Вероятность распределения данного спроса отражает геометрическая прогрессия. Затраты на хранение единицы товара в единицу времени Сх составляют 20 у. е. / т · дн. затраты, связанные с неудовлетворением спросанаединицутоваравединицувремениСу – 50 у. е.
Необходимо определить размер товара S, закладываемого на склад, с тем условием, чтобы суммарные издержки на хранение и издержки, связанные снеудовлетворением спроса, были минимальны.
Так как спрос на товар имеет случайный характер, необходимо рассматривать два случая:
1. Спрос за время τ не превысит размера хранимого на складе товара (рисунок 6.1).
Рисунок 6.1 – Графическая интерпретация определения параметров складской системы в случае, если спрос не превышает размера хранимого на кладе товара
51
S q;
S S S q S q . 2 2
2. Спрос за время τ превысит размер хранимого на складе товара (рисунок 6.2).
Рисунок 6.2 – Графическая интерпретация определения параметров складской системы в случае, если спрос не превышает размера хранимого на кладе товара
|
|
|
|
|
|
|
|
|
|
|
S q; |
|
|
|
|
S |
|
|
|
S 2 |
|
|
|
||
|
|
|
|
|
|
|
|
|||||
так как |
1 |
|
|
, то S |
|
|
|
. |
|
|
||
q |
2 q |
|
|
|||||||||
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
(q S) 2 |
; |
|
|
|
|
|
|
|
|
S |
|||||
|
|
|
|
|
|
|
|
|
|
|
2 |
|
52
|
|
|
|
|
2 |
q S |
|
|
|
|
|
|
(q S)2 |
|
|
|
|
|
|
|
|
||||||||
так как |
, то |
SH |
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
q |
|
|
2 q |
|
. |
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
Найдем математическое ожидание суммарных издержек: |
|||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
q S |
|
|
|
|
|
q |
|
|
|
|
|
||||
|
|
|
|
|
|
|
М С(S) |
Cx · S |
2 |
·P q |
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
q 0 |
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
q |
|
|
S |
2 |
|
Cy · q S |
2 |
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
Cx · |
|
|
|
|
|
·P q |
min; |
|||||||||||||||
|
|
|
|
|
2q |
|
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2q |
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
q S 1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
M |
|
C(S) |
|
|
q S |
|
|
|
|
|
q |
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
Cx · |
P |
q |
|
|
|
Cx ·S |
Cy · q S |
·P q 0; |
||||||||||||||||
|
|
S |
|
|
|
|
|||||||||||||||||||||||
|
|
|
|
|
|
q 0 |
|
|
|
|
q S 1 |
q |
|
|
|
|
q |
|
|||||||||||
|
Cx |
· P q S Cx Cy · |
P q |
Cy · |
P q 0; |
||||||||||||||||||||||||
|
|
|
|
q S |
|
|
|
|
|
|
|
|
|
|
|
q |
|
|
|
|
|
|
q |
|
|||||
|
|
|
|
q 0 |
|
|
|
|
|
|
|
|
|
|
|
q S 1 |
q |
|
|
|
q S 1 |
|
|||||||
|
|
|
|
|
|
|
|
|
q |
P q |
q S |
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
P q 1; |
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
q S 1 |
|
|
|
|
q 0 |
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
q |
|
|
|
|
|
|
|
|
|
|
|
|
|
q S |
|
|
|
|||
|
|
|
|
|
|
|
Cy |
P q |
Cy Cy P q ; |
|
|||||||||||||||||||
|
|
|
|
Cx |
|
|
|
q S 1 |
|
|
|
|
|
|
|
|
|
|
q 0 |
q Cy ; |
|||||||||
|
|
|
|
Cy P q S Cx Cy |
P |
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
q S |
|
|
|
|
|
|
|
|
|
|
|
|
|
q |
|
|
|
||||
|
|
|
|
|
|
|
|
q 0 |
|
|
|
|
|
|
|
|
|
|
|
|
q S 1 |
q |
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Cy |
|
|
; |
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Cx Cy |
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
q S |
|
|
|
|
|
|
q |
P |
q |
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
P |
q S |
|
|
|
|
|
|
. |
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
q |
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
q 0 |
|
|
|
|
|
|
q S 1 |
|
|
|
|
|
|
|
|
||||||
53
Далее необходимо произвести расчет, когда
q S |
q |
P q |
|
|
P q S |
|
|
. |
|
q |
||||
q 0 |
q S 1 |
|
По итогу расчетов неравенство будет иметь решение при S ≥ 7.
6.3. Определение параметров складской системы на основе массового обслуживания
6.3.1. Расчет системы массового обслуживания с потерями
Пусть имеется простейший поток требований, который удовлетворяет следующим условиям:
– ординарности (в малый отрезок времени не может поступить сразу два требования):
P1( t) P 1( t);
P0 ( t) P1( t) 1 – за время t может поступить одно или ни
одного требования
P1( t) – вероятность того, что за время t поступит ровно одно требование;
–стационарности (вероятности не зависят от расположения данного интервала на временной оси – λ = const);
–отсутствию последействия (вероятность времени обслуживания канала не зависит от того, сколько он обслуживался до этого).
Простейший поток требований описывается законом Пуассона:
Pk (t) ( kt!)k e t
(вероятность поступления k требований за время t;
P0 ( t) e t P(T t) 1 P(T t)
(не поступило ни одного требования);
54
P(T t) 1 e t
(вероятность того, что интервал между поступлениями двух смежных требований будет меньше t),
где T – интервал времени между поступлениями двух смежных требований.
Пусть имеется n-канальная система. Пусть время обслуживания одного требования имеет показательный закон распределения:
F (t) 1 e t ;
M1T
где v – интенсивность обслуживания требования (сколько требований имеется возможность обслужить за единицу времени).
Будем рассматривать систему массового обслуживания в момент времени t и дадим бесконечно малое приращение t. Состояние системы массового обслуживания определяется количеством требований, находящихся в системе.
Обозначим состояние системы: X – некоторое состояние системы;
Xk – в системе находится k требований.
На рисунке 6.3 представлен граф состояний.
Рисунок 6.3 – Граф состояний системы массового обслуживания в момент времени t
t + t – система будет находится в состоянии А.
55
Событие А: в момент времени t система находится в состоянии Xk–1. Событие В: ни одно из требований k-1 не будет обслужено за время t. Событие С: за это же время поступило ровно одно требование.
1. Рассчитаем вероятность события А:
Пусть Pk(t) – вероятность того, что в системе находиться k требований в момент времени t.
P ( t) 1 P ( t) 1 e t 1 (1 t) t; |
|
|
|
|||||||||||||||||||||||
1 |
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
k n f k (x |
) |
|
( x)k 0( x)k 1 |
|
|
|
|
|
||||||||||||
f (x |
|
x) |
|
0 |
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
0 |
|
|
|
k 0 |
k! |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
(ряд Тейлора). |
|
|
|
|
|
|
|
а x x, тогда |
|
|
|
|
|
|
|
|
|
|
|
|||||||
Предположим, что x0 0, |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
n xk |
|
|
n e0 |
|
xk |
n 1 xk |
|
|
x |
|
|
x2 |
|
|
|
|
|
|||||||||
ex |
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
... |
|
|
|
||||
|
k! |
k! |
k! |
1! |
2! |
|
|
|
||||||||||||||||||
n 0 |
|
n 0 |
n 0 |
|
|
|
|
|
|
|
|
|||||||||||||||
P(T t) 1 P(T t) e t – |
вероятность |
того, что |
время |
|||||||||||||||||||||||
обслуживания одного требования будет больше |
|
|
t. |
|
|
|
|
|
|
|
||||||||||||||||
Вероятность того, что за время |
t не будет обслужено ни одного |
|||||||||||||||||||||||||
k–1 требования, будет определяться по следующей формуле: |
|
|
|
|||||||||||||||||||||||
P |
1 |
( t) e (k l) t |
1 (k 1) t. |
|
|
|
|
|
||||||||||||||||||
k |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
В формулу Тейлора подставляем вместо x (–λ· |
t), тогда |
|
x2 |
0, |
||||||||||||||||||||||
|
2! |
|||||||||||||||||||||||||
т. е. 2 ( t)2 0 ex 1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
x |
1 x, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
2! |
|
|
|
|
|
1! |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
т. е.
e t 1 ( t) 1 t,
P(A) Pk 1(t) t 1 (k 1) t Pk 1(t) t.
56
2. Рассчитаем вероятность события B:
P(B) P1(t) e t e t Pk (t) 1 ( k ) t .
3. Рассчитаем вероятность события С, вероятность того, что ровно одно изk–1 находящихсяв системетребованийбудет обслужено:
P(T t) 1 e t t;
Ck1 1 t e k t (k 1) t(1 k t) (k 1) t;
P(C) Pk 1(t)(1 t)(k 1) t (k 1) t Pk 1(t);
где P(T t) – вероятность того, что время обслуживания будет меньше t;
Ck1 1 – число сочетаний из k + 1 по одному.
Таким образом, вероятность того, что система находиться в состоянии Xk в момент времени t + t равна:
Pk (t t) P(A) P(B) P(C),
Pk (t t) Pk 1(t) t Pk (t) 1 ( k ) t (k 1) t Pk 1(t),
Pk (t t) Pk (t) Pk 1(t) ( k ) Pk (t) (k 1) Pk 1(t),
t
дPдkt(t) Pk 1(t) ( k ) Pk (t) (k 1) Pk 1(t),
дPд0t(t) P0 (t) P1(t),
57
дPдnt(t) Pn 1(t) n Pn (t),
где n – число каналов обслуживания.
Решив систему уравнений, можно найти вероятности для любого состояния системы. Однако на практике используют установившийся режим функционирования системы массового обслуживания,
т. е. такой режим, когда t .
Таким образом, для установившегося режима Pk (t) const т. е. не зависит от времени t, и производная Pk (t) равна нулю, система преобразуется к следующему виду:
0 P0 P1; |
|
|
|
|
0 Pn 1 n Pn ; |
|
|
|
k 1 P |
|
|
|
0 Pk 1 k P |
. |
|
|
k |
k 1 |
|
Сделаем следующую подстановку:
Zk Pk 1 k Pk ;Z0 0;
Zk Zk 1 0;
Zk 0;
Pk k Pk 1;
P1 P0 ;
P |
|
P |
1 |
2 |
P . |
|
|
|
|
|
|||
2 |
2 |
1 |
2 |
0 |
||
|
|
|
|
|
||
58
Пусть , тогда
|
|
|
|
P |
|
k |
P ; |
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
k |
|
k! |
0 |
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
||
k n |
|
|
|
k n |
k |
|
|
|
|
|
|
1 |
|
|
P |
1 |
|
|
|
P |
1 |
P |
|
|
|
, |
|||
k! |
|
k |
||||||||||||
k 0 |
k |
|
|
k 0 |
0 |
|
0 |
|
k n |
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
k! |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
k 0 |
|
|
где Pn – вероятность того, что все каналы будут заняты, т. е. веро-
ятность отказа поступившего требования на обслуживание; А – вероятность обслуживания;
Pобслуж 1 Pотказа – вероятность того, что одно требование, по-
ступившее в систему будет обслужено. Абсолютная пропускная способность системы:
NPобслуж (1 Pn ).
6.3.2.Определение параметров складской системы
Q 1530 т – годовой грузооборот склада;
T 365 дн;
q 24,23 т – средний размер партий груза;
1 4,28 дн;
d 0,55 т/м2 – средняя нагрузка груза на площадку склада.
Необходимо определить площадь склада, практически обеспечивающую хранение грузооборота. Вероятность обслуживания долж-
на быть не менее 0,95(Pобслуж 0,95) .
|
Q |
|
1530 |
0,173 партий/день; |
|
q T |
24,23 365 |
||||
|
|
|
59
