最小二乘法（多项式为基地）MATLAB实现

function A= zxecf(x,y,n)
% x,y为对应的插值点，为行向量；n为多项式阶数
n=n+1;
G=ones(n,n);
B=ones(n,1);
for i=1:n
    a=x.^(i-1);
    for j=1:n
        b=x.^(j-1);
        G(i,j)=sum(a.*b);
    end
    B(i,1)=sum(y.*a);
end
A=(G^(-1))*B;
end