NumPy 教程（第 9 章）：数组广播

广播(Broadcast)是 numpy 对不同形状(shape)的数组进行数值计算的方式， 对数组的算术运算通常在相应的元素上进行

如果两个数组 a 和 b 形状相同，即满足 a.shape == b.shape，那么 a*b 的结果就是 a 与 b 数组对应位相乘。这要求维数相同，且各维度的长度相同

In [1]: import numpy as np

In [2]: a = np.array([1,2,3,4])

In [3]: b = np.array([10,10,10,10])

In [4]: a
Out[4]: array([1, 2, 3, 4])

In [5]: b
Out[5]: array([10, 10, 10, 10])

In [6]: num = a * b

In [7]: num
Out[7]: array([10, 20, 30, 40])

当运算中的 2 个数组的形状不同时，numpy 将自动触发广播机制

In [1]: import numpy as np

In [2]: a = np.array([[ 0, 0, 0], [10,10,10], [20,20,20], [30,30,30]])

In [3]: b = np.array([1,2,3])

In [4]: a
Out[4]:
array([[ 0,  0,  0],
       [10, 10, 10],
       [20, 20, 20],
       [30, 30, 30]])

In [5]: b
Out[5]: array([1, 2, 3])

In [6]: sum = a + b

In [7]: sum
Out[7]:
array([[ 1,  2,  3],
       [11, 12, 13],
       [21, 22, 23],
       [31, 32, 33]])

In [8]: sum = a * b

In [9]: sum
Out[9]:
array([[ 0,  0,  0],
       [10, 20, 30],
       [20, 40, 60],
       [30, 60, 90]])

4x3 的二维数组与长为 3 的一维数组相加，等效于把数组 b 在二维上重复 4 次再运算

In [10]: a
Out[10]:
array([[ 0,  0,  0],
       [10, 10, 10],
       [20, 20, 20],
       [30, 30, 30]])

In [11]: b
Out[11]: array([1, 2, 3])

In [12]: b = np.tile(b, (4, 1))

In [13]: b
Out[13]:
array([[1, 2, 3],
       [1, 2, 3],
       [1, 2, 3],
       [1, 2, 3]])

In [14]: sum = a + b

In [15]: sum
Out[15]:
array([[ 1,  2,  3],
       [11, 12, 13],
       [21, 22, 23],
       [31, 32, 33]])

广播的规则:

• 让所有输入数组都向其中形状最长的数组看齐，形状中不足的部分都通过在前面加 1 补齐

• 输出数组的形状是输入数组形状的各个维度上的最大值

• 如果输入数组的某个维度和输出数组的对应维度的长度相同或者其长度为 1 时，这个数组能够用来计算，否则出错

• 当输入数组的某个维度的长度为 1 时，沿着此维度运算时都用此维度上的第一组值

简单理解：

对两个数组，分别比较他们的每一个维度（若其中一个数组没有当前维度则忽略），满足：

• 数组拥有相同形状

• 当前维度的值相等

• 当前维度的值有一个是 1

若条件不满足，抛出 ValueError: frames are not aligned 异常