如何在Python中利用NumPy进行主成分分析(PCA)降维?

冬丽大大_2704

冬丽大大_2704

2026-08-30

350人浏览

原创

必须先标准化再调用np.cov(),否则量纲大的特征会主导协方差矩阵,导致主成分方向偏移;标准化需减均值再除标准差(ddof=1),并确保中心化后均值为0、标准差为1。

如何在python中利用numpy进行主成分分析(pca)降维?

为什么必须先标准化数据再调用 np.cov()

不标准化直接算协方差矩阵,结果会被量纲大的特征主导。比如一个特征单位是「毫米」,另一个是「万元」,前者数值天然大几个数量级,np.cov() 输出的协方差矩阵就严重失真,主成分方向完全偏移。

标准化要做的只有两步:减均值(X - np.mean(X, axis=0)),再除标准差(/ np.std(X, axis=0, ddof=1))。注意 ddof=1 才符合样本协方差定义,否则 np.cov() 默认用 ddof=0,和 sklearn 的行为不一致。

  • 别用 sklearn.preprocessing.StandardScaler 后再自己算协方差——它返回的是 ndarray,但内部用了 ddof=0,和纯 NumPy 手动实现有细微偏差
  • 如果原始数据含缺失值,np.mean() 和 np.std() 默认会返回 nan,得提前用 np.nanmean() / np.nanstd() 并设 axis=0
  • 标准化后务必检查:np.allclose(np.mean(X_std, axis=0), 0) 和 np.allclose(np.std(X_std, axis=0, ddof=1), 1) 应该都为 True

手算协方差矩阵后,该用 np.linalg.eig() 还是 np.linalg.svd()

np.linalg.eig() 直接作用于协方差矩阵,返回特征值和特征向量,逻辑清晰;np.linalg.svd() 对中心化后的数据矩阵做奇异值分解,更数值稳定,尤其当特征数远大于样本数时(比如 1000 维、50 个样本),np.cov() 会得到病态矩阵,eig() 可能崩溃或返回复数特征向量。

实际选哪个,看数据形状:

  • 若 X_centered.shape[1] (特征数 ≤ 样本数),用 <code>np.linalg.eig(np.cov(X_centered.T)) 更直观
  • 若特征数远大于样本数,改用 U, s, Vt = np.linalg.svd(X_centered, full_matrices=False),主成分就是 Vt.T 的列,前 k 列即降维投影矩阵
  • 两种方法得到的主成分方向可能符号相反(特征向量可反向),不影响投影结果,但别拿它们直接比对 pca.components_

pca.explained_variance_ratio_ 和手动计算的累计方差比对不上怎么办

这是因为 sklearn 的 explained_variance_ratio_ 默认基于 SVD 的奇异值平方(即特征值)归一化,而手动用 eig() 算出的特征值若没排序、没归一化,或者用了不同自由度的协方差,就会不一致。

361SCAN security check your skill python scripts
361SCAN security check your skill python scripts

扫描并分析已安装技能,用于查看指定技能目录的名称、描述和详情,或列出所有已安装的技能。

下载

要对齐结果,手动计算时必须:

  • 对特征值数组 eigvals 用 np.sort(eigvals)[::-1] 降序排列
  • 用 eigvals / eigvals.sum() 计算比例,不是除以 len(eigvals) 或其他常数
  • 确认协方差矩阵是 np.cov(X_centered.T, ddof=1),不是 .T 写反或漏掉 ddof=1
  • 如果用 SVD,奇异值 s 要平方再归一:s**2 / np.sum(s**2)

降维投影时,X @ components.T 和 X_centered @ components.T 哪个对

必须用中心化后的数据 X_centered。PCA 的投影本质是把数据点往主成分方向做正交投影,这个操作只对去均值后的坐标有意义。原始 X 没减均值,投影结果整体偏移,后续任何距离、聚类、可视化都会出错。

典型错误写法:X_pca = X @ components.T —— 这相当于在原点未移动的情况下硬投,第一主成分坐标均值不为 0,pca.explained_variance_ 也会虚高。

  • 正确流程固定三步:标准化 → 得 X_centered → 投影:X_pca = X_centered @ components.T
  • 如果后续要还原(inverse_transform),也必须用同一套均值和标准差,不能只存 components
  • 用 sklearn 时,fit_transform() 内部自动完成中心化,但手写时这一步绝不能跳过

手算 PCA 最容易被忽略的,不是公式推导,而是标准化与中心化的数值细节——ddof、nan 处理、矩阵转置方向、投影输入是否已中心化,每一步错一点,结果就漂移。跑通一个 Iris 数据集不难,但在真实业务数据上复现 sklearn 结果,往往卡在这几步。

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

相关文章

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

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

下载

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

相关专题

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

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

2023.07.20

1751

4

python能做什么
python能做什么

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

2023.07.25

4544

7

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

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

2023.07.31

1729

3

python教程
python教程

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

2023.08.03

26737

23

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

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

2023.08.04

3247

5

python eval
python eval

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

2023.08.04

3247

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

2523

5

热门下载

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

精品课程

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

共0课时 | 0人学习

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

共0课时 | 0人学习

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

共0课时 | 0人学习