工程计算基础第8章-1.pptVIP

? 复化 Simpson 公式： 4 4 4 4 4 = Sn 注：为方便编程，可采用另一记法：令 n’ = 2n 为偶数， 这时 ，有 ? 收敛速度与误差估计： 定义 　　 若一个复化积分公式的误差满足 且C ? 0，则称该公式是 p 阶收敛的。 复化梯形公式： /*中值定理*/ 2阶收敛 4阶收敛 6阶收敛 例1：计算 解： 其中 = 3.138988494 其中 = 3.141592502 用8等分的梯形公式和4等分的Simpson公式计算 运算量基本相同，都用了9个点 Q: 给定精度 ?，如何取 n ? 例如：要求 ，如何判断 n = ？ 复化求积公式之间的关系 …… 复化梯形公式的细分关系 复化Simpson公式的细分关系 …… 复化Cotes公式的细分关系 …… 余项 练习 8.1 8.3：（1）， （4） 工程计算基础 第9讲 为什么要数值积分？ 在微积分里，按Newton-Leibniz公式求定积分 要求被积函数f(x) ? 有解析表达式； ? f(x)的原函数F(x)为初等函数．  问题 ? f(x)没有解析表达式，只有数表形式 e.g. ? f(x)有表达式，但原函数不是初等函数 e.g.　　　 ，　　　　　　　　 它们的原函数都不是初等函数. 还有被积函数f(x)的原函数能用初等函数表示，但表达式太复杂，例如函数 并不复杂，但积分后其表达式却很复杂，积分后其原函数F(x)为： x 1 2 3 4 5 f(x) 4 4.5 6 8 8.5 积分近似（数值积分）思想 Aj为求积系数，与求积节点xj有关，与f(x)的具体表达形式无关！ E(f)为余项 数值求积问题 构造或确定一个求积公式，要讨论解决的问题有 　 (i) 确定求积系数Ak和求积节点n； (ii)　求积公式的误差估计和收敛性． 　　　用什么标准来判定两个节点数相同的求积公式的“好”与“差”呢？通常用“代数精确度”的高低作为求积公式“好”与“差”的一个标准． 代数精度 　　　定义１　若对任意的r次多项式，求积公式中有E(f)=0，而至少有一个r+1次多项式 ，则求积公式具有r+1次代数精度。 为了保证精度，我们希望求积公式能对“尽可能多”的函数准确地成立，这就提出了所谓代数精度的概念． 　　 定义2 若对函数　　　　　　　　， 求积公式(2)精确成立，即 而　　　　　　 ， 则称其具有n次代数精确度． 　　　 一般地，欲使求积公式 具有m次代数精度，只要令它对于f (x) = 1，x，…，xm 都能准确成立，这就要求 例1. 试确定下面积分公式中的参数使其代数精确度尽量高. 解: 因此 所以该积分公式具有3次代数精确度 给定[a, b]上的点列 (xi, f(xi)), i = 0,1,2,…,n, 构造Ln(x), 则 余项： 插值型求积公式 注 对于一般的函数 f(x)，En( f )≠0, 但若 f(x) 为次数小于 n 的多项式时，因为 f (n+1) (x) = 0, ∴ En( f ) = 0。因此 n 次插值多项式形式的数值积分公式至少有 n 阶代数精 度。 1、Newton-Cotes数值求积公式 Newton-Cotes公式是指等距节点下使用Lagrange插值 多项式建立的数值求积公式 各节点为 其中 而 因此对于定积分 有 令 即有 n阶Newton-Cotes求积公式 Newton-Cotes公式的余项(误差) 注意是等距节点 所以Newton-Cotes公式化为 Nowton-Cotes型求积公式的误差分析 定理 Newton-Cotes求积公式的余项可表示为： 其中 其中 2、低阶Newton-Cotes公式及其余项 在Newton-Cotes公式中,n=1,2,4时的公式是最常用也 最重要三个公式,称为低阶公式 (1). 梯形公式及其余项 Cotes系数为 求积公式为 上式称为梯形求积公式,也称两点公式，记为 梯形公式的余