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

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,