数学建模-SARS疫情对某些经济指标影响例题

1047521767

            数学建模-SARS疫情对某些经济指标影响例题目录

• 一、问题
• 二、对例题如何评估分析思路
• • 1.第一种
  • 2.第二种
  • 3.第三种
    9 R9 K% @' Q& ?5 i  M# v- a' N
• 三、模型的分析和预测
• • 1.给出假设
  • 2.建立灰色预测模型GM(1,1)
  • 3.模型的求解
  • • a为什么取0.4？
    • z(1)(k)是怎么得出的？
      , I# D1 k) Y, T2 K0 s( B7 Z/ p9 w# w
    
    ' z$ ~: ]7 x$ ?
• 四、编程实现
• 一、问题
• 
  
• 表1：
  
• 
  + H6 |( g7 S7 k2 f9 I

• 2003年的SARS病毒，我们只需要根据几年前的数据预测出2003年在正常情况下的数据，然后和实际2003年的数据做出比较，就能进行评估

  表3：

  
  2 z' e2 z$ |( m% S

• 
  

• 二、对例题如何评估分析思路

  1.第一种

  我们只需要找到商品零售业、旅游业和综合服务业的历年的的基本水平，然后分析一下如果按照正常的模式，如果没有疫情，在未来的某一年应该达到怎样的水平，然后再和实际的水平做一个比较，就可以得到定量的分析

  2.第二种

  我们也可以根据每年的某个月来预测2003年的某个月的数据，比如用每年2月的数据来预测2003年2月的数据

  3.第三种

  我们可以对每年每个月的和进行预测，然后如果要预测某个月，把这个月的数据每年占总和的比例进行预测，然后用2003年的某个月的比例乘以总和就是这个月的数据

  
  

  三、模型的分析和预测

• 1.给出假设
• 
  
• 2.建立灰色预测模型GM(1,1)
  
• 
  
• 
  
• 
  7 O& J; n9 Q: U
• 3.模型的求解
  
• 
    g6 N0 j0 t7 o0 n

• 这里平均值是每年12个月的平均值，总共六年，六个数据

  a为什么取0.4？

  a先取0.5，如果发现经不起检验就适当的调高或者调低，但是范围在0~1里面

  z(1)(k)是怎么得出的？

  z(1)(k)=0.4*x(1)(k)+(1-0.4)*x(1)(k-1)

  
  
• 四、编程实现
• 
  . b! j& B$ @1 t% P$ u
• 
  3 y$ Y" {6 d! c
• 代码演示：
  
• 
  clc,clear
  han1=[83.0,79.8,78.1,85.1,86.6,88.2,90.3,86.7,93.3,92.5,90.9,96.9;
      101.7,85.1,87.8,91.6,93.4,94.5,97.4,99.5,104.2,102.3,101.0,123.5;
  3 z5 l8 A: _/ ?7 O; k    92.2,114.0,93.3,101.0,103.5,105.2,109.5,109.2,109.6,111.2,121.7,131.3;
  2 y* V$ ~0 A- [    105.0,125.7,106.6,116.0,117.6,118.0,121.7,118.7,120.2,127.8,121.8,121.9;
      139.3,129.5,122.5,124.5,135.7,130.8,138.7,133.7,136.8,138.9,129.6,133.7;
  9 ~" [7 w9 ~$ }; k    137.5,135.3,133.0,133.4,142.8,141.6,142.9,147.3,159.6,162.1,153.5,155.9;
  " X9 D5 W: m0 n; Q. `! y    163.2,159.7,158.4,145.2,124,144.1,157,162.6,171.8,180.7,173.5,176.5]
  4 K& K$ W8 d  ?0 Y2 D  ^han1(end,=[];%相当于han1=han1(1:6;
  m=size(han1,2);%把月份提取出来
  : j; e' D. S/ }6 T8 I  [) px0=mean(han1,2);%返回x矩阵每行的平均值，其中的2代表返回行
  $ \1 r9 w% ~. C; S5 \7 e* Rx1=cumsum(x0)%一次累加
  alpha=0.4;n=length(x0);
  z1=alpha*x1(2:n)+(1-alpha)*x1(1:n-1)%求邻域生成数
  : h- N/ @; C* H% m* _9 t. A- `" k* qY=x0(2:n);B=[-z1,ones(n-1,1)];
  ab=B\Y
  5 _  F  s8 `! x8 u' ^k=6;
  $ ^" }1 y: y% m7 kx7hat=(x0(1)-ab(2)/ab(1))*(exp(-ab(1)*k)-exp(-ab(1)*(k-1)))
  + s2 O! z6 U2 j, Nz=m*x7hat
  $ {8 I1 c. ~8 @6 N$ S" au=sum(han1)/sum(sum(han1))  %sum(han1)是每一列的数据，sum(sum(han1))是每一列加每一行总的数据，相当于每一月的均值
  v=z*u
  4 f$ P# g6 J8 }
  2 h1 n5 e9 X4 J7 p5 J7 N
  

• 输出：
  
• han1 =
  
  0 q9 e( r4 m2 Z( f
    列 1 至 11
  7 S) x# h0 L. s& {; ^4 Y$ J8 @: X) o
  5 \) o4 H+ n6 |  S6 i7 c# u: d
  ( X  Z) l+ T' E+ U2 G9 V% D; z   83.0000   79.8000   78.1000   85.1000   86.6000   88.2000   90.3000   86.7000   93.3000   92.5000   90.9000
  7 X+ }9 e' q! l+ ?  101.7000   85.1000   87.8000   91.6000   93.4000   94.5000   97.4000   99.5000  104.2000  102.3000  101.0000
  1 W1 Y) |. E+ P  C5 E   92.2000  114.0000   93.3000  101.0000  103.5000  105.2000  109.5000  109.2000  109.6000  111.2000  121.7000
    105.0000  125.7000  106.6000  116.0000  117.6000  118.0000  121.7000  118.7000  120.2000  127.8000  121.8000
    139.3000  129.5000  122.5000  124.5000  135.7000  130.8000  138.7000  133.7000  136.8000  138.9000  129.6000
  & o' t6 J/ j5 K  137.5000  135.3000  133.0000  133.4000  142.8000  141.6000  142.9000  147.3000  159.6000  162.1000  153.5000
    163.2000  159.7000  158.4000  145.2000  124.0000  144.1000  157.0000  162.6000  171.8000  180.7000  173.5000
  
  
    列 12
  ' J! d, M/ U, o& ?; Y8 A' C8 |
  # I; B8 X9 i, b% J# ]6 p3 m0 Y
  ' o9 ]& W8 X: F* v1 y   96.9000
  3 K( ^, x  y) o! i' y: R* x  123.5000
  + |  [1 f4 y% l% K5 N  131.3000
    121.9000
  - k7 s/ @8 W9 P  133.7000
    155.9000
  + K* M7 W7 @" T3 m8 c: x/ O5 _  176.5000
  ) t1 j' R) X, a0 C  h6 F
  
  1 s/ }" i7 O4 P1 v+ c: P
  & ?" s% Z. |/ f3 Z; C2 G
  x1 =
  / b! Y2 `1 l9 E% O
  
     87.6167
    186.1167
  ( w3 S( ?) ?9 F6 q" a  294.5917
    413.0083
    545.8167
  8 E0 q$ _, O+ |+ P  691.2250
  
  0 Z& _4 P" e6 t: D
  ( G# a4 R2 i6 _1 J9 c
  
  ) Y+ ^8 d# [: d% h7 o& Zz1 =
  
  - v3 M" u6 n7 P4 b1 L* c
    127.0167
    229.5067
  . g) {  z- K6 P1 E% p- Z  |5 w/ O) i  341.9583
    466.1317
    603.9800
  
  ; H; E5 p/ r/ ~/ k0 T) j2 ]
  " k( Y$ ?$ T9 p1 O7 y
  : R+ J% @6 l& x7 V' @
  ab =
  
  
  3 D& D. o7 U+ g* h- e   -0.0993
  % Z' G- x9 ~$ W  Z. d) j% _- o   85.5985
  # h  K# _* A- t/ c& Z
  
  
  
  x7hat =
  
  
    {7 {- {' H, X- d9 a  162.8793
  ' l5 s) w! n* @6 L) @
  
  
  
  z =
  
  . W) Q' @% ~+ P$ }6 I0 h* I; f
  6 _  j6 ~  }7 Q9 h   1.9546e+03
  
  
  
  ' n5 n$ C6 u5 Q
  - D, j7 I& q2 S& }u =
  
  ! e+ x' f7 ]: H2 q, B6 q
  + D* B9 _3 a4 k5 p  a  列 1 至 11
  
  ( V+ }/ d1 _) G' @! d: P- N" r
  - k" D- `- _8 @- j    0.0794    0.0807    0.0749    0.0786    0.0819    0.0818    0.0845    0.0838    0.0872    0.0886    0.0866
  / F; J# J; e8 R  N
  
    @& g6 A9 m+ x' `4 p: i  列 12
  
  - n, A8 k% h: o* J8 f
  - u+ P7 V4 Q5 O% t    0.0920
  
  , Q+ W$ t- ^! N: M* x2 L
  
  ! h2 K8 d. V+ C/ @  D9 c( U
  $ A' [* n, Y" t7 S* n$ y1 G# V9 rv =
  
  0 |3 }7 Y/ E  C2 U1 j! P2 C. ^- ~& a
  ) I) j7 c4 @7 P" N  列 1 至 11
  
    h# ~, U/ k4 u4 p- y
  8 w1 y6 x5 h! M) m' F0 W# z  155.2152  157.7365  146.4023  153.5421  160.1400  159.8337  165.0649  163.7924  170.5317  173.1473  169.3064
  , ?, J/ l, @& {5 O4 ^$ c/ Z2 A3 N/ F
  
  0 Y' b, Q" {& X1 h  列 12
  
  
  1 Q& y5 H  v! c; r( {* j' L1 }# c  179.8394
  ) n: z* j2 r1 T, t+ e0 j; j
  
  
• 
  6 ?, G0 ?6 P" j2 [1 l1 b
  + u" |5 E# w" G
• 我们可以得出
  0 d' v; ~7 k6 h- z3 @
• 
  
  6 h8 j# q% l6 c) \1 U0 z4 `

v就是2003年的实际数据
+ m# W7 r7 Q9 H9 b( f$ A6 E3 V" {% G

• 
  ( I& X( Z: f% E% L; L9 Y% H2 [- J
  " U2 U7 J, A& @

) m( e; R  O8 z6 e& }; e3 v

sjlxdn

11111111111111
* {. l  s3 Q' x/ Y; C3 n! g