计算方法教案详解.docVIP

计 算 方 法 教 案 姓名：魏哲哲 学号：120444072 班级：12级金融班 院系：数学科学学院 数值积分 数值积分： 求定积分的近似值的数值方法。即用被积函数的有限个抽样值的离散或加权平均近似值代替定积分的值。 本节内容： 基础知识 牛顿—柯特斯公式 一、基础知识 1、问题的提出 计算连续函数在区间[a,b]上的定积分，可以用牛顿—莱布尼兹（Newton-Leibnitz）公式： =F(b)-F(a) 其中F(x)是f(x)在区间[a,b]上的一个原函数。 在某些情况下，不适合用牛顿—莱布尼兹公式求定积分： 通过观测等方法得到f(x)上某些离散的点，不知道或不存在f(x)的具体解析表达式。 f(x)的原函数F(x)不能用有限形式表示。 求原函数F(x)较麻烦，但利用计算机等工具，用数值积分的方法求解较方便。数值积分方法求得的结果虽然有误差，但可以控制在要求的误差范围之内。实际工程问题一般不需要理论上无误差的精确解，能得到误差满足要求的近似解就可以了。 数值积分的主要思想，是避开求f(x)的原函数F(x)，只由f(x)上若干离散的点，用对应的数值积分算法，来求f(x)在[a,b]上定积分的近似值。 2、数值积分公式 =(b-a)f(ξ) f(ξ)为f(x)在区间[a,b]上的平均高度。 如果知道了平均高度f(ξ)，就能够求出积分 。现在的问题是，如何由f(x)上若干个点f(xi)，(i=0,1,2,……,n)，构造1个表达式，替代定积分中值定理中的平均高度f(ξ) 。 数值积分方法就是在积分区间[a,b]上选取若干节点x0,x1,……,xn，用f(x0),f(x1),……,f(xn)的线性组合作为积分的近似值。数值积分公式的一般形式为： ≈ f(xi) R[f]= － f(xi) 其中xi称为求积节点，Ai称为求积系数，R[f]为数值积分公式的余项。 Ai是与求积节点对应的常数，它不依赖被积函数f(x)的具体形式。这一类求积方法避开了求f(x)的原函数F(x)，易于用计算机实现，因此又称为机械求积。 1.中矩形公式 I=(b-a)f((a+b)/2) 取积分区间[a,b]的中点x=(a+b)/2近似地代替ξ，用f((a+b)/2)取代平均高度f(ξ)。中矩形公式是1点公式，又是开公式（求积节点不包含积分区间的端点）。1点公式还有左矩形公式和右矩形公式。 2.左矩形公式 I=(b-a)f(a) 取积分区间[a,b]的端点x=a近似地代替ξ，用f(a)取代平均高 度f(ξ)。 3. 左矩形公式 I=(b-a)f(b) 取积分区间[a,b]的端点x=b近似地代替ξ，用f(b)取代平均高度f(ξ)。 左矩形公式和右矩形公式是半开半闭公式（求积节点仅包含积分区间1侧的端点）。 4. 梯形公式 I=(b-a)(f(a)+f(b))/2) 取f(a)、f(b)的平均值f(a)+f(b))/2近似地代替平均高度f(ξ) 。梯形公式是2点公式，又是闭公式（求积节点包含积分区间的端点）。 3、代数精度 代数精度是间接地反应某一数值积分公式对被积函数逼近能力的一个参数。 定义：如果求积公式对所有不超过m次的代数多项式函数都能精确成立，而对m＋1次的代数多项式函数至少有1个不能精确成立，则称此求积公式具有m次代数精度。 若被积函数为n次代数多项式函数，求积公式有m次代数精度，且m≥n，则求得的积分值误差为0。 若被积函数不是代数多项式函数，那么不能保证有限次代数精度的数值积分公式的余项为0。过高阶次代数精度的数值积分有时会出现很大误差，因此建议不使用代数精度超过6次的数值积分公式。当积分区间[a,b]较大时，可以采用复化求积公式等方法来减小误差。 定理：求积公式具有m次代数精度的充要条件是当被积函数f(x)=xr，r=0,1,2,……,m时精确成立，当被积函数f(x)=xm+1时不能精确成立。 中矩形公式和梯形公式具有1次代数精度，左矩形公式和右矩形公式具有0次代数精度。 4、插值型求积公式 插值与数值积分有密切联系。我们可以先对求积节点做插值，再对插值函数求积分，如果插值函数与被积函数足够接近，就可以把插值函数的积分近似为被积函数的积分。代数插值函数具有存在唯一性，易构造，易求积分，易用计算机实现，因此本章用第5章的方法构造代数插值函数。 用拉格郎日插值函数代替被积函数求积分： ≈= = f(xi) 式中 是与求积节点xi对应的常数，不依赖被积函数的具体形式，象这样构造出的求积公式称为插值型求积公式，它的余项函数为 R[f]= - = 当x∈[a,b]时