最小二乘法矩阵求导过程的推导

关于最小二乘问题的求解,之前已有梯度下降法,还有比较快速的牛顿迭代。今天来介绍一种方法,是基于矩阵求导来计算的,它的计算方式更加简洁高效,不需要大量迭代,只需解一个正规方程组。在开始之前,首先来认识一个概念和一些用到的定理。矩阵的迹定义如下 

一个的矩阵的迹是指的主对角线上各元素的总和,记作。即      

 

 

                        

             

 好了,有了上述7个定理,就可以来求最小二乘解了。设

 

  

那么进一步得到 

    

 接下来会涉及到矩阵求导,因为 

    

那么进一步利用矩阵求导并利用上述定理,得到

 

    

 我们知道在极值点处梯度值为零,即 :

    

上述得到的方程组叫做正规方程组,那么最终得到 

    

这样最小二乘问题只需解一个线性方程组即可,不再需要像梯度下降那样迭代了。 

既然说到了正规方程组,在介绍一种方程组,叫做超定方程组,它的定义为:把方程个数大于未知量个数的方程组叫做超定方程组。通常来说,对于一个超定方程组来说,求最小二乘解只需要两边同时乘的转置,然后得到正规方程组,然后解这个方程就得到了最小二乘解。

  • 5
    点赞
  • 20
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
最小二乘法拟合圆是一种常用的曲线拟合方法,可以用于将一组点拟合成一个圆。其原理是将点到圆心的距离的平方作为误差,通过最小化误差来求解圆心坐标和半径。 以下是最小二乘法拟合圆的公式推导及其在C语言中的实现: 1. 公式推导 设圆的方程为:(x-a)^2 + (y-b)^2 = r^2 将每个点代入圆的方程中,得到如下方程组: (x1-a)^2 + (y1-b)^2 = r^2 (x2-a)^2 + (y2-b)^2 = r^2 ... (xn-a)^2 + (yn-b)^2 = r^2 将圆心坐标和半径表示为未知数,即: A = [a,b,r] 则上述方程组可以表示为: f(A) = [ (x1-a)^2 + (y1-b)^2 - r^2, (x2-a)^2 + (y2-b)^2 - r^2, ... (xn-a)^2 + (yn-b)^2 - r^2 ] 误差函数为: E(A) = f(A)^T * f(A) 对误差函数求导,得到: dE/dA = -2 * J^T * f(A) 其中,J为雅可比矩阵,表示误差函数对未知数的偏导数,即: J = [ -2*(x1-a), -2*(y1-b), -2*r; -2*(x2-a), -2*(y2-b), -2*r; ... -2*(xn-a), -2*(yn-b), -2*r ] 将雅可比矩阵和方程组代入上式,得到: dE/dA = -2 * ( J^T * J * A + J^T * f(A) ) 令dE/dA = 0,解得: A = ( J^T * J )^-1 * J^T * f(A) 即可求出圆心坐标和半径。 2. C语言实现 在C语言中,可以使用矩阵运算库(如Eigen)来实现矩阵的乘法、转置和求逆等操作。以下是一个简单的C语言实现: ```c #include <Eigen/Dense> #include <vector> using namespace Eigen; using namespace std; // 最小二乘法拟合圆 // 输入:points-点集,格式为vector<Vector2d> // 输出:圆心坐标和半径,格式为Vector3d Vector3d fitCircle(const vector<Vector2d>& points) { int n = points.size(); // 构建J和f(A)矩阵 MatrixXd J(n, 3); VectorXd f(n); for (int i = 0; i < n; i++) { double x = points[i](0); double y = points[i](1); J(i, 0) = -2 * x; J(i, 1) = -2 * y; J(i, 2) = -2; f(i) = -(x * x + y * y); } // 求解圆心坐标和半径 Vector3d A = (J.transpose() * J).inverse() * J.transpose() * f; A(0) = -A(0) / 2; A(1) = -A(1) / 2; A(2) = sqrt(A(0) * A(0) + A(1) * A(1) - A(2)); return A; } ``` 使用示例: ```c vector<Vector2d> points; points.push_back(Vector2d(0, 1)); points.push_back(Vector2d(1, 0)); points.push_back(Vector2d(0, -1)); points.push_back(Vector2d(-1, 0)); Vector3d A = fitCircle(points); cout << "圆心坐标:" << A.head<2>().transpose() << endl; cout << "半径:" << A(2) << endl; ```
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值