Python怎么处理稀疏矩阵_SciPy库稀疏矩阵四种存储格式对比

酷瑶吖_4234

酷瑶吖_4234

2026-03-21

980人浏览

原创

选csc_matrix还是csr_matrix取决于操作方向:列切片用csc_matrix,行切片用csr_matrix;coo_matrix仅适用于构造阶段,需转为csr/csc才能运算;lil_matrix仅适合动态构建,完成后必须转换;拼接矩阵前须统一格式并校验形状。

python怎么处理稀疏矩阵_scipy库稀疏矩阵四种存储格式对比

稀疏矩阵该选 csc_matrix 还是 csr_matrix?

取决于你主要做「列切片」还是「行切片」——csc_matrix(Compressed Sparse Column)对列操作快,csr_matrix(Compressed Sparse Row)对行操作快。底层是用三个数组存非零值、行索引、列指针(或反之),结构差异直接决定访问效率。

常见错误现象:mat[:, 5] 在 csr_matrix 上会触发全矩阵转置再切片,慢得离谱;反过来,mat[3, :] 在 csc_matrix 上也一样。

  • 做特征提取、按列标准化(比如每列减均值)→ 优先用 csc_matrix
  • 做样本采样、按行计算 L2 范数、批量预测 → 优先用 csr_matrix
  • 从 scipy.sparse.load_npz() 加载后默认是 csr,别直接拿去列操作
  • 转换开销不可忽略:mat.tocsc() 或 mat.tocsr() 是 O(nnz) 时间,不是常数

coo_matrix 为什么不能直接做乘法?

coo_matrix 是最原始的三元组格式(行索引、列索引、值),不支持算术运算和切片,因为没建索引结构——它连「某行有哪些非零元」都得遍历找。

使用场景:只适合构造阶段。比如从文件逐行读取 (i, j, val),append 到三个列表里,最后一次性转成 coo_matrix,再立刻转成 csr 或 csc。

  • 误用 coo_matrix 做 @ 矩阵乘 → 报错 TypeError: unsupported operand type(s) for @: 'coo_matrix' and 'coo_matrix'
  • coo_matrix 允许重复 (i,j),但转成 csr/csc 时会自动求和,这是隐式行为,容易漏掉数据异常
  • 构造完必须转:用 .tocsr() 或 .tocsc(),别留着 coo 做后续计算

lil_matrix 适合动态插入,但千万别用它存最终结果

lil_matrix(List of Lists)内部用 Python 列表存每行的列索引和值,支持 mat[i, j] = x 这种随机写入,但内存占用高、访问慢、不支持广播运算。

python-parallelization
python-parallelization

将顺序 Python 代码转换为并行/并发实现,用于在需要并行化 Python 代码、提升性能时使用。

下载

性能影响明显:10 万行 × 10 万列的稀疏矩阵,用 lil 存可能吃掉 2–3 倍内存,且 mat.sum(axis=1) 比 csr 慢一个数量级。

  • 只在「边生成边填」场景用:比如迭代算法中逐行构造矩阵
  • 填完立刻转:mat.tocsr() 是标准收尾动作,不是可选项
  • 别用 lil_matrix 做 np.dot 或传给 sklearn 模型——很多模型内部会检查格式,遇到 lil 直接报 ValueError: Expected CSR or CSC matrix

用 scipy.sparse.bmat 拼接矩阵时维度对不上怎么办?

bmat 要求子块在逻辑上能排成矩形,比如 [[A, B], [C, D]] 中,A.shape[0] 必须等于 C.shape[0],A.shape[1] 等于 B.shape[1],否则报 ValueError: blocks must be 2-D 或更模糊的 inconsistent shapes。

容易踩的坑是混用不同格式:比如 A 是 csr,B 是 coo,bmat 不会自动统一,可能在拼接中途因格式不兼容崩掉。

  • 拼之前先统一格式:A.tocsr(), B.tocsr(), C.tocsr(), D.tocsr()
  • 手动校验形状:assert A.shape[0] == C.shape[0] 和 assert A.shape[1] == B.shape[1],比看报错信息快得多
  • 如果只是水平拼([A, B]),用 scipy.sparse.hstack([A, B]) 更稳;垂直拼用 vstack,它们对形状要求更明确,报错也更直白

稀疏格式转换不是“无感”的,每次 .tocsr() 都在重排数据;而不同格式对同一操作的时间差可以到百倍。真正卡住你的往往不是算法,而是矩阵刚加载完就用错了格式。

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

相关文章

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

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

下载

相关标签:

python

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

相关专题

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

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

2023.07.20

1731

4

python能做什么
python能做什么

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

2023.07.25

4384

7

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

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

2023.07.31

1709

3

python教程
python教程

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

2023.08.03

25697

23

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

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

2023.08.04

3127

5

python eval
python eval

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

2023.08.04

3147

5

scratch和python区别
scratch和python区别

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

2023.08.11

1183

5

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

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

2023.08.10

596

4

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

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

2023.08.11

2443

5

热门下载

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

精品课程

更多
相关推荐
/
热门推荐
/
最新课程
Uvicorn官方安装与入门手册
Uvicorn官方安装与入门手册

共0课时 | 0人学习

Python uv官方项目管理指南
Python uv官方项目管理指南

共0课时 | 0人学习

Python uv官方安装与入门文档
Python uv官方安装与入门文档

共0课时 | 0人学习