其中
Ra=y
r1 ( x1 ) rm ( x1 )
R
,
r1 ( xn ) rm ( xn )
a1
a
,
am
(3)
y1
y
yn
定理: 当RTR可逆时,超定方程组(3)存在最小二乘解, 且即为方程组
RTRa=RTy
解: a=(RTR)-1RTy
第13页
设 y*=a+bxi , 令 δi=yi-y*i=yi-a-bxi, 依据最 小二乘原理, 即使误差平方和达到最小, 也就是令
第22页
用MATLAB作非线性最小二乘拟合 Matlab(优化工具箱)提供了两个求非线性最小二 乘拟合函数: lsqcurvefit和lsqnonlin。 两个命令都要先建立M-文献定义函数, 然后编写 主程序M-文献求解非线性模型未知参数
注意: 二者定义f(x)方式是不同.
第23页
1. lsqcurvefit 已知数据点: xdata=(xdata1,xdata2,…,xdatan),
1.线性最小二乘拟合 2.非线性最小二乘拟合 3.cftool(自学)
ployfit \
lsqcurvefit lsqnonlin
第18页
多项式在x处值y可用下列命令计算: y=polyval(a,x)
用MATLAB作线性最小二乘拟合
1.作多项式f(x)=a1xm+ …+amx+am+1拟合,可利用已有
第28页
课堂小结
x0赋值:注意单位/依据结果重新赋值 ployfit(x,y,m)中x,y赋值: y换行要加… 画图: 非线性模型画出直线
非线性线性化后没有还原
第29页