Метод двойного сглаживания Брауна
Онлайн калькулятор позволит провести прогнозирование временного ряда yt методом двойного сглаживания Брауна.
Инструкция. Укажите количество данных (количество строк), нажмите Далее. Полученное решение сохраняется в файле Word.
Здесь будет отображаться решение.
1. Первое сглаживание: S1t=αYt + (1-α)S1t-1
2. Второе сглаживание: S2t=αS1t + (1-α)S2t-1
3. Уровень: at=2S1t-S2t
4. Тренд: bt=
5. Прогноз на период t+k: yt+k=at+btk
Пример №1.
Угол наклона прямой, образованной 3 точками первой реализации, можно найти с помощью метода наименьших квадратов.
Линейное уравнение тренда имеет вид y = bt + a
Система уравнений МНК:
an + b∑t = ∑y
a∑t + b∑t2 = ∑y*t
| t | y | t2 | y2 | t y |
| 1 | 7266 | 1 | 52794756 | 7266 |
| 2 | 6961 | 4 | 48455521 | 13922 |
| 3 | 6334 | 9 | 40119556 | 19002 |
| 6 | 20561 | 14 | 141369833 | 40190 |
| Ср.знач. | 6853.667 | 4.667 | 47123277.667 | 13396.667 |
3a + 6b = 20561
6a + 14b = 40190
Из первого уравнения выражаем a и подставим во второе уравнение. Получаем a = 7785.667, b = -466
Начальные условия: S1=a+b=7785.6666666667-466=7319.667
S2=
=
=8407
Покажем расчет первых четырех уровней прогноза.
S11=7319.667; S21=7319.667; a1=2*7319.667-8407=6232.333; b1=0.3/(1-0.3)*(7319.667-8407)=-466
S12=0.3*6961+(1-0.3)*7319.667=7212.067
S22=0.3*7212.067+(1-0.3)*8407=8048.52
a2=2*7212.067-8048.52=6375.613
b2=0.3/(1-0.3)*(7212.067-8048.52)=-358.48
S13=0.3*6334+(1-0.3)*7212.067=6948.647
S23=0.3*6948.647+(1-0.3)*8048.52=7718.558
a3=2*6948.647-7718.558=6178.735
b3=0.3/(1-0.3)*(6948.647-7718.558)=-329.962
S14=0.3*6150+(1-0.3)*6948.647=6709.053
S24=0.3*6709.053+(1-0.3)*7718.558=7415.7064
a4=2*6709.053-7415.706=6002.4
b4=0.3/(1-0.3)*(6709.053-7415.706)=-302.8516
| t | yt | Сглаживание, S1t | Сглаживание, S2t | Уровень, at | Тренд, bt | Прогноз, yt* (Ft) | (yt-yt*2) |
| 1 | 7266 | 7319.667 | 8407 | 6232.333 | -466 | 5766.333 | - |
| 2 | 6961 | 7212.067 | 8048.52 | 6375.613 | -358.48 | 6017.133 | 342677.55 |
| 3 | 6334 | 6948.647 | 7718.558 | 6178.735 | -329.962 | 5848.773 | 24107.117 |
| 4 | 6150 | 6709.053 | 7415.706 | 6002.399 | -302.852 | 5699.547 | 21786.075 |
| 5 | 5178 | 6249.737 | 7065.916 | 5433.558 | -349.791 | 5083.767 | 65309.99 |
| 6 | 4871 | 5836.116 | 6696.976 | 4975.256 | -368.94 | 4606.316 | 10869.312 |
| 7 | 4791 | 5522.581 | 6344.657 | 4700.505 | -352.318 | 4348.187 | 8189.367 |
| 8 | 4294 | 5154.007 | 5987.462 | 4320.551 | -357.195 | 3963.356 | 704.976 |
| 9 | 4270 | 4888.805 | 5657.865 | 4119.745 | -329.597 | 3790.147 | 22576.698 |
| 10 | 3611 | 4505.463 | 5312.144 | 3698.782 | -345.72 | 3353.062 | 7705.714 |
| 11 | 3357 | 4160.924 | 4966.778 | 3355.07 | -345.366 | 3009.704 | 3.724 |
| 12 | 3259 | 3890.347 | 4643.849 | 3136.845 | -322.929 | 2813.916 | 14921.828 |
| 13 | 2983 | 3618.143 | 4336.137 | 2900.149 | -307.712 | 2592.437 | 6864.343 |
| 14 | 2830 | 3381.7 | 4049.806 | 2713.594 | -286.331 | 2427.263 | 13550.343 |
| 15 | 2732 | 3186.79 | 3790.901 | 2582.679 | -258.905 | 2323.774 | 22296.811 |
| 16 | 2725 | 3048.253 | 3568.107 | 2528.399 | -222.794 | 2305.605 | 38651.844 |
| 17 | 2700 | 2943.777 | 3380.808 | 2506.746 | -187.299 | 2319.447 | 37346.969 |
| 637562.663 |
Прогнозирование
Первое значение суммы остатков не учитываем, т.к. при однократном экспоненциальном сглаживании первый прогноз возможен только со второго шага (k=16).y(17+1) = 2506.746-187.299 = 2319.447
Ошибка прогноза определяется по формуле:
Определим значение критерия Стьюдента для числа степеней свободы k = n-m = 16-2 = 14 и уровня значимости 0.05: t(14;0.05) = 2.145
Нижняя граница прогноза:
y1=2319.447-2.145*213.401 = 1861.701
Нижняя граница прогноза:
y2=2319.447+2.145*213.401 = 2777.193
Анализ и выводы
Ряд демонстрирует устойчивый отрицательный тренд (убывание). Начальное значение Y1=7266, конечное Y17=2700. Трендовая компонента Tt изначально резко отрицательна, но постепенно становится менее крутой. Это указывает на замедление темпов снижения — ряд продолжает убывать, но с меньшей скоростью. При прогнозировании следует учитывать ошибку прогноза s=457.75.Пример №2. Имеются данные о стоимости произведенной продукции (О) за десять месяцев, а также стоимости основных производственных фондов (Ф) за двенадцать месяцев текущего года.
- Выбрать факторный и результативный признаки. Произвести графический анализ данных. Выбрать приемлемую модель, произвести её спецификацию.
- Определить МНК-оценку параметров модели, выяснить их значимость, а также уравнения в целом.
- Методом экстраполяции линейного периода спрогнозировать стоимость произведенной продукции за ноябрь декабрь месяцы.
- Методами корреляционно-регрессионного анализа, а также адаптивного сглаживания (методом Брауна) спрогнозировать стоимость произведенной продукции за ноябрь и декабрь.
- Проанализировать результаты прогнозирования
| t | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 | 12 |
| О | 2 | 6 | 10 | 11 | 14 | 18 | 20 | 25 | 23 | 24 | ||
| Ф | 15 | 19 | 22 | 20 | 25 | 23 | 24 | 28 | 30 | 31 | 34 | 35 |
Алгоритм адаптации
- Вычислить ошибку прогноза: et=Yt−Ft
- Обновить оценку дисперсии ошибок (например, экспоненциально взвешенную): Qt=γe2t+(1−γ)Qt-1, где γ∈(0,1) — параметр адаптации.
- Вычислить текущий αₜ: αt=|et|/sqrt(Qt) (или другая функция, например, αt=min(0.9,max(0.1,k⋅|et|))
- Использовать αₜ для обновления S1t и S2t вместо фиксированного α.