第七讲-定积的近似计算.pptVIP

数学实验 定积分的近似计算 问题背景和实验目的 主要内容 矩形法 定积分的定义： 矩形法 定积分的近似： 左点法、右点法和中点法 矩形法举例 矩形法举例 定积分几何意义 梯形法 梯形法 梯形法举例 抛物线法 抛物线法 抛物线法 抛物线法 抛物线法 Matlab 函数 总结 Matlab 计算定积分函数介绍 梯形法：trapz trapz 举例 抛物线法 quad 举例 dblquad 抛物线法计算二重积分： dblquad dblquad 举例 dblquad 举例 int 符号积分 int 举例 int 举例 其它相关函数 数值实验 数值实验 triplequad 抛物线法计算三重积分： triplequad triplequad 举例 triplequad 举例 上机作业 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 函数近似计算二重积分 triplequad(f,a,b,c,d,e,f,tol) tol 为计算精度，若不指定，则缺省精度为 10-6 f 可以是： 字符串；inline 定义的内联函数；函数句柄 [a,b] 是 第一积分变量 的积分区间， [c,d] 是 第二积分变量 的积分区间 [e,f]是第二积分变量 的积分区间 f(x,y,z) 中关于前两个自变量的运算是数组运算， 即把 x， y 看成是向量，z 看成是标量。 也可以全部采用数组运算 例2：计算三重积分 triplequad(inline(4*x.*y+3*x.^2+z^2),-1,1,0,2,0,1) syms x y z int(int(int(4*x*y+3*x^2+z^2,x,-1,1),y,0,2),z,0,1) triplequad(@(x,y,z)4*x.*y+3*x.^2+z^2 , -1,1, 0, 2,0,1) 指定 x、y, z 分别是第一、二、三积分变量 triplequad(inline(4*x.*y+3*x.^2+z^2), -1,1, 0, 2,0,1) 被积函数 f (x,y,z) 的另一种定义方法：匿名函数 例2：计算三重积分 梯形数值积分命令 trapz( ) clear x=0:pi/100:pi; y=sin(x); trapz(x,y) plot(x,y,b*) clear x=sort(rand(1,101)*pi); y=sin(x); trapz(x,y) plot(x,y,rd) rand(1,101)产生101个均匀随机数,每个数都介于0-1之间 辛卜生求数值积分命令quad() clear fun=inline(1./(x.^3-2*x-5)) ezplot(fun,[0,2]) [q,n]=quad(fun,0,2) 二重积分dblquad()与三重积分 fun=inline(y*sin(x)+x*cos(y)) Q=dblquad(fun,pi,2*pi,0,pi) [x,y]=meshgrid(pi:.1:2*pi,0:.1:pi); z=fun(x,y); mesh(x,y,z) 1.分别用梯形法与抛物线法，计算 取n = 120，并常识直接使用函数trapz()、quad()进行计算求解，比较结果的差异。 2.试计算定积分 注意：可以运用trapz()、quad()或附录程序求解吗？ 3.学习fuluBsum.m的程序设计方法，尝试用函数sum改写附录A和附录C的程序，避免for循环。 用Matlab演示 Matlab演示 用Matlab演示 Matlab演示 用Matlab演示 * * 定积分的近似计算 定积分计算的基本公式是牛顿－莱布尼兹公式。但当被积函数的原函数不知道时，如何计算？这时就需要利用近似计算。特别是在许多实际应用中，被积函数甚至没有解析表达式，而是一条实验记录曲线，或一组离散的采样值，此时只能用近似方法计算定积分。 本实验主要研究定积分的三种近似计算算法：矩形法、梯形法和抛物线法。同时介绍 Matlab 计算定积分的相关函数。 矩形法 梯形法 抛物线法 数值积分的常见算法 Matlab 求积分函数 数值积分函数：trapz、quad、d