Birkhoff意义下Hojman-Urrutia方程的离散变分计算.pdfVIP

第 11卷第3期2013年9月 动 力 学 与控 制 学报 Vo1．11No．3 1672-6553／2013／11(31／1994 JOURNALOFDYNAMICSANDCONTROL Sep 2013 Birkhoff意义下Hojman．Urrutia方程的离散变分计算术 宋端 刘世兴 (1．辽东学院影像物理教研室，丹东 l18001)(2．辽宁大学物理学院，沈阳 110036) 摘要 在Birkhoff框架下，采用离散变分方法研究了非Hamilton系统一Hojman—Urrutia方程的数值解法，并通 过和传统的Runge—Kutta方法进行比较，说明了在Birkhoff 架下研究这类不具有简单辛结构的非 Hamilton 系统可以得到更可靠和精确的数值结果． 关键词 Birkhoff方程 ，Hojman—Urrutia方程， 非Hamilton系统， 离散变分计算 DOI： 10．6052／1672—6553—2013-()60 上发表了一篇文章，在该文中，他们给出了如下的 引言 方程 ]： Hamilton系统具有简单的辛结构，自动满足 自 +岁=0， +Y=0 (1) 伴随性质，可以用来描述耗散可忽略的保守动力学 称为Hojman—Urrutia方程．该方程在关于Lagrange 系统，并在动力学系统的保结构算法研究中具有重 逆问题和 Hamilton逆 问题的研究以及 Birkhoff力 要意义 ]．但是对于本质非 自伴随的动力学系 学的发展过程 中起着重要的作用 ． 统 ，在保持实验室可观测量或动力学函数物理意义 取如下初始条件 ： 不变的情况下，则不能表示为 Hamilton系统．虽然 (0)=0， 戈(0)=2， 通过 Darboux变换，可以实现简单的辛结构表示， Y(0)=1， (0)=1 (2) 但是此时的实验室可观测量或动力学函数已经失 则方程(1)有如下的解析解 ： 去了直接的物理意义，我们称这类不能表示为简单 (t)=sin(t)一COS()+t+1 Hamilton方程 的动力学系统为非Hamilton系统 ． Y(t)=sin(t)+COS() (3) 对于非Hamilton系统，由于不具有简单的辛结构， 方程(1)看似简单，却不能表示成 Lagrange方程或 Hamilton系统的保辛算法已经不再适用，因此需要 Hamilton正则方程的形式，因此其是一个典型的非 找到一种较理想的数值算法来数值求解这类非 Hamilton系统问题，其一阶形式没有辛结构．但是 Hamilton系统的运动方程问题． 在Birkhoff力学框架下，该方程可 以表示为满足 自 本文以Hojman．Urrutia(H．U)方程 为例，在 伴随条件的Birkhoff方程的形式，具有简单的辛结 Birkhoff意义下，将 H—u方程表示为 自治Birkhoff方 构．女口果取 口= ，。=Y，。。=戈，0=夕，币0用 Hojman 程的形式，并采用离散变分方法，给出研究这类非 方法可以构造该系统的Birkhoff表示 ，] R1=Ⅱ +0。 Hamilton系统的数值积分子．通过数值实验说明了 ， 在 Birkhoff框架下研究非Hamil