GM（11）模型代码

灰色理论的微分方程型模型称为GM模型，G表示grey(灰)，M表示Model(模型)．GM(1，N)表示1阶的，N个变量的微分 方程型模型．而GM(1，1)则是1阶的， 1个变量的微分方程型模型。

function []=greymodel（y）
% 本程序主要用来计算根据灰色理论建立的模型的预测值。
% 应用的数学模型是 GM（11）。
% 原始数据的处理方法是一次累加法。
y=input（‘请输入数据 ‘）;
n=length（y）;
yy=ones（n1）;
yy（1）=y（1）;
for i=2:n
    yy（i）=yy（i-1）+y（i）;
end
B=ones（n-12）;
for i=1:（n-1）
    B（i1）=-（yy（i）+yy（i+1））/2;
    B（i2）=1;
end
BT=B‘;
for j=1:n-1
    YN（j）=y（j+1）;
end
YN=YN‘;
A=inv（BT*B）*BT*YN;
a=A（1）;
u=A（2）;
t=u/a;
i=1:n+2;