DGM.rar_DGM_DGM模型_raysrfx_分数阶累加_灰色模型

function [error]=cap_dgm_ga(r) global y0; x0=y0; n=length(x0); t=1:n; xr=leijia_ga(x0,r); Y=[xr(2:n)']; %x为一次累加后数据 A=ones(1,n-1); B=[xr(1:n-1)',A']; value=(pinv(B'*B))*B'*Y; %value=（ -0.1617 6.2706） b1=value(1); b2=value(2); x1_y(1)=x0(1); %digits(20);%规定了运算精度是6位有效数 for k=2:1:n x1_y(k)=(x0(1)-b2/(1-b1))*b1^(k-1)+b2/(1-b1); end x0_y=leijian_ga(x1_y,r); for m=2:n err(m)=abs(x0_y(m)-x0(m))/x0(m); end error=sum(err)/(n-1);