function?[z?xnew]=GM11（x0t）
%z=[a?u]?xnew為預(yù)測(cè)矩陣?t為預(yù)測(cè)的年份
n=length（x0）;
for?i=1:n
????x1（i）=sum（x0（1:i））;
end
Z（1）=0;
for?i=2:n
????Z（i）=0.5*（x1（i）+x1（i-1））;
end
A=[-Z（2:n）‘ones（n-11）];
b=x0（2:n）‘;
z=A\b;
xnew（1）=x0（1）;
a=z（1）;
u=z（2）;
for?i=2:n+t
????xnew（i）=（x0（1）-u/a）*（exp（-a*（i-1））-exp（-a*（i-2）））;
end

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件????????325??2009-08-13?17:44??GM11_2.m

?????文件????????335??2009-08-13?17:23??GM11_1.m

-----------?---------??----------?-----??----

??????????????????660????????????????????2