sublspoly.m

直线拟合与多项式拟合 涉及：《数值计算》等相关课程的基础内容 共两个例子

function [C,F] = sublspoly(X,Y,M)

% 调用形式：
%       [C,F] = sublspoly(X,Y,M)
% 功能：
%       用最小二乘法多项式拟合非线性数据
% 输入：
%       -- X   行向量，包含拟合数据点的横坐标
%       -- Y   行向量，包含拟合数据点的纵坐标
%       -- M   标量，给定的多项式的阶数
% 输出：
%       -- F N x (M+1)的矩阵，每一行表示一个数据的x值，同一行的元素从左到右分别表示x的0次方到M次方
%      --C  列向量，从低阶到高阶排列多项式的系数
% 袁怡圃 2003/04/30

n=length(X);
B=zeros(1:M+1);
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;