GM(1,1)新陈代谢

function[caby]=GM_1_1oxcdx(x0,x)%师兄的原始GM新陈代谢x0=[4.445.898.068.749.6310.74];formatlongn=length(x0);m=length(x0);y=0;while(x=m)fork=1:nx1(k)=sum(x0(1:k));%生成x1序列end%x1;fork=2:nz1(k-1)=(x1(k)+x1(k-1))/2;%生成背景值z1序列endY=x0(2:n)';B=[-z1'ones(n-1,1)];N=inv(B'*B)*B'*Y;a=roundn(N(1,1),-6)b=roundn(N(2,1),-6)%a,b正确说明z1构造成功过b_chu_a=roundn(b/a,-6);fork=0:n-1X1(k+1)=b/a+(x0(1)-b_chu_a)*exp(-a*k);%改正x（k）的生成end%X1fork=1:n-1X0(k+1)=X1(k+1)-X1(k);%还原值x0（k）的生成endX0;y=b/a+(x0(1)-b_chu_a)*exp(-a*n);X1(n+1)=y;y-X1(n)%第一个预测值，不是模拟x0（n+1）的改正，下面要替换掉x0（1）fork=1:length(x0)-1x0(k)=x0(k+1);endx0(length(x0))=y-X1(n);x0m=m+1;endend