灰色预测工具箱matlab,Matlab灰色预测工具箱——走过数模

2009-07-02 23:05

灰色预测几乎是每年数模培训必不可少的内容，相对来说也是比较简单，这里写了四个函数，方便在Matlab里面调用，分别是GM(1，1)，残差GM(1，1)，新陈代谢GM(1，1)，Verhust

自己写得难免有所疏忽，需要的朋友自己找本书本来试验一下。。

Gm(1，1)

function [px0,ab,rel]=gm11(x0,number)

%[px0,ab,rel]=gm11(x0,number)

%px0为预测数列，rel为平均相对误差,rel为平均相对误差(为百分比)

%默认的number参数为原数组大小

if nargin==1

number=max(size(x0));

end

n=max(size(x0));

x1=zeros(size(x0));

for k=1:n

for i=1:k

x1(k)=x1(k)+x0(i);

end

end

z=zeros(size(x0));

for k=2:n

z(k)=0.5*(x1(k)+x1(k-1));

end

y=x0";

y(1)=[];

b(:,1)=-z";

b(:,2)=1;

b(1,:)=[];

ab=inv(b"*b)*b"*y;