最小二乘拟合三次曲面

最新推荐文章于 2024-05-03 21:23:11 发布

Dana_ranmo

最新推荐文章于 2024-05-03 21:23:11 发布

阅读量752

点赞数 15

文章标签：线性代数算法最小二乘法

本文链接：https://blog.csdn.net/Dana_ranmo/article/details/135820891

版权

本文介绍了如何在最小二乘法的基础上扩展至三次曲面拟合，对比了与二次曲面的区别，强调了额外系数的处理，并提及了使用OpenCV库的实际应用，通过线性方程组求解10个系数的三次曲面方程。

摘要由CSDN通过智能技术生成

最小二乘法实现三次曲面拟合

最小二乘拟合二次曲面的方法有很多博主有创作，比如
https://blog.csdn.net/weixin_43319685/article/details/103167458
我在此基础上推导最小二乘法拟合三次曲面的方法，大体思路差不多，如有问题请指正，使用Opencv4 cpp实现的代码我将放在另一篇博客中。
假设三次曲面方程为：
在这里插入图片描述相比二次曲面方程，三次曲面方程多了4个系数，从6个系数增加到10个系数，采用最小二乘法求解，使得下面的式子误差最小：
$\sigma {2}=\sum_{i=1}{n}[I(x_{i},y_{i})-z_{i}]{2}=\sum_{i=1}{n}[(ax_{i}{3}+by_{i}{3}+cx_{i}{2}y_{i}+dx_{i}y_{i}{2}+ex_{i}{2}+fy_{i}{2}+gx_{i}y_{i}+hx_{i}+iy_{i}+j)-z_{i}]^{2}$
应满足：
$\frac{\partial \sigma ^{2}}{\partial a}=\frac{\partial \sum_{i=1}{n}[I(x_{i},y_{i})-z_{i}]{2}}{\partial a}=\sum_{i=1}{n}2x_{i}{3}[(ax_{i}{3}+by_{i}{3}+cx_{i}{2}y_{i}+dx_{i}y_{i}{2}+ex_{i}{2}+fy_{i}{2}+gx_{i}y_{i}+hx_{i}+iy_{i}+j)-z_{i}]=0$
在这里插入图片描述

写成矩阵形式

其中各项为
$A=\sum_{}^{}\begin{bmatrix} x{6}&x{3}y{3}&x{5}y&x{4}y{2}&x{5}&x{3}y{2}&x{4}y&x{4}&x{3}y&x^{3}\ x{3}y{3} & y^{6} & x{2}y{4} & xy^{5} & x{2}y{3} & y^{5} & xy^{4} & xy^{3} & y^{4} & y^{3}\ x^{5}y & x{2}y{4} & x{4}y{2} & x{3}y{3} & x^{4}y & x{2}y{3} & x{3}y{2} & x^{3}y & x{2}y{2} & x^{2}y \ x{4}y{2} & xy^{5} & x{3}y{3} & x{2}y{4} & x{3}y{2} & xy^{4} & x{2}y{3} & x{2}y{2} & xy^{3} & xy^{2}\ x^{5} & x{2}y{3} & x^{4}y & x{3}y{2} & x^{4} & x{2}y{2} & x^{3}y & x^{3} & x^{2}y & x^{2}\ x{3}y{2} & y^{5} & x{2}y{3} & xy^{4} & x{2}y{2} & y^{4} & xy^{3} & xy^{2} & y^{3} & y^{2}\ x^{4}y & xy^{4} & x{3}y{2} & x{2}y{3} & x^{3}y & xy^{3} & x{2}y{2} & x^{2}y & xy^{2} & xy\ x^{4} & xy^{3} & x^{3}y & x{2}y{2} & x^{3} & xy^{2} & x^{2}y & x^{2} & xy & x\ x^{3}y & y^{4} & x{2}y{2} & xy^{3} & x^{2}y & y^{3} & xy^{2} & xy & y^{2} & y\ x^{3} & y^{3} & x^{2}y & xy^{2} & x^{2} & y^{2} & xy & x & y & 1 \end{bmatrix}$