数值分析李庆杨数值积分与数值微分.pptxVIP

来求定积分。前面公式虽然在理论上或在解决实际问题中都起了很大的作用,但它并不能完全解决定积分的计算问题。因为定积分的计算常常会碰到以下三种情况: (1)被积函数f(x)的原函数F(x)不易找到。许多很简单的函数,例如 等,其原函数都不能用初等函数表示成有限形式。 ; (2)被积函数f(x)没有具体的解析表达式。其函数关系由表格或图形表示,无法求出原函数。 (3)尽管f(x)的原函数能表示成有限形式但其表达式相当复杂。例如定积分 的被积函数 的原函数就比较复杂,从数值计算角度来看,计算量太大。 ; 如图4.1,若用左矩形近似地代替曲边梯形,则得到左矩形公式 ? ; 图 4.1 ; 如图4.2,若用梯形的面积近似地代替曲边梯形的面积,则得到计算定积分的梯形公式 ;图 4.2;二、代数精度的概念;(1.4);这就是梯形公式(1.1),它表明利用线性方程组(1.4)推出的求积公式，与用通过两点(a,f(a))与(b,f(b))的直线近似曲线y=f(x)得到的结果一致。当f(x)=x2时(1.4)式的第三个式子不成立，因为; 在方程组(1.4)中如果节点xi及系数Ai都不确定，那么方程组(1.4)是关于xi及Ai(i=0,1,…,n)的2n+2个参数的非线性方程组。此方程组当n1时求解是很困难的，但当n=0及n=1的情形还可通过求解方程组(1.4)得到相应的求积公式。;三、插值型求积公式;四、求积公式的余项;（1.9）;例.给定形如 的求积公式，试确定系数 ，使公式具有尽可能高的代数精确度，并求此求积公式的余项。;解：根据题意可令 分别代入求积公式使它精确成立得 解得 于是有 当 时 . 而上式右端为1/3，故公式对 不精确成立，其代数精确度为2. 故余项表达式为 令 得 ，于是有 故得 ;五、求积公式的收敛性和稳定性;§4.2 牛顿 ― 柯特斯(Newton―Cotes) 公式 建立数值积分公式最基本的思想是选取一个既简单又有足够精度的函数φ(x),用φ(x)代替被积函数f(x),于是有 现用第2章介绍的插值多项式Pn(x)来代替被积函数f(x),即有 ; 利用拉格朗日插值多项式 ; 这里yi=f(xi),对式(2―1)两边积分得 ; 为牛顿 ― 柯特斯求积公式,Rn(f)为牛顿 ― 柯特斯求积公式的余项。 令 x=x0+sh , 0≤s≤n;(2―6) ; 称C(n)i为柯特斯求积系数。 很显然,当n=1时,可算得 ; 这是梯形公式。 当n=2时,可得 ; 这是抛物线(辛普森Simpson)公式。 当n=3时, ; 代入(2―5)式得到求积公式 ; 表 4―1 ; 柯特斯系数C(n)i仅与n和i有关,与被积函数f(x)无关,且满足 ; 例1 试分别用梯形公式和抛物线公式计算积分 解 利用梯形公式 ; 原积分的准确值 ; 2.2 误差估计 现对牛顿 ― 柯特斯求积公式所产生的误差作一个分析。由式(2―4),牛顿 ― 柯特斯求积公式的余项为 易知,牛顿 ― 柯特斯求积公式(2―5)对任何不高于n次的多项式是准确成立的。这是因为 f(n+1)(ξ)≡0 故 Rn(f)≡0; 一般说来,若某个求积公式对于次数不高于m的多项式都准确成立(即Rn(f)≡0),而对于某一次数为m+1的多项式并不准确成立(即Rn(