求解非线性刚性初值问题隐显线性多步方法误差分析.pdfVIP

2015年 2月 计 算 数 学 第 37卷第 1期 Feb．，2015 MATHEMATICA NUMERICA SINICA Vo1．37．No．1 求解非线性刚性初值问题的隐显线性多步方 法的误差分析冰1) 张根根 易 星 肖爱国 (湘潭大学数学与计算科学学院，湖南湘潭411105) 摘 要 隐显线性多步方法由隐式线性多步方法和显式线性多步法组合而成．本文主要讨论求解满足 单边Lipschitz条件的非线性刚性初值问题和一类奇异摄动初值问题的隐显线性多步方法的误差 分析．最后，由数值例子验证了所获的理论结果的正确性及方法处理这两类问题的有效性． 关键词：刚性问题；奇异摄动问题；隐显多步方法；误差分析 MR (2000)主题分类：65L04，65L06，65Lll，65L20 — 1．引 言 刚性初值问题 [。，n]常出现在控制系统、电子网络、生物学、化学动力学等许多科学工程 领域中，它不适宜用显式方法求解；若考虑用隐式方法求解，则对于刚性微分方程大系统而言 其计算量又很大．为了提高计算效率，常采用隐显方法求解刚性初值问题，即将其右端函数分 裂为刚性部分和非刚性部分，分别采用隐式方法和显式方法离散 [卜，M，]．因此，将隐显方法 用于刚性问题的求解具有重要的理论与应用价值． 关于隐显线性多步方法，文献中主要讨论了它们的阶条件、线性稳定性、强稳定性、高效 算法构造及应用等 I卜。，，。，14,¨]．Ascher等 3【]给出了隐显线性 k一步方法的阶条件，具体构造 了一至四阶隐显线性多步方法，并利用 Schur定理获得了一些线性稳定性结果．in’tHoutIsJ 针对线性 自治系统给出了在一定条件下，隐显线性多步方法的稳定域包含所用隐式方法的稳 定域和所用显式方法稳定域的直积．Hundsdorfer和 Ruuth 【】针对带松弛项的抛物型系统的 隐显线性多步方法，研究了方法的单调性和有界性．Koto_9_研究了隐显线性多步法求解标量 试验延迟微分方程的稳定性．Wang和 Ruuth[14J构造了一至四阶变步长的隐显线性多步方 法，并获得 了一些线性稳定性结果．而隐显线性多步方法包括 了重要的隐显 BDF方法 (即对 刚性部分采用 BDF方法离散)，在文献 (如 1【3，5，9，10，13，14])中已有大量讨论和成功应 用。但是，很少有文献研究隐显线性多步方法求解非线性刚性初值问题的收敛性．本文主要讨 论求解两类典型非线性刚性初值问题的隐显线性多步法的误差分析，并通过数值例子证实所 获理论结果的有效性． 本文在第二节提出了所求解的两类典型刚性问题及由隐式线性多步方法和显式线性多步 法组合而成的一类隐显线性多步法．在第三节，给出了隐显多步方法关于两类典型刚性问题 的误差分析结果及其证明．在第四节，数值算例验证了所获理论结果的正确性及方法处理这 两类问题的有效性． 2013年 12月 28日收到． )基金项 目：国家 自然科学基金项 目(No，教育部创新团队项 目(No．IRTl179)及湖南省教育厅重 点项 目(No．14A146)． 2 计 算 数 学 2．问题类与方法 考虑刚性常微分方程初值问题 』I(0())=厂．(，())+g(，())，∈0[，】， (2．1) = 叩， 叼 ∈R ， 一 这里 为 m维欧氏空问，f：D= 0【，T] 11~ 具有刚性，g：D × _÷ 满足经典 Lipschitz条件 g(t，Y)一g(t，)ll Lily—zll，