[2018年最新整理]6数值积分与微分.docVIP

第六章 数值积分与数值微分 教学内容： 1．机械求积公式与代数精度 2．牛顿—柯特斯求积公式 3．龙贝格求积公式 4．高斯求积公式 5．数值微分 教学重点： 代数精度，各种求积公式的构造、收敛性、收敛速度及误差估计 教学难点： 复化求积公式的误差估计，龙贝格算法，高斯求积公式 计划学时：5-7学时 授课提纲： 机械求积公式与代数精度 数值求积的必要性 在工程及实际工作中经常会遇到求定积分的问题。如求一块石棉瓦曲边的长度，假定石棉瓦曲边服从函数f(x)=sinx，从x=0到x=b，则曲边弧长 这是一个第二类椭圆积分。又如在微分方程初值问题 的特解时，可以把解表示为 这是一个积分方程，其中需要求积分 由微积分理论可知，只要被积函数在区间[a,b]连续，就可以使用牛顿—莱布尼兹(Newton-Leibniz)公式 计算。然而，在许多实际问题中，这种解析方法是无能为力的，或者是非常麻烦的。原因有三： 找不到原函数F(x)，如第一个例子中的第二类椭圆积分，又如 等，有些被积函数是以表格形式给出，没有有限的解析表达式。 (2) 虽然可以找到原函数F(x)，但它比被积函数复杂得多，而且有时难以给出最后的数值结果。 (3) 除一些特殊的无穷积分外，通常很难求无穷积分的值。 故采用数值方法求定积是很有必要的，同时也可以证明后面给出的数值求积公式是有效的。 数值求积的基本思想 对于积分 由积分中值定理，若，则有 这样，如何近似的代替平均高度是数值求积的根本问题。如 用代替，得到左矩形公式 用代替，得到右矩形公式 用代替，得到中矩形公式 用代替，得到梯形公式 用代替，得到辛卜生公式 一般地，我们可以在区间[a,b]上适当选取若干个节点，然后用的加权平均值作为平均高度的近似值，构造求积公式 去近似。称为求积节点，为求积系数，为近似求积公式。这类数值求积方法常称作机械求积，（6-6）称为机械求积公式。该方法避开了牛顿—莱布尼兹公式需要求原函数的困难。 一但、确定，近似求积公式也就确定了。即、的不同选取，得到不同的求积公式。之后需要讨论的问题就是求积公式的代数精度、收敛性及稳定性。 求积公式的代数精度 不同的数值求积公式好坏的度量，是将求积公式应用在多项式上进行比较的，希望求积公式对次数尽可能高的多项式能精确成立，这就是代数精度的概念。 定义6-1 若某个求积公式对任意次数不超过m的多项式能精确成立，而对m+1多项式不能精确成立，则称该求积公式具有m次代数精度。 定理6-1 求积公式具有m次代数精度的充分必要条件是公式对精确成立，而对不精确。 证明 必要性 若求积公式具有m次代数精度，即对任意的数 由积分性质，有 故 由设求积公式对某个m+1次多项式不精确，即 由于 即有 即 充分性 若求积公式对精确成立，即 设为任一不超过m次的多项式，则有 即 又由于 故 。 该定理给出了验证求积公式的代数精度的有效手段，即验证求积公式对精确成立，而对不精确成立，则公式就具有m次代数精度。 不难验证梯形公式和辛卜生公式具有1次和3次代数精度。 试确定一个至少具有2次代数精度的求积公式 解 为使求积公式至少具有2次代数精度，公式必须对成立，即 解之得 ，， 即有 将代入求积公式两端，由于 ，故该求积公式 仅具有2次代数精度。 试确定求积公式中的参数，使公式具有尽可能高的代数精度。 解 令分别代入求积公式，并使等号成立，得 解之得 ，， 由于，而， 故求积公式仅有3次代数精度。 插值型求积公式 给定一组节点，且已知函数在这些节点处的函数值，容易得到以上述节点为插值基点的拉格朗日插值函数，令 记 得到插值型求积公式 由插值公式的余项定理容易得到插值型求积公式的余项 式中与变量x有关，。 推论6-1 对于插值型公式的求积系数，有 定理6-2 形如（6-7）的求积公式至少具有n次代数精度的充分必要条件是该公式是插值型求积公式。 证明 充分性 若公式是插值型求积公