GM(1,1)预测模型的MATLAB程序求助，急！！！

xinzhiyong

GM(1,1)灰色模型的程序实现function GM1=fungry1(x0) %输入原始数据x0 T=input('T=');%从键盘输入从最后一个历史数据算起的第T时点 x1=zeros(1,length(x0));B=zeros(length(x0)-1,2); + o& [9 R2 b4 s5 I3 ~yn=zeros(length(x0)-1,1);Hatx0=zeros(1,length(x0)+T); ) U9 _: ~) Z8 S! oHatx00=zeros(1,length(x0));Hatx1=zeros(1,length(x0)+T); epsilon=zeros(length(x0),1);omega=zeros(length(x0),1); 4 C: Z& C- E) r1 w# ^for i=1:length(x0)     for j=1:i         x1(i)=x1(i)+x0(j); 0 S7 l: Z) T2 H- x    end 2 |, n( l. T: s8 e" }: g: |" Pend for i=1:length(x0)-1 $ U/ r3 s2 @' P% u/ i9 T- v  q6 N& C    B(i,1)=(-1/2)*(x1(i)+x1(i+1));     B(i,2)=1;     yn(i)=x0(i+1); end : j) l% J8 }  H, ]1 P; DHatA=(inv(B'*B))*B'*yn % GM(1,1)模型参数估计 1 D; Y0 B# b: D6 r1 q2 L0 Ofor k=1:length(x0)+T & j# G6 _' z3 L0 b' {/ j7 `    Hatx1(k)=(x0(1)-HatA(2)/HatA(1))*exp(-HatA(1)*(k-1))+HatA(2)/HatA(1); 3 ~6 S3 j% O" d+ Wend 5 k8 N, r" |3 w5 _% r: u  x- @Hatx0(1)=Hatx1(1); 2 {( v. @0 s+ ^+ _/ tfor k=2:length(x0)+T     Hatx0(k)=Hatx1(k)-Hatx1(k-1);%累计还原得到历史数据的模拟值 end for i=1:length(x0) %开始模型检验 $ t  }' J9 \* i% [. M    epsilon(i)=x0(i)-Hatx0(i);     omega(i)=(epsilon(i)/x0(i))*100; end 6 t, U. z0 [/ ^' O: O' `  S% x0;Hatx0;epsilon;omega;  %必要时去掉％得到各种数据 c=std(epsilon)/std(x0);p=0; for i=1:length(x0)     if abs(epsilon(i)-mean(epsilon))<0.6745*std(x0) . T/ e+ J& t1 A8 m" h9 N  n        p=p+1;     end 0 u5 l" O9 D" o) ^0 eend p=p/length(x0) if p>0.95 & c<0.35 . |* o2 c$ ^& I+ p% _' l    disp('The model is good,and the forecast is:'),     disp(Hatx0(length(x0)+T)) elseif p>0.85 & c<0.5     disp('The model is eligibility,and the forecast is:'),     disp(Hatx0(length(x0)+T)) / X/ [' m2 b$ h  I4 telseif p>0.7 & c>0.65 5 W0 L, D, k+ I" q; P    disp('The model is not good,and the forecast is:'), - l9 N" M3 E4 ]) X$ G7 i) Q, A2 K7 ]- ?    disp(Hatx0(length(x0)+T)) $ L- O# \6 i/ x. o" z& Y6 celse p<=0.7 & c>0.65     disp('The model is bad and try again') " ~4 p* ]% j& @  jend for i=1:length(x0) 6 B9 c- _7 ]# j; U3 I( x    Hatx00(i)=Hatx0(i); end $ h& P; W9 q7 q0 Z! C% \* s1 @z=1:length(x0); plot(z,x0,'-',z,Hatx00,':') %将原始数据和模拟值画在一个图上帮助观察 text(2,x0(2),'History data: real line')   X; u: p! p' S% ^  z" w4 V2 a- jtext(length(x0)/2,Hatx00(length(x0))/2,'Simulation data:broken line') - u: o+ u, M  x( ^8 T- o- XendT=input('T=');%从键盘输入从最后一个历史数据算起的第T时点????是指什么啊，请大哥们，大姐们教一下，我急用，请快，谢谢我的初始值x0=[1.620938526 : S; F) ]5 r; z1 \0.07925621 0.052318818 . D3 f5 p7 Q; y  M0.041252502 0.021800479 0.053132975 8 C$ A, w6 ?0 J  ^( C4 R0.089908836 ) C4 j- i' @6 E. @9 r- [- \+ p( H$ P3 `0.109153219 0 `6 T4 H. @- P/ X0.079331832 0.342192598 0.099718142 0.135194823 5 Y8 k: D2 P1 h/ X' z0.109274037 0.08152013 ' A& L* S6 B3 c4 C4 d0.067876355 0.064706843 - ]( K0 ~: F! J2 M6 I( A1 U0.055562197 0.050848544 % X! A( `2 K/ y8 o! u' a( F9 X]';

daiqiang5566

楼主莫急..

yysclshi

a= -0.0080
u= 0.0713
2 u5 o$ g( X* h/ o* \8 d预测值
) J; o- Y/ |$ \- k: x! A    1.6209    0.0846    0.0853    0.0859    0.0866    0.0873
    0.0880    0.0887    0.0895    0.0902    0.0909    0.0916
    0.0924    0.0931    0.0938    0.0946    0.0954    0.0961
初始值
    1.6209    0.0793    0.0523    0.0413    0.0218    0.0531
2 q, u. \* p& `# h' h/ U. i    0.0899    0.1092    0.0793    0.3422    0.0997    0.1352
& T# j; k, S: ?5 j    0.1093    0.0815    0.0679    0.0647    0.0556    0.0508
  J9 T' q+ v% F: Q# q( R残差
1 L  @% F' F- H: m- B7 }         0   -0.0053   -0.0329   -0.0447   -0.0648   -0.0342
- I) M' I- I; G. h& [; h1 j( g    0.0019    0.0204   -0.0101    0.2520    0.0088    0.0436
- {3 p% S( n; c  \    0.0169   -0.0116   -0.0260   -0.0299   -0.0398   -0.0453
相对误差
2 e' s) O7 m; R: e7 `; l1 ]         0    0.0672    0.6297    1.0835    2.9741    0.6437
$ j4 |, Q; u3 x! U. V4 ]    0.0209    0.1870    0.1276    0.7365    0.0885    0.3223
    0.1548    0.1420    0.3826    0.4619    0.7162    0.8903
5 K$ |6 @+ N( O) V+ k/ G- x! L* H. e方差比
  A  v. L6 M+ C( g, C8 X4 M3 `    0.1869
p =
* h* f; }3 ~4 I' G     1

yysclshi

这个图片我不会发，真难为情

杨晓敬

不好意思，不会啊！帮不上忙啊

gxj820

很好啊！！顶顶！

一剑卡卡

不懂。。。。。。。。。。。。。。

jshzncd

不懂啊，楼主~不好意思是哈