clear all;
x0=input('请输入原始数据列');
m=input('请输入预测数据个数');
n=length(x0);
for i=1:n;
x1(i)=sum(x0(1:i));
end
for k=2:n;
x2(k-1)=x0(k)-x0(k-1);
end
B=ones(n-1,2);
Y=ones(n-1,1);
for j=1:n-1;
B(j,1)=-x0(j+1);
Y(j,1)=x2(j);
end
zxec=inv(B'*B)*B'*Y;
a=zxec(1);
b=zxec(2);
syms k;
sjxys=(b/a^2-x0(1)/a)*exp(-a*k)+(b/a)*(k+1)+((1+a)/a)*(x0(1)-b/a);
for k=1:n-1;
moni(k+1)=(b/a^2-x0(1)/a)*exp((-a)*(k-1))*(exp(-a)-1)+b/a;
end
moni(1)=x0(1);
disp('模拟数据');
