抱歉,您的浏览器无法访问本站
本页面需要浏览器支持(启用)JavaScript
了解详情 >

Hello World

人よ、幸福に生きろ!

矩阵的指数运算 eBe^{\mathbf{B}} 定义为以下形式的无穷级数:

eB=k=0Bkk!,e^{\mathbf{B}} = \sum_{k=0}^\infty \frac{\mathbf{B}^k}{k!},

证明:

det(eB)=etr(B)\det(e^{\mathbf{B}}) = e^{\operatorname{tr}(\mathbf{B})}


解答如下:
任何方阵 B\mathbf{B} 都可以通过相似变换化为其 Jordan 标准型:

B=PJP1,\mathbf{B} = \mathbf{P} \mathbf{J} \mathbf{P}^{-1},

且矩阵的相似变换不改变行列式的值

单个 Jordan 块 Ji\mathbf{J}_i 的形式为:

Ji=[λi1000λi1000λi1000λi].\mathbf{J}_i = \begin{bmatrix} \lambda_i & 1 & 0 & \cdots & 0 \\ 0 & \lambda_i & 1 & \cdots & 0 \\ \vdots & \vdots & \ddots & \ddots & \vdots \\ 0 & 0 & \cdots & \lambda_i & 1 \\ 0 & 0 & \cdots & 0 & \lambda_i \end{bmatrix}.

我们可以通过矩阵指数的定义和 Jordan 块的特殊结构来推导出 eJie^{\mathbf{J}_i} 的表达式。以下是详细的推导过程:

当计算 Jik\mathbf{J}_i^k 时,超对角线的结构会引入系数(类似于组合数),而对角线上的 λik\lambda_i^k 将按照普通的标量次幂增长。

N=JiλiI\mathbf{N} = \mathbf{J}_i - \lambda_i \mathbf{I},则:

Ji=λiI+N,\mathbf{J}_i = \lambda_i \mathbf{I} + \mathbf{N},

其中 N\mathbf{N} 是严格上三角矩阵,即其对角线为 0,超对角线是 1。

由于严格上三角矩阵 N\mathbf{N} 的性质,当其阶数为 n 时,Nn=0\mathbf{N}^n = 0(零矩阵)。

展开 (λiI+N)k(\lambda_i \mathbf{I} + \mathbf{N})^k

Jik=(λiI+N)k=j=0k(kj)(λiI)kjNj.\mathbf{J}_i^k = (\lambda_i \mathbf{I} + \mathbf{N})^k = \sum_{j=0}^k \binom{k}{j} (\lambda_i \mathbf{I})^{k-j} \mathbf{N}^j.

由于 Nj=0\mathbf{N}^j = 0jnj \geq n,我们只需考虑 j<nj < n的项。

Jik\mathbf{J}_i^k 带入矩阵指数的定义:

eJi=k=0Jikk!=k=01k!j=0k(kj)λikjNj.e^{\mathbf{J}_i} = \sum_{k=0}^\infty \frac{\mathbf{J}_i^k}{k!} = \sum_{k=0}^\infty \frac{1}{k!} \sum_{j=0}^k \binom{k}{j} \lambda_i^{k-j} \mathbf{N}^j.

交换求和顺序:

eJi=j=0Njj!k=jλikj(kj)!.e^{\mathbf{J}_i} = \sum_{j=0}^\infty \frac{\mathbf{N}^j}{j!} \sum_{k=j}^\infty \frac{\lambda_i^{k-j}}{(k-j)!}.

观察内层和式:

k=jλikj(kj)!=m=0λimm!=eλi.\sum_{k=j}^\infty \frac{\lambda_i^{k-j}}{(k-j)!} = \sum_{m=0}^\infty \frac{\lambda_i^m}{m!} = e^{\lambda_i}.

因此:

eJi=eλij=0Njj!.e^{\mathbf{J}_i} = e^{\lambda_i} \sum_{j=0}^\infty \frac{\mathbf{N}^j}{j!}.

由于 N\mathbf{N} 是严格上三角矩阵,其幂次具有明确结构:

  • Nj\mathbf{N}^j 的对角线为零;
  • 超对角线第 jj 条上的元素是 1j!\frac{1}{j!}

因此:

j=0Njj!=[1112!1(n1)!0111(n2)!00110001].\sum_{j=0}^\infty \frac{\mathbf{N}^j}{j!} = \begin{bmatrix} 1 & 1 & \frac{1}{2!} & \cdots & \frac{1}{(n-1)!} \\ 0 & 1 & 1 & \cdots & \frac{1}{(n-2)!} \\ \vdots & \vdots & \ddots & \ddots & \vdots \\ 0 & 0 & \cdots & 1 & 1 \\ 0 & 0 & \cdots & 0 & 1 \end{bmatrix}.

eJi=eλi[1112!1(n1)!0111(n2)!00110001].e^{\mathbf{J}_i} = e^{\lambda_i} \begin{bmatrix} 1 & 1 & \frac{1}{2!} & \cdots & \frac{1}{(n-1)!} \\ 0 & 1 & 1 & \cdots & \frac{1}{(n-2)!} \\ \vdots & \vdots & \ddots & \ddots & \vdots \\ 0 & 0 & \cdots & 1 & 1 \\ 0 & 0 & \cdots & 0 & 1 \end{bmatrix}.

从结果中可以看出,eJie^{\mathbf{J}_i} 的对角线元素都是 eλie^{\lambda_i},而上三角部分为某些幂级数系数,但这些非对角元素不会影响行列式。

因此:

det(eJi)=(eλi)n,\det(e^{\mathbf{J}_i}) = (e^{\lambda_i})^n,

其中 nn 是该 Jordan 块的大小。

对于整个 Jordan 标准型矩阵 J\mathbf{J}eJe^{\mathbf{J}} 由所有 Jordan 块的指数组成:

eJ=[eJ1eJ2].e^{\mathbf{J}} = \begin{bmatrix} e^{\mathbf{J}_1} & & \\ & e^{\mathbf{J}_2} & \\ & & \ddots \end{bmatrix}.

行列式的性质告诉我们:

det(eJ)=idet(eJi).\det(e^{\mathbf{J}}) = \prod_i \det(e^{\mathbf{J}_i}).

每个 eJie^{\mathbf{J}_i} 的行列式为 eλie^{\lambda_i} 的积,因此:

det(eJ)=eλi.\det(e^{\mathbf{J}}) = e^{\sum \lambda_i}.

由于 tr(B)=λi\operatorname{tr}(\mathbf{B}) = \sum \lambda_i,我们得到:

det(eB)=etr(B).\det(e^{\mathbf{B}}) = e^{\operatorname{tr}(\mathbf{B})}.

评论