灰色预测的matlab代码

% 灰色预测 GM（11）模型
% 预测精度等级 
% 好 P>0.95 C<0.35 
% 合格 P>0.80 C<0.45 
% 勉强合格 P>0.70 C<0.50 
% 不合格 P<=0.70 C>=0.65

clcclear;
syms a b;
c=[a b]‘;
A=[1240	1243	1850	2641	2692]; %输入数据;
B=cumsum（A）; %原始数据累加
n=length（A）;
for i=1:（n-1）
    C（i）=（B（i）+B（i+1））/2; %生成累加矩阵
end
%计算待定参数的值
D=A;D（1）=[];
D=D‘;
E=[-C;ones（1n-1）];
c=inv（E