惯性聚合 高效追踪和阅读你感兴趣的博客、新闻、科技资讯
阅读原文 在惯性聚合中打开

推荐订阅源

Martin Fowler
Martin Fowler
WordPress大学
WordPress大学
月光博客
月光博客
OSCHINA 社区最新新闻
OSCHINA 社区最新新闻
大猫的无限游戏
大猫的无限游戏
让小产品的独立变现更简单 - ezindie.com
让小产品的独立变现更简单 - ezindie.com
博客园 - 聂微东
Apple Machine Learning Research
Apple Machine Learning Research
钛媒体:引领未来商业与生活新知
钛媒体:引领未来商业与生活新知
雷峰网
雷峰网
小众软件
小众软件
酷 壳 – CoolShell
酷 壳 – CoolShell
博客园 - 叶小钗
美团技术团队
宝玉的分享
宝玉的分享
Hugging Face - Blog
Hugging Face - Blog
阮一峰的网络日志
阮一峰的网络日志
A
About on SuperTechFans
Jina AI
Jina AI
D
Docker
Last Week in AI
Last Week in AI
MongoDB | Blog
MongoDB | Blog
Stack Overflow Blog
Stack Overflow Blog
Microsoft Azure Blog
Microsoft Azure Blog

又见苍岚

COLMAP PatchMatch Stereo 算法详解 事件驱动的状态机框架:从理论到工程实践 Git 在国内网络环境下无法 Push 的排查与修复 —— 配置 Clash 代理 分段五次多项式插值原理详解 路径插值方法深度对比研究 Claude Code 使用指南 OpenClaw 记忆管理与技能创建指南 CBS(Conflict-Based Search)算法详解 A* 算法及其变种详解 OpenClaw 配置多 Agents Windows Powershell 无法加载文件,因为在此系统上禁止运行脚本问题的解决方案 MaxClaw 安装流程 大模型 AI 名词介绍 AList 网盘聚合工具简介 Protobuf 简介与测试 Claude Code 简介以及 GLM 4.7 模型接入 Github 歌词下载工具 163MusicLyrics Python __getattr__ 懒加载 Python TypedDict 机器人仿真平台 Gazebo 安装记录 机器人仿真平台 Gazebo 简介 多机器人路径规划问题(Multi-Agent Path Finding, MAPF)简介 Python exifread 读取修改过的 jpeg 信息错误问题修复 3D 坐标系变换的理解 3D 旋转矩阵基本概念 MongoDB Compass 介绍 Python 环境管理工具 uv Flutter 开发指南 Snipaste 安装下载与黑屏问题解决方案 全局路径规划算法记录
矩阵分解 -2- 特征值分解
Yiwei Zhang · 2022-09-27 · via 又见苍岚

线性代数中,特征分解(Eigendecomposition),又称谱分解(Spectral decomposition)是将矩阵分解为由其特征值和特征向量表示的矩阵之积的方法。

定义

  • 线性代数中,特征分解(Eigendecomposition),又称谱分解(Spectral decomposition)是将矩阵分解为由其特征值和特征向量表示的矩阵之积的方法。

基础理论

  • $N$ 维非零向量 $v$ 是 $N×N$ 的矩阵 $A$ 的特征向量,当且仅当下式成立:

$$
{\displaystyle \mathbf {A} \mathbf {v} =\lambda \mathbf {v} }
$$

  • 其中 $λ$ 为一标量,称为 $v$ 对应的特征值。也称 $v$ 为特征值 $λ$ 对应的特征向量。也即特征向量被施以线性变换 $A$ 只会使向量伸长或缩短而其方向不被改变。

  • 由上式可得

$$
{\displaystyle p\left(\lambda \right):=\det \left(\mathbf {A} -\lambda \mathbf {I} \right)=0 }
$$

  • 称多项式 $p(λ)$ 为矩阵 $A$ 的特征多项式。上式亦称为矩阵 $A$ 的特征方程。特征多项式是关于未知数 $λ$ 的 $N$ 次多项式。由代数基本定理,特征方程有 $N$ 个解。这些解的解集也就是特征值的集合,有时也称为“谱”(Spectrum)。
  • 我们可以对多项式 p 进行因式分解,而得到

$$ {\displaystyle p\left(\lambda \right)=(\lambda -\lambda _{1})^{n_{1}}(\lambda -\lambda _{2})^{n_{2}}\cdots (\lambda -\lambda _{k})^{n_{k}}=0 } $$

  • 其中

