matlab定积分矩形法实验报告,定积分的近似计算(数学实验报告matlab版).doc
数学实验报告
实验序号:2日期:2013年 12月 5日
班级
2011应数一班姓名孙婉婉学号1101114143
实验名称
定积分的近似计算加深理解积分理论中分割、近似、求和、取极限的思想方法了解定积分近似计算的矩形法、梯形法与抛物线法会用MATLAB语言编写求定积分近似值的程序,会用MALAB中的命令求定积分。,取n=120.要求使用函数trapz( )、quad( )进行计算求解,并比较结果的差异;
2 试计算定积分.(注意:可以运用trapz( )、quad( )或附录程序求解吗?为什么?);
3 学习fuluBsum.m的程序设计方法,尝试用函数sum改写附录C的程序,避免for循环。实验原理与数学模型:
实验原理与数学模型:
1.? 矩形法:根据定积分的定义,每一个积分和都可以看作是定积分的一个近似值,即在几何意义上,这是用一系列小矩形面积近似小曲边梯形的结果,所以把这个近似计算方法称为矩形法.不过,只有当积分区间被分割得很细时,矩形法才有一定的精确度.
针对不同的取法,计算结果会有不同。
(1) 左点法:对等分区间
,在区间上取左端点,即取。
(2)右点法:同(1)中划分区间,在区间上取右端点,即取。
(3)中点法:同(1)中划分区间,在区间上取中点,即取。
2.? 梯形法
等分区间
,相应函数值为 ().曲线上相应的点为 ()将曲线的每一段弧用过点,的弦(线性函数)
来代替,这使得每个上的曲边梯形成为真正的梯形,其面积为
,.于是各个小梯形面积之和就是曲边梯形面积的近似值,
,
即 ,称此式为梯形公式。?
3.? 抛物线法
将积分区间作等分,分点依次为
,,对应函数值为
(),曲线上相应点为
().
现把区间上的曲线段用通过三点,,的抛物线
来近似代替,然后求函数从到的定积分:
由于,代入上式整理后得
同样也有
……
将这个积分相加即得原来所要计算的定积分的近似值:
,
即
这就是抛物线法公式,也称为辛卜生(Simpson)公式.?
实验所用软件及版本:2012B主要内容:
1,分别用梯形法与抛物线法,计算,将积分区间[1,2]作120等分。并尝试用函数trapz(),quad()进行算求解,比较结果的差异。
2,试计算定积分.
(注意:可以运用trapz()、quad()、或附录程序求解吗?为什么?)
3,学习fuluBsum.m的程序设计方法,尝试用函数sum改写矩形法和抛物线法的程序,避免for循环。
实验过程记录(含基本步骤、主要程序清单及异常情况记录等):
1、(1)梯形法:
①、format long
n=120;a=1;b=2;inum=0;
syms x fx
fx=1/x;
for i=1:n
xj=a+(i-1)*(b-a)/n;
xi=a+i*(b-a)/n;
fxj=subs(fx,'x',xj);
fxi=subs(fx,'x',xi);
inum=inum+(fxj+fxi)*(b-a)/(2*n);
end
inum
②x=1:1/120:2;
y=1./x;
trapz(x,y)
(2)抛物线法:
① format long
n=120;a=1;b=2;inum=0;
syms x fx
fx=1/x;
for i=1:n
xj=a+(i-1)*(b-a)/n;
xi=a+i*(b-a)/n;
xk=(xi+xj)/2;
fxj=subs(fx,'x',xj);
fxi=subs(fx,'x',xi);
fxk=subs(fx,'x',xk);
inum=inum+(fxj+4*fxk+fxi)*(b-a)/(6*n);
end
inum
②quad('1./x',1,2)
2、(1)符号求积分:
int('sin(x)/x','x',0,inf)
(2)quad('sin(x)./x',0,inf)
3 ex2zhi2tixingfa.m
format long
n=120;a=1;b=2;inum=0;
syms x fx
fx=1./x;
for i=1:n
xj=a+(i-1)*(b-a)/n; xi=a+i*(b-a)/n;
fxj=sub
matlab定积分矩形法实验报告,定积分的近似计算(数学实验报告matlab版).doc相关推荐
- c语言中求定积分矩形法,C语言使用矩形法求定积分的通用函数
本文为大家分享了C语言使用矩形法求定积分的通用函数,供大家参考,具体内容如下 要求: 写一个用矩形法求定积分的通用函数,分别求:sin(x),cos(x),e^x . 分析: 矩形法,学过高等数学就知 ...
- 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 ...
- matlab某分子由25个原子组成,清华大学数学实验-实验7无约束优化1
数学实验 实验 7 无约束优化 实验 7 无约束优化 生医 0 王言 2010013212 实验目的 1 掌握用 MATLAB 优化工具箱的基本用法 对不同算法进行初步分析 比较 2 练习用无约束优化 ...
- matlab某分子由25个原子组成,清华大学数学实验实验7无约束优化1
清华大学数学实验实验7无约束优化1 (8页) 本资源提供全文预览,点击全文预览即可全文预览,如果喜欢文档就下载吧,查找使用更方便哦! 11.9 积分 实验7无约束优化生医0 王言 2010013212 ...
- matlab内将数组内所有数累乘,数学实验上机汇总【未完成】
数学实验上机作业整理∈hyd 实验一 1. 计算球体体积(半径r=5) r=5;v=(4/3)*pi*r^3 v =523.5988 2.设矩阵1234567891023416A ?? ? = ? ? ...
- 数学建模还原纸片matlab程序,碎纸片的拼接复原问题大学生数学建模全国一等奖论文.doc...
碎纸片的拼接复原问题 摘要 为解决碎纸片的拼接复原问题,我们通过定义差异度指数.高度差,建立0-1规划模型,使用聚类分析.MATLAB搜索算法和人工干预等相结合,得到了所有附件复原序号和复原图片. 针 ...
- 矩形法求定积分的原理和实现
矩形法求解定积分 做个小笔记,如表述有不对的地方,欢迎指正. 原理 积分在图形上表示就是面积,把被积分区域划分成一个个小矩形,求解这些矩形的面积和即可. 被积函数:f(x)f(x)f(x); 积分区间 ...
- 求矩形的最小值c语言,C语言复习---矩形法求定积分函数
#include #include using namespace std; int main(){ float fun(float x); float a, b; cout << &qu ...
- c 语言 用矩形法求定积分,如何用矩形法(梯形法)求定积分
分析: 高中的时候,我们学习过,可以通过矩形法或者矩形法来求定积分. 思路就是将积分区间划分成n等份,然后将这n等份近似看成矩形(或梯形),然后对所有的矩形(或梯形)的面积进行求和. 简单的例子: 求 ...
最新文章
- torch.log函数
- SAP Spartacus org unit list点击item之后的页面跳转实现
- Android之华为手机打开app奔溃提示java.io.FileNotFoundException: res/drawable/abc_vector_test.xml
- spring jaxb_自定义Spring命名空间使JAXB变得更容易
- 少儿编程100讲轻松学python(二)-python cv2模块怎么安装
- UE4中的字符串转换
- 禅道能连接数据库无法打开页面
- JavaScript学习(五十九)—原型、原型链、闭包以及闭包的不足
- 电线电缆行业MES解决方案
- 关于使用WS2811点灯
- 指数函数为例的超越函数的逼近拟合误差分析
- EasyUI - panel 高度自适应
- MAC系统中快速切换输入法的方法
- 第三课 排序算法LowB三人组
- 查看计算机ip地址配置信息的命令是,查看电脑配置指令大全
- ROS工业机器人和工业自动化竞赛ARIAC 2021即将启动
- 显示屏信号超过输入范围1600x900@60HZ ------解决办法
- dict.get / dict.setdefault / defaultdict
- 神奇的cyclegan!
- with torch.no_grad()
热门文章
- 常用软件简简单单变成绿色版
- php pdo 1064,php-即使有错误,PDO错误代码也总是00000
- 人工智能研究中心快递柜——源码部署及分析综述
- 对于室内定位的认识和看法?
- Windows升级日志文件可以删除吗?
- window.print()打印时出现网页的链接等信息的页眉页脚,去除方法
- 计算机基础知识一级分值,计算机一级考试内容及分值
- Map.putall()方法
- Google卫星地图下的港珠澳大桥
- Android getResources()。getDrawable()已弃用API 22