函授 数值 课件 chap3：数值积分1(牛顿——柯特斯公式).pptVIP

华长生制作 * 第四章 数值积分与数值微分 第四章 数值积分与数值微分 4.1 Newton-Cotes公式 4.2 复合求积法 4.3 Romberg算法 4.4* Gauss求积法 4.5 数值微分 本章要点 (1) 等距节点下的:Newton-Cotes公式和Romberg公式 (2) 数值微分公式 本章应用题: 为了计算瑞士国土的面积,首先对地图作了如下测量:以西 向东方向为x轴,由南向北方向为y轴,选择方便的原点,并将 从最西边界到最东边界在x轴上的区间适当地划分为若干 段,在每个分点的y方向测出南边界点和北边界点的y坐标, 数据如表(单位mm): x 7.0 10.5 13.0 17.5 34.0 40.5 44.5 48.0 56.0 y1 44 45 47 50 50 38 30 30 34 y2 44 59 70 72 93 100 110 110 110 x 61.0 68.5 76.5 80.5 91.0 96.0 101.0 104.0 106.5 y1 36 34 41 45 46 43 37 33 28 y2 117 118 116 118 118 121 124 121 121 x 111.5 118.0 123.5 136.5 142.0 146.0 150.0 157.0 158.0 y1 32 65 55 54 52 50 66 66 68 y2 121 122 116 83 81 82 86 85 68 瑞士地图的外形如图(比例尺18mm:40km) 试由测量数据计算瑞士国土的近似面积,并与其精确值 41288平方公里比较 4.1 Newton-Cotes公式 对于积分 但是在工程技术和科学研究中,常会见到以下现象: 以上这些现象,Newton-Leibniz很难发挥作用 只能建立积分的近似计算方法 这类方法很多,但为方便起见,最常用的一种方法是利用 插值多项式来构造数值求积公式,具体步骤如下: 不同的 插值方法 有不同的 基函数 这就是数值求积公式 为了使一个求积公式能对更多的积分具有较好的实际计 算意义,就要求它对尽可能多的被积函数都准确地成立 因此定义代数精度的概念: 定义1. 若求积公式 则称该求积公式具有m次的代数精度 代数精度也称 代数精确度 例1. 试确定下面积分公式中的参数使其代数精确度尽量高. 解: 因此 所以该积分公式具有3次代数精确度 一、Newton-Cotes数值求积公式 Newton-Cotes公式是指等距节点下使用Lagrange插值 多项式建立的数值求积公式 各节点为 其中 而 因此对于定积分 有 令 即有 n阶Newton-Cotes求积公式 Newton-Cotes公式的余项(误差) 注意是等距节点 所以Newton-Cotes公式化为 思考 使用n次Lagrange插值多项式的Newton-Cotes 公式至少具有n次代数精度,并且n为偶数时至 少具有n+1次代数精度,试以n=1,2,4为例说明 该结果 二、低阶Newton-Cotes公式及其余项 在Newton-Cotes公式中,n=1,2,4时的公式是最常用也 最重要三个公式,称为低阶公式 1.梯形(trapezia)公式及其余项 Cotes系数为 求积公式为 上式称为梯形求积公式,也称两点公式，记为 梯形公式的余项为 第二积分 中值定理 梯形(trapezia)公式具有1次代数精度 故 2.Simpson公式及其余项 Cotes系数为 求积公式为 上式称为Simpson求积公式，也称三点公式或抛物线公式 记为 Simpson公式的余项为 Simpson公式具有3次代数精度 3.Cotes公式及其余项 Cotes系数为 求积公式为 上式称为Cotes求积公式，也称五点公式 记为 Cotes公式的余项为 Cotes公式具有5次代数精度 三、Newton-Cotes公式的稳定性(舍入误差) 考察Cotes系数 因此用Newton-Cotes公式计算积分的舍入误差主要由 其值可以精确给定 记 而理论值为 即 Newton-Cotes公式的舍入误差只是函数值误差的