灰色预测GM11带后验差检验

灰色预测GM(1,1)代码，带后验差检验，可计算C和p。附带数据，简单易懂。

%建立符号变量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）代码\