有效地减去不同形状的numpy数组

Question

使用 numpy 的优秀广播规则，您可以使用

从形状 (5,3) 数组 X 中减去形状 (3,) 数组 v

X - v

结果是一个形状 (5,3) 数组，其中每一行 i 是差 X[i] - v。

有没有办法从 X 中减去形状 (n,3) 数组 w，以便在不显式使用循环的情况下从整个数组 X 中减去 w 的每一行？

Answer 1

您需要将X 的维度扩展为None/np.newaxis 以形成一个3D 数组，然后通过w 进行减法运算。这将在3D 操作中引入broadcasting，并产生形状为(5,n,3) 的输出。实现看起来像这样 -

X[:,None] - w  # or X[:,np.newaxis] - w

相反，如果所需的排序是(n,5,3)，那么您需要改为扩展w 的维度，就像这样 -

X - w[:,None] # or X - w[:,np.newaxis] 

示例运行 -

In [39]: X
Out[39]: 
array([[5, 5, 4],
       [8, 1, 8],
       [0, 1, 5],
       [0, 3, 1],
       [6, 2, 5]])

In [40]: w
Out[40]: 
array([[8, 5, 1],
       [7, 8, 6]])

In [41]: (X[:,None] - w).shape
Out[41]: (5, 2, 3)

In [42]: (X - w[:,None]).shape
Out[42]: (2, 5, 3)