где yt прогнозные значения временного ряда; среднее значение прогноза (тренд); vt составляющая прогноза, отражающая сезонные колебания (сезонная волна); εt случайная величина отклонения прогноза.
В частных случаях количество составляющих модели меньше, например, только и vt. Подробно вопросы прогнозирования с использованием методов экстраполяции изложены в ряде работ, но ввиду отсутствия общепринятого алгоритма обработки временных рядов может быть предложена следующая последовательность расчета:
3.4.2. Примеры прогноза
|
1-й цикл | 2-й цикл | 3-й цикл | ||||||
---|---|---|---|---|---|---|---|---|
День | Спрос, ед. |
Всего с начала цикла |
День | Спрос, ед. |
Всего с начала цикла |
День | Спрос, ед. |
Всего с начала цикла |
1 | 9 | 9 | 11 | 0 | 0 | 21 | 5 | 5 |
2 | 2 | 11 | 12 | 6 | 6 | 22 | 5 | 10 |
3 | 1 | 12 | 13 | 5 | 11 | 23 | 4 | 14 |
4 | 3 | 15 | 14 | 7 | 18 | 24 | 3 | 17 |
5 | 7 | 22 | 15 | 10 | 28 | 25 | 4 | 21 |
6 | 5 | 27 | 16 | 7 | 35 | 26 | 1 | 22 |
7 | 4 | 31 | 17 | 6 | 41 | 27 | 2 | 24 |
8 | 8 | 39 | 18 | 9 | 50 | 28 | 8 | 32 |
9 | 6 | 45 | 19 | * | 50 | 29 | 3 | 35 |
10 | 5 | 50 | 20 | * | 50 | 30 | 4 | 39 |
Проиллюстрируем возможные варианты прогнозов для одной реализации.
Пример 3.4.1.
Воспользуемся первой реализацией. Допустим, что нам известны значения расхода деталей со склада за пять дней работы (табл. 3.4.2).
Таблица 3.4.2
ti, дн. | yi, ед. | yiti | Прогноз yi* | (yt–yi)2 | |
---|---|---|---|---|---|
1 | 41 | 1 | 41 | 42 | 1 |
2 | 39 | 4 | 78 | 39 | 0 |
3 | 38 | 9 | 114 | 36 | 4 |
4 | 35 | 16 | 140 | 33 | 4 |
5 | 28 | 25 | 140 | 30 | 4 |
Суммы |
* Значения округлены
Выберем уравнение тренда yt в виде линейной зависимости:
(3.4.2) |
Расчет коэффициентов уравнения и производится по формулам, полученных на основе метода наименьших квадратов:
(3.4.3) |
(3.4.4) |
Находим: a0 = 45,2, a1 = 3,0. Таким образом, уравнение прогноза пишется в виде:
(3.4.5) |
Для оценки границ интервального прогноза необходимо рассчитать среднее квадратичное отклонение σt:
(3.4.6) |
Подставляя значения в формулу, находим σt:
(3.4.7) |
На основании полученных зависимостей yt и σt рассчитываются прогнозные оценки:
Расчет прогнозной величины среднего времени расхода производится по формуле
(3.4.8) |
Приняв yt = 0, находим :
Для расчета страхового запаса воспользуемся формулой:
(3.4.9) |
где σt среднее квадратичное отклонение,
tβ параметр нормального закона распределения, соответствующий доверительной вероятности β.
Параметр tβ определяет для нормального закона число средних квадратических отклонений, которые нужно отложить от центра рассеивания (влево и вправо) для того, чтобы вероятность попадания в полученный участок была равна β.
В нашем случае доверительные интервалы откладывают вверх и вниз от среднего значения уt. .
В табл. 3.4.3 приведены наиболее часто встречающиеся в практических расчетах значения вероятности β и параметра tβ для нормального закона распределения.
Таблица 3.4.3
β | tβ | β | tβ |
---|---|---|---|
0,80 | 1,282 | 0,92 | 1,750 |
0,82 | 1,340 | 0,94 | 1,880 |
0,84 | 1,404 | 0,95 | 1,960 |
0,86 | 1,475 | 0,96 | 2,053 |
0,88 | 1,554 | 0,98 | 2,325 |
0,90 | 1,643 | 0,99 | 2,576 |
0,91 | 1,694 | 0,999 | 3,290 |
Страховой запас рассчитывается так же, как и границы интервального прогноза. Для рассматриваемого примера при доверительной вероятности β=0,9 находим по табл. 3.4.3 tβ = 1,643. Тогда величина страхового запаса составит:
Примем yc=3,0.
На рис. 3.4.2 приведены границы интервального прогноза при β = 0,9.
Рис. 3.4.2. Прогноз текущего расхода деталей на складе (N = 5):
1 исходные данные; 2 уравнение тренда;
3, 3' границы интервального прогноза; 4 время расхода запаса
Рассчитанное значение страхового запаса соответствует только одному дню наступления дефицита, а именно согласно прогнозу T = 15. Для учета возможных нарушений срока поставки необходимо также при расчете страхового запаса оценить влияние задержки, связанной с выполнением заказа, в частности с транспортировкой.
К сожалению, по одной реализации невозможно оценить вероятностный характер длительности функциональных циклов поставки. Однако можно предположить, что выявленная тенденция расхода запаса сохранится. В этом случае для оценки прогнозной величины страхового запаса можно воспользоваться формулой
(3.4.10) |
где τ параметр, характеризующий количество дней задержки поставки заказа.
Рассчитаем величину страхового запаса при условии задержки на один день по сравнению с прогнозной оценкой T = 15 дней, т. е. на 16-й день:
Аналогично, при τ = 2 (17 день)
Для оценки вероятности отсутствия дефицита допускается, что отклонения ежедневного расхода деталей от среднего значения (тренда) подчиняются нормальному закону распределения. Тогда, пользуясь уравнением функции нормального закона, определяют вероятность отсутствия дефицита:
(3.4.11) |
где yt уравнение тренда;
σ среднее квадратическое отклонение.
В табл. 3.4.4 приведен ряд значений функции Ф(х) и Р(х).
x | Ф(х) | Р(х) | x | Ф(х) | Р(х) |
---|---|---|---|---|---|
0,00 | 0,50 | 0,50 | -1,280 | 0,10 | 0,90 |
-0,125 | 0,45 | 0,55 | -1,405 | 0,08 | 0,92 |
-0,253 | 0,40 | 0,60 | -1,555 | 0,06 | 0,94 |
-0385 | 0,35 | 0,65 | -1,645 | 0,05 | 0,95 |
-0,525 | 0,30 | 0,70 | -1,75 | 0,04 | 0,96 |
-0,675 | 0,25 | 0,75 | -2,05 | 0,02 | 0,98 |
-0,842 | 0,20 | 0,80 | -2,30 | 0,01 | 0,99 |
-1,037 | 0,15 | 0,85 | -3,10 | 0,001 | 0,999 |
Появление дефицита означает, что текущая величина запаса на складе равна нулю, т. е. у = 0.
Для определения вероятности отсутствия дефицита необходимо:
Для рассматриваемого примера рассчитаем вероятности отсутствия дефицита деталей на складе на 13-й, 14-й и 15-й дни. Так, для T = 13 получаем:
По табл. 3.4.4 находим РТ=13 > 0,999, т. е, вероятность отсутствия дефицита ничтожно мала.
Аналогично, для T = 14 получим yТ=14 = 3,2, x = 1,78, и вероятность отсутствия дефицита РТ=14 = 0,95.
Наконец, для T = 15 вероятность отсутствия дефицита Р = 0,5.
Следует подчеркнуть, что так же, как при оценке прогнозной величины страхового запаса, определение вероятности отсутствия дефицита по одной реализации справедливо только при строгом соблюдении сроков поставки. Если они не соблюдаются, то расчет должен проводиться с учетом рассеивания длительности функциональных циклов поставки.
В заключение определим ошибку прогноза среднего времени Т:
(3.4.12) |
где Tф, Tп соответственно фактическая и прогнозная продолжительность цикла, дн.
Получим:
Ошибка прогноза велика, но это закономерно, так как нарушено одно из эмпирических правил экстраполяционного прогнозирования: между предпрогнозным периодом t и периодом упреждения (прогноза) τ = T t должно соблюдаться соотношение:
(3.4.13) |
При T = 5 допустимая величина времени прогноза:
Следовательно, величина надежного прогноза соответствует T ≈ 7 дн. и период упреждения составляет τ = 2 дн.
Пример 3.4.2.
Считается, что средняя длина функционального цикла расхода запасов составляет T = 10 дн. Тогда t = 7,5 дн.
Увеличим длину динамического ряда до N = 7 (рис. 3.4.3).
Рис. 3.4.3. Прогноз текущего расхода деталей на складе (N = 7):
1 исходные данные; 2 уравнение тренда;
3, 3' границы интервального прогноза; 4 время расхода запаса
Выполним расчеты аналогично примеру 3.4.1, полученные данные занесем в табл. 3.4.5.
Таблица 3.4.5
ti | yi | yiti | yi | (yt–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 |
Суммы |
Получим уравнение тренда:
Соответственно,
Рассчитаем среднее прогнозное время расхода запаса со склада
Рассчитаем величину страхового запаса yc для 12-го, 13-го и 14-го дней. Примем β = 0,95, т. е, tβ = 1,96. Тогда:
Определим вероятность дефицита на складе на 10-й день. Находим ; по табл. 3.4.4 , т.е. наличие дефицита маловероятно. Аналогично, для для .
В заключение можно сделать следующее замечание: рассчитанные величины среднего запаса получены при условии, что наблюдающая величина дефицита и вариация ежедневного расхода независимые величины. Несомненно, это допущение требует проверки.
На формирование стратегии автотранспортного предприятия (АТП) на рынке влияют факторы как внешней, так и внутренней среды, в том числе определяющие состояние спроса на услуги. Основным является вопрос о потенциальных возможностях предприятия, определяемых технико-технологическими и организационно-финансовыми факторами среды. Принципиальное различие между предъявляемыми к перевозке грузами (или спросом) и провозными возможностями АТП состоит в том, что первое следует отнести к условиям внешней среды, т. е. «природе», состояние которой формируется под действием большого количества факторов и в подавляющем большинстве случаев не зависит от транспортной политики конкретного АТП (если рассматриваемое предприятие не является монополистом в данном сегменте рынка транспортных услуг), а второе определяется политикой и тактикой действий предприятия, не имея случайного характера, а, скорее, подчиняясь неким внутренним закономерностям. Таким образом, под влиянием случайных факторов объем перевозок представляет собой случайную величину, подчиняющуюся определенному закону или функции распределения F(Q). Введение функции распределения для описания состояния «природы» позволяет, согласно теории статистических решений, использовать вероятностные критерии принятия решений в условиях риска.
Что касается состояния АТП, то оно может быть представлено в виде различных стратегий Ai, каждая из которых количественно характеризуется числом автомобилей Ni и их провозными возможностями Wi.
Указанные стратегии Ai, являются дискретными величинами, если используется число автомобилей N, или непрерывными за счет варьирования показателей, входящих в расчет производительности автомобиля Wi.
Связь между Ai стратегией и объемом перевозок Qi определяется в виде матрицы (табл. 3.4.6), элементы которой (aij) отражают «выигрыш», получаемый АТП при выборе i-й стратегии.
Таблица 3.4.6
Стратегия АТП | Объëм перевозок | |||||
---|---|---|---|---|---|---|
Q1 | Q2 | ... | Qj | ... | Qn | |
A1 | a11 | a12 | ... | a1j | ... | a1n |
A2 | a21 | a22 | ... | a2j | ... | a2n |
... | ... | ... | ... | ... | ... | ... |
Ai | ai1 | ai2 | ... | aij | ... | ain |
... | ... | ... | ... | ... | ... | ... |
Am | am1 | am2 | ... | amj | ... | amn |
В ряде работ, где предпринимались попытки использования теории статистических решений для конкретных хозяйственных объектов, в качестве элемента матрицы aij «выигрыша» использовались условные величины. В качестве «выигрыша» могут быть использованы различные экономические показатели: доход, прибыль и другие, а также показатели, способствующие усилению конкурентных или рыночных позиций, усилению влияния на клиентуру и укреплению имиджа предприятия, улучшению качества производимых услуг.
Возможны три соотношения между объемом перевозок Qi и стратегией предприятия А. Первое состояние «выигрыша»; второе нейтральное состояние; третье состояние «проигрыша». Величина вероятностное отклонение за счет случайного характера величин, определяющих значение элементов матрицы. Теоретически возможен вариант, когда области значений aij будут расположены иным образом, чем это показано на рис. 3.4.4.
Рис. 3.4.4. Распределение на различные области матрицы стратегий АТП:
1 «выигрыш», 2 нейтральное состояние, 3 «проигрыш»
В частности, введение оценки «упущенной выгоды» может изменить границы областей 1-3. Считается, что наилучшей стратегией A = Ai является та, при которой показатель Ai обращается в максимум:
(3.4.14) |
где Qj = F(Qj) вероятность j-го состояния «природы».
Таким образом, оптимальная стратегия АТП может быть определена при наличии F(Qj) и матрицы стратегий aij.
Рассмотрим возможные варианты расчета F(Qj). Традиционно для количественной оценки прогноза Qj, используется метод экстраполяции по динамическим рядам с использованием полиномов различной степени.
Результаты прогноза представляются в виде среднего значения Q и дисперсии DQ, по которым определяется вид функции распределения F(Qj). Далее c использованием условия максимизации Ai выбираем стратегию АТП.
Основная трудность использования вышеописанной методики это невысокая точность прогноза, Повышение точности может быть достигнуто за счет комбинированных методов прогноза, предусматривающих синтез двух и более прогнозных вариантов.
Каждый метод прогнозирования обладает определенной достоверностью, имеет свои преимущества и недостатки, Считается, что комбинированные методы прогнозирования (синтез прогнозов) позволяют компенсировать недостатки одних способов достоинствами других. На рис. 3.4.5 представлена блок-схема комбинированного прогноза для двух вариантов прогноза, один из которых прогноз, выполненный эвристическим методом, основанным на статистической обработке мнений экспертов.
Рис. 3.4.5. Блок-схема выбора стратегии АТП
в целевом сегменте рынка транспортных услуг
Процедура получения экспертных оценок может быть формализована и представлена в виде блок-схемы (рис. 3.4.6).
Рис. 3.4.6. Блок-схема прогноза на основе экспертных опросов
Рассмотрим некоторые блоки подробнее.
Формирование группы экспертов важнейшая составляющая экспертного метода. Не останавливаясь подробно на вопросах персонального подбора, затронем только количественную сторону, а именно число экспертов. Известно, что при прогнозировании в целях минимизации расходов на прогноз стремятся привлекать минимальное число экспертов при условии обеспечения ошибки результата прогнозирования не более Е, где 0 < Е < 1. Поэтому рекомендуемое число экспертов может быть определено по формуле:
(3.4.15) |
При подстановке предельных значений Е находим:
|
Таким образом, минимальное количество экспертов равно 4. Для определения максимальной численности экспертной группы используется неравенство:
(3.4.16) |
где Ki компетентность i-го эксперта, рассчитываемая на основе анкеты самооценки;
Kmax максимально возможная компетентность по используемой шкале компетентности экспертов.
Статистический анализ результатов опроса предусматривает проведение двух взаимосвязанных процедур: традиционной статистической обработки в виде средних значений, дисперсий и т. п., а также оценки всей экспертной группы степени согласованности, взаимосвязи и других показателей мнений экспертов. Оценка группы экспертов проводится с использованием части полученных статистических оценок. Если последние не удовлетворяют соответствующим критериям, то в блок-схеме предусмотрена корректировка, которая приводит, в частности, к изменению состава экспертов и повторной процедуре опроса.
Методика статистической обработки данных включает следующие этапы:
(3.4.17) |
где aij ранг, присвоенный j-м экспертом i -му фактору;
m число экспертов.
(3.4.18) |
где k число факторов.
(3.4.19) |
(3.4.20) |
Если W существенно отличается от нуля, то можно полагать, что между оценками экспертов существует определенное согласие.
(3.4.21) |
где табличное значение.
В случае соблюдения неравенства с доверительной вероятностью можно утверждать, что мнения экспертов относительно вероятности факторов согласуются неслучайно.
Представленный вариант получения прогноза на основе экспертных оценок является универсальным и в случае использования баллов заканчивается построением ранжированной диаграммы рангов.
Для перехода к конкретному прогнозу, в частности объема перевозок, последовательности расчета сводятся к следующему:
Полагаем, что интервалу Qj соответствует наиболее правдоподобная гипотеза (П1), затем вторая ( (П2) и т.д.
(3.4.22) |
Значения весовых коэффициентов для определения комбинированных оценок вероятностей каждого интервала находим по формулам:
(3.4.23) |
где μ1 и Dq весовой коэффициент и дисперсия экстраполяционного прогноза;
μ2 и Dэq весовой коэффициент и дисперсия экспертного прогноза.
(3.4.24) |
|
|