%建立符號變量a（發展系數）和b（灰作用量）?
syms?a?b;?
c?=?[a?b]‘;?
%原始數列?A?
load?x0.mat
A?=?x0;?
n?=?length（A）;??
%對原始數列?A?做累加得到數列?B?
B?=?cumsum（A）;?
%對數列?B?做緊鄰均值生成?
for?i?=?2:n?
????C（i）?=?（B（i）?+?B（i?-?1））/2;?
end
C（1）?=?[];?
%構造數據矩陣?
B?=?[-C;ones（1n-1）];?
Y?=?A;?Y（1）?=?[];?Y?=?Y‘;?
%使用最小二乘法計算參數?a（發展系數）和b（灰作用量）
c?=?inv（B*B‘）*B*Y;
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
%對數列?F?累減還原得到預測出的數據?
G?=?[];?G（1）?=?A（1）;
for?i?=?2:（n+4）
????G（i）?=?F（i）?-?F（i-1）;%得到預測出來的數據?
end
disp（‘預測數據為：‘）;
G?
%模型檢驗?
H?=?G（1:104）;?
%計算殘差序列?
epsilon?=?A?-?H;?
%法一：相對殘差Q檢驗?
%計算相對誤差序列?
delta?=?abs（epsilon./A）;
%計算相對誤差Q?
disp（‘相對殘差Q檢驗：‘）
Q?=?mean（delta）
%法二：方差比C檢驗
disp（‘方差比C檢驗：‘）
C?=?std（epsilon?1）/std（A?1）

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????1262??2019-04-13?15:10??灰色預測GM（11）代碼\GM.m
?????文件?????????533??2019-04-13?09:07??灰色預測GM（11）代碼\x0.mat
?????目錄???????????0??2019-04-13?15:29??灰色預測GM（11）代碼\