canchagm.m

用matlab实现灰色序列残差预测模型GM（1

e0=[31.0615   40.4060   77.1353   21.1143]; %原始序列
[m,n]=size(e0);
a_l=exp(-2/(n+1))            %计算级比覆盖范围
a_h=exp(2/(n+1))
e0=sqrt(e0);                 %对原始数据作两次方根变换   
e0=sqrt(e0);
for i=2:n
    q(i)=e0(i-1)/e0(i);      %计算实际级比
end
q
for i=1:n                     %计算一次累加序列
    e1(i)=sum(e0(1:i));       
end
for k=1:n-1                     
    z1(k)=1/2*e1(k)+1/2*e1(k+1);
end
Y=e0(2:n)';                    %得出Y矩阵
B=[];
B(:,1)=-z1';
B(:,2)=ones(1);                 %得出B矩阵
aU1=(inv(B'*B))*B'*Y;            %计算a,U 
a1=aU1(1);
U1=aU1(2);
for i=1:n 
    ee1(i)=(e0(1)-U1/a1)*exp(-a1*(i-1))+U1/a1; %解偏微分方程
end
for i=1:n-1
    ee0(i+1)=ee1(i+1)-ee1(i);           %一次累减，得出预测值
end
ee0=ee0.^4;                              %作四次方还原