nihe.m

来自「就是一件小小的源码」· M 代码 · 共 13 行

% 最小二乘法多项式拟合 
function p = nihe(x,y,m)    
% 用途：多项式拟合
% 格式：p=nihe(x,y,m),x,y是数据向量，m为拟合次数，p为返回多项式系数降幂排列
A=zeros(m+1,m+1);
for i =0 : m
    for j =0 : m
        A(i+1,j+1)=sum(x.^(i+j));
    end
    b(i+1)=sum(x.^i.*y);
end
a=A\b';
p=fliplr(a');