Deep Learning(花书)教材笔记-Math and Machine Learning Basics(线性代数拾遗)

I. Linear Algebra

1. 基础概念回顾

  • scalar: 标量
  • vector: 矢量,an array of numbers.
  • matrix: 矩阵, 2-D array of numbers.
  • tensor: 张量, 更高维的一组数据集合。
  • identity Matricx:单位矩阵
  • inverse Matrix:逆矩阵,也称非奇异函数。当矩阵A的行列式\(|A|≠0\)时,则存在\(A^{-1}\).

2. Span

3. Norm

\(L^p\) norm 定义如右: \(||x||_p=(\sum_i|x_i|^p)^{\frac{1}{p}}\) for \(p∈R,p≥1\).

任何满足如下条件的函数都可视为norm:

  • \(f(x)=0 \, \Rightarrow x=0\)
  • \(f(x+y)≤f(x)+f(y)\) (三角不等式)
  • \(\forall α ∈R,f(αx)=|α|f(x)\)

1) \(L^2\) Norm

最常用的是二范式,即\(L^2\) norm,也称为Euclidean norm(欧几里得范数)。因为在机器学习中常用到求导,二范式求导之后只与输入数据本身有关,所以比较实用。

2) \(L^1\) Norm

但是二范式在零点附近增长很慢,而且有的机器学习应用需要在零点和非零点之间进行区分,此时二范式显得力不从心,所以我们可以选择一范式,即\(L^1\) norm,其表达式为:\(||x||_1=\sum_i|x_i|\).

3) \(L^0\) Norm

0范式表示矢量中非0的元素的个数。其实0范式这个说法是不严谨的,因为它不满足第三个条件,but whatever~

4) \(L^∞\) Norm

无穷大范式,也叫max norm,它表示矢量中所有元素绝对值的最大值,即

\[||x||_∞=max |x_i| \]

5) F norm

F norm全称是Frobenius Norm,其表达式如下:

\[||A||_F=\sqrt{\sum_{i,j}A_{i,j}^2} \]

4.特殊矩阵和向量

1) Diagonal matrix(对角矩阵)

定义: a matrix \(D\) is diagonal if and only if \(D_{i,j}=0\) for all \(i≠j\).

仔细看定义!!!这里并没有说必须是squre matrix(方阵),所以对角矩阵不一定是方阵,rectangle matrix也有可能是对角矩阵(只要对角线上不为0,其余部分都为0)。

2) Orthogonal Matrix(正交矩阵)

定义: 若\(A^TA=AA^T=I\),那么n阶实矩阵A则为正交矩阵。

注意矩阵A必须为方阵,另外有定义可知 \(A^{-1}=A^T\)

3) Orthonomal Matrix(标准正交矩阵)

定义: 满足正交矩阵的要求,且为x和y均为unit vector(单位矢量)。

5. Eigendecomposition(特征分解)

很多数学概念其实都可以分解成很小的组成部分,然后通过观察这些组成进而找出它们可能存在的通用的性质。例如对于一个整数12,我们会试着把它分解成12=2×2×3,由这个表达式我们可以得到一些有用的结论,例如12不能被5整除,任何数乘以12后都能被3整除等等。

很自然地,对于矩阵,我们也想看看他是否也能被拆分呢,所以就引入了特征分解的概念,通过特征分解我们会得到矩阵\(A\)的(一组)eigenvector(特征向量): \(v\)eigenvalue(特征值): \(λ\),它们满足如下等式:

\[Av=λv \]

(特征向量当然也可以在右边,但是通常更习惯于放在右边。)

假设矩阵\(A\)有n个线性独立的特征向量\(\{v^{(1)}, ..., v^{(n)}\}\)以及对应的特征值\(\{ λ_1, ...,λ_n \}\)。记
\(V=[v^{(1)}, ..., v^{(n)}],λ=[λ_1, ...,λ_n ]\),则矩阵A的特征分解如下:

\[A=Vdiag(λ)V^{-1} \]

另外实对称矩阵的特征分解用得比较多,表达式为\(A=Q\Lambda Q^{-1}\),\(Q\)表示由特征向量组成的正交矩阵,\(\Lambda\)表示对角矩阵,注意\(Q\)\(\Lambda\)的值是一一对应的。

  • 当一个矩阵的特征值都为正时,该矩阵则为positive definite(正定矩阵).
  • 当一个矩阵的特征值都大于等于0时,该矩阵则为positive semidefinite(半正定矩阵).
  • 当一个矩阵的特征值都为负时,该矩阵则为negative definite(负定矩阵).
  • 当一个矩阵的特征值都小于等于0时,该矩阵则为negative semidefinite(半负定矩阵).

6. Singular Value Decomposition(奇异值分解)

Singular Value Decomposition (SVD) 可以把一个矩阵分解得到 singular vectors和singular values。SVD可以像特征值分解一样帮助我们对一个矩阵进行分析,并且SVD适用性更广。每个实矩阵都能做SVD,但是不一定能做特征值分解。比如说如果一个矩阵不是方阵,那么就不能做特征分解,但是我们可以做SVD。

SVD分解后的矩阵表达式如下:

\[A=UDV^T \]

假设A是一个m×n矩阵,那么U定义为m×m矩阵,D是m×n矩阵,V是n×n矩阵。
除此以外

  • 矩阵U和V都是orthogonal matrix,其中矩阵U的列向量是left-singular vectors,矩阵V的列向量是right-singular vectors。矩阵A的left-singular vectors是矩阵\(A^TA\)的特征向量,right-singular vectors是矩阵\(AA^T\)的特征向量。矩阵A的非零奇异值是矩阵\(AA^T\)或者\(A^TA\)的平方根。
  • 矩阵D是diagonal matrix,注意不一定是方阵。D对角线上的即为矩阵A的奇异值(singular value)。

讲这么多,肯定对SVD还没有一个直观的理解,下面一节会介绍SVD的应用。

7. Moore-Penrose Pseudoinverse

我们在求一个矩阵的逆(matrix inverse)的时候,一般都需要规定这个矩阵是方阵。

假设有一个线性方程\(Ax=y\),为了解出这个方程,我们很直观地希望能够造出一个left-inverse矩阵B和A相乘,从而求出x,即\(x=By\)

如果A是一个非方阵的矩阵,当它的row大于column时,很有可能此时无解;而当row小于column时,可能有多解。

Moore-Penrose Pseudoinverse就是为了解决这个问题的,矩阵A的伪逆定义如下:

\[A^+=lim_{α\searrow{0}}(A^TA+αI)^{-1}A^T$$。 上面的公式实际很少用,一般都是使用SVD的公式,即 $$A^+=VD^+U^T\]

U,D,V是上节中提到的矩阵A的奇异分解。\(D^+\)是矩阵D的伪逆,它是首先将D的非零元素取倒数得到一个矩阵,然后将这个矩阵转置之后就得到了\(D^+\)

当矩阵A的row比column少时,使用伪逆可以得到很多解。但是,\(x=A^+y\)这个解是所有解中有最小Euclidean norm(\(||x||_2\))的。

当矩阵A的row比column多时,可能无解。但是使用伪逆求得的解x ,能使得\(Ax\)尽可能的接近\(y\),也就是说能使得\(||Ax-y||_2\)最小。

8. Trace Operator(迹)

trace运算符是将矩阵对角线上的所有元素求和,即\(Tr(A)=\sum_iA_{i,i}\)



MARSGGBO原创





2018-12-01



posted @ 2018-12-01 16:44  marsggbo  阅读(1237)  评论(0编辑  收藏  举报