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

定积分的近似计算 hjTang@xidian.edu.cn 问题背景和实验目的 矩形法 定积分的定义： 矩形法 定积分的近似： 左点法、右点法和中点法 矩形法举例 矩形法举例 定积分几何意义 梯形法 梯形法 梯形法举例 抛物线法 抛物线法 抛物线法 抛物线法 抛物线法 Matlab 计算定积分函数介绍 梯形法：trapz trapz 举例 抛物线法 quad 举例 二重积分的计算 抛物线法计算二重积分： dblquad dblquad 举例 dblquad 举例 三重积分的计算 抛物线法计算二重积分： triplequad triplequad 举例 triplequad 举例 int 符号积分 int 举例 int 举例 其它相关函数 数值实验 数值实验 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演示 用Matlab演示 * * 定积分的近似计算 定积分计算的基本公式是牛顿－莱布尼兹公式。但当被积函数的原函数不知道时，如何计算？这时就需要利用近似计算。特别是在许多实际应用中，被积函数甚至没有解析表达式，而是一条实验记录曲线，或一组离散的采样值，此时只能用近似方法计算定积分。 本实验主要研究定积分的三种近似计算算法：矩形法、梯形法和抛物线法。同时介绍 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) 。 设过以上三点的抛物线方程为： 则在区间 [x0, x2] 上，有 y = ? x2 + ? x + ? = p1(x) 同理可得： 相加即得： 整理后可得： 或 辛普森 (Simpson) 公式 抛物线法公式 == 例：用抛物线法计算下面定积分 ( 取 n=100 )， 并计算相对误差 解： a=0, b=1, n=100, yi = f (xi) = 1/( 1+xi2 ) 相对误差： trapz(x,y)x 为分割点（节点）组成的向量，y 为被积函数在节点上的函数值组成的向量。 Matlab 近似计算定积分的相关函数 前面的做法 例：用梯形法计算下面定积分 ( 取 n=100) 解： a=0, b=1, n=100, yi = f (xi) = 1/( 1+xi2 ) x=0:1/100:1; y=1./(1+x.^2); trapz(x, y) trapz函数 trapz(x,1./(1+x.^2)) quad(f,a