Matlab怎么實現灰色預測

這篇文章主要介紹“Matlab怎么實現灰色預測”的相關知識，小編通過實際案例向大家展示操作過程，操作方法簡單快捷，實用性強，希望這篇“Matlab怎么實現灰色預測”文章能幫助大家解決問題。

模型介紹

略微帶過一下原理：

灰色預測對于趨勢不強的數據，將其原始數據進行累加后得到具有明顯趨勢的新數據進行擬合，假設原數據為： 

則新數據集X⁽¹⁾中數據為：

通常認為累加數據服從指數分布，欸那我們直接假設：

直接進行一個非線性擬合不就完事了，但是，從小學二年級開始，老師就教導我們，像這用指數函數擬合啊，各個參數的變化對于整體曲線的影響效果差別很大啊，直接擬合誤差會很大啊，怎么能給他整成一個線性擬合??？(直接擬合我有試過，能夠大體描述趨勢并做出預測，但是誤差會比灰色預測大一點)。

指數函數一般符合一個微分方程： 

解常微分方程易得： 

這樣我們只需要將a,u這倆常數求出來就能得到x⁽¹⁾(t) ,就能得到X⁽¹⁾序列然后逐項做差就能得到X⁽⁰⁾序列。

當k≤t≤k+1時，有：

帶入回之前的微分方程就有： 

實際上直接進行線性擬合就好了，當然我們也可以拿最小二乘法裝模作樣的分析一下,令: 

基礎代碼

核心代碼非常短，只有六行，代碼中給出了詳細的注釋，基礎繪圖還是MATLAB風比較簡陋，后面會給出圖像修飾代碼。

X0=[15 16.1 17.3 18.4 18.7 19.1 19.9 21.3 22.5];
t=1:9;    % 原始數據自變量范圍
pt=10:12; % 預測數據自變量范圍

X1=cumsum(X0);                    % 累加生成趨勢明顯新序列
Z=X1(1:end-1)+diff(X1)./2;        % 均值，即(X1(1:end-1)+X1(2:end))./2
a_u=polyfit(-Z,X0(2:end),1);      % a_u=(B'*B)\B'*Y,B=[-Z,ones]
a=a_u(1);u=a_u(2);
P=(X1(1)-u/a)./exp(a.*([t,pt]-1))+u/a; % 求X1擬合值
P=[P(1),diff(P)];                      % X1逐項做差求預測值

% 繪圖
plot([t,pt],P,'*-');
hold on
plot(t,X0,'s-')  
legend('預測值','真實值');

基礎代碼+修飾

X0=[15 16.1 17.3 18.4 18.7 19.1 19.9 21.3 22.5];
t=1:9;    % 原始數據自變量范圍
pt=10:12; % 預測數據自變量范圍

X1=cumsum(X0);                    % 累加生成趨勢明顯新序列
Z=X1(1:end-1)+diff(X1)./2;        % 均值，即(X1(1:end-1)+X1(2:end))./2
a_u=polyfit(-Z,X0(2:end),1);      % a_u=(B'*B)\B'*Y,B=[-Z,ones]
a=a_u(1);u=a_u(2);
P=(X1(1)-u/a)./exp(a.*([t,pt]-1))+u/a; % 求X1擬合值
P=[P(1),diff(P)];                      % X1逐項做差求預測值

% 繪圖
plot([t,pt],P,'s-','Color',[82,124,179]./255,'MarkerFaceColor',[82,124,179]./255,...
    'MarkerEdgeColor',[1,1,1],'LineWidth',2,'MarkerSize',16)  
hold on
plot(t,X0,'d-.','Color',[169,64,71]./255,'MarkerFaceColor',[169,64,71]./255,...
    'MarkerEdgeColor',[1,1,1],'LineWidth',2,'MarkerSize',16)  

% 增添圖例
lgd=legend('fitting result','original data');
lgd.Location='best';
lgd.FontSize=16;

% 坐標區域修飾
ax=gca;grid on;box off
ax.LineWidth=2;
ax.Color=[249,250,245]./255;
ax.XMinorTick='on';
ax.YMinorTick='on';
ax.GridLineStyle='-.';
ax.XColor=[1,1,1].*.2;
ax.YColor=[1,1,1].*.2;
ax.FontName='Cambria';
ax.FontSize=14;

關于“Matlab怎么實現灰色預測”的內容就介紹到這里了，感謝大家的閱讀。如果想了解更多行業相關的知識，可以關注億速云行業資訊頻道，小編每天都會為大家更新不同的知識點。