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); " G7 _* R+ B+ q; W1 u# wyn=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) - R% f% S. k8 z+ w    for j=1:i 2 X; n) b; r2 B6 {        x1(i)=x1(i)+x0(j);     end 8 `4 R  `' B/ ^- i; uend : X, z( `! P" l- ?: K: hfor i=1:length(x0)-1 7 h& w) H0 Z, p) F    B(i,1)=(-1/2)*(x1(i)+x1(i+1)); 3 Z' ]3 L. @* M5 v& S    B(i,2)=1;     yn(i)=x0(i+1); end HatA=(inv(B'*B))*B'*yn % GM(1,1)模型参数估计 for k=1:length(x0)+T , q; ]# E! P& B    Hatx1(k)=(x0(1)-HatA(2)/HatA(1))*exp(-HatA(1)*(k-1))+HatA(2)/HatA(1); + E# L7 Z2 N' Z5 N, l$ M$ V5 _8 Dend Hatx0(1)=Hatx1(1); % q- u7 h7 @* F6 i: H) tfor k=2:length(x0)+T 9 q4 {' J- c1 K. G6 K3 Y    Hatx0(k)=Hatx1(k)-Hatx1(k-1);%累计还原得到历史数据的模拟值 end 1 |8 ^7 C, k7 e) B6 wfor i=1:length(x0) %开始模型检验     epsilon(i)=x0(i)-Hatx0(i);     omega(i)=(epsilon(i)/x0(i))*100; 1 \' b; v% L" b; @end 8 N( g3 O) k' A1 R+ F% x0;Hatx0;epsilon;omega;  %必要时去掉％得到各种数据 c=std(epsilon)/std(x0);p=0; ( E+ Z; C3 G4 {for i=1:length(x0)   m/ r6 u- V5 T    if abs(epsilon(i)-mean(epsilon))<0.6745*std(x0) : |3 @) D, A. J7 j5 ~8 {        p=p+1; # q$ @4 i0 t/ ], t# o    end ' |; |! Z) ]' y3 s  E0 X) G- t: |end p=p/length(x0) ' w$ i, D8 o* }! Z6 Qif p>0.95 & c<0.35     disp('The model is good,and the forecast is:'),   W2 }# G8 x. j' S0 t$ c& _    disp(Hatx0(length(x0)+T)) 8 X; b6 y% e9 b2 ]- A0 ^elseif p>0.85 & c<0.5     disp('The model is eligibility,and the forecast is:'),     disp(Hatx0(length(x0)+T)) ; V% X" b& g; s5 v4 Y/ C2 melseif p>0.7 & c>0.65     disp('The model is not good,and the forecast is:'),     disp(Hatx0(length(x0)+T)) else p<=0.7 & c>0.65 - r7 F$ A! h: l7 Q% g    disp('The model is bad and try again') end % l" s8 I3 p- g1 i9 Yfor i=1:length(x0)     Hatx00(i)=Hatx0(i); end z=1:length(x0); plot(z,x0,'-',z,Hatx00,':') %将原始数据和模拟值画在一个图上帮助观察 text(2,x0(2),'History data: real line') text(length(x0)/2,Hatx00(length(x0))/2,'Simulation data:broken line') " R, m. \* o, T) u* IendT=input('T=');%从键盘输入从最后一个历史数据算起的第T时点????是指什么啊，请大哥们，大姐们教一下，我急用，请快，谢谢我的初始值x0=[1.620938526 ; n& b# l% n- {0.07925621 0.052318818 0.041252502 0.021800479 5 |2 b3 e5 x/ e" x" k, p0.053132975 0.089908836 * `, O- k) N: H) o4 F0.109153219 0.079331832 0.342192598 & T& d) j; \# b* i. i5 t- h0.099718142 0.135194823 9 }4 a% }: A0 E. v0.109274037 6 Q5 L/ I  f: ]0.08152013 % O: ~8 }+ a) q$ E4 a# J0.067876355 0.064706843 0.055562197 + p, W" G; @2 e9 e0.050848544 * c7 _: V9 g$ J]';

daiqiang5566

楼主莫急..

yysclshi

a= -0.0080
u= 0.0713
4 V: o+ e( M2 p预测值
( c8 E  z6 H% k4 \, M# ~% l5 \    1.6209    0.0846    0.0853    0.0859    0.0866    0.0873
1 y+ V& F# Z2 ?    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
6 v5 i9 z: V  F9 ?& T    0.0899    0.1092    0.0793    0.3422    0.0997    0.1352
6 X% |" G) |) ?+ K7 |8 @    0.1093    0.0815    0.0679    0.0647    0.0556    0.0508
残差
         0   -0.0053   -0.0329   -0.0447   -0.0648   -0.0342
& O: p: \, h1 t6 _2 m4 h    0.0019    0.0204   -0.0101    0.2520    0.0088    0.0436
    0.0169   -0.0116   -0.0260   -0.0299   -0.0398   -0.0453
相对误差
4 O3 A4 F4 c$ |" n# |. ~         0    0.0672    0.6297    1.0835    2.9741    0.6437
    0.0209    0.1870    0.1276    0.7365    0.0885    0.3223
: P6 @+ v- K5 Q  J, `4 j    0.1548    0.1420    0.3826    0.4619    0.7162    0.8903
方差比
    0.1869
p =
$ M6 h* Y$ L& w6 v     1

yysclshi

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

杨晓敬

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

gxj820

很好啊！！顶顶！

一剑卡卡

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

jshzncd

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