函數pre = GM11(x)
%x，建模數據序列x = [x(1)，x(2)，...，x(n)]
% pre，GM11模型得到的預測值，即x估計的x(n+1)。
n =長度(x)；
x 1 = cumsum(x)；
z 1 =(x 1(1:n-1)+x 1(2:n))/2；
Y = x(2:長度(x))’；
B = [-z1 '，個位(n-1，1)]；
a =(B ' * B)\ B ' * Y；%GM(1，1)模型的A、B參數組成的向量。
beta = A(2)/(1+. 5 * A(1))；
alpha = A(1)/(1+. 5 * A(1))；
%初始序列x0(k)模擬序列隨時間t變化
f = @(t)(beta-(alpha)* x(1))* exp(-(A(1)*(t-2)))；
pre = f(n+1)；
轉移
x =[1879 2208 2948 2852 2774 2926 4680 3691]；
GM11(x)