QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 3383|回复: 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 m: O* A0 L+ j3 }6 _

    , n  I# T7 l9 x" {  {移动平均法7 g; Z. R2 S3 L7 q$ N5 x/ S
    移动平均法是根据时间序列资料逐渐推移,依次计算包含一定项数的时序平均数, 以反映长期趋势的方法。当时间序列的数值由于受周期变动和不规则变动的影响,起伏 较大,不易显示出发展趋势时,可用移动平均法,消除这些因素的影响,分析、预测序列的长期趋势。  移动平均法有简单移动平均法,加权移动平均法,趋势移动平均法等。 - p9 @# x0 g: B6 z! f; ^

    . l* T* h* W% N& B0 m. e简单移动平均法
    : d& x5 U' p+ `
    / X7 {9 E! M# L: \; O; H! Z0 C, y" @4 `( |. B; B: \5 t
    0 B$ S. ?- v; q) O5 j, X* V
    近N 期序列值的平均值作为未来各期的预测结果。一般 N 的取值范围: 5≤N≤ 200。当历史序列的基本趋势变化不大且序列中随机变动成分较多时,N 的 取值应较大一些。否则 N 的取值应小一些。在有确定的季节变动周期的资料中,移动平均的项数应取周期长度。选择佳 N 值的一个有效方法是,比较若干模型的预测误 差。预测标准误差小者为好。
    ! D" q6 u( ^# ]* _+ N* D" ~# [* j4 t! o
    简单移动平均法只适合做近期预测,而且是预测目标的发展趋势变化不大的情况。 如果目标的发展趋势存在其它的变化,采用简单移动平均法就会产生较大的预测偏差和滞后。 2 p+ [9 b) C1 y$ h: q

    4 u7 R5 z9 N: L2 I例 1  某企业 1 月~11 月份的销售收入时间序列如表 1 示。试用一次简单滑动平 均法预测第 12 月份的销售收入。
    : {3 Z$ T; I) Y+ ~9 {9 g# e1 W4 ?
    % P) c( u+ Z+ u6 [0 i" w7 {' f; |5 Y7 b# u
    & a/ O1 a; q5 w/ ]  W
    : z+ K' L  z( f/ T! ]" b& ]

    7 J( b- N$ h# h% P' x( n6 J( g# c计算的 Matlab 程序如下:
    4 O1 b% M6 I  n' u' o% Z: e" S! h& J( A. U! n/ H; h; \
    clc,clear . I+ ?$ q8 g6 m6 i
    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);   
    5 |6 |( V/ M; h4 F; ?0 k n=[4,5];   %n 为移动平均的项数
    * r( i3 z! {6 s5 W( G& }1 kfor i=1:length(n)   
    - y! r. O5 f4 O. T% ]+ { %由于 n 的取值不同,yhat 的长度不一致,下面使用了细胞数组   
    4 p6 Q; w9 [1 a5 S$ L    for j=1:m-n(i)+1         2 P5 i5 _* t7 G! V
            yhat{i}(j)=sum(y(j:j+n(i)-1))/n(i);     
      `$ s! H$ F& n" v- u, d    end   
    % \& C# ^# ~- l' _+ @- u* j    y12(i)=yhat{i}(end);     5 q7 s, j! F6 z9 q) f
        s(i)=sqrt(mean((y(n(i)+1:m)-yhat{i}(1:end-1)).^2));
    : ~6 A! d8 P2 s+ |3 c- d" wend
    5 {0 e8 Q8 h9 qy12,s
    / }; n9 t, B+ v* [5 U' _# u) `3 u8 e7 C$ s4 B  ~
    加权移动平均法) }# S$ N2 M7 ]/ u5 M2 t# q0 m* X
    在简单移动平均公式中,每期数据在求平均时的作用是等同的。但是,每期数据所包含的信息量不一样,近期数据包含着更多关于未来情况的信心。因此,把各期数据等同看待是不尽合理的,应考虑各期数据的重要性,对近期数据给予较大的权重,这就 是加权移动平均法的基本思想。 ; B: n1 |" @8 S. ?& R
    ' [  ~, n# H% U2 k% p

    # F- j* d/ h% d3 H* w. Q# l. e! i: r: A2 [7 x4 E
    例 2  我国 1979~1988 年原煤产量如表 2 所示,试用加权移动平均法预测 1989 年 的产量5 H1 k! r: ~: p; Y) I

    8 l6 O5 Z. H2 K: u- o8 _3 d7 [
    8 v  Q# f% q$ A+ A
    ; _# F+ U4 Y' \1 t6 f
    ' p5 g+ l. Q& V# f/ c# s' k
    & `9 T% X1 c* U( V4 E
    9 N7 ?2 I' u) T- O- H+ p0 J# f
    ' n0 ?/ Q3 @0 N1 @2 q* B+ G0 c2 g! X6 C
    计算的 MATLAB 程序如下: ! h& W- S2 u& A& M

    0 [5 S# W9 V2 Q7 Y0 R: _y=[6.35 6.20    6.22    6.66    7.15    7.89    8.72    8.94    9.28    9.8]; 7 Q! L6 \6 s9 V$ R+ h* |4 [
    w=[1/6;2/6;3/6];
    7 I; E0 F/ M" Jm=length(y);n=3; " x! N/ @6 _/ L; g% w, t
    for i=1:m-n+1     
    . V3 T5 ]( ], @    yhat(i)=y(i:i+n-1)*w; 6 w( {" Y/ U" z+ d1 ]* d* {* N: e' ^6 ]
    end
    * ]. {0 v: a+ X5 T2 E! ~$ B) ayhat
    ' m' n8 B+ W' K9 `0 Z1 S7 Yerr=abs(y(n+1:m)-yhat(1:end-1))./y(n+1:m)
    5 N8 [* q! w' A7 J- qT_err=1-sum(yhat(1:end-1))/sum(y(n+1:m))
    8 E2 s; h6 F6 q9 Y- ^0 iy1989=yhat(end)/(1-T_err)
    4 F  j% F+ N: R& G: B2 P: c- t, ]4 r
    在加权移动平均法中,   的选择,同样具有一定的经验性。一般的原则是:近期 数据的权数大,远期数据的权数小。至于大到什么程度和小到什么程度,则需要按照预 测者对序列的了解和分析来确定。
    % g2 E. ^% Q, ~7 g0 O0 G* A
      A& i1 N8 x1 q趋势移动平均法
    ; b/ J4 y1 G: r简单移动平均法和加权移动平均法,在时间序列没有明显的趋势变动时,能够准确 反映实际情况。但当时间序列出现直线增加或减少的变动趋势时,用简单移动平均法和 加权移动平均法来预测就会出现滞后偏差。因此,需要进行修正,修正的方法是作二次 移动平均,利用移动平均滞后偏差的规律来建立直线趋势的预测模型。这就是趋势移动平均法。  一次移动的平均数为 " L( e4 D7 J# R4 m! Z! H! }
    9 m" c/ Y5 H0 S
    : d, U/ Y8 b- Z, K$ }. w; u4 J+ |

    ) ]" o" r, z$ v1 D) t
    : |8 S1 Q  M# i5 V2 T# |! ^9 n5 P, Q# H2 |

    8 T- N$ B% `: H4 q
    ' w% t8 ]$ F$ N$ `2 Y例 3  我国 1965~1985 年的发电总量如表 3 所示,试预测 1986 年和 1987 年的发 电总量。
    9 J2 D! C/ e3 q$ j6 {- n" w. X
    & ^' X$ O2 o1 N1 ~; o9 S
    $ R# j+ m# o/ j6 I5 e
    $ @0 G1 G* h% L+ ^. S% u* E解  由散点图 1 可以看出,发电总量基本呈直线上升趋势,可用趋势移动平均法 来预测。
    ! Y' a9 j. O5 L: q2 |) K& D2 b4 B, W8 n/ i+ `& y* q  Q6 v: _
    4 c) X5 J# O- k$ f  X( O
    5 B2 g' u- B+ z+ ^; b. I) l
    计算的 MATLAB 程序如下:
    & t* @5 _0 B  O$ L  r# O# D9 T1 `8 V1 W9 g
    clc,clear - a9 p9 {5 Q9 S6 u6 r  n1 G
    load y.txt   %把原始数据保存在纯文本文件 y.txt 中
    9 G4 j2 C$ I7 ?/ I/ rm1=length(y);   
    4 C& S7 t. I! V% I* l8 q5 ~n=6;   %n 为移动平均的项数
    . L0 Z" H  i" B- efor i=1:m1-n+1     # l) B8 g  V+ ]
        yhat1(i)=sum(y(i:i+n-1))/n;
    2 d) O; O- C- B; S+ l) v/ Mend # B) y( U( M5 c8 }* f
    yhat1 ; o* _. Y* R6 r& \
    m2=length(yhat1);
    & r/ A5 Q' @4 L% d8 K, {for i=1:m2-n+1    0 T! ?) \3 _* r, @
        yhat2(i)=sum(yhat1(i:i+n-1))/n;
    # {" f: S: _: vend 3 [% q: |) A& F1 N/ u: r  _
    yhat2   ' }; W" ~; E: _
    plot(1:21,y,'*')
    . W: T; Z1 I: V, e9 Ea21=2*yhat1(end)-yhat2(end)
    9 Y: F; e) o; J9 Eb21=2*(yhat1(end)-yhat2(end))/(n-1)
    / O+ U' ?$ n/ r  g$ i" Xy1986=a21+b21
    6 A4 s6 f0 p( k1 N3 [& u& Oy1987=a21+2*b21 * i+ d& C+ }) d6 o! c% s$ Y# b/ u
    $ ^1 E- I$ O' t5 m
    9 M5 l$ k' a" y3 @3 X$ ^/ j7 a. _
    ' q- ^% C* K' `& _/ G" p
    趋势移动平均法对于同时存在直线趋势与周期波动的序列,是一种既能反映趋势变 化,又可以有效地分离出来周期变动的方法。
    4 {' k. O1 J8 M) D3 G1 c0 s0 A! B% G) c% D
    # ]  }+ w  Q; I/ J
    * l; x$ _5 d3 o- S: S# S+ i
    ————————————————
    , Z. e! z; V  g5 W8 v% |7 B- v. |7 ?版权声明:本文为CSDN博主「wamg潇潇」的原创文章,遵循CC 4.0 BY-SA版权协议,转载请附上原文出处链接及本声明。
    ; W  [5 q& @$ ~$ V5 g0 m. U4 [原文链接:https://blog.csdn.net/qq_29831163/article/details/89440426
    2 [' `# b5 d2 D! p9 b+ R3 h3 W9 N* P+ D% s

    * E4 y3 E; S5 K' @
    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-7-29 01:28 , Processed in 0.450791 second(s), 51 queries .

    回顶部