最小二乘法的矩阵化表述

对于矩阵方程: A x = [ ( row ⁡ 1 ) ( row ⁡ 2 ) . . . ( row ⁡ m ) ] [ x 1 x 2 . . . x n ] = [ ( row ⁡ 1 ) ⋅ x ( row ⁡ 2 ) ⋅ x . . . . ( row ⁡ m ) ⋅ x ] = b A x=\left[\begin{array}{ll}(\operatorname{row} 1) \\ (\operatorname{row} 2) \\ ...\\(\operatorname{row} m)\end{array}\right]\left[\begin{array}{ll} x_{1} \\ x_{2} \\ ...\\ x_{n} \end{array}\right]=\left[\begin{array}{ll}(\operatorname{row} 1) \cdot x \\ (\operatorname{row} 2)\cdot x \\ ....\\(\operatorname{row} m) \cdot x\end{array}\right]=b Ax=(row1)(row2)...(rowm)x1x2...xn=(row1)x(row2)x....(rowm)x=b
当行数m>列数n时,方程可能无解。

而最小二乘法的目的就是让求取 x ^ \widehat{x} x 使得 E = ∥ A x − b ∥ 2 E=\|A x-b\|^{2} E=Axb2 最小。

根据矩阵的导数公式: ∂ ∂ X ( X b + c ) T D ( X b + c ) = ( D + D T ) ( X b + c ) b T \frac{\partial}{\partial \mathbf{X}}(\mathbf{X b}+\mathbf{c})^{T} \mathbf{D}(\mathbf{X b}+\mathbf{c})=\left(\mathbf{D}+\mathbf{D}^{T}\right)(\mathbf{X b}+\mathbf{c}) \mathbf{b}^{T} X(Xb+c)TD(Xb+c)=(D+DT)(Xb+c)bT

∂ E ∂ x = ∂ ∂ x ( A x − b ) T ( A x − b ) = 2 A T ( A x − b ) \frac{\partial E}{\partial x}=\frac{\partial}{\partial x}(A x-b)^{T}(A x-b)=2 A^{T}(A x-b) xE=x(Axb)T(Axb)=2AT(Axb)=0
因此: A T A x ^ = A T b A^{\mathrm{T}} A \hat{x}=A^{\mathrm{T}} b ATAx^=ATb
其中 x ^ \hat{x} x^ 是最小二乘解!

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

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值