?數(shù)學建模可以直接帶數(shù)據(jù)使用的灰色預測的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]';%需替換的數(shù)據(jù) 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

總結

以上是生活随笔為你收集整理的灰色预测的MATLAB代码的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內容還不錯，歡迎將生活随笔推薦給好友。