实验二、定积分的近似计算.pptVIP

* 第1页，共30页，编辑于2022年，星期六 定积分计算的基本公式是牛顿－莱布尼兹公式。但当被积函数的原函数不知道时，如何计算？这时就需要利用近似计算。特别是在许多实际应用中，被积函数甚至没有解析表达式，而是一条实验记录曲线，或一组离散的采样值，此时只能用近似方法计算定积分。 本实验主要研究定积分的三种近似计算算法：矩形法、梯形法和抛物线法。同时介绍 Matlab 计算定积分的相关函数。 问题背景和实验目的 定积分的近似计算 * 第2页，共30页，编辑于2022年，星期六 矩形法 梯形法 抛物线法 数值积分的常见算法 主要内容 Matlab 求积分函数 数值积分函数：trapz、quad、dblquad 符号积分函数：int * 第3页，共30页，编辑于2022年，星期六 定积分的定义 定积分的近似 * 第4页，共30页，编辑于2022年，星期六 矩形法 n 充分大，?x 充分小 通常我们取 左点法 右点法 中点法 点 可以任意选取，常见的取法有： 左端点 , 右端点 和中点 。 定积分的近似： * 第5页，共30页，编辑于2022年，星期六 步长 节点 矩形法 左点法 右点法 中点法 * 第6页，共30页，编辑于2022年，星期六 矩形法举例 例：用不同的矩形法计算下面的定积分 ( 取 n=100 )， 并比较这三种方法的相对误差。 左点法： 右点法： 中点法： 解： h =1/n=0.01, xi = i*h, a=0, b=1, n=100 (i = 0, 1, 2, ..., 100) * 第7页，共30页，编辑于2022年，星期六 理论值： 左点法相对误差： 相对误差分析 矩形法举例 右点法相对误差： 中点法相对误差： 不同的算法有不同的计算精度 有没有更好的近似计算定积分的方法 ? * 第8页，共30页，编辑于2022年，星期六 定积分几何意义 * 第9页，共30页，编辑于2022年，星期六 曲边小梯形的面积可以由直边小梯形的面积来近似 整个曲边梯形的面积： 梯形法 * 第10页，共30页，编辑于2022年，星期六 如果我们 n 等分区间 [a,b]，即令： 则 == 梯形公式 梯形法 梯形公式与中点公式有什么区别 ? * 第11页，共30页，编辑于2022年，星期六 解： == 例：用梯形法计算下面定积分 ( 取 n=100 )，并计算相对误差 梯形法举例 a=0, b=1, n=100, f (x) = 1/( 1+x2 ) == h =1/100=0.01, xi = i*h, yi = f (xi) 相对误差： * 第12页，共30页，编辑于2022年，星期六 2n 等分区间 [a,b] ，得 用抛物线代替该直线，计算精度是否会更好？ 计算每个节点上的函数值： 抛物线法 在区间 [x0, x2] 上，用过以下三点 的抛物线来近似原函数 f (x) 。 * 第13页，共30页，编辑于2022年，星期六 设过以上三点的抛物线方程为： 则在区间 [x0, x2] 上，有 y = ? x2 + ? x + ? = p1(x) 抛物线法 * 第14页，共30页，编辑于2022年，星期六 同理可得： 相加即得： 抛物线法 * 第15页，共30页，编辑于2022年，星期六 整理后可得： 或 辛卜生 (Simpson) 公式 抛物线法公式 抛物线法 * 第16页，共30页，编辑于2022年，星期六 == 例：用抛物线法计算下面定积分 ( 取 n=100 )，并计算相对误差 解： a=0, b=1, n=100, yi = f (xi) = 1/( 1+xi2 ) 抛物线法 相对误差： * 第17页，共30页，编辑于2022年，星期六 矩形法 梯形法 抛物线法 数值积分的常见算法 Matlab 函数 Matlab 求积分函数 数值积分函数：trapz、quad、dblquad 符号积分函数：int * 第18页，共30页，编辑于2022年，星期六 trapz(x,y)x 为分割点（节点）组成的向量，y 为被积函数在节点上的函数值组成的向量。 trapz trapz 梯形法 * 第19页，共30页，编辑于2022年，星期六 前面的做法 例：用梯形法计算下面定积分 ( 取 n=100 ) 解： a=0, b=1, n=100, yi = f (xi) = 1/( 1+xi2 ) x=0:1/100:1; y=1./(1+x.