数学建模社区-数学中国

标题: ARIMA模型 指数平滑法 [打印本页]

作者: 2744557306    时间: 2024-3-20 10:30
标题: ARIMA模型 指数平滑法
当涉及到时间序列数据分析和预测时,指数平滑法是一种常用的方法之一。它是一种简单而有效的平滑技术,可以用于处理具有趋势和季节性的时间序列数据。
8 v% x4 K% ^& U5 d指数平滑法基于以下假设:8 N! s# r' z; G: L
1.过去的观测值对于预测未来的值更为重要。
. P3 R3 ~. I9 o4 v3 A3 j- s' L2.近期的观测值对于预测的影响程度比过去的观测值更大。4 _8 E8 _) E" R* H8 e
3.观测值中的随机波动应该被消除,使得预测值更稳定。
) T, S0 \  |5 ~3 a5 [4 a7 \' D* |2 G: Y: l# y1 h
指数平滑法有几种不同的变体,其中最常见的是简单指数平滑、二次指数平滑和三次指数平滑。
! ?0 k, T6 x, c+ \' \7 Y+ F! Y0 W% v1 |( v9 x* o
4.简单指数平滑(Simple Exponential Smoothing):
- A& Y5 N# z' Q/ i! Q; U简单指数平滑方法将所有观测值视为等权重,并使用一个平滑参数(通常表示为α)来确定各个观测值的权重。具体来说,简单指数平滑法的计算公式如下:. z/ `( S4 R. W6 A9 d, M/ \, E6 Y
* Y: c% _. B' f' E
5.初始值:S₁ = X₁! U7 G" b- W' N) q+ P+ H" R$ e; E/ a& s
6.平滑公式:Sₜ = α * Xₜ + (1 - α) * Sₜ₋₁
' d  X! t+ h3 t8 }2 K其中,Sₜ表示时刻 t 的平滑值,Xₜ表示时刻 t 的观测值。, P: x' |% J: E3 }7 ?
7.二次指数平滑(Double Exponential Smoothing):' l* V* B, V' |
二次指数平滑法相对于简单指数平滑法,将趋势也纳入考虑。它不仅考虑过去观测值的权重,还考虑了过去观测值的趋势。具体来说,二次指数平滑法的计算公式如下:
5 E7 ]6 ?+ y0 t& e; u& u. O
- D2 ]. V. g' T% b$ D6 Z8.初始水平值:L₁ = X₁
5 w9 }9 n0 B- T, O" R9.初始趋势值:T₁ = X₂ - X₁
9 Q! y9 F0 Y- \4 [5 V# L10.平滑公式(水平值):Lₜ = α * Xₜ + (1 - α) * (Lₜ₋₁ + Tₜ₋₁)
. @3 g9 j5 G+ E$ X+ M5 f# k11.平滑公式(趋势值):Tₜ = β * (Lₜ - Lₜ₋₁) + (1 - β) * Tₜ₋₁
* p6 z; i( ^; }3 ^其中,Lₜ表示时刻 t 的水平值,Tₜ表示时刻 t 的趋势值,Xₜ表示时刻 t 的观测值,α和β是平滑参数。
" m4 N+ T$ H0 N  o1 W, C* R6 L12.三次指数平滑(Triple Exponential Smoothing):" n3 z% C6 v& I- ]" `
三次指数平滑法是在二次指数平滑法的基础上进一步考虑了季节性的影响。它适用于具有明显季节性模式的时间序列数据。具体来说,三次指数平滑法的计算公式如下:
( d6 U0 y  c* ~* e9 ]9 [$ p3 G
8 d/ ?% N- Z+ g, F1 t13.初始水平值:L₁ = X₁- h" U* g" Q$ b6 Z+ V( T7 H
14.初始趋势值:T₁ = X₂ - X₁
) y7 m1 ~4 B3 k& }3 m5 }$ N15.初始季节性值:S₁ = X₃ - 2T₁ - L₁) {! M0 D% a4 G- t# _4 o
16.平滑公式(水平值):Lₜ = α * (Xₜ - Sₜ₋₃) + (1 - α) * (Lₜ₋₁ + Tₜ₋₁)
% \1 \# ~5 k9 G17.平滑公式(趋势值):Tₜ = β * (Lₜ - Lₜ₋₁) + (1 - β) * Tₜ₋₁. y- a) T* b/ x6 n  h
18.平滑公式(季节性值):Sₜ = γ * (Xₜ - Lₜ - Tₜ) + (1 - γ) * Sₜ₋₃- u9 h' I' w2 E; M% P) x* N
其中,Lₜ表示时刻 t 的水平值,Tₜ表示时刻 t 的趋势值,Sₜ表示时刻 t 的季节性值,Xₜ表示时刻 t 的观测值,α、β和γ是平滑参数。8 [/ _. n- z/ i& W! w1 X5 }