$$ {\displaystyle \sum \limits _{i=1}^{k}{n_{i}}=N} $$

  • 对每一个特征值 λi ,我们都有下式成立:

$$
{\displaystyle \left(\mathbf {A} -\lambda _{i}\mathbf {I} \right)\mathbf {v} =0 }
$$

  • 对每一个特征方程,都会有 ${\displaystyle m_{i}}( {\displaystyle 1\leq m_{i}\leq n_{i}})$个线性无关的解。这 $m_i$ 个向量与一个特征值 $λ_i$ 相对应。这里,整数 $m_i$ 称为特征值 $λ_i$ 的几何重数,而 $n_i$ 称为代数重数
  • 这里需要注意的是几何重数与代数重数可以相等,但也可以不相等。一种最简单的情况是 $m_i = n_i = 1$。特征向量的极大线性无关向量组中向量的个数可以由所有特征值的几何重数之和来确定。

特征值分解

  • 令 $A$ 是一个 $N×N$ 的方阵,且有 $N$ 个线性独立的特征向量 ${\displaystyle q_{i},,(i=1,\dots ,N)}$ 。这样, $A$ 可以被分解为

$$
{\displaystyle \mathbf {A} =\mathbf {Q} \mathbf {\Lambda } \mathbf {Q} ^{-1}}
$$

  • 其中 $Q$ 是 $N×N$ 方阵,且其第 $i$ 列为 $A$ 的特征向量 $q_i$。

  • $Λ$ 是对角矩阵,其对角线上的元素为对应的特征值,也即 $\Lambda_{ii}=\lambda_i$。这里需要注意只有可对角化矩阵才可以作特征分解

    比如 ${\displaystyle {\begin{bmatrix}1&1\\0&1\\\end{bmatrix}}}$ 不能被对角化,也就不能特征分解。

  • 一般来说,特征向量 $q_i ,, (i = 1, \dots, N)$ 一般被单位化(但这不是必须的)。未被单位化的特征向量组 $v_i ,, (i = 1, \dots, N)$, 也可以作为 $Q$ 的列向量。

对称矩阵

  • 任意的 $N×N$ 实对称矩阵的特征值都是实数且都有 $N$ 个线性无关的特征向量。并且这些特征向量都可以正交单位化而得到一组正交且模为 1 的向量。故实对称矩阵 $A$ 可被分解成

$$
{\displaystyle \mathbf {A} =\mathbf {Q} \mathbf {\Lambda } \mathbf {Q} ^{-1}=\mathbf {Q} \mathbf {\Lambda } \mathbf {Q} ^{T}}
$$

  • 其中 $Q$ 为 正交矩阵, $Λ$ 为实对角矩阵。

证明 $\mathbf Q$ 正交

  • $A$ 为 $N×N$ 实对称矩阵,特征值都是实数且都有 $N$ 个线性无关的特征向量。
  • 则 $A$ 可分解为:

$$
\mathbf {A} =\mathbf {Q} \mathbf {\Lambda } \mathbf {Q} ^{-1}
$$

  • 转置上式:

$$
\mathbf {A}^T =(\mathbf {Q} {-1})T \mathbf {\Lambda } \mathbf {Q}^T
$$

  • 由于 $A$ 为对称矩阵:

$$
A^T=A
$$

  • 因此对应位置的分解矩阵完全相同,有:

$$
\mathbf {Q} ^{-1}=\mathbf {Q} ^{T}
$$

  • 因此 $\mathbf Q$ 为 正交矩阵

通过特征分解求反(逆)矩阵

  • 若矩阵 $A$ 可被特征分解并特征值中不含零,则矩阵 $A$ 为非奇异矩阵,且其逆矩阵可以由下式给出:

$$
{\displaystyle \mathbf {A} ^{-1}=\mathbf {Q} \mathbf {\Lambda } ^{-1}\mathbf {Q} ^{-1}}
$$

  • 因为 $ Λ$ 为对角矩阵,其逆矩阵容易计算出:

$$
{\displaystyle \left[\Lambda ^{-1}\right]_{ii}={\frac {1}{\lambda _{i}}}}
$$

参考资料

文章链接:
https://www.zywvvd.com/notes/study/linear-algebra/matrix-decomp/eigenvalue-deco/eigendecomposition/