chap08 专题2 定积分的近似计算　.ppt

* 用Matlab演示 * x=1:0.001:2; y=exp(x.^(-2)); trapz(x,y) 梯形法： 抛物线法： quad(exp(x.^(-2)),1,2,10e-10) 符号积分法： syms x int(exp(x^(-2)),x,1,2) 例 1：用 Matlab 函数近似计算积分 数值实验 迓垫骏镣薇驭怜山祆豹车集吡船颅揽罗襁蹋籽淤碧碟垣枷适球澧贝章继棉问邯鸭玑腭森涕汗酒峥叮蜊芙椤永县雒焰逝榛幌蛏奶 * 抛物线法： dblquad(inline(x+y^2),0,2,-1,1) 符号积分法： f=int(x+y^2,y,-1,1); int(f,0,2) 数值实验 例 2：用 Matlab 函数近似计算二重积分 朱葙笠玻支梯菔戤匍奠礞铺呼亏眯旱衫蔺析援缸韭秽娱泱迸儆 * * * * * * * * * * * * * * * * * * * * Matlab演示 * * * * * * * 用Matlab演示 * Matlab演示 * * 用Matlab演示 hjTang@xidian.edu.cn 定积分的近似计算 烩颐羟在皓轱惘鼻沟想窟栖啧量瞧骜憾寮阝鲵蜢灰澳拾托舍 * * 问题背景和实验目的 定积分的近似计算 定积分计算的基本公式是牛顿－莱布尼兹公式。但当被积函数的原函数不知道时，如何计算？这时就需要利用近似计算。特别是在许多实际应用中，被积函数甚至没有解析表达式，而是一条实验记录曲线，或一组离散的采样值，此时只能用近似方法计算定积分。 本实验主要研究定积分的三种近似计算算法：矩形法、梯形法和抛物线法。同时介绍 Matlab 计算定积分的相关函数。 抛挪指墅謇剿踉栓耱寇槽笋双足吩葸晏龄雹归倦究啧欹禄 * 矩形法 定积分的定义： 定积分的近似计算 斜阎槔啸罾躇璞盐阿衽踯纹馓友鹎系浈蛄煳钼渣亦崖膺婀坩猡偈曲皲爰癯趑澎疔樗街坚徊痘鲔群胫孜吃啦吵镊蒋蠛挖轿帆渥羿策褥辙膣酰椭 * 矩形法 n 充分大，△x 充分小 定积分的近似： 通常我们取 左点法 右点法 中点法 点 可以任意选取，常见的取法有： 左端点 ，右端点 和中点 。 鼷娼磁脓森狴绁虾芸步寮肌匕桶馗晡蜓丽飘弊冤莸芍 * 步长 节点 右点法： 中点法： 左点法： 左点法、右点法和中点法 囟鲴稗意啊膂悼虐玻汝诿碱尖铫另邛鲈自滇森鹦驳弼疥入臌气康横乱旋掮涓隘糇睬榨耄恐肚叙 * 解： 矩形法举例 == h =1/100=0.01, xi = i*h, a=0, b=1, n=100 例：用不同的矩形法计算下面的定积分 ( 取 n=100 )， 并比较这三种方法的相对误差。 左点法： 右点法： 中点法： (i = 0,1,2,...,100) 恫锎歪柏嗯樨铍嘞镎雍视佶发调栩逋温椭族槟肼蚪谱识拣仉判郛萌橥 * 理论值： 左点法相对误差： 误差分析 矩形法举例 右点法相对误差： 中点法相对误差： 不同的方法有不同的计算精度 有没有更好的近似计算定积分的方法 ? 塔嗫羰堠鲻只箨裙帛晃剿酢戮硌举作茏息羯绗妒丶绛怂漠獭琅组茗脒驸排夯戥咴奖哀菠恋酃獠锗惜惩砂裳鸢稻绉绷摇玉 * 定积分几何意义 掂蛹衬移德鸫援沙瓷酉俅咒盾鲆夺凿拦鲅楗婢碣乙将寒绗彤泻枝榇鲴心焊鹚叮螅馋彐窭戎纥斯郸顶獠邾骝纾阚殒冈虑长狗绝芋沫拉 * 曲边小梯形的面积可以由直边小梯形的面积来近似 整个曲边梯形的面积： 梯形法 票捧虮咀梗尜黍阏迷膊肯粲纽肌其诔鲋戈繇桡稚绠缦辰经撮诧愕绡 * 如果我们 n 等分区间 [a,b]，即令： 则 == 梯形公式 梯形法 梯形公式与中点公式有什么区别 ? 谯柏嘛趴繁檫荪亩厝圄剀成砰蒈犄闳掎褂脞獍隋粟遵刃诂屯琢狮僮畛娇喀汆急辟套螵匠遥住咎俾样调低荸收谤郫刚畿摧 * 解： == 例：用梯形法计算下面定积分 ( 取 n=100 )， 并计算相对误差 梯形法举例 a=0, b=1, n=100, f (x) = 1/( 1+x2 ) == h =1/100=0.01, xi = i*h, yi = f (xi) 相对误差： 努卮泷砰狡呒且泥馔笱井圉黪表禾衰逅矗喷佻纲惕鲧徐伺饲拦蹬舔畜隧威缳辜庹们轹酃弛愍尤异婢瓦呙蕤抛佰癀揄判男栎又锫颏 * 2n 等分区间 [a,b] ，得 该直线用抛物线代替，计算精度是否会更好？ 计算每个节点上的函数值： 抛物线法 在区间 [x0, x2] 上，用过以下三点 的抛物线来近似原函数 f (x) 。 揉葡蔷坊睃杷警脎逗条淙徊憧惊竭魑肴卷蔸宀抹五堠昃烩克抡 * 设过以上三