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

Matlab灰色预测模型GM(1,1)代码function c7fun73 X0=[2.874 3.278 3.307 3.39 3.679]; AU=c7fun73(X0);a=AU(1);u=AU

Matlab灰色预测模型GM(1,1)代码 function c7fun73 X0=[2.874 3.278 3.307 3.39 3.679]; AU=c7fun73(X0);a=AU(1);u=AU(2);m2=length(X0); for k=1:1:m2-1 xx1(k+1)=(X0(1)-u/a)*exp(-a*k)+u/a; end s=0;xx0(1)=X0(1); for jj=2:1:m2; xx0(jj)=xx1(jj)-xx1(jj-1); end disp('GM(1,1)对数列进行预测结果');xx0 disp('数列1原始观测数据');X0 disp('a');AU(1) disp('u');AU(2) function au=c7fun73(X0) m=length(X0); s1=0; for jj=1:1:m;X1(jj)=s1+X0(jj); s1=X1(jj); end for ii=1:1:m-1; B(ii)=-(X1(ii)+X1(ii+1))/2; end B=[B(:),ones(m-1,1)]; y=X0([2:m])'; au=inv((B'*B))*B'*y;