线性代数知识点总结
1. 基础概念和符号
线性代数提供了一种紧凑地表示和操作线性方程组的方法。 例如，以下方程组：
4
x
1
−
5
x
2
=
−
13
4x_1 − 5x_2 = −13
4x1−5x2=−13
− 2 x 1 + 3 x 2 = 9 −2x_1 + 3x_2 = 9 −2x1+3x2=9
这是两个方程和两个变量，正如你从高中代数中所知，你可以找到 x 1 x_1 x1 和 x 2 x_2 x2 的唯一解（除非方程以某种方式退化，例如，如果第二个方程只是第一个的倍数，但在上面的情况下，实际上只有一个唯一解）。 在矩阵表示法中，我们可以更紧凑地表达：
A x = b Ax= b Ax=b
with A = [ 4 − 5 − 2 3 ] , b = [ − 13 9 ] \text { with } A=\left[\begin{array}{cc}{4} & {-5} \\ {-2} & {3}\end{array}\right], b=\left[\begin{array}{c}{-13} \\ {9}\end{array}\right] with A=[4−2−53],b=[−139]
我们可以看到，这种形式的线性方程有许多优点（比如明显地节省空间）。
1.1 基本符号
我们使用以下符号：
-
A ∈ R m × n A \in \mathbb{R}^{m \times n} A∈Rm×n，表示 A A A 为由实数组成具有 m m m行和 n n n列的矩阵。
-
x ∈ R n x \in \mathbb{R}^{ n} x∈Rn，表示具有 n n n个元素的向量。 通常，向量 x x x将表示列向量: 即，具有 n n n行和 1 1 1列的矩阵。 如果我们想要明确地表示行向量: 具有 1 1 1 行和 n n n列的矩阵 - 我们通常写 x T x^T xT（这里 x T x^T xT x x x的转置）。
-
x i x_i xi表示向量 x x x的第 i i i个元素
x = [ x 1 x 2 ⋮ x n ] x=\left[\begin{array}{c}{x_{1}} \\ {x_{2}} \\ {\vdots} \\ {x_{n}}\end{array}\right] x=⎣⎢⎢⎢⎡x1x2⋮xn⎦⎥⎥⎥⎤
- 我们使用符号 a i j a_{ij} aij（或 A i j A_{ij} Aij, A i , j A_{i,j} Ai,j等）来表示第 i i i 行和第 j j j列中的 A A A 的元素：
A = [ a 11 a 12 ⋯ a 1 n a 21 a 22 ⋯ a 2 n ⋮ ⋮ ⋱ ⋮ a m 1 a m 2 ⋯ a m n ] A=\left[\begin{array}{cccc}{a_{11}} & {a_{12}} & {\cdots} & {a_{1 n}} \\ {a_{21}} & {a_{22}} & {\cdots} & {a_{2 n}} \\ {\vdots} & {\vdots} & {\ddots} & {\vdots} \\ {a_{m 1}} & {a_{m 2}} & {\cdots} & {a_{m n}}\end{array}\right] A=⎣⎢⎢⎢⎡a11a21⋮am1a12a22⋮am2⋯⋯⋱⋯a1na2n⋮amn⎦⎥⎥⎥⎤
- 我们用 a j a^j aj或者 A : , j A_{:,j} A:,j表示矩阵 A A A的第 j j j列：
A = [ ∣ ∣ ∣ a 1 a 2 ⋯ a n ∣ ∣ ∣ ] A=\left[\begin{array}{llll}{ |} & { |} & {} & { |} \\ {a^{1}} & {a^{2}} & {\cdots} & {a^{n}} \\ { |} & { |} & {} & { |}\end{array}\right] A=⎣⎡∣a1∣∣a2∣⋯∣an∣⎦⎤
- 我们用 a i T a^T_i aiT或者 A i , : A_{i,:} Ai,:表示矩阵 A A A的第 i i i行：
A = [ − a 1 T − − a 2 T − ⋮ − a m T − ] A=\left[\begin{array}{c}{-a_{1}^{T}-} \\ {-a_{2}^{T}-} \\ {\vdots} \\ {-a_{m}^{T}-}\end{array}\right] A=⎣⎢⎢⎢⎡−a1T−−a2T−⋮−amT−⎦⎥⎥⎥⎤
- 在许多情况下，将矩阵视为列向量或行向量的集合非常重要且方便。 通常，在向量而不是标量上操作在数学上（和概念上）更清晰。只要明确定义了符号，用于矩阵的列或行的表示方式并没有通用约定。
2.矩阵乘法
两个矩阵相乘，其中 A ∈ R m × n A \in \mathbb{R}^{m \times n} A∈Rm×n and B ∈ R n × p B \in \mathbb{R}^{n \times p} B∈Rn×p ，则：
C = A B ∈ R m × p C = AB \in \mathbb{R}^{m \times p} C=AB∈Rm×p
其中：
C i j = ∑ k = 1 n A i k B k j C_{i j}=\sum_{k=1}^{n} A_{i k} B_{k j} Cij=k=1∑nAikBkj
请注意，为了使矩阵乘积存在， A A A中的列数必须等于 B B B中的行数。有很多方法可以查看矩阵乘法，我们将从检查一些特殊情况开始。
2.1 向量-向量乘法
给定两个向量 x , y ∈ R n x, y \in \mathbb{R}^{n} x,y∈Rn, x T y x^T y xTy通常称为向量内积或者点积，结果是个实数。
x T y ∈ R = [ x 1 x 2 ⋯ x n ] [ y 1 y 2 ⋮ y n ] = ∑ i = 1 n x i y i x^{T} y \in \mathbb{R}=\left[\begin{array}{llll}{x_{1}} & {x_{2}} & {\cdots} & {x_{n}}\end{array}\right]\left[\begin{array}{c}{y_{1}} \\ {y_{2}} \\ {\vdots} \\ {y_{n}}\end{array}\right]=\sum_{i=1}^{n} x_{i} y_{i} xTy∈R=[x1x2⋯xn]⎣⎢⎢⎢⎡y1y2⋮yn⎦⎥⎥⎥⎤=i=1∑nxiyi
注意： x T y = y T x x^T y = y^Tx xTy=yTx 始终成立。
给定向量 x ∈ R m x \in \mathbb{R}^{m} x∈Rm, y ∈ R n y \in \mathbb{R}^{n} y∈Rn (他们的维度是否相同都没关系)， x y T ∈ R m × n xy^T \in \mathbb{R}^{m \times n} xyT∈Rm×n叫做向量外积 , 当 ( x y T ) i j = x i y j (xy^T)_{ij} = x_iy_j (xyT)ij=xiyj 的时候，它是一个矩阵。
x
y
T
∈
R
m
×
n
=
[
x
1
x
2
⋮
x
m
]
[
y
1
y
2
⋯
y
n
]
=
[
x
1
y
1
x
1
y
2
⋯
x
1
y
n
x
2
y
1
x
2
y
2
⋯
x
2
y
n
⋮
⋮
⋱
⋮
x
m
y
1
x
m
y
2
⋯
x
m
y
n
]
x y^{T} \in \mathbb{R}^{m \times n}=\left[\begin{array}{c}{x_{1}} \\ {x_{2}} \\ {\vdots} \\ {x_{m}}\end{array}\right]\left[\begin{array}{llll}{y_{1}} & {y_{2}} & {\cdots} & {y_{n}}\end{array}\right]=\left[\begin{array}{cccc}{x_{1} y_{1}} & {x_{1} y_{2}} & {\cdots} & {x_{1} y_{n}} \\ {x_{2} y_{1}} & {x_{2} y_{2}} & {\cdots} & {x_{2} y_{n}} \\ {\vdots} & {\vdots} & {\ddots} & {\vdots} \\ {x_{m} y_{1}} & {x_{m} y_{2}} & {\cdots} & {x_{m} y_{n}}\end{array}\right]
xyT∈Rm×n=⎣⎢⎢⎢⎡x1x2⋮xm⎦⎥⎥⎥⎤[y1y2⋯yn]=⎣⎢⎢⎢⎡x1y1x2y1⋮xmy1x1y2x2y2⋮xmy2⋯⋯⋱⋯x1ynx2yn⋮xmyn⎦⎥⎥⎥⎤
举一个外积如何使用的一个例子：让
1
∈
R
n
1\in R^{n}
1∈Rn表示一个
n
n
n维向量，其元素都等于1，此外，考虑矩阵
A
∈
R
m
×
n
A \in R^{m \times n}
A∈Rm×n，其列全部等于某个向量
x
∈
R
m
x \in R^{m}
x∈Rm。 我们可以使用外积紧凑地表示矩阵
A
A
A:
A = [ ∣ ∣ ∣ x x ⋯ x ∣ ∣ ∣ ] = [ x 1 x 1 ⋯ x 1 x 2 x 2 ⋯ x 2 ⋮ ⋮ ⋱ ⋮ x m x m ⋯ x m ] = [ x 1 x 2 ⋮ x m ] [ 1 1 ⋯ 1 ] = x 1 T A=\left[\begin{array}{llll}{ |} & { |} & {} & { |} \\ {x} & {x} & {\cdots} & {x} \\ { |} & { |} & {} & { |}\end{array}\right]=\left[\begin{array}{cccc}{x_{1}} & {x_{1}} & {\cdots} & {x_{1}} \\ {x_{2}} & {x_{2}} & {\cdots} & {x_{2}} \\ {\vdots} & {\vdots} & {\ddots} & {\vdots} \\ {x_{m}} & {x_{m}} & {\cdots} & {x_{m}}\end{array}\right]=\left[\begin{array}{c}{x_{1}} \\ {x_{2}} \\ {\vdots} \\ {x_{m}}\end{array}\right]\left[\begin{array}{lll}{1} & {1} & {\cdots} & {1}\end{array}\right]=x \mathbf{1}^{T} A=⎣⎡∣x∣∣x∣⋯∣x∣⎦⎤=⎣⎢⎢⎢⎡x1x2⋮xmx1x2⋮xm⋯⋯⋱⋯x1x2⋮xm⎦⎥⎥⎥⎤=⎣⎢⎢⎢⎡x1x2⋮xm⎦⎥⎥⎥⎤[11⋯1]=x1T
2.2 矩阵-向量乘法
给定矩阵 A ∈ R m × n A \in \mathbb{R}^{m \times n} A∈Rm×n，向量 x ∈ R n x \in \mathbb{R}^{n} x∈Rn , 它们的积是一个向量 y = A x ∈ R m y = Ax \in R^{m} y=Ax∈Rm。 有几种方法可以查看矩阵向量乘法，我们将依次查看它们中的每一种。
如果我们按行写 A A A，那么我们可以表示 A x Ax Ax为：
y = A x = [ − a 1 T − − a 2 T − ⋮ − a m T − ] x = [ a 1 T x a 2 T x ⋮ a m T x ] y=A x=\left[\begin{array}{ccc}{-} & {a_{1}^{T}} & {-} \\ {-} & {a_{2}^{T}} & {-} \\ {} & {\vdots} & {} \\ {-} & {a_{m}^{T}} & {-}\end{array}\right] x=\left[\begin{array}{c}{a_{1}^{T} x} \\ {a_{2}^{T} x} \\ {\vdots} \\ {a_{m}^{T} x}\end{array}\right] y=Ax=⎣⎢⎢⎢⎡−−−a1Ta2T⋮amT−−−⎦⎥⎥⎥⎤x=⎣⎢⎢⎢⎡a1Txa2Tx⋮amTx⎦⎥⎥⎥⎤
换句话说，第 i i i个 y y y是 A A A的第 i i i行和 x x x的内积，即： y i = y i = a i T x y_i = y_{i}=a_{i}^{T} x yi=yi=aiTx。
同样的， 可以把 A A A 写成列的方式，则公式如下：
y = A x = [ ∣ ∣ ∣ a 1 a 2 ⋯ a n ∣ ∣ ∣ ] [ x 1 x 2 ⋮ x n ] = [ a 1 x 1 ] + [ a 2 ] x 2 + ⋯ + [ a n ] x n y=A x=\left[\begin{array}{cccc} \mid & \mid & & \mid \\ a^{1} & a^{2} & \cdots & a^{n} \\ \mid & \mid & & \mid \end{array}\right]\left[\begin{array}{c} x_{1} \\ x_{2} \\ \vdots \\ x_{n} \end{array}\right]=\left[\begin{array}{l} a^{1} \\ x_{1} \end{array}\right]+\left[\begin{array}{l} a^{2} \end{array}\right] x_{2}+\cdots+\left[a^{n}\right] x_{n} y=Ax=⎣⎡∣a1∣∣a2∣⋯∣an∣⎦⎤⎣⎢⎢⎢⎡x1x2⋮xn⎦⎥⎥⎥⎤=[a1x1]+[a2]x2+⋯+[an]xn
换句话说， y y y是 A A A的列的线性组合，其中线性组合的系数由 x x x的元素给出。
到目前为止，我们一直在右侧乘以列向量，但也可以在左侧乘以行向量。 这是写的， y T = x T A y^T = x^TA yT=xTA 表示 A ∈ R m × n A \in \mathbb{R}^{m \times n} A∈Rm×n， x ∈ R m x \in \mathbb{R}^{m} x∈Rm， y ∈ R n y \in \mathbb{R}^{n} y∈Rn。 和以前一样，我们可以用两种可行的方式表达 y T y^T yT，这取决于我们是否根据行或列表达 A A A.
第一种情况，我们把 A A A用列表示：
y T = x T A = x T [ ∣ ∣ ∣ a 1 a 2 ⋯ a n ∣ ∣ ∣ ] = [ x T a 1 x T a 2 … x T a n ] y^{T}=x^{T} A=x^{T}\left[\begin{array}{cccc}{ |} & { |} & {} & { |} \\ {a^{1}} & {a^{2}} & {\cdots} & {a^{n}} \\ { |} & { |} & {} & { |}\end{array}\right]=\left[\begin{array}{cccc}{x^{T} a^{1}} & {x^{T} a^{2}} & {\dots} & {x^{T} a^{n}}\end{array}\right] yT=xTA=xT⎣⎡∣a1∣∣a2∣⋯∣an∣⎦⎤=[xTa1xTa2…xTan]
这表明 y T y^T yT的第 i i i个元素等于 x x x和 A A A的第 i i i列的内积。
最后，根据行表示 A A A，我们得到了向量-矩阵乘积的最终表示:
y
T
=
x
T
A
=
[
x
1
x
2
⋯
x
n
]
[
−
a
1
T
−
−
a
2
T
−
⋮
−
a
m
T
−
]
=
x
1
[
−
a
1
T
−
]
+
x
2
[
−
a
2
T
−
]
+
…
+
x
n
[
−
a
n
T
−
]
y^{T}=x^{T} A=\left[\begin{array}{llll} x_{1} & x_{2} & \cdots & x_{n} \end{array}\right]\left[\begin{array}{c} -a_{1}^{T}- \\ -a_{2}^{T}- \\ \vdots \\ -a_{m}^{T}- \end{array}\right]=x_{1}\left[-a_{1}^{T}-\right]+x_{2}\left[-a_{2}^{T}-\right]+\ldots+x_{n}\left[-a_{n}^{T}-\right]
yT=xTA=[x1x2⋯xn]⎣⎢⎢⎢⎡−a1T−−a2T−⋮−amT−⎦⎥⎥⎥⎤=x1[−a1T−]+x2[−a2T−]+…+xn[−anT−]
所以我们看到
y
T
y^T
yT是
A
A
A的行的线性组合，其中线性组合的系数由
x
x
x的元素给出。
2.3 矩阵-矩阵乘法
有了这些知识，我们现在可以看看四种不同的（形式不同，但结果是相同的）矩阵-矩阵乘法：也就是本节开头所定义的 C = A B C=AB C=AB的乘法。
首先，我们可以将矩阵 - 矩阵乘法视为一组向量-向量乘积。 从定义中可以得出：最明显的观点是$C
的
的
的( i，j )
元
素
等
于
元素等于
元素等于A
的
第
的第
的第i
行
和
行和
行和B
的
的
的的
的的j$列的内积。如下面的公式所示：
C
=
A
B
=
[
−
a
1
T
−
−
a
2
T
−
⋮
−
a
m
T
−
]
[
∣
∣
∣
b
1
b
2
⋯
b
p
∣
∣
∣
]
=
[
a
1
T
b
1
a
1
T
b
2
⋯
a
1
T
b
p
a
2
T
b
1
a
2
T
b
2
⋯
a
2
T
b
p
⋮
⋮
⋱
⋮
a
m
T
b
1
a
m
T
b
2
⋯
a
m
T
b
p
]
C=A B=\left[\begin{array}{cc}{-} & {a_{1}^{T}} &{-} \\ {-} & {a_{2}^{T}} &{-} \\ {} & {\vdots} \\ {-} & {a_{m}^{T}} &{-} \end{array}\right]\left[\begin{array}{cccc}{ |} & { |} & {} & { |} \\ {b_{1}} & {b_{2}} & {\cdots} & {b_{p}} \\ { |} & { |} & {} & { |}\end{array}\right]=\left[\begin{array}{cccc}{a_{1}^{T} b_{1}} & {a_{1}^{T} b_{2}} & {\cdots} & {a_{1}^{T} b_{p}} \\ {a_{2}^{T} b_{1}} & {a_{2}^{T} b_{2}} & {\cdots} & {a_{2}^{T} b_{p}} \\ {\vdots} & {\vdots} & {\ddots} & {\vdots} \\ {a_{m}^{T} b_{1}} & {a_{m}^{T} b_{2}} & {\cdots} & {a_{m}^{T} b_{p}}\end{array}\right]
C=AB=⎣⎢⎢⎢⎡−−−a1Ta2T⋮amT−−−⎦⎥⎥⎥⎤⎣⎡∣b1∣∣b2∣⋯∣bp∣⎦⎤=⎣⎢⎢⎢⎡a1Tb1a2Tb1⋮amTb1a1Tb2a2Tb2⋮amTb2⋯⋯⋱⋯a1Tbpa2Tbp⋮amTbp⎦⎥⎥⎥⎤
这里的
A
∈
R
m
×
n
A \in \mathbb{R}^{m\times n}
A∈Rm×n ，
B
∈
R
n
×
p
B \in \mathbb{R}^{n \times p}
B∈Rn×p，
a
i
∈
R
n
a_i \in \mathbb{R}^n
ai∈Rn ，
b
j
∈
R
n
×
p
b^j \in \mathbb{R}^{n \times p}
bj∈Rn×p， 这里的
A
∈
R
m
×
n
A \in \mathbb{R}^ {m \times n}
A∈Rm×n，
B
∈
R
n
×
p
B \in \mathbb{R}^ {n \times p}
B∈Rn×p，
a
i
∈
R
n
a_i \in \mathbb{R} ^ n
ai∈Rn，
b
j
∈
R
n
×
p
b ^ j \in \mathbb{R} ^ {n \times p}
bj∈Rn×p，所以它们可以计算内积。 我们用通常用行表示
A
A
A而用列表示
B
B
B。
或者，我们可以用列表示
A
A
A，用行表示
B
B
B，这时
A
B
AB
AB是求外积的和。公式如下：
C
=
A
B
=
[
∣
∣
∣
a
1
a
2
⋯
a
n
∣
∣
∣
]
[
−
b
1
T
−
−
b
2
T
−
⋮
−
b
n
T
−
]
=
∑
i
=
1
n
a
i
b
i
T
C=A B=\left[\begin{array}{cccc}{ |} & { |} & {} & { |} \\ {a_{1}} & {a_{2}} & {\cdots} & {a_{n}} \\ { |} & { |} & {} & { |}\end{array}\right]\left[\begin{array}{c}{-}& {b_{1}^{T}}&{-} \\ {-}& {b_{2}^{T}}&{-} \\ {\vdots} \\{-}& {b_{n}^{T}}&{-}\end{array}\right]=\sum_{i=1}^{n} a_{i} b_{i}^{T}
C=AB=⎣⎡∣a1∣∣a2∣⋯∣an∣⎦⎤⎣⎢⎢⎢⎡−−⋮−b1Tb2TbnT−−−⎦⎥⎥⎥⎤=i=1∑naibiT
换句话说，
A
B
AB
AB等于所有的
A
A
A的第
i
i
i列和
B
B
B第
i
i
i行的外积的和。因此，在这种情况下，
a
i
∈
R
m
a_i \in \mathbb{R}^ m
ai∈Rm和
b
i
∈
R
p
b_i \in \mathbb{R}^p
bi∈Rp， 外积
a
i
b
i
T
a^ib_i^T
aibiT的维度是
m
×
p
m×p
m×p，与
C
C
C的维度一致。
其次，我们还可以将矩阵 - 矩阵乘法视为一组矩阵向量积。如果我们把 B B B用列表示，我们可以将 C C C的列视为 A A A和 B B B的列的矩阵向量积。公式如下：
C
=
A
B
=
A
[
∣
∣
∣
b
1
b
2
⋯
b
p
∣
∣
∣
]
=
[
∣
∣
∣
A
b
1
A
b
2
⋯
A
b
p
∣
∣
∣
]
C=A B=A\left[\begin{array}{cccc}{ |} & { |} & {} & { |} \\ {b_{1}} & {b_{2}} & {\cdots} & {b_{p}} \\ { |} & { |} & {} & { |}\end{array}\right]=\left[\begin{array}{cccc}{ |} & { |} & {} & { |} \\ {A b_{1}} & {A b_{2}} & {\cdots} & {A b_{p}} \\ { |} & { |} & {} & { |}\end{array}\right]
C=AB=A⎣⎡∣b1∣∣b2∣⋯∣bp∣⎦⎤=⎣⎡∣Ab1∣∣Ab2∣⋯∣Abp∣⎦⎤
这里
C
C
C的第
i
i
i列由矩阵向量乘积给出，右边的向量为
c
i
=
A
b
i
c_i = Ab_i
ci=Abi。 这些矩阵向量乘积可以使用前一小节中给出的两个观点来解释。
最后，我们有类似的观点，我们用行表示
A
A
A，
C
C
C的行作为
A
A
A和
C
C
C行之间的矩阵向量积。公式如下：
C
=
A
B
=
[
−
a
1
T
−
−
a
2
T
−
⋮
−
a
m
T
−
]
B
=
[
−
a
1
T
B
−
−
a
2
T
B
−
⋮
−
a
m
T
B
−
]
C=A B=\left[\begin{array}{ccc}{-} & {a_{1}^{T}} & {-} \\ {-} & {a_{2}^{T}} & {-} \\ {} & {\vdots} & {} \\ {-} & {a_{m}^{T}} & {-}\end{array}\right] B=\left[\begin{array}{c} {-} & {a_{1}^{T} B} & {-}\\ {-} & {a_{2}^{T} B} & {-} \\ {\vdots} \\ {-} & {a_{m}^{T} B}& {-}\end{array}\right]
C=AB=⎣⎢⎢⎢⎡−−−a1Ta2T⋮amT−−−⎦⎥⎥⎥⎤B=⎣⎢⎢⎢⎡−−⋮−a1TBa2TBamTB−−−⎦⎥⎥⎥⎤
这里第 i i i行的 C C C由左边的向量的矩阵向量乘积给出： c i T = a i T B c_i^T = a_i^T B ciT=aiTB
将矩阵乘法剖析到如此大的程度似乎有点过分，特别是当所有这些观点都紧跟在我们在本节开头给出的初始定义（在一行数学中）之后。
这些不同方法的直接优势在于它们允许您在向量的级别/单位而不是标量上进行操作。 为了完全理解线性代数而不会迷失在复杂的索引操作中，关键是要用尽可能多的概念进行操作。
实际上所有的线性代数都处理某种矩阵乘法，花一些时间对这里提出的观点进行直观的理解是非常必要的。
除此之外，了解一些更高级别的矩阵乘法的基本属性是很有必要的：
-
矩阵乘法结合律: ( A B ) C = A ( B C ) (AB)C = A(BC) (AB)C=A(BC)
-
矩阵乘法分配律: A ( B + C ) = A B + A C A(B + C) = AB + AC A(B+C)=AB+AC
-
矩阵乘法通常不是可交换的; 也就是说，通常 A B ≠ B A AB \ne BA AB=BA。 （例如，假设 A ∈ R m × n A \in \mathbb{R}^ {m \times n} A∈Rm×n, B ∈ R n × p B \in \mathbb{R}^ {n \times p} B∈Rn×p，如果 m m m和 q q q不相等，矩阵乘积 B A BA BA甚至不存在！）
如果您不熟悉这些属性，请花点时间自己验证它们。 例如，为了检查矩阵乘法的相关性，假设 A ∈ R m × n A \in \mathbb{R}^ {m \times n} A∈Rm×n, B ∈ R n × p B \in \mathbb{R}^ {n \times p} B∈Rn×p， C ∈ R p × q C \in \mathbb{R}^ {p \times q} C∈Rp×q。 注意 A B ∈ R m × p AB \in \mathbb{R}^ {m \times p} AB∈Rm×p，所以 ( A B ) C ∈ R m × q (AB)C \in \mathbb{R}^ {m \times q} (AB)C∈Rm×q。 类似地， B C ∈ R n × q BC \in \mathbb{R}^ {n \times q} BC∈Rn×q，所以 A ( B C ) ∈ R m × q A(BC) \in \mathbb{R}^ {m \times q} A(BC)∈Rm×q。 因此，所得矩阵的维度一致。 为了表明矩阵乘法是相关的，足以检查 ( A B ) C (AB)C (AB)C的第 ( i , j ) (i,j) (i,j)个元素是否等于 A ( B C ) A(BC) A(BC)的第 ( i , j ) (i,j) (i,j)个元素。 我们可以使用矩阵乘法的定义直接验证这一点：
( ( A B ) C ) i j = ∑ k = 1 p ( A B ) i k C k j = ∑ k = 1 p ( ∑ l = 1 n A i l B l k ) C k j = ∑ k = 1 p ( ∑ l = 1 n A i l B l k C k j ) = ∑ l = 1 n ( ∑ k = 1 p A i l B l k C k j ) = ∑ l = 1 n A i l ( ∑ k = 1 p B l k C k j ) = ∑ l = 1 n A i l ( B C ) l j = ( A ( B C ) ) i j \begin{aligned}((A B) C)_{i j} &=\sum_{k=1}^{p}(A B)_{i k} C_{k j}=\sum_{k=1}^{p}\left(\sum_{l=1}^{n} A_{i l} B_{l k}\right) C_{k j} \\ &=\sum_{k=1}^{p}\left(\sum_{l=1}^{n} A_{i l} B_{l k} C_{k j}\right)=\sum_{l=1}^{n}\left(\sum_{k=1}^{p} A_{i l} B_{l k} C_{k j}\right) \\ &=\sum_{l=1}^{n} A_{i l}\left(\sum_{k=1}^{p} B_{l k} C_{k j}\right)=\sum_{l=1}^{n} A_{i l}(B C)_{l j}=(A(B C))_{i j} \end{aligned} ((AB)C)ij=k=1∑p(AB)ikCkj=k=1∑p(l=1∑nAilBlk)Ckj=k=1∑p(l=1∑nAilBlkCkj)=l=1∑n(k=1∑pAilBlkCkj)=l=1∑nAil(k=1∑pBlkCkj)=l=1∑nAil(BC)lj=(A(BC))ij
3 运算和属性
在本节中，我们介绍矩阵和向量的几种运算和属性。 希望能够为您复习大量此类内容，这些笔记可以作为这些主题的参考。
3.1 单位矩阵和对角矩阵
单位矩阵,
I
∈
R
n
×
n
I \in \mathbb{R}^{n \times n}
I∈Rn×n，它是一个方阵，对角线的元素是1，其余元素都是0：
I
i
j
=
{
1
i
=
j
0
i
≠
j
I_{i j}=\left\{\begin{array}{ll}{1} & {i=j} \\ {0} & {i \neq j}\end{array}\right.
Iij={10i=ji=j
对于所有
A
∈
R
m
×
n
A \in \mathbb{R}^ {m \times n}
A∈Rm×n，有：
A
I
=
A
=
I
A
AI = A = IA
AI=A=IA
注意，在某种意义上，单位矩阵的表示法是不明确的，因为它没有指定
I
I
I的维数。通常，
I
I
I的维数是从上下文推断出来的，以便使矩阵乘法成为可能。 例如，在上面的等式中，
A
I
=
A
AI = A
AI=A中的I是
n
×
n
n\times n
n×n矩阵，而
A
=
I
A
A = IA
A=IA中的
I
I
I是
m
×
m
m\times m
m×m矩阵。
对角矩阵是一种这样的矩阵：对角线之外的元素全为0。对角阵通常表示为：
D
=
d
i
a
g
(
d
1
,
d
2
,
.
.
.
,
d
n
)
D= diag(d_1, d_2, . . . , d_n)
D=diag(d1,d2,...,dn)，其中：
D
i
j
=
{
d
i
i
=
j
0
i
≠
j
D_{i j}=\left\{\begin{array}{ll}{d_{i}} & {i=j} \\ {0} & {i \neq j}\end{array}\right.
Dij={di0i=ji=j
很明显：单位矩阵
I
=
d
i
a
g
(
1
,
1
,
.
.
.
,
1
)
I = diag(1, 1, . . . , 1)
I=diag(1,1,...,1)。
3.2 转置
矩阵的转置是指翻转矩阵的行和列。
给定一个矩阵：
A
∈
R
m
×
n
A \in \mathbb{R}^ {m \times n}
A∈Rm×n, 它的转置为
n
×
m
n \times m
n×m的矩阵
A
T
∈
R
n
×
m
A^T \in \mathbb{R}^ {n \times m}
AT∈Rn×m ，其中的元素为：
(
A
T
)
i
j
=
A
j
i
(A^T)_{ij} = A_{ji}
(AT)ij=Aji
事实上，我们在描述行向量时已经使用了转置，因为列向量的转置自然是行向量。
转置的以下属性很容易验证：
- ( A T ) T = A (A^T )^T = A (AT)T=A
- ( A B ) T = B T A T (AB)^T = B^T A^T (AB)T=BTAT
- ( A + B ) T = A T + B T (A + B)^T = A^T + B^T (A+B)T=AT+BT
3.3 对称矩阵
如果
A
=
A
T
,
A=A^{T},
A=AT, 则矩阵
A
∈
R
n
×
n
A \in \mathbb{R}^{n \times n}
A∈Rn×n 是对称矩阵。如果
A
=
−
A
T
,
A=-A^{T},
A=−AT, 它是反对称的。很容易证明, 对于任何矩阵
A
∈
R
n
×
n
,
A \in \mathbb{R}^{n \times n},
A∈Rn×n, 矩阵
A
+
A
T
A+A^{T}
A+AT 是对称的, 矩车
A
−
A
T
A-A^{T}
A−AT 是反对称的。由此得出, 任何方短阵
A
∈
R
n
×
n
A \in \mathbb{R}^{n \times n}
A∈Rn×n 可以表示为对称矩阵和反对称矩阵的和, 所以:
A
=
1
2
(
A
+
A
T
)
+
1
2
(
A
−
A
T
)
A=\frac{1}{2}(A+A^T)+\frac{1}{2}(A-A^T)
A=21(A+AT)+21(A−AT)
上面公式的右边的第一个矩阵是对称矩阵，而第二个矩阵是反对称矩阵。 事实证明，对称矩阵在实践中用到很多，它们有很多很好的属性，我们很快就会看到它们。
通常将大小为
n
n
n的所有对称矩阵的集合表示为
S
n
\mathbb{S}^n
Sn，因此
A
∈
S
n
A \in \mathbb{S}^n
A∈Sn意味着
A
A
A是对称的
n
×
n
n\times n
n×n矩阵。
3.4 矩阵的迹
方矩阵
A
∈
R
n
×
n
A \in \mathbb{R}^ {n \times n}
A∈Rn×n的迹，表示为
tr
(
A
)
\operatorname{tr} (A)
tr(A)（或者只是
tr
A
\operatorname{tr} A
trA，如果括号显然是隐含的），是矩阵中对角元素的总和：
tr
A
=
∑
i
=
1
n
A
i
i
\operatorname{tr} A=\sum_{i=1}^{n} A_{i i}
trA=i=1∑nAii
迹具有以下属性（如下所示）：
-
对于矩阵 A ∈ R n × n A \in \mathbb{R}^ {n \times n} A∈Rn×n，则： tr A = tr A T \operatorname{tr}A =\operatorname{tr}A^T trA=trAT
-
对于矩阵 A , B ∈ R n × n A,B \in \mathbb{R}^ {n \times n} A,B∈Rn×n，则： tr ( A + B ) = tr A + tr B \operatorname{tr}(A + B) = \operatorname{tr}A + \operatorname{tr}B tr(A+B)=trA+trB
-
对于矩阵 A ∈ R n × n A \in \mathbb{R}^ {n \times n} A∈Rn×n， t ∈ R t \in \mathbb{R} t∈R，则： tr ( t A ) = t tr A \operatorname{tr}(tA) = t\operatorname{tr}A tr(tA)=ttrA.
-
对于矩阵 A A A, B B B， A B AB AB 为方阵, 则： tr A B = tr B A \operatorname{tr}AB = \operatorname{tr}BA trAB=trBA
-
对于矩阵 A A A, B B B, C C C, A B C ABC ABC为方阵, 则： tr A B C = tr B C A = tr C A B \operatorname{tr}ABC = \operatorname{tr}BCA=\operatorname{tr}CAB trABC=trBCA=trCAB, 同理，更多矩阵的积也是有这个性质。
作为如何证明这些属性的示例，我们将考虑上面给出的第四个属性。 假设 A ∈ R m × n A \in \mathbb{R}^ {m \times n} A∈Rm×n和 B ∈ R n × m B \in \mathbb{R}^ {n \times m} B∈Rn×m（因此 A B ∈ R m × m AB \in \mathbb{R}^ {m \times m} AB∈Rm×m是方阵）。 观察到 B A ∈ R n × n BA \in \mathbb{R}^ {n \times n} BA∈Rn×n也是一个方阵，因此对它们进行迹的运算是有意义的。 要证明 tr A B = tr B A \operatorname{tr}AB = \operatorname{tr}BA trAB=trBA，请注意：
tr A B = ∑ i = 1 m ( A B ) i i = ∑ i = 1 m ( ∑ j = 1 n A i j B j i ) = ∑ i = 1 m ∑ j = 1 n A i j B j i = ∑ j = 1 n ∑ i = 1 m B j i A i j = ∑ j = 1 n ( ∑ i = 1 m B j i A i j ) = ∑ j = 1 n ( B A ) j j = tr B A \begin{aligned} \operatorname{tr} A B &=\sum_{i=1}^{m}(A B)_{i i}=\sum_{i=1}^{m}\left(\sum_{j=1}^{n} A_{i j} B_{j i}\right) \\ &=\sum_{i=1}^{m} \sum_{j=1}^{n} A_{i j} B_{j i}=\sum_{j=1}^{n} \sum_{i=1}^{m} B_{j i} A_{i j} \\ &=\sum_{j=1}^{n}\left(\sum_{i=1}^{m} B_{j i} A_{i j}\right)=\sum_{j=1}^{n}(B A)_{j j}=\operatorname{tr} B A \end{aligned} trAB=i=1∑m(AB)ii=i=1∑m(j=1∑nAijBji)=i=1∑mj=1∑nAijBji=j=1∑ni=1∑mBjiAij=j=1∑n(i=1∑mBjiAij)=j=1∑n(BA)jj=trBA
这里，第一个和最后两个等式使用迹运算符和矩阵乘法的定义，重点在第四个等式，使用标量乘法的可交换性来反转每个乘积中的项的顺序，以及标量加法的可交换性和相关性，以便重新排列求和的顺序。
3.5 范数
向量的范数
∥
x
∥
\|x\|
∥x∥是非正式度量的向量的“长度” 。 例如，我们有常用的欧几里德或
ℓ
2
\ell_{2}
ℓ2范数，
∥
x
∥
2
=
∑
i
=
1
n
x
i
2
\|x\|_{2}=\sqrt{\sum_{i=1}^{n} x_{i}^{2}}
∥x∥2=i=1∑nxi2
注意：
∥
x
∥
2
2
=
x
T
x
\|x\|_{2}^{2}=x^{T} x
∥x∥22=xTx
更正式地，范数是满足4个属性的函数（ f : R n → R f : \mathbb{R}^{n} \rightarrow \mathbb{R} f:Rn→R）：
- 对于所有的 x ∈ R n x \in \mathbb{R}^ {n} x∈Rn, f ( x ) ≥ 0 f(x) \geq 0 f(x)≥0(非负).
- 当且仅当 x = 0 x = 0 x=0 时， f ( x ) = 0 f(x) = 0 f(x)=0 (明确性).
- 对于所有 x ∈ R n x \in \mathbb{R}^ {n} x∈Rn, t ∈ R t\in \mathbb{R} t∈R，则 f ( t x ) = ∣ t ∣ f ( x ) f(tx) = \left| t \right|f(x) f(tx)=∣t∣f(x) (正齐次性).
- 对于所有 x , y ∈ R n x,y \in \mathbb{R}^ {n} x,y∈Rn, f ( x + y ) ≤ f ( x ) + f ( y ) f(x + y) \leq f(x) + f(y) f(x+y)≤f(x)+f(y) (三角不等式)
其他范数的例子是
ℓ
1
\ell_1
ℓ1范数:
∥
x
∥
1
=
∑
i
=
1
n
∣
x
i
∣
\|x\|_{1}=\sum_{i=1}^{n}\left|x_{i}\right|
∥x∥1=i=1∑n∣xi∣
和
ℓ
∞
\ell_{\infty }
ℓ∞范数：
∥
x
∥
∞
=
max
i
∣
x
i
∣
\|x\|_{\infty}=\max _{i}\left|x_{i}\right|
∥x∥∞=imax∣xi∣
事实上，到目前为止所提出的所有三个范数都是
ℓ
p
\ell_p
ℓp范数族的例子，它们由实数
p
≥
1
p \geq 1
p≥1参数化，并定义为：
∥
x
∥
p
=
(
∑
i
=
1
n
∣
x
i
∣
p
)
1
/
p
\|x\|_{p}=\left(\sum_{i=1}^{n}\left|x_{i}\right|^{p}\right)^{1 / p}
∥x∥p=(i=1∑n∣xi∣p)1/p
也可以为矩阵定义范数，例如Frobenius范数:
∥
A
∥
F
=
∑
i
=
1
m
∑
j
=
1
n
A
i
j
2
=
tr
(
A
T
A
)
\|A\|_{F}=\sqrt{\sum_{i=1}^{m} \sum_{j=1}^{n} A_{i j}^{2}}=\sqrt{\operatorname{tr}\left(A^{T} A\right)}
∥A∥F=i=1∑mj=1∑nAij2
=tr(ATA)
许多其他更多的范数，但它们超出了这个复习材料的范围。
3.6 线性相关性和秩
一组向量
x
1
,
x
2
,
⋯
x
n
∈
R
{x_1,x_2, \cdots x_n} \in \mathbb{R}
x1,x2,⋯xn∈R， 如果没有向量可以表示为其余向量的线性组合，则称称该向量是线性无相关的。 相反，如果属于该组的一个向量可以表示为其余向量的线性组合，则称该向量是线性相关的。 也就是说，如果：
x
n
=
∑
i
=
1
n
−
1
α
i
x
i
x_{n}=\sum_{i=1}^{n-1} \alpha_{i} x_{i}
xn=i=1∑n−1αixi
对于某些标量值
α
1
,
⋯
α
n
−
1
∈
R
\alpha_1,\cdots \alpha_n-1 \in \mathbb{R}
α1,⋯αn−1∈R，要么向量
x
1
,
x
2
,
⋯
x
n
x_1,x_2, \cdots x_n
x1,x2,⋯xn是线性相关的; 否则，向量是线性无关的。 例如，向量：
x
1
=
[
1
2
3
]
x
2
=
[
4
1
5
]
x
3
=
[
2
−
3
−
1
]
x_{1}=\left[\begin{array}{l}{1} \\ {2} \\ {3}\end{array}\right] \quad x_{2}=\left[\begin{array}{c}{4} \\ {1} \\ {5}\end{array}\right] \quad x_{3}=\left[\begin{array}{c}{2} \\ {-3} \\ {-1}\end{array}\right]
x1=⎣⎡123⎦⎤x2=⎣⎡415⎦⎤x3=⎣⎡2−3−1⎦⎤
是线性相关的，因为：
x
3
=
−
2
x
1
+
x
2
x_3=-2x_1+x_2
x3=−2x1+x2。
矩阵 A ∈ R m × n A \in \mathbb{R}^{m \times n} A∈Rm×n的列秩是构成线性无关集合的 A A A的最大列子集的大小。 由于术语的多样性，这通常简称为 A A A的线性无关列的数量。同样，行秩是构成线性无关集合的 A A A的最大行数。 对于任何矩阵 A ∈ R m × n A \in \mathbb{R}^{m \times n} A∈Rm×n，事实证明 A A A的列秩等于 A A A的行秩（尽管我们不会证明这一点），因此两个量统称为 A A A的秩，用 rank ( A ) \text{rank}(A) rank(A)表示。 以下是秩的一些基本属性：
- 对于 A ∈ R m × n A \in \mathbb{R}^{m \times n} A∈Rm×n， rank ( A ) ≤ m i n ( m , n ) \text{rank}(A) \leq min(m,n) rank(A)≤min(m,n)，如果 ( A ) = min ( m , n ) \text(A) = \text{min} (m, n) (A)=min(m,n)，则 A A A 被称作满秩。
- 对于 A ∈ R m × n A \in \mathbb{R}^{m \times n} A∈Rm×n， rank ( A ) = rank ( A T ) \text{rank}(A) = \text{rank}(A^T) rank(A)=rank(AT)
- 对于 A ∈ R m × n A \in \mathbb{R}^{m \times n} A∈Rm×n, B ∈ R n × p B \in \mathbb{R}^{n \times p} B∈Rn×p , rank ( A B ) ≤ min ( rank ( A ) , rank ( B ) ) \text{rank}(AB) \leq \text{min} ( \text{rank}(A), \text{rank}(B)) rank(AB)≤min(rank(A),rank(B))
- 对于 A , B ∈ R m × n A,B \in \mathbb{R}^{m \times n} A,B∈Rm×n， rank ( A + B ) ≤ rank ( A ) + rank ( B ) \text{rank}(A + B) \leq \text{rank}(A) + \text{rank}(B) rank(A+B)≤rank(A)+rank(B)
3.7 方阵的逆
方阵
A
∈
R
n
×
n
A \in \mathbb{R}^{n \times n}
A∈Rn×n的倒数表示为
A
−
1
A^{-1}
A−1，并且是这样的独特矩阵:
A
−
1
A
=
I
=
A
A
−
1
A^{-1}A=I=AA^{-1}
A−1A=I=AA−1
请注意，并非所有矩阵都具有逆。 例如，非方形矩阵根据定义没有逆。 然而，对于一些方形矩阵
A
A
A，可能仍然存在
A
−
1
A^{-1}
A−1可能不存在的情况。 特别是，如果
A
−
1
A^{-1}
A−1存在，我们说
A
A
A是可逆的或非奇异的，否则就是不可逆或奇异的。
为了使方阵A具有逆
A
−
1
A^{-1}
A−1，则
A
A
A必须是满秩。 我们很快就会发现，除了满秩之外，还有许多其它的充分必要条件。
以下是逆的属性; 假设
A
,
B
∈
R
n
×
n
A,B \in \mathbb{R}^{n \times n}
A,B∈Rn×n，而且是非奇异的：
- ( A − 1 ) − 1 = A (A^{-1})^{-1} = A (A−1)−1=A
- ( A B ) − 1 = B − 1 A − 1 (AB)^{-1} = B^{-1}A^{-1} (AB)−1=B−1A−1
-
(
A
−
1
)
T
=
(
A
T
)
−
1
(A^{-1})^{T} =(A^{T})^{-1}
(A−1)T=(AT)−1因此，该矩阵通常表示为
A
−
T
A^{-T}
A−T。
作为如何使用逆的示例，考虑线性方程组， A x = b Ax = b Ax=b，其中 A ∈ R n × n A \in \mathbb{R}^{n \times n} A∈Rn×n， x , b ∈ R x,b\in \mathbb{R} x,b∈R， 如果 A A A是非奇异的（即可逆的），那么 x = A − 1 b x = A^{-1}b x=A−1b。 （如果 A ∈ R m × n A \in \mathbb{R}^{m \times n} A∈Rm×n不是方阵，这公式还有用吗？）
3.8 正交阵
如果 x T y = 0 x^Ty=0 xTy=0，则两个向量 x , y ∈ R n x,y\in \mathbb{R}^{n} x,y∈Rn 是正交的。如果 ∥ x ∥ 2 = 1 \|x\|_2=1 ∥x∥2=1，则向量 x ∈ R n x\in \mathbb{R}^{n} x∈Rn 被归一化。如果一个方阵 U ∈ R n × n U\in \mathbb{R}^{n \times n} U∈Rn×n的所有列彼此正交并被归一化（这些列然后被称为正交），则方阵 U U U是正交阵（注意在讨论向量时的意义不一样）。
它可以从正交性和正态性的定义中得出:
U
T
U
=
I
=
U
U
T
U^ TU = I = U U^T
UTU=I=UUT
换句话说，正交矩阵的逆是其转置。 注意，如果
U
U
U不是方阵 :即，
U
∈
R
m
×
n
U\in \mathbb{R}^{m \times n}
U∈Rm×n，
n
<
m
n <m
n<m ，但其列仍然是正交的，则
U
T
U
=
I
U^TU = I
UTU=I，但是
U
U
T
≠
I
UU^T \neq I
UUT=I。我们通常只使用术语"正交"来描述先前的情况 ，其中
U
U
U是方阵。
正交矩阵的另一个好的特性是在具有正交矩阵的向量上操作不会改变其欧几里德范数，即:
∥
U
x
∥
2
=
∥
x
∥
2
\|U x\|_{2}=\|x\|_{2}
∥Ux∥2=∥x∥2
对于任何
x
∈
R
x\in \mathbb{R}
x∈R ,
U
∈
R
n
U\in \mathbb{R}^{n}
U∈Rn是正交的。
3.9 矩阵的值域和零空间
一组向量
{
x
1
,
…
x
n
}
\{x_{1}, \ldots x_{n}\}
{x1,…xn}是可以表示为
{
x
1
,
…
x
n
}
\{x_{1}, \ldots x_{n}\}
{x1,…xn}的线性组合的所有向量的集合。 即：
span
(
{
x
1
,
…
x
n
}
)
=
{
v
:
v
=
∑
i
=
1
n
α
i
x
i
,
α
i
∈
R
}
\operatorname{span}\left(\left\{x_{1}, \ldots x_{n}\right\}\right)=\left\{v : v=\sum_{i=1}^{n} \alpha_{i} x_{i}, \quad \alpha_{i} \in \mathbb{R}\right\}
span({x1,…xn})={v:v=i=1∑nαixi,αi∈R}
可以证明，如果
{
x
1
,
…
x
n
}
\{x_{1}, \ldots x_{n}\}
{x1,…xn}是一组
n
n
n个线性无关的向量，其中每个
x
i
∈
R
n
x_i \in \mathbb{R}^{n}
xi∈Rn，则
span
(
{
x
1
,
…
x
n
}
)
=
R
n
\text{span}(\{x_{1}, \ldots x_{n}\})=\mathbb{R}^{n}
span({x1,…xn})=Rn。 换句话说，任何向量
v
∈
R
n
v\in \mathbb{R}^{n}
v∈Rn都可以写成
x
1
x_1
x1到
x
n
x_n
xn的线性组合。
向量 y ∈ R m y\in \mathbb{R}^{m} y∈Rm投影到 { x 1 , … x n } \{x_{1}, \ldots x_{n}\} {x1,…xn}（这里我们假设 x i ∈ R m x_i \in \mathbb{R}^{m} xi∈Rm）得到向量 v ∈ span ( { x 1 , … , x n } ) v \in \operatorname{span}(\{x_{1}, \ldots, x_{n}\}) v∈span({x1,…,xn})，由欧几里德范数 ∥ v − y ∥ 2 \|v - y\|_2 ∥v−y∥2可以得知，这样 v v v尽可能接近 y y y。
我们将投影表示为
Proj
(
y
;
{
x
1
,
…
x
n
}
)
\operatorname{Proj}\left(y ;\left\{x_{1}, \ldots x_{n}\right\}\right)
Proj(y;{x1,…xn})，并且可以将其正式定义为:
Proj
(
y
;
{
x
1
,
…
x
n
}
)
=
argmin
v
∈
span
(
{
x
1
,
…
,
x
n
}
)
∥
y
−
v
∥
2
\operatorname{Proj}\left(y ;\left\{x_{1}, \ldots x_{n}\right\}\right)=\operatorname{argmin}_{v \in \operatorname{span}\left(\left\{x_{1}, \ldots, x_{n}\right\}\right)}\|y-v\|_{2}
Proj(y;{x1,…xn})=argminv∈span({x1,…,xn})∥y−v∥2
矩阵
A
∈
R
m
×
n
A\in \mathbb{R}^{m \times n}
A∈Rm×n的值域（有时也称为列空间），表示为
R
(
A
)
\mathcal{R}(A)
R(A)，是
A
A
A列的跨度。换句话说，
R
(
A
)
=
{
v
∈
R
m
:
v
=
A
x
,
x
∈
R
n
}
\mathcal{R}(A)=\left\{v \in \mathbb{R}^{m} : v=A x, x \in \mathbb{R}^{n}\right\}
R(A)={v∈Rm:v=Ax,x∈Rn}
做一些技术性的假设（即
A
A
A是满秩且
n
<
m
n <m
n<m），向量
y
∈
R
m
y \in \mathbb{R}^{m}
y∈Rm到
A
A
A的范围的投影由下式给出:
Proj
(
y
;
A
)
=
argmin
v
∈
R
(
A
)
∥
v
−
y
∥
2
=
A
(
A
T
A
)
−
1
A
T
y
\operatorname{Proj}(y ; A)=\operatorname{argmin}_{v \in \mathcal{R}(A)}\|v-y\|_{2}=A\left(A^{T} A\right)^{-1} A^{T} y
Proj(y;A)=argminv∈R(A)∥v−y∥2=A(ATA)−1ATy
这个最后的方程应该看起来非常熟悉，因为它几乎与我们在课程中（我们将很快再次得出）得到的公式：用于参数的最小二乘估计一样。 看一下投影的定义，显而易见，这实际上是我们在最小二乘问题中最小化的目标（除了范数的平方这里有点不一样，这不会影响找到最优解），所以这些问题自然是非常相关的。
当
A
A
A只包含一列时，
a
∈
R
m
a \in \mathbb{R}^{m}
a∈Rm，这给出了向量投影到一条线上的特殊情况：
Proj
(
y
;
a
)
=
a
a
T
a
T
a
y
\operatorname{Proj}(y ; a)=\frac{a a^{T}}{a^{T} a} y
Proj(y;a)=aTaaaTy
一个矩阵
A
∈
R
m
×
n
A\in \mathbb{R}^{m \times n}
A∈Rm×n的零空间
N
(
A
)
\mathcal{N}(A)
N(A) 是所有乘以
A
A
A时等于0向量的集合，即：
N
(
A
)
=
{
x
∈
R
n
:
A
x
=
0
}
\mathcal{N}(A)=\left\{x \in \mathbb{R}^{n} : A x=0\right\}
N(A)={x∈Rn:Ax=0}
注意，
R
(
A
)
\mathcal{R}(A)
R(A)中的向量的大小为
m
m
m，而
N
(
A
)
\mathcal{N}(A)
N(A) 中的向量的大小为
n
n
n，因此
R
(
A
T
)
\mathcal{R}(A^T)
R(AT)和
N
(
A
)
\mathcal{N}(A)
N(A) 中的向量的大小均为
R
n
\mathbb{R}^{n}
Rn。 事实上，还有很多例子。 证明：
{
w
:
w
=
u
+
v
,
u
∈
R
(
A
T
)
,
v
∈
N
(
A
)
}
=
R
n
and
R
(
A
T
)
∩
N
(
A
)
=
{
0
}
\left\{w : w=u+v, u \in \mathcal{R}\left(A^{T}\right), v \in \mathcal{N}(A)\right\}=\mathbb{R}^{n} \text { and } \mathcal{R}\left(A^{T}\right) \cap \mathcal{N}(A)=\{\mathbf{0}\}
{w:w=u+v,u∈R(AT),v∈N(A)}=Rn and R(AT)∩N(A)={0}
换句话说，
R
(
A
T
)
\mathcal{R}(A^T)
R(AT)和
N
(
A
)
\mathcal{N}(A)
N(A) 是不相交的子集，它们一起跨越
R
n
\mathbb{R}^{n}
Rn的整个空间。 这种类型的集合称为正交补，我们用
R
(
A
T
)
=
N
(
A
)
⊥
\mathcal{R}(A^T)= \mathcal{N}(A)^{\perp}
R(AT)=N(A)⊥表示。
3.10 行列式
一个方阵 A ∈ R n × n A \in \mathbb{R}^{n \times n} A∈Rn×n的行列式是函数 det \text {det} det： R n × n → R n \mathbb{R}^{n \times n} \rightarrow \mathbb{R}^{n} Rn×n→Rn，并且表示为 ∣ A ∣ \left| A \right| ∣A∣。 或者 det A \text{det} A detA（有点像迹运算符，我们通常省略括号）。 从代数的角度来说，我们可以写出一个关于 A A A行列式的显式公式。 因此，我们首先提供行列式的几何解释，然后探讨它的一些特定的代数性质。
给定一个矩阵：
[
−
a
1
T
−
−
a
2
T
−
⋮
−
a
n
T
−
]
\left[\begin{array}{cccc}{-} & {a_{1}^{T}} & {-} \\ {-} & {a_{2}^{T}} & {-} \\ {} & {\vdots} & {} \\ {-} & {a_{n}^{T}} & {-}\end{array}\right]
⎣⎢⎢⎢⎡−−−a1Ta2T⋮anT−−−⎦⎥⎥⎥⎤
考虑通过采用
A
A
A行向量
a
1
,
…
a
n
∈
R
n
a_{1}, \ldots a_{n}\in \mathbb{R}^{n}
a1,…an∈Rn的所有可能线性组合形成的点
S
⊂
R
n
S \subset \mathbb{R}^{n}
S⊂Rn的集合，其中线性组合的系数都在0和1之间; 也就是说，集合
S
S
S是
span
(
{
a
1
,
…
a
n
}
)
\text{span}(\{a_{1}, \ldots a_{n}\})
span({a1,…an})受到系数
a
1
,
…
a
n
a_{1}, \ldots a_{n}
a1,…an的限制的线性组合，
α
1
,
⋯
,
α
n
\alpha_1, \cdots ,\alpha_n
α1,⋯,αn满足
0
≤
α
i
≤
1
,
i
=
1
,
…
,
n
0 \leq \alpha_{i} \leq 1, i=1, \ldots, n
0≤αi≤1,i=1,…,n。从形式上看，
S
=
{
v
∈
R
n
:
v
=
∑
i
=
1
n
α
i
a
i
where
0
≤
α
i
≤
1
,
i
=
1
,
…
,
n
}
S=\left\{v \in \mathbb{R}^{n} : v=\sum_{i=1}^{n} \alpha_{i} a_{i} \text { where } 0 \leq \alpha_{i} \leq 1, i=1, \ldots, n\right\}
S={v∈Rn:v=i=1∑nαiai where 0≤αi≤1,i=1,…,n}
事实证明，
A
A
A的行列式的绝对值是对集合
S
S
S的“体积”的度量。
比方说：一个
2
×
2
2 \times2
2×2的矩阵(4)：
A
=
[
1
3
3
2
]
A=\left[\begin{array}{ll}{1} & {3} \\ {3} & {2}\end{array}\right]
A=[1332]
它的矩阵的行是：
a
1
=
[
1
3
]
a
2
=
[
3
2
]
a_{1}=\left[\begin{array}{l}{1} \\ {3}\end{array}\right] \quad a_{2}=\left[\begin{array}{l}{3} \\ {2}\end{array}\right]
a1=[13]a2=[32]
对应于这些行对应的集合
S
S
S如图1所示。对于二维矩阵，
S
S
S通常具有平行四边形的形状。 在我们的例子中，行列式的值是
∣
A
∣
=
−
7
\left| A \right| = -7
∣A∣=−7（可以使用本节后面显示的公式计算），因此平行四边形的面积为7。
在三维中，集合 S S S对应于一个称为平行六面体的对象（一个有倾斜边的三维框，这样每个面都有一个平行四边形）。行定义 S S S的 3 × 3 3×3 3×3矩阵S的行列式的绝对值给出了平行六面体的三维体积。在更高的维度中，集合 S S S是一个称为 n n n维平行切的对象。
在代数上，行列式满足以下三个属性（所有其他属性都遵循这些属性，包括通用公式）：
-
恒等式的行列式为1, ∣ I ∣ = 1 \left| I \right|= 1 ∣I∣=1（几何上，单位超立方体的体积为1）。
-
给定一个矩阵 A ∈ R n × n A \in \mathbb{R}^{n \times n} A∈Rn×n, 如果我们将 A A A中的一行乘上一个标量 t ∈ R t \in \mathbb{R} t∈R，那么新矩阵的行列式是 t ∣ A ∣ t\left| A \right| t∣A∣
∣ [ − t a 1 T − − a 2 T − ⋮ a m T − ] ∣ = t ∣ A ∣ \left|\left[\begin{array}{ccc}{-} & {t a_{1}^{T}} & {-} \\ {-} & {a_{2}^{T}} & {-} \\ {} & {\vdots} & {} \\ {} & {a_{m}^{T}} & {-}\end{array}\right]\right|=t|A| ∣∣∣∣∣∣∣∣∣⎣⎢⎢⎢⎡−−ta1Ta2T⋮amT−−−⎦⎥⎥⎥⎤∣∣∣∣∣∣∣∣∣=t∣A∣
几何上，将集合 S S S的一个边乘以系数 t t t，体积也会增加一个系数 t t t。
- 如果我们交换任意两行在 a i T a_i^T aiT和 a j T a_j^T ajT，那么新矩阵的行列式是 − ∣ A ∣ -\left| A \right| −∣A∣，例如：
∣ [ − a 2 T − − a 1 T − ⋮ − a m T − ] ∣ = − ∣ A ∣ \left|\left[\begin{array}{ccc}{-} & {a_{2}^{T}} & {-} \\ {-} & {a_{1}^{T}} & {-} \\ {} & {\vdots} & {} \\ {-} & {a_{m}^{T}} & {-}\end{array}\right]\right|=-|A| ∣∣∣∣∣∣∣∣∣⎣⎢⎢⎢⎡−−−a2Ta1T⋮amT−−−⎦⎥⎥⎥⎤∣∣∣∣∣∣∣∣∣=−∣A∣
你一定很奇怪，满足上述三个属性的函数的存在并不多。事实上，这样的函数确实存在，而且是唯一的（我们在这里不再证明了）。
从上述三个属性中得出的几个属性包括：
- 对于 A ∈ R n × n A \in \mathbb{R}^{n \times n} A∈Rn×n, ∣ A ∣ = ∣ A T ∣ \left| A \right| = \left| A^T \right| ∣A∣=∣∣AT∣∣
- 对于 A , B ∈ R n × n A,B \in \mathbb{R}^{n \times n} A,B∈Rn×n, ∣ A B ∣ = ∣ A ∣ ∣ B ∣ \left| AB \right|= \left| A \right|\left| B \right| ∣AB∣=∣A∣∣B∣
- 对于 A ∈ R n × n A \in \mathbb{R}^{n \times n} A∈Rn×n, 有且只有当 A A A是奇异的（比如不可逆） ，则： ∣ A ∣ = 0 \left| A \right|= 0 ∣A∣=0
- 对于 A ∈ R n × n A \in \mathbb{R}^{n \times n} A∈Rn×n 同时， A A A为非奇异的，则： ∣ A − 1 ∣ = 1 / ∣ A ∣ \left| A ^{−1}\right| = 1/\left| A \right| ∣∣A−1∣∣=1/∣A∣
在给出行列式的一般定义之前，我们定义，对于
A
∈
R
n
×
n
A \in \mathbb{R}^{n \times n}
A∈Rn×n，
A
\
i
,
\
j
∈
R
(
n
−
1
)
×
(
n
−
1
)
A_{\backslash i, \backslash j}\in \mathbb{R}^{(n-1) \times (n-1)}
A\i,\j∈R(n−1)×(n−1)是由于删除第
i
i
i行和第
j
j
j列而产生的矩阵。 行列式的一般（递归）公式是：
∣
A
∣
=
∑
i
=
1
n
(
−
1
)
i
+
j
a
i
j
∣
A
\
i
,
\
j
∣
(
for any
j
∈
1
,
…
,
n
)
=
∑
j
=
1
n
(
−
1
)
i
+
j
a
i
j
∣
A
\
i
,
\
j
∣
(
for any
i
∈
1
,
…
,
n
)
\begin{aligned}|A| &=\sum_{i=1}^{n}(-1)^{i+j} a_{i j}\left|A_{\backslash i, \backslash j}\right| \quad(\text { for any } j \in 1, \ldots, n) \\ &=\sum_{j=1}^{n}(-1)^{i+j} a_{i j}\left|A_{\backslash i, \backslash j}\right| \quad(\text { for any } i \in 1, \ldots, n) \end{aligned}
∣A∣=i=1∑n(−1)i+jaij∣∣A\i,\j∣∣( for any j∈1,…,n)=j=1∑n(−1)i+jaij∣∣A\i,\j∣∣( for any i∈1,…,n)
对于
A
∈
R
1
×
1
A \in \mathbb{R}^{1 \times 1}
A∈R1×1，初始情况为
∣
A
∣
=
a
11
\left| A \right|= a_{11}
∣A∣=a11。如果我们把这个公式完全展开为
A
∈
R
n
×
n
A \in \mathbb{R}^{n \times n}
A∈Rn×n，就等于
n
!
n!
n!（
n
n
n阶乘）不同的项。因此，对于大于
3
×
3
3×3
3×3的矩阵，我们几乎没有明确地写出完整的行列式方程。然而，
3
×
3
3×3
3×3大小的矩阵的行列式方程是相当常见的，建议好好地了解它们：
∣
[
a
11
]
∣
=
a
11
\left|\left[a_{11}\right]\right|=a_{11}
∣[a11]∣=a11
∣ [ a 11 a 12 a 21 a 22 ] ∣ = a 11 a 22 − a 12 a 21 \left|\left[\begin{array}{ll}{a_{11}} & {a_{12}} \\ {a_{21}} & {a_{22}}\end{array}\right]\right|=a_{11} a_{22}-a_{12} a_{21} ∣∣∣∣[a11a21a12a22]∣∣∣∣=a11a22−a12a21
∣ [ a 11 a 12 a 13 a 21 a 22 a 23 a 31 a 32 a 33 ] ∣ = a 11 a 22 a 33 + a 12 a 23 a 31 + a 13 a 21 a 32 − a 11 a 23 a 32 − a 12 a 21 a 33 − a 13 a 22 a 31 \left|\left[\begin{array}{l}{a_{11}} & {a_{12}} & {a_{13}} \\ {a_{21}} & {a_{22}} & {a_{23}} \\ {a_{31}} & {a_{32}} & {a_{33}}\end{array}\right]\right|=\quad \begin{array}{c}{a_{11} a_{22} a_{33}+a_{12} a_{23} a_{31}+a_{13} a_{21} a_{32}} \\\quad \quad {-a_{11} a_{23} a_{32}-a_{12} a_{21} a_{33}-a_{13} a_{22} a_{31}} \\ {}\end{array} ∣∣∣∣∣∣⎣⎡a11a21a31a12a22a32a13a23a33⎦⎤∣∣∣∣∣∣=a11a22a33+a12a23a31+a13a21a32−a11a23a32−a12a21a33−a13a22a31
矩阵
A
∈
R
n
×
n
A \in \mathbb{R}^{n \times n}
A∈Rn×n的经典伴随矩阵（通常称为伴随矩阵）表示为
adj
(
A
)
\operatorname{adj}(A)
adj(A)，并定义为：
adj
(
A
)
∈
R
n
×
n
,
(
adj
(
A
)
)
i
j
=
(
−
1
)
i
+
j
∣
A
\
j
,
\
i
∣
\operatorname{adj}(A) \in \mathbb{R}^{n \times n}, \quad(\operatorname{adj}(A))_{i j}=(-1)^{i+j}\left|A_{\backslash j, \backslash i}\right|
adj(A)∈Rn×n,(adj(A))ij=(−1)i+j∣∣A\j,\i∣∣
（注意索引
A
\
j
,
\
i
A_{\backslash j, \backslash i}
A\j,\i中的变化）。可以看出，对于任何非奇异
A
∈
R
n
×
n
A \in \mathbb{R}^{n \times n}
A∈Rn×n，
A
−
1
=
1
∣
A
∣
adj
(
A
)
A^{-1}=\frac{1}{|A|} \operatorname{adj}(A)
A−1=∣A∣1adj(A)
虽然这是一个很好的“显式”的逆矩阵公式，但我们应该注意，从数字上讲，有很多更有效的方法来计算逆矩阵。
3.11 二次型和半正定矩阵
给定方矩阵
A
∈
R
n
×
n
A \in \mathbb{R}^{n \times n}
A∈Rn×n和向量
x
∈
R
n
x \in \mathbb{R}^{n}
x∈Rn，标量值
x
T
A
x
x^T Ax
xTAx被称为二次型。 写得清楚些，我们可以看到：
x
T
A
x
=
∑
i
=
1
n
x
i
(
A
x
)
i
=
∑
i
=
1
n
x
i
(
∑
j
=
1
n
A
i
j
x
j
)
=
∑
i
=
1
n
∑
j
=
1
n
A
i
j
x
i
x
j
x^{T} A x=\sum_{i=1}^{n} x_{i}(A x)_{i}=\sum_{i=1}^{n} x_{i}\left(\sum_{j=1}^{n} A_{i j} x_{j}\right)=\sum_{i=1}^{n} \sum_{j=1}^{n} A_{i j} x_{i} x_{j}
xTAx=i=1∑nxi(Ax)i=i=1∑nxi(j=1∑nAijxj)=i=1∑nj=1∑nAijxixj
注意：
x
T
A
x
=
(
x
T
A
x
)
T
=
x
T
A
T
x
=
x
T
(
1
2
A
+
1
2
A
T
)
x
x^{T} A x=\left(x^{T} A x\right)^{T}=x^{T} A^{T} x=x^{T}\left(\frac{1}{2} A+\frac{1}{2} A^{T}\right) x
xTAx=(xTAx)T=xTATx=xT(21A+21AT)x
第一个等号的是因为是标量的转置与自身相等，而第二个等号是因为是我们平均两个本身相等的量。 由此，我们可以得出结论，只有
A
A
A的对称部分有助于形成二次型。 出于这个原因，我们经常隐含地假设以二次型出现的矩阵是对称阵。
我们给出以下定义：
-
对于所有非零向量 x ∈ R n x \in \mathbb{R}^n x∈Rn， x T A x > 0 x^TAx>0 xTAx>0，对称阵 A ∈ S n A \in \mathbb{S}^n A∈Sn为正定（positive definite,PD）。这通常表示为 A ≻ 0 A\succ0 A≻0（或 A > 0 A>0 A>0），并且通常将所有正定矩阵的集合表示为 S + + n \mathbb{S}_{++}^n S++n。
-
对于所有向量 x T A x ≥ 0 x^TAx\geq 0 xTAx≥0，对称矩阵 A ∈ S n A \in \mathbb{S}^n A∈Sn是半正定(positive semidefinite ,PSD)。 这写为（或 A ⪰ 0 A \succeq 0 A⪰0仅 A ≥ 0 A≥0 A≥0），并且所有半正定矩阵的集合通常表示为 S + n \mathbb{S}_+^n S+n。
-
同样，对称矩阵 A ∈ S n A \in \mathbb{S}^n A∈Sn是负定（negative definite,ND），如果对于所有非零 x ∈ R n x \in \mathbb{R}^n x∈Rn，则 x T A x < 0 x^TAx <0 xTAx<0表示为 A ≺ 0 A\prec0 A≺0（或 A < 0 A <0 A<0）。
-
类似地，对称矩阵 A ∈ S n A \in \mathbb{S}^n A∈Sn是半负定(negative semidefinite,NSD），如果对于所有 x ∈ R n x \in \mathbb{R}^n x∈Rn，则 x T A x ≤ 0 x^TAx \leq 0 xTAx≤0表示为 A ⪯ 0 A\preceq 0 A⪯0（或 A ≤ 0 A≤0 A≤0）。
-
最后，对称矩阵 A ∈ S n A \in \mathbb{S}^n A∈Sn是不定的，如果它既不是正半定也不是负半定，即，如果存在 x 1 , x 2 ∈ R n x_1,x_2 \in \mathbb{R}^n x1,x2∈Rn，那么 x 1 T A x 1 > 0 x_1^TAx_1>0 x1TAx1>0且 x 2 T A x 2 < 0 x_2^TAx_2<0 x2TAx2<0。
很明显，如果 A A A是正定的，那么 − A −A −A是负定的，反之亦然。同样，如果 A A A是半正定的，那么 − A −A −A是是半负定的，反之亦然。如果果 A A A是不定的，那么 − A −A −A是也是不定的。
正定矩阵和负定矩阵的一个重要性质是它们总是满秩，因此是可逆的。为了了解这是为什么，假设某个矩阵
A
∈
S
n
A \in \mathbb{S}^n
A∈Sn不是满秩。然后，假设
A
A
A的第
j
j
j列可以表示为其他
n
−
1
n-1
n−1列的线性组合：
a
j
=
∑
i
≠
j
x
i
a
i
a_{j}=\sum_{i \neq j} x_{i} a_{i}
aj=i=j∑xiai
对于某些
x
1
,
⋯
x
j
−
1
,
x
j
+
1
,
⋯
,
x
n
∈
R
x_1,\cdots x_{j-1},x_{j + 1} ,\cdots ,x_n\in \mathbb{R}
x1,⋯xj−1,xj+1,⋯,xn∈R。设
x
j
=
−
1
x_j = -1
xj=−1，则：
A
x
=
∑
i
≠
j
x
i
a
i
=
0
Ax=\sum_{i \neq j} x_{i} a_{i}=0
Ax=i=j∑xiai=0
但这意味着对于某些非零向量
x
x
x，
x
T
A
x
=
0
x^T Ax = 0
xTAx=0，因此
A
A
A必须既不是正定也不是负定。如果
A
A
A是正定或负定，则必须是满秩。
最后，有一种类型的正定矩阵经常出现，因此值得特别提及。 给定矩阵
A
∈
R
m
×
n
A \in \mathbb{R}^{m \times n}
A∈Rm×n（不一定是对称或偶数平方），矩阵
G
=
A
T
A
G = A^T A
G=ATA（有时称为Gram矩阵）总是半正定的。 此外，如果
m
≥
n
m\geq n
m≥n（同时为了方便起见，我们假设
A
A
A是满秩），则
G
=
A
T
A
G = A^T A
G=ATA是正定的。
3.12 特征值和特征向量
给定一个方阵 A ∈ R n × n A \in\mathbb{R}^{n\times n} A∈Rn×n，我们认为在以下条件下， λ ∈ C \lambda \in\mathbb{C} λ∈C是 A A A的特征值， x ∈ C n x\in\mathbb{C}^n x∈Cn是相应的特征向量：
A x = λ x , x ≠ 0 Ax=\lambda x,x \ne 0 Ax=λx,x=0
直观地说，这个定义意味着将 A A A乘以向量 x x x会得到一个新的向量，该向量指向与 x x x相同的方向，但按系数 λ \lambda λ缩放。值得注意的是，对于任何特征向量 x ∈ C n x\in\mathbb{C}^n x∈Cn和标量 t ∈ C t\in\mathbb{C} t∈C， A ( c x ) = c A x = c λ x = λ ( c x ) A(cx)=cAx=c\lambda x=\lambda(cx) A(cx)=cAx=cλx=λ(cx)， c x cx cx也是一个特征向量。因此，当我们讨论与 λ \lambda λ相关的特征向量时，我们通常假设特征向量被标准化为长度为1（这仍然会造成一些歧义，因为 x x x和 − x −x −x都是特征向量，但我们必须接受这一点）。
我们可以重写上面的等式来说明
(
λ
,
x
)
(\lambda,x)
(λ,x)是
A
A
A的特征值和特征向量的组合：
(
λ
I
−
A
)
x
=
0
,
x
≠
0
(\lambda I-A)x=0,x \ne 0
(λI−A)x=0,x=0
但是
(
λ
I
−
A
)
x
=
0
(\lambda I-A)x=0
(λI−A)x=0只有当
(
λ
I
−
A
)
(\lambda I-A)
(λI−A)有一个非空零空间时，同时
(
λ
I
−
A
)
(\lambda I-A)
(λI−A)是奇异的，
x
x
x才具有非零解，即：
∣
(
λ
I
−
A
)
∣
=
0
|(\lambda I-A)|=0
∣(λI−A)∣=0
现在，我们可以使用行列式的先前定义将表达式
∣
(
λ
I
−
A
)
∣
|(\lambda I-A)|
∣(λI−A)∣扩展为
λ
\lambda
λ中的（非常大的）多项式，其中，
λ
\lambda
λ的度为
n
n
n。它通常被称为矩阵
A
A
A的特征多项式。
然后我们找到这个特征多项式的 n n n（可能是复数）根，并用 λ 1 , ⋯ , λ n \lambda_1,\cdots,\lambda_n λ1,⋯,λn表示。这些都是矩阵 A A A的特征值，但我们注意到它们可能不明显。为了找到特征值 λ i \lambda_i λi对应的特征向量，我们只需解线性方程 ( λ I − A ) x = 0 (\lambda I-A)x=0 (λI−A)x=0，因为 ( λ I − A ) (\lambda I-A) (λI−A)是奇异的，所以保证有一个非零解（但也可能有多个或无穷多个解）。
应该注意的是，这不是实际用于数值计算特征值和特征向量的方法（记住行列式的完全展开式有 n ! n! n!项），这是一个数学上的争议。
以下是特征值和特征向量的属性（所有假设在 A ∈ R n × n A \in\mathbb{R}^{n\times n} A∈Rn×n具有特征值 λ 1 , ⋯ , λ n \lambda_1,\cdots,\lambda_n λ1,⋯,λn的前提下）：
-
A A A的迹等于其特征值之和
tr A = ∑ i = 1 n λ i \operatorname{tr} A=\sum_{i=1}^{n} \lambda_{i} trA=i=1∑nλi -
A A A的行列式等于其特征值的乘积
∣ A ∣ = ∏ i = 1 n λ i |A|=\prod_{i=1}^{n} \lambda_{i} ∣A∣=i=1∏nλi -
A A A的秩等于 A A A的非零特征值的个数
-
假设 A A A非奇异，其特征值为 λ \lambda λ和特征向量为 x x x。那么 1 / λ 1/\lambda 1/λ是具有相关特征向量 x x x的 A − 1 A^{-1} A−1的特征值，即 A − 1 x = ( 1 / λ ) x A^{-1}x=(1/\lambda)x A−1x=(1/λ)x。（要证明这一点，取特征向量方程， A x = λ x Ax=\lambda x Ax=λx，两边都左乘 A − 1 A^{-1} A−1）
-
对角阵的特征值 d = d i a g ( d 1 ， ⋯ , d n ) d=diag(d_1，\cdots,d_n) d=diag(d1，⋯,dn)实际上就是对角元素 d 1 ， ⋯ , d n d_1，\cdots,d_n d1，⋯,dn
3.13 对称矩阵的特征值和特征向量
通常情况下，一般的方阵的特征值和特征向量的结构可以很细微地表示出来。
值得庆幸的是，在机器学习的大多数场景下，处理对称实矩阵就足够了，其处理的对称实矩阵的特征值和特征向量具有显着的特性。
在本节中，我们假设 A A A是实对称矩阵, 具有以下属性：
-
A A A的所有特征值都是实数。 我们用用 λ 1 , ⋯ , λ n \lambda_1,\cdots,\lambda_n λ1,⋯,λn表示。
-
存在一组特征向量 u 1 ， ⋯ u n u_1，\cdots u_n u1，⋯un，对于所有 i i i， u i u_i ui是具有特征值 λ i \lambda_{i} λi和 b b b的特征向量。 u 1 ， ⋯ u n u_1，\cdots u_n u1，⋯un是单位向量并且彼此正交。
设
U
U
U是包含
u
i
u_i
ui作为列的正交矩阵：
U
=
[
∣
∣
∣
u
1
u
2
⋯
u
n
∣
∣
∣
]
U=\left[\begin{array}{cccc}{ |} & { |} & {} & { |} \\ {u_{1}} & {u_{2}} & {\cdots} & {u_{n}} \\ { |} & { |} & {} & { |}\end{array}\right]
U=⎣⎡∣u1∣∣u2∣⋯∣un∣⎦⎤
设
Λ
=
d
i
a
g
(
λ
1
,
⋯
,
λ
n
)
\Lambda= diag(\lambda_1,\cdots,\lambda_n)
Λ=diag(λ1,⋯,λn)是包含
λ
1
,
⋯
,
λ
n
\lambda_1,\cdots,\lambda_n
λ1,⋯,λn作为对角线上的元素的对角矩阵。 使用2.3节的方程（2）中的矩阵 - 矩阵向量乘法的方法，我们可以验证：
A
U
=
[
∣
∣
∣
A
u
1
A
u
2
⋯
A
u
n
∣
∣
∣
]
=
[
∣
∣
∣
∣
λ
1
u
1
λ
2
u
2
⋯
λ
n
u
n
∣
∣
∣
∣
]
=
U
diag
(
λ
1
,
…
,
λ
n
)
=
U
Λ
A U=\left[\begin{array}{cccc}{ |} & { |} & {} & { |} \\ {A u_{1}} & {A u_{2}} & {\cdots} & {A u_{n}} \\ { |} & { |} & {} & { |}\end{array}\right]=\left[\begin{array}{ccc}{ |} & { |} & { |} & { |}\\ {\lambda_{1} u_{1}} & {\lambda_{2} u_{2}} & {\cdots} & {\lambda_{n} u_{n}} \\ { |} & { |} & {|} & { |}\end{array}\right]=U \operatorname{diag}\left(\lambda_{1}, \ldots, \lambda_{n}\right)=U \Lambda
AU=⎣⎡∣Au1∣∣Au2∣⋯∣Aun∣⎦⎤=⎣⎡∣λ1u1∣∣λ2u2∣∣⋯∣∣λnun∣⎦⎤=Udiag(λ1,…,λn)=UΛ
考虑到正交矩阵
U
U
U满足
U
U
T
=
I
UU^T=I
UUT=I，利用上面的方程，我们得到：
A
=
A
U
U
T
=
U
Λ
U
T
A=AUU^T=U\Lambda U^T
A=AUUT=UΛUT
这种 A A A的新的表示形式为 U Λ U T U\Lambda U^T UΛUT，通常称为矩阵 A A A的对角化。术语对角化是这样来的：通过这种表示，我们通常可以有效地将对称矩阵 A A A视为对角矩阵 , 这更容易理解。关于由特征向量 U U U定义的基础， 我们将通过几个例子详细说明。
背景知识：代表另一个基的向量。
任何正交矩阵 U = [ ∣ ∣ ∣ u 1 u 2 ⋯ u n ∣ ∣ ∣ ] U=\left[\begin{array}{cccc}{ |} & { |} & {} & { |} \\ {u_{1}} & {u_{2}} & {\cdots} & {u_{n}} \\ { |} & { |} & {} & { |}\end{array}\right] U=⎣⎡∣u1∣∣u2∣⋯∣un∣⎦⎤定义了一个新的属于 R n \mathbb {R}^{n} Rn的基（坐标系），意义如下：对于任何向量 x ∈ R n x \in\mathbb{R}^{n} x∈Rn都可以表示为 u 1 ， ⋯ u n u_1，\cdots u_n u1，⋯un的线性组合，其系数为 x 1 , ⋯ x n x_1,\cdots x_n x1,⋯xn：
x = x ^ 1 u 1 + ⋯ + ⋯ x ^ n u n = U x ^ x=\hat x_1u_1+\cdots +\cdots \hat x_nu_n=U\hat x x=x^1u1+⋯+⋯x^nun=Ux^
在第二个等式中，我们使用矩阵和向量相乘的方法。 实际上，这种
x
^
\hat x
x^是唯一存在的:
x
=
U
x
^
⇔
U
T
x
=
x
^
x=U \hat{x} \Leftrightarrow U^{T} x=\hat{x}
x=Ux^⇔UTx=x^
换句话说，向量
x
^
=
U
T
x
\hat x=U^Tx
x^=UTx可以作为向量
x
x
x的另一种表示，与
U
U
U定义的基有关。
“对角化”矩阵向量乘法。 通过上面的设置，我们将看到左乘矩阵
A
A
A可以被视为左乘以对角矩阵关于特征向量的基。 假设
x
x
x是一个向量，
x
^
\hat x
x^表示
U
U
U的基。设
z
=
A
x
z=Ax
z=Ax为矩阵向量积。现在让我们计算关于
U
U
U的基
z
z
z：
然后，再利用
U
U
T
=
U
T
=
I
UU^T=U^T=I
UUT=UT=I和方程
A
=
A
U
U
T
=
U
Λ
U
T
A=AUU^T=U\Lambda U^T
A=AUUT=UΛUT，我们得到：
z
^
=
U
T
z
=
U
T
A
x
=
U
T
U
Λ
U
T
x
=
Λ
x
^
=
[
λ
1
x
^
1
λ
2
x
^
2
⋮
λ
n
x
^
n
]
\hat{z}=U^{T} z=U^{T} A x=U^{T} U \Lambda U^{T} x=\Lambda \hat{x}=\left[\begin{array}{c}{\lambda_{1} \hat{x}_{1}} \\ {\lambda_{2} \hat{x}_{2}} \\ {\vdots} \\ {\lambda_{n} \hat{x}_{n}}\end{array}\right]
z^=UTz=UTAx=UTUΛUTx=Λx^=⎣⎢⎢⎢⎡λ1x^1λ2x^2⋮λnx^n⎦⎥⎥⎥⎤
我们可以看到，原始空间中的左乘矩阵
A
A
A等于左乘对角矩阵
Λ
\Lambda
Λ相对于新的基，即仅将每个坐标缩放相应的特征值。
在新的基上，矩阵多次相乘也变得简单多了。例如，假设
q
=
A
A
A
x
q=AAAx
q=AAAx。根据
A
A
A的元素导出
q
q
q的分析形式，使用原始的基可能是一场噩梦，但使用新的基就容易多了：
q
^
=
U
T
q
=
U
T
A
A
A
x
=
U
T
U
Λ
U
T
U
Λ
U
T
U
Λ
U
T
x
=
Λ
3
x
^
=
[
λ
1
3
x
^
1
λ
2
3
x
^
2
⋮
λ
n
3
x
^
n
]
\hat{q}=U^{T} q=U^{T} AAA x=U^{T} U \Lambda U^{T} U \Lambda U^{T} U \Lambda U^{T} x=\Lambda^{3} \hat{x}=\left[\begin{array}{c}{\lambda_{1}^{3} \hat{x}_{1}} \\ {\lambda_{2}^{3} \hat{x}_{2}} \\ {\vdots} \\ {\lambda_{n}^{3} \hat{x}_{n}}\end{array}\right]
q^=UTq=UTAAAx=UTUΛUTUΛUTUΛUTx=Λ3x^=⎣⎢⎢⎢⎡λ13x^1λ23x^2⋮λn3x^n⎦⎥⎥⎥⎤
“对角化”二次型。作为直接的推论，二次型
x
T
A
x
x^TAx
xTAx也可以在新的基上简化。
x
T
A
x
=
x
T
U
Λ
U
T
x
=
x
^
Λ
x
^
=
∑
i
=
1
n
λ
i
x
^
i
2
x^{T} A x=x^{T} U \Lambda U^{T} x=\hat{x} \Lambda \hat{x}=\sum_{i=1}^{n} \lambda_{i} \hat{x}_{i}^{2}
xTAx=xTUΛUTx=x^Λx^=i=1∑nλix^i2
(回想一下，在旧的表示法中，
x
T
A
x
=
∑
i
=
1
,
j
=
1
n
x
i
x
j
A
i
j
x^{T} A x=\sum_{i=1, j=1}^{n} x_{i} x_{j} A_{i j}
xTAx=∑i=1,j=1nxixjAij涉及一个
n
2
n^2
n2项的和，而不是上面等式中的
n
n
n项。)利用这个观点，我们还可以证明矩阵
A
A
A的正定性完全取决于其特征值的符号：
- 如果所有的 λ i > 0 \lambda_i>0 λi>0，则矩阵 A A A正定的，因为对于任意的 x ^ ≠ 0 \hat x \ne 0 x^=0, x T A x = ∑ i = 1 n λ i x ^ i 2 > 0 x^{T} A x=\sum_{i=1}^{n} \lambda_{i} \hat{x}_{i}^{2}>0 xTAx=∑i=1nλix^i2>0
- 如果所有的 λ i ≥ 0 \lambda_i\geq 0 λi≥0，则矩阵 A A A是为正半定，因为对于任意的$\hat x , , ,x^{T} A x=\sum_{i=1}^{n} \lambda_{i} \hat{x}_{i}^{2} \geq 0$
- 同样，如果所有 λ i < 0 \lambda_i<0 λi<0或 λ i ≤ 0 \lambda_i\leq 0 λi≤0，则矩阵 A A A分别为负定或半负定。
- 最后，如果 A A A同时具有正特征值和负特征值，比如λ λ i > 0 \lambda_i>0 λi>0和 λ j < 0 \lambda_j<0 λj<0，那么它是不定的。这是因为如果我们让 x ^ \hat x x^满足 x ^ i = 1 \hat x_i=1 x^i=1和 x ^ k = 0 \hat x_k=0 x^k=0，同时所有的 k ≠ i k\ne i k=i，那么 x T A x = ∑ i = 1 n λ i x ^ i 2 > 0 x^{T} A x=\sum_{i=1}^{n} \lambda_{i} \hat{x}_{i}^{2}>0 xTAx=∑i=1nλix^i2>0 ,我们让 x ^ \hat x x^满足 x ^ i = 1 \hat x_i=1 x^i=1和 x ^ k = 0 \hat x_k=0 x^k=0，同时所有的 k ≠ i k\ne i k=i，那么 x T A x = ∑ i = 1 n λ i x ^ i 2 < 0 x^{T} A x=\sum_{i=1}^{n} \lambda_{i} \hat{x}_{i}^{2}<0 xTAx=∑i=1nλix^i2<0
特征值和特征向量经常出现的应用是最大化矩阵的某些函数。特别是对于矩阵
A
∈
S
n
A \in \mathbb{S}^{n}
A∈Sn，考虑以下最大化问题：
max
x
∈
R
n
x
T
A
x
=
∑
i
=
1
n
λ
i
x
^
i
2
subject to
∥
x
∥
2
2
=
1
\max _{x \in \mathbb{R}^{n}} \ x^{T} A x=\sum_{i=1}^{n} \lambda_{i} \hat{x}_{i}^{2} \quad \text { subject to }\|x\|_{2}^{2}=1
x∈Rnmax xTAx=i=1∑nλix^i2 subject to ∥x∥22=1
也就是说，我们要找到（范数1）的向量，它使二次型最大化。假设特征值的阶数为
λ
1
≥
λ
2
≥
⋯
λ
n
\lambda_1 \geq \lambda _2 \geq \cdots \lambda_n
λ1≥λ2≥⋯λn，此优化问题的最优值为
λ
1
\lambda_1
λ1，且与
λ
1
\lambda_1
λ1对应的任何特征向量
u
1
u_1
u1都是最大值之一。（如果
λ
1
>
λ
2
\lambda_1 > \lambda_2
λ1>λ2，那么有一个与特征值
λ
1
\lambda_1
λ1对应的唯一特征向量，它是上面那个优化问题的唯一最大值。）
我们可以通过使用对角化技术来证明这一点：注意，通过公式
∥
U
x
∥
2
=
∥
x
∥
2
\|U x\|_{2}=\|x\|_{2}
∥Ux∥2=∥x∥2推出
∥
x
∥
2
=
∥
x
^
∥
2
\|x\|_{2}=\|\hat{x}\|_{2}
∥x∥2=∥x^∥2，并利用公式：
x
T
A
x
=
x
T
U
Λ
U
T
x
=
x
^
Λ
x
^
=
∑
i
=
1
n
λ
i
x
^
i
2
x^{T} A x=x^{T} U \Lambda U^{T} x=\hat{x} \Lambda \hat{x}=\sum_{i=1}^{n} \lambda_{i} \hat{x}_{i}^{2}
xTAx=xTUΛUTx=x^Λx^=∑i=1nλix^i2，我们可以将上面那个优化问题改写为：
max
x
^
∈
R
n
x
^
T
Λ
x
^
=
∑
i
=
1
n
λ
i
x
^
i
2
subject to
∥
x
^
∥
2
2
=
1
\max _{\hat{x} \in \mathbb{R}^{n}}\ \hat{x}^{T} \Lambda \hat{x}=\sum_{i=1}^{n} \lambda_{i} \hat{x}_{i}^{2} \quad \text { subject to }\|\hat{x}\|_{2}^{2}=1
x^∈Rnmax x^TΛx^=i=1∑nλix^i2 subject to ∥x^∥22=1
然后，我们得到目标的上界为
λ
1
\lambda_1
λ1：
x
^
T
Λ
x
^
=
∑
i
=
1
n
λ
i
x
^
i
2
≤
∑
i
=
1
n
λ
1
x
^
i
2
=
λ
1
\hat{x}^{T} \Lambda \hat{x}=\sum_{i=1}^{n} \lambda_{i} \hat{x}_{i}^{2} \leq \sum_{i=1}^{n} \lambda_{1} \hat{x}_{i}^{2}=\lambda_{1}
x^TΛx^=i=1∑nλix^i2≤i=1∑nλ1x^i2=λ1
此外，设置
x
^
=
[
1
0
⋮
0
]
\hat{x}=\left[\begin{array}{c}{1} \\ {0} \\ {\vdots} \\ {0}\end{array}\right]
x^=⎣⎢⎢⎢⎡10⋮0⎦⎥⎥⎥⎤可让上述等式成立，这与设置
x
=
u
1
x=u_1
x=u1相对应。
4.矩阵微积分
虽然前面章节中的主题通常包含在线性代数的标准课程中，但似乎很少涉及（我们将广泛使用）的一个主题是微积分扩展到向量设置展。尽管我们使用的所有实际微积分都是相对微不足道的，但是符号通常会使事情看起来比实际困难得多。 在本节中，我们将介绍矩阵微积分的一些基本定义，并提供一些示例。
4.1 梯度
假设
f
:
R
m
×
n
→
R
f: \mathbb{R}^{m \times n} \rightarrow \mathbb{R}
f:Rm×n→R是将维度为
m
×
n
m \times n
m×n的矩阵
A
∈
R
m
×
n
A\in \mathbb{R}^{m \times n}
A∈Rm×n作为输入并返回实数值的函数。 然后
f
f
f的梯度（相对于
A
∈
R
m
×
n
A\in \mathbb{R}^{m \times n}
A∈Rm×n）是偏导数矩阵，定义如下：
∇
A
f
(
A
)
∈
R
m
×
n
=
[
∂
f
(
A
)
∂
A
11
∂
f
(
A
)
∂
A
12
⋯
∂
f
(
A
)
∂
A
1
n
∂
f
(
A
)
∂
A
21
∂
f
(
A
)
∂
A
22
⋯
∂
f
(
A
)
∂
A
2
n
⋮
⋮
⋱
⋮
∂
f
(
A
)
∂
A
m
1
∂
f
(
A
)
∂
A
m
2
⋯
∂
f
(
A
)
∂
A
m
n
]
\nabla_{A} f(A) \in \mathbb{R}^{m \times n}=\left[\begin{array}{cccc}{\frac{\partial f(A)}{\partial A_{11}}} & {\frac{\partial f(A)}{\partial A_{12}}} & {\cdots} & {\frac{\partial f(A)}{\partial A_{1n}}} \\ {\frac{\partial f(A)}{\partial A_{21}}} & {\frac{\partial f(A)}{\partial A_{22}}} & {\cdots} & {\frac{\partial f(A)}{\partial A_{2 n}}} \\ {\vdots} & {\vdots} & {\ddots} & {\vdots} \\ {\frac{\partial f(A)}{\partial A_{m 1}}} & {\frac{\partial f(A)}{\partial A_{m 2}}} & {\cdots} & {\frac{\partial f(A)}{\partial A_{m n}}}\end{array}\right]
∇Af(A)∈Rm×n=⎣⎢⎢⎢⎢⎡∂A11∂f(A)∂A21∂f(A)⋮∂Am1∂f(A)∂A12∂f(A)∂A22∂f(A)⋮∂Am2∂f(A)⋯⋯⋱⋯∂A1n∂f(A)∂A2n∂f(A)⋮∂Amn∂f(A)⎦⎥⎥⎥⎥⎤
即，
m
×
n
m \times n
m×n矩阵:
(
∇
A
f
(
A
)
)
i
j
=
∂
f
(
A
)
∂
A
i
j
\left(\nabla_{A} f(A)\right)_{i j}=\frac{\partial f(A)}{\partial A_{i j}}
(∇Af(A))ij=∂Aij∂f(A)
请注意，
∇
A
f
(
A
)
\nabla_{A} f(A)
∇Af(A)的维度始终与
A
A
A的维度相同。特殊情况，如果
A
A
A只是向量
A
∈
R
n
A\in \mathbb{R}^{n}
A∈Rn，则
∇
x
f
(
x
)
=
[
∂
f
(
x
)
∂
x
1
∂
f
(
x
)
∂
x
2
⋮
∂
f
(
x
)
∂
x
n
]
\nabla_{x} f(x)=\left[\begin{array}{c}{\frac{\partial f(x)}{\partial x_{1}}} \\ {\frac{\partial f(x)}{\partial x_{2}}} \\ {\vdots} \\ {\frac{\partial f(x)}{\partial x_{n}}}\end{array}\right]
∇xf(x)=⎣⎢⎢⎢⎢⎡∂x1∂f(x)∂x2∂f(x)⋮∂xn∂f(x)⎦⎥⎥⎥⎥⎤
重要的是要记住，只有当函数是实值时，即如果函数返回标量值，才定义函数的梯度。例如，
A
∈
R
m
×
n
A\in \mathbb{R}^{m \times n}
A∈Rm×n相对于
x
x
x，我们不能取
A
x
Ax
Ax的梯度，因为这个量是向量值。
它直接从偏导数的等价性质得出：
-
∇ x ( f ( x ) + g ( x ) ) = ∇ x f ( x ) + ∇ x g ( x ) \nabla_{x}(f(x)+g(x))=\nabla_{x} f(x)+\nabla_{x} g(x) ∇x(f(x)+g(x))=∇xf(x)+∇xg(x)
-
对于 t ∈ R t \in \mathbb{R} t∈R ， ∇ x ( t f ( x ) ) = t ∇ x f ( x ) \nabla_{x}(t f(x))=t \nabla_{x} f(x) ∇x(tf(x))=t∇xf(x)
原则上，梯度是偏导数对多变量函数的自然延伸。然而，在实践中，由于符号的原因，使用梯度有时是很困难的。例如，假设
A
∈
R
m
×
n
A\in \mathbb{R}^{m \times n}
A∈Rm×n是一个固定系数矩阵，假设
b
∈
R
m
b\in \mathbb{R}^{m}
b∈Rm是一个固定系数向量。设
f
:
R
m
×
n
→
R
f: \mathbb{R}^{m \times n} \rightarrow \mathbb{R}
f:Rm×n→R为
f
(
z
)
=
z
T
z
f(z)=z^Tz
f(z)=zTz定义的函数，因此
∇
z
f
(
z
)
=
2
z
\nabla_{z}f(z)=2z
∇zf(z)=2z。但现在考虑表达式，
∇
f
(
A
x
)
\nabla f(Ax)
∇f(Ax)
该表达式应该如何解释？ 至少有两种可能性：
1.在第一个解释中，回想起
∇
z
f
(
z
)
=
2
z
\nabla_{z}f(z)=2z
∇zf(z)=2z。 在这里，我们将
∇
f
(
A
x
)
\nabla f(Ax)
∇f(Ax)解释为评估点
A
x
Ax
Ax处的梯度，因此:
∇
f
(
A
x
)
=
2
(
A
x
)
=
2
A
x
∈
R
m
\nabla f(A x)=2(A x)=2 A x \in \mathbb{R}^{m}
∇f(Ax)=2(Ax)=2Ax∈Rm
2.在第二种解释中，我们将数量
f
(
A
x
)
f(Ax)
f(Ax)视为输入变量
x
x
x的函数。 更正式地说，设
g
(
x
)
=
f
(
A
x
)
g(x) =f(Ax)
g(x)=f(Ax)。 然后在这个解释中:
∇
f
(
A
x
)
=
∇
x
g
(
x
)
∈
R
n
\nabla f(A x)=\nabla_{x} g(x) \in \mathbb{R}^{n}
∇f(Ax)=∇xg(x)∈Rn
在这里，我们可以看到这两种解释确实不同。 一种解释产生 m m m维向量作为结果，而另一种解释产生 n n n维向量作为结果！ 我们怎么解决这个问题？
这里，关键是要明确我们要区分的变量。
在第一种情况下，我们将函数
f
f
f与其参数
z
z
z进行区分，然后替换参数
A
x
Ax
Ax。
在第二种情况下，我们将复合函数
g
(
x
)
=
f
(
A
x
)
g(x)=f(Ax)
g(x)=f(Ax)直接与
x
x
x进行微分。
我们将第一种情况表示为 ∇ z f ( A x ) \nabla zf(Ax) ∇zf(Ax)，第二种情况表示为 ∇ x f ( A x ) \nabla xf(Ax) ∇xf(Ax)。
保持符号清晰是非常重要的，以后完成课程作业时候你就会发现。
4.2 黑塞矩阵
假设
f
:
R
n
→
R
f: \mathbb{R}^{n} \rightarrow \mathbb{R}
f:Rn→R是一个函数，它接受
R
n
\mathbb{R}^{n}
Rn中的向量并返回实数。那么关于
x
x
x的黑塞矩阵（也有翻译作海森矩阵），写做：
∇
x
2
f
(
A
x
)
\nabla_x ^2 f(A x)
∇x2f(Ax)，或者简单地说，
H
H
H是
n
×
n
n \times n
n×n矩阵的偏导数：
∇
x
2
f
(
x
)
∈
R
n
×
n
=
[
∂
2
f
(
x
)
∂
x
1
2
∂
2
f
(
x
)
∂
x
1
∂
x
2
⋯
∂
2
f
(
x
)
∂
x
1
∂
x
n
∂
2
f
(
x
)
∂
x
2
∂
x
1
∂
2
f
(
x
)
∂
x
2
2
⋯
∂
2
f
(
x
)
∂
x
2
∂
x
n
⋮
⋮
⋱
⋮
∂
2
f
(
x
)
∂
x
n
∂
x
1
∂
2
f
(
x
)
∂
x
n
∂
x
2
⋯
∂
2
f
(
x
)
∂
x
n
2
]
\nabla_{x}^{2} f(x) \in \mathbb{R}^{n \times n}=\left[\begin{array}{cccc}{\frac{\partial^{2} f(x)}{\partial x_{1}^{2}}} & {\frac{\partial^{2} f(x)}{\partial x_{1} \partial x_{2}}} & {\cdots} & {\frac{\partial^{2} f(x)}{\partial x_{1} \partial x_{n}}} \\ {\frac{\partial^{2} f(x)}{\partial x_{2} \partial x_{1}}} & {\frac{\partial^{2} f(x)}{\partial x_{2}^{2}}} & {\cdots} & {\frac{\partial^{2} f(x)}{\partial x_{2} \partial x_{n}}} \\ {\vdots} & {\vdots} & {\ddots} & {\vdots} \\ {\frac{\partial^{2} f(x)}{\partial x_{n} \partial x_{1}}} & {\frac{\partial^{2} f(x)}{\partial x_{n} \partial x_{2}}} & {\cdots} & {\frac{\partial^{2} f(x)}{\partial x_{n}^{2}}}\end{array}\right]
∇x2f(x)∈Rn×n=⎣⎢⎢⎢⎢⎢⎡∂x12∂2f(x)∂x2∂x1∂2f(x)⋮∂xn∂x1∂2f(x)∂x1∂x2∂2f(x)∂x22∂2f(x)⋮∂xn∂x2∂2f(x)⋯⋯⋱⋯∂x1∂xn∂2f(x)∂x2∂xn∂2f(x)⋮∂xn2∂2f(x)⎦⎥⎥⎥⎥⎥⎤
换句话说，
∇
x
2
f
(
x
)
∈
R
n
×
n
\nabla_{x}^{2} f(x) \in \mathbb{R}^{n \times n}
∇x2f(x)∈Rn×n，其：
(
∇
x
2
f
(
x
)
)
i
j
=
∂
2
f
(
x
)
∂
x
i
∂
x
j
\left(\nabla_{x}^{2} f(x)\right)_{i j}=\frac{\partial^{2} f(x)}{\partial x_{i} \partial x_{j}}
(∇x2f(x))ij=∂xi∂xj∂2f(x)
注意：黑塞矩阵通常是对称阵：
∂
2
f
(
x
)
∂
x
i
∂
x
j
=
∂
2
f
(
x
)
∂
x
j
∂
x
i
\frac{\partial^{2} f(x)}{\partial x_{i} \partial x_{j}}=\frac{\partial^{2} f(x)}{\partial x_{j} \partial x_{i}}
∂xi∂xj∂2f(x)=∂xj∂xi∂2f(x)
与梯度相似，只有当
f
(
x
)
f(x)
f(x)为实值时才定义黑塞矩阵。
很自然地认为梯度与向量函数的一阶导数的相似，而黑塞矩阵与二阶导数的相似（我们使用的符号也暗示了这种关系）。 这种直觉通常是正确的，但需要记住以下几个注意事项。
首先，对于一个变量
f
:
R
→
R
f: \mathbb{R} \rightarrow \mathbb{R}
f:R→R的实值函数，它的基本定义：二阶导数是一阶导数的导数，即：
∂
2
f
(
x
)
∂
x
2
=
∂
∂
x
∂
∂
x
f
(
x
)
\frac{\partial^{2} f(x)}{\partial x^{2}}=\frac{\partial}{\partial x} \frac{\partial}{\partial x} f(x)
∂x2∂2f(x)=∂x∂∂x∂f(x)
然而，对于向量的函数，函数的梯度是一个向量，我们不能取向量的梯度，即:
∇
x
∇
x
f
(
x
)
=
∇
x
[
∂
f
(
x
)
∂
x
1
∂
f
(
x
)
∂
x
2
⋮
∂
f
(
x
)
∂
x
n
]
\nabla_{x} \nabla_{x} f(x)=\nabla_{x}\left[\begin{array}{c}{\frac{\partial f(x)}{\partial x_{1}}} \\ {\frac{\partial f(x)}{\partial x_{2}}} \\ {\vdots} \\ {\frac{\partial f(x)}{\partial x_{n}}}\end{array}\right]
∇x∇xf(x)=∇x⎣⎢⎢⎢⎢⎡∂x1∂f(x)∂x2∂f(x)⋮∂xn∂f(x)⎦⎥⎥⎥⎥⎤
上面这个表达式没有意义。 因此，黑塞矩阵不是梯度的梯度。 然而，下面这种情况却这几乎是正确的：如果我们看一下梯度
(
∇
x
f
(
x
)
)
i
=
∂
f
(
x
)
/
∂
x
i
\left(\nabla_{x} f(x)\right)_{i}=\partial f(x) / \partial x_{i}
(∇xf(x))i=∂f(x)/∂xi的第
i
i
i个元素，并取关于于
x
x
x的梯度我们得到：
∇
x
∂
f
(
x
)
∂
x
i
=
[
∂
2
f
(
x
)
∂
x
i
∂
x
1
∂
2
f
(
x
)
∂
x
2
∂
x
2
⋮
∂
f
(
x
)
∂
x
i
∂
x
n
]
\nabla_{x} \frac{\partial f(x)}{\partial x_{i}}=\left[\begin{array}{c}{\frac{\partial^{2} f(x)}{\partial x_{i} \partial x_{1}}} \\ {\frac{\partial^{2} f(x)}{\partial x_{2} \partial x_{2}}} \\ {\vdots} \\ {\frac{\partial f(x)}{\partial x_{i} \partial x_{n}}}\end{array}\right]
∇x∂xi∂f(x)=⎣⎢⎢⎢⎢⎡∂xi∂x1∂2f(x)∂x2∂x2∂2f(x)⋮∂xi∂xn∂f(x)⎦⎥⎥⎥⎥⎤
这是黑塞矩阵第
i
i
i行（列）,所以：
∇
x
2
f
(
x
)
=
[
∇
x
(
∇
x
f
(
x
)
)
1
∇
x
(
∇
x
f
(
x
)
)
2
⋯
∇
x
(
∇
x
f
(
x
)
)
n
]
\nabla_{x}^{2} f(x)=\left[\nabla_{x}\left(\nabla_{x} f(x)\right)_{1} \quad \nabla_{x}\left(\nabla_{x} f(x)\right)_{2} \quad \cdots \quad \nabla_{x}\left(\nabla_{x} f(x)\right)_{n}\right]
∇x2f(x)=[∇x(∇xf(x))1∇x(∇xf(x))2⋯∇x(∇xf(x))n]
简单地说：我们可以说由于：
∇
x
2
f
(
x
)
=
∇
x
(
∇
x
f
(
x
)
)
T
\nabla_{x}^{2} f(x)=\nabla_{x}\left(\nabla_{x} f(x)\right)^{T}
∇x2f(x)=∇x(∇xf(x))T，只要我们理解，这实际上是取
∇
x
f
(
x
)
\nabla_{x} f(x)
∇xf(x)的每个元素的梯度，而不是整个向量的梯度。
最后，请注意，虽然我们可以对矩阵
A
∈
R
n
A\in \mathbb{R}^{n}
A∈Rn取梯度，但对于这门课，我们只考虑对向量
x
∈
R
n
x \in \mathbb{R}^{n}
x∈Rn取黑塞矩阵。
这会方便很多（事实上，我们所做的任何计算都不要求我们找到关于矩阵的黑森方程），因为关于矩阵的黑塞方程就必须对矩阵所有元素求偏导数
∂
2
f
(
A
)
/
(
∂
A
i
j
∂
A
k
ℓ
)
\partial^{2} f(A) /\left(\partial A_{i j} \partial A_{k \ell}\right)
∂2f(A)/(∂Aij∂Akℓ)，将其表示为矩阵相当麻烦。
4.3 二次函数和线性函数的梯度和黑塞矩阵
现在让我们尝试确定几个简单函数的梯度和黑塞矩阵。 应该注意的是，这里给出的所有梯度都是梯度的特殊情况。
对于 x ∈ R n x \in \mathbb{R}^{n} x∈Rn, 设 f ( x ) = b T x f(x)=b^Tx f(x)=bTx 的某些已知向量 b ∈ R n b \in \mathbb{R}^{n} b∈Rn ，则：
f
(
x
)
=
∑
i
=
1
n
b
i
x
i
f(x)=\sum_{i=1}^{n} b_{i} x_{i}
f(x)=i=1∑nbixi
所以：
∂
f
(
x
)
∂
x
k
=
∂
∂
x
k
∑
i
=
1
n
b
i
x
i
=
b
k
\frac{\partial f(x)}{\partial x_{k}}=\frac{\partial}{\partial x_{k}} \sum_{i=1}^{n} b_{i} x_{i}=b_{k}
∂xk∂f(x)=∂xk∂i=1∑nbixi=bk
由此我们可以很容易地看出
∇
x
b
T
x
=
b
\nabla_{x} b^{T} x=b
∇xbTx=b。 这应该与单变量微积分中的类似情况进行比较，其中
∂
/
(
∂
x
)
a
x
=
a
\partial /(\partial x) a x=a
∂/(∂x)ax=a。
现在考虑
A
∈
S
n
A\in \mathbb{S}^{n}
A∈Sn的二次函数
f
(
x
)
=
x
T
A
x
f(x)=x^TAx
f(x)=xTAx。 记住这一点：
f
(
x
)
=
∑
i
=
1
n
∑
j
=
1
n
A
i
j
x
i
x
j
f(x)=\sum_{i=1}^{n} \sum_{j=1}^{n} A_{i j} x_{i} x_{j}
f(x)=i=1∑nj=1∑nAijxixj
为了取偏导数，我们将分别考虑包括
x
k
x_k
xk和
x
2
k
x_2^k
x2k因子的项：
∂
f
(
x
)
∂
x
k
=
∂
∂
x
k
∑
i
=
1
n
∑
j
=
1
n
A
i
j
x
i
x
j
=
∂
∂
x
k
[
∑
i
≠
k
∑
j
≠
k
A
i
j
x
i
x
j
+
∑
i
≠
k
A
i
k
x
i
x
k
+
∑
j
≠
k
A
k
j
x
k
x
j
+
A
k
k
x
k
2
]
=
∑
i
≠
k
A
i
k
x
i
+
∑
j
≠
k
A
k
j
x
j
+
2
A
k
k
x
k
=
∑
i
=
1
n
A
i
k
x
i
+
∑
j
=
1
n
A
k
j
x
j
=
2
∑
i
=
1
n
A
k
i
x
i
\begin{aligned} \frac{\partial f(x)}{\partial x_{k}} &=\frac{\partial}{\partial x_{k}} \sum_{i=1}^{n} \sum_{j=1}^{n} A_{i j} x_{i} x_{j} \\ &=\frac{\partial}{\partial x_{k}}\left[\sum_{i \neq k} \sum_{j \neq k} A_{i j} x_{i} x_{j}+\sum_{i \neq k} A_{i k} x_{i} x_{k}+\sum_{j \neq k} A_{k j} x_{k} x_{j}+A_{k k} x_{k}^{2}\right] \\ &=\sum_{i \neq k} A_{i k} x_{i}+\sum_{j \neq k} A_{k j} x_{j}+2 A_{k k} x_{k} \\ &=\sum_{i=1}^{n} A_{i k} x_{i}+\sum_{j=1}^{n} A_{k j} x_{j}=2 \sum_{i=1}^{n} A_{k i} x_{i} \end{aligned}
∂xk∂f(x)=∂xk∂i=1∑nj=1∑nAijxixj=∂xk∂⎣⎡i=k∑j=k∑Aijxixj+i=k∑Aikxixk+j=k∑Akjxkxj+Akkxk2⎦⎤=i=k∑Aikxi+j=k∑Akjxj+2Akkxk=i=1∑nAikxi+j=1∑nAkjxj=2i=1∑nAkixi
最后一个等式，是因为
A
A
A是对称的（我们可以安全地假设，因为它以二次形式出现）。 注意，
∇
x
f
(
x
)
\nabla_{x} f(x)
∇xf(x)的第
k
k
k个元素是
A
A
A和
x
x
x的第
k
k
k行的内积。 因此，
∇
x
x
T
A
x
=
2
A
x
\nabla_{x} x^{T} A x=2 A x
∇xxTAx=2Ax。 同样，这应该提醒你单变量微积分中的类似事实，即
∂
/
(
∂
x
)
a
x
2
=
2
a
x
\partial /(\partial x) a x^{2}=2 a x
∂/(∂x)ax2=2ax。
最后，让我们来看看二次函数
f
(
x
)
=
x
T
A
x
f(x)=x^TAx
f(x)=xTAx黑塞矩阵（显然，线性函数
b
T
x
b^Tx
bTx的黑塞矩阵为零）。在这种情况下:
∂
2
f
(
x
)
∂
x
k
∂
x
ℓ
=
∂
∂
x
k
[
∂
f
(
x
)
∂
x
ℓ
]
=
∂
∂
x
k
[
2
∑
i
=
1
n
A
ℓ
i
x
i
]
=
2
A
ℓ
k
=
2
A
k
ℓ
\frac{\partial^{2} f(x)}{\partial x_{k} \partial x_{\ell}}=\frac{\partial}{\partial x_{k}}\left[\frac{\partial f(x)}{\partial x_{\ell}}\right]=\frac{\partial}{\partial x_{k}}\left[2 \sum_{i=1}^{n} A_{\ell i} x_{i}\right]=2 A_{\ell k}=2 A_{k \ell}
∂xk∂xℓ∂2f(x)=∂xk∂[∂xℓ∂f(x)]=∂xk∂[2i=1∑nAℓixi]=2Aℓk=2Akℓ
因此，应该很清楚
∇
x
2
x
T
A
x
=
2
A
\nabla_{x}^2 x^{T} A x=2 A
∇x2xTAx=2A，这应该是完全可以理解的（同样类似于
∂
2
/
(
∂
x
2
)
a
x
2
=
2
a
\partial^2 /(\partial x^2) a x^{2}=2a
∂2/(∂x2)ax2=2a的单变量事实）。
简要概括起来：
-
∇ x b T x = b \nabla_{x} b^{T} x=b ∇xbTx=b
-
∇ x x T A x = 2 A x \nabla_{x} x^{T} A x=2 A x ∇xxTAx=2Ax (如果 A A A是对称阵)
-
∇ x 2 x T A x = 2 A \nabla_{x}^2 x^{T} A x=2 A ∇x2xTAx=2A (如果 A A A是对称阵)
4.4 最小二乘法
让我们应用上一节中得到的方程来推导最小二乘方程。假设我们得到矩阵 A ∈ R m × n A\in \mathbb{R}^{m \times n} A∈Rm×n（为了简单起见，我们假设 A A A是满秩）和向量 b ∈ R m b\in \mathbb{R}^{m} b∈Rm，从而使 b ∉ R ( A ) b \notin \mathcal{R}(A) b∈/R(A)。在这种情况下，我们将无法找到向量 x ∈ R n x\in \mathbb{R}^{n} x∈Rn，由于 A x = b Ax = b Ax=b，因此我们想要找到一个向量 x x x，使得 A x Ax Ax尽可能接近 b b b，用欧几里德范数的平方$|A x-b|_{2}^{2} $来衡量。
使用公式 ∥ x ∥ 2 = x T x \|x\|^{2}=x^Tx ∥x∥2=xTx，我们可以得到：
∥
A
x
−
b
∥
2
2
=
(
A
x
−
b
)
T
(
A
x
−
b
)
=
x
T
A
T
A
x
−
2
b
T
A
x
+
b
T
b
\begin{aligned}\|A x-b\|_{2}^{2} &=(A x-b)^{T}(A x-b) \\ &=x^{T} A^{T} A x-2 b^{T} A x+b^{T} b \end{aligned}
∥Ax−b∥22=(Ax−b)T(Ax−b)=xTATAx−2bTAx+bTb
根据
x
x
x的梯度，并利用上一节中推导的性质：
∇
x
(
x
T
A
T
A
x
−
2
b
T
A
x
+
b
T
b
)
=
∇
x
x
T
A
T
A
x
−
∇
x
2
b
T
A
x
+
∇
x
b
T
b
=
2
A
T
A
x
−
2
A
T
b
\begin{aligned} \nabla_{x}\left(x^{T} A^{T} A x-2 b^{T} A x+b^{T} b\right) &=\nabla_{x} x^{T} A^{T} A x-\nabla_{x} 2 b^{T} A x+\nabla_{x} b^{T} b \\ &=2 A^{T} A x-2 A^{T} b \end{aligned}
∇x(xTATAx−2bTAx+bTb)=∇xxTATAx−∇x2bTAx+∇xbTb=2ATAx−2ATb
将最后一个表达式设置为零，然后解出
x
x
x，得到了正规方程：
x
=
(
A
T
A
)
−
1
A
T
b
x = (A^TA)^{-1}A^Tb
x=(ATA)−1ATb
这和我们在课堂上得到的相同。
4.5 行列式的梯度
现在让我们考虑一种情况，我们找到一个函数相对于矩阵的梯度，也就是说，对于
A
∈
R
n
×
n
A\in \mathbb{R}^{n \times n}
A∈Rn×n，我们要找到
∇
A
∣
A
∣
\nabla_{A}|A|
∇A∣A∣。回想一下我们对行列式的讨论：
∣
A
∣
=
∑
i
=
1
n
(
−
1
)
i
+
j
A
i
j
∣
A
\
i
,
\
j
∣
(
for any
j
∈
1
,
…
,
n
)
|A|=\sum_{i=1}^{n}(-1)^{i+j} A_{i j}\left|A_{\backslash i, \backslash j}\right| \quad(\text { for any } j \in 1, \ldots, n)
∣A∣=i=1∑n(−1)i+jAij∣∣A\i,\j∣∣( for any j∈1,…,n)
所以：
∂
∂
A
k
ℓ
∣
A
∣
=
∂
∂
A
k
ℓ
∑
i
=
1
n
(
−
1
)
i
+
j
A
i
j
∣
A
\
i
,
\
j
∣
=
(
−
1
)
k
+
ℓ
∣
A
\
k
,
\
ℓ
∣
=
(
adj
(
A
)
)
ℓ
k
\frac{\partial}{\partial A_{k \ell}}|A|=\frac{\partial}{\partial A_{k \ell}} \sum_{i=1}^{n}(-1)^{i+j} A_{i j}\left|A_{\backslash i, \backslash j}\right|=(-1)^{k+\ell}\left|A_{\backslash k,\backslash \ell}\right|=(\operatorname{adj}(A))_{\ell k}
∂Akℓ∂∣A∣=∂Akℓ∂i=1∑n(−1)i+jAij∣∣A\i,\j∣∣=(−1)k+ℓ∣∣A\k,\ℓ∣∣=(adj(A))ℓk
从这里可以知道，它直接从伴随矩阵的性质得出：
∇
A
∣
A
∣
=
(
adj
(
A
)
)
T
=
∣
A
∣
A
−
T
\nabla_{A}|A|=(\operatorname{adj}(A))^{T}=|A| A^{-T}
∇A∣A∣=(adj(A))T=∣A∣A−T
现在我们来考虑函数
f
:
S
+
+
n
→
R
f : \mathbb{S}_{++}^{n} \rightarrow \mathbb{R}
f:S++n→R，
f
(
A
)
=
log
∣
A
∣
f(A)=\log |A|
f(A)=log∣A∣。注意，我们必须将
f
f
f的域限制为正定矩阵，因为这确保了
∣
A
∣
>
0
|A|>0
∣A∣>0，因此
∣
A
∣
|A|
∣A∣的对数是实数。在这种情况下，我们可以使用链式法则（没什么奇怪的，只是单变量演算中的普通链式法则）来看看：
∂
log
∣
A
∣
∂
A
i
j
=
∂
log
∣
A
∣
∂
∣
A
∣
∂
∣
A
∣
∂
A
i
j
=
1
∣
A
∣
∂
∣
A
∣
∂
A
i
j
\frac{\partial \log |A|}{\partial A_{i j}}=\frac{\partial \log |A|}{\partial|A|} \frac{\partial|A|}{\partial A_{i j}}=\frac{1}{|A|} \frac{\partial|A|}{\partial A_{i j}}
∂Aij∂log∣A∣=∂∣A∣∂log∣A∣∂Aij∂∣A∣=∣A∣1∂Aij∂∣A∣
从这一点可以明显看出：
∇
A
log
∣
A
∣
=
1
∣
A
∣
∇
A
∣
A
∣
=
A
−
1
\nabla_{A} \log |A|=\frac{1}{|A|} \nabla_{A}|A|=A^{-1}
∇Alog∣A∣=∣A∣1∇A∣A∣=A−1
我们可以在最后一个表达式中删除转置，因为
A
A
A是对称的。注意与单值情况的相似性，其中
∂
/
(
∂
x
)
log
x
=
1
/
x
\partial /(\partial x) \log x=1 / x
∂/(∂x)logx=1/x。
4.6 特征值优化
最后，我们使用矩阵演算以直接导致特征值/特征向量分析的方式求解优化问题。 考虑以下等式约束优化问题：
max
x
∈
R
n
x
T
A
x
subject to
∥
x
∥
2
2
=
1
\max _{x \in \mathbb{R}^{n}} x^{T} A x \quad \text { subject to }\|x\|_{2}^{2}=1
x∈RnmaxxTAx subject to ∥x∥22=1
对于对称矩阵
A
∈
S
n
A\in \mathbb{S}^{n}
A∈Sn。求解等式约束优化问题的标准方法是采用拉格朗日形式，一种包含等式约束的目标函数，在这种情况下，拉格朗日函数可由以下公式给出：
L
(
x
,
λ
)
=
x
T
A
x
−
λ
x
T
x
\mathcal{L}(x, \lambda)=x^{T} A x-\lambda x^{T} x
L(x,λ)=xTAx−λxTx
其中，
λ
\lambda
λ被称为与等式约束关联的拉格朗日乘子。可以确定，要使
x
∗
x^*
x∗成为问题的最佳点，拉格朗日的梯度必须在
x
∗
x^*
x∗处为零（这不是唯一的条件，但它是必需的）。也就是说，
∇
x
L
(
x
,
λ
)
=
∇
x
(
x
T
A
x
−
λ
x
T
x
)
=
2
A
T
x
−
2
λ
x
=
0
\nabla_{x} \mathcal{L}(x, \lambda)=\nabla_{x}\left(x^{T} A x-\lambda x^{T} x\right)=2 A^{T} x-2 \lambda x=0
∇xL(x,λ)=∇x(xTAx−λxTx)=2ATx−2λx=0
请注意，这只是线性方程
A
x
=
λ
x
Ax =\lambda x
Ax=λx。 这表明假设
x
T
x
=
1
x^T x = 1
xTx=1，可能最大化（或最小化）
x
T
A
x
x^T Ax
xTAx的唯一点是
A
A
A的特征向量。