灰色模型检验的MATLAB程序 （好像有三种检验方法：残差 关联度 方差比和小概率）

司马空觞

x=   ;%随便一列数
n=length(x);
x0=x(:,2);
for i=1:n
    x1(i)=sum(x0(1:i));
end
for k=2:n %紧邻均生成Z
  z(k)=0.5*x1(k)+0.5*x1(k-1);
end
for i=1:n-1
b(i,1)=-z(i+1);
Y(i)=x0(i+1);
end
b(:,2)=1;
au=inv(b'*b)*b'*Y';%计算au
yc1(1)=x0(1);
for k=1:n
c=x0(1)-au(2)/au(1);
yc1(k+1)=c*exp(-au(1)*k)+au(2)/au(1);
end
yc0(1)=x0(1);
for k=1:n
yc0(k+1)=yc1(k+1)-yc1(k);
end
disp(uint16(yc0(2:1:n+1)));

袁海亮

楼主 这是遇到什么事情了呢 ？

司马空觞

袁海亮 发表于 2013-8-27 07:48
楼主 这是遇到什么事情了呢 ？

跟自己过不去，做几个模型玩一下

magic2728

这，。。。