第五章 数值积分及数值微分2014.ppt

高斯－勒让德公式 此为高斯－勒让德公式，区间为[-1,1],勒让德正交多项式Pn+1(x)的零点就是其高斯点。 例：取 其两个零点为 。求积公式为 令它对f(x)=1,x成立，有 高斯－勒让德公式(续） 验证：分别令f(x)=x2,x3,x4则 两点高斯－勒让德公式 高斯－勒让德公式(续） 区间[-1,1]上权函数W(x)=1的Gauss型求积公式,称为Gauss-Legendre求积公式,其Gauss点为Legendre多项式的零点. Gauss-Legendre求积公式 公式的Gauss点和求积系数可在数学用表中查到 . 由 因此,[a,b]上权函数W(x)=1的Gauss型求积公式为 高斯－勒让德公式(续） n xk Ak n xk Ak 1 0 2 6 ±0.9324695142 ±0.6612093865 ±0.2386191861 0.1713244924 0.3607615730 0.4679139346 2 ±0.5773502692 1 3 ±0.7745966692 0 0.5555555556 0.8888888889 7 ±0.9491079123 ±0.7415311856 ±0.4058451514 0 0.1294849662 0.2797053915 0.3818300505 0.4179591837 4 ±0.8611363116 ±0.3399810436 0.3478548451 0.6521451549 8 ±0.9602898565 ±0.7966664774 ±0.5255324099 ±0.1834346425 0.1012285363 0.2223810345 0.3137066459 0.3626837834 5 ±0.9061798459 ±0.5384693101 0 0.2369268851 0.4786286705 0.5688888889 例题 运用高斯――勒让德公式计算积分 解：两点公式 两点梯形公式 例题（续） 三点公式: 三点辛普森公式: 例题（续） 例 用二点高斯-勒让德公式计算积分 解 作变量代换 则 记 ,因为节点 得 所以,由二点高斯公式 带权的高斯公式 对于任意次数不超过2n+1的多项式均能准确成立 称其为带权的高斯公式。其中ρ(x)为权函数 称为高斯－切比雪夫公式。高斯点为n+1次切比雪夫多项式的零点 当 时，所建立的高斯公式 高斯公式例题 例：构造下列形式的高斯公式 解：令它对于f(x)=1,x,x2,x3准确成立，得 由于x0A0+x1A1=x0(A0+A1)+(x1-x0)A1 利用第1式，可将2式化为 同样，利用2式化3式，利用3式化4式，分别得 高斯公式例题(续) 高斯公式是 5.6 数值微分 自然而又简单的方法就是，取极限的近似值，即差商 1.当函数f(x)以离散点列给出时，要求我们给出导数值， 2.函数f(x)过于复杂 这两种情况都要求我们用数值的方法求函数的导数值 微积分中，关于导数的定义如下： 5.6 数值微分 5.6 数值微分 由Taylor展开 因此，有误差 5.6 数值微分 5.6 数值微分 由Taylor展开 因此，有误差 5.6 数值微分 5.6 数值微分 由Taylor展开 因此，有误差 5.6 数值微分 由误差表达式，h越小，误差越小，但同时舍入误差增大， 所以，有个最佳步长。我们可以用事后误差估计的方法来确定。设D(h)D(h/2)分别为步长为h，h/2的差商公式。则 时的步长h/2 就是合适的步长 5.6 数值微分 h f’(1.15) R(x) h f’(1.15) R(x) 0.10 3.1630 -0.0048 0.05 3.1590 -0.0008 0.09 3.1622 -0.0040 0.04 3.1588 -0.0006 0.08 3.1613 -0.0031 0.03 3.1583 -0.0001 0.07 3.1607 -0.0025 0.02 3.1575 -0.0007 0.06 3.1600 -0.0018 0.01 3.1550 -0.0032 5.6 数值微分 插值是建立逼近函数的手段，用以研究原函数的性质。因此，可以用插值函数的导数近似为原函数的导数插值型数值微分 插值型数值微分 误差 5.6 数值微分 给定点列 且, 求 解： 5.6 数值微分 称为三点公式Taylor展开分析，可以知道，它们都是 ,称为三点公式 5.6 数值微分 用三次样条函数求数值微分设s(x)是f(x)在各区间［xk,xk+1］上的三次样条插值函数,则在区间［xk, xk+1］上可通过三次