实验三定积分的近似计算

一、问题背景与实验目的

利用牛顿—莱布尼兹公式虽然可以精确地计算定积分的值,但它仅适用于被积函数的原函数能用初等函数表达出来的情形.如果这点办不到或者不容易办到,这就有必要考虑近似计算的方法.在定积分的很多应用问题中,被积函数甚至没有解析表达式,可能只是一条实验记录曲线,或者是一组离散的采样值,这时只能应用近似方法去计算相应的定积分.

本实验将主要研究定积分的三种近似计算算法:矩形法、梯形法、抛物线法.对于定积分的近似数值计算,Matlab有专门函数可用.

二、相关函数(命令)及简介

1.sum(a):求数组a的和.

2.format long:长格式,即屏幕显示15位有效数字.

(注:由于本实验要比较近似解法和精确求解间的误差,需要更高的精度).3.double():若输入的是字符则转化为相应的ASCII码;若输入的是整型数值则转化为相应的实型数值.

4.quad():抛物线法求数值积分.

格式: quad(fun,a,b) ,注意此处的fun是函数,并且为数值形式的,所以使用*、/、^等运算时要在其前加上小数点,即 .*、./、.^等.例:Q = quad('1./(x.^3-2*x-5)',0,2);

5.trapz():梯形法求数值积分.

格式:trapz(x,y)

其中x为带有步长的积分区间;y为数值形式的运算(相当于上面介绍的函数fun)

例:计算

0sin()d

x x

π

x=0:pi/100:pi;y=sin(x);

trapz(x,y)

6.dblquad():抛物线法求二重数值积分.

格式:dblquad(fun,xmin,xmax,ymin,ymax),fun可以用inline定义,也可以通过某个函数文件的句柄传递.

例1:Q1 = dblquad(inline('y*sin(x)'), pi, 2*pi, 0, pi)

顺便计算下面的Q2,通过计算,比较Q1 与Q2结果(或加上手工验算),找出积分变量x、y的上下限的函数代入方法.

Q2 = dblquad(inline('y*sin(x)'), 0, pi, pi, 2*pi)

例2:Q3 = dblquad(@integrnd, pi, 2*pi, 0, pi)

这时必须存在一个函数文件:

53

matlab定积分矩形法实验,MATLAB实验三 定积分的近似计算相关推荐

  1. matlab定积分矩形法实验报告,定积分的近似计算(数学实验报告matlab版).doc

    数学实验报告 实验序号:2日期:2013年 12月 5日 班级 2011应数一班姓名孙婉婉学号1101114143 实验名称 定积分的近似计算加深理解积分理论中分割.近似.求和.取极限的思想方法了解定 ...

  2. matlab用矩形法算积分,编写程序,使用矩形法、梯形法、辛普生法计算一重积分...

    function ss=exam3_3(t,a,b,n) %format long %exam3_3(1,0,1,200) %exam3_3(2,0,1,200) %exam3_3(3,0,1,200 ...

  3. 如何用matlab做判别分析法,判别分析matlab流程

    matlab的判别分析_信息与通信_工程科技_专业资料.广西某锰矿床已知两种不同... 在判别分析中称这 m 个变量为判别因子.现有一个新的样本点 x = ( x1 , ? , xm )T 要 判断此 ...

  4. c语言中求定积分矩形法,C语言使用矩形法求定积分的通用函数

    本文为大家分享了C语言使用矩形法求定积分的通用函数,供大家参考,具体内容如下 要求: 写一个用矩形法求定积分的通用函数,分别求:sin(x),cos(x),e^x . 分析: 矩形法,学过高等数学就知 ...

  5. 有限差分方法的matlab编程,有限差分方法的MATLAB编程

    有限差分方法的MATLAB编程 实践中应用的电磁场,其场域的边界大多数比较复杂,使解析法难以应用.虽然有些电磁问题经简化后能解析求解,但解的 (本文共6页) 阅读全文>> 以二叉树方法和有 ...

  6. matlab大津法函数,matlab 大津阈值分割【相关词_ 大津法阈值分割matlab】

    阈值法 阈值分割程序 Otsu Thresholding 赞(0) 踩(0) 收藏(0) 说明:大津法实现图像阈值分割的matlab源程序代码 (Otsu Thresholding Image sou ...

  7. matlab重叠保留法,【matlab实现】重叠相加法与重叠保留法

    两个有限长的序列,当一个序列的长度远大于另外一个序列的时候,如果仍采用常规的处理方法,效率会很低,所以引入了重叠相加法和重叠保留法.下面我们就利用MATLAB来实现这两种方法. 在下面的代码中,我使用 ...

  8. matlab 赋权法,基于MATLAB的非线性曲线拟合赋权法

    成都航空职业技术学院学报 Journal of Chengdu Aeronautic Polytechnic 2018 年 12 月第 4 期(总第 117 期) Vol.34 No.4(Serial ...

  9. matlab的面积法,基于MATLAB的面积计算方法

    胥泽银 XU Ze- yin 第26 卷 第 2 期 2004 年 5 月 物探化探计算技术 Vol126 No. 2 M ay 2004 COM PU T IN G TECHN IQU ES FOR ...

  10. 有限差分方法的matlab编程,有限差分方法的MATLAB编程(1)

    第22卷第4期 <新疆师范大学学报>(自然科学版)JournalofXinjiangNormalUniversity (NaturalSciencesEdition) VoI.22,No. ...

最新文章

  1. 反编译使用yield关键字的方法
  2. max flow value 是网络流里的什么_为什么你这么努力,还是没有通过投行面试
  3. 典型问题分析(十五)
  4. gitee使用svn_Gitee SVN支持
  5. [转]Linux配置防火墙
  6. nanodet学习笔记 tensorrt
  7. FESCAR:阿里重磅开源分布式事务解决方案
  8. 洛谷P1852:跳跳棋(LCA,树形结构)
  9. 数据库表连接总结:等值连接, 自然连接,左外连接,右外连接,内连接,全外连接;
  10. 20111228_给窗体设置默认输入法
  11. wire 和 reg_Wire和Launchkit开源,水流监控系统以及更多新闻
  12. 笨办法学 Python · 续 练习 33:解析器
  13. 牛客网暑期ACM多校训练营(第三场) E Sort String 哈希处理字符串(模板)
  14. Android粒子爆炸特效[可用于任意控件]
  15. 使用jQuery Uploadify在ASP.NET 上传附件
  16. class在java中_class在java里是什么意思
  17. 企业微信扫码登录流程
  18. 印象笔记如何与微信连接到服务器,如何保存微信和微博到印象笔记?
  19. 惠斯通电桥你知道怎么用吗?
  20. win7系统电脑蓝屏怎么解决,如何解决win7电脑蓝屏

热门文章

  1. 学生管理系统 c语言版
  2. 如何写一个NB的商业计划书
  3. Java、前端页面中文乱码解决方式
  4. 【MFC开发(15)】进度条控件Progress Control
  5. 蓝牙耳机连接电脑,提示无法安装驱动程序
  6. 小米/VIVO/OPPO全系列救砖+解锁+工具+教程+激活账户技术
  7. java图书管理系统开发环境_应用Java开发图书管理系统
  8. Qi v1.2.4协议 之 9 Stand-by Power 【英文翻译】
  9. svn删除文件出错的经验总结
  10. 哈里斯鹰优化算法(HHO)附代码