QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 3312|回复: 0
打印 上一主题 下一主题

[建模教程] 时间序列模型 (二):移动平均法

[复制链接]
字体大小: 正常 放大
浅夏110 实名认证       

542

主题

15

听众

1万

积分

  • TA的每日心情
    开心
    2020-11-14 17:15
  • 签到天数: 74 天

    [LV.6]常住居民II

    邮箱绑定达人

    群组2019美赛冲刺课程

    群组站长地区赛培训

    群组2019考研数学 桃子老师

    群组2018教师培训(呼伦贝

    群组2019考研数学 站长系列

    跳转到指定楼层
    1#
    发表于 2020-5-30 15:04 |只看该作者 |倒序浏览
    |招呼Ta 关注Ta |邮箱已经成功绑定
    移动平均法 可以作为一种数据平滑的方式 ,以每天的气温数据为例,今天的天气可能与过去的十天的气温有线性关系;或者有的人对食物有一种节俭的美德,他们做的饭菜能看出有些是上一顿的,当然也有一部分是今天的做的,再假设隔两顿的都被倒掉了,并且每天都是这样的,那么这碗饭菜可能就是一部分上一顿的再加上一部分今天现做的,这就是一个一阶的移动平均。
    7 I8 A( c. ]6 E6 @8 V7 w) \( o: }3 E: Y# ]9 u1 b2 P& U9 \9 @
    移动平均法; t" {% F" }+ |5 [5 H/ \  P
    移动平均法是根据时间序列资料逐渐推移,依次计算包含一定项数的时序平均数, 以反映长期趋势的方法。当时间序列的数值由于受周期变动和不规则变动的影响,起伏 较大,不易显示出发展趋势时,可用移动平均法,消除这些因素的影响,分析、预测序列的长期趋势。  移动平均法有简单移动平均法,加权移动平均法,趋势移动平均法等。
      h. g8 ?0 {1 ~: V  j) Q9 c& r8 A
    ! [2 s1 h  t  {8 g简单移动平均法
    * X' ^5 F! ]4 G' ^3 t* I8 f1 Z" Z" q# c
    $ E" V( u9 Z  f

    6 H( d% P) m; Y近N 期序列值的平均值作为未来各期的预测结果。一般 N 的取值范围: 5≤N≤ 200。当历史序列的基本趋势变化不大且序列中随机变动成分较多时,N 的 取值应较大一些。否则 N 的取值应小一些。在有确定的季节变动周期的资料中,移动平均的项数应取周期长度。选择佳 N 值的一个有效方法是,比较若干模型的预测误 差。预测标准误差小者为好。
    1 H! W2 p9 ^/ F2 E8 |9 q5 a( @+ f* q' \% |% c
    简单移动平均法只适合做近期预测,而且是预测目标的发展趋势变化不大的情况。 如果目标的发展趋势存在其它的变化,采用简单移动平均法就会产生较大的预测偏差和滞后。 0 E, `9 i" ]; u6 a9 x6 }

    7 p$ }2 b! c9 |! m2 d0 ^9 ~例 1  某企业 1 月~11 月份的销售收入时间序列如表 1 示。试用一次简单滑动平 均法预测第 12 月份的销售收入。
    / R/ ]  x2 U3 ]3 ^
    1 G: S% c7 p0 S: n& Y1 v. b: K3 p* e; h6 E/ n8 g6 e! O* w$ x. w
    . F( Z1 Q/ p9 h
    3 o# ~+ m% f5 ]- Z# V
    ' q: \1 B: T) a4 Q$ _/ ~" P9 |
    计算的 Matlab 程序如下: 0 r* ?% R/ Y& @: s$ ~1 c8 W, v" G

    * i" R# F5 s. y& Y( Qclc,clear 4 g3 }6 Z4 @7 R. J9 F
    y=[533.8  574.6  606.9  649.8   705.1  772.0  816.4  892.7  963.9  1015.1  1102.7]; m=length(y);   9 `7 V( I, U) p- R) {( I
    n=[4,5];   %n 为移动平均的项数 6 I- e$ U# B) X3 W8 `+ }" D
    for i=1:length(n)   
    ! ^2 A. X- W7 C6 h1 v %由于 n 的取值不同,yhat 的长度不一致,下面使用了细胞数组    / R0 k9 K$ T" d2 a9 b1 Z
        for j=1:m-n(i)+1         8 h8 K2 J- @) g7 H4 N- \( u
            yhat{i}(j)=sum(y(j:j+n(i)-1))/n(i);     
    - }3 w* \. e/ c6 z4 [+ Q* W    end   
    - E5 Z0 g) Z9 B9 u    y12(i)=yhat{i}(end);     1 ^. X6 j2 Q/ r$ D. M" o/ |8 W
        s(i)=sqrt(mean((y(n(i)+1:m)-yhat{i}(1:end-1)).^2));
    # W& Q( I) U* g' w* K: Zend
    ; ?/ w' {9 M- U& ^4 Q4 W) my12,s
    $ x$ a0 G$ g  R7 F
    $ I- a  p3 g4 X4 D3 \3 V加权移动平均法7 d/ u3 w7 g  _% b
    在简单移动平均公式中,每期数据在求平均时的作用是等同的。但是,每期数据所包含的信息量不一样,近期数据包含着更多关于未来情况的信心。因此,把各期数据等同看待是不尽合理的,应考虑各期数据的重要性,对近期数据给予较大的权重,这就 是加权移动平均法的基本思想。 / z5 Z1 N* K! ^- ~4 F+ w
    ( K$ E8 b; A5 t* F! _1 a( |  e
    ( c, ?. U  B. ]1 o5 k8 B( ^

    # r3 i8 i( X9 I5 m! R2 G- P例 2  我国 1979~1988 年原煤产量如表 2 所示,试用加权移动平均法预测 1989 年 的产量
    , `7 q* y4 P& v! B4 _% r
    ; }( W7 b, D1 i. f% p1 I! _" e
    " _* |' y9 s6 }; `- b
    9 n# a# I  u3 n, K, w  ?# I% ~- E. y
    & @% k% O6 w2 k* U
    : s  k- m. J6 |: p) a+ |6 h1 u  B* [0 j* A4 D
    : c: O+ D$ ^7 i! k

    1 s1 H, R2 U. L! K2 k% X计算的 MATLAB 程序如下:
    7 R( \* E( i' g) o* A  f* Y$ j, M
    9 H& G: V* @! }8 S) |y=[6.35 6.20    6.22    6.66    7.15    7.89    8.72    8.94    9.28    9.8]; 0 F6 T: h8 w0 t
    w=[1/6;2/6;3/6];
    1 {' {3 s8 k; W* G$ o/ u7 [$ km=length(y);n=3; / Z0 Q+ z" X% s" ?/ D9 O2 l
    for i=1:m-n+1     # U& ?' i& T6 \% W' L9 Z
        yhat(i)=y(i:i+n-1)*w; % Y' L# }' h( i6 x9 Q
    end 8 @" Y$ s) G) ?
    yhat ; ~( B; ~; W! ?% ^1 W: R
    err=abs(y(n+1:m)-yhat(1:end-1))./y(n+1:m) ! _, H# q! a4 t4 p
    T_err=1-sum(yhat(1:end-1))/sum(y(n+1:m)) ( ^* y. A( `  I) c. M; B0 T0 v
    y1989=yhat(end)/(1-T_err)
    " j7 @1 \3 c0 h1 a, t+ W  V; `, Y2 G! S4 N& V3 K
    在加权移动平均法中,   的选择,同样具有一定的经验性。一般的原则是:近期 数据的权数大,远期数据的权数小。至于大到什么程度和小到什么程度,则需要按照预 测者对序列的了解和分析来确定。
    . M( g' X! |" a: i
    & c. c+ ^$ t  [! h7 y趋势移动平均法2 N5 c4 d, V; j% o' a6 i* h
    简单移动平均法和加权移动平均法,在时间序列没有明显的趋势变动时,能够准确 反映实际情况。但当时间序列出现直线增加或减少的变动趋势时,用简单移动平均法和 加权移动平均法来预测就会出现滞后偏差。因此,需要进行修正,修正的方法是作二次 移动平均,利用移动平均滞后偏差的规律来建立直线趋势的预测模型。这就是趋势移动平均法。  一次移动的平均数为
    ; |- \% ^: B1 h2 U$ g! x2 i  A7 ]' R$ O( U' J/ G- W& Z9 F: O6 e$ Y5 |9 o

    & x) P' G* m# Z" g+ C7 |
    . p% ~7 Z8 k. D+ l% g9 \! X5 E! {, T9 q8 f, Z& s$ Z- Q
    8 y) K5 u3 T  e

    9 ^6 @5 L" X$ E8 z3 g' y1 V1 {/ D% b: m0 h  O, G3 A
    例 3  我国 1965~1985 年的发电总量如表 3 所示,试预测 1986 年和 1987 年的发 电总量。 1 l0 a! z4 f( X# B( _

    - E) j, ?5 Q8 }- G9 A: Z  \8 r9 B5 u2 W; t8 d5 O

    * f. f# y& y/ D3 D. U8 K解  由散点图 1 可以看出,发电总量基本呈直线上升趋势,可用趋势移动平均法 来预测。 : J* g* G6 C# R  h, I: N
    8 G' u7 m% Y1 r0 ?! m( H7 H
    % G' @# U. S0 ]+ i% F  q7 h
    1 Y4 V+ I* x, l  P; [& o! _3 J0 I
    计算的 MATLAB 程序如下: 9 J) Q& J& P4 y. f
      E5 W& D) Q8 S2 v
    clc,clear
    . G& E! X3 ~7 n% @" dload y.txt   %把原始数据保存在纯文本文件 y.txt 中 % k0 n( u. Y8 I! Q& ^% d
    m1=length(y);    & D  K$ S# D: n! F/ M; \4 l
    n=6;   %n 为移动平均的项数
    3 r4 @( A: d! H! Q! K' qfor i=1:m1-n+1     
    : I- N8 ^) i# F0 [! p# _2 t    yhat1(i)=sum(y(i:i+n-1))/n; 7 d$ o, C- }7 [: F
    end 0 M8 }9 i$ V5 a2 V
    yhat1
    / c* A# p  [4 _& x. T8 pm2=length(yhat1);
    ; r! d$ f5 U+ J( nfor i=1:m2-n+1   
    9 r. G( T! q. K5 h: d. A) ?    yhat2(i)=sum(yhat1(i:i+n-1))/n; 6 A2 p$ F2 z& u  \: Y
    end
    ; m& E/ Y# ~+ y. X4 H/ i. @( Y: Iyhat2   2 w( B4 h: `& h. X5 e, H
    plot(1:21,y,'*')
    & B, v6 u+ m3 b1 D( \! D. q' \# aa21=2*yhat1(end)-yhat2(end) # r. R* j* d0 l
    b21=2*(yhat1(end)-yhat2(end))/(n-1) % i, r1 |+ ~5 M) g# b1 l% T
    y1986=a21+b21
    ( K& Q. X9 V7 `5 n0 By1987=a21+2*b21 5 m4 A7 T" y+ F8 t& \

    ! M: @7 w! Y' J7 u
    7 d9 K8 r+ [4 B( M( Z2 O
    * N. m7 p& Z& |趋势移动平均法对于同时存在直线趋势与周期波动的序列,是一种既能反映趋势变 化,又可以有效地分离出来周期变动的方法。( p' c7 \' Y$ L, T$ V

    % Z" w* F$ }) S6 T* p
    $ o, E' ]8 ]- |, o' o6 O
    / b/ P% |9 E1 c2 ^! e2 N  e1 a' R1 O, {————————————————
    / F% Q9 N  {' `0 R- L0 d4 t版权声明:本文为CSDN博主「wamg潇潇」的原创文章,遵循CC 4.0 BY-SA版权协议,转载请附上原文出处链接及本声明。
    8 I" ^6 C4 \% k% ^原文链接:https://blog.csdn.net/qq_29831163/article/details/89440426
    & J9 i8 m2 u+ a! P' a" t* _" k( Z# W2 m7 y0 L

    / Y5 G1 m0 L5 w% D% S
    zan
    转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信
    您需要登录后才可以回帖 登录 | 注册地址

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

    关于我们| 联系我们| 诚征英才| 对外合作| 产品服务| QQ

    手机版|Archiver| |繁體中文 手机客户端  

    蒙公网安备 15010502000194号

    Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问:王兆丰

    GMT+8, 2026-4-14 17:58 , Processed in 0.295635 second(s), 51 queries .

    回顶部