第1个回答 推荐于2018-04-11
%灰色预测模型程序
clear
syms ab;
c=[a b]';
A=[27.4 29.061 31.851 34.223 40.521 46.54 47.065 52.085];%这里填已知的前1900个数(前面已填写的是8个我之前做的,直接删掉填你的就行了)
B=cumsum(A);%累加
n=length(A);
for i=1:(n-1)
C(i)=(B(i)+B(i+1))/2;
end
%计算待定参数
D=A;
D(1)=[];
D=D';
E=[-C; ones(1,n-1)];
c=inv(E*E')*E*D;
c=c';
a=c(1);
b=c(2);
%预测往后预测100个数据
F=[];F(1)=A(1);
for i=2:(n+100)
F(i)=(A(1)-b/a)/exp(a*(i-1))+b/a;
end
G=[];G(1)=A(1);
for i=2:(n+100)
G(i)=F(i)-F(i-1);
end
t1=1:1900;% 原始数据长度
t2=1:2000;%预测后的总长度
G
plot(t1,A,'o',t2,G)本回答被提问者和网友采纳