x=load('x.txt'); %加载观测数据文件
n=length(x)
lamd=x(1:n-1)./x(2:n) %计算级比
exp(-2/(n+1))
exp(2/(n+1))
xx(1)=x(1:1);
length(x)
for k1=2:length(x)
xx(k1)=xx(k1-1)+x(k1); %运用 AGO 累加得到新的数列
end
xx
for k2=1:length(xx)-1
z(k2)=(xx(k2)+xx(k2+1))/2;
end
z
B=[-z',ones(length(xx)-1,1)] %构建矩阵 B
Y=x(2:end,:) %构建矩阵 Y
a=(inv(B'*B)*(B'*Y)); %求矩阵 a
a
for k3=1:length(x)
m1(k3)=(x(1)-(a(2)/a(1)))*(exp(-a(1)*(k3-1)))+(a(2)/a(1));%计算预测累加值
end
m1
m2(1)=m1(1); %对 m2(1)进行赋值