基于抛物线逼近法方法的常微分方程数值解法研究.pdfVIP

第 29卷 第2期 河 南 科 学 V01．29 No．2 2011年 2月 HENAN SCIENCE Feb．2011 文章编号：1004—3918(2011)02—0127—06 基于抛物线逼近法方法的常微分方程数值解法研究 杨小远 ， 张 凯2， 唐文琦27 魏 华2 (1．北京航空航天大学 数学与系统科学学院，数学、信息、行为教育部重点实验室，北京 100191； 2．北京航空航天大学 计算机学院，北京 100191) 摘 要：通过对欧拉法的近似方法进行分析，采用抛物线与坐标轴所围面积作为积分值的近似，从而获得 了比欧拉 法精确度更高并且可以涵盖梯形法以及龙格一库塔法的一种常微分方程数值解的方法——抛物线逼近法，并就其局 部误差，收敛性及稳定性进行 了分析，最后给出了该方法的数值例子． 关键词：抛物线；辛 卜生公式；泰勒展式 中图分类号：0177．2 文献标识码：A 对于常微分方程 问题 』 ()， (1) {u(0)：U。． (2) 其解 U()在 区间 [O，]上存在 、唯一且足够光滑 ，则微分方程 (1)和下面积分方程 M(+)=u(￡)+ff(v+u()dr (3) 等价 (其中步长 T)· 对于初值 问题 (1)和 (2)，目前数值算法有很多，比如欧拉法 、龙格一库塔法等 ，文献 [1—3]研究了数值方 法 ．本文提出一种新的数值方法，并给 出了详细的理论分析和数值试验 ． 1 抛物线逼近法 初值问题 (1)和 (2)，考虑用抛物线与坐标轴所围成的面积作为 (3)式右端定积分的近似值，利用辛 卜生 求积公式2[1得 ／(r+“()dr一-h -[f(t,u(c))+4f(￡十h) 。 ，u(+争))(f+h，M(+^))]+尺(+)， 若令 t．=t，则得 +】) h [f(t j +4f( ) _))_]+RJ．．， (4) 其 中 t j+l ： “( 等[ 十4f( ) tj+l~t／())]． 进一步可以得到 =一 · u㈥()=一 ( ∈( +1)． 为得到积分方程 的解 ，我们可以采用隐式法或者显式法进行计算，详细讨论如下 ． 收稿 日期：2O10—12一l5 基金项 目：国家 自然基金项 目：基于框架理论的遥感 图像压缩算法研究 ：北京市精品课程建设项 目；学校重点 教改项 目(工科数学分析)开放式教学研究与实践 作者简介 杨小远 (1964一)，女，辽 沈阳人 ，教授 ，博士，博导，主要研究方 向为微分方程数值解法． 一 128一 河 南 科 学 第29卷 第2期 1．1 隐 式法 h 若采用新步长 - - I~]t川：+ ，则(4)式可化为 2 u()()+}[／(，，())+4f(t,，“())(+，u())]+R川I，