代码之家 › 专栏 › 技术社区 › b-fg

二维元素的Numpy矩阵乘法

matrix-multiplication numpy python

b-fg · 技术社区 · 6 年前

我有一个 a 努比 ndarray

a =  ([[ uu, uv, uw],
       [ uv, vv, vw],
       [ uw, vw, ww]])

每个组件本身就是一个大小为的二维数组 (N,M) ,所以 一 矩阵有一个 (3,3,N,M)

我怎样才能执行 a*a 以蟒蛇的方式? a@a 抛出以下错误(对于N=1218和M=540):

(尺寸3)!=1218(尺寸2)

我希望能够执行这个操作,就像其中只有标量值不会抛出与其形状相关的错误,因为它是一个简单的3x3矩阵乘法。

谢谢。

1 回复 | 直到 6 年前

Divakar 6 年前

假设您希望沿最后两个轴对每个元素执行矩阵乘法,我们可以使用 np.einsum -

np.einsum('ijkl,jmkl->imkl',a,a)

样品运行验证-

In [43]: np.random.seed(0)

In [44]: a = np.random.rand(3,3,4,5)

In [45]: a[:,:,0,0].dot(a[:,:,0,0])
Out[45]: 
array([[0.71750146, 1.17057872, 1.11135764],
       [0.62938365, 0.86437796, 0.74541383],
       [1.04636618, 1.62011127, 1.35483565]])

In [46]: np.einsum('ijkl,jmkl->imkl',a,a)[:,:,0,0]
Out[46]: 
array([[0.71750146, 1.17057872, 1.11135764],
       [0.62938365, 0.86437796, 0.74541383],
       [1.04636618, 1.62011127, 1.35483565]])

推荐文章

Google User · Django管理员在`list_display中未显示`creation_date`字段`

5 月前

user29747013 · 如何创建一个新的数据框架,其中包含原始数据框架中列的聚合列?

5 月前

ÎÎÎ½Î· ÎÎ®Î¹Î½Î¿Ï · Python lxml.html语法错误:使用lxml find时XPATH的谓词无效

5 月前

user29715306 · from_users=和chats=电视节目中的差异

5 月前

Redshoe · 当执行numpy.genfromtxt()时,python是否会读取文件的所有行?

5 月前

RASEL MAHMUD · 为什么以及如何在is_even()函数内的IF条件中递归X变量在满足0后递增?[副本]

5 月前

prayner · 更新嵌套字典包含列表中的项

5 月前

Bringo Jr · 我可以在O(n)中解决这个问题吗?

5 月前

Dave · 如何在for循环中修改列表值

5 月前

Shukurullox Komiljonov · 从记录中获得相互和解。使用SQL

5 月前