二阶常微分方程初值问题数值方法的研究综述.docVIP

二阶常微分方程初值问题数值方法的研究综述 摘 要 ： 二阶常微分方程周期初值问题数值方法近些年来倍受人们的关注。我们将对近些年来二阶常微分方程周期初值问题数值解的研究做一个简要的综述, 并提出进一步研究的设想。 关键词 ： 周期初值问题 ； 数值方法； 代数精度； 周期性区间； P 一稳定性； 相延迟阶 1 、 引言 考虑二阶常微分方程初值问题 里假设 有足够的 滑性， 且不含 Y的一阶导数项， 并假设这一问有周期解或振荡解。这一问题的数值方法近些来备受人们的关注，因为在科学与工程的许多域中都能碰到， 比如天体力学， 量子力学【 1 ， 2 】 论物理与化学， 电子学以及偏微分方程的半离等等。例如， 考虑一维S c h r ~ d i n g e r 方程 其中l 是一个非负整数， W( x )是势能函数，E是常数。这一方程的数值解近二十年来引起人广泛的兴趣，出现了许许多多的数值方法[ 3 ] 。易见， 这是二阶周期问题数值解的一个特款。本文我们将对近些年来二阶常微分方程周期初值问题数值解的研究做一个简要的综述， 共分以下几个部分。在第二节中， 我们给出这一问题的一些基本概念与基本理论；在第三节中， 我们将结合自己的研究工作对解这一问题的数值方法进行一些综述； 最后在评 论中， 我们提 出进一步研究的设想。 2 、 基本理论 考虑有周期解或振荡解的二阶常微分方程初值问题 ( 1 )。在设计数值方法时， 方法的代数 精度通常是考虑的首要因素， 这是因为代数精度 愈高， 我们得到的数值解就愈精确。 定义 1 .[ 4 ] 考虑解问题 ( 1 )的线性 l 步方法 其中 是系数且 作为 的 近似值，如果方法 ( 2 ) 对充分可微函数 有下式满足 则称方法 ( 2 ) 有 q阶代数精度。下面考虑单步方法的代数精度。将( 1 )化成与其等价的一阶系统 定义2. [ 5 ] 考虑解问题 ( 1 )的单步方法， 如果对充分可微函数 有下式满足 其中) 为精确解， 为数值解， 则称此单步方法有r 阶代数精度。同时，由于问题的精确解是周期的或振 荡 的， 故我们希望数值解也应有类似的性质， 这样就产 生 了数值 方 法 的稳 定性 问题 。 1 9 7 6年 ，L a mb e r t ， Wa t s o n在讨论解问题 ( 1 )的线性对称多步法的稳定性时， 基于纯量模型方程。 提出了两个非常重要的概念： 周期性区间与 P _ 稳定性f 4 ] 。 当某个线性对称多步法用于 ( 4 ) 时，我们可以得该方法的特征多项式，如果存在一个实区问使得对所有的 h 是时问步长，特征多项式有两个共轭复根 ( 称它们为主根 ( p r i n - c i p a l r o o t s )) 位于单位圆周上， 而其余的根 ( 称它们为伪根 ( s p u r i o u s r o o t s )) 均在单位圆内， 则这个区间就是该方法的周期性区问。如果周期性区问是 则称该方法是 P 一 稳定的。这意味着， 如果选取一个步长 h使得 ( h h ) 位于周期性区问中， 则由主根引起的扰动对解析解是一个周期逼近， 而相对而言， 由伪根引起的扰动则可以忽略不计。从几何上讲， 此时数值解会停 留在周期轨道上， 而不会内旋或外旋， 有时称这种现象为轨道稳定的 ( o r b i t a l l y s t a b l e )。P 一 稳定性则保证了无论步长取多大， 数值解的周期性均能保持。此外 ， 方法 的相延迟 ( p h a s e — l a g ) 性质也是一个考 虑的重要因素。相延迟是指精确解与数值解之 间 的夹角， 它反映了两者的某种同步关系。这一概念最早 由 B r u s a ， N i g r o在讨论线性结构动力系统 的数值方法时提出的【 6 ] 。后来， G a l d we l l ， T h o ma s ，d e r Ho u w e n ， S o mme r i j e r ， C o l e ma n ， S i mo s等 人 详细讨论一些方法的相延迟性质， 提出了相延迟阶 的概念[ 7 - 1 1 ， 3 ] 。为了便于说明， 我们下面给出单步方法与两步方法的周期性区间， P 一 稳定性以及 相延迟阶的定义与一些性质。考虑单步方法。将某一单步法运用至模型方 程 ( 4 )， 则我们得到它的特征方程