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

Использование экономико-математических методов и моделей в логистике

.pdf
Скачиваний:
0
Добавлен:
24.11.2025
Размер:
2.19 Mб
Скачать

 

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

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]