Python如何计算NumPy数组的协方差矩阵_调用cov函数进行特征分析

星敏君_4200

星敏君_4200

2026-03-30

937人浏览

原创

np.cov 默认按行计算协方差,即每行视为一个变量、每列视为一个样本;若数据为标准格式 (n_samples, n_features),需显式设置 rowvar=false 才能得到正确的 n_features × n_features 协方差矩阵。

python如何计算numpy数组的协方差矩阵_调用cov函数进行特征分析

np.cov 默认按行还是按列计算协方差?

默认按**行**把每行当作一个变量(即观测维度),列是样本点——这和大多数机器学习场景的直觉相反,容易导致协方差矩阵形状错乱、特征相关性算反。

  • 如果你的数据是 shape=(n_samples, n_features)(标准格式),直接传给 np.cov 会把每个样本当变量,结果矩阵是 n_samples × n_samples,完全没用
  • 正确做法是加参数 rowvar=False,告诉它“变量在列上”,此时才输出 n_features × n_features 协方差矩阵
  • 不设 rowvar 时,np.cov(X) 等价于 np.cov(X, rowvar=True),别依赖默认值

为什么 np.cov(X.T) 和 np.cov(X, rowvar=False) 结果一样但要慎用?

两者数学等价,但 X.T 在高维或内存受限时可能触发复制开销,尤其对大数组;而 rowvar=False 是原地逻辑切换,更轻量。

  • np.cov(X.T) 强制生成转置副本,若 X 是 100 万 × 100 的 float64,转置后多占约 800MB 内存
  • np.cov(X, rowvar=False) 不额外分配,适合生产环境或 notebook 中反复调试
  • 注意:如果 X 是 Fortran-order(列优先)数组,X.T 是 view,但 rowvar=False 依然更明确、可读性更强

NaN 值导致 np.cov 返回全 nan 怎么办?

np.cov 遇到任何 nan 就整行/整列失效,默认返回全 nan 矩阵,不是跳过或报错。

Python Packaging
Python Packaging

深度Python打包工作流——pyproject元数据、依赖与可选额外项、构建后端、wheel、版本控制、发布及CI发布规范……

下载
  • 用 nan_policy='omit' 参数可启用剔除含 nan 的样本对(pairwise deletion),这是最常用解法
  • 不要用 nan_policy='propagate'(默认)或 'raise',前者静默失败,后者直接中断
  • 示例:np.cov(X, rowvar=False, nan_policy='omit') —— 它对每一对特征独立丢弃该对中任一为 nan 的样本行
  • 注意:pairwise 删除可能导致协方差矩阵非正定,后续做 PCA 或 Cholesky 分解前需检查

协方差矩阵数值不准?检查数据是否已中心化

np.cov 内部会自动减均值,但前提是它能正确算出均值;若数据含 inf、极端离群值或 dtype 不匹配(如 int64 与 float32 混用),均值计算可能溢出或截断,导致协方差失真。

  • 先用 np.isfinite(X).all() 确保无 inf/nan
  • 显式中心化更可控:X_centered = X - X.mean(axis=0),再用 np.dot(X_centered.T, X_centered) / (X.shape[0] - 1)
  • dtype 推荐统一为 float64,避免 int32 数组在求和时悄悄溢出
  • 小样本(n_samples )下,样本协方差估计方差大,别急着下结论

协方差矩阵本身不标准化,数值大小严重依赖原始特征量纲;做特征分析前,通常得先标准化或直接用相关系数矩阵(np.corrcoef)。这点常被跳过,但直接影响后续解释。

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

相关专题

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

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

2023.07.20

1631

4

python能做什么
python能做什么

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

2023.07.25

4024

7

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

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

2023.07.31

1629

3

python教程
python教程

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

2023.08.03

23157

23

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

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

2023.08.04

2847

5

python eval
python eval

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

2023.08.04

2887

5

scratch和python区别
scratch和python区别

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

2023.08.11

1123

5

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

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

2023.08.10

596

4

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

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

2023.08.11

2223

5

热门下载

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

精品课程

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