单刚性SPDIDE系统线性多当方法的数值收敛性.pdfVIP

2010年2月 湖北第二师范学院学报 Feb．2010 第27卷第2期 JournalofHubeiUniversityofEducation V01．27 No．2 单刚性 SPDIDE系统线性麦 方法的数值收敛性 周丹丹 赵 然 1湖北第二师范学院数学与数量经济学院，武汉430205；2华中科技大学数学系，武汉430074 摘 要：讨论数值求解单刚性奇异摄动延迟积分微分系数的误差分析和数值收敛性。利用离散变分方法获得 了A 稳定的线性多步法关于这类时滞奇异摄动问题的收敛性结果。数值例子进一步证实了理论结果的正确性。 关键词：时滞奇异摄动；线性多步方法；收敛性 中图分类号：o175．6 文献标识码 ：A 文章编号：1674．344X 2010 2-0012-03 作者简介：周丹丹 1982一 ，女，湖北钟祥人，助教，硕士，研究方向为分布参数系统的控制ntCe-~数值算法的稳定性。 奇异摄动初值问题是一类特殊刚性问题，其特殊 令r ， p 一 一 · 一 一 ·一 的结构特点使得这类问题不能被 B理论满意地覆盖。 r ， 5 这一问题主要包括单刚性和多刚性两大类。目前，对 其中rj xi 是离散变分式 5 中的各项系数。 于最简单形式的单参数奇异摄动初值问题，国内外已 rm 由 2 4 式可得：e - i ‘ 6 有一些数值方法收敛性研究。Hairer[1】，Lubichp[引， ． Strehmel ，Schneider引 Kadalbajoo ，肖爱国 】，侯 故若系数 rm-j xi 存在恰当的界，则 6 式可以用 ， 文星¨0]，甘四清和孙耿 ¨-13]等人基于直接法或8渐 来估计误差项em的界。特别的，系数 实际上就 进展开法分别给出了关于线性多步法，Runge—Kutta 是用线性多步方法 p， 以步长h 1求解齐次线性 方法，Rosenbrock型方法，以及一般线性方法的收敛性 方程Y m · ，， 0得到的各个数值解 ，，，其中取定 结果。 初始启动值为Y一+l … Y一， 0，Yo a一 I ～． 而对于积分延迟奇异摄动问题 【1数值方法的收 则在假设 一 ≠0的情况下， 有恰当的界存 敛性性质尚少见相关报道。本文在前人工作的基础 在。故通过复合叠加等运算可得到整体误差e的界。 上，研究了线性多少方法救解一类单刚性奇异摄动延 2 线性多步算法的收敛性 迟积分微分 SPDIDE 系统的数值收敛性。 考虑SPDIDE系统 I 离散变分公式方法 ， ·Y ，‘Y t ，Y t一下 ， 离散变分公式方法被普遍用于解决有关奇异摄动 』：一k ￡，s，y s ds ，tE[，]， 强刚性系统数值方法的收敛性研究。 I．y ￡ ￡ ，￡≤￡0，占 O 考虑具有如下形式的误差方程： 其中奇异摄量 远小于 1．f，k均为已知光滑函 Otf一 i e dm m≥0 1 数，满足经典Lipschitz条件： 为了找出恰当的条件使得整体误差 e有界，首 lIf t