向量 与Ｐｙｔｈｏｎ

前言

在 上文 中我简单概括了矩阵的基本运算，并给出了两个应用实例。这篇文章我们继续谈谈向量。

向量是线性代数中的基本概念，也是机器学习的基础数据表示形式。例如计算机阅读文本的过程首先就会将文本分词，然后用向量表示。这是因为向量很适合在高维空间中表达和处理。在机器学习中会接触到的诸如投影、降维的概念，都是在向量的基础上做的。

在 Rn 空间中定义的向量 V⃗  ，可以用一个包含 n 个实数的有序集来表示，即 V⃗ =⎡⎣⎢⎢⎢v1v2…vn⎤⎦⎥⎥⎥ ，这个有序集里的每个元素称为向量的 分量 。例如一个 R2 空间中的向量 [21] ，有些地方也会用 (2,1) 或 <2,1> 这样的形式来表示。

绘图表示这个变量：

向量的长度被定义为

∥∥v⃗ ∥∥=v21+v22+…+v2n−−−−−−−−−−−−−−√

，和我们以往所接触的距离公式一模一样。长度为 1 的向量称为  单位向量 。

基本运算

加

向量 a 与向量 b 的加法定义为:

a+b=⎡⎣⎢⎢⎢a1+b1a2+b2…an+bn⎤⎦⎥⎥⎥

绘图示意向量 a=[−12] 与 b=[31] 的相加，值为 [23] ：

在 Python 中，可以直接用 Numpy 的 ndarray 来表示向量。

import numpy as np
a = np.array([-1, 2])
b = np.array([3, 1])
print a + b # [2 3]

减

a−b=⎡⎣⎢⎢⎢a1−b1a2−b2…an−bn⎤⎦⎥⎥⎥

从几何角度讲，向量减相当于加上一个反向的向量。

import numpy as np
a = np.array([-1, 2])
b = np.array([3, 1])
print a - b  # [-4, 1]

乘

标量乘向量

标量 c 乘以向量 a 定义为：

c⋅a=⎡⎣⎢⎢⎢c⋅a1c⋅a2…c⋅an⎤⎦⎥⎥⎥=⎡⎣⎢⎢⎢a1⋅ca2⋅c…an⋅c⎤⎦⎥⎥⎥

绘图示意向量 a=[−12] 乘以一个标量 3 得到 [−36] ：

Python 实现：

import numpy as np
a = np.array([-1, 2])
print a * 3 #[-3, 6]

向量点积

向量的点积（又叫点乘）定义如下：

a⃗ ⋅b⃗ =⎡⎣⎢⎢⎢a1a2…an⎤⎦⎥⎥⎥⋅⎡⎣⎢⎢⎢b1b2…bn⎤⎦⎥⎥⎥=a1b1+a2b2+…+anbn

可见点积得到的是一个标量。

例如：

⎡⎣⎢352⎤⎦⎥⋅⎡⎣⎢147⎤⎦⎥=3⋅1+5⋅4+2⋅7=37

Python 示例：

import numpy as np
a = np.array([3, 5, 2])
b = np.array([1, 4, 7])
print a.dot(b)  # 37
print np.dot(a, b)  # 37（另一种等价写法）

容易证明点积满足乘法交换律、分配律和结合律。

我们知道向量的长度定义为 ∥∥v⃗ ∥∥=v21+v22+…+v2n−−−−−−−−−−−−−−√ ，联立点积的定义，可以得出：

eq: 1 »

∥∥v⃗ ∥∥=v21+v22+…+v2n−−−−−−−−−−−−−−√=v⃗ ⋅v⃗ −−−−√

关于点积还有一个非常重要的性质，称为 柯西不等式 ：

• 对两个非 0 向量 x⃗ ,y⃗ ∈Rn ， |x⃗ ⋅y⃗ |≤∥∥x⃗ ∥∥∥∥y⃗ ∥∥ 。
• 当且仅当 x⃗ =cy⃗  时，等式成立。

虽然受限于篇幅不去证明它，但这个性质非常重要，后面会有很多向量的理论都建立在它的基础之上。例如，对一个向量 (x⃗ +y⃗ ) ,利用这个性质，结合公式 1，我们可以得到

∥∥x⃗ +y⃗ ∥∥2  =(x⃗ +y⃗ )⋅(x⃗ +y⃗ )=∥∥x⃗ ∥∥2+2x⃗ y⃗ +∥∥y⃗ ∥∥2≤∥∥x⃗ ∥∥2+2∥∥x⃗ ∥∥∥∥y⃗ ∥∥+∥∥y⃗ ∥∥2

所以：

∥∥x⃗ +y⃗ ∥∥2≤(∥∥x⃗ ∥∥+∥∥y⃗ ∥∥)2

两边开平方得到：

