清华第五版数值分析第4章课件.pptVIP

为什么要数值积分？ 解决办法 1) 我们用不同的办法近似 可得到不同的积分公式。 2）用简单曲线的积分代替复杂曲线的积分。 求积公式举例 1）梯形公式 2）中矩形公式 3）一般公式 求积节点， 求积系数，也称为节点 的权，权仅与节点的选取有关，不依赖于被积函数的具体形式。 这种数值积分方法称为机械求积，特点是积分问题转变为被积函数值的计算。 例 设有求积公式 求A0，A1，A2，使其代数精度尽量高，并问此时求积公式的代数精度 解：（3个未知系数需三个方程） 令求积公式分别对f(x) = 1、x、x2准确成立。即 ? 解之得A0 = A2 = 1/3，A1 = 4/3，即有 又易知求积公式对f(x) = x3也准确成立： 但 所以该求积公式具有3次代数精度。 容易验证: 梯形公式 1 次代数精度 中矩形公式 1 次代数精度 求积公式的构造 1）若已经选定求积节点 则解线性方程组 2）若系数和节点都不确定，则解关于 和 的非线性方程组。 3）用简单函数的积分代替被积函数的积分。 定理 求积公式至少具有 n 次代数精度的充分必要条件是它是插值型的。 证明：1） 2） 求积公式的余项，收敛性，稳定性 若求积公式代数精度为 ，则可设 求出K即可。K不依赖于函数f。令 得到 梯形公式余项 中矩形公式余项 例：求 余项 解：1）确定代数精度是2. 2）设 3）令 求得 收敛性定义 在 中，若 则称求积公式是收敛的。 稳定性定义 设 对任给 若存在 ，只要 就有 则称求积公式是稳定的。 定理：若 ， 则求积公式稳定。 证明 当n=1时， Newton-Cotes公式(2)为梯形公式 当n=２时， Newton-Cotes公式(2)为辛普森（Simpson）公式 证明：令 代人求积公式俩边得到。 3）当 柯特思系数出现负值。 对于n阶的Newton-Cotes求积公式 当n为奇数时，至少具有n次代数精度；当n为偶数时，至少具有n+1次代数精度. 因为被积函数为奇函数。 梯形公式的余项 定理　若f(x)? C2[a,b] ,则梯形公式余项为 Simpson公式的余项 定理　若f(x)? C4[a,b] ,则Simpson公式余项为 证明：1）确定代数精度为3. 2）令 3） 柯特思公式余项 例 用梯形公式和Simpson公式计算积分 ,并估计误差. 解 记a=0, b=1, f(x)=e-x ,则f (x)=-e-x f(x)=e-x , f(x)=-e-x , f(4)(x)=e-x 复合求积法 将区间[a,b]适当分割成若干个子区间，对每个子区间使用低阶求积公式，构成所谓的复合求积公式，这是提高积分精度的一个常用的方法。 复合梯形公式 将区间[a, b]分成n等份，其分点为xi=a+ih (i=0,1,2,…n), 步长h=(b-a)/n. 在每个小区间[xk-1, xk] (k=1,2,…n)上利 用梯形公式 则 称 为复合梯形公式 复合梯形公式的误差 定理：若 ,则 证 ( ) 由于f(x)?C2[a,b] ，利用闭区间上连续函数的介值定理知存在一点??[a,b],使 因为求积系数为正，所以复合梯形公式稳定。 证明 ( ) (