9 B  k& g3 c) j% K- X4 Z请注意,指数平滑法适用于简单的时间序列模式,例如渐变性趋势或没有明显季节性模式的数据。对于更复杂的模式,如长期趋势、季节性和其他周期性模式,可以考虑使用更高级的模型,例如ARIMA模型。
+ C8 F8 ~) k5 h) K' o+ a3 BARIMA(Autoregressive Integrated Moving Average)模型是一种常用于时间序列分析和预测的统计模型,它将自回归(AR)、差分(I)和移动平均(MA)三种方法结合在一起。ARIMA模型根据时间序列数据的自相关性和移动平均性来拟合模型,并用于预测未来的数值。ARIMA模型通常用于处理具有复杂模式的时间序列数据,如具有长期趋势、季节性和周期性模式的数据。2 l6 ]- J0 p- V; I5 W% h
2 Y5 B) Z; ]: N& Y; c; |! _% X) L
# 导入所需的库) ]5 N4 A: ~  C0 y9 s, x
import numpy as np
- T% B1 \2 E% W2 dimport pandas as pd/ p/ R% S* ^8 ~% q" W: g( n* T

" h* x) c0 u$ K& a2 X, a# 源数据
+ w: P; h* R$ C/ xdf = pd.DataFrame({. ^+ U9 [! G. c  r% W
    't': [i for i in range(1, 11)],1 n% _, g4 {/ }4 I
    'production': [2031, 2234, 2566, 2820, 3006,
3 U. d# W5 x0 @# P) n                   3093, 3277, 3514, 3770, 4107],
' V5 I# A0 L& G' w: P9 C& h+ q( u4 d})1 h5 Y8 z6 O) i+ ~) [% P
6 a  i6 u0 B1 R0 r  D7 L/ j
这段代码创建了一个DataFrame df,包含了时间(t)和产量(production)的数据。5 I' i! L* E& n' T- g* S9 M5 f
# 设 alpha=.3,计算一次、二次指数平滑
3 v  g/ W, W& h" r7 e5 Q' d- q" @$ Lalpha = .3! e5 F! l/ f1 }, E5 P& b2 G, T
s1, s2 = [int(df['production'][0])],\' p3 K0 V* j4 A9 x7 g& Y+ H% k
    [int(df['production'][0])]+ E# |% A7 b$ ?" w
/ d; Y: m; w$ T: v! U
for i in range(1, len(df['t'])):
" I# K( F" F' p( u, g# v+ r    s1.append(alpha*df['production'][i] + (1-alpha)*s1[i-1])
( @% L3 L% j/ A! F    s2.append(alpha*s1[i] + (1-alpha)*s2[i-1])
; `( I$ E! q. adf['s1'] = s12 m5 v7 V* R4 Z; ^4 [) a+ h+ R
df['s2'] = s20 l3 I' t4 X6 e9 E& U0 M7 Y

2 v, P# r9 A/ B  W4 b这段代码使用指数平滑法计算了一次指数平滑和二次指数平滑的值,并将结果添加到DataFrame df中的's1'和's2'列中。$ R. L2 o; I; o: d
# 计算过去年的预测值,以及未来年的线性表达式0 j% Z5 A& R6 P
predict_list = [None]
! J' j- _6 t2 T1 U6 c& Z, xfor i in range(len(df['t'])-1):
. t# e& w1 x: C& ^) J6 R9 O    a = 2*df['s1'][i] - df['s2'][i], O6 o6 W& P- K
    b = (alpha / (1-alpha)) * (df['s1'][i] - df['s2'][i])/ k, E& L0 z, t- p! _( J- B+ ^& v
    predict_list.append(a + b)
( t: k+ s1 ^3 R/ h4 ]
, s- l* }8 c, \% x3 At = 107 Y0 o3 c" `# g9 `
a = 2*df['s1'][t-1] - df['s2'][t-1]: z8 u( Q: |1 C
b = (alpha / (1-alpha)) * (df['s1'][t-1] - df['s2'][t-1])
% j6 f- c+ n1 X0 T( P3 e/ Kdf['predict'] = predict_list
" _( S2 B) r* F8 ]1 Mprint('at =', a)
% W+ `* N# L% c5 `- ?, c7 |" \print('bt =', b)
' V8 N; l' F  \- b* |$ s/ U
+ ?& Z7 a5 y# \" V这段代码计算了过去年份的预测值,并给出了未来年份的线性表达式。其中a和b是线性表达式的参数,对应于二次指数平滑法的结果。将过去年份的预测值添加到DataFrame df的'predict'列中,并打印出a和b的值。: E0 r  @% Q  a8 D2 l4 S
# 计算未来年的预测值' u! t# J7 h* y' }3 ~) _  V8 `- w
pred_df = df.copy()
' ?5 J, d* T$ @/ }# l" Y+ ~for i in range(5):
  L- g1 J- w& N0 m$ L, W% E    pred_df = pd.concat([pred_df, pd.DataFrame({7 s' H$ _( h* a$ b2 i% z, M9 h0 f# z
        't': [10+i+1],
" {1 z3 X0 w# G) Q        'predict': a + b*i,. g* j2 ~: R; j' _+ y) v, z+ ]$ c/ u$ b3 E
    })])
: x/ ?( _% }$ C- u. d% ?% ?& hpred_df
7 @: T& h7 i7 U9 p4 K' U: r( j
) f! i) N) n7 h  b% y7 E+ o这段代码使用线性表达式计算了未来5年的预测值,并将结果添加到新的DataFrame pred_df中。
. c9 y* c/ e  T- S1 u# 画图* \/ A& J9 m/ D! S
import matplotlib.pyplot as plt
9 v" d# A& S) F/ v, q, h/ \- p- C) ]
plt.plot(pred_df['t'], pred_df['production'],label='production')
$ p4 V9 N& k, uplt.scatter(pred_df['t'], pred_df['production'])
. v7 b. [0 f: c1 Wplt.plot(pred_df['t'], pred_df['s1'],label='s1')  Z2 Y1 g; Y. U& h$ |
plt.scatter(pred_df['t'], pred_df['s1'])
* ?2 s* {! K5 y% k0 T# nplt.plot(pred_df['t'], pred_df['s2'],label='s2')
# k  m1 \  T1 Y- w* W2 zplt.scatter(pred_df['t'], pred_df['s2'])
  p! c/ q/ X+ q$ J7 F- hplt.plot(pred_df['t'], pred_df['predict'],label='predict'): ^' h- ?1 l" l' k- u+ H! u
plt.scatter(pred_df['t'], pred_df['predict'])
) I' m) L7 f' g  ~7 u7 \plt.legend()
* ]( N' o8 g% F! \$ j/ \3 {+ d0 `
0 M- p0 v- c" P- e6 V5 o8 q, @; Z这段代码使用matplotlib库将实际产量、一次指数平滑、二次指数平滑和预测值的曲线图绘制出来。+ |( `& ]2 k4 y1 t! O) {( q
希望以上解释对你有帮助。如果你还有任何问题,请随时提问。
- u3 f  o) ~- j' l5 o/ M5 S+ z: Q( S: r8 V& _$ c6 v  I/ e

" v2 t/ C/ u/ M0 U! H5 L

24.exponential_smoothing.py

1.64 KB, 下载次数: 1, 下载积分: 体力 -2 点

售价: 2 点体力  [记录]






欢迎光临 数学建模社区-数学中国 (http://www.madio.net/) Powered by Discuz! X2.5