基于MATLAB实现最小二乘法的曲线拟合

function [a] = fitting(t,x,order)
k=zeros(length(x),order+1);
for i=1:length(x)
    for j=1:order+1
        k(i,j)=t(i)^(j-1);
    end
end
y=k'*x';
k=k'*k;
a=sorm(k,y,1.8);%超松弛迭代法求解正规方程
a=fliplr(a);
end

输出的a向量即为拟合方程曲线的多项式系数向量