超定方程的最小二乘解的三维几何解释

  • 原始方程 Ax=bAx = b,解为 x=A1bx = A^{-1}b,matlab描述 x = A\b
  • 超定方程乘以 ATA^T 变为方阵 ATAx=ATbA^TAx = A^Tb
  • 列向量的形式 ATAA^TA 直接是一个数,简化计算
  • 再把 ATAA^TA 作为一个整体除过去
    x=(ATA)1ATbx = (A^TA)^{-1} A^Tb

  • 最小二乘解,向量 bb' 在张成的平面之外,解的满足误差最小,合成的向量是 bb' 在张成的投影
    AT(bAx)=0A^T(b'-Ax')=0
  • v 和 w 均为列向量
    A=[v,w]=[x1x2y1y2]A = [v,w]=\begin{bmatrix} x_1&x_2\\ y_1&y_2 \end{bmatrix}
    AT=[vTwT]=[x1y1x2y2]A^T = \begin{bmatrix} v^T\\ w^T \end{bmatrix}=\begin{bmatrix} x_1&y_1\\ x_2&y_2 \end{bmatrix}
  • 向量 vv 和 向量 ee 垂直
    x1x3+y1y3=0[x1,y1][x3y3]=0vT[x3y3]=vTe=0x_1x_3+y_1y_3=0\Rightarrow[x_1,y_1]\begin{bmatrix}x_3\\ y_3\end{bmatrix}=0\Rightarrow v^T\begin{bmatrix}x_3\\ y_3\end{bmatrix}=v^Te=0
  • 向量 ee 与张成的平面垂直
    ATe=0A^Te=0

  • bb 位于张成的向量内
    k1v+k2w=bk_1v+k_2w=b
    [v,w][k1k2]=b[v,w]\begin{bmatrix} k_1\\ k_2 \end{bmatrix}=b
    超定方程的最小二乘解的三维几何解释
  • 在matlab里有子程序对应 (ATA)1AT(A^TA)^{-1} A^T 伪逆 pinv(A)