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

第4章 数值积分与数值微分 Newton-Cotes公式 数值微分 §4.6 高斯型积分 /* Gaussian Quadrature */ 构造具有2n+1次代数精度的求积公式 将节点 x0 … xn 以及系数 A0 … An 都作为待定系数。令 f (x) = 1, x, x2, …, x2n+1 代入可求解，得到的公式具有2n+1 次代数精度。这样的节点称为Gauss 点，公式称为Gauss 型求积公式。 例：求 的 2 点 Gauss 公式。 解：设 ，应有 3 次代数精度。 ? + ? 1 0 1 1 0 0 ) ( ) ( ) ( x f A x f A dx x f x 代入 f (x) = 1, x, x2, x3 不是线性方程组，不易求解。 正交多项式 证明： “?” 对任意次数不大于n 的多项式 Pm(x)， Pm(x) w(x)的次数不大于2n+1，则代入公式应精确成立： 0 = 0 ? x0 … xn 为 Gauss 点, 则公式 至少有2n+1次代数精度。 与任意次数不大于n 的多项式 P(x) （带权）正交。 　 x0 … xn 为 Gauss 点 ? “?” 要证明 x0 … xn 为 Gauss 点，即要证公式对任意次数不大于2n+1 的多项式 Pm(x) 精确成立，即证明： 设 0 ? n次多项式 TH5表明，在[a,b]上带权的n+1次正交多项式的零点就是求积公式的高斯点。 如何通过正交多项式求高斯点？ 正交多项式族{ ?0, ?1, …, ?n, … }有性质： 任意次数不大于n 的多项式 P(x) 必与?n+1 正交。 若取 w(x) 为其中的?n+1则?n+1的根就是 Gauss 点。 再解上例： ? + ? 1 0 1 1 0 0 ) ( ) ( ) ( x f A x f A dx x f x Step 1：构造正交多项式?2 设 c bx x x a x x x + + = + = = 2 2 1 0 ) ( , ) ( , 1 ) ( j j j ? ? 5 3 - = a 0 ) ( 1 0 = + ? dx a x x 0 ) , ( 1 0 = j j ? ? = + + - ? = = + + ? = 1 0 2 1 1 0 2 1 0 0 ) )( 5 3 ( 0 ) , ( 0 ) ( 0 ) , ( dx c bx x x x dx c bx x x j j j j 21 5 9 10 = - = c b 即： Step 2：求?2 = 0 的 2 个根，即为 Gauss 点 x0 ，x1 Step 3：代入 f (x) = 1, x 以求解 A0 ，A1 解线性方程组，简单。 结果与前一方法相同： ? 利用此公式计算 的值 注：构造正交多项式也可以利用 L-S 拟合中介绍过的递推式进行。 Gauss型求积公式具有数值结果精度高,收敛得以保证、计算简便、易于在计算机上实现等优点,并且在积分区间[a,b]有限时便于推广到高维数值积分.不足之处是公式的构造比较困难,另一个是由于相邻次数的正交多项式的根,从而造成增加求积节点以提高计算结果的精度时,原先所有求积节点上的函数值全部无用.所以在具体应用Gauss求积公式计算数值积分时,n取得都较小. ? Gauss 公式的余项： /* 设P为f 的过x0 … xn的插值多项式 */ /*只要P 的阶数不大于2n+1，则下一步等式成立*/ 插值多项式的余项 Q：什么样的插值多项式在 x0 … xn 上有 2n+1 阶？ A：Hermite 多项式！ 满足 Th6 高斯求积公式的求积系数Ak(k=0,1,2…,n)全是正的。 是n次多项式。 2n次多项式 因而，高斯公式能够准确成立。 推论：由定理6得知，高斯求积公式是稳定的。 ? 特殊正交多项式族： ① Legendre 多项式族： 1 ) ( ? x r 定义在[?1, 1]上， 满足： 由 有递推 以 Pn+1 的根为节点的求积公式称为Gauss-Legendre 公式。 ② Chebyshev 多项式族： 2 1 1 ) ( x x - = r 定义在[?1, 1]上， Tn+1 的根为 k = 0, …, n 以此为节点构造公式 称为 Gauss-Chebyshev 公式。