lspoly.m

来自「余向阳博士Matlab教程的源代码。 包括：1.非线性方程的解法；2.线性方程」· M 代码 · 共 14 行

function c=lspoly(x,y,m)
%  x是数据点的横坐标组成的向量，y是纵坐标组成的向量。
%  m是要构造的多项式的次数，c是多项式由高次到低次的系数所组成的向量。
n=length(x);
b=zeros(1:m+1);
f=zeros(n,m+1);
for k=1:m+1
    f(:,k)=x'.^(k-1);
end
a=f'*f;
b=f'*y';
c=a\b;
c=flipud(c);