二阶偏微分方程组 龙格库塔法_数值积分公式的推广:龙格库塔法(一)

1. 写在前面的话

龙格库塔法是求解常微分方程的数值计算方法。它始于1895年,历史悠久,后人以龙格(Runge)和库塔(Kutta)两位计算数学家的名字为它命名。关于它的研究文献有很多,至今仍有新的论文不断涌现,卷帙之浩繁,恐怕难以完全统计。

先后在此领域内作出重要贡献的科学家除了龙格和库塔两位外,还有霍伊恩(Heun),奈斯特龙(Nyström),胡塔(Huta),吉尔(Gill),柯蒂斯(Curtis),布切尔(Butcher),默森(Merson),库珀(Cooper),海尔(Hairer),万纳(Wanner),桑兹塞尔纳(Sanz-Serna)等。

其中新西兰数学家布切尔是该问题的当代大师. 他创造性地将抽象代数和图论的方法引入到龙格库塔法的研究当中,奠定了龙格库塔法根树表示论的理论基础,进而改变了整个数值方法的面貌。以他的名字命名的数学名词有:布切尔群,布切尔表,布切尔级数(简称为B级数)。关于龙格库塔法,除了一系列重要的研究性论文,他还写过好几篇值得一读的综述类文章,并有专著传世,感兴趣的读者可以查阅其著作或文献(比如【1,2】)。

本文仅就龙格库塔法背后隐藏的机理作一番初步探讨,有兴趣的读者可以进一步思考或加以研究。

2. 数值积分公式

在介绍龙格库塔法之前,我们先来回顾一下数值积分公式。随便翻开一本计算方法的教科书,大家都能看到下面的数值求积公式:5a2bd5da97da47e5f264ac5399636095.png其中 941dc6063a1da534c53d0dcdbe2f4021.png 称为权系数,而 72a9361142d94d1308b2830b53cf9482.png 称为求积节点。这里假定被积函数比较光滑。通常情况下,节点在积分区间 (0, 1)内部或边界上选取。

其中最常用的求积公式是所谓的插值型求积公式,即用函数 f 的拉格朗日插值多项式来逼近被积函数, 从而得到a2eca231362995a612c0953596272a61.png再将积分号与求和符号交换, 得到5a2bd5da97da47e5f264ac5399636095.png其中873ea6111623c9e56afe759c9a6b2ca7.png注意, 当节点个数只有一个时 (即 s = 1),可定义 6af6d4768a919ee7eb771aa5638dd4c5.png

求积节点的选取可以多种多样,不同的选取方式将导致不同的逼近误差,最终影响到数值求积的精度。比如,若分别选取高斯(Gauss)节点、拉道(Radau)节点、洛巴托(Lobatto)节点,则相应得到高斯型、拉道型及洛巴托型数值求积公式,其代数精确度分别为 fdc9c4e5514a1a2ad4e23daa397211dc.png。所谓数值积分公式的代数精确度,是指当被积函数 为多项式函数时,使得数值求积公式精确成立的多项式的最高次

  • 1
    点赞
  • 4
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值