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); yn=zeros(length(x0)-1,1);Hatx0=zeros(1,length(x0)+T); Hatx00=zeros(1,length(x0));Hatx1=zeros(1,length(x0)+T); 2 W$ {* b/ g. l/ B4 }epsilon=zeros(length(x0),1);omega=zeros(length(x0),1); & V  g- D6 F& Y! Mfor i=1:length(x0) , b! Q* {4 A# I# u, f    for j=1:i ' w+ P) W# e3 J. C- F        x1(i)=x1(i)+x0(j);     end end for i=1:length(x0)-1 " X9 u2 A6 |+ j: B8 h. A* v    B(i,1)=(-1/2)*(x1(i)+x1(i+1));     B(i,2)=1;     yn(i)=x0(i+1); 0 }+ _+ j7 g  h' r+ Fend ! [: N% u: K% L: D. J# R/ RHatA=(inv(B'*B))*B'*yn % GM(1,1)模型参数估计 for k=1:length(x0)+T     Hatx1(k)=(x0(1)-HatA(2)/HatA(1))*exp(-HatA(1)*(k-1))+HatA(2)/HatA(1); end Hatx0(1)=Hatx1(1); for k=2:length(x0)+T     Hatx0(k)=Hatx1(k)-Hatx1(k-1);%累计还原得到历史数据的模拟值 ' F7 N5 u' c: R! A  Pend for i=1:length(x0) %开始模型检验 8 ~  |& T5 A+ i9 F2 F    epsilon(i)=x0(i)-Hatx0(i);     omega(i)=(epsilon(i)/x0(i))*100; 1 t# N+ q& j9 a6 |4 J; eend % x0;Hatx0;epsilon;omega;  %必要时去掉％得到各种数据 7 u9 f" c3 p% F' lc=std(epsilon)/std(x0);p=0; ! j! [4 f* ]; B! c* q! qfor i=1:length(x0) . c: h0 A) n! i6 l4 e    if abs(epsilon(i)-mean(epsilon))<0.6745*std(x0) 1 H. o* H6 r+ [' g: v        p=p+1; 4 a  Z; w% U  @    end 3 j) a( [& a- E( V$ z8 K7 ]8 aend 6 E; M6 L" e1 Q/ t1 o' wp=p/length(x0) if p>0.95 & c<0.35     disp('The model is good,and the forecast is:'),     disp(Hatx0(length(x0)+T)) - w- z7 v- t6 _" \elseif p>0.85 & c<0.5     disp('The model is eligibility,and the forecast is:'), - d  K8 o. m2 T3 y    disp(Hatx0(length(x0)+T)) 5 e1 X& u8 M0 k9 G9 Qelseif p>0.7 & c>0.65     disp('The model is not good,and the forecast is:'),     disp(Hatx0(length(x0)+T)) ; F' U  o  g- ?9 s0 L1 h& ^6 V+ \else p<=0.7 & c>0.65 6 Z' F* Q* k( I% p    disp('The model is bad and try again') end for i=1:length(x0)     Hatx00(i)=Hatx0(i); end ; f0 N" Y7 H, ^# F/ x; Zz=1:length(x0); plot(z,x0,'-',z,Hatx00,':') %将原始数据和模拟值画在一个图上帮助观察 " x) }! [, w/ Y; D' wtext(2,x0(2),'History data: real line') text(length(x0)/2,Hatx00(length(x0))/2,'Simulation data:broken line') 9 g8 u( k; d3 H2 MendT=input('T=');%从键盘输入从最后一个历史数据算起的第T时点????是指什么啊，请大哥们，大姐们教一下，我急用，请快，谢谢我的初始值x0=[1.620938526 1 D8 l" z) k; s. d, I: c3 q/ t* I; v0.07925621 0.052318818 0.041252502   g* I! C1 w9 k0 j9 B1 x, s0.021800479 ' ?$ O/ u, P5 g& I+ O0.053132975 0.089908836 0.109153219 % r5 Z! x6 M6 t0.079331832 0.342192598 0.099718142 0 ]2 n9 p' r8 P: r; k6 V8 J4 C; g0.135194823 0.109274037 0.08152013 3 p0 C1 r# H: q0 Q: o0.067876355 0.064706843 0.055562197 0.050848544 ! f+ X+ e) m! z2 R9 s8 e]';

daiqiang5566

楼主莫急..

yysclshi

a= -0.0080
u= 0.0713
3 Z# Y+ z- X/ }, \1 d. d4 }# v! a预测值
: I& v3 l1 [4 U7 z    1.6209    0.0846    0.0853    0.0859    0.0866    0.0873
! e8 N  f8 x, `9 d1 i$ {$ e    0.0880    0.0887    0.0895    0.0902    0.0909    0.0916
' K/ P: h- q6 ~8 Y4 \% ^' S/ c$ i    0.0924    0.0931    0.0938    0.0946    0.0954    0.0961
% f9 p/ {- Z6 q- Q+ Y: v+ ?# r初始值
" E8 t, c+ u$ Q    1.6209    0.0793    0.0523    0.0413    0.0218    0.0531
    0.0899    0.1092    0.0793    0.3422    0.0997    0.1352
    0.1093    0.0815    0.0679    0.0647    0.0556    0.0508
残差
( `1 H" ]9 s+ ^) i/ Q         0   -0.0053   -0.0329   -0.0447   -0.0648   -0.0342
( R& f8 a4 M+ w* N, x. P    0.0019    0.0204   -0.0101    0.2520    0.0088    0.0436
; V' k5 a0 W, d- p9 v    0.0169   -0.0116   -0.0260   -0.0299   -0.0398   -0.0453
3 S- m( }( o9 A+ e相对误差
         0    0.0672    0.6297    1.0835    2.9741    0.6437
    0.0209    0.1870    0.1276    0.7365    0.0885    0.3223
9 ]( X* J( P; m( D# Z% L) `0 K  X3 C- Q% a    0.1548    0.1420    0.3826    0.4619    0.7162    0.8903
: ^$ M1 [% B! l方差比
    0.1869
p =
/ _  ^0 o+ e* a! |' ^     1

yysclshi

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

杨晓敬

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

gxj820

很好啊！！顶顶！

一剑卡卡

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

jshzncd

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