QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 3418|回复: 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 |邮箱已经成功绑定
    移动平均法 可以作为一种数据平滑的方式 ,以每天的气温数据为例,今天的天气可能与过去的十天的气温有线性关系;或者有的人对食物有一种节俭的美德,他们做的饭菜能看出有些是上一顿的,当然也有一部分是今天的做的,再假设隔两顿的都被倒掉了,并且每天都是这样的,那么这碗饭菜可能就是一部分上一顿的再加上一部分今天现做的,这就是一个一阶的移动平均。0 x8 ^& ^7 C$ j9 z8 f
    + t5 Z& B0 Q- W4 j/ @- X+ z) u
    移动平均法) A' O. ]. d7 V" w+ l! B6 }
    移动平均法是根据时间序列资料逐渐推移,依次计算包含一定项数的时序平均数, 以反映长期趋势的方法。当时间序列的数值由于受周期变动和不规则变动的影响,起伏 较大,不易显示出发展趋势时,可用移动平均法,消除这些因素的影响,分析、预测序列的长期趋势。  移动平均法有简单移动平均法,加权移动平均法,趋势移动平均法等。
    ( b7 B! _' f+ |+ T5 p2 U8 h$ ^* b$ o3 @6 i! G) _: E
    简单移动平均法
    & ]9 i+ {5 P& u  h
    5 }& k9 a' j& v) o- d: F( O, n/ g3 N

    9 b6 W( r! G0 l) K近N 期序列值的平均值作为未来各期的预测结果。一般 N 的取值范围: 5≤N≤ 200。当历史序列的基本趋势变化不大且序列中随机变动成分较多时,N 的 取值应较大一些。否则 N 的取值应小一些。在有确定的季节变动周期的资料中,移动平均的项数应取周期长度。选择佳 N 值的一个有效方法是,比较若干模型的预测误 差。预测标准误差小者为好。 ! ~6 P5 x( H2 j
    & r: C3 p# i6 f5 k7 k0 X7 Q
    简单移动平均法只适合做近期预测,而且是预测目标的发展趋势变化不大的情况。 如果目标的发展趋势存在其它的变化,采用简单移动平均法就会产生较大的预测偏差和滞后。 " Q& g4 P7 ]8 [; y2 c# T8 x
    9 D- \7 L: f, Q
    例 1  某企业 1 月~11 月份的销售收入时间序列如表 1 示。试用一次简单滑动平 均法预测第 12 月份的销售收入。
    0 [% T7 ~. M7 C+ ?- a5 H4 o2 l) b& c+ s  _, I

    6 Y. n# B7 H4 @8 w# }: ^
    ; `% u( @: U$ p% j  h; y" t& R# g/ m9 o/ R" E& c

    / _) B- R6 E( P' Y计算的 Matlab 程序如下:
    , s. R' q1 ]6 R/ w6 _. j/ c+ g; A' K! T) N. ?5 D9 ~8 j
    clc,clear
    ) t9 m: ?1 {  z  i( xy=[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);   
    : h6 m3 i6 M, B- w n=[4,5];   %n 为移动平均的项数 ; \, J2 I' Q  _, @+ U7 O- A; ~" _+ w
    for i=1:length(n)    ) M+ f3 B& g9 Y0 c4 i5 T& O* k
    %由于 n 的取值不同,yhat 的长度不一致,下面使用了细胞数组    7 u6 |6 U9 j, L' S/ u* X: p
        for j=1:m-n(i)+1         
    - {" V) N( l  T/ [        yhat{i}(j)=sum(y(j:j+n(i)-1))/n(i);     
    ; f* h/ I+ A  y* O    end   
    / e) C& C8 I, e, u- Z. A  s! _  S7 t: v    y12(i)=yhat{i}(end);     ! I% q! c, `2 t3 ^$ I
        s(i)=sqrt(mean((y(n(i)+1:m)-yhat{i}(1:end-1)).^2));
      A' d$ @- N) @end 9 @6 ^( A) |  ]: w+ H
    y12,s
    9 H2 A  U% V2 k5 k( }) h- a5 D% q! C1 n6 \/ M$ P$ X0 t; c! ?5 T# |
    加权移动平均法
    " q  w, @# s" L/ J在简单移动平均公式中,每期数据在求平均时的作用是等同的。但是,每期数据所包含的信息量不一样,近期数据包含着更多关于未来情况的信心。因此,把各期数据等同看待是不尽合理的,应考虑各期数据的重要性,对近期数据给予较大的权重,这就 是加权移动平均法的基本思想。
    ' I  h. L# \+ I' t# J7 G4 B! `( m$ s( Q3 W

    * }6 L# [/ N- I" |* b0 [7 z/ Q+ r3 A% q% H" {" r
    例 2  我国 1979~1988 年原煤产量如表 2 所示,试用加权移动平均法预测 1989 年 的产量
    ' [3 i7 i. M1 i0 G; h6 N8 P6 H9 I
    ) y; _" U7 U$ }1 u# z% _
    ( Y: r$ S* }& F" [+ C' I, [5 W4 E5 v2 G! N8 k* W8 b* K# ?' m5 V
    # x% B' J6 A" @# M" o

    - H& ?/ h# N, V' F9 {: i; j2 g$ ]" s+ F* t% K; r. y

    9 C% l! D0 L0 Z! P& o7 f) N" p2 o
    5 l6 V1 B' T3 N% L1 ?计算的 MATLAB 程序如下: $ d( H1 w, Z8 R1 G( i: R! ]

    * {9 U5 j- z; x$ c( ]! Uy=[6.35 6.20    6.22    6.66    7.15    7.89    8.72    8.94    9.28    9.8]; * q4 x! D$ j* I2 e  k2 r  D
    w=[1/6;2/6;3/6];
    ! c2 {; C: t0 r" f2 ym=length(y);n=3;
    8 l. v3 C: `" n& i# p4 \for i=1:m-n+1     ! w( V6 K  n7 T0 v3 C  T+ t
        yhat(i)=y(i:i+n-1)*w; / ?' h" k" E! H2 b6 Q' R2 M0 N
    end % j$ F) [- G- ?& x/ j) ~' v
    yhat
    + E5 H, f' [/ T! s* E, perr=abs(y(n+1:m)-yhat(1:end-1))./y(n+1:m) 0 _. r6 e4 p, k" k2 V- J% V
    T_err=1-sum(yhat(1:end-1))/sum(y(n+1:m))
    ; s9 M! Y$ k: ?; N+ |  }$ Yy1989=yhat(end)/(1-T_err)
    3 K( u5 t. M- n  Z' j3 N+ n8 ?+ e6 u0 B
    在加权移动平均法中,   的选择,同样具有一定的经验性。一般的原则是:近期 数据的权数大,远期数据的权数小。至于大到什么程度和小到什么程度,则需要按照预 测者对序列的了解和分析来确定。 9 _3 j% H5 S7 ?) x# y* g4 I' E

    0 m- [5 B) W$ e; R趋势移动平均法
    * C# c+ H' s  [- U) I简单移动平均法和加权移动平均法,在时间序列没有明显的趋势变动时,能够准确 反映实际情况。但当时间序列出现直线增加或减少的变动趋势时,用简单移动平均法和 加权移动平均法来预测就会出现滞后偏差。因此,需要进行修正,修正的方法是作二次 移动平均,利用移动平均滞后偏差的规律来建立直线趋势的预测模型。这就是趋势移动平均法。  一次移动的平均数为 . U! k3 J' }: u
    - I1 H( J% j! |2 ]  l% h, b$ N

    : b! l, l7 N+ u' L) L: i8 t" S. k
    * o- s; G+ G: B
    5 N5 ]6 ]1 G6 w) ~% [
    # ~0 ^/ y: U, }9 {" Q6 v. r- G( J+ G
    3 ^/ V" R) B$ e% b6 ?6 f9 a$ V
    例 3  我国 1965~1985 年的发电总量如表 3 所示,试预测 1986 年和 1987 年的发 电总量。
    ( o1 T$ n8 b7 U& ~# Z* R8 A  T" k2 M, i$ j. J$ b, O
    1 j1 U3 T$ M: U5 l* I; p% o
    ! ^- a0 ~9 }/ r6 l& X3 U
    解  由散点图 1 可以看出,发电总量基本呈直线上升趋势,可用趋势移动平均法 来预测。 % c, k6 H$ z: P3 r$ U

    ! {  t# I. o- Y/ d4 j
    # @, B; c' J/ \: R" p; m5 l8 D8 r3 l+ m, l3 n
    计算的 MATLAB 程序如下: + w6 K7 Z% a# l: k) s2 U: J

    . q4 e8 N& @7 Z0 h- q1 Wclc,clear
    1 ?( t% H: _0 E, ?load y.txt   %把原始数据保存在纯文本文件 y.txt 中
    # i6 r$ \9 \5 f4 I) km1=length(y);    ! F; Q. J4 T$ g. n3 G% I, W
    n=6;   %n 为移动平均的项数 * j" r1 d# B3 ]9 {9 f$ f
    for i=1:m1-n+1     ( R7 C/ H. `2 m! z" h6 g2 t
        yhat1(i)=sum(y(i:i+n-1))/n;
    $ K5 _: Y+ |! ]  q0 |5 Z% Gend
    # O2 m, n( f  L9 [6 d- f" qyhat1 3 ^+ S% F( N! S# j
    m2=length(yhat1); ; S* }# H2 {4 [% \
    for i=1:m2-n+1   
    8 P3 @5 w6 c4 `9 b    yhat2(i)=sum(yhat1(i:i+n-1))/n; # m) Q% W' o' h& {/ ?
    end
    5 Q, q1 X2 g0 h. I+ b. s) ^yhat2   - V3 D2 W& S# f7 j. X
    plot(1:21,y,'*') 7 p( z% k  f" E1 B1 p# r, F: O
    a21=2*yhat1(end)-yhat2(end)
    7 z  W, Z$ b  @( `4 w% e! Pb21=2*(yhat1(end)-yhat2(end))/(n-1) & l& d( g( R; |; `* v" g- T$ W" R
    y1986=a21+b21 1 E4 F, c' z! r# M3 O
    y1987=a21+2*b21
    / u3 d, g0 `1 {( \# y, N* u. T$ F
      {& B7 F  u/ p
    ( ?, {5 K# O" A. ^# X4 r  i1 N  O4 k( U7 G- M) X1 u' P
    趋势移动平均法对于同时存在直线趋势与周期波动的序列,是一种既能反映趋势变 化,又可以有效地分离出来周期变动的方法。2 C' h5 e, o. W$ c2 w" S- S& z

    0 x8 M5 z2 A" k1 @! G, c# A6 K2 q. Q; O3 O, y

    6 A& L4 V5 L$ U8 }/ s————————————————
    ' t0 M0 ~+ R0 d# X" {, M3 B版权声明:本文为CSDN博主「wamg潇潇」的原创文章,遵循CC 4.0 BY-SA版权协议,转载请附上原文出处链接及本声明。
    5 v, d" Q: e/ [# i: X原文链接:https://blog.csdn.net/qq_29831163/article/details/894404264 e0 L: N6 C" c5 @% j, _& J

    . i- \0 q; `% h  t& {
      E* X% s- c) l  C8 N; q
    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-9-13 18:30 , Processed in 0.318029 second(s), 51 queries .

    回顶部