跳转到内容

马格努斯展开

维基百科,自由的百科全书

数学物理学中,马格努斯展开(英语:Magnus expansion)为线性算子的一阶齐次线性微分方程的解提供了指数表示,得名于数学家威廉·马格努斯。特别地,这种方法提供了变系数n阶线性常微分方程组的基矩阵。指数是无穷级数,其项涉及多重积分和嵌套换元。

确定情形

[编辑]

马格努斯方法及其解释

[编辑]

给定n × n系数矩阵A(t),我们希望求解与线性常微分方程相关的初值问题

其中Y(t)是未知n维向量函数。

n = 1时,解为

A(t)对任意一组t, t1t2仍满足A(t1) A(t2) = A(t2) A(t1),则此式可推广到n > 1情形。At无关时尤为如此。在一般情形下,上述表达式不再是问题的解。

马格努斯提出的解矩阵初值问题的方法是,用某个n阶方阵函数Ω(t, t0)的指数来表示解:

稍后可将其构造为级数展开式:

为简单起见习惯将Ω(t, t0)写作Ω(t),并取t0 = 0.

马格努斯意识到,由于d/dt (eΩ) e−Ω = A(t),可利用庞加莱-豪斯多夫矩阵恒等式将Ω的时间导数和伯努利数及 Ω伴随自同态联系起来

并以“BCH展开的连续类似物”递归求解Ω,下详。

上式构成了矩阵线性初值问题求解的马格努斯展开式马格努斯级数。前4项:

其中[A, B] ≡ A BB AAB的矩阵交换子

这些方程可解释如下:Ω1(t)与标量(n = 1)情形下的指数完全重合,但这方程无法给出整个解。若坚持要用指数表示(李群),则要对指数进行修正。马格努斯级数的剩余部分系统地提供了修正:Ω或其部分在解的李群李代数中。

在应用中,很少能对马格努斯级数精确求和,而要截断才能得到近似解。马格努斯方法的主要优势在于,中截级数通常和精确解具有相同的重要性质,这异于传统摄动理论。例如,经典力学中,时间演化辛几何特征在每阶近似都得到保留。同样,量子力学时间演化算子的幺正性也得到保留(例如,与解决同一问题的戴森级数相反)。

扩展的收敛性

[编辑]

从数学角度看,收敛问题如下:给定矩阵A(t),何时可得作为马格努斯级数和的指数Ω(t)

t ∈ [0,T)时,级数收敛的充分条件是

其中表示矩阵范数。这个结果是通用的,因为可构造特定矩阵,t > T时级数都发散。

马格努斯生成器

[编辑]

生成马格努斯展开式中所有项的递归过程利用了下面的递归定义的矩阵 Sn(k)

然后得到

此处adkΩ是迭代交换子的简写(参见伴随自同态):

其中Bj伯努利数,而B1 = −1/2

最后,明确算得这一递归后,就可将Ωn(t)表为涉及n个矩阵A的n-1个嵌套换元的n重积分的线性组合:

随着n增加,这个式子会变复杂。

随机情形

[编辑]

推广到随机常微分方程

[编辑]

要推广到随机常微分方程,令布朗运动概率空间 上,有限时间区间和自然过滤。现在,考虑线性矩阵值随机伊藤积分方程(索引j采用爱因斯坦求和约定)

其中是逐步可测的值有界随机过程单位矩阵。参考确定情形,并依随机情形做修改[1],相应的矩阵对数将变为伊藤过程,其展开的前两项为, 其中ij根据爱因斯坦求和约定

推广的收敛性

[编辑]

随机情形下,收敛将受制于停止时间,第一个收敛结果如下:[2] 在前面关于系数的假设下,存在强解及严格为正的停止时间,使得:

  1. 在时间之前有实数对数,即
  2. 以下表示有把握成立:

    其中是随机马格努斯展开的第n项,定义见下文马格努斯展开式小节;
  3. 存在正常数C,仅取决于,其中,于是

马格努斯展开式

[编辑]

随机马格努斯展开的推广形式:

其中通用项是形式为下式的伊藤过程:

项可递归定义为

其中

算子S定义为

应用

[编辑]

1960年代以来,马格努斯展开作为一种摄动理论工具,已成功应用于物理学化学的许多领域,从原子物理学分子物理学核磁共振[3]量子电动力学。自1998年以来,它还用于构建矩阵线性微分方程数值积分的实用算法。由于它们集成了马格努斯展开的特性,保留了问题的定性特征,因此相应方案也是几何积分的典型例子。

另见

[编辑]

注释

[编辑]
  1. ^ Kamm, Pagliarani & Pascucci 2021
  2. ^ Kamm, Pagliarani & Pascucci 2021,Theorem 1.1
  3. ^ Haeberlen, U.; Waugh, J.S. Coherent Averaging Effects in Magnetic Resonance. Phys. Rev. 1968, 175 (2): 453–467. Bibcode:1968PhRv..175..453H. doi:10.1103/PhysRev.175.453. 

参考文献

[编辑]