最小二乘多项式拟合.m

来自「拉格朗日插值多项式拟合,牛顿插值多项式,欧拉方程解偏微分方程,使用极限微分求解导」· M 代码 · 共 25 行

%__________________________________
%最小二乘多项式拟合
%格式：Pm(X)=c1+c2X+c3X^2+......+cM+1X^M, 拟合N个数据点（Xk,Yk)(k=1...N)
%__________________________________
function C=lspoly(X,Y,M)
%X和Y为输入的两个对应的向量
%M为构造的M阶最小二乘多项式
n=length(X);
B=zeros(1:n);
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([1 2;3 4;5 6;7 8])
 %  ans =
 %  7     8
 %  5     6
 %  3     4
 %  1     2
C=flipud(C);