灰色预测.txt

clear symsab; c=[ab]'; A=[489498546572606608]; B=cumsum(A);%yuanshishujuleijia n=length(A); fori=1:(n-1) C(i)=(B(i)+B(i+1))/2%shengchengleijiajuzhen end %jisuandaidingcanshudezhi 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); %yucehoux