ercheng.m

来自「数值分析常用几种方法」· M 代码 · 共 15 行

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