笔记 | numpy-09.广播(Broadcast)

#%%
# 广播(Broadcast)是 numpy 对不同形状(shape)的数组进行数值计算的方式， 对数组的算术运算通常在相应的元素上进行
# 如果两个数组 a 和 b 形状相同，即满足 a.shape == b.shape，那么 a*b 的结果就是 a 与 b 数组对应位相乘。这要求维数相同，且各维度的长度相同
#%%
import numpy as np

a = np.array([1, 2, 3, 4]) # 4r1c
a
#%%
b = np.array([10, 20, 30, 40]) # 4r1c
b
#%%
c = a * b
c
#%%
# 当运算中的 2 个数组的形状不同时，numpy 将自动触发广播机制
import numpy as np

a = np.array([
    [0, 0, 0],
    [10, 10, 10],
    [20, 20, 20],
    [30, 30, 30]
])
a # 4r3c
#%%
b = np.array([1, 2, 3])
b # 1r3c
#%%
"""
a行索引0 0+1=1 0+2=2 0+3=3
a行索引1 10+1=11 10+2=12 10+3=13
a行索引2 20+1=21 20+2=22 20+3=23
a行索引3 30+1=31 30+2=32 30+3=33
"""
a + b # 将b依次带入每行的a | a的列对应b的行 | 计算时会与找边个数相同的进行依次计算
#%%
bb = np.tile(b, (4, 1)) # 重复 b的各维度 | 1 垂直方向重复 2 水平方向重复 | 4 重复4次
bb
#%%
a + bb # 等同于 a + b
#%%
"""
广播的规则:

让所有输入数组都向其中形状最长的数组看齐，形状中不足的部分都通过在前面加 1 补齐。
输出数组的形状是输入数组形状的各个维度上的最大值。
如果输入数组的某个维度和输出数组的对应维度的长度相同或者其长度为 1 时，这个数组能够用来计算，否则出错。
当输入数组的某个维度的长度为 1 时，沿着此维度运算时都用此维度上的第一组值

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

数组拥有相同形状。
当前维度的值相等。
当前维度的值有一个是 1。
若条件不满足，抛出 "ValueError: frames are not aligned" 异常
"""