矩阵指数

✍ dations ◷ 2024-12-27 00:07:49 #矩阵论,李群,指数

矩阵指数是方块矩阵的一种矩阵函数,与指数函数类似。矩阵指数给出了矩阵李代数与对应的李群之间的关系。

设为×的实数或复数矩阵。的指数,用或exp()来表示,是由以下幂级数所给出的×矩阵:

以上的级数总是收敛的,因此的指数是定义良好的。注意,如果是1×1的矩阵,则的矩阵指数就是由的元素的指数所组成的1×1矩阵。

设和为×的复数矩阵,并设和为任意的复数。我们把×的单位矩阵记为,把零矩阵记为0。

我们可以从指数级数的定义直接得到矩阵指数的如下性质:

接下来是一个关键性质:

由此导出的推论有:

矩阵指数的一个重要性,是它可以用来解微分方程。从(1)可知,以下微分方程

其中是矩阵,具有解

矩阵指数也可以用来解非齐次方程:

参见以下的例子。

当不是常数时,以下形式的微分方程没有闭式解:

但马格努斯级数可以给出无穷级数形式的解。

根据雅可比公式,对任意复矩阵,下列迹等式成立:

det ( e A ) = e tr ( A )   {\displaystyle \det(e^{A})=e^{\operatorname {tr} (A)}~} ) ≠ 0,从而必可逆。

我们知道,对于任何实数(标量)和,指数函数都满足公式 + = 。类似的等式对于可交换矩阵也成立:如果矩阵和是可交换的(即 = ),则:

但是,如果它们不是可交换的,则以上的等式不一定成立。

这个命题反过来不成立:+=并不一定就意味着和是可交换的。但是,如果和只含有代数数,而且它们的大小至少为2×2,则反过来也成立。

和不可交换的情况可以用以下方法计算:

即使 X {\displaystyle X} 的逆矩阵由−给出。这与复数的指数总是非零的事实类似。这样,矩阵指数就给出了一个映射:

这是从所有×矩阵的空间到一般线性群(所有非奇异方阵所组成的群)的映射。实际上,这个映射是满射,就是说每一个非奇异方阵都可以写成某个矩阵的指数。矩阵对数就是这个映射的逆映射。

对于任何两个矩阵和,我们有:

其中|| · ||表示任意的矩阵范数。从中可以推出,指数映射在(C)的紧子集内是连续和利普希茨连续的。

以下的映射

定义了一般线性群中的一条光滑曲线,当 = 0时穿过单位元。实际上,这给出了一般线性群的一个单参数子群,这是由于:

这条曲线在点的导数(或切向量)由以下等式给出:

= 0时的导数就是矩阵,所以我们可以说,是这个单参数子群的推广。

更加一般地:

寻找计算矩阵指数的可靠和准确的方法是困难的,目前在数学和数值分析领域中仍然是一个正在研究的话题。有些方法列举如下。

如果矩阵是对角的:

则把主对角线上的所有元素取指数,就是原矩阵的指数:

这也允许了我们计算可对角化矩阵的指数。如果 A = U D U 1 {\displaystyle A=UDU^{-1}} 是对角矩阵,则 e A = U e D U 1 {\displaystyle e^{A}=Ue^{D}U^{-1}} ,有 = 0,则矩阵称为幂零矩阵。在这种情况下,矩阵指数可以直接从级数展开式来计算,这是因为级数在有限个项后就终止了:

当矩阵的最小多项式可以分解为一次多项式的积时,它就可以表示为以下的和:

其中:

这称为Dunford分解。

这就是说,我们可以通过化为前两种情况,来计算的指数:

注意为了让最后一步成立, 和必须是可交换的。

另外一个密切相关的方法,是利用的若尔当标准型。假设 =  −1,其中是的若尔当标准型。那么:

另外,由于

因此,我们只需要知道怎样计算若尔当块的矩阵指数。但是,每一个若尔当块都具有形式

