[生活]第一讲定积分的数值计算.docx

［生活］第一讲 定积分的数值计算 第一讲定积分的数值计算 【主要目的】围绕定积分的概念与数值计算方法这一大家非常熟悉的主题，突 出数值实验、几何观察、数值分析等实验特性，学生通过实验与理论的对照，加深 对数学思想和数学知识的理解和掌握，学习如何从实验角度创新知识、发现知识, 并上升到理论分析的角度。 【主要内容】定积分的数值计算方法，包括：矩形法、梯形法与辛普森法；对误 差的了解:精度与收敛速度 引言 /（町巾丄］［為轲 首先回忆一下函数在区间上的定积分概念的建立过程。考虑在区间 内任意插 入个分点的分法： 入个分点的分法： 討一1丁总=呵 5 S-i 把分割成个小区间，第个子区间的长度为 臥W【S丙］/匾）g （出二12…「M） 任取数，做乘积，把所有这些乘积相加得到和 日⑺=聽伽吒［心 如果无论区间怎样划分及分点怎样选取，当时，该和式 都趋于同一常数，则 称函数在区间上可积，且称此常数为在区间 上的定积分，即 川町仏切八）［算慕”磁=龈。言心g 0称和式为积分和或黎曼和。 在定积分的概念中包含了两个任意性，即对区间的分割和点的选取都是任意 的。 显然，对于区间的不同分割或者点的选取不同，得到的和式一般不同。定积分 的定 ⑺一＞0 义中要求在对区间无限细分（）的条件下，所有这些和式都趋于同一数值。这 一点初学者较难理解。我们将通过数值实验来加以理解。 /（町和上］/^（刃切 当在区间上连续，为在区间上的原函数时，我们可以 用牛顿-莱布尼兹公式 方便地求得。但是有些函数其原函数不能用初等函数表示出来，这样对应的定 积分通常也不能用牛顿-莱布尼兹公式算出其精确值。而且，在自然科学与工程技 术中有许多问题，被积函数并不是用具体函数表达式解析表示的，而经常是通过实 验或测量方法用表格或图形给出的，这就导出了定积分的数值计算问题。我们将利 用“分割取近似，作和求极限”这一定积分思想方法，来构造一些数值计算方法, 并进行数值实验。 实验一定积分概念的深化 达布和 设函数在区间上有界。考虑将将区间任意分割成个子区间 Jt-J r （）的分法，设在子区间上的上、下确界分别 上莎Jt二加尤-叫［耳?心1川力 为，称为在子区间上的振幅，和式 F a ⑺二評％ KO = 丫如/W/⑴T 分别称为关于该分割的达布（Darboux）大和与达布小和。由定义可知，函数 对 应于同一分割的积分和有无穷多个，但达布大和与达布小和却都各只有一个。当 在区间上不连续时，达布和不一定是积分和，但它们都与积分和有着密切的联系, 易知道对于同一分割，有 N Jt-L .可以证明在区间上可积的充分必要条件是 宀“ lim ［S（r）-s（r）］= lim Y 恋, = 0 /（町忆切伸丸 垃屮n 现在假定在区间上非负连续，那么达布大和在几何上就表示在子区 珀円1见目⑺心 间上以为高所做的个小矩形构成的阶梯形的面积；达布小和表示在 子区间上以 为高所做的个小矩形构成的阶梯形的面积，它们的差 ［珀円1 g ⑺■切二詁加八）［以］奶T 0 就是这两个阶梯形面积之差。由于函数在区间上可积，所以当，即 当区间被 无限细分时，这两个阶梯形面积都趋于该曲边梯形的面积，从而这两个阶梯 形面积之差为零，即 当考虑对区间进行等分时，我们有 相应地将、分别记作和. 特别，如果在区间上单调增加，那么达布小和就是左和 H亠1 A -— 0、 A — e 如=L如=少+上——)—— ，达布大和就是右和 左 3)二 Eg竝 W)= 斗上Gk i a /(刃巾上］ 如果在区间上单调减少，那么达布大和就是左和，达布小和就是右和，即 ￡㈤=血⑻.附)=殛俎⑷㈣=汽盟］I［丽阴 数值实验1对区间上作等分，观察在上的达布大和与达布小 和之差随增加时 的变化趋势。 ￡3)-(对)閔 Mathematica 程序(ch1-ex1.nb) 实验过程: 改变分割次数，观察；(2)改变被积函数观察 实验结果分析与理解： 5璋寸一卫辰方总谕描 从实验看出，对于函数，它在上的达布大和与达布小和之差 随增加而趋于0.达布大和与达布小和分别趋于曲边三角形的面积。 实验二定积分数值计算方法一一近似计算 如果在区间上可积，那么我们已经知道用它的左和或右和来逼近它，我们 称之为矩形求积公式。当越大，逼近的精度越高。根据上述求积公式 I矽』y仗）必￡馬筋￡（呵（/ （打 当为增函数时，， 用g触（科）￡ E』（力必疙 砂炒［/ （町 当为减函数时，， 我们甚至知道什么时候左和及右和给出的是过剩的近似值还是不足的近似值。 从上面的数值实验例子可以看到，当=2时，左和给出了一个相当差的不足近似 值, 而右和也只给出了一个相当差的过剩近似值。当然，当充分大时，它们都能给 出好的近似值。但是，在给定的条件下，我们如何修改计算求积公式，使本例