实验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 functi