Bootstrap

最小二乘解(Least-squares Minimization )

对于线性方程组,解的判别条件如下:
1. Ax=0 总有解,至少有零解
2. Am×nx=0
 当 r(A)=n ,只有零解
 当 r(A)<n ,有无穷多解
3. Am×nx=b
 当 r(A)r(A|b) ,无解
 当 r(A)=r(A|b)=n ,有唯一解
 当 r(A)=r(A|b)=r<n ,有无穷多解

  我一般我们会面临形如 Am×nx=b 的方程。我们考虑测量数据和我们需要的解的参数之间的关系,该方程的解可以分为以下几种情况:

  1. 如果 m<n ,未知数大于方程数。那么解不唯一,存在一个解矢量空间。
  2. 如果 m=n ,那么只要 A 可逆(非奇异,也就是满秩)就有唯一解,解为 x=A1b
  3. 如果 m>n ,方程数大于未知数。方程一般没有解,除非 b 属于 A 的列向量组成的子空间。

  我们考虑 mn 并且 r(A)=n 的情况。如果解不存在,我们找一个最接近 Am×nx=b 的解矢量仍然有意义,这个方程成为超定方程(方程大于未知数)。也就是说,我们寻找一个向量 x 使得 Axb 最小,这里的 表示矢量范数。这样的 x 称为该超定方程组的最小二乘解。接下来讨论三种解最小二乘的方法,分别用奇异值分解正规方程QR分解


奇异值分解

  奇异值分解(SVD)是最有用的矩阵分解方法中的一种。给定一个矩阵 Am×n(mn) ,存在一个正交矩阵 Um×m Vn×n ,有

A=UDVT=U[Σ0]VT

  其中 Σ=diag(σ1,σ2,,σn)Rn×n ,且 σ1σ2σn0 。上述分解就称为 奇异值分解 σ1,σ2,,σn 称为 A 奇异值,矩阵 U V 满足 UTU=I VTV=I 。由上式可知
ATA=VDUTUDVT=VΣ2VTAAT=U[
;