灰色预测gm(1,1)MATLAB代码

clear

%原始数据输入；

x0=[3.936;4.575;4.968;5.063;5.968;5.507];

n=length(x0);

%由数列x0生成数列x1;

s=0;

for i=1:n

s=s+x0(i);

x1(i)=s;

end

%由数列x1生成矩阵C和A；

for i=1 : (n-1);

A(i)=-(x1(i)+x1(i+1))/2;

end

A=[A' ones(n-1,1)];

for i=1 : (n-1);

C(i)=x0(i+1);

end

C=C';

beta=inv(A'*A)*A'*C;

%计算出a和u;

a=beta(1);

u=beta(2);

%求出数列x1的预测值数列x2;

m=input('请输入预测的数目:')+length(x0);

for k=0 : (m-1)

x2(k+1)=(x0(1)-u/a)*exp(-a*k)+u/a;

end

%求出原始数列x0的预测值数列x3;

x3(1)=x0(1);

for k=1 : (m-1)

x3(k+1)=(1-exp(a))*(x0(1)-u/a)*exp(-a*k); end

x3

t0=1:length(x0);

t3=1:m;

plot(t0,x0,'*',t3,x3)