第六章 数值积分和数值微分.ppt

Newton-Cotes 积分 若积分节点可自由选取，则一个简单的办法就是取等距节点，即对区间做等距剖分。 这种数值积分称为Newton-Cotes积分 * a=x0 x1 x2 xn-2 xn-1 x4 xn=b x3 设节点步长 (b-a) 与积分点坐标、积分区间长度无关，故可以预先求出。 Newton-Cotes系数 * 积分系数 n＝1时 梯形公式 * n＝2 时 Simpson公式 * * 1、梯形公式 此处用了积分中值定理 误差 线性（1次）插值余项 * 不难看出，梯形公式具有1阶代数精度 2、Simpson 公式 注意到 Simpson 公式有3 阶代数精度，为了推导该数值积分公式的误差，我们构造一个3 次多项式 为0,why? * 则通常有 因此，在 N-C 积分中，n为偶数时有n+1 阶代数精度，而n为奇数时有n 阶代数精度。 例如：n=1时， 梯形公式，1阶代数精度； n=2时，Simpson公式，2+1=3阶代数精度。 令 * 6.2 复化数值积分 数值积分公式与多项式插值关系非常密切！因此Runge现象的存在，使得我们不能用太多的积分点计算。与插值时的情况类似，可采用分段、低阶的方法。 * 取等距节点， 复化梯形公式 * 积分误差 a=x0 x1 x2 xn-2 xn-1 x4 xn=b x3 由均值定理, 可以看出，复化梯形公式是收敛的。 思考：如果节点不等距，还可以做复化积分吗？ * 积分误差 取等距节点， 复化Simpson公式 * 积分误差 由均值定理知， 可以看出，复化Simpson公式是收敛的。 * a=x0 x1 x2 xn-2 xn-1 x4 xn=b x3 问题：一共做了多少次Simpson？ n/2=m次 积分误差 a=x0 x1 x2 x3 x4 f(x) x5 b=x6 复化梯形与复化Simpson 设节点等距， 定义 　　 若一个积分公式的误差满足 且C ? 0，则称该公式是 p 阶收敛的。 ~ ~ ~ 例：计算 解： 其中 = 3.138988494 其中 = 3.141592502 运算量基本相同 * Lab06 复化积分 1.分别编写用复化Simpson积分公式和复化梯形积  分公式计算积分的通用程序。 2. 用如上程序计算积分 取等距节点，记节点{xi , i=0,…N}，其中N为 {2k,k=0,1,…,10} ，并计算误差(用科学计数形式) ，同时给出误差阶（用浮点形式，比如1.8789）。 3. 简要分析两种方法的优劣。 * 误差阶： 则，相应的误差阶为： * 记步长为h时的误差为 ,步长为h/n时的误差为 （这里n=2）, Sample Output (? represents a space) 复化梯形积分，误差(科学计数形式)和误差阶为 k=0?,?e0=0. ############e00 k=1?,?e1= 0.############e-1, d1=？ 比如 d1= 1.1111 k=2?,?e2= 0.############, d2=？ ... 复化Simpson积分，误差和误差阶为 k=1?,? e0= 0. ############e00 k=2?,? e1 = 0.############e-1, d1=？ ... * 数 学 系 University of Science and Technology of China DEPARTMENT OF MATHEMATICS 数 学 系 University of Science and Technology of China DEPARTMENT OF MATHEMATICS 第六章 数值积分和数值微分 * 6.1 数值积分 回顾：在微积分中，Newton-Leibniz 公式 但在实际中，该公式往往没法用；为什么？ 数值积分：即是有限离散节点上的函数值的某种线性组合 * 1、只知道目标函数 的一些离散数据值； 2、原函数F(x)不存在； 3、原函数F(x)存在但非常复杂，不易求出或求不出； 其中，称 为积分系数；与f(x)无关，与积分区间和积分点有关。 例： * 为数值积分， 为精确积分， 则称数值积分 有 k 阶代数精度是指： 问题：如果判断好坏？ 代数精度 对任意次数不高于k次的多项式 f(x)， 数值积分 都没有误差(但k+1次的有误差)。 * 记 用插值函数的积分，作为数值积分 代数精度 由Lagra