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

xinzhiyong

GM(1,1)灰色模型的程序实现function GM1=fungry1(x0) %输入原始数据x0 ) w) D9 Y3 w- p% BT=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); ! I; K5 C* _4 M, N5 t0 \Hatx00=zeros(1,length(x0));Hatx1=zeros(1,length(x0)+T); * v, Z3 J2 O! M0 i. V" {* \epsilon=zeros(length(x0),1);omega=zeros(length(x0),1); for i=1:length(x0)     for j=1:i / I  Q$ X5 C$ n        x1(i)=x1(i)+x0(j);     end end for i=1:length(x0)-1 . Q6 q3 P3 n0 k/ H    B(i,1)=(-1/2)*(x1(i)+x1(i+1)); " C7 M- k$ B' r. z% G/ G% T    B(i,2)=1; ; o1 I* Z/ j4 O0 O    yn(i)=x0(i+1); 4 u+ s, t; s: O% Send HatA=(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); 7 S1 `' c$ r; `0 h- ?  oend ) Q# Z) k1 @" bHatx0(1)=Hatx1(1); + v% D. t. D! Mfor k=2:length(x0)+T 1 K& a, X! D' s( @1 [/ U    Hatx0(k)=Hatx1(k)-Hatx1(k-1);%累计还原得到历史数据的模拟值 ; ~4 D8 b' ?; Y3 _! ]) {& E! fend for i=1:length(x0) %开始模型检验     epsilon(i)=x0(i)-Hatx0(i); / E4 ~6 u2 N2 x9 H3 h$ _    omega(i)=(epsilon(i)/x0(i))*100; end + M7 S) d, O. T) Y$ Y3 j6 _2 S% x0;Hatx0;epsilon;omega;  %必要时去掉％得到各种数据 6 @2 J, l8 q* _& Gc=std(epsilon)/std(x0);p=0; ; |' W4 u+ D- _. }; M$ r7 Z7 Cfor i=1:length(x0) 1 l8 A. e/ v% _4 l( o4 T6 I    if abs(epsilon(i)-mean(epsilon))<0.6745*std(x0)         p=p+1; ' Z4 L* Q0 N- N8 \2 |& z% y    end end ' Y1 o1 b- F, ?* o" b$ \p=p/length(x0) & a# r6 U, z/ G7 ]8 d1 Aif p>0.95 & c<0.35     disp('The model is good,and the forecast is:'),     disp(Hatx0(length(x0)+T)) ! O6 |: d9 ]8 d( V: m( K" ]elseif p>0.85 & c<0.5     disp('The model is eligibility,and the forecast is:'), ! R. v7 p0 @/ L# }    disp(Hatx0(length(x0)+T)) elseif 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 : V# ]1 e6 S, W- W7 {; ~0 [    disp('The model is bad and try again') end for 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') endT=input('T=');%从键盘输入从最后一个历史数据算起的第T时点????是指什么啊，请大哥们，大姐们教一下，我急用，请快，谢谢我的初始值x0=[1.620938526 6 H5 Y, U* W' w8 C6 ?0 B0.07925621 - d3 N. E; N/ r5 T" D4 k, Q2 W# s0.052318818 ( m( r/ g; w% ?0.041252502 & x$ d( q& m5 z  {5 y3 V) K: a4 y0.021800479 ' I- |! Z$ p1 W0.053132975 ! z: {# G$ r- a9 M0.089908836 1 T/ G6 J5 k- O* g2 ?2 G; y/ r0.109153219 0.079331832 0 G- }8 ]8 G" ~# u0.342192598 5 q: @# R! O$ Z0.099718142 + h  f$ C/ Z  w0.135194823 ( {0 N1 m! G2 U& g5 {; t& `0.109274037 8 }" n) k# m5 ?4 n8 F3 a* M2 t0.08152013 5 o3 G: _% b7 |* `' n7 H0.067876355 0.064706843 0.055562197 0 b, e: Q, \2 J; {0.050848544 ]';

jshzncd

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

一剑卡卡

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

gxj820

很好啊！！顶顶！

杨晓敬

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

yysclshi

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

yysclshi

a= -0.0080
u= 0.0713
预测值
& c% m4 }$ w# E+ ?% Q    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
初始值
! \5 ^' f8 O& A: ]3 M    1.6209    0.0793    0.0523    0.0413    0.0218    0.0531
# @8 \4 w! i' k: ]( j6 h1 q    0.0899    0.1092    0.0793    0.3422    0.0997    0.1352
    0.1093    0.0815    0.0679    0.0647    0.0556    0.0508
残差
         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
相对误差
( M+ [# f- n9 k# b7 f! X) K! X         0    0.0672    0.6297    1.0835    2.9741    0.6437
, U2 \3 b! ^# X: h$ a& y    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 =
" m) W. r  x$ ]- E* z/ V# ?* w     1

daiqiang5566

楼主莫急..