GM（1，1）灰色預測模型的運算和檢驗代碼

function []=greymodel(y)
% 本程序主要用來計算根據灰色理論建立的模型的預測值。
% 應用的數學模型是 GM(1,1)。
% 原始數據的處理方法是一次累加法。
y=input(‘請輸入數據 ‘);
n=length(y);
yy=ones(n,1);
yy(1)=y(1);
for i=2:n
yy(i)=yy(i-1)+y(i);
end
B=ones(n-1,2);
for i=1:(n-1)
B(i,1)=-(yy(i)+yy(i+1))/2;
B(i,2)=1;
end
BT=B’;
for j=1:n-1
YN(j)=y(j+1);
end
YN=YN’;
A=inv(BT*B)BTYN;
a=A(1);
u=A(2);
t=u/a;
i=1:n+2;
yys(i+1)=(y(1)-t).*exp(-a.i)+t;
yys(1)=y(1);
for j=n+2:-12
ys(j)=yys(j)-yys(j-1);
end
x=1:n;
xs=2:n+2;
yn=ys(2:n+2);
plot(x,y,’^r’,xs,yn,’-b’);
det=0;

sum1=0;
sumpe=0;
for i=1:n
sumpe=sumpe+y(i);
end
pe=sumpe/n;
for i=1:n;
sum1=sum1+(y(i)-pe).^2;
end
s1=sqrt(sum1/n);
sumce=0;
for i=2:n
sumce=sumce+(y(i)-yn(i));
end
ce=sumce/(n-1);
sum2=0;
for i=2:n;
sum2=sum2+(y(i)-yn(i)-ce).^2;
end
s2=sqrt(sum2/(n-1));
c=(s2)/(s1);
disp([‘后驗差比值為：’,num2str?]);
if c<0.35
disp(‘系統預測精度好’)
else if c<0.5
disp(‘系統預測精度合格’)
else if c<0.65
disp(‘系統預測精度勉強’)
else
disp(‘系統預測精度不合格’)
end
end
end

disp(['下個擬合值為 ',num2str(ys(n+1))]);
disp([‘再下個擬合值為’,num2str(ys(n+2))]);

總結

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

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