《数值积分方法.pptx》由会员分享,可在线阅读,更多相关《数值积分方法.pptx(28页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、取取左左端点端点矩形矩形近似近似 数值积分的数值积分的思想:思想:分割分割、近似、近似、求和求和取取右右端点端点矩形矩形近似近似 定积分定积分几何几何意义:意义:曲边梯形的面积曲边梯形的面积第1页/共28页 数值积分公式的数值积分公式的一般形式一般形式:其中其中求积求积节点节点求积求积系数系数仅与仅与求积节点求积节点有关有关求积公式的求积公式的截断误差截断误差或或余项余项:第2页/共28页5.1 插值型求积公式插值型求积公式思思想想用被积函数用被积函数 在区间在区间 上的上的插值多项式插值多项式近似代替计算近似代替计算作作n次次Lagrange插值多项式插值多项式:设已知函数设已知函数 在节点
2、在节点上的函数值上的函数值第3页/共28页其中其中余项余项第4页/共28页则有数值积分公式则有数值积分公式 这是用插值函数代替被积函数导出的定积分近似这是用插值函数代替被积函数导出的定积分近似计算公式,称为计算公式,称为插值型数值积分公式插值型数值积分公式。第5页/共28页n=1时的求积公式时的求积公式一、梯形一、梯形公式公式第6页/共28页用用梯形梯形面积近似面积近似 这是用线性插值函数代替被积函数导出的定积分近这是用线性插值函数代替被积函数导出的定积分近似计算公式,称为似计算公式,称为梯形数值积分公式梯形数值积分公式。几何意义几何意义第7页/共28页截断误差:已知线性插值的截断误差为 积分
3、中值定理:积分中值定理:连连续、不变号续、不变号第8页/共28页n=2时的求积公式时的求积公式二、二、Simpson公式公式将将 a,b 二二 等分,等分节点等分,等分节点 x0=a,x1=(a+b)/2,x2=b 作为积分节点,构造二次作为积分节点,构造二次Lagrange插值多插值多项式项式L2(x):第9页/共28页 这是用二次插值函数代替被积函数导出的定积分近似计这是用二次插值函数代替被积函数导出的定积分近似计算公式,称为算公式,称为辛普森数值积分公式辛普森数值积分公式。几何意义:几何意义:几何意义:几何意义:第10页/共28页SimpsonSimpson积分公式的截断误差(定理):积
4、分中值定理:积分中值定理:连续、不变号连续、不变号第11页/共28页复合求积法 通常把积分区间等分成若干个子区间,在每个子区间上用低阶的求积公式(如梯形积分公式Simpson积分公式),对所有的子区间求和即得整个区间a,b上的积分公式,这种方法称为复合求积法。5.2 复合复合求积求积公式公式第12页/共28页5.2.1 5.2.1 复化梯形积分复化梯形积分 将将a,b分成若干小区间,在每个区间分成若干小区间,在每个区间xi,xi+1上用梯形积分公式,再将这些小区间上的数上用梯形积分公式,再将这些小区间上的数值积分累加起来,就得到区间值积分累加起来,就得到区间a,b上的数值积上的数值积分。这种方
5、法称为分。这种方法称为复化梯形积分复化梯形积分。计算公式计算公式 将将a,b n等分等分,h=xi+1-xi=(b-a)/n,xi=a+ih,i=0,1,2,n,第13页/共28页记为记为记为记为 T T(h h)或或或或 T Tn n(f f):):复化梯形复化梯形公式的几何意义公式的几何意义小梯形小梯形面积面积之和之和近似近似复化梯形复化梯形公式公式第14页/共28页复化梯形复化梯形公式的余项公式的余项设设由由介值介值定理定理余项估计式余项估计式第15页/共28页 计算公式计算公式 将将a,b 2m 等分等分,m 为积分子区间数,记为积分子区间数,记 n=2m,n+1为节点总数为节点总数,
6、h=xi+1-xi=(b-a)/n,xi=a+ih,i=0,1,2,n,5.2.2 复化复化Simpson公式:公式:第16页/共28页复化复化Simpson公式公式复化复化Simpson公式的几何意义公式的几何意义小抛物小抛物面积面积之和之和近似近似系数首尾为系数首尾为1,奇数点为,奇数点为4,偶数点为,偶数点为2第17页/共28页复化复化Simpson公式的余项公式的余项设设由由介值介值定理定理余项估计式余项估计式第18页/共28页例:例:分别利用复化分别利用复化梯形梯形公式、公式、复化复化Simpson公式公式计算计算积分积分 的近似值,要求按复化的近似值,要求按复化Simpson公公式
7、计算时误差不超过式计算时误差不超过 。解:解:首先来确定首先来确定步长步长复化复化Simpson公式的余项:公式的余项:其中其中第19页/共28页本题本题 的求法:的求法:由由归纳法归纳法知知第20页/共28页解不等式得解不等式得将区间将区间 8等分,分别采用复化等分,分别采用复化Simpson、梯形梯形公式公式 0 1/8 1/4 3/8 10.9973980.9896880.976727 1/2 5/8 6/8 7/8 10.9588510.9361560.9088580.8771930.841471第21页/共28页复化复化梯形梯形公式公式(n=8)复化复化Simpson公式公式(n=4
8、)第22页/共28页代数精度代数精度的判别方的判别方法法 如果求积公式如果求积公式对一切不高于对一切不高于m次的多项式都次的多项式都恒成立恒成立,而对于某个,而对于某个m+1次多项式次多项式不能精确成立不能精确成立,则称该求积公式具有则称该求积公式具有m次代数精度。次代数精度。定理定理 求积公式求积公式具有次具有次m代数精度的充要条件是代数精度的充要条件是 为为 时求积公式时求积公式精确成立精确成立,而,而 为为 时求积公式不能成为等式。时求积公式不能成为等式。5.3 数值积分公式的代数精度和数值积分公式的代数精度和 Gauss求积求积公式公式第23页/共28页例例2 见见p73的例的例5.5
9、第24页/共28页 Gauss求积求积公式公式一、一、Gauss积分问题的提法积分问题的提法 前述前述的的求积公式中求积节点是取求积公式中求积节点是取等距节点等距节点,求积系数计算方便,但,求积系数计算方便,但代数精度代数精度要受到要受到限制;限制;为了提高为了提高代数精度代数精度,需要适当选择求积节点,需要适当选择求积节点:当求积节点个数确定后,不管这些求积节点如何选当求积节点个数确定后,不管这些求积节点如何选 取,取,求积公式的求积公式的代数精代数精度度最高最高能达到多少?能达到多少?具有具有最高最高代数精度代数精度的求积公式的求积公式中求积节点如何选取?中求积节点如何选取?积分公式的积分公式的一般形式一般形式:第25页/共28页 形如形如 的的插值型插值型求积公式的代数精度最高不超过求积公式的代数精度最高不超过2n+1次。次。定理定理第26页/共28页 这样由方程组的4个方程就能求出4个未知数,得根据定理知三点插值型求积公式的代数精度为5,同理可以去验证三点高斯求积公式二、二、Gauss求积公式的应用求积公式的应用积分积分见书上的见书上的P75.第27页/共28页谢谢您的观看!第28页/共28页
限制150内