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

xinzhiyong

GM(1,1)灰色模型的程序实现function GM1=fungry1(x0) %输入原始数据x0 T=input('T=');%从键盘输入从最后一个历史数据算起的第T时点 ( n1 @7 u4 Z9 x/ z1 qx1=zeros(1,length(x0));B=zeros(length(x0)-1,2); yn=zeros(length(x0)-1,1);Hatx0=zeros(1,length(x0)+T); / B6 P0 n. a% l8 k+ K, w1 fHatx00=zeros(1,length(x0));Hatx1=zeros(1,length(x0)+T); 3 e0 f- {: O4 i8 f5 B0 k/ `( |epsilon=zeros(length(x0),1);omega=zeros(length(x0),1); : T: I& ?7 ?7 P+ W. Y) Bfor i=1:length(x0)     for j=1:i         x1(i)=x1(i)+x0(j);     end end * c7 v$ q7 G/ ^6 [% M# Z% Mfor i=1:length(x0)-1 % m. A/ w8 S+ s4 V    B(i,1)=(-1/2)*(x1(i)+x1(i+1));   D0 r9 W  a1 z7 w3 S! z) |    B(i,2)=1; 2 N' x1 R( Z: u7 T5 S: D" {    yn(i)=x0(i+1); ! \8 e* a3 x9 r8 x" W( u3 Uend 2 i5 q& y8 r2 hHatA=(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 * z0 ~% A7 Z' A. C1 @* t    Hatx0(k)=Hatx1(k)-Hatx1(k-1);%累计还原得到历史数据的模拟值 end 0 z( {+ t% u; e9 {for i=1:length(x0) %开始模型检验 ; i& A& e# D( m2 h/ C    epsilon(i)=x0(i)-Hatx0(i); 3 S( D+ D' W% t6 T4 d3 n/ V& L; X    omega(i)=(epsilon(i)/x0(i))*100; ! g' \3 v/ k8 B7 J. L& _end % }8 t8 U; h+ r. j4 }% x0;Hatx0;epsilon;omega;  %必要时去掉％得到各种数据 c=std(epsilon)/std(x0);p=0; for i=1:length(x0) 3 L. O; |4 A, \% }    if abs(epsilon(i)-mean(epsilon))<0.6745*std(x0)         p=p+1; ( I0 c9 \. K! i1 e7 z; q    end - k: u' \: W- Q. l! |end / }1 n5 Y1 D7 ]; n1 ^& t- qp=p/length(x0) if p>0.95 & c<0.35     disp('The model is good,and the forecast is:'),     disp(Hatx0(length(x0)+T)) 6 y  ]* |& y$ S# \elseif p>0.85 & c<0.5     disp('The model is eligibility,and the forecast is:'),     disp(Hatx0(length(x0)+T)) 9 P- v- W- R0 ^elseif p>0.7 & c>0.65 8 ~% W% @, @! V" f9 A/ r    disp('The model is not good,and the forecast is:'),     disp(Hatx0(length(x0)+T)) % {. k  W! l- E% E% V. k2 G- V9 `else p<=0.7 & c>0.65     disp('The model is bad and try again') end & B$ `9 k: ~# X; L0 afor i=1:length(x0) % m- E+ m& v, ?! V, Z! P    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') ' h, B- w3 M) d  v8 z/ jendT=input('T=');%从键盘输入从最后一个历史数据算起的第T时点????是指什么啊，请大哥们，大姐们教一下，我急用，请快，谢谢我的初始值x0=[1.620938526 0.07925621 0.052318818 0.041252502 0.021800479 6 o4 Z! s* r* f4 h$ R  h% ^0.053132975 1 ~& m: w, l2 j/ ~& R  ~0.089908836 0.109153219 * f3 {( C! d* T0.079331832 ' L3 T  u; s4 b+ q0.342192598 0.099718142 0.135194823 0.109274037 0.08152013 6 L+ Z; l, V7 }, X. d0.067876355 2 D: Z2 {) e7 y: N0.064706843 6 C8 a, w1 j6 d( ]. T" w% Z( X6 g0.055562197 0.050848544 1 D5 ^* i& J% \( j  p/ ]]';

daiqiang5566

楼主莫急..

yysclshi

a= -0.0080
$ U$ g6 y: ^5 w- ^. Q& D$ T6 N# Qu= 0.0713
预测值
$ q7 \, j( w4 q7 x    1.6209    0.0846    0.0853    0.0859    0.0866    0.0873
    0.0880    0.0887    0.0895    0.0902    0.0909    0.0916
: n' [! o5 F! x! B: b4 J& D. \: f    0.0924    0.0931    0.0938    0.0946    0.0954    0.0961
初始值
2 ]; A6 I% S0 P- Z7 [1 N6 o    1.6209    0.0793    0.0523    0.0413    0.0218    0.0531
' D" t/ [2 a! h$ E8 F! K    0.0899    0.1092    0.0793    0.3422    0.0997    0.1352
  ?6 r( X7 M' e5 q3 l    0.1093    0.0815    0.0679    0.0647    0.0556    0.0508
残差
' S6 g% f+ c$ J& _: l         0   -0.0053   -0.0329   -0.0447   -0.0648   -0.0342
4 w, M4 d5 {2 X, o1 z    0.0019    0.0204   -0.0101    0.2520    0.0088    0.0436
/ L; s2 w5 `' ?: h    0.0169   -0.0116   -0.0260   -0.0299   -0.0398   -0.0453
8 `, ^% r: a$ s% `* Y相对误差
' _/ Y8 p) [3 q         0    0.0672    0.6297    1.0835    2.9741    0.6437
    0.0209    0.1870    0.1276    0.7365    0.0885    0.3223
: Q/ ~% z; Y# e1 T: y. M& I# v    0.1548    0.1420    0.3826    0.4619    0.7162    0.8903
方差比
    0.1869
- T! W/ W  {% E( j3 N) w6 S2 Ap =
     1

yysclshi

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

杨晓敬

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

gxj820

很好啊！！顶顶！

一剑卡卡

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

jshzncd

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