- 1、本文档共8页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
定积分的运用
实验四 定积分及应用
实验的目的
1、掌握利用Matlab进行积分运算;
2、掌握积分在计算面积、体积等问题中的应用;
3、掌握各种积分指令的区别与特点。
实验的基本理论与方法
1、定积分定义:
函数在区间上的定积分定义为:设函数在上有界,在区间上任取个分点:,把分成个小区间, 。这些分点构成对区间的一个分割,用T表示。小区间的长度为。记,称为分割T的模。在区间上取点,做函数值与小区间长度的乘积,并作和。当时,和总趋于确定的极限,这时这个极限为函数在区间上的定积分,记作。即
。
2、定积分的应用
①计算平面图形的面积:由连续曲线,直线及x轴所围成的曲边梯形面积为;
②计算旋转体的体积:由连续曲线,直线及x轴所围成的曲边梯形绕x轴旋转一周所成立体的体积为;
③计算平面曲线的弧长:设曲线弧由直线坐标方程给出,其中在上具有一阶连续导数,则曲线弧长;
设曲线弧由参数方程给出,其中在上具有连续导数,则曲线弧长;
设曲线弧由极坐标方程给出,其中在上具有连续导数,则曲线弧长。
④定积分的近似计算:
用分点将区间分成个长度相等的小区间,每个小区间长度为。
矩形法公式:
梯形法公式:
抛物线法公式:
实验使用的函数与命令
1、符号积分命令int
int(fun):求函数fun的不定积分;
int(fun,var):求函数fun关于变量var的不定积分;
int(fun, var, a,b,):求函数fun的在[a,b]间的定积分或广义积分;
2、数值计算定积分quad, quad1,trapz
①quad(fun,a,b,tol):fun为被积的函数名,a, b微积分上下限,tol为精度,若缺省,其缺省值为1.0e-6;
②quad(fun,a,b,tol,trace):参数fun,a,b,tol用法与上面相同,而输入第五个非零参数trace,是对积分过程通过被积函数上的图像进行跟踪。
对于quad使用自适应步长Simpson法,而quadl的调用格式与quad一致,但它使用Lobbato算法,其精度和速度要比quad的精度高。
③Z = trapz(X,Y):用于进行梯形积分,精度低,适用于数值函数和光滑性不好的函数。 其中X表示积分区间的离散化变量,Y与X同维的向量,表示被积函数,Z返回积分的近似值。
quad, quad1,trapz都不能用于广义积分,此外由于数值积分的特点,对一些假奇异积分也不能直接求解。
3、Inline
实验指导
例1 用符号积分命令int计算下列各题
1)不定积分与;2)定积分;3)广义积分
解:1)由于该题比较简单,我们直接在命令窗口输入Matlab代码
clear;syms x y z;
int(sin(x*y+z),z)
得到的计算结果
ans =
-cos(x*y+z)
再继续输入
int(x^2*sin(x))
得到的计算结果
ans =
-x^2*cos(x)+2*cos(x)+2*x*sin(x)
在以上的结果基础上,可以利用diff命令验证计算结果的正确性。
2)在命令窗口输入Matlab代码
syms x;
int(1/(3+2*x+x^2),x,0,1)
得到定积分的值
ans =
1/2*atan(2^(1/2))*2^(1/2)-1/2*atan(1/2*2^(1/2))*2^(1/2)
3)在2)的基础上继续输入
int(1/(3+2*x+x^2),x,-inf,inf)
得到广义积分的值
ans =
1/2*pi*2^(1/2)
例2 计算数值积分。
解:我们先用trapz(X,Y)来求,在命令窗口输入Matlab代码
x=-1:0.1:2; y=x.^3;%取积分步长为0.1
trapz(x,y)
输出结果为
ans =
3.7575
我们现在把步长取到0.01
x=-1:0.01:2; y=x.^3;
trapz(x,y)
输出结果为
ans =
3.7501
通过上面的运算结果可以看到,利用不同的步长进行计算,可以看到步长和精度之间的关系。如果取步长为0.001则输出结果为3.7500与步长取到0.01值接近。
可以利用符号积分命令计算
syms x;
int(x^3,x,-1,2)
ans = 15/4
下面我们利用quad,quad1来计算,首先要建立名为jifen.m的文件。
function y=jifen(x)
y=x.^3;
然后进行在命令窗口进行调用并输出结果
quad(jifen,-1,2)
ans =
3.7500:
利用quad1同样得到上面结果。
例3 计算数值积分。并用符号积分指令求解,观察输出结果。
解:利用trapz计算数值积分
x=0:0.1:1; y=sin(x.^2)./(1+x);
tra
您可能关注的文档
- 孙长才-加强应收帐款的管理.doc
- 孟子与企业管理.doc
- 学习《教学教育过程最优化》的详细笔记-巴班斯基.doc
- 学习中医之复习精华.doc
- 孝直小学环境设计方案.doc
- 学分制管理办法(修订).doc
- 学前专业手工基础教程平时作业.doc
- 学前儿童古诗90首.doc
- 学校“创先争优”活动方案.doc
- 学校安全管理员工作职责.doc
- 英语人教PEP版八年级(上册)Unit4+writing+写作.pptx
- 人美版美术四年级(上册)8 笔的世界 课件 (1).pptx
- 人美版美术七年级(上册)龙的制作.pptx
- 英语人教PEP版六年级(上册)Unit 2 第一课时.pptx
- 数学苏教版三年级(上册)3.3 长方形和正方形周长的计算 苏教版(共12张PPT).pptx
- 音乐人教版八年级(上册)青春舞曲 课件2.pptx
- 音乐人教版四年级(上册) 第一单元 音乐知识 附点四分音符|人教版.pptx
- 英语人教PEP版四年级(上册)Unit 6 Part B let's learn 1.pptx
- 道德与法治人教版二年级(上册)课件-3.11大家排好队部编版(共18张PPT).pptx
- 人美版美术七年级(上册)《黄山天下奇》课件1.pptx
文档评论(0)