【Matlab 041期】【數學建模】Matlab 電力預測預測之灰度預測組合預測指數平滑回歸分析

clc;
clear all
data=[43.83
46.58
48.29
51.77
53.30
57.22
58.94
63.37
65.36
70.48
72.57
78.30
86.61
87.42
96.37
97.38
107.46
100.80
109.88
118.03
114.39
118.17
139.01
186.48
198.45];
t=1995:2019;
T=1995:2025;
%% 灰度預測
 
%% 迴歸預測
 
disp(['迴歸預測均值方差=',num2str(wucha1)]);
disp(['迴歸預測百分絕對方差=',num2str(wucha2*100),'%']);
figure(2)
plot(t,data,'bo-',T,G_huigui,'r*-');  %原始資料與預測資料的比較
xlabel('年份')
ylabel('用電量');
legend('實際值','迴歸預測值')
title('迴歸預測用電量')
%% %% 二次指數平滑
 
 
S(3)=std(G_pinghua);%標準差
format short
[wucha1,wucha2]=mse(G_pinghua(21:25),data(21:25));
wucha1=roundn(wucha1,-2);
wucha2=roundn(wucha2,-4);
disp(['二次指數平滑預測均值方差=',num2str(wucha1)]);
disp(['二次指數平滑預測百分絕對方差=',num2str(wucha2*100),'%']);
figure(3)
plot(t,data,'bo-',T,G_pinghua,'r*-');  %原始資料與預測資料的比較
xlabel('年份')
ylabel('用電量');
legend('實際值','二次指數平滑預測值')
title('二次指數平滑預測用電量')
%% 組合預測
 
disp(['變權組合預測均值方差=',num2str(wucha1)]);
disp(['變權組合預測預測百分絕對方差=',num2str(wucha2*100),'%']);
figure(4)
plot(t,data,'bo-',T,G_zuhe,'r*-');  %原始資料與預測資料的比較
xlabel('年份')
ylabel('用電量');
legend('實際值','變權組合預測值')
title('變權組合預測用電量')

注：完整程式碼或者代寫新增QQ2449341593