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

xinzhiyong

GM(1,1)灰色模型的程序实现function GM1=fungry1(x0) %输入原始数据x0 T=input('T=');%从键盘输入从最后一个历史数据算起的第T时点 * d! O3 G4 l8 zx1=zeros(1,length(x0));B=zeros(length(x0)-1,2); yn=zeros(length(x0)-1,1);Hatx0=zeros(1,length(x0)+T); 4 u0 b3 D) v& Q0 M/ `, [6 uHatx00=zeros(1,length(x0));Hatx1=zeros(1,length(x0)+T); # S) c7 H# g6 F$ b* p4 V/ J( c7 ^: Mepsilon=zeros(length(x0),1);omega=zeros(length(x0),1); for i=1:length(x0)     for j=1:i 1 k% Y* V9 p) y; o$ A        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; + {' n1 V4 t9 O- U* A8 p! x    yn(i)=x0(i+1); 6 L- r: o, [! B1 C8 X& j+ Rend . r  R. s: ]) j1 M( JHatA=(inv(B'*B))*B'*yn % GM(1,1)模型参数估计 for k=1:length(x0)+T 5 E6 [2 ]$ K# [2 {    Hatx1(k)=(x0(1)-HatA(2)/HatA(1))*exp(-HatA(1)*(k-1))+HatA(2)/HatA(1); 6 W& F3 I5 y6 P8 iend - W6 m$ c, M; H# t7 E$ zHatx0(1)=Hatx1(1); for k=2:length(x0)+T     Hatx0(k)=Hatx1(k)-Hatx1(k-1);%累计还原得到历史数据的模拟值 " H: o' C: N6 ~# p" ~' G; F' q2 l6 send / @! C7 _' E$ J$ q% g$ U; Lfor i=1:length(x0) %开始模型检验 - W. {7 N. h/ R, Q    epsilon(i)=x0(i)-Hatx0(i);     omega(i)=(epsilon(i)/x0(i))*100; end - M5 o* E- u% G* W6 x: Q3 s+ Z% \% x0;Hatx0;epsilon;omega;  %必要时去掉％得到各种数据 c=std(epsilon)/std(x0);p=0; : u! f6 I! t) D% E. x* Lfor i=1:length(x0) 5 {: z6 B* v3 t0 R& z" p8 p    if abs(epsilon(i)-mean(epsilon))<0.6745*std(x0)         p=p+1; ( X* }! h: r2 W' p- p    end end 0 q4 N7 @; M$ C' }p=p/length(x0) 2 @6 r. N2 t, f) eif p>0.95 & c<0.35     disp('The model is good,and the forecast is:'), : N* ?9 A% [. U8 o* c    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)) elseif p>0.7 & c>0.65 3 \7 ^8 M/ |0 ^, j( Z    disp('The model is not good,and the forecast is:'), : {8 o- ?$ W: D" }& b    disp(Hatx0(length(x0)+T)) ; L$ h& A- x( Q- delse p<=0.7 & c>0.65 , D8 N9 R! u) u# b& r. i! N! }5 F    disp('The model is bad and try again') 2 ?$ Z9 W- L. `, oend ; X% D$ a2 W- y1 ^5 O& ffor i=1:length(x0)     Hatx00(i)=Hatx0(i); end z=1:length(x0); : _+ k6 _8 i3 O! I: O+ kplot(z,x0,'-',z,Hatx00,':') %将原始数据和模拟值画在一个图上帮助观察 % C& b4 y! i9 q2 `' a9 D/ r5 Ntext(2,x0(2),'History data: real line') text(length(x0)/2,Hatx00(length(x0))/2,'Simulation data:broken line') . r6 L' o8 L4 p9 v( c& L5 nendT=input('T=');%从键盘输入从最后一个历史数据算起的第T时点????是指什么啊，请大哥们，大姐们教一下，我急用，请快，谢谢我的初始值x0=[1.620938526 0.07925621 , ^% c# q+ [2 T2 |2 Y, L0.052318818 , {- w: J  U, h4 Q4 R0.041252502 0.021800479 0.053132975   ]1 K+ L9 V* n5 N* @* v! b0.089908836 3 g' ], t9 K1 S0.109153219 ( M% h" \- u5 ?0.079331832 2 z# |6 A5 j# j& ?# H0 ~5 |% o0.342192598 0.099718142 0 m1 B; H/ F8 O+ y0 J0.135194823 0.109274037 0 }5 [' X# @+ z% Q0.08152013 0.067876355 0.064706843 & F& g; @: x0 E# N1 H$ s: S  @, K0.055562197 0.050848544 ]';

daiqiang5566

楼主莫急..

yysclshi

a= -0.0080
u= 0.0713
& [/ i4 C& p" g  M1 z预测值
    1.6209    0.0846    0.0853    0.0859    0.0866    0.0873
. N- Q5 B0 R9 Z  r1 t: v1 B  b    0.0880    0.0887    0.0895    0.0902    0.0909    0.0916
% }$ O; v+ `4 @  ~7 {    0.0924    0.0931    0.0938    0.0946    0.0954    0.0961
初始值
    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
残差
& f* d( ]+ T: R' @- f. q$ q         0   -0.0053   -0.0329   -0.0447   -0.0648   -0.0342
    0.0019    0.0204   -0.0101    0.2520    0.0088    0.0436
    0.0169   -0.0116   -0.0260   -0.0299   -0.0398   -0.0453
相对误差
         0    0.0672    0.6297    1.0835    2.9741    0.6437
6 Q- ?! X6 _8 u+ H    0.0209    0.1870    0.1276    0.7365    0.0885    0.3223
! q! ?( m: ?( W+ z/ {    0.1548    0.1420    0.3826    0.4619    0.7162    0.8903
- F8 n8 w4 S( k9 X* J3 A方差比
; {8 i. a  s# @2 l7 t  `4 M    0.1869
+ R" \3 b+ x& K$ @0 y6 Xp =
     1

yysclshi

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

杨晓敬

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

gxj820

很好啊！！顶顶！

一剑卡卡

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

jshzncd

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