数学建模社区-数学中国
标题:
时间序列模型 （四）：差分指数平滑法、 自适应滤波法
[打印本页]
作者:
浅夏110
时间:
2020-5-31 14:31
标题:
时间序列模型 （四）：差分指数平滑法、 自适应滤波法
差分指数平滑法
$ v0 }+ o$ G. h/ p* _9 ~
在上节我们已经讲过，当时间序列的变动具有直线趋势时，用一次指数平滑法会出现滞后偏差，其原因在于数据不满足模型要求。因此，我们也可以从数据变换的角度来考虑改进措施，即在运用指数平滑法以前先对数据作一些技术上的处理，使之能适合于 一次指数平滑模型，以后再对输出结果作技术上的返回处理，使之恢复为原变量的形态。 差分方法是改变数据变动趋势的简易方法。下面我们讨论如何用差分方法来改进指数平 滑法。
' }- C: K0 S; Y7 p) t
4 c1 r# }# H, Q, {& c3 G
一阶差分指数平滑法
`8 R( L/ Y/ C' G2 Y4 ]( b
4 U: S ]) K1 U1 M: `6 O
! y/ i) q5 a+ g+ P$ K2 B5 c5 G' r2 {
( p/ Q @# D# }1 f% O
在前面我们已分析过，指数平滑值实际上是一种加权平均数。因此把序列中逐期增 量的加权平均数（指数平滑值）加上当前值的实际数进行预测，比一次指数平滑法只用 变量以往取值的加权平均数作为下一期的预测更合理。从而使预测值始终围绕实际值上 下波动，从根本上解决了在有直线增长趋势的情况下，用一次指数平滑法所得出的结果 始终落后于实际值的问题。
l* H: r$ p% `5 N+ U( L* ~$ |
4 J' S1 R6 j7 K1 g! [( R0 t
3 M# J1 b9 a. Y; m' f* b9 l) t# W
! U, A0 A- _ S: W; Y
例 7 某工业企业 1977～1986 年锅炉燃料消耗量资料如表 8 所示，试预测 1987 年 的燃料消耗量。
& w+ G( M7 q" k
G1 r, f; v& T1 X
! V/ O! f1 V- [2 o; ~! z2 K& J
8 I( J: M+ E/ l/ ~
解 由资料可以看出，燃料消耗量，除个别年份外，逐期增长量大体在 200 吨左右，即呈直线增长，因此可用一阶差分指数平滑模型来预测。我们取 α = 0.4，初始值为新序列首项值，计算结果列于表 8 中。预测 1987 年燃料消耗量为
( n2 a Q% h: H1 p) F2 I+ J
: ?) |. z& q! ]5 L; J, \$ k
( 百吨 )
i4 J9 Y& V* I) i1 W
* S6 R- o& g: _6 t0 t2 O
0 R7 B( B- b' {* D" n) Z
7 |/ h1 X# k: o) O
二阶差分指数平滑模型
5 E" q9 H/ P" Y0 y& N
# b' ^8 L+ {! b; G7 E$ R: w9 J3 C
. _' N n& l% t8 n! r' T, T
2 _: a- y# M( H; N
2 e( |$ m2 S( c/ d& F6 f* i
差分方法和指数平滑法的联合运用，除了能克服一次指数平滑法的滞后偏差之外， 对初始值的问题也有显著的改进。因为数据经过差分处理后，所产生的新序列基本上是平稳的。这时，初始值取新序列的第一期数据对于未来预测值不会有多大影响。其次， 它拓展了指数平滑法的适用范围，使一些原来需要运用配合直线趋势模型处理的情况可用这种组合模型来取代。但是，于指数平滑法存在的加权系数 α 的选择问题，以及只能逐期预测问题，差分指数平滑模型也没有改进。
5 J, M% h# J; U6 A E
?! z3 z7 I* M }; T
自适应滤波法
4 W" s5 K, t; K; g( c s2 K
自适应滤波法的基本过程
# a2 [4 Q9 i5 }+ }. z
自适应滤波法与移动平均法、指数平滑法一样，也是以时间序列的历史观测值进行 某种加权平均来预测的，它要寻找一组“佳”的权数，其办法是先用一组给定的权数 来计算一个预测值，然后计算预测误差，再根据预测误差调整权数以减少误差。这样反 复进行，直至找出一组“佳”权数，使误差减少到低限度。由于这种调整权数的过 程与通讯工程中的传输噪声过滤过程极为接近，故称为自适应滤波法。
4 x3 \8 D: D, ]
2 f2 ^1 g9 w' [0 Y4 t6 ]8 O
; j- U+ \) k* I+ D. y
3 ^0 x, m( Z7 j: q
6 z- Z- g7 g6 B5 h
2 X. u/ U3 ^$ Z; p
下面举一个简单的例子来说明此法的全过程。设有一个时间序列包括 10 个观测值， 如表 9 所示。试用自适应滤波法，以两个权数来求第 11 期的预测值。
- @5 H; J' f# R
8 {* o1 U1 t3 y* Z6 h% c% L' r
: g& l) R( z+ e" t% i* C
3 @- B- n; r! j i8 ~! J' r
O1 H9 g# h9 ^! R$ q* V0 }
2 u. H k* n+ q+ m9 \$ L# m
在实际应用中，权数调整计算工作量可能很大，必须借助于计算机才能实现。
8 X9 r; ?5 h. E: R
; f, n+ ~0 L+ R* E4 n6 l8 n) z
计算的 MATLAB 程序如下：
$ u: a! N8 f" B, B
0 t$ a* s ~7 f5 H. C' H& I
clc,clear
( Z; z$ }3 v) M( p0 q& y1 O
yt=0.1:0.1:1;
6 X9 s. p8 O5 V/ T: [3 V
m=length(yt);
6 ~* V: f9 B6 E5 `6 o6 ~
k=0.9; N=2; Terr=10000;
5 C' p0 Z9 Y6 v) u
w=ones(1,N)/N;
/ L: F7 P9 A2 M; t3 w
while abs(Terr)>0.00001
6 p' W. [1 P. \
Terr=[];
* y6 B+ _7 Z* a. l6 y
for j=N+1:m-1
# g. M, M& |, h7 F- g( M& d
yhat(j)=w*yt(j-1:-1:j-N)';
$ H/ u- j7 B3 Z7 |
err=yt(j)-yhat(j);
- t* Z/ W$ D) y+ `
Terr=[Terr,abs(err)];
% g8 ?" n4 j; |; F
w=w+2*k*err*yt(j-1:-1:j-N);
6 Y' ] Q$ L- H" A# f
end
# s0 }0 L- {6 Z
Terr=max(Terr);
6 H0 I+ G0 c5 e& g% F* ]; ^4 e
end
( }: P* p3 E0 _& X$ x
w, yhat
7 [% B9 R7 i; @* y5 ]( ~9 r! o+ ?
% M1 o6 f( k% L! @% O
N, k 值和初始权数的确定
; j% d/ k8 x$ p/ N) Y5 b
在开始调整权数时，首先要确定权数个数 N 和学习常数k 。一般说来，当时间序 列的观测值呈季节变动时， N 应取季节性长度值。如序列以一年为周期进行季节变动 时，若数据是月度的，则取 N =12 。若季节是季度的，则取N=4 。如果时间序列无 明显的周期变动，则可用自相关系数法来确定，即取 N 为高自相关系数的滞后时期。
& k9 x2 |4 R7 f
O+ C% Y5 W9 [7 {5 {& b
k 的取值一般可定为1 / N，也可以用不同的k 值来进行计算，以确定一个能使S 小的k 值。
0 g9 h' Z; m! ]3 B: R
) X2 j4 N. u! N
, q8 A' L. I2 q3 G6 x Q2 W2 o$ h
0 O% n, H- N- V! E% b5 V% f
自适应滤波法有两个明显的优点：一是技术比较简单，可根据预测意图来选择权数 的个数和学习常数，以控制预测。也可以由计算机自动选定。二是它使用了全部历史数 据来寻求佳权系数，并随数据轨迹的变化而不断更新权数，从而不断改进预测。 由于自适应滤波法的预测模型简单，又可以在计算机上对数据进行处理，所以这种 预测方法应用较为广泛.
) u! \$ R$ r& b9 j; G: ^
9 D- n9 L9 M2 V c+ [% r3 N
: F) {7 t( @5 \9 H& ~: a/ s
————————————————
5 f* O' O& m- O3 n, |. ?/ p- ~/ d' M
版权声明：本文为CSDN博主「wamg潇潇」的原创文章，遵循CC 4.0 BY-SA版权协议，转载请附上原文出处链接及本声明。
) X7 a3 f2 \ h5 K4 Y) E- P5 ^
原文链接：https://blog.csdn.net/qq_29831163/article/details/89447556
) r; j' }- F, M/ G; V0 \
0 ~& J1 }6 Z2 c+ G" m$ G& ]4 z) ?
( ^- v. k" Z, \% V' B1 m* f
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5
; j- U+ \) k* I+ D. y
: g& l) R( z+ e" t% i* C
, q8 A' L. I2 q3 G6 x Q2 W2 o$ h