b12积分(下).ppt

第七章 数值积分与微 分（下） §6 高斯型求积公式 一般理论(续1） 一般理论(续2） 一般理论（续3） 证明插值型求积公式的代数精度最高不超过2n－1次 一般理论举例 一般理论举例（续1） 一般理论举例（续2） 一般理论举例（续3） 带权函数的Gauss型求积公式 举例说明Gauss点与正交多项式的关系 例11（续1） 例11（续2） 例11（续3） 求解例11方法小结 求Gauss点的一般方法 Gauss点的充要条件 定理7.4（续1） 定理7.4（续2） 定理7.4证明 定理7.4证明（充分性） 定理7.4证明（充分性）（续） 6.2 Gauss型求积公式的误差 求积公式系数Ak 0的证明： 6.3 几种常用的Gauss型求积公式 Gauss—Legendre（勒让德）求积公式（续） n个Gauss点xk可取作pn(x)的零点 求相应的系数Ak 求相应的系数Ak （续） 求相应的系数Ak（续1） 求具体的零点(Gauss点)及对应的系数 N=3的零点(Gauss点)及对应的系数 从[-1,1]推广到[a,b] 表7-2 Gauss-Legendre求积公式计算举例 例12（续） Gauss-Legendre求积公式计算举例（续） 2. 高斯-切比雪夫（Gauss-Chebyshev）求积公式 3. 高斯-拉盖尔（Gauss-Laguerre）求积公式 高斯-拉盖尔（Gauss-Laguerre）求积公式(续) 4. 高斯-埃尔米特（Gauss-Hermite）求积公式 高斯-埃尔米特（Gauss-Hermite）求积公式（续） Gauss-Hermite求积公式计算积分举例 用Gauss型求积公式计算积分近似值的步骤 解非线性方程组例子 正交多项式求解例子 §7 数值微分 中点法 中点法的Taylor公式推导 中点法公式计算导数f ? (a) 差商代替微商举例 7.2 插值型的微分公式 插值型的微分公式（续1） 插值型的微分公式（续2） 等距节点在节点xi处的导数值的求法 2. 三点公式 ： 误 差 插值型的微分公式举例 例14 (续） 7.3 外推法 外推法（续1） 外推法（续2） 外推法（续3） 外推法（续4） 外推法举例 外推法举例（续） 第七章 结 束 下面求具体的零点(Gauss点)及对应的系数 ,由它们即可构成 求积分的Gauus型公式: 三次代数精度 截断误差随n的增大， 越来越小? 0，非常 精确，? 更高阶的 （n更大的） Gauus型求积公式不用。 代数精度为最高5次 (三个点最低为2次)。 上面是在[－1 , 1]上，对任意区间[ a , b ]可作变量代换： 点tk为n次Legender多项式pn ( t )的零点，Ak同前一样。 给出了部分Gauss-Legendre求积公式的节点与 求积系数值，以便查用。 n xk Ak n xk Ak 1 0 2 5 ±0.9061798 ±0.5384693 ±0.2369268 ±0.4786286 ±0.5688888 2 ±0.5773502 1 3 ±0.5773502 ±0±06 ±0.9324695 ±0.6612093 ±0.2386191 ±0.1713244 ±0.3607615 ±0.4679139 4 ±0±0±0例12 若用n=2的Gauss-Legendre求积公式计算，由表7-2， x1=-0 x2=0A1,2=1，则 ： （紧接下屏） 若用n=3的Gauss-Legendre求积公式，查表7-2得： 其误差为： 如果用三点的Simpson公式求积分近似值，则有 I 的精确值为1。由例12结果可见，当被积函数充 分光滑，在节点数相同时，Gauss型求积公式比 N-C公式的精度高。 Chebyshev多项式是[－1, 1] 上关于权函数: 的正交多 项式系， 其表达式为： n次Chebysher多项式Tn(x)的零点为 以xk为节点，求积系数为： 的带有 权函数 的Gauss型 求积公式: 称为Gauss-Chebyshev求积公式。其余项为： Laguerre多项式: 是区间[0,+?]上关于权函数 ? (x)= e - x的正交多项式系，以 n次Laguerre多项式的零点{xk}为Gauss点导出的Gauss型 求积公式: 称为Gauss-Laguerre求积公式。求积系数的计算公式为: 求积公式（7-38） 的截