灰色GM(1,1)模型预测代码

预测未来10的代码：

N=size(X);m1=N(1);m2=N(2);

for i=1:m1

n=i

x0=X(i,:)

E=triu(ones(m2));x1=x0*E

b1=x1;b1(1)=[];b2=x1;b2(m2)=[];b=-0.5*(b1+b2);

B=[b;ones(1,m2-1)];B=B';

y0=x0;y0(1)=[];y0=y0';

A=((inv(B'*B))*B')*y0;

a=A(1),u=A(2),u_a=u/a

k=0:m2+9;x2=(x0(1)-u_a)*exp(-k*a)+u_a;

x3=x2;x3(m2+10)=[];x4=[0 x3];

x5=x2-x4

x6=x5(1:m2);Q=x0-x6;r=Q./x0;s1=std(x0);Qmean=mean(Q);s2=std(Q);C=s2/s1

D=abs(Q-Qmean);p0=0.6745*s1

t=0;

for j=1:m2

d=D(j);

if(dt=t+1;

end

end

P=t/m2

figure(i)

t1=1:m2;

t2=1:m2+10;

plot(t1,x0,'^r',t2,x5,'*-b');

end

预测未来3的代码：

N=size(X);m1=N(1);m2=N(2);

for i=1:m1

n=i

x0=X(i,:)

E=triu(ones(m2));x1=x0*E

b1=x1;b1(1)=[];b2=x1;b2(m2)=[];b=-0.5*(b1+b2);

B=[b;ones(1,m2-1)];B=B';

y0=x0;y0(1)=[];y0=y0';

A=((inv(B'*B))*B')*y0;

a=A(1),u=A(2),u_a=u/a

k=0:m2+2;x2=(x0(1)-u_a)*exp(-k*a)+u_a;

x3=x2;x3(m2+3)=[];x4=[0 x3];

x5=x2-x4

x6=x5(1:m2);Q=x0-x6;r=Q./x0;s1=std(x0);Qmean=mean(Q);s2=std(Q);C=s2/s1

D=abs(Q-Qmean);p0=0.6745*s1

t=0;

for j=1:m2

d=D(j);

if(dt=t+1;

end

end

P=t/m2

figure(i)

t1=1:m2;

t2=1:m2+3;

plot(t1,x0,'^r',t2,x5,'*-b');

end