QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 3316|回复: 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 |邮箱已经成功绑定
    移动平均法 可以作为一种数据平滑的方式 ,以每天的气温数据为例,今天的天气可能与过去的十天的气温有线性关系;或者有的人对食物有一种节俭的美德,他们做的饭菜能看出有些是上一顿的,当然也有一部分是今天的做的,再假设隔两顿的都被倒掉了,并且每天都是这样的,那么这碗饭菜可能就是一部分上一顿的再加上一部分今天现做的,这就是一个一阶的移动平均。
    + t3 n( P/ f. i- k- [
      U! n$ I- M, X$ Y移动平均法+ G! n1 P: Z5 S3 T3 n6 _+ ~) r
    移动平均法是根据时间序列资料逐渐推移,依次计算包含一定项数的时序平均数, 以反映长期趋势的方法。当时间序列的数值由于受周期变动和不规则变动的影响,起伏 较大,不易显示出发展趋势时,可用移动平均法,消除这些因素的影响,分析、预测序列的长期趋势。  移动平均法有简单移动平均法,加权移动平均法,趋势移动平均法等。
      m8 m" V, p1 u7 k/ Y+ c% C2 j3 y6 M5 w7 ~( A0 Q( m# v, N1 l9 Y
    简单移动平均法  z. d" h+ e7 t1 D

    2 {( k+ H/ H8 ?) |7 R
    % G$ d( Q. q5 @: a
    9 M4 }; w/ a" o. L# w3 `( a) P0 ~近N 期序列值的平均值作为未来各期的预测结果。一般 N 的取值范围: 5≤N≤ 200。当历史序列的基本趋势变化不大且序列中随机变动成分较多时,N 的 取值应较大一些。否则 N 的取值应小一些。在有确定的季节变动周期的资料中,移动平均的项数应取周期长度。选择佳 N 值的一个有效方法是,比较若干模型的预测误 差。预测标准误差小者为好。 0 B& l; c" F& n/ b3 l) y7 a
    : B" b0 T0 c: o
    简单移动平均法只适合做近期预测,而且是预测目标的发展趋势变化不大的情况。 如果目标的发展趋势存在其它的变化,采用简单移动平均法就会产生较大的预测偏差和滞后。 1 Z/ J( X! e  h1 B& g4 m2 s

    1 C. m' L" o. I$ o/ O2 @4 r9 B. j例 1  某企业 1 月~11 月份的销售收入时间序列如表 1 示。试用一次简单滑动平 均法预测第 12 月份的销售收入。
    1 f' G0 v. S& Q' ]' [
    ( B* P% F4 b' g; [& R) V4 K9 w7 V5 y/ s+ t- A: S# ]& x
    + _0 d' O- {/ U
    5 \" X. k4 ]" d: L, }% ~8 J3 p: V

    * ^1 G" e) [! T计算的 Matlab 程序如下:
    + j/ @1 X. T/ d7 g7 `
    7 m+ k, @) W; R( ~4 X* H: V, g; ]! sclc,clear   @& W6 T5 j' d) T
    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);   
    , o: n2 y" S: f: X- s0 ^5 i6 I* b n=[4,5];   %n 为移动平均的项数 - l1 h, y1 P' I  G3 I
    for i=1:length(n)   
    " q* Z& m* h* e2 _6 r# p1 Z %由于 n 的取值不同,yhat 的长度不一致,下面使用了细胞数组    ; Z5 |8 ?/ K6 K+ G- Z, x! G$ x
        for j=1:m-n(i)+1         
    . }4 ^# y, X+ {' o. `( h- `        yhat{i}(j)=sum(y(j:j+n(i)-1))/n(i);     / N" B) i# e9 j9 t/ E
        end   
    ( A+ K* [( n' Z! F7 f    y12(i)=yhat{i}(end);     6 t) L3 g% p5 k4 Z4 j( \. c8 f
        s(i)=sqrt(mean((y(n(i)+1:m)-yhat{i}(1:end-1)).^2));
    , o6 ?' F4 y' C( C2 D- N# r' qend ! w" _4 X3 L3 h6 n4 h) o
    y12,s ' |$ K1 o' g; F8 V

    9 v  n/ m( {% ~0 n3 d5 ?+ t4 l; c加权移动平均法3 k+ g' K8 w. X+ e
    在简单移动平均公式中,每期数据在求平均时的作用是等同的。但是,每期数据所包含的信息量不一样,近期数据包含着更多关于未来情况的信心。因此,把各期数据等同看待是不尽合理的,应考虑各期数据的重要性,对近期数据给予较大的权重,这就 是加权移动平均法的基本思想。
    ' G/ |4 x1 U- V: E' z3 ~9 _6 Q: P; y. c4 Z9 i
    & I4 N  K# R8 ~* S" M$ n. i" y& F
    1 ?3 E% @* V! d# ~4 H8 W) X
    例 2  我国 1979~1988 年原煤产量如表 2 所示,试用加权移动平均法预测 1989 年 的产量
    : H0 N9 I( _: Y/ a* N$ H" q& u6 I$ M7 u5 [% `1 n& Y

    : ?3 n! J% p6 [: b  j# a+ y$ F: U: D

    6 T5 p% k6 V7 V" @  b3 |
    / E, ?& @% U- T+ {1 [) P
    . F3 p9 v: p0 ?" f+ }* D/ E" q3 r# E/ y- Q
    . _+ L( o/ f! F
    计算的 MATLAB 程序如下:
    7 A1 R2 l4 l( l% ]. L$ x8 T9 S4 [- h
    y=[6.35 6.20    6.22    6.66    7.15    7.89    8.72    8.94    9.28    9.8]; / ^+ }2 A; L; k7 p7 d- Q
    w=[1/6;2/6;3/6];
    7 \- t( |8 ^: Y- ?m=length(y);n=3;
      D6 Q  c- r( o  ?for i=1:m-n+1     
    / t7 b7 m  O. X. ^5 S8 j, `    yhat(i)=y(i:i+n-1)*w; # Y+ t# a. D9 Q% K
    end ) \: Q  |5 B- t( [
    yhat
    + Q2 l3 y; q9 O6 Z9 kerr=abs(y(n+1:m)-yhat(1:end-1))./y(n+1:m)
    + M. ]' y# A9 U+ g% ]1 PT_err=1-sum(yhat(1:end-1))/sum(y(n+1:m)) / w: q6 ~$ q- w$ t" v( p
    y1989=yhat(end)/(1-T_err)
    1 H  w1 S; |& e* O8 W  t$ f& \! p
    ( k  Y4 o8 |  _% m$ f 在加权移动平均法中,   的选择,同样具有一定的经验性。一般的原则是:近期 数据的权数大,远期数据的权数小。至于大到什么程度和小到什么程度,则需要按照预 测者对序列的了解和分析来确定。
    0 V$ X% F7 K4 f6 S+ m8 o7 B& c3 k
    趋势移动平均法( P1 c+ `4 j, F* n- o9 M
    简单移动平均法和加权移动平均法,在时间序列没有明显的趋势变动时,能够准确 反映实际情况。但当时间序列出现直线增加或减少的变动趋势时,用简单移动平均法和 加权移动平均法来预测就会出现滞后偏差。因此,需要进行修正,修正的方法是作二次 移动平均,利用移动平均滞后偏差的规律来建立直线趋势的预测模型。这就是趋势移动平均法。  一次移动的平均数为 * E% \  n$ Y+ Q) w
    6 @, i# c( M: h" c5 V# m- O: k
    : j* u8 N6 z( u; M! Y5 b, Z( i

    8 c5 k9 o& m1 e' x. b/ A+ G) m
    0 z( R7 N" E" Q+ P  ]0 P+ E- ]6 {* X2 Z+ u6 Q) C
      x' e- P- G- W/ s3 O

    9 M0 t) B- Q. W: w! c# q9 O( j) U例 3  我国 1965~1985 年的发电总量如表 3 所示,试预测 1986 年和 1987 年的发 电总量。
    ; n% [; W" y2 R/ M
    " S! y* l6 J5 m6 ~8 Q% D8 ]2 o; P. R, Z
    7 d+ b3 J# t+ x' X
    解  由散点图 1 可以看出,发电总量基本呈直线上升趋势,可用趋势移动平均法 来预测。 , H+ l& _7 q$ J$ q3 b: e4 M% G
    2 `/ o4 }, v# e+ I# D1 V; [! |
    # ^2 J7 e. P0 J

    7 B5 V* c' [- [% p0 c/ F* S计算的 MATLAB 程序如下:
    # [$ M2 t; E! s2 W" y* P
    8 [# |6 h$ R+ L1 [clc,clear 9 ?- P, B! X# e9 K
    load y.txt   %把原始数据保存在纯文本文件 y.txt 中
    ( l0 J( O  v1 ~/ |8 B. N1 Um1=length(y);    : J# F! \8 M; q3 h" Z1 H
    n=6;   %n 为移动平均的项数 3 v% Z* P# B8 O
    for i=1:m1-n+1     - F: b6 c$ [9 F1 V% M5 O( Q
        yhat1(i)=sum(y(i:i+n-1))/n; & Q9 S" E' J5 y$ z5 p
    end & L  X1 i4 f# Z% D
    yhat1   `& O4 O: e) O+ C- m
    m2=length(yhat1); # e2 s3 L+ V5 n8 n9 F4 T
    for i=1:m2-n+1    3 j  Z' }0 u, _6 F+ s
        yhat2(i)=sum(yhat1(i:i+n-1))/n;
    0 j, J& [% i; S4 Hend % S6 V0 M9 [8 ]7 f, x9 d, S% A
    yhat2   
    1 \3 L5 j+ d' G1 E2 c& m) gplot(1:21,y,'*') ! v2 J* y2 }" P( K1 _& [5 @% _& a
    a21=2*yhat1(end)-yhat2(end)
    : z1 G- i$ B! {b21=2*(yhat1(end)-yhat2(end))/(n-1)
    . d4 r/ I- m# E. ]2 Ly1986=a21+b21
    3 l3 N  M0 p, g7 @* u& ny1987=a21+2*b21
    ( u% N" g. |2 O, j; w) K( P1 n9 F: Y

    % E6 C; k( ], x3 f  P3 d! L$ V+ C+ Q
    1 O  u8 p0 C6 J6 ~趋势移动平均法对于同时存在直线趋势与周期波动的序列,是一种既能反映趋势变 化,又可以有效地分离出来周期变动的方法。
      ~) U: t6 i7 R1 u/ ^4 z
    & Z! j  @: Q  d7 w9 E' b. @/ V* L! J5 r6 ]9 y! p. I6 j9 @

    # p, m% i$ X( z# x* ]5 N* P* a————————————————! g+ W! S3 B! C% ~/ |! i) p: e& p7 a
    版权声明:本文为CSDN博主「wamg潇潇」的原创文章,遵循CC 4.0 BY-SA版权协议,转载请附上原文出处链接及本声明。
    5 ^$ F, l( M1 f1 h! y4 i原文链接:https://blog.csdn.net/qq_29831163/article/details/89440426
    2 Y, A- ^5 Y  `+ V; z7 `% h8 q5 W1 ~+ [3 W: R- }1 a# ~

    " T  G9 ?! r& j; }' ~  @0 A: n% U
    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-17 07:38 , Processed in 0.304879 second(s), 50 queries .

    回顶部