第5章 节 数值积分 计算方法 《代码优化》课件.pptVIP

第5章 数值积分 ;1.1 牛顿 ― 柯特斯(Newton―Cotes) 公式 ; 来求定积分。公式(5―1)虽然在理论上或在解决实际问题中都起了很大的作用,?但它并不能完全解决定积分的计算问题。因为定积分的计算常常会碰到以下三种?情况: (1)被积函数f(x)的原函数F(x)不易找到。许多很简单的函?数,例如 等,其原函数都不能用初等函数表示成有限形式。 ; (2)被积函数f(x)没有具体的解析表达式。其函数关系由表格或图形表?示,无法求出原函数。 (3)尽管f(x)的原函数能表示成有限形式但其表达式相当复杂。例如定?积分 的被积函数 的原函数就比较复杂,从数值计算角度来?看,计算量太大。 ; 如图5.1,若用?左矩形近似地代替曲边梯形,则得到左矩形公式 ? ;图 5.2; 1.1 牛顿 ― 柯特斯(Newton―Cotes) 公式 建立数值积分公式最基本的思想是选取一个既简单又有足够精度的函数φ(x),?用φ(x)代替被积函数f(x),于是有 现用第四章介绍的插值多项式Pn(x)来代替被积函数f(x),即有 ; 利用拉格朗日插值多项式 ; 这里yi=f(xi),对式(5―6)两边积分得 ; 为牛顿 ― 柯特斯求积公式,Rn(f)为牛顿 ― 柯特斯求积公式的余项。 令 x=x0+sh , 0≤s≤n;(5―11) ; 称C(n)i为柯特斯求积系数。 很显然,当n=1时,可算得 ; 这是梯形公式。 当n=2时,可得 ; 这是抛物线公式。 当n=3时, ; 代入(5―10)式得到求积公式 ; 柯特斯系数C(n)i仅与n和i有关,与被积函数f(x)无关,且满足 ; 例1 试分别用梯形公式和抛物线公式计算积分 解 利用梯形公式 ; 原积分的准确值 ; 表 5―1 ; 1.2 误差估计 现对牛顿 ― 柯特斯求积公式所产生的误差作一个分析。由式(5―9),牛顿 ― 柯特斯求积公式的余项为 易知,牛顿 ― 柯特斯求积公式(5―10)对任何不高于n次的多项式是准确成立的。这是因为 f(n+1)(ξ)≡0 故 Rn(f)≡0; 一般说来,若某个求积公式对于次数不高于m的多??式都准确成立(即Rn(f)≡0),而对于某一次数为m+1的多项式并不准确成立(即Rn(f)0),则称这一求积公式的代数精确度为m。 ? 牛顿 ― 柯特斯求积公式的代数精确度至少为n。通常在基点个数相等的情况下,代数精确度愈高,求积公式就愈精确。 定理1 (梯形公式的误差)设f(x)在区间［a, b］上具有连续的二阶导数,则梯形求积公式的误差为 ; 由于 ω1(x)=(x-a)(x-b) ; 在区间(a, b)内不变号,f″(ξ)是x的函数且在［a,b］上连续,故根据积分第二中值定理参见有关《数学分析》教材中“一元函数积分学第二中值定理”。 知,存在某一η∈(a, b)使 ; 定理2 (抛物线公式的误差)设f(x)在［a, b］上有连续的四阶导数,则抛物线公式的误差为 ;§2 复合求积公式 ; 相加后得 ; 因而 ; 例2 若用复合梯形公式计算积分 问积分区间要等分多少才能保证有五位有效数字 解 由余项(5―21)式 ; 由于原积分的准确值具有一位整数,因此要使近似积分值有五位有效数字,只需取n满足 ; 2.2 复合抛物线公式 类似复合梯形公式的做法,把区间［a,b］分成n个相等的子区间［x2i,x2i+2］(i=0,1,…,n