高等数学定积分的近似计算.pptVIP

精品课件 前面的做法 例：用梯形法计算下面定积分 ( 取 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)) trapz 举例 精品课件 quad(f,a,b,tol) f = f(x) 为被积函数，[a,b] 为积分区间，tol 为计算精度 将自变量看成是向量 不用自己分割积分区间 可以指定计算精度，若不指定，缺省精度是 10-6 精度越高，函数运行的时间越长 此处的函数 f 是数值形式，应该使用数组运算，即： .* ./ .\ .^ quad quad 抛物线法 用Matlab演示 用Matlab演示 精品课件 实验二 定积分的近似计算 数学实验 精品课件 定积分计算的基本公式是牛顿－莱布尼兹公式。但当被积函数的原函数不知道时，如何计算？这时就需要利用近似计算。特别是在许多实际应用中，被积函数甚至没有解析表达式，而是一条实验记录曲线，或一组离散的采样值，此时只能用近似方法计算定积分。 本实验主要研究定积分的三种近似计算算法：矩形法、梯形法和抛物线法。同时介绍 Matlab 计算定积分的相关函数。 问题背景和实验目的 定积分的近似计算 精品课件 矩形法 梯形法 抛物线法 数值积分的常见算法 主要内容 Matlab 求积分函数 数值积分函数：trapz、quad、dblquad 符号积分函数：int 精品课件 定积分的定义 定积分的近似 精品课件 矩形法 n 充分大，?x 充分小 通常我们取 左点法 右点法 中点法 点 可以任意选取，常见的取法有： 左端点 , 右端点 和中点 。 定积分的近似： 精品课件 步长 节点 矩形法 左点法 右点法 中点法 fuluA.m 精品课件 矩形法举例 例：用不同的矩形法计算下面的定积分 ( 取 n=100 )， 并比较这三种方法的相对误差。 左点法： 右点法： 中点法： 解： h =1/n=0.01, xi = i*h, a=0, b=1, n=100 (i = 0, 1, 2, ..., 100) 精品课件 理论值： 左点法相对误差： 相对误差分析 矩形法举例 右点法相对误差： 中点法相对误差： 不同的算法有不同的计算精度 有没有更好的近似计算定积分的方法 ? 精品课件 定积分几何意义 精品课件 曲边小梯形的面积可以由直边小梯形的面积来近似 整个曲边梯形的面积： 梯形法 精品课件 如果我们 n 等分区间 [a,b]，即令： 则 == 梯形公式 梯形法 梯形公式与中点公式有什么区别 ? fuluB.m 精品课件 解： == 例：用梯形法计算下面定积分 ( 取 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) 公式 抛物线法公式 抛物线法 fuluC.m 精品课件 == 例：用抛物线法计算下面定积分 ( 取 n=100 )，并计算相对误差 解： a=0, b=1, n=100, yi = f (xi) = 1/( 1+xi2 ) 抛物线法 相对误差： 精品课件 矩形法 梯形法 抛物线法 数值积分的常见算法 Matlab 函数 Matlab 求积分函数 数值积分函数：trapz、quad、dblquad 符号积分函数：int 精品课件 trapz(x,y)x 为分割点（节点）组成的向量，y 为被积函数在节点上的函数值组成的向量。 trapz trapz 梯形法 用Matlab演示 用Matlab演示 *