数学建模源码

灰色系统理论建模源代码 function GM1_1(X0) %format long  [m,n]=size(X0);   X1=cumsum(X0);   %累加  X2=[]; for i=1:n-1      X2(i,:)=X1(i)+X1(i+1); end  B=-0.5.*X2  t=ones(n-1,1);  B=[B,t]  ;      % 求B矩阵 YN=X0(2:end)  ;  P_t=YN./X1(1:(length(X0)-1)) %对原始数据序列X0进行准光滑性检验，                             %序列x0的光滑比P(t)=X0(t)/X1(t-1) A=inv(B."*B)*B."*YN."  a=A(1)  u=A(2)  c=u/a  ; b=X0(1)-c    X=[num2str(b),"exp","(",num2str(-a),"k",")",num2str(c)];  strcat("X(k+1)=",X)  %syms k;   for t=1:length(X0)      k(1,t)=t-1;  end  k  Y_k_1=b*exp(-a*k)+c; for j=1:length(k)-1     Y(1,j)=Y_k_1(j+1)-Y_k_1(j); end  XY=[Y_k_1(1),Y]    %预测值 CA=abs(XY-X0)     %残差数列  Theta=CA       %残差检验 绝对误差序列  XD_Theta= CA ./ X0   %残差检验 相对误差序列 AV=mean(CA); &