最小二乘椭球拟合 知乎_6.7 最小二乘拟合问题

本文探讨了最小二乘拟合问题,特别是在椭球拟合中的应用。介绍了数据拟合的一般形式,强调拟合函数不局限于多项式。通过设置残差平方和最小化,形成线性最小二乘问题,并利用偏导数为零找到极值点。讨论了点集函数的内积及其性质,以及如何将其转换为矩阵形式,即Gram矩阵。文章还提醒读者,仅函数系的线性无关性不足以保证Gram矩阵非奇异,需附加Haar条件,当函数系满足此条件时,Gram矩阵可逆。
摘要由CSDN通过智能技术生成

数据拟合问题的一般形式:

任给一组离散数据

equation?tex=%5C%7B%28x_%7Bi%7D%2Cy_%7Bi%7D%29%5C%7D%280%5Cleq+i%5Cleq+m%29 ,设拟合函数为:

equation?tex=P%28x%29%3Da_%7B0%7D%5Cvarphi_%7B0%7D%28x%29%2Ba_%7B1%7D%5Cvarphi_%7B1%7D%28x%29%2Ba_%7B2%7D%5Cvarphi_%7B2%7D%28x%29%2B...%2Ba_%7Bn%7D%5Cvarphi_%7Bn%7D%28x%29 ,其中
equation?tex=%5C%7B%5Cvarphi_%7Bi%7D%28x%29%5C%7D%280%5Cleq+i+%5Cleq+n%29 为给定的线性无关函数系,
equation?tex=a_%7B0%7D%2Ca_%7B1%7D%2C...a_%7Bn%7D 为要求得的未知参数。

(注:这里的拟合函数不一定为多项式函数)

记残量的平方和为:

equation?tex=r%28a_%7B0%7D%2Ca_%7B1%7D%2C...a_%7Bn%7D%29%3D%5Csum_%7Bj%3D0%7D%5E%7Bm%7D%7B%28P%28x_%7Bj%7D%29-y_%7Bj%7D%29%7D%5E%7B2%7D

求使得残量平方和最小得一组系数就是线性最小二乘问题,

equation?tex=%5C%7B%5Cvarphi_%7Bi%7D%28x%29%5C%7D%280%5Cleq+i+%5Cleq+n%29 为最小二乘问题得基函数,求得的拟合函数为数据的最小二乘拟合。

求解:利用偏导数为零得到极值点的原理可以得到最小二乘问题满足的方程组,求解方程组中未知系数(也就是拟合函数中的参数)即可。

点集函数的内积:

91b952da2140121059a21f46e7851693.png

内积的三条基本性质:1.非负性 2.对称性 3.线性性

通过引入点集函数的内积,便可将上面的方程组写成矩阵的形式:

65769abc73434e6e893268ce5e7a04a4.png

上面的系数矩阵被称为Gram矩阵

cec37576c1004ee7b7594872e6da70b2.png

558ca349023ae7cebfe9355a0fb3e25a.png

注意:仅凭函数系的线性无关性推不出Gram矩阵非奇异

如:

equation?tex=%5Cvarphi_%7B0%7D%3Dsinx%2C%5Cvarphi_%7B1%7D%3Dsin2x%2Csinx%E4%B8%8Esin2x%E5%9C%A8%5B0%2C2%5Cpi%5D%E4%B8%8A%E6%98%8E%E6%98%BE%E7%BA%BF%E6%80%A7%E6%97%A0%E5%85%B3%2C%E5%8F%96x_%7B0%7D%3D0%2C+x_%7B1%7D%3D%5Cpi%2Cx_%7B2%7D%3D2%5Cpi%5C%5C%28%5Cvarphi_%7B0%7D%2C%5Cvarphi_%7B0%7D%29%3D0%2C%28%5Cvarphi_%7B1%7D%2C%5Cvarphi_%7B1%7D%29%3D0%2C%28%5Cvarphi_%7B0%7D%2C%5Cvarphi_%7B1%7D%29%3D0%2C%E6%95%85G%3D0%EF%BC%8CGram%E7%9F%A9%E9%98%B5%E4%B8%BA%E4%B8%8D%E5%8F%AF%E9%80%86%E7%9F%A9%E9%98%B5%E3%80%82

故必须加上下面的Haar(哈尔)条件:

函数系

equation?tex=%5C%7B%5Cvarphi_%7Bi%7D%28x%29%5C%7D%280%5Cleq+i+%5Cleq+n%29 的任意线性组合在点集
equation?tex=%5C%7Bx_%7Bi%7D%5C%7D%280%5Cleq+i+%5Cleq+m%29%28m%5Cgeq+n%29 上至多有n个不同的零点,则称函数系在该点集上满足Harr条件。

定理:若函数系在该点集上满足Harr条件,则得到的Gram矩阵是可逆的。

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值