%灰色模型預測是在數據不呈現一定規律下可以采取的一種建模和預測方法，其預測數據與原始數據存在一定的規律相似性

%下面程序是灰色模型GM（11）程序二次擬合和等維新陳代謝改進預測程序matlab6.5?使用本程序請注明，程序存儲為gm1.m

%x?=?[59995903584857007884];gm1（x）;??測試數據?

%二次擬合預測GM（11）模型
function??gmcal=gm1（x）
sizexd2?=?size（x2）;
%求數組長度

k=0;
for?y1=x
????k=k+1;
????if?k>1
????????x1（k）=x1（k-1）+x（k）;
????????%累加生成
????????z1（k-1）=-0.5*（x1（k）+x1（k-1））;???
????????%z1維數減1，用于計算B
????????yn1（k-1）=x（k）;
????else
????????x1（k）=x（k）;
????end
end
%x1z1kyn1

sizez1=size（z12）;
%size（yn1）;
z2?=?z1‘;
z3?=?ones（1sizez1）‘;

YN?=?yn1‘;???%轉置
%YN

B=[z2?z3];
au0=inv（B‘*B）*B‘*YN;
au?=?au0‘;
%Bau0au

afor?=?au（1）;
ufor?=?au（2）;
ua?=?au（2）./au（1）;
%aforuforua?
%輸出預測的??a?u?和?u/a的值

constant1?=?x（1）-ua;
afor1?=?-afor;
x1t1?=?‘x1（t+1）‘;
estr?=?‘exp‘;