getab.m

来自「这是一个有预测功能的程序」· M 代码 · 共 9 行

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';