clear all

x=[0.8600    0.8244    0.8347    0.8451    0.8557    0.8664				
];
n=size(x,2);
%求級數比
for i=1:n-1
    y(i)=x(i)/x(i+1);
end
y；
%判斷級數是否范圍內

x
%求x1
sum=0;
a=0.5;
for i=1:n
    x1(i)=x(i)+sum;
    sum=sum+x(i);
end
x1;
%求z1
for i=1:n-1
    z1(i)=a*x1(i)+(1-a)*x1(i+1);
end
   z1
   %生成B
   B=[-z1',ones(n-1,1)]
   %生成Y1
   for i=1:n-1
       y1(i)=x(i+1);
   end
   Y1=y1';
   Y1
   
   x;
   %求(a,b)
   ab=inv(B'*B)*B'*Y1
   %求預測值
   for i=1:n
       yux(i)=(x(1)-ab(2)/ab(1))*exp(-ab(1)*i)+ab(2)/ab(1);
       yux1(i)=(x(1)-ab(2)/ab(1))*(exp(-ab(1)*i)-exp(-ab(1)*(i-1)))
   end
   yux
   x1
   x