怎么在Python中使用NumPy实现矩阵高效相乘

秋浩小哥_3700

秋浩小哥_3700

2026-10-08

539人浏览

原创

优先用 @ 运算符——语义专一、可读性强、报错及时;np.dot() 语义复杂,一维时为内积,高维时易因自动轴对齐导致意外结果。

怎么在python中使用numpy实现矩阵高效相乘

用 np.dot() 还是 @ 运算符?选哪个更稳妥

两者在绝大多数场景下行为一致,但语义和可读性有差别。@ 是 Python 3.5+ 引入的矩阵乘法专用运算符,专为二维数组(或最后两维为矩阵)设计;np.dot() 历史更久,但语义稍杂:对一维数组做内积,对二维数组才做矩阵乘,高维时按规则收缩最后/倒数第二维,容易误用。

实操建议:

  • 明确做矩阵乘时,优先用 @ —— 更直观、不易歧义,比如 A @ B
  • 兼容老版本 Python(np.dot(A, B)
  • 避免用 np.multiply() 或 *,那是逐元素乘,不是矩阵乘

输入维度不匹配时,ValueError: matmul: Input operand X has a mismatch in its core dimension 怎么快速定位

这是 @ 和 np.matmul() 最常见的报错,本质是「左矩阵列数 ≠ 右矩阵行数」。但实际中常因数组维度隐式变化而难以察觉,比如把一维向量当行向量用却没 reshape。

排查要点:

  • 用 .shape 显式检查:确保 A.shape[-1] == B.shape[-2](对高维也适用)
  • 一维数组参与矩阵乘时,默认被当作行向量((n,) 视为 (1, n))或列向量((n,) 视为 (n, 1)),取决于位置:v @ A 要求 v.shape[0] == A.shape[0],而 A @ v 要求 A.shape[-1] == v.shape[0]
  • 不确定时,显式 reshape:用 v.reshape(-1, 1) 强制列向量,v.reshape(1, -1) 强制行向量

np.matmul() 和 np.dot() 在三维及以上数组中的行为差异

这是最容易踩坑的地方:二者对高维数组的广播逻辑不同。np.matmul() 严格按「批处理矩阵乘」设计,只对最后两维做矩阵乘,前面维度广播对齐;np.dot() 则会对倒数第一维和倒数第二维做求和,导致结果维度可能意外变化。

示例对比:

python-pro
python-pro

高级 Python 特性、异步编程、性能调优、静态类型、内存管理、Python 内部机制及生态库方面的专家。

下载
import numpy as np
A = np.ones((2, 3, 4))  # batch=2, 3x4 矩阵
B = np.ones((2, 4, 5))  # batch=2, 4x5 矩阵
<p>print(np.matmul(A, B).shape)  # (2, 3, 5) ✅ 批处理乘
print(np.dot(A, B).shape)    # (2, 3, 2, 5) ❌ 对 A 的 -1 维和 B 的 -2 维求和,引入新轴</p>

实操建议:

  • 只要涉及三维及以上(如 batched image features、Transformer attention weights),一律用 np.matmul() 或 @
  • 除非你明确需要 np.dot() 的旧式高维收缩行为(极少)

为什么有时候用 np.einsum('ij,jk->ik', A, B) 反而更快?

np.einsum() 在某些硬件或特定尺寸组合下,能绕过 NumPy 默认 BLAS 后端的调度开销,或触发更优的底层实现(如 Intel MKL 的特殊优化路径)。但它不是“总是更快”,而是有条件优势。

适用场景:

  • 小矩阵(如 16x16 或 32x32),BLAS 启动成本占比高,einsum 可省掉这部分
  • 需要融合多个操作(如乘完立刻加偏置、再激活),用 einsum 可减少中间数组内存分配
  • 非标准收缩模式(如 'ijk,ilk->ijl'),此时 matmul 无法直接表达

注意:einsum 默认不使用多线程,若想压满 CPU,需配合 optimize=True 并确认后端支持。

实际项目里,95% 的矩阵乘用 @ 就够了;但一旦出现维度报错或性能瓶颈,别只盯着算法,先看形状是否真对齐、再确认用的是不是 matmul 而非 dot——这两个细节,比调参更容易卡住进度。

Python免费学习笔记(深入):立即使用
在学习笔记中,你将探索 Python 的核心概念和高级技巧!

相关文章

PHP速学视频免费教程(入门到精通)
PHP速学视频免费教程(入门到精通)

PHP怎么学习?PHP怎么入门?PHP在哪学?PHP怎么学才快?不用担心,这里为大家提供了PHP速学教程(入门到精通),有需要的小伙伴保存下载就能学习啦!

下载

相关标签:

python

本站声明:本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系admin@php.cn

相关专题

更多
python打包成可执行文件
python打包成可执行文件

本专题为大家带来python打包成可执行文件相关的文章,大家可以免费的下载体验。

2023.07.20

1671

4

python能做什么
python能做什么

python能做的有:可用于开发基于控制台的应用程序、多媒体部分开发、用于开发基于Web的应用程序、使用python处理数据、系统编程等等。本专题为大家提供python相关的各种文章、以及下载和课程。

2023.07.25

4184

7

format在python中的用法
format在python中的用法

Python中的format是一种字符串格式化方法,用于将变量或值插入到字符串中的占位符位置。通过format方法,我们可以动态地构建字符串,使其包含不同值。php中文网给大家带来了相关的教程以及文章,欢迎大家前来阅读学习。

2023.07.31

1669

3

python教程
python教程

Python已成为一门网红语言,即使是在非编程开发者当中,也掀起了一股学习的热潮。本专题为大家带来python教程的相关文章,大家可以免费体验学习。

2023.08.03

24237

23

python环境变量的配置
python环境变量的配置

Python是一种流行的编程语言,被广泛用于软件开发、数据分析和科学计算等领域。在安装Python之后,我们需要配置环境变量,以便在任何位置都能够访问Python的可执行文件。php中文网给大家带来了相关的教程以及文章,欢迎大家前来学习阅读。

2023.08.04

2967

5

python eval
python eval

eval函数是Python中一个非常强大的函数,它可以将字符串作为Python代码进行执行,实现动态编程的效果。然而,由于其潜在的安全风险和性能问题,需要谨慎使用。php中文网给大家带来了相关的教程以及文章,欢迎大家前来学习阅读。

2023.08.04

2987

5

scratch和python区别
scratch和python区别

scratch和python的区别:1、scratch是一种专为初学者设计的图形化编程语言,python是一种文本编程语言;2、scratch使用的是基于积木的编程语法,python采用更加传统的文本编程语法等等。本专题为大家提供scratch和python相关的文章、下载、课程内容,供大家免费下载体验。

2023.08.11

1163

5

python合并两个列表
python合并两个列表

Python是一种强大的编程语言,具有许多方便的功能和工具。在Python中,有多种方法可以合并两个列表。php中文网给大家带来了相关的教程以及文章,欢迎大家前来学习阅读。

2023.08.10

596

4

python是前端还是后端
python是前端还是后端

Python属于前端也属于后端,其灵活性和丰富的生态系统使得开发人员能够在不同的领域中灵活运用。本专题为大家提供python相关的文章、下载、课程内容,供大家免费下载体验。

2023.08.11

2323

5

热门下载

更多
网站特效
/
网站源码
/
网站素材
/
前端模板

精品课程

更多
相关推荐
/
热门推荐
/
最新课程