Ch[].数值分析方法.pptVIP

* 第四章 数值积分与数值微分 本章主要内容： §1 插值型数值求积公式 1.1 一般求积公式及其代数精度 1、提出问题： 一、一般求积公式 2、解决方法： 与f(x)无关的常数,称为积分系数 求积节点 数值积分 介绍一些常用的数值积分公式 数值微分 介绍一些常用的数值微分公式 外推方法 （加速收敛的方法） (1.2)和(1.3)都称之为数值求积公式或机械求积公式 。 为求积公式的截断误差（方法误差）。 3、衡量某种方法好坏的标准： （1）代数精度 （2）数值稳定性 （3）收敛性 对多少次多项式该方法无误差，即计算值是准确的。 或者说成灵敏度如何，也就是看舍入误差对计算结果影响的 即是截断误差的大小。 大小。比如病态方程组，当系数矩阵中的元素有微小变化时，引起 方程组无解。这实际上是由舍入误差或者说成舍入误差的传递造成 的。 二、代数精度 1、代数精度的定义 定义1 若求积公式（1.2）对任意不高于m次的代数多项式都 注：讨论具体问题时，不可能把所有次数小于或等于m的多 等价定义1′：若(1.2)中对于1，x,…,xm都精确成立，对xm+1不精 另外，若代数精度为m，也就是对xm(1.2)或(1.3)精确成立。则 等价定义1′′： 若(1.3)中 R[xi]=0,(i=0,1,…,m)，而R[xm+1]不为0， 精确成立，而对 xm+1不能精确成立，则称该求积公式具有m次代 数精度。 因此有等价定义。 项式列出来验证，因此只要验证对1，x,…,xm 精确成立即可。 确成立，则称(1.2)的代数精度为m。 (1.3)中若f(x)是x的m次多项式，有R(f)=0,因此定义1也可写成： 则称(1.2)的代数精度为m。 2、例子 分析：由等价定义, 解： 所以该求积公式的代数精度m=3。 求代数精度，只对最简单的函数xm来验证。 问题： 1、方法： 1.2 插值型求积公式 插值基函数 与f 无关,记为Ai 插值多项式 插值型求积公式的定义 定义2 对给定互异求积节点 ,若求积系数 2、优点 定理1 数值求积公式(1.2)或(1.3)是插值型的当且仅当它的代数 证明：（必要性 ） 设求积公式（1.2）是插值型的，则 是由(1.4)给出的，则称该求积公式是插值型的。此 时数值求积公式(1.2)称为插值型求积公式。 精度 。 推论1 对给定求积节点 ，代数精度最高的 说明：不研究一般的求积公式。 推论2 若 ，(1.3)是插值型求积公式，则有余项公式 求积公式是插值型求积公式。 3、举例 例2：求插值型求积公式 并确定其代数精度。 分析：实际上该题目是求A0，A1，并确定其代数精度。 从而求积公式为 且n=1,因而代数精度大于等于1，以下验证代数精度从m=2开始 解 (一) 解法（二） 因为是插值型的, 所以代数精度大于或等于 1, 因而对x0 =1, x 该 公式精确成立，即有方程组 以下介绍常用的求积公式 1.3 Newton-Cotes 求积公式 则插值型求积公式称为N-C求积公式。 插值型求积公式 成立。 因此,Newton-Cotes公式为 1、N－C求积系数及公式 2、Cotes系数特点： 表4-1 可查表4-1 首先, （2） (-1)n+k=(-1)n-k 特点： 3、常用的N－C公式及名称 （m=1） －中矩形公式 －左矩形公式 －右矩形公式 几何意义(中矩形公式为例): －梯形公式 （m=1） －梯形公式 （m=1） 注：实际上是用1次Lagrange插值公式近似f(x)计算。 几何意义： 用直边梯形的面积近似。（如图） 以f(x)为曲边的曲边梯形的面积， Newton-Cotes公式称为Simpson公式 （m=3） 几何意义: 为曲边的曲边梯形的面积来近似，因此 用插值抛物线 （英.1710－1761））或抛物线公式。 该公式也称为抛物线公式。(如图) 以f(x)为曲边的曲边梯形的面积， 当n=3,4时，Newton—Cotes公式分别称为Simpson 法则 （m=5） （m=3） （公式）和Cotes公式。 *