QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 3310|回复: 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 |邮箱已经成功绑定
    移动平均法 可以作为一种数据平滑的方式 ,以每天的气温数据为例,今天的天气可能与过去的十天的气温有线性关系;或者有的人对食物有一种节俭的美德,他们做的饭菜能看出有些是上一顿的,当然也有一部分是今天的做的,再假设隔两顿的都被倒掉了,并且每天都是这样的,那么这碗饭菜可能就是一部分上一顿的再加上一部分今天现做的,这就是一个一阶的移动平均。; `. n5 u4 ?- v2 C% e* |* I

    1 Q  H: P8 E6 P移动平均法2 F8 v/ c9 L% Y* h' N
    移动平均法是根据时间序列资料逐渐推移,依次计算包含一定项数的时序平均数, 以反映长期趋势的方法。当时间序列的数值由于受周期变动和不规则变动的影响,起伏 较大,不易显示出发展趋势时,可用移动平均法,消除这些因素的影响,分析、预测序列的长期趋势。  移动平均法有简单移动平均法,加权移动平均法,趋势移动平均法等。   R7 n0 x8 |+ y/ R8 F/ \
    9 U& G- ^  A# h$ R( I: g
    简单移动平均法9 a% E! _% r- v. S" s: Z( S6 ^
    * K# E7 \% E, w7 ^

    : l0 A" `7 Y7 _* A& C9 g* Q0 R3 y+ P, U# [, A
    近N 期序列值的平均值作为未来各期的预测结果。一般 N 的取值范围: 5≤N≤ 200。当历史序列的基本趋势变化不大且序列中随机变动成分较多时,N 的 取值应较大一些。否则 N 的取值应小一些。在有确定的季节变动周期的资料中,移动平均的项数应取周期长度。选择佳 N 值的一个有效方法是,比较若干模型的预测误 差。预测标准误差小者为好。 ( Y. X9 Q( A$ J- C8 c) j
    6 d2 K! k( s1 `- N
    简单移动平均法只适合做近期预测,而且是预测目标的发展趋势变化不大的情况。 如果目标的发展趋势存在其它的变化,采用简单移动平均法就会产生较大的预测偏差和滞后。 . k: P5 k! B& n2 \

    8 X  S# r) ^: J: k1 [例 1  某企业 1 月~11 月份的销售收入时间序列如表 1 示。试用一次简单滑动平 均法预测第 12 月份的销售收入。
      \+ n5 x9 y% d; H$ j
    / ]! C$ c3 e& T' v+ K1 W1 y4 ~
    0 j; t$ N: F+ L& L3 k" b5 N/ |( u: L. D

    3 P; y* O1 B6 D( p! Z
    + v' k& ?" N1 a7 [( r5 h计算的 Matlab 程序如下:
    9 f2 s* [! }( S6 t- }% @: V
    5 z: X- O% Z/ D, V) L+ d! _clc,clear # O1 M4 z; k: C( ]1 N  P( c' m2 ^/ F8 A
    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);   0 G0 i$ W7 S/ \" G+ t
    n=[4,5];   %n 为移动平均的项数
    - A* @2 J. x9 v8 o  y( z' _for i=1:length(n)    $ Y0 h. M6 U, r/ I
    %由于 n 的取值不同,yhat 的长度不一致,下面使用了细胞数组   
    . j' B; M2 N8 k% v    for j=1:m-n(i)+1         
    * y$ b4 t4 \7 L6 e  o/ ]/ b% c        yhat{i}(j)=sum(y(j:j+n(i)-1))/n(i);     
    * d8 R7 U6 J. r# y    end   
    / c7 k& Z1 O" p, w( m4 v8 f2 ]    y12(i)=yhat{i}(end);     
    * U  l: H3 C; S' P# p8 X; Q    s(i)=sqrt(mean((y(n(i)+1:m)-yhat{i}(1:end-1)).^2)); & n/ B( v8 P' D" S2 B
    end
    + Y$ `4 Z3 V. O. S" q. D; Zy12,s 1 ~7 \! E+ f, m
    3 k5 @" Y# M. {
    加权移动平均法" v( C5 j! Y" v; c+ x
    在简单移动平均公式中,每期数据在求平均时的作用是等同的。但是,每期数据所包含的信息量不一样,近期数据包含着更多关于未来情况的信心。因此,把各期数据等同看待是不尽合理的,应考虑各期数据的重要性,对近期数据给予较大的权重,这就 是加权移动平均法的基本思想。
    $ D8 a: t4 V! g: R5 i0 Z
    9 k# i& {9 n8 n: k, C: v+ O, ^, L

    7 ~; S9 q. E8 G& D例 2  我国 1979~1988 年原煤产量如表 2 所示,试用加权移动平均法预测 1989 年 的产量7 e3 Q: m( y6 {9 ?' J( ?3 D% h+ `
    $ `- y6 s( g& z$ \9 q  y9 ~$ s, g

    % A: c5 @) Q* i; W. ]
    ( M0 V, h9 C$ D  U* h) S* {5 A3 L$ Y/ W" D! r
    2 k; @3 ]5 I+ K7 P

    0 D& E% L7 s/ i- K
    " Z8 ~0 s) v6 ]" K0 V, n8 S4 B# w2 O# A& Y3 B9 c. l+ {
    计算的 MATLAB 程序如下: ( g+ [8 `4 y8 m1 l

    % q/ y7 q8 V% ny=[6.35 6.20    6.22    6.66    7.15    7.89    8.72    8.94    9.28    9.8];
    7 V' P5 m5 n9 P! J' A' g' k! p' Ew=[1/6;2/6;3/6]; $ N6 R4 ]7 B- U  w# a
    m=length(y);n=3; ) y( y% Y, ^% \
    for i=1:m-n+1     0 z5 R) t' \4 Z3 V5 {1 t
        yhat(i)=y(i:i+n-1)*w;
    7 G$ Y4 X3 k+ H% f: z6 ]" d2 Nend
    ( L8 `# J9 S, Gyhat 0 h" ^. o1 g8 @( N9 O
    err=abs(y(n+1:m)-yhat(1:end-1))./y(n+1:m) " }0 `- z, R- `. c0 m4 A! D/ \  F
    T_err=1-sum(yhat(1:end-1))/sum(y(n+1:m)) 7 }2 i' v- x- _+ }/ t0 Y* w
    y1989=yhat(end)/(1-T_err)6 C8 O4 P& l0 y! a

    9 T, V! i- }3 T% p  J- k# j 在加权移动平均法中,   的选择,同样具有一定的经验性。一般的原则是:近期 数据的权数大,远期数据的权数小。至于大到什么程度和小到什么程度,则需要按照预 测者对序列的了解和分析来确定。 8 T* l. Z' H- y% [+ \, k

    * I" w. z& e2 k; T. y趋势移动平均法
    2 H& @; @% K. \/ W  k简单移动平均法和加权移动平均法,在时间序列没有明显的趋势变动时,能够准确 反映实际情况。但当时间序列出现直线增加或减少的变动趋势时,用简单移动平均法和 加权移动平均法来预测就会出现滞后偏差。因此,需要进行修正,修正的方法是作二次 移动平均,利用移动平均滞后偏差的规律来建立直线趋势的预测模型。这就是趋势移动平均法。  一次移动的平均数为 6 U+ X5 W% j3 @) S

    . f. K" {: I* @1 i" B' b: e3 v
    9 M0 p3 Y7 i! ?1 R  K1 u- p/ N( a7 t5 `; W

    $ x# w' Y5 Z3 j$ D1 q) b4 v( q2 P8 D* G7 _3 P

    & W" Q" m/ I, s/ K! l& a
    $ P" q; m- A3 z8 j例 3  我国 1965~1985 年的发电总量如表 3 所示,试预测 1986 年和 1987 年的发 电总量。 ; I/ `" F  z1 i7 j) X

    ; _5 f* u& \) N8 }' M2 b8 W, q
    2 z1 Z, t: K* W4 c
    ; ~+ _, i# l5 m4 `解  由散点图 1 可以看出,发电总量基本呈直线上升趋势,可用趋势移动平均法 来预测。
      ~* ^7 D% y5 e) ^& y: w" g( N4 m7 A) J0 B
    3 K  x5 t! g5 N. T0 U7 p0 T
    * p% x( A  _, F# o+ ^; u1 i# Z
    计算的 MATLAB 程序如下: * `7 I9 R9 @2 Q5 |& |" N+ N! n

    8 _& ]/ q/ {# ]2 d: ?clc,clear
    3 y" x2 K  A* p: y9 Q: U. ~$ }7 vload y.txt   %把原始数据保存在纯文本文件 y.txt 中 " q$ x2 a: ]3 x, Z( K- ]
    m1=length(y);    2 R0 S( s+ ?8 S, q. j% n
    n=6;   %n 为移动平均的项数 8 H, K: i$ W' W: `; H
    for i=1:m1-n+1     5 J: D" R2 T& i/ N" h
        yhat1(i)=sum(y(i:i+n-1))/n; 3 o- \4 ^) n4 y, e+ v4 N  _, s
    end 5 c  Z! ]/ `* s9 f- o5 t
    yhat1 ' x* E( p- b* _0 A4 T: Q" t
    m2=length(yhat1);
    / H8 |' Q% c/ z6 z- m0 Q" }for i=1:m2-n+1   
    , `- d( H. |% {5 p    yhat2(i)=sum(yhat1(i:i+n-1))/n;
    4 M4 N; S& U% M0 y8 o  h: jend
    & C# B* a) B3 g) u4 V) T: V9 eyhat2   # s/ X) I& h/ h: B
    plot(1:21,y,'*') ; k& u8 r; B) A! p5 _$ A
    a21=2*yhat1(end)-yhat2(end)
    * V9 f9 s: g7 Cb21=2*(yhat1(end)-yhat2(end))/(n-1) $ i0 \- e& s$ V3 ^9 T0 X. Z, c7 T
    y1986=a21+b21 7 Y( _* K' T. w2 b& I. Q9 B
    y1987=a21+2*b21
    8 K0 Z/ Q" s) Q, {2 B, a2 N( D! s: f, f/ d1 p5 P, p# y

    ) V( r" F% V# r/ z) E6 ^$ r) ^) N! a  J) c
    趋势移动平均法对于同时存在直线趋势与周期波动的序列,是一种既能反映趋势变 化,又可以有效地分离出来周期变动的方法。9 E! ~) v7 F, r5 L) y3 W6 I& p5 A
    3 `6 Z' n$ S& b1 t$ X; H' s
    5 m5 h- n7 H6 s1 t2 X7 C' ~* A
    " U; J/ c% R* F$ [+ ]4 K7 N
    ————————————————' B- Q/ I' ?/ b7 @4 P
    版权声明:本文为CSDN博主「wamg潇潇」的原创文章,遵循CC 4.0 BY-SA版权协议,转载请附上原文出处链接及本声明。5 J( M' I/ ~" t1 T  ^' _0 J5 P
    原文链接:https://blog.csdn.net/qq_29831163/article/details/89440426( J4 \9 F* U3 `8 ?. e, q

    7 A$ j3 d, u9 P$ b. N1 J
    8 B. c" ^( J% x, s
    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-12 14:22 , Processed in 0.579029 second(s), 52 queries .

    回顶部