数值分析07-数值积分及微分(下).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） 外推法举例 外推法举例（续） 第七章 结 束 仿上可知：若以T(h)近似f ?(x)，误差阶为o(h2)将h缩小一半： 将此式乘 4减去式 （7-58） 式得： 回到这里： 所以Gauus - Legender求积公式的系数为： Ak还有其它形式。 并且当f ( x ) ?1代入可得： 而余项： 因为： 所以有： 下面求具体的零点(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