∥∥x⃗ +y⃗ ∥∥≤∥∥x⃗ ∥∥+∥∥y⃗ ∥∥

这就得到了三角不等式。

从几何的角度来说，向量的点积与向量间夹角 θ 的余弦有关：

a⃗ ⋅b⃗ =∥∥a⃗ ∥∥∥∥b⃗ ∥∥cosθ

，这意味着向量的点积其实反映了向量 a⃗  在向量 b⃗  上的  投影 ，即两个向量在同个方向上的相同程度。当两向量正交时， cosθ 的值为0，点积的值为0，投影最小。当两向量平行时， cosθ 的值为1，点积值最大，投影也最大。

观察上图， L 是 v⃗  向量两端延伸出来的直线，即 L=cv⃗ |c∈R 。记向量 x⃗  在 L 上的投影为 ProjL(x⃗ ) 。根据点积的性质，可得：

(x⃗ −cv⃗ ProjL(x⃗ ))⋅v⃗  x⃗ ⋅v⃗ −cv⃗ ⋅v⃗  c⋅v⃗ ⋅v⃗  c=0=0=x⃗ ⋅v⃗ =x⃗ ⋅v⃗ v⃗ ⋅v⃗ 

有了 c ， 我们就可以求出投影 ProjL(x⃗ ) 为：

ProjL(x⃗ )=cv⃗ =(x⃗ ⋅v⃗ v⃗ ⋅v⃗ )v⃗ 

例如，向量 a⃗ =[12] ，向量 b⃗ =[11] ，那么 a⃗  在 b⃗  方向 L 上的投影为：

ProjL(a⃗ )=cb⃗ =(a⃗ ⋅b⃗ b⃗ ⋅b⃗ )b⃗ =32b⃗ 

Python 示例：

defget_projection(a, b):
    return a.dot(b)*1.0*b/b.dot(b)

a = np.array([1, 2])
b = np.array([2, 2])
print get_projection(a, b)  # [1.5 1.5]

向量外积

向量的（又叫叉乘、向量积、叉积）只在 R2 和 R3 中定义：

R2 的向量外积：

[a1a2]×[b1b2]=[a1⋅b2−a2⋅b1]

例如：

[12]×[34]=[1⋅4−3⋅2]=[−2]

R3 的向量外积：

⎡⎣⎢a1a2a3⎤⎦⎥×⎡⎣⎢b1b2b3⎤⎦⎥=⎡⎣⎢a2⋅b3−a3⋅b2a3⋅b1−a1⋅b3a1⋅b2−a2⋅b1⎤⎦⎥

例如：

⎡⎣⎢352⎤⎦⎥×⎡⎣⎢147⎤⎦⎥=⎡⎣⎢5⋅7−2⋅42⋅1−3⋅73⋅4−5⋅1⎤⎦⎥=⎡⎣⎢27−197⎤⎦⎥

可见向量间外积的结果会得到一个新的向量。

Python 示例：

import numpy as np
a = np.array([3, 5, 2])
b = np.array([1, 4, 7])
print np.cross(a, b)  # [27, -19, 7]

外积的一个重要作用是可以得到一个和 a⃗  、 b⃗  两个原向量正交的新向量 c⃗  ，且可以通过右手法则来确定新向量的方向（一个简单的确定满足“右手定则”的结果向量的方向的方法是这样的：若坐标系是满足右手定则的，当右手的四指从 a⃗  以不超过180度的转角转向 b⃗  时，竖起的大拇指指向是 c⃗  的方向）。

从几何的角度来说，向量的外积与向量间夹角 θ 的正弦有关：

∥∥a⃗ ×b⃗ ∥∥=∥∥a⃗ ∥∥∥∥b⃗ ∥∥sinθ

，这意味着向量的外积反映了向量 a⃗  与向量 b⃗  的正交程度。当两向量平行时， sinθ 的值为0，外积的值为0，正交程度最小。当两向量正交时， sinθ 的值为1，外积值最大，正交程度最大。

矩阵向量积

当矩阵 A 的列数与向量 x⃗  的分量数相同时，矩阵和向量的积有定义：

Am×nx⃗ =⎡⎣⎢⎢⎢a11a21…am1a12a22am2………a1na2namn⎤⎦⎥⎥⎥⎡⎣⎢⎢⎢x1x2…xn⎤⎦⎥⎥⎥=⎡⎣⎢⎢⎢a11x1+a12x2+…+a1nxna21x1+a22x2+…+a2nxn…am1x1+am2x2+…+amnxn⎤⎦⎥⎥⎥

例如矩阵 A=[413215] 乘以向量 x⃗ =⎡⎣⎢527⎤⎦⎥ 的结果为：

[4⋅5+3⋅2+1⋅71⋅5+2⋅2+5⋅7]=[3344]

