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

根据引用和引用提供的信息，GM(2,1)的matlab代码如下：

clc,clear
x0=[41,49,61,78,96,104];
n=length(x0);
x1=cumsum(x0);
a_x0=diff(x0);
a_x0=[0,a_x0];
for i=2:n
    z(i)=0.5*(x1(i)+x1(i-1));
end
B=[-x0(2:end)',-z(2:end)',ones(n-1,1)];
Y=a_x0(2:end)';
u=B\Y;
x=dsolve('D2x+a1*Dx+a2*x=b','x(0)=c1,x(5)=c2');
x=subs(x,{'a1','a2','b','c1','c2'},{u(1),u(2),u(3),x1(1),x1(6)});
yuce=subs(x,'t',0:n-1);
digits(6)
x=vpa(x);
x0_hat=[yuce(1),diff(yuce(1:n-1))];
epsilon=x0-x0_hat;
delta=abs(epsilon./x0);

该代码实现了GM(2,1)模型的计算和预测，其中给定的输入数据为x0=[41,49,61,78,96,104]。通过GM(2,1)模型的计算，可以得到预测结果x0_hat，以及预测误差epsilon和误差率delta。