- 1、本文档共30页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
实验二定积分的近似计算 数学软件与数学实验 教学课件
数学实验 实验二 定积分的近似计算 问题背景和实验目的 矩形法 定积分的定义: 矩形法 定积分的近似: 左点法、右点法和中点法 矩形法举例 矩形法举例 定积分几何意义 梯形法 梯形法 梯形法举例 抛物线法 抛物线法 抛物线法 抛物线法 抛物线法 Matlab 计算定积分函数介绍 梯形法:trapz trapz 举例 抛物线法 quad 举例 二重积分的计算 抛物线法计算二重积分: dblquad dblquad 举例 dblquad 举例 符号积分 int 举例 int 举例 其它相关函数 数值实验 数值实验 上机作业 * * 实验二、定积分的近似计算 定积分计算的基本公式是牛顿-莱布尼兹公式。但当被积函数的原函数不知道时,如何计算?这时就需要利用近似计算。特别是在许多实际应用中,被积函数甚至没有解析表达式,而是一条实验记录曲线,或一组离散的采样值,此时只能用近似方法计算定积分。 本实验主要研究定积分的三种近似计算算法:矩形法、梯形法和抛物线法。同时介绍 Matlab 计算定积分的相关函数。 实验二、定积分的近似计算 n 充分大,?x 充分小 通常我们取 左点法 右点法 中点法 点 可以任意选取,常见的取法有: 左端点 ,右端点 和中点 。 步长 节点 右点法: 中点法: 左点法: 相应的 Matlab 程序见 fuluA.m 解: == h =1/100=0.01, xi = i*h, a=0, b=1, n=100 例:用不同的矩形法计算下面的定积分 ( 取 n=100 ), 并比较这三种方法的相对误差。 左点法: 右点法: 中点法: (i = 0,1,2,...,100) 理论值: 左点法相对误差: 误差分析 右点法相对误差: 中点法相对误差: 不同的方法有不同的计算精度 有没有更好的近似计算定积分的方法 ? 曲边小梯形的面积可以由直边小梯形的面积来近似 整个曲边梯形的面积: 如果我们 n 等分区间 [a,b],即令: 则 == 梯形公式 梯形公式与中点公式有什么区别 ? Matlab 程序见 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) 公式 抛物线法公式 Matlab 程序见 fuluC.m == 例:用抛物线法计算下面定积分 ( 取 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,b,tol) f = f(x) 为被积函数,[a,b] 为积分区间,tol 为计算精度 将自变量看成是向量 抛物线法:quad 不用自己分割积分区间 可以指定计算精度,若不指定,缺省精度是 10-6 精度越高,函数运行的时间越长 此处的函数 f 是数值形式,应该使用数组运算,即 点运算:.*,./ ,.\ ,.^ 注: 解: quad(1./(1+x.^2),0,1) quad(1./(1+x.^2),0,1,1e-10) quad(1./(1+x.^2),0,1,1e-16) 函数表达式一定要用 单引号 括起来! 涉及的运算一定要用 数组运算! 例:用 quad 计算
您可能关注的文档
- 实习3 多进程(线程)快速排序 操作系统高级 教学课件.ppt
- 实习4 内存管理 操作系统高级 教学课件.ppt
- 实习指导书 水质工程学.doc
- 宝贝成长记录相册PPT模板.pptx
- 实体模型 [ppt]Mastercam+V10基础教程.ppt
- 实例_动物种群的相互竞争与相互依存的模型.doc
- 宝贝baby成长纪念册可爱宝宝相册PPT模板.pptx
- 实功率有限信号相关函数的定义 信号与线性系统分析(4版)电子教案.ppt
- 实例分析 [ppt]Mastercam V10基础教程.ppt
- 实例与练习 栅格数据的空间分析 ArcGIS地理信息系统基础实验练习.pdf
- 2024年中国钽材市场调查研究报告.docx
- 2024年中国不锈钢清洗车市场调查研究报告.docx
- 2024年中国分类垃圾箱市场调查研究报告.docx
- 2024年中国水气电磁阀市场调查研究报告.docx
- 2024年中国绿藻片市场调查研究报告.docx
- 2010-2023历年初中毕业升学考试(青海西宁卷)数学(带解析).docx
- 2010-2023历年福建厦门高一下学期质量检测地理卷.docx
- 2010-2023历年初中数学单元提优测试卷公式法(带解析).docx
- 2010-2023历年初中毕业升学考试(山东德州卷)化学(带解析).docx
- 2010-2023历年初中毕业升学考试(四川省泸州卷)化学(带解析).docx
最近下载
- 内燃机车钳工(中级)职业鉴定理论考试题及答案.doc VIP
- 聚焦贯彻落实新时代党的组织路线存在不足.docx VIP
- 研究生学术英语读写教程(于华)课后习题答案.pdf
- 跨学科实践活动2 自制净水器(同步课件)-【上好课】2024-2025学年九年级化学上册同步精品课堂(沪教版2024).pptx VIP
- 抗凝剂皮下注射相关并发症及处理对策.docx
- 第十章 甲状腺功能减退症PPT课件.ppt
- 《红楼梦的饮食文化》.docx VIP
- PICMG_COMe R3.0原版标准完整.docx
- ERAS加速康复外科在产科应用进展.pptx VIP
- 河南省南阳市六校2021-2022学年高一上学期第一次联合测试地理试题.pdf VIP
文档评论(0)