Python 示例：

a = np.matrix('4 3 1;1 2 5')
x = np.array([[5], [2], [7]])
print a*x  # [[33] [44]]

矩阵的向量积可以当成是矩阵的所有列向量的线性组合：

Am×nx⃗ =⎡⎣⎢⎢⎢⎢⎢⎢⎢⎡⎣⎢⎢⎢a11a21…am1⎤⎦⎥⎥⎥V1→⎡⎣⎢⎢⎢a12a22…am2⎤⎦⎥⎥⎥V2→…⎡⎣⎢⎢⎢a1na2n…amn⎤⎦⎥⎥⎥Vn→⎤⎦⎥⎥⎥⎥⎥⎥⎥⎡⎣⎢⎢⎢x1x2…xn⎤⎦⎥⎥⎥=x1V1→+x2V2→+…+xnVn→

而向量 x⃗  的每一个分量可以看成是 A 的每个列向量的加权。

向量的转置

向量 V⃗ =⎡⎣⎢⎢⎢v1v2…vn⎤⎦⎥⎥⎥n×1 的转置定义为 V⃗ T=[v1v2…vn]1×n

例如向量 A⃗ =[24] 的转置就是 A⃗ T=[24] 。

Python 示例：

>>> a = np.array([[2, 4]])
>>> a.T
array([[2],
       [4]])

注意上面声明 a 时用了两对 [] ，以生成一个二维向量。一维的向量转置结果是不会变化的：

>>> b = np.array([2, 4])
>>> b.T
array([2, 4])

向量的转置有两个性质：一个向量 v⃗  点乘另一个向量 w⃗  ，其结果和向量 v⃗  转置后和向量 w⃗  做矩阵乘法相同。即 v⃗ ⋅w⃗ =v⃗ Tw⃗  。

线性无关

张成空间

一组向量的张成空间说白了就是指这些向量随便线性组合后能够表示多少个向量。记为 span(a⃗ ,b⃗ ) 。

例如，对于 R2 空间中两个不平行的非0向量 a⃗ =[21] 和向量 b⃗ =[03] ，不难发现这两个向量能够表示二维空间中任一其他向量，即 span(a⃗ ,b⃗ )=R2 。证明如下：

对于 R2 中任一向量 [xy] ，假设可以由 a⃗  和 b⃗  线性组合而成，那么有：

c1[21]+c2[03]=[xy]

即：

{c1⋅2 c1⋅1+c2⋅0+c2⋅3=x=y

求解该方程得：

⎧⎩⎨⎪⎪c1c2=x2=y3−x6

由于 x 、 y 的值已确定，所以 c1 、 c2 的值也必然唯一。

线性相关和线性无关

当一个向量集合里的每个向量都对张成的空间有贡献时，称这个向量集合线性无关。反之称为线性无关。能够表示一个空间的最少向量组合称为空间的基。

听起来有点难理解，其实就是非常简单的道理：假如一个向量集合中存在某个向量能由集合里的其他向量线性组合而成，那这个集合对于张成空间而言就存在多余的向量。此时就是线性相关；反之，假如集合里每一个元素都没法由其他元素组合而成，那么这个集合每个元素都对张成空间有贡献，这个集合就是线性无关的。

例如，对于上述的例子，如果再增加一个向量 c⃗ =[52] ，由于 c⃗  可以由 a⃗  和 b⃗  线性组合而成，由 a 、 b 和 c 共同张成的空间并没有变化，仍然是 R2 ，因此称集合 {a⃗ ,b⃗ ,c⃗ } 线性相关。

判断是否线性相关

一个向量集合 s=v1,v2,…,vn 线性相关的充分必要条件是存在一部分非0系数使得 c1v1+c2v2+…+cnvn=0=⎡⎣⎢⎢⎢00…0⎤⎦⎥⎥⎥ 。

例如有向量 [21] 和 [32] ，则可以先写出如下的等式：

c1[21]+c2[32]=[00]

容易求解得 [c1c2]=[00] ，说明两个向量线性无关。也说明这两个向量可以张成 R2 。

类似地，对于三个 R3 中的向量 ⎡⎣⎢200⎤⎦⎥ 、 ⎡⎣⎢010⎤⎦⎥ 和 ⎡⎣⎢007⎤⎦⎥ ，不难判断这三个向量是线性无关的，他们共同张成了 R3 空间。

而对于向量集合 {[21],[32],[12]} ，不难算出存在非 0 的系数 ⎡⎣⎢c1c2c3⎤⎦⎥=⎡⎣⎢−43−1⎤⎦⎥ 使得 c1[21]+c2[32]+c3[12]=[00] 。因此集合 {[21],[32],[12]} 线性相关。

下篇文章将进阶讨论线性子空间和特征向量。