QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 3363|回复: 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 |邮箱已经成功绑定
    移动平均法 可以作为一种数据平滑的方式 ,以每天的气温数据为例,今天的天气可能与过去的十天的气温有线性关系;或者有的人对食物有一种节俭的美德,他们做的饭菜能看出有些是上一顿的,当然也有一部分是今天的做的,再假设隔两顿的都被倒掉了,并且每天都是这样的,那么这碗饭菜可能就是一部分上一顿的再加上一部分今天现做的,这就是一个一阶的移动平均。
    ' r, e7 d7 S0 ?7 a4 K1 s! n- p
    6 A% p5 L! @% m; ?; k7 g移动平均法
    $ s) F4 S# w' O移动平均法是根据时间序列资料逐渐推移,依次计算包含一定项数的时序平均数, 以反映长期趋势的方法。当时间序列的数值由于受周期变动和不规则变动的影响,起伏 较大,不易显示出发展趋势时,可用移动平均法,消除这些因素的影响,分析、预测序列的长期趋势。  移动平均法有简单移动平均法,加权移动平均法,趋势移动平均法等。 5 T- D* }' N0 o7 }7 `; }
    + f% A: U; k  s7 U& z0 ~
    简单移动平均法
    4 `) o; |  i! ^. ]: f( S" T/ M
    1 F. t5 @, {7 O$ W: M0 D, l
    ; A! f. y, A6 l$ [1 D3 N% r3 ?3 p" @: H* H! T8 z/ _, `
    近N 期序列值的平均值作为未来各期的预测结果。一般 N 的取值范围: 5≤N≤ 200。当历史序列的基本趋势变化不大且序列中随机变动成分较多时,N 的 取值应较大一些。否则 N 的取值应小一些。在有确定的季节变动周期的资料中,移动平均的项数应取周期长度。选择佳 N 值的一个有效方法是,比较若干模型的预测误 差。预测标准误差小者为好。
    . C) Y; k, b1 ^( j* u; |! k6 R' @# `3 n& f; H" N0 T0 f
    简单移动平均法只适合做近期预测,而且是预测目标的发展趋势变化不大的情况。 如果目标的发展趋势存在其它的变化,采用简单移动平均法就会产生较大的预测偏差和滞后。 & [; C- R2 i$ y, O+ z6 G8 E, W, ~

    % y: e: Y& U1 \- Z$ X- Z例 1  某企业 1 月~11 月份的销售收入时间序列如表 1 示。试用一次简单滑动平 均法预测第 12 月份的销售收入。 / D. p  Q) E: y! X" v, u, m

    + \6 C) A7 |$ T8 ~3 _
    4 y# o( W) K% s4 p. R
    1 @2 J# |! X* I9 w0 _
    ! Q6 L$ a! w, o* g. f  n/ d$ _2 E
    计算的 Matlab 程序如下: " c* _$ f$ f- F
    , h- T5 b; H. A2 A! w% Y1 H+ d
    clc,clear
    # ]! H$ m' G1 y6 W$ J4 q: V: Ry=[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);   
    , h" c  \6 A2 J+ Y0 M  h( _ n=[4,5];   %n 为移动平均的项数 8 L7 V! q" y) ?8 E$ D. M
    for i=1:length(n)   
    , W: L% T. e( U9 v1 O7 K %由于 n 的取值不同,yhat 的长度不一致,下面使用了细胞数组    " C1 t3 P1 g0 Q+ i7 f
        for j=1:m-n(i)+1         
    ! h# z, @0 u, g  f        yhat{i}(j)=sum(y(j:j+n(i)-1))/n(i);     
    5 ~0 t  A4 b3 {    end   
    ( F" o9 y1 Q8 `# y! [. H    y12(i)=yhat{i}(end);     7 S% S6 s7 n# `& d% j
        s(i)=sqrt(mean((y(n(i)+1:m)-yhat{i}(1:end-1)).^2));
    * W4 c- }: r  G  n- I0 \+ Aend
    % {- z+ u  C$ c: Dy12,s
    2 j, h8 U7 f7 W% x: t- l8 [, j  ~# s% d* x8 d0 h
    加权移动平均法
      U  `) Y3 }* V- x在简单移动平均公式中,每期数据在求平均时的作用是等同的。但是,每期数据所包含的信息量不一样,近期数据包含着更多关于未来情况的信心。因此,把各期数据等同看待是不尽合理的,应考虑各期数据的重要性,对近期数据给予较大的权重,这就 是加权移动平均法的基本思想。 7 Z5 A& `4 D2 W+ O8 T
    / F$ `6 W8 x8 d1 x, G( S

    9 F% m( o* w* {% M$ D  ]
    8 E3 V. }% S7 K. @* z4 s1 ?/ u例 2  我国 1979~1988 年原煤产量如表 2 所示,试用加权移动平均法预测 1989 年 的产量
    ! O/ W- W" J8 g* m. O! \0 d6 z9 @
    3 m6 U3 ~) K! T4 [1 t' s, U
    % ]! _: v1 D4 ^$ Q8 \2 p! g
    7 v% Y/ P6 J; X8 j. w& E" E! [* N3 I/ T

    6 Q8 U4 c/ q% `/ u  k0 F
    6 Z. z0 t0 Y: z* g
    ( _" j! o9 G6 j% B2 _2 i' D& T8 N' A, y7 c( y
    计算的 MATLAB 程序如下: " w) @0 C5 d! o: B0 A: T# H

    ; y1 O9 K  J- l. p; A8 S9 P* `y=[6.35 6.20    6.22    6.66    7.15    7.89    8.72    8.94    9.28    9.8];
    ) u# k- c5 U: L/ e) @& Ww=[1/6;2/6;3/6]; 3 T9 y4 W6 \+ j( K2 d, H' P9 k  o9 `
    m=length(y);n=3; , V, n+ R& ?7 }5 v
    for i=1:m-n+1     $ |% {6 [+ R) |6 G$ E1 @2 ^, N
        yhat(i)=y(i:i+n-1)*w; : U) b8 s) e5 K$ s0 ?& t
    end # g% `) @6 D! ~  K8 K; u2 }' k! @0 j
    yhat
    / |. d) h5 _* S9 H+ ^4 g, jerr=abs(y(n+1:m)-yhat(1:end-1))./y(n+1:m)
    7 V0 ~9 z& j& R. C% qT_err=1-sum(yhat(1:end-1))/sum(y(n+1:m)) - D1 g) J' ?1 `( Y3 V$ x% }
    y1989=yhat(end)/(1-T_err)
    8 o8 ~, h9 n  t; s" U
      _, `, G; V; m, r 在加权移动平均法中,   的选择,同样具有一定的经验性。一般的原则是:近期 数据的权数大,远期数据的权数小。至于大到什么程度和小到什么程度,则需要按照预 测者对序列的了解和分析来确定。
    ; B9 t8 L( O3 S' y' d3 Q$ g, F$ b) U7 b/ Z6 Q3 W
    趋势移动平均法
    4 y' }2 t: ?" {' l: ~9 ^" W简单移动平均法和加权移动平均法,在时间序列没有明显的趋势变动时,能够准确 反映实际情况。但当时间序列出现直线增加或减少的变动趋势时,用简单移动平均法和 加权移动平均法来预测就会出现滞后偏差。因此,需要进行修正,修正的方法是作二次 移动平均,利用移动平均滞后偏差的规律来建立直线趋势的预测模型。这就是趋势移动平均法。  一次移动的平均数为 5 p5 j9 C! e$ y6 |; ]3 S- I

    * z, L1 }2 N; n/ r7 _2 s! R$ ]9 g9 Y  H+ n9 ^1 I- ^, T! y9 x8 H" z

    , D! A, S# W& B/ T: l! c; P' o( D$ O, t: J. g& w( M1 @

    3 u, |4 T( k9 I9 J/ E$ ^2 `9 L& x! V8 l# e

    4 J6 K' L* i1 L- S, L5 I! U; X例 3  我国 1965~1985 年的发电总量如表 3 所示,试预测 1986 年和 1987 年的发 电总量。
    7 E& u  Q) B  ?8 }- Q, c  w" |( H- @6 r6 r" S7 Z1 b* f2 ^

    / a" o& }# y, Z! x4 w" S! a  S" v1 h! C  q& V! Z4 V& G
    解  由散点图 1 可以看出,发电总量基本呈直线上升趋势,可用趋势移动平均法 来预测。 % {  Z& k  m7 j/ \
    & `) v2 b$ r8 W8 Q- l
      h1 c' n# y% V+ F8 O

    2 \% G1 k4 ?, ~2 O计算的 MATLAB 程序如下:
    4 s6 [8 w* H0 T* Y- s
    ! g& {& i, q4 ]" r8 {6 Qclc,clear
    6 X) x5 h9 _1 w: mload y.txt   %把原始数据保存在纯文本文件 y.txt 中
    4 g+ I  U, [' I) Fm1=length(y);   
    : a/ D6 f- z" q6 on=6;   %n 为移动平均的项数
      p% x" q; O- M2 U1 Z% Vfor i=1:m1-n+1     
    5 X; @: l+ k& @7 @5 P    yhat1(i)=sum(y(i:i+n-1))/n;
    - C5 B" x8 `; @7 H( G% N$ Kend
    3 r$ ?0 `. W+ Y8 dyhat1 $ f+ h' ?5 k6 ]
    m2=length(yhat1); 2 r. n: O( P$ e  r
    for i=1:m2-n+1   
    % |4 i  E/ Y  S. d* D9 D3 S    yhat2(i)=sum(yhat1(i:i+n-1))/n; 2 v" v3 O4 M2 Z9 f8 s
    end
    " j2 j. {" |/ E  ^, o1 gyhat2   ( Z9 G' u: X( j. v
    plot(1:21,y,'*')
    # p, O$ P5 C5 _; ta21=2*yhat1(end)-yhat2(end) : \5 r: K- _3 q: \8 A
    b21=2*(yhat1(end)-yhat2(end))/(n-1) 7 O+ q+ A1 g. h: ]6 M4 z
    y1986=a21+b21
    1 u; S3 [3 c) z0 h* ay1987=a21+2*b21 ) L  m5 `8 p' J) a$ Q4 Z/ A
    0 m$ ]; z0 |/ f/ p% [
    * Q! l$ r) c2 B  V" \

    ! I- V# d/ Y; X# @. F趋势移动平均法对于同时存在直线趋势与周期波动的序列,是一种既能反映趋势变 化,又可以有效地分离出来周期变动的方法。. R1 \% Z6 e* r6 C+ V) K- N
    9 L& v+ A- ?! o" v: Y" C+ E
    $ N$ [7 P# j9 n& h3 m' _
    5 u7 W0 d5 T5 f. S0 i2 S
    ————————————————
    * m" G6 H  O' a6 h: B# z+ Y6 m版权声明:本文为CSDN博主「wamg潇潇」的原创文章,遵循CC 4.0 BY-SA版权协议,转载请附上原文出处链接及本声明。: v9 F6 P1 l) c8 W" M, w
    原文链接:https://blog.csdn.net/qq_29831163/article/details/89440426
    3 L3 b3 d  ^* j: t. o' h! t) w0 a4 m: q* W  x: Y4 `9 Q  M; Q
    * I$ r3 C2 ^0 z6 I$ `  y" `
    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-6-15 07:25 , Processed in 7.215201 second(s), 51 queries .

    回顶部