数值计算方法-最小二乘拟合

要求：
编写一个函数，输入为节点矩阵，基函数，输出图像，显示节点和拟合曲线；
理论：
最小二乘是要求对于给定数据列 ，要求存在某个函数类

中寻求一个函数：
，使得 满足
根据以上条件可知，点

是多元函数

的极小点，从而 满足方程组

即 ，

记
，则上述方程组可表示成 ， (k=0,1, …， n)
写成矩阵形式为

，这个方程组成为法方程组，可以证明，当
线性无关时，它有唯一解。
特别地，曲线拟合的一种常用情况为代数多项式，即取

，则

(k=0,1, …， n)
故相应的法方程组变为

，这就是最小二乘法的原理。
Matlab代码：

①主函数
function u=OLSE_t(x,y,w,a)
syms b;
p=zeros(a,a);
f=zeros(a,1);
fai=zeros(a,length(x));
for m = 1:length(x)
    fai(1,m)=log(x(m));
    fai(2,m)=cos(x(m));
    fai(3,m)=exp(x(m)