快捷導(dǎo)航

Matlab實(shí)現(xiàn)灰色預(yù)測(cè)的示例代碼

更新時(shí)間：2022年05月17日 16:10:32 作者：slandarer

這篇文章主要為大家詳細(xì)介紹了如何利用Matlab實(shí)現(xiàn)灰色預(yù)測(cè)，文中的示例代碼講解詳細(xì)，具有一定的參考價(jià)值，感興趣的小伙伴可以了解一下

模型介紹

略微帶過(guò)一下原理：

灰色預(yù)測(cè)對(duì)于趨勢(shì)不強(qiáng)的數(shù)據(jù)，將其原始數(shù)據(jù)進(jìn)行累加后得到具有明顯趨勢(shì)的新數(shù)據(jù)進(jìn)行擬合，假設(shè)原數(shù)據(jù)為：

則新數(shù)據(jù)集X⁽¹⁾中數(shù)據(jù)為：

通常認(rèn)為累加數(shù)據(jù)服從指數(shù)分布，欸那我們直接假設(shè)：

直接進(jìn)行一個(gè)非線性擬合不就完事了，但是，從小學(xué)二年級(jí)開(kāi)始，老師就教導(dǎo)我們，像這用指數(shù)函數(shù)擬合啊，各個(gè)參數(shù)的變化對(duì)于整體曲線的影響效果差別很大啊，直接擬合誤差會(huì)很大啊，怎么能給他整成一個(gè)線性擬合??？(直接擬合我有試過(guò)，能夠大體描述趨勢(shì)并做出預(yù)測(cè)，但是誤差會(huì)比灰色預(yù)測(cè)大一點(diǎn))。

指數(shù)函數(shù)一般符合一個(gè)微分方程：

解常微分方程易得：

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

當(dāng)k≤t≤k+1時(shí)，有：

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

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

基礎(chǔ)代碼

核心代碼非常短，只有六行，代碼中給出了詳細(xì)的注釋?zhuān)A(chǔ)繪圖還是MATLAB風(fēng)比較簡(jiǎn)陋，后面會(huì)給出圖像修飾代碼。

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

X1=cumsum(X0);                    % 累加生成趨勢(shì)明顯新序列
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逐項(xiàng)做差求預(yù)測(cè)值

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

基礎(chǔ)代碼+修飾

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

X1=cumsum(X0);                    % 累加生成趨勢(shì)明顯新序列
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逐項(xiàng)做差求預(yù)測(cè)值

% 繪圖
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;

% 坐標(biāo)區(qū)域修飾
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;