版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、数学实验报告实验序号:2日期:2013年 12月 5日班级2011应数一班姓名孙婉婉学号1101114143实验名称定积分的近似计算问题背景描述: 牛顿莱布尼兹公式仅使用于被积函数的原函数能用初等函数表达出来的情形。如果这点办不到或者不容易办到,就有必要考虑近似计算的方法。在定积分的许多应用问题中,被积函数甚至没有解析表达式,可能只是一条实验记录曲线,或者是一组离散的采样值,这时应用近似方法去计算相应的定积分,以便更好地熟练地掌握定积分的近似计算. 实验目的:1 加深理解积分理论中分割、近似、求和、取极限的思想方法;2 了解定积分近似计算的矩形法、梯形法与抛物线法;3 会用MATLAB语言编写
2、求定积分近似值的程序,会用MALAB中的命令求定积分。实验内容1 分别用梯形法与抛物线法,计算,取n=120.要求使用函数trapz( )、quad( )进行计算求解,并比较结果的差异;2 试计算定积分.(注意:可以运用trapz( )、quad( )或附录程序求解吗?为什么?);3 学习fuluBsum.m的程序设计方法,尝试用函数sum改写附录C的程序,避免for循环。实验原理与数学模型:实验原理与数学模型: 1 矩形法:根据定积分的定义,每一个积分和都可以看作是定积分的一个近似值,即在几何意义上,这是用一系列小矩形面积近似小曲边梯形的结果,所以把这个近似计算方法称为矩形法不过
3、,只有当积分区间被分割得很细时,矩形法才有一定的精确度针对不同的取法,计算结果会有不同。(1) 左点法:对等分区间,在区间上取左端点,即取。(2)右点法:同(1)中划分区间,在区间上取右端点,即取。(3)中点法:同(1)中划分区间,在区间上取中点,即取。2 梯形法等分区间,相应函数值为 ()曲线上相应的点为 ()将曲线的每一段弧用过点,的弦(线性函数)来代替,这使得每个上的曲边梯形成为真正的梯形,其面积为,于是各个小梯形面积之和就是曲边梯形面积的近似值,即 ,称此式为梯形公式。 3 抛物线法将积分区间作等分,分点依次为,对应函数值为(),曲线上相应点为()现把
4、区间上的曲线段用通过三点,的抛物线来近似代替,然后求函数从到的定积分:由于,代入上式整理后得同样也有将这个积分相加即得原来所要计算的定积分的近似值:,即这就是抛物线法公式,也称为辛卜生(Simpson)公式 实验所用软件及版本:2012B主要内容: 1,分别用梯形法与抛物线法,计算,将积分区间1,2作120等分。并尝试用函数trapz(),quad()进行算求解,比较结果的差异。2,试计算定积分.(注意:可以运用trapz()、quad()、或附录程序求解吗?为什么?) 3,学习fuluBsum.m的程序设计方法,尝试用函数sum改写矩形法和抛物线法的程序,避免for循环。实验过程
5、记录(含基本步骤、主要程序清单及异常情况记录等):1、(1)梯形法:、format longn=120;a=1;b=2;inum=0;syms x fxfx=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);endinumx=1:1/120:2;y=1./x;trapz(x,y)(2)抛物线法: format longn=120;a=1;b=2;inum=0;sym
6、s x fxfx=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);endinum quad('1./x',1,2)2、(1)符号求积分:int('sin(x)/x','x',0,inf)(2)quad('sin(x
7、)./x',0,inf)3 ex2zhi2tixingfa.m format long n=120;a=1;b=2;inum=0;syms x fxfx=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 integrate=int(fx,1,2)integrate=double(integrate)fprintf('The rela
8、tive error between inum and real-value is about:%gnn',.abs(inum-integrate)/integrate)ex2zhi2paowuxianfa.mformat long n=120;a=1;b=2;m=0;p=0;syms x fxfx=1./x;for k=1:n i=2*k-1; j=2*k; xi=a+(b-a)*i/(2*n); xj=a+(b-a)*j/(2*n); fxi=subs(fx,'x',xi); fxj=subs(fx,'x',xj); m=m+fxi;p=p+fxj;
9、end fxa=subs(fx,'x',a); fxb=subs(fx,'x',b); inum=(b-a)*(fxa+fxb+4*m+2*p)/(6*n)integrate=int(fx,1,2)integrate=double(integrate)fprintf('The relative error between inum and real-value is about:%gnn',.abs(inum-integrate)/integrate)ex2zhi2paowuxianfa.mformat long n=120;a=1;b=2;m=
10、0;p=0;syms x fxfx=1./x;for k=1:n i=2*k-1; j=2*k; xi=a+(b-a)*i/(2*n); xj=a+(b-a)*j/(2*n); fxi=subs(fx,'x',xi); fxj=subs(fx,'x',xj); m=m+fxi;p=p+fxj;end fxa=subs(fx,'x',a); fxb=subs(fx,'x',b); inum=(b-a)*(fxa+fxb+4*m+2*p)/(6*n)integrate=int(fx,1,2)integrate=double(integ
11、rate)fprintf('The relative error between inum and real-value is about:%gnn',.abs(inum-integrate)/integrate)改写矩形法ex6gaixiejuxingfa.mformat long n=100;a=0;b=1; syms x fxfx=1/(1+x2);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);fxij=(fxi+fxj)/2;
12、m=fxj*(b-a)/n;p=fxi*(b-a)/n;k=fxij*(b-a)/n;inum1=sum(m)inum2=sum(p)inum3=sum(k)改写抛物线法 ex6gaixiepaowuxianfa.mformat long n=100;a=0;b=1; syms x fxfx=1/(1+x2);i=0:(n-1);xj=a+(2*i)*(b-a)/(2*n);xi=a+(2*i+1)*(b-a)/(2*n);xk=a+(2*i+2)*(b-a)/(2*n);fxj=subs(fx,'x',xj);fxi=subs(fx,'x',xi);fxk=s
13、ubs(fx,'x',xk);m=(fxj+4*fxi+fxk)*(b-a)/(6*n);inum=sum(m)实验结果报告与实验总结:run('E:matlabex2zhi2tixingfa.m')inum = 0.69315152080005integrate =log(2)integrate =0.69314718055995The relative error between inum and real-value is about:6.26164e-006>> 0.69315152080005integrate =log(2)integra
14、te =0.69314718055995The relative error between inum and real-value is about:6.26164e-006run('E:matlabex2zhi2paowuxianfa.m')inum = 0.69453606945825integrate =log(2)integrate = 0.69314718055995The relative error between inum and real-value is about:0.00200374run('E:matlabex2zhi2trapz.m')inum = 0.69315152080005integrate =log(2) integrate = 0.69314718055995The relative error between inum and real-value is about:6.26164e-006、(1)梯形法 inum = 0.693151520800048 ans = 0.693151520800048(2)抛物线法 inum = 0.693147180569364 ans = 0.6931471998629702、(1)ans =pi/21, 求定积分时用
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 部编本二年级上册语文第二至七单元(内容含课文口语交际及语文园地)全部教案
- 城市规划学徒指导手册
- 油气勘探钻探施工合同
- 劳务派遣员工健康检查
- 汽车制造锅炉房施工合同
- 环保项目严禁参与虚假环保承诺
- 硫酸厂宿舍楼施工协议
- 科技园区研发创新车库改造协议
- 石油公司出纳人员聘用合同
- 室内运动场地坪施工协议
- 平凡的世界英文简介ppt
- 海康威视视频车位诱导与反向寻车系统与解决与方案
- 四年级湘版美术知识点(复习提纲)
- 牙、牙合、颌的生长发育 颅面的生长发育
- 幕墙预埋件工程专项施工方案
- 公文管理中的错误
- 2023年11月北京地区成人本科学士学位英语真题及答案
- 2020年城市燃气服务企业组织结构及部门职责
- 《计算机操作系统》汤小丹
- GB/T 19682-2005翻译服务译文质量要求
- GB/T 12755-1991建筑用压型钢板
评论
0/150
提交评论