其中是幂零矩阵。则这个区块的矩阵指数由下式给出:

假设我们想要计算以下矩阵的指数。

它的若尔当型为:

其中矩阵由下式给出:

我们首先来计算exp()。我们有:

1×1矩阵的指数仅仅是该矩阵的元素的指数,因此exp(1(4)) = 。 J 2 ( 16 ) {\displaystyle J_{2}(16)} ) = λ exp()来算出:

因此,原矩阵的指数为:

矩阵指数在解线性微分方程时十分有用。前面曾提到,以下形式的微分方程

具有解C。如果我们考虑以下向量

我们就可以把线性微分方程表示为:

如果我们作一个猜想,把两边乘以一个积分因子 −,便得到:

如果我们可以计算,那么就得到了微分方程的解。

假设我们有以下的微分方程组:

相关的矩阵为:

在以上的例子中,我们计算了矩阵指数

因此微分方程组的通解为:

也就是说,

对于非齐次的情况,我们可以用积分因子的方法(类似于参数变换的方法)。我们找到形为p() = exp()()一个特解:

为了让p为方程的解,必须有:

因此,

其中由问题的初始条件决定。

假设我们有以下的微分方程组:

那么我们有

以及

用前面的方法,我们可以得出齐次微分方程的解。由于齐次方程的通解与非齐次方程的特解的和就是非齐次方程的通解,因此我们只需要找到一个特解(用参数变换法)。

我们有:

进一步简化,就可以得到原方程的特解。

相关

  • 负回授负反馈(英语:negative feedback),是反馈的一种。是指系统的输出会影响系统的输入,在输出变动时,所造成的影响恰和原来变动的趋势相反;反之,就称为正反馈。另一种说法是系统在一个条
  • 楚瓦什共和国楚瓦什共和国(俄语:Чува́шская Респу́блика — Чува́шия, Chuvashskaya Respublika — Chuvashiya; 楚瓦什语:Чăваш Республик
  • 行星适居性行星适居性是天文学里对星体上生命的出现与繁衍潜力的评估指标,其可以适用于行星及行星的天然卫星。生命的必要条件是能量来源(通常是太阳能但并不全然)。但通常是当其他众多条
  • 1157年重要事件及趋势重要人物
  • 秃头脱发,又称为脱发、掉发,指的是毛发从头部或是身体上掉落,一般来说至少会包含头发的脱落,其严重程度可以从一个小范围,乃至全身受影响。一般来说不会伴随发炎或出现疤痕。正常人的
  • span class=nowrapAgFsub2/sub/span二氟化银(化学式:AgF2),IUPAC名称氟化银(II),是少见的银(II)化合物之一。二氟化银是无色或白色结晶粉末,对光敏感,有杂质氧化银和氟化银存在时呈棕色或黄色。它是强氟化剂,需保存在
  • 保罗·塞尔瓦保罗·约瑟夫·塞尔瓦(英语:Paul Joseph Selva),美国空军上将,现为第10任美国参谋长联席会议副主席,之前曾任美国运输司令部司令,2015年7月31日就任现职。
  • 卢博米尔·胡萨尔卢博米尔·胡萨尔(乌克兰语:Любомир Гузар、英语:Lubomyr Husar:1933年2月26日-2017年5月31日),是乌克兰籍天主教司铎级枢机。也是乌克兰希腊礼天主教会基辅-加利奇荣
  • 埃利普索伊德山坐标:77°48′S 163°50′E / 77.800°S 163.833°E / -77.800; 163.833埃利普索伊德山(英语:Ellipsoid Hill)是南极洲的山峰,位于维多利亚地的斯科特海岸,在布卢冰川以北,处于斯菲
  • 哥斯拉 (1954年)《哥斯拉》(日语:ゴジラ)是1954年11月3日首映的日本电影,为哥斯拉系列电影的第一部作品,片长97分钟,为黑白摄制。当时并映作品为古装喜剧电影《仇讨珍剑法》(日文原题:仇討珍剣法)。