灰色模型的MATLAB代碼——數學建模必備

先簡單介紹一下灰色模型（grey models）：

灰色模型（grey models）是通過少量的、不完全的信息，建立灰色微分預測模型，對事物發展規律作出模糊性的長期描述。

常見適用範圍：

描述系統具有層次、結構關係的模糊性

描述對象具有動態變化的隨機性

指標數據具有不完備或不確定性

前段時間參加了學校的數模建模大賽，第一次了解到灰色模型，在預測指標或者評估發展趨勢的用處還是挺大的

話不多說，源代碼如下，需要的同學複製粘貼運行就好了（會彈出預測圖的）~

clc,clear; syms a b; c=[a b]; A=[99215,109655,120333,135823,159878,182321,209407];%輸入你的原始數據值，我這裡用了7個B=cumsum(A); %對原始數據進行累加n=length(A); for i=1:(n-1) C(i)=(B(i)+B(i+1))/2; %生成累加矩陣end %計算待定參數的值D=A;D(1)=[]; D=D; E=[-C;ones(1,n-1)]; c=inv(E*E)*E*D; c=c; a=c(1);b=c(2); %對後續數據進行預測 F=[];F(1)=A(1); for i=2:(n+4) F(i)=(A(1)-b/a)/exp(a*(i-1))+b/a; end G=[];G(1)=A(1); for i=2:(n+5) G(i)=F(i)-F(i-1); %n加的數是指預測的數據個數，可以修改，我這裡預測了五個數據end t1=2011:2017; %假設數據按年份的時間軸變化，一開始用了七個數據，所以這裡從2011-2017，共七年t2=2011:2022; %原始的七個數據加上預測的五個數據共12個，故從2011-2022G plot(t1,A,o,t2,G,-); %這是繪圖函數grid onxlabel(年份) %橫軸名稱ylabel(預測值) %縱軸名稱title(基於2011-2017年的數據對未來五年的預測圖) %預測圖名稱

推薦閱讀：