灰色预测系统GM（1，1）模型程序

firmsoul

可用于一般时间序列的预测

1. function GM=huise(data,N)
   
2. T=length(data);
   
3. X0=data;
   
4. for i=2:T
   
5.    X1(1)=X0(1);
   
6.    X1(i)=X1(i-1)+X0(i);                %用AGO生成一阶累加生成模块
   
7. end
   
8. for i=1:T-1
   
9.    M(i)=-(0.5*(X1(i)+X1(i+1)));
   
10. end
    
11. B=zeros(T-1,2);                       %构造累加矩阵B
    
12. for i=1:T-1
    
13.      for j=1:2
    
14.          if j<2
    
15.             B(i,j)=M(i);
    
16.          elseif j>1
    
17.             B(i,j)=1;
    
18.          end
    
19.      end
    
20. end
    
21. for i=2:T                          %构造常数项向量Y
    
22.     Y(i-1)=X0(i);
    
23. end
    
24. HCS=inv(B'*B)*B'*Y'             %用最小二乘法求灰参数HCS
    
25. H=HCS';
    
26. for i=2:T+N                        %计算出累加序列
    
27. XR1(i)=(X0(1)-H(2)/H(1))*exp(-1*H(1)*(i-1))+H(2)/H(1);
    
28. end
    
29. for i=11:T+N                      %还原计算出预测值
    
30.         K(i-10)=XR1(i)-XR1(i-1);
    
31. end
    
32. GM=K;

复制代码

将函数保存为huise.m文件，然后在命令行中输入以下格式
输入格式：

1. data＝[x1,x1,...,xn]
   
2. huise(data,N)

复制代码

然后回车
哈哈，就可以了，没事干写来练练手的，残差分析部分就不写了，要是觉得精度不够就自己补上去：）写得不好，大家将就用吧

[ 本帖最后由 suffer 于 2006-10-9 20:08 编辑 ]

zsh1228

执行不了呀!哪个地方弄错了吧!
我是新手!!请多指点!!
\谢了!

caichengtao

楼主，你的K是做什么用的？？？

mcx8305

程序没经过调试请别拿出来！:@@

firmsoul

一年多前写的小程序了调试过了，只是输入格式搞错了。那天喝多了没注意
输入数据其实应该是data=[x1 x2 x3....xn]
是偶不对，不过逗号那么明显的错误，你们也不细心看看，汗

h_zf00

好象还是有错误

okor

写的还是不错的