Olverae Assets

1 把 eˣ 里的 x 换成一个矩阵

你熟悉的 eˣ = 1 + x + x²/2! + x³/3! + … 是一个无穷级数。神奇的是:把这个 x 整个换成一个方阵 A,级数照样收敛,得到的就是矩阵指数 e^A——它也是一个同样大小的矩阵。这不是文字游戏,而是解微分方程组的钥匙。

矩阵指数 e^A = I + A + A²/2! + A³/3! + …。注意开头是单位矩阵 I(对应数里的 1),每一项都是矩阵,无穷个加起来仍是一个矩阵。

这个级数对任何方阵都收敛(因为 n! 增长压过矩阵幂)。但别按定义硬算无穷项——实际计算靠对角化,见下文。

2 先看标量:级数怎么"长成"指数曲线

理解矩阵版之前,先把标量 eˣ 的级数看个明白。下面拖项数 k:蓝线是真正的 eˣ,橙线是只取前 k 项的多项式近似。项数越多,橙线在越宽的范围内贴住蓝线——这正是"无穷级数定义指数"的含义。

泰勒级数逼近 eˣ

取前 k 项3

在 x=2 处真值 e²≈7.389前 k 项近似 = —误差 = —

只取 1 项就是常数 1,2 项是直线 1+x……项数越多,多项式越能"掰弯"去拟合指数曲线。矩阵指数用的是一模一样的级数,只是把数换成矩阵。

3 实战算法:借对角化一步到位

硬加无穷项不现实。但若 A 能对角化 A = P D P⁻¹,有个漂亮捷径:e^A = P e^D P⁻¹,而 e^D 只是把对角线上每个特征值各取 e^λ!整个无穷级数瞬间塌缩成"对每个特征值取一次指数"。

若 A = PDP⁻¹,则 e^A = P e^D P⁻¹,其中 e^D = diag(e^{λ₁}, e^{λ₂}, …)。

道理:A^k = P D^k P⁻¹,代进级数,中间的 P⁻¹P 成对抵消,P 和 P⁻¹ 提到级数外,里头每条对角线就变成标量级数 Σλ^k/k! = e^λ。例:A=diag(2,−1) → e^A=diag(e², e⁻¹)≈diag(7.39, 0.37)。

对角矩阵最省事

若 A 本就对角,e^A 就是把每个对角元素取指数,完事。这是所有计算的基石。

可对角化:换基再取指

换到特征基(P⁻¹)→ 各轴取 e^λ(e^D)→ 换回来(P)。复杂矩阵的指数被拆成几个标量指数。

缺陷矩阵:用 Jordan

不可对角化时改用 Jordan 型,e^{J} 会冒出 t·eᵗ 这类"多项式×指数"项——正是重根微分方程解的来源。

4 它到底有什么用:解 x′ = Ax

最大用途:线性微分方程组 x′ = Ax 的解就是 x(t) = e^{At} x(0)——和标量方程 x′=ax 解为 x=eᵃᵗx₀ 完全平行!而 e^{At} 的行为由特征值 λ 主宰:看标量 e^{λt} 就懂全局。拖 λ 和 t 体会增长/衰减。

e^{λt}:特征值符号决定命运

特征值 λ0.6

行为 = 增长

一个大坑:e^{A+B} = e^A e^B 只在 AB=BA(可交换)时成立!矩阵一般不交换,所以这个熟悉的指数律对矩阵通常失效。但 e^0=I、d/dt e^{At}=A e^{At}、(e^A)⁻¹=e^{−A} 这些始终成立。

5 试三道(先想,再点)

三道题,点选项即时判;答题进度显示在右侧目录。

一句话带走:

矩阵指数 e^A = I + A + A²/2! + …(把 eˣ 级数里的 x 换成方阵)。实算靠对角化:A=PDP⁻¹ ⇒ e^A=P e^D P⁻¹,e^D 对每个特征值取 e^λ;缺陷矩阵用 Jordan。核心用途:x′=Ax 的解是 x(t)=e^{At}x(0),行为由 e^{λt} 决定(λ>0 爆炸、λ<0 衰减)。注意 e^{A+B}=e^Ae^B 仅当 AB=BA。