灰色预测的MATLAB代码

 数学建模可以直接带数据使用的灰色预测的MATLAB代码

clc,clear,close all %1.色预测程序
GDP=[108.3;109.08;110.03;110.08;111.31;112.68;114.16;109.63;109.11]';%需替换的数据
t=[1; 2; 3;4; 5;6;7;8;9]';
workers=[5441.4;5520.7;5748.6;6622.1;6896;7358.4;7875.2;8837.6;8831.2]';
n=length(GDP);
lamda=GDP(1:n-1)./GDP(2:n);
range=minmax(lamda) ;%最小最大值
x1=cumsum(GDP);%求和
for i=2:n
z(i)=0.5*(x1(i)+x1(i-1));%均值
end
B=[-z(2:n)',ones(n-1,1)];
 Y=GDP(2:n)';
u=B\Y
x=dsolve('Dx+a*x=b','x(0)=x0');
x=subs(x,{'a','b','x0'},{u(1),u(2),x1(1)});
yuce1=subs(x,'t',[0:n-1]);
digits(6),y=vpa(x)
%为提高预测精度，先计算预测值，再显示微分方程的解
yuce=diff(yuce1);
yuce=[x1(1),yuce]
epsilon=GDP-yuce%计算相对误差
delta=abs(epsilon./GDP)%计算残差
rho=1-(1-0.5*u (1))/(1+0.5*u(1))*lamda %计算级比偏差值

运行的结果：
u = -0.0013

110.0673

y =87600.2*exp(0.00125803*t) - 87491.9
yuce =[ 1083/10, 0, 0, 0, 0, 0, 0, 0, 0, 0]

级比偏差值：
rho = 0.0059    0.0074   -0.0008    0.0098    0.0109    0.0117   -0.0426   -0.0060