zxec.m

来自「数值线性代数的Matlab应用程序包 共13个程序函数」· M 代码 · 共 21 行

function c=ZXEC(x,y,n)

%x,y原始数据点
%n多项式次数

X=zeros(n+1,n+1);
U=zeros(n+1,1);
m=length(x);
for i=1:n+1
    for j=1:n+1
        for k=1:m
            X(i,n+2-j)=X(i,n+2-j)+x(k).^(i+j-2);
        end
    end
end
for i=1:n+1
    for k=1:m
        U(i)=U(i)+y(k)*x(k).^(i-1);
    end
end
c=X\U;