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

第三章 数值积分与数值微分（二）;一、梯形公式的递推公式及事后估计法; 注意到每个子区间[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)右端积分用辛普森求积公式，则有;四、三次样条求导;五、数值微分的外推算法;(6.8