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

xinzhiyong

GM(1,1)灰色模型的程序实现function GM1=fungry1(x0) %输入原始数据x0 6 A0 M/ Z( l7 o2 t9 P! S8 \) KT=input('T=');%从键盘输入从最后一个历史数据算起的第T时点 ! z. F2 p6 b; b  Z5 [x1=zeros(1,length(x0));B=zeros(length(x0)-1,2); yn=zeros(length(x0)-1,1);Hatx0=zeros(1,length(x0)+T); Hatx00=zeros(1,length(x0));Hatx1=zeros(1,length(x0)+T); epsilon=zeros(length(x0),1);omega=zeros(length(x0),1); for i=1:length(x0) 2 u4 D: o" p2 L; f    for j=1:i 8 P" u& [# f3 k9 @% K* P" ^/ ]6 N        x1(i)=x1(i)+x0(j);     end end for i=1:length(x0)-1     B(i,1)=(-1/2)*(x1(i)+x1(i+1));     B(i,2)=1; : Q$ J8 z3 V8 |3 k) _    yn(i)=x0(i+1); end % f9 ?4 f3 J' `HatA=(inv(B'*B))*B'*yn % GM(1,1)模型参数估计 2 j' h! d  d* j% p, a! i* Yfor k=1:length(x0)+T 5 r2 B8 R) Y* d, F    Hatx1(k)=(x0(1)-HatA(2)/HatA(1))*exp(-HatA(1)*(k-1))+HatA(2)/HatA(1); end / K: @# `+ {; A6 qHatx0(1)=Hatx1(1); / L' e. V% T  _7 [) a9 Jfor k=2:length(x0)+T 5 e& X4 q9 K# I' i; U    Hatx0(k)=Hatx1(k)-Hatx1(k-1);%累计还原得到历史数据的模拟值 ( ]& n1 [+ r' W9 b. kend 0 M/ t  c! K4 D" {; X0 D9 T2 }( wfor i=1:length(x0) %开始模型检验     epsilon(i)=x0(i)-Hatx0(i);     omega(i)=(epsilon(i)/x0(i))*100; end % x0;Hatx0;epsilon;omega;  %必要时去掉％得到各种数据 c=std(epsilon)/std(x0);p=0; ( {6 U7 p) w) ]& C8 W" o0 Nfor i=1:length(x0) 2 C' r; W( Q% T; j3 p' N. L    if abs(epsilon(i)-mean(epsilon))<0.6745*std(x0) & Q. c( d9 F; k  \- l. j        p=p+1; 3 i2 {( O7 a7 I; I5 `- @    end ' ^2 U1 @6 J$ J8 B0 X3 ]end . _% e' ~: T0 F" p. n( jp=p/length(x0) % h% H, |  j7 x7 ~8 q  I) cif p>0.95 & c<0.35 3 z2 W9 m  S( k8 C% _    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:'), $ ]9 l7 s% t" i$ o! y1 S, \3 [% H    disp(Hatx0(length(x0)+T)) 6 d: H* X$ s5 d9 Oelseif p>0.7 & c>0.65     disp('The model is not good,and the forecast is:'), . i, g+ c% x8 b0 D+ @4 |    disp(Hatx0(length(x0)+T)) else p<=0.7 & c>0.65     disp('The model is bad and try again') end & V) E5 F* T8 sfor i=1:length(x0)     Hatx00(i)=Hatx0(i); + j* Y) u- u: x: q- z/ ^& [: {end # A- V6 e8 D. jz=1:length(x0); plot(z,x0,'-',z,Hatx00,':') %将原始数据和模拟值画在一个图上帮助观察 : p  b& x5 w" z1 ^5 Ltext(2,x0(2),'History data: real line') 7 W7 a# M1 M7 q4 D" v' A% Ctext(length(x0)/2,Hatx00(length(x0))/2,'Simulation data:broken line') * Y$ P9 u" I! `: B( b. j! RendT=input('T=');%从键盘输入从最后一个历史数据算起的第T时点????是指什么啊，请大哥们，大姐们教一下，我急用，请快，谢谢我的初始值x0=[1.620938526 8 B* v* ^0 x  _6 ~. P- j0.07925621 # G* X: l9 T) N/ C, M0.052318818 % _) A0 R+ x0 \8 `0 P3 f$ E6 G' O0.041252502 0.021800479 6 S9 h' f- s# X  d! ^8 p# K0.053132975 0.089908836 8 Q3 a, `0 R8 q, V: k0.109153219 0.079331832 0.342192598 0.099718142 4 n- d2 L$ x. p* V- I  ]1 B3 p0.135194823 0.109274037 0.08152013 2 A! ^9 T9 y, A6 y0.067876355 0.064706843 0.055562197 0.050848544 0 L  O* B% b4 u: C; K]';

daiqiang5566

楼主莫急..

yysclshi

a= -0.0080
' p( O& t4 Z/ Z3 W1 T& o2 Nu= 0.0713
8 x0 t) x. M8 ]; U5 o$ B预测值
    1.6209    0.0846    0.0853    0.0859    0.0866    0.0873
/ s$ K' i9 ^+ M6 m4 Q$ ?    0.0880    0.0887    0.0895    0.0902    0.0909    0.0916
' e) l; U# A) i    0.0924    0.0931    0.0938    0.0946    0.0954    0.0961
- U1 @1 F: ^4 Z( u  W初始值
5 R% u/ m  [5 H! ^+ N# j    1.6209    0.0793    0.0523    0.0413    0.0218    0.0531
* X0 [5 E$ z; ?7 D: {    0.0899    0.1092    0.0793    0.3422    0.0997    0.1352
8 g; y& Y0 _2 @" P3 y1 i) G5 S6 ~    0.1093    0.0815    0.0679    0.0647    0.0556    0.0508
残差
         0   -0.0053   -0.0329   -0.0447   -0.0648   -0.0342
6 i7 {% B) F, [2 W    0.0019    0.0204   -0.0101    0.2520    0.0088    0.0436
    0.0169   -0.0116   -0.0260   -0.0299   -0.0398   -0.0453
相对误差
/ \. n1 Z3 J- l" D; Z, `- U         0    0.0672    0.6297    1.0835    2.9741    0.6437
    0.0209    0.1870    0.1276    0.7365    0.0885    0.3223
    0.1548    0.1420    0.3826    0.4619    0.7162    0.8903
方差比
    0.1869
p =
9 h1 f  B3 ]- {  W6 ^     1

yysclshi

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

杨晓敬

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

gxj820

很好啊！！顶顶！

一剑卡卡

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

jshzncd

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