QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 3390|回复: 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 |邮箱已经成功绑定
    移动平均法 可以作为一种数据平滑的方式 ,以每天的气温数据为例,今天的天气可能与过去的十天的气温有线性关系;或者有的人对食物有一种节俭的美德,他们做的饭菜能看出有些是上一顿的,当然也有一部分是今天的做的,再假设隔两顿的都被倒掉了,并且每天都是这样的,那么这碗饭菜可能就是一部分上一顿的再加上一部分今天现做的,这就是一个一阶的移动平均。
    9 K$ ?* G0 [6 `8 D$ _. U
    ; A: ^7 V& `8 k& v移动平均法
    7 p0 `% x, s% `$ x5 K3 B移动平均法是根据时间序列资料逐渐推移,依次计算包含一定项数的时序平均数, 以反映长期趋势的方法。当时间序列的数值由于受周期变动和不规则变动的影响,起伏 较大,不易显示出发展趋势时,可用移动平均法,消除这些因素的影响,分析、预测序列的长期趋势。  移动平均法有简单移动平均法,加权移动平均法,趋势移动平均法等。
      t/ h) p  _6 n. v0 M2 [& C1 b; r) z# ]8 @& U% ?% W
    简单移动平均法
    # d* w8 h' L1 V- F
    + x5 U' v! z; b; X- P4 @* |; _* E
    3 U4 a, H4 I" u) P6 h
    1 H: E, j% i+ J- v& i; ^近N 期序列值的平均值作为未来各期的预测结果。一般 N 的取值范围: 5≤N≤ 200。当历史序列的基本趋势变化不大且序列中随机变动成分较多时,N 的 取值应较大一些。否则 N 的取值应小一些。在有确定的季节变动周期的资料中,移动平均的项数应取周期长度。选择佳 N 值的一个有效方法是,比较若干模型的预测误 差。预测标准误差小者为好。 6 N0 N2 \. K& y0 X2 z+ L
    1 y& f1 L; f. D+ j9 P
    简单移动平均法只适合做近期预测,而且是预测目标的发展趋势变化不大的情况。 如果目标的发展趋势存在其它的变化,采用简单移动平均法就会产生较大的预测偏差和滞后。 ' d5 t  R6 y" f% W. K0 a) A6 }
    4 {: d! q6 m9 U: y, b- C
    例 1  某企业 1 月~11 月份的销售收入时间序列如表 1 示。试用一次简单滑动平 均法预测第 12 月份的销售收入。 ; k" L1 J: u, O/ g0 J. K5 l! c/ J
    ' n$ G3 c4 k; Q9 _* f0 ?9 a

    - p1 n0 `9 i, T7 s6 j7 I
    7 n! E( |' h: j. k/ X2 r
    7 N4 v1 R, U' q+ q! Q3 I" C7 Z8 f3 e/ P" x
    计算的 Matlab 程序如下: + f! g' \6 B; A, Q

    # _; @/ t6 z' x3 V6 Y& |; Jclc,clear + ?9 w, J3 W2 e& X
    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);   
    & \7 E4 U; l/ O* s& {5 {0 H1 t n=[4,5];   %n 为移动平均的项数
    0 t& E/ m* M) m- r: Hfor i=1:length(n)   
    7 t. S- ]3 G/ z; _6 O %由于 n 的取值不同,yhat 的长度不一致,下面使用了细胞数组   
    2 ?* W' `/ l0 [# g& T: @+ I- O    for j=1:m-n(i)+1         
    0 l) u0 v& U  y4 O6 ~  r/ ~        yhat{i}(j)=sum(y(j:j+n(i)-1))/n(i);     2 `$ p8 L( f1 d6 V5 Z
        end    8 ^7 K% e6 B! [6 S
        y12(i)=yhat{i}(end);     # G7 J$ N% X7 G5 o; N1 N
        s(i)=sqrt(mean((y(n(i)+1:m)-yhat{i}(1:end-1)).^2));
    3 J! v. s! @: q* w; b4 I& {" o' ^& G# {% Zend . q! ?8 G3 m! j/ n' ?, N: A
    y12,s
    0 |, q* j! S  D% f* }# k
    0 K# M! l* A# S加权移动平均法
    ) t( c6 f! r. E在简单移动平均公式中,每期数据在求平均时的作用是等同的。但是,每期数据所包含的信息量不一样,近期数据包含着更多关于未来情况的信心。因此,把各期数据等同看待是不尽合理的,应考虑各期数据的重要性,对近期数据给予较大的权重,这就 是加权移动平均法的基本思想。 2 b$ ?! r5 H( p% h. q, @

    ( y4 k- o- h4 g$ t9 @& k: r! x) S6 @+ j9 j4 t* [1 U- X

    2 m; r& K& h+ `( S5 `& W. y" J例 2  我国 1979~1988 年原煤产量如表 2 所示,试用加权移动平均法预测 1989 年 的产量9 g' r6 ^+ U( L$ Q* W1 ?

    - d; _: F* q4 R3 K7 G1 Y/ x/ r' i; I8 @; U/ Y; T. w* t& U

    " }5 _5 r6 Y3 h: d6 i
    ; u1 N' G: j. B, o8 _" b& @! D5 d" {6 q+ d

    . ?+ x% P' K0 V& A- v0 d; a4 @( y) ?3 ^( g2 H8 A- E4 Y/ G

    5 N) h) ?" Z: h! l( K9 C$ m7 `计算的 MATLAB 程序如下:
    - o7 E5 q9 |7 n* Q8 A7 g
    & N8 l/ ]# |! M7 }/ W+ wy=[6.35 6.20    6.22    6.66    7.15    7.89    8.72    8.94    9.28    9.8];   M( H4 L3 z, X$ ^, ^4 c5 p
    w=[1/6;2/6;3/6];
    # f* _" h& ]9 S( _+ Z3 Fm=length(y);n=3; ( m: z  [* n' Z7 S; s
    for i=1:m-n+1     
    $ A* b! q$ r5 J5 S2 J' J- |2 u    yhat(i)=y(i:i+n-1)*w;
    5 ~: z/ x, G, D5 eend % z6 p' f) y. r& I
    yhat 7 F* j% B% t4 ]* ?9 O1 M- R8 w$ T
    err=abs(y(n+1:m)-yhat(1:end-1))./y(n+1:m)
    . u2 h7 t5 [1 |& X+ H3 |T_err=1-sum(yhat(1:end-1))/sum(y(n+1:m))
    ( E% P; ]/ O+ Q8 \& k" By1989=yhat(end)/(1-T_err)" n4 G  m2 p- ~+ I6 k

    # J% H4 N  T7 K, }5 k5 Z' r+ ~: z 在加权移动平均法中,   的选择,同样具有一定的经验性。一般的原则是:近期 数据的权数大,远期数据的权数小。至于大到什么程度和小到什么程度,则需要按照预 测者对序列的了解和分析来确定。 $ F4 j5 K/ L- N, ?& D
    - ]% F, ]- N4 ~3 a0 E- `" m( ?/ ?
    趋势移动平均法8 ?# g0 Z" S! {1 X6 T/ {1 b
    简单移动平均法和加权移动平均法,在时间序列没有明显的趋势变动时,能够准确 反映实际情况。但当时间序列出现直线增加或减少的变动趋势时,用简单移动平均法和 加权移动平均法来预测就会出现滞后偏差。因此,需要进行修正,修正的方法是作二次 移动平均,利用移动平均滞后偏差的规律来建立直线趋势的预测模型。这就是趋势移动平均法。  一次移动的平均数为
    # k  X$ |$ ~+ Y/ _/ s: U$ b* c' v% z3 l; M
    8 c! ]. X% X1 e& M
    3 m/ s# e, ~6 Z- E# t1 E* d& U
    6 S# I. F: A  T$ Z6 h$ h

    . ]5 c, P* S- w' ~' V1 ]8 q
    4 a% R1 U/ g1 e3 E: Y3 n
    ( w/ n3 ~6 N1 {$ V7 _例 3  我国 1965~1985 年的发电总量如表 3 所示,试预测 1986 年和 1987 年的发 电总量。 4 N2 A# O: y( B3 G$ W. q- L$ z
    + C0 S+ u: v% {8 A" `; m% T7 }; ~9 D

    " [2 L* _& o, k
    ! k& `2 N% T' Z5 X9 m6 g9 M, I& h解  由散点图 1 可以看出,发电总量基本呈直线上升趋势,可用趋势移动平均法 来预测。 7 }4 B+ a5 y6 U7 A7 J
    & g% d) C$ [! z; _) q! \/ e1 K" a0 k

    " n, o' ], @' W  J+ c4 j$ W, J' A! d( b
    计算的 MATLAB 程序如下: 8 K( x) K( z' k1 v

    9 |, @; V6 O4 }7 l' @- hclc,clear 6 s/ o( ^, {& Q  J! e: K
    load y.txt   %把原始数据保存在纯文本文件 y.txt 中
    . i( S. {' L8 M( Mm1=length(y);   
    8 s, {- {" U% u/ @8 _9 i; r7 m3 y6 Jn=6;   %n 为移动平均的项数 + O6 H+ G$ A; F( {! Q% B
    for i=1:m1-n+1     6 J, `% K# k+ G* A
        yhat1(i)=sum(y(i:i+n-1))/n; 0 H7 ^- g# \# s; g; h. V! W$ Z
    end
    * F. m' Q5 D& ]; Wyhat1
    4 n! R( h+ S# R" d) H' R+ f7 L! lm2=length(yhat1); - F) S2 b! \8 Q/ R4 B, g  o# K; F- s: M
    for i=1:m2-n+1   
    & C! d! D1 x, p2 x  W    yhat2(i)=sum(yhat1(i:i+n-1))/n;
    4 Q% p- o$ G" j0 x, q  Wend
    * ?  l: g% ^" @1 M* B2 z; K; Ayhat2   2 N9 l* O, P% h' {
    plot(1:21,y,'*')
    2 @$ i+ V6 H( L5 O$ q; {. K, }a21=2*yhat1(end)-yhat2(end) ; F2 K, I; v- `. X5 m7 [7 ?
    b21=2*(yhat1(end)-yhat2(end))/(n-1) 0 B+ z# T9 b$ K+ {( x6 E
    y1986=a21+b21 - O: ?9 s  q" ^1 z
    y1987=a21+2*b21 ! @  l) ]7 O& g2 V7 _

    & ~. ^5 c8 ~/ U: S$ u
    " r. V  Z4 R# w1 }; p! L. z6 P
    0 G  E- p( P" l+ D9 e趋势移动平均法对于同时存在直线趋势与周期波动的序列,是一种既能反映趋势变 化,又可以有效地分离出来周期变动的方法。, i% J' ?* U, u; n3 w$ A
    ' D, i! S, A# y$ r, M7 v6 M: {
    # p, |4 p- n/ ]) g; W

    ' X4 t) [4 x0 Z1 g" A! Y————————————————
    6 P* Y8 d; V0 u5 v, E- L+ Y版权声明:本文为CSDN博主「wamg潇潇」的原创文章,遵循CC 4.0 BY-SA版权协议,转载请附上原文出处链接及本声明。
    # r/ t4 h: s  A% {, D. x1 Z原文链接:https://blog.csdn.net/qq_29831163/article/details/89440426. g: w4 l" j( L7 G( p3 l
    - d7 v' A0 E/ k! _! R2 h0 d2 z
    ' y* V  @7 d& E6 S% ?- _! |8 a
    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-8-4 11:57 , Processed in 0.485736 second(s), 50 queries .

    回顶部