<< Пред. стр. 5 (из 10) След. >>
1112
15
22
27
31
39
45
50 11
12
13
14
15
16
17
18
19
20 0
6
5
7
10
7
6
9
*
* 0
6
11
18
28
35
41
50
50
50 21
22
23
24
25
26
27
28
29
30 5
5
4
3
4
1
2
8
3
4 5
10
14
17
21
22
24
32
35
39 Примечание: * - дефицит
Для оценки границ интервального прогноза необходимо рассчитать среднее квадратичное отклонение
(7.5)
Вспомогательные расчеты приведены в табл. 7.2. Подставляя значения в формулу (7.5), находим ?t
?t =
На основании полученных зависимостей yi и ?t раcсчитываются прогнозные оценки
* среднего времени расхода текущего запаса ;
* страхового запаса ус с заданной доверительной вероятностью Р;
* вероятности отсутствия дефицита деталей на складе в течение прогнозируемого периода.
Расчет прогнозной величины среднего времени расхода производится по формуле (7.2). Приняв yt = 0, находим
дней
Для расчета страхового запаса воспользуемся формулой
(7.6)
где ?t - среднее квадратичное отклонение, формула (7.5);
t? - параметр нормального закона распределения, соответствующий доверительной вероятности ?.
Параметр t? определяет для нормального закона число средних квадратических отклонений, которые нужно отложить от центра рассеивания (влево и вправо) для того, чтобы вероятность попадания в полученный участок была равна ?.
В нашем случае доверительные интервалы откладывают вверх и вниз от среднего значения yt.
На рис. 7.2. приведены границы, соответствующие yt ± ?.
В табл. 7.3. приведены наиболее часто встречающиеся в практических расчетах значения вероятности ? и параметра t? для нормального закона распределения.
Таким образом, страховой запас рассчитывается также, как и границы интервального прогноза, т.е. по формуле (7.6).
Таблица 7.2.
Исходные данные и результаты расчета
коэффициентов уравнения (7.2) при N = 5
ti, дни yi, ед. ti2 yiti прогноз (yt - yi)2 1
2
3
4
5 41
39
38
35
28 1
4
9
16
25 41
78
114
140
140 42
39
36
33
30 1
0
4
4
4 Суммы
?Примечание: значения округлены
Для рассматриваемого примера для доверительной вероятности ?=0,9 находим по табл. 7.3 t? = 1,643 и по формуле (7.5) величину страхового запаса
yc = 1,8?1,643 = 2,96.
Примем yc = 3,0.
На рис. 7.2. приведены границы интервального прогноза при ?=0,9.
Рассчитанное значение страхового запаса соответствует только одному дню наступления дефицита, а именно согласно прогнозу Т=15. Для учета возможных нарушений срока поставки, необходимо также при расчете страхового запаса оценить влияние задержки, связанной с выполнением заказа, в частности, с транспортировкой.
Рис. 7.2. Прогноз текущего расхода деталей на складе (N = 5)
1-исходные данные; 2-уравнение тренда; 3, 3'-границы интервального прогноза; 4-время расхода запаса Т.
К сожалению, по одной реализации невозможно оценить вероятностный характер длительности функциональных циклов поставки. Однако можно предположить, что выявленная тенденция расхода запаса, формула (7.2), сохраниться. В этом случае, для оценки прогнозной величины страхового запаса можно воспользоваться формулой
yc? = |a1|?+t??t (7.7)
где ? - параметр, характеризующий количество дней задержки поставки заказа.
Таблица 7.3.
Доверительная вероятность ? и параметр t?
нормального закона распределения
? t? ? t? 0,80
0,82
0,84
0,86
0,88
0,90
0,91 1,282
1,340
1,404
1,475
1,554
1,643
1,694 0,92
0,94
0,95
0,96
0,98
0,99
0,999 1,750
1,880
1,960
2,053
2,325
2,576
3,290
Рассчитаем величину страхового запаса при условии задержки на один день по сравнению с прогнозной оценкой Т=15 дней, т.е. на 16 день. По формуле (7.7) находим
yc? = |-3,0|1,0+1,643?1,8=6,0 ед.
Аналогично, при ?=2 (17 день) yc?=9,0 ед.
Для оценки вероятности отсутствия дефицита допустим, что отклонения ежедневного расхода деталей от среднего значения (тренда) подчиняются нормальному закону распределения. Тогда, воспользовавшись уравнением функции нормального закона, определим вероятность отсутствия дефицита по формуле
Р(y) = 1 - F(y) = 1 - , (7.8)
где yt - уравнение тренда, формула (7.2)
? - среднее квадратическое отклонение, формула (7.5).
Таблица 7.4.
Значения? нормальной функции
распределения Ф (x), вероятности Р(x) и параметра x
x Ф (x) Р(x) x Ф (x) Р(x) 0,00
-0,125
-0,253
-0,385
-0,525
-0,675
-0,842
-1,037 0,50
0,45
0,40
0,35
0,30
0,25
0,20
0,15 0,50
0,55
0,60
0,65
0,70
0,75
0,80
0,85 -1,280
-1,405
-1,555
-1,645
-1,75
-2,05
-2,30
-3,10 0,10
0,08
0,06
0,05
0,04
0,02
0,01
0,001 0,90
0,92
0,94
0,95
0,96
0,98
0,99
0,999 ? некоторые значения округлены
Сделаем в интеграле замену переменной
(7.9)
и приведем его к виду
F(x) = , (7.10)
Интеграл (7.10) не выражается через элементарные функции, поэтому для расчетов можно воспользоваться численными методами и ЭВМ или специальными таблицами. Для нормальной функции распределения с параметрами среднее значение mx=0 и ?x=1.
Ф(x) = (7.11)
Очевидно, что F(y) = Ф().
В табл. 7.4. приведен ряд значений функции Ф(x) и Р(x).
Между параметрами ? и x, а также ? и Ф(x) существует соотношение
2Ф(x) - 1 = ? (7.12)
На рис. 7.3 приведены графики нормальной функции распределения и плотности нормального распределения.
Появление дефицита означает, что текущая величина запаса на складе равна нулю, т. е. y = 0.
Следовательно, для определения вероятности отсутствия дефицита необходимо:
* по формуле (7.9) рассчитать x=;
* по табл. 7.4 с помощью x найти Р(x).
Для рассматриваемого примера рассчитаем вероятности отсутствия дефицита деталей на складе на 13, 14 и 15 день. Так, для t=13 получаем
,
и
x = .
По табл. 7.4 находим Р(Т=13)?0,999, т.е. вероятность дефицита ничтожно мала.
Аналогично, для Т=14 получим ; x = -1,78, и вероятность отсутствия дефицита РТ=14 ? 0,95.
Наконец, для Т=15, вероятность отсутствия дефицита Р ? 0,5.
Следует подчеркнуть, что также как при оценке прогнозной величины страхового запаса, определение вероятности отсутствия дефицита по одной реализации справедливо только при строгом соблюдении сроков поставки. Если они не соблюдаются, то расчет должен проводиться с учетом рассеивания длительности функциональных циклов поставки.
а) б)
Рис. 7.3. Нормальный закон распределения:
а) плотность распределения; б) функция распределения.
В заключении определим ошибку прогноза среднего времени Т, поскольку имеются реальные данные о текущем расходе в табл. 7.1.
?Т = ?100% (7.13)
где Тф, Тп - соответственно фактическая и прогнозная продолжительность цикла, дн.
Подставив значения в (7.13), находим
?Т = ?100%=50%
Ошибка прогноза велика, но это закономерно, так как нарушено одно из эмпирических правил экстраполяционного прогнозирования: между предпрогнозным периодом t и периодом упреждения (прогноза) ? =Т- t должно соблюдаться соотношение:
(7.14)
Если следовать соотношению (7.14), то при t = 5 допустимая величина времени прогноза
Т= (7.15)
Следовательно, величина надежного прогноза соответствует Т?7 дней и период упреждения составляет ?=2 дня.
Пример 2. В работе [2] указывается, что средняя длина функционального цикла расхода запасов составляет =10 дней. Тогда, по формуле (7.15) находим t=7,5 дней.
Увеличим длину динамического ряда до N=7 и, выполнив аналогичные расчеты (табл. 7.5), получим уравнение тренда
y?t=47-3,9t
Соответственно, ?t?=2,1.
Рассчитаем среднее прогнозное время расхода запаса со склада
= дней,
и ошибку прогноза, формула (7.13)
?Т = ?100%=20%
Таблица 7.5
Исходные данные и результаты расчета
коэффициентов уравнения (7.2) при N=7
ti yi ti2 yiti (-yi)2 1 41 1 41 43,1 4,41 2 39 4 78 39,2 0,04 3 38 9 114 35,3 7,29 4 35 16 140 31,4 12,96 5 28 25 140 27,6 0,25 6 23 36 138 23,6 0,36 7 19 49 133 19,7 0,49 Суммы