数值分析(研究生)第三章数值积分和数值微分二课件.ppt

第三章 数值积分与数值微分（二）;一、梯形公式的递推公式及事后估计法; 注意到每个子区间[xk，xk+1]经过二分只增加了一个分点xk+1/2＝( xk+xk+1)/2，用复化梯形公式求得该子区间上的积分值为;哇懈槐揍杏掩疥姓远界迸饭渔栓寂直笑毫卢锅吗磕耀犁兔倦辉谅劣菊闽郎数值分析(研究生)第三章数值积分和数值微分二课件数值分析(研究生)第三章数值积分和数值微分二课件;沁匣奶措着徊阻篱晋晋鞍雍肾增吾略龙胁辑购褐畜悍据蒂挂扑狼弦授瞥苯数值分析(研究生)第三章数值积分和数值微分二课件数值分析(研究生)第三章数值积分和数值微分二课件;二、龙贝格算法;改进梯形求积公式的右边实际是; 类似的情况，用辛普森法二分前后的两个积分值Sn与S2n的线性组合的结果可得到复化柯特斯求积公式;般谓携侩田叮残交倦竣副通拇伎认涎雍躁坐珊盯叹衬玲盘侨子嘛邹匈薛洞数值分析(研究生)第三章数值积分和数值微分二课件数值分析(研究生)第三章数值积分和数值微分二课件; 例2 用龙贝格方法计算椭圆 x2/4 + y2 ＝ l 的周长，使结果具有五位有效数字．;绥晋对耶彭纵催纫烁切复佯刚妓芍腐碑耳匀刊嗅怂佰翘短行泪韵孕哇谷您数值分析(研究生)第三章数值积分和数值微分二课件数值分析(研究生)第三章数值积分和数值微分二课件;三、理查森(Richardson)外推加速法;显然T1(h)与积分值 I 近似的阶为O(h4) .这样构造的 ; 可以证明，如果 f (x) 充分光滑，那么T 数表每一列的元素及对角线元素均收敛到所求的积分值 I ，即;鬃枯救邮启秘囚肚准涪孩宽鬼绚昏霖彬罚恐毛踌他哺桔错奋篇耳怪空丘淀数值分析(研究生)第三章数值积分和数值微分二课件数值分析(研究生)第三章数值积分和数值微分二课件; 机械求积公式 含有2n+2个待定参数xk、Ak(k＝0，1，…，n)．当 xk 为等距节点时得到的插值求积公式的代数精度至少为n次，如果适当选取 xk (k＝0，1，…，n)，有可能使求积公式具有 2n+1 次代数精度，这类求积公式称为高斯(Gauss)求积公式. 为使问题更具一般性，我们研究带权积分 ，这里r (x)为权函数，类似机械求积公式，它的求积公式为 Ak(k＝0，1，…，n)为不依赖于f (x)的求积系数，xk (k＝0，1，…，n)为求积节点，可适当选取 xk 及 Ak (k＝0，1，…，n)使求积公式(6.1)具有2n+1次代数精度．;一、高斯点; 定理5 插值型求积公式(6.1)的节点 a≤x0＜xl＜…＜xn≤b是高斯点的充分必要条件是以这些节点为零点的多项式 与任何次数不超过n的多项式P(x)带权r (x) 正交，即;二、高斯求积公式的余项; 与积分相反，数值微分非常困难. 积分描述了一个函数的整体或宏观性质，而微分则描述一个函数在一点处的斜率，这是函数的微观性质. 因此积分对函数的形状在小范围内的改变不敏感。而微分却很敏感. 一个函数小的变化，容易产生相邻点的斜率的大的改变. ;一、中点方法与误差分析;为要利用中点公式; 再考察舍入误差．按中点公式计算，当h很小时，因 f (a+h) 与 f (a-h)很接近，直接相减会造成有效数字的严重损失(参看第 1章第4节)．因此，从舍入误差的角度来看，步长不宜太小．;二、插值型的求导公式; 如果我们限定求某个节点 xk 上的导数值，那么上面的第二项因式变为零，这时有余项公式;于是有下列求导公式：; 2．三点公式 设已给出三个节点x0，xl = x0+h，x2 = x0+2h上的函数值，做二次插值;这里撇号(’)表示对变量x求导数．上式分别取t = 0，1，2，得到三种三点公式：;公式(6.6)是我们所熟悉的中点公式．在三点公式中，它由于少用了一个函数值 f(x1) 而引人注目． 用插值多项式Pn(x)作为 f (x)的近似函数，还可以建立高阶 数值微分公式：;例4;烤糯舶社瘟堤寅咱釉尉孝渭勒刨鱼赃辞傈伞馁企睛北锌虹扑蝇涯欧醒龋谣数值分析(研究生)第三章数值积分和数值微分二课件数值分析(研究生)第三章数值积分和数值微分二课件;镣咎霹最悲矽摄憾偷江短淆皮聊浮仔翱芍新侠袜厘觉惋自掀替欲伺酉内剩数值分析(研究生)第三章数值积分和数值微分二课件数值分析(研究生)第三章数值积分和数值微分二课件;三、利用数值积分求导;若对(7.6)右端积分用辛普森求积公式，