脉冲微分方程的一个修正block-by-block数值格式-论文.pdfVIP

第33卷 第5期 工 程 数 学 学 报 Vol. 33 No. 5 2016年10月 Oct. 2016 文章编号 脉冲微分方程的一个修正 数值格式∗ 曹俊英 , 马群长 , 王自强 贵州民族大学理学院，贵阳 北京大学数学科学学院，北京 摘 要: 本文利用修正的block-by-block 方法针对脉冲微分方程构造了高阶数值格式．修正 的block-by-block 方法是传统的block-by-block 方法的改进，其优点是除第一块外其余 每块都能够解耦求解积分方程的高阶数值方法．首先，把脉冲微分方程等价转化为脉 冲型积分方程，并利用修正的block-by-block 方法进行离散，得到在两个相邻脉冲点 中除第一块外其余每块都解耦的高阶数值格式．其次，利用离散的Grownwall 不等式 证明了数值解逼近精度为四阶．最后，一系列的数值算例验证了理论分析的正确性． 关键词: 脉冲微分方程；block-by-block 算法；高阶数值格式；收敛性分析 分类号: (2000) 34K45; 65N15 中图分类号: O241.82 文献标识码: A 1 引引引言言言 随着科学技术的不断发展，许多反映客观现实的物理模型都具有这样的特征：在发展 的某些阶段，由于受到外部的作用或系统内部自身的原因，使得系统瞬间改变原有的状 态，这种瞬间改变原有状态的现象通常称为脉冲现象．脉冲现象在现代科技各领域的实际 问题中是广泛存在的，其数学模型往往可归结为脉冲微分系统．脉冲微分系统能够更深 刻、更精确地反映事物的变化规律．随着科学技术的突飞猛进，人们越来越认识到脉冲微 分系统的重要性以及在实践中的应用价值．例如，具有脉冲影响系统应用在物理、生物科 技、工业、机器人、辐射能、药物动力学、种群动态、生态、最优控制、微生物的研究、 繁殖、经济学、生产理论等方面．正是由于脉冲微分系统能够更好地反映瞬间突变事物的 变化规律．因此，它已经吸引了国内外众多的专家和学者对脉冲微分方程的理论进行了研 - 究 ． 由于脉冲微分方程起源较早，研究脉冲微分方程的专家和学者已经从各个方面对 脉冲微分方程进行了研究，从而使得脉冲微分方程理论取得了丰富的成果．但是对于 脉冲微分方程数值方法的研究还刚刚起步，关于脉冲微分方程数值方法的成果还较 少．2000 年，Ranbelovi´c 等人 基于著名的Runge-Kutta 数值方法提出了一种计算脉冲 微分方程的新算法．2005 年，Hosseini 利用修正的谱方法求解具有脉冲项的常微分 方程．2007 年，Liu 等人 研究了一个线性脉冲微分方程，并指出若将Runge-Kutta 方 法应用于该系统，那么该方法是数值稳定的．2008 年，Ran 等人 主要研究了线性 收稿日期: 2015-07-30. 作者简介: 曹俊英(1981年5月生)，女，博士，教授. 研究方向：微分方程数值解. 基金项目: 国家自然科学基金；贵州省科学技术基金([2014]2098; [2013]2144). 第5期 曹俊英，等：脉冲微分方程的一个修正block-by-block 数值格式 507 脉冲微分系统的Runge-Kutta 方法的数值解的稳定性．2008 年，Liang 和Liu 主要 针对一类非线性的脉冲微分方程进行研究，在 分别满足单边Lipschitz 条件和经 典Lipschitz 条件下得到了其解析解及数值解的稳定