Hamilton力学的辛算法

哈密尔顿系统的辛几何算法哈密尔顿系统是一类具有特殊的物理意义的动力系统，其在物理学、力学、动力学和计算力学等领域有着广泛应用。

哈密尔顿系统通常具有一组关于位置和动量的相变量，其演化满足哈密尔顿方程。

由于哈密尔顿系统具有良好的保持量和结构稳定性，因此在数值模拟中的算法设计尤其重要。

辛几何算法是一类特殊的数值演化方法，其以保持哈密尔顿系统相变量守恒和辛结构稳定性为目标，常常用于哈密尔顿系统的数值积分。

辛几何算法最早由李约瑟于 1988 年提出，其不仅能够在数值计算中保持相变量的守恒，还能够在哈密尔顿系统的长期演化中保持辛结构稳定性。

辛几何算法主要由两个部分组成，即辛映射和辛算子。

辛映射指的是从一个相变量向下一个相变量的映射，它通常满足“保相量”和辛结构不变性的特点。

保相量指的是相变量在变化过程中的守恒，而辛结构不变性则指的是哈密尔顿系统在演化过程中的辛不变性。

而辛算子则是这个辛映射的数值逼近，常常采用辛波发方法、显式和隐式辛算法等方法。

在演化哈密尔顿系统时，辛几何算法通常采用显式辛算法进行数值模拟。

显式辛算法的主要思路是采用辛映射和辛算子的组合，来实现对哈密尔顿系统的数值模拟。

在模拟过程中，辛几何算法需要保证每一步的演化都是辛的，这样系统才能保持哈密尔顿量以及其他相变量不变。

因此，辛几何算法在数值模拟中的应用非常广泛。

然而，辛几何算法的实现却比较困难。

在数值模拟时，辛几何算法需要考虑一系列问题，如相变量的数值守恒、哈密尔顿量的捕获和重构、快速演化、长时间演化、难以计算的高维效应等等。

这些问题都需要采用一些特殊的技巧和策略来解决。

总的来说，辛几何算法是一种特殊的数值演化方法，其以保持哈密尔顿系统相变量守恒和辛结构稳定性为目标，在计算力学、物理学、动力学等领域有着广泛应用。

辛算法研究冯康教授于1984年提出Hamilton系统的辛几何算法，首次将保持Hamilton 系统几何结构的思想引入数值分析，随后引来了国内外在这方面的极大兴趣。

冯康教授开辟了一个新的研究领域，并在近十年的时间里带领他的研究小组在辛几何算法等保结构算法及其在数值分析中的应用等方面进行了广泛深入的研究，取得了丰硕的成果。

现在人们已越来越意识到保结构算法的重要性，事实上保结构算法已在很多领域包括天体力学，量子化学，非线性波，不可压流体，大气物理和地物勘探数据处理等，找到很好的应用。

另一方面，在这些领域中的应用反过来又必将促进辛几何算法等保结构算法本身的不断完善何不断发展。

大量数值实验结果已证明，较之传统算法，保结构算法对保结构动力系统的计算具有令人信服的优势。

在辛几何算法理论方面，仍有一些悬而未决的难点问题和最新提出的一些重点问题有待解决；另一方面，在辛算法的基础上，现在国际上又兴起了专门解保守型偏微分方程的多辛算法。

孤立子波动方程是广泛应用于物理领域的非常重要的守恒型偏微分方程，对其复杂孤立子波的数值模拟二十多年来一直是国际上一个热门课题，数值结果已经显示辛算法，多辛算法的优势和可行性，但至尽仍然存在许多有待解决的问题。

我们主要研究的是多辛算法，现阶段的研究水平已完全与国际科研水平持平，并保持同步发展。

我们把多辛几何算法应用到孤子方程中去，数值研究表明多辛算法和常微分方程情形一样具有长期跟踪能力, 不会带来人为污染, 能正确反应孤子碰撞问题.我们发表在Phys.A Mathematical gen(2000), 33:18, 3613--3626上的论文给出了一个三层12点格式. 此文发表后不久, 两美国学者Schultz M ,Trimper S 在此刊同卷41期上发表论文:“动力运动产生孤子”, 文中称我们的方法是著名的方法。

李群算法是最近才发展起来的一种很有发展前途的数值算法。

它用来求解齐次流形上的常微分方程，使得所求的数值解仍在同一流形上李群算法是冯康先生辛几何算法的拓扩，把保几何结构的思想推广到一般李群上。

