getab.m

这是一个有预测功能的程序

function [coe]=getab(x,x1)
%得到系数a，b
n=length(x);
for i=1:n-1
  b(i)=-(x1(i)+x1(i+1))/2;
  yn(i)=x(i+1);
end
B=[b;ones(size(b))]';
coe=pinv(B)*yn';