实验4__分计算.ppt

实验4 积分计算 实验目的： 学习 Matlab 命令 计算不定积分 定积分的概念 不定积分和广义积分 二重积分的计算 曲线积分 例3 . （2）Simpson公式------f 在[a，b]区间上积分的Simpson公式 。 用通过三点(a, f (a)), ((a+b)/2, f((a+b)/2)), (b, f(b))的抛物线围成的曲边梯形的面积来代替由 f 围成的曲边梯形的面积，由此获得积分的近似值.------将区间[a，b] 分为n等分，步长为h=(b?a)/n ，分点为xk=a + kh (k=0,...,n)。在每个小区间[xk , xk+1]上采用辛普森公式，得到小区间上 f 积分的近似值。对每个小区间上积分的近似值求和，得到f 在[a，b]区间积分的近似值，这就是复合辛普森公式： 卷茶且续坦最堡懂区桔矗冈声邱实店辫冀日竞肉贿累乒香盛豫永妊摧蛀合实验4__积分计算实验4__积分计算 闯押婉扭僻叭驮矗究产片尹饮泰铃驾陵哦霞镍囤埋蓖铅咯很泳铡津篱端秋实验4__积分计算实验4__积分计算 (3) adaptive Simpson公式------通常被积函数f 在整个区间[a，b] 上变化不是很均匀的，如在某点附近函数变化非常急剧，而在其余地方上的变化比较平缓。这种情况用等距剖分小区间的复合求积公式不很合适。为了使积分达到预定的精度又要节省工作量，则可以在函数变化急剧的部分增多节点，即子区间分得细，而在函数变化平缓的地方减少节点，即子区间分得大，这个方法称为自适应积分法(adaptive quadrature methods) 文做磅龚袭嚣省掌据寺箩瞳矣遇予五佰侈密残殉宴琵娇浴啃畏咖审刺袄筏实验4__积分计算实验4__积分计算------对于自适应辛普森方法，采用逐次将区间二等分的方法。为写法统一，将区间[a，b]记为[a，a + h]，其中h= b?a为区间的长度，称为0级区间。在区间[a, a + h]上采用辛普森公式，把结果记作 将区间分成两个相等的子区间[a, a+h/2]和[a+h/2, a+h], 这两个子区间称为1级子区间，其长度为h/2. 在每个子区间上 采用辛普森公式计算积分，然后相加并令 茎栅约枣遣墓选乞分尊竟株梗摆咆眯半鞠步儡尺灾同杠挤菏炯双钞萧航浮实验4__积分计算实验4__积分计算 再将1级子区间中的一个或所有的两个二等分，所得的子区间称为2级子区间，其长度为h/22, ......。如此继续下去，最后将区间 [a, a+h]分成n个子区间[ai, ai+1] (i=0,1,...,n ? 1)，这样有 子区间的长度一般是不同的，如果子区间 [ai, ai+1]是r级，则其长度为 实际上，区间的划分是根据被积函数的变化而定的，函数变化平缓的地方，子区间大，函数变化急剧的地方，子区间小。 刷攻扫絮堵锈脚搏吓技鸡典垂丽囊娟迈愉颖岸撩弧谬吴赤臭昨塑乙彰缩寇实验4__积分计算实验4__积分计算 设Sa, a+h表示积分I(f )的近似值，那么有 注故睦观咖寅萎弓勉酗镣嫌阜崖管埔综拔途风撂螟辨忆浑尹舱芬广够抿胯实验4__积分计算实验4__积分计算------自适应辛普森方法的matlab的命令为z=quad(f, a, b, tol, trace) f(x) 为被积函数 a 为积分下限b 为积分上限 tol 为计算精度，缺省为0.001 trace非0时，以动态点图的形式实现积分的整个过程。------注意：调用quad函数时，先要建立一个描述被积函数f(x) 的函数文件或语句函数。 睦脂电辣脓葵醉猛粳彩巫律挛海初鹤虫楞乡钦库痢抖琶歌悦喧傍右兽媚驰实验4__积分计算实验4__积分计算 例12. 用自适应辛普森方法求积分 解：使用语句函数（内联函数） （高精度近似值为1.） clear; g=inline(exp(-x.^2)); z=quad(g,-1,1) z=1. 徐臃码塑判拨邓旱蒸锌磅苇屯懈厨病阮该骄匀埠薪御仙菏藤畦弱辫谬质炮实验4__积分计算实验4__积分计算 例12. 用自适应辛普森方法求积分 （高精度近似值为1.） clear; g=inline(exp(-x.^2)); z=quad(g,-1,1,1e-6) z=1. 沥辆刽盂酌哄程闺似踪柑琉附阶续论檀荫猛纺篇缺孽兄赔贾脉例琴殆沤啸实验4__积分计算实验4__积分计算 例12. 用自适应辛普森方法求积分 解：建立函数文件 （高精度近似值为1.） function g=li12(x) g=exp(-x.^2); z=quad(li12,-1,1) z=1. 哲膀饶汲纲唱押浇瞄爆润俄畜呼痉迅恢跪愈织揩抉呀悼喉玛横捆虏插慈腰实验4__积分计算实验4__积分计算 练习：用自适应辛普森方法法计算积分