zuixiaoerchengfa.m

来自「牛顿迭代法」· M 代码 · 共 23 行

function S=zuixiaoerchengfa(x,y,n,w)
% 数据的最小二乘拟合，其中
%   x,y --- 数据的(x,y)坐标
%   n  ---  数据拟合的次数，缺省时n=1
%   w  ---  权值，缺省时w=1
%   S  ---  数据拟合的系数
% 需要另写的函数
%   phi(x) --- 基函数，通常为多项式1, x, x^2, ...

global i; global j;
if nargin <4  w=1; end
if nargin <3  n=1; end
Phi2=zeros(n+1);
for i=0:n
    for j=0:n        
        Phi2(i+1,j+1)=sum((w.*phi_k(x,i)).*phi_k(x,j));
    end
end
PhiF=zeros(n+1,1);
for i=0:n
    PhiF(i+1)=sum((w.*phi_k(x,i)).*y);
end
S=Phi2\PhiF;