数值分析课件-(第4章).ppt

《数值分析课件-(第4章).ppt》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、第四章数值积分和数值微分内容提要4.1引言4.2牛顿-柯特斯公式4.3复化求积公式4.4龙贝格求积公式4.5高斯求积公式4.6数值微分4.1引言一、数值求积的基本思想对定义在区间［a,b］上的定积分但有时原函数不能用初等函数表示，有时原函数又十分复杂，难于求出或计算；另外如被积函数是由测量或数值计算给出的一张数据表示时，上述方法也不能直接运用。因此有必要研究积分的数值计算问题。积分中值定理告诉我们：平均高度f(ζ)aζbyxy=f(x)0af((a+b)/2)byxy=f(x)0abyxy=f(x)0梯形公式平均高度中矩形公式平均高度更一般地，我们构造具有下列形式的求积公式求积节点求积系数这

2、类数值方法通常称为机械求积，其特点是将积分求值问题归结为函数值的计算，这就避开了牛顿-莱布尼兹公式需要寻求原函数的困难。二、代数精度的概念利用代数精度的概念构造求积公式三、插值型的求积公式4.2牛顿-柯特斯公式一、牛顿-柯特斯公式的导出柯特斯系数牛顿-柯特斯公式的代数精度4.3复合求积公式一、问题与基本思想在使用牛顿-柯特斯公式时将导致求积系数出现负数(当n≥8时,牛顿.柯特斯求积系数会出现负数)，因而不可能通过提高阶的方法来提高求积精度。为了提高精度通常采用将积分区间划分成若干个小区间，在各小区间上采用低次的求积公式(梯形公式或辛普森公式)，然后再利用积分的可加性，把各区间上的积分加起来，

3、便得到新的求积公式，这就是复化求积公式的基本思想。本节只讨论复化的梯形公式和复化的辛普森公式。二、复合梯形公式三、复合辛普森公式xi01/81/43/81/2f(xi)10.99739780.98961580.97672670.9588510xi5/83/47/81f(xi)0.93615560.90885160.87719250.84147094.4龙贝格求积公式一、梯形法的递推化(变步长求积法)于是可以逐次对分形成一个序列{T1,T2,T4,T8,…},此序列收敛于积分真值I。当

4、T2n-Tn

5、<ε时，取T2n为I的近似值。以上算法称为变步长求积法。但由于此序列收敛太慢。下节我们将其改造

6、成为收敛快的序列。二、龙贝格算法如何提高收敛速度以节省计算量是龙贝格算法要讨论的中心问题。这样我们从收敛较慢的{Tn}序列推出了收敛较快的{Sn}序列。可以证明{Sn}序列实际上就是逐次分半的复化辛普森公式序列。这样我们从{Cn}序列又推出了收敛更快的{Rn}序列.{Rn}序列也称为龙贝格序列。我们从收敛较慢的{Tn}序列只用了一些四则运算，便推出了收敛更快的{Sn}序列,{Cn}序列和{Rn}序列。T1T2S1T4S2C1T8S4C2R1T16S8C4R2﹕﹕﹕﹕运算顺序表kT2kS2k-1C2k-2R2k-300.920735510.93979330.946145920.94451350

7、.94608690.946083030.94569090.94608330.94608310.9460831这里利用二分3次的数据（它们的精度都很差，只有两三位有效数字）通过三次加速求得R1=0.9460831,这个结果的每一位数字都是有效数字，可见加速效果是十分显著的。4.5高斯求积公式一、一般理论4.6数值微分一、中点方法与误差分析数值微分就是要用函数值的线性组合近似函数在某点的导数值。由导数定义差商近似导数得到数值微分公式。hG(h)hG(h)hG(h)10.36600.050.35300.0010.35000.50.35640.010.35000.00050.30000.10.353

8、50.0050.35000.00010.3000二、插值型的求导公式知识结构图四数值积分与数值微分数值积分基本概念牛顿-柯特斯公式复合求积公式数值微分中点方法插值型求导公式龙贝格求积公式高斯求积公式End!