最小二乘算法Matlab程序算法

1．最小二乘算法Matlab程序算法：
(1) 输入数据和权系数选择函数逼近系。
(2) 求解法方程组。
(3)求出最小二乘逼近。
数据的最小二乘拟和假定某天的气温变化记录如表所示 （表没了）

function S=square_least(x,y,n,w)
%数据的最小二乘拟合,其中
%x,y为数据的(x,y)坐标;
%n为数据拟合的次数，缺省时n=1;
%w为权值，缺省值w=1;
%S为数据拟合的系数;
%需要另写的外部函数.
%phi(x)为基函数,通常为多项式1,x,x^2,...
    global;global;
if nargin < 4w = 1; end
if nargin < 3n = 1; end
phi2=zeros(n+1);
for i=0:n
    for i=0:n
        phi2(i+1, j+1)=sum(w.* phi_k(x,i)).*phi_k(x,j);
    end
end
phiF= zeros(n+1,1);
end
S=phi\phiF