求下面这个改进的灰色预测模型的matlab实现，用到了kolmogorov方程，做得好可加价

Question

求下面这个改进的灰色预测模型的matlab实现，用到了kolmogorov方程，残差修正预测模型，做得好可加价

Answer 1

可以借鉴下‘

function [XR]=GM1_1(x0)
[a,b]=size(x0);
x1=zeros(a,b);
z1=zeros(a,b-1);
for i =1:b %建立数列x1
    for j = 1:i
        x1(1,i)=x1(1,i)+x0(1,j);
    end
end
for i=2:b
    z1(1,i-1)=(x1(1,i)+x1(1,i-1))/2;
end
B=[-1.*z1',ones(b-1,1)];
Y=x0(1,2:b)';
U=(B.'*B)^(-1)*B.'*Y;
p=U(1,1);
q=U(2,1);
m=q/p;
XR=zeros(a,b);
for i =1:b+1
    XR(1,i)=m+(x0(1,1)-m)*exp(-p*(i-1));
end
for i=1:b
    XR(1,b+2-i)=XR(1,b+2-i)-XR(1,b+1-i);
end
end