数值积分数值微分21599.ppt

设 表示二分次后求得的梯形值： 计算步骤： k表示二分次数，m表示加速次数 表示表示序列 的m次加速值。 则： 例如:地球卫星轨道是一个椭圆，椭圆周长的计算公式是 这里a是椭圆的半径轴，c是地球中心与轨道中心（椭圆中心） 的距离，记h为近地点距离，H为远地点距离，R=6371（km）为 地球半径，则 我国第一颗地球卫星近地点距离h=439(km)，远地点距离 H=2384(km）。试求卫星轨道的周长误差不超过10^-5。 1.564646 1.564646 1.564646 2 1.564648 1.564646 1 1.564640 0 即人造卫星轨道的周长为48708km 从而有 解： §5 高斯求积公式 定义：若一组节点 x0 … xn ∈[a,b],是使插值型求积公式 具有2n+1次代数精度。这样的节点称为Gauss 点，Ak称为Gauss系数，求积公式称为Gauss 型求积公式。 一、高斯求积的一般理论 节点 x0 … xn 以及系数 A0 … An 都作为待定系数。 要使求积公式具有2n+1 次代数精度，令 f (x) = 1, x, x2, …, x2n+1 代入求积公式精确成立，解出xk和 Ak . 例：求 的 2 点 Gauss 公式。 解：设 ，应有 3 次代数精度。 ? + ? 1 0 1 1 0 0 ) ( ) ( ) ( x f A x f A dx x f x 代入 f (x) = 1, x, x2, x3 不是线性方程组，不易求解。 从求积过程知需求解非线性方程组，可以利用正交多项式的特性来构造求积公式。 *数计学院《数值计算》课程建设组QAB 第五章 数值积分与数值微分 5.1 引言 5.2 Newton-Cotes 公式 5.3 复化求积公式 5.4 龙贝格求积公式 5.5 高斯型积分 5.6 数值微分 【本章重点】1. 求积公式代数精确度定义，应用此定义建立求积公式。2. 梯形公式，Simpson公式及它们的复合公式及余项表达式和误差估计。3. 外推原理及Romberg求积公式。4. Gauss型求积公式及求积节点（即Gauss点）的充要条件，Gauss-Legendre求积公式与Gauss-chebyshev求积公式。5. 求积公式收敛性与稳定性概念和基本结论。 【课前思考】1．什么是代数精确度？如何确定和验证求积公式的代数精确 度和验证求积公式的代数精确度次数？如何利用代数精确 度概念确定求积公式系数与节点？2．写出梯形公式、Simpson公式及复合梯形公式及复合Simpson 公式和它们的截断误差。复合梯形公式和复合Simpson公式的 误差是步长h的几阶小量？若给定误差 ，如何根据复合求积 公式的误差估计求积区间等分数3 ．什么是Gauss型求积公式？什么是Gauss点？如何求Gauss 点？ §1 引言 一、数值求积的基本思想 1、牛顿-莱布尼兹公式 但是求函数f(x)的原函数F(x)不一定比计算积分容易， 例如函数 找不到用初等函数表示的原函数。 另外若给出的函数f(x)是数据表，也不好求函数的积分。 计算定积分的方法： 2、积分中值定理 但是点 的具体位置一般不知道，故难以准确算出 的值。 有两种近似方法： 得到 梯形公式 另一种是矩形法， 左矩形公式 一种是梯形法，用 代替 中矩形公式 右矩形公式 简称矩形公式 定义：公式 叫做数值求积公式（机械求积） 其中xk称为求积节点， Ak称为求积系数,亦称伴随节点xk的权。 二、 代数精度的概念 定义 　　　若某个求积公式对于次数不超过m的多项式均能准确成立，但对于m+1次的多项式就不能准确成立，则称此求积公式的代数精度为 m 。 一般要使求积公式 具有m次 精度，只要令它使 都能准确成立，即： 解：逐次检查公式是否精确成立 当 f(x)=1： = 当 f(x)=x ： = 当 f(x)= x2 ： ? 故：代数精度 = 1 例1 考察 有几次代数精度。 三、 插值型的求积公式 设给定一组节点: 且已知函数在这些节点处的函数值f(xi)(i=0,1,…n),由第二章知可以作插值函数Ln(x),由于Ln(x)为多项式，所以其积分可以很容易求得： 记 则上式= 称公式： 其中