最小二乘多项式拟合.m

拉格朗日插值多项式拟合,牛顿插值多项式,欧拉方程解偏微分方程,使用极限微分求解导数（微分）,微分方程组的N=4龙格库塔解法,雅可比爹迭代法解方程AX=B,最小二乘多项式拟合,组合辛普生公式求解积分,用

%__________________________________
%最小二乘多项式拟合
%格式：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);