辛算法在电磁计算中的应用摘要近几年，随着计算机性能的飞速发展和计算物理中各种新型算法的出现，各种电磁场数值方法层出不穷，但很多算法面临着计算时间长、储存空间不足及计算精度低等方面的困难。

Hamilton系统理论是当代数学物理中的一个重要的工具。

一切守恒的物理过程，总能表示成适当的Hamilton系统。

辛算法正是保持Hamilton系统内在性质的一种新型数值方法，该算法在长时间的数值计算中，具有一般数值方法无可比拟的计算优势。

本文首先介绍了电磁学的基本背景和电磁计算的研究，然后介绍了辛算法。

接着，介绍了辛算法在Maxwell方程中的应用，然后在无耗煤质和散射存在时的情况下分析了辛时域有限差分法的计算式。

最后，以真空中一维的高斯脉冲电磁波为例用辛算法进行了数值运算。

关键词：电磁计算；辛算法；Hamilton系统；Maxwell方程一．引言电磁场理论的应用遍及地理学、生命科学、医学、材料科学和信息科学等几乎所有技术学科领域。

计算电磁学是以电磁场理论为基础，以高性能的计算技术为手段，运用计算数学提供的各种方法，解决复杂电磁场理论和工程问题的应用科学。

因此，开展计算电磁学的研究不仅可以产生国际水平的基础研究成果，更重要的是可以促进我国民用和军用电磁学相关领域的发展。

早在1864年，Maxwell在前人理论和实验的基础上建立了统一的电磁场理论，并用数学模型揭示了自然界一切宏观电磁现象所遵循的普遍规律，这就是Maxwell方程组，它包括微分形式和积分形式。

简单地说，所有的宏观电磁问题都可以归结为Maxwell方程组在各种边界条件下的求解问题。

计算电磁学自20世纪60年代兴起，至今40余年。

纵观整个电磁理论发展的过程，电磁学的发展可以分为两个阶段。

以20世纪60年代为分界点，之前可以称为经典电磁学阶段，在这个时期，电磁场理论和工程中的许多问题大多采用解析或渐进的方法进行处理，即在几种可分离变量的坐标系中求解Maxwell方程组或其退化形式，最后得到解析解。

哈密尔顿系统有限元的守恒性和辛性质哈密尔顿系统是最重要的动力系统。

冯康院士曾指出，一切真实的无耗散的物理过程都可表示为这样或那样的Hamilton形式，它们都是常微分或偏微分方程组。

Hamilton系统有两个最重要的特性：守恒性和辛结构。

在数值计算中能否保持这些特性具有重要意义。

1983年冯康研究此问题，他惊讶地发现，此前这里竟是一片空白，许多经典的算法都不适应，计算几万步后有时已面目全非。

他1984年首创性提出辛差分算法，并作了深入系统的研究，开辟了一大片研究新领域。

以后国内外许多学者作了多方面推广和广泛应用。

冯康的这项首创工作得到了国际一致公认。

但是任何离散算法，一般不可能同时保辛又保能量（Ge-Masden定理）。

辛差分算法很好地保辛，但只在格式精度意义下保能量。

而许多学者认为，有时保能量更重要。

因此我们转向有限元法，却发现至今有关研究极少。

而我们的研究表明，有限元总是保能量的，对线性系统也是辛的，对非线性系统是高精度保辛的，而且长时间计算稳定且精度高，效果非常好。

这些结论已刻划了有限元的基本特征。

因此有限元法是与辛差分算法完全不同的另一种算法，从另一方面弥补了辛差分算法的不足。

本研究是对辛算法的一次重要推进。

本文主要创新点如下：（1）．首次系统深入研究任意m次有限元解非线性Hamilton系统，证明了在任何节点上能量总是守恒的，因此在相平面上轨道总是稳定的。

并首次提出用超收敛分析方法研究有限元的辛性质；（2）．对线性Hamilton系统的任意m次有限元，得到了一个深刻的高阶超收敛O(h^2m+1)新估计，首次证明m次有限元的节点值是2m阶对角Páde 逼近，因而是辛格式。

此结果与冯康等研究的辛差分格式结论一致；（3）．对非线性Hamilton系统的任意m次有限元法，构造了新的辅助问题，并得到误差估计和负范数估计，首次证明m次有限元对每一次步进是高精度O(h^2m+1)意义下近似保辛的。