《南通大学实验报告》由会员分享，可在线阅读，更多相关《南通大学实验报告(21页珍藏版)》请在人人文库网上搜索。

1、南通大学实验报告定积分与定积分的近似计算学院： 理学院 班级： 数师153班 学号： 1502012072 姓名： 顾阳 第一部分实验报告书解读一、实验目的实验主要是分析用矩阵公式，梯形公式，辛普森公式求定积分的近似值，并比较它们与定积分的近似情况。可以先学习定积分的数值计算方法，理解定积分的定义，掌握牛顿-莱布尼茨公式。二、实验材料1.1定积分的数值计算计算定积分的近似值，可将积分区间等分而得矩形公式程序为或也可用梯形公式近似计算如果要准确些，可用辛普森公式对于，矩形公式、梯形公式、辛普森公式的Mathematica程序为a=0;b=1;k=10;fx_:=Sinx;d=NIntegrate。

2、fx,x,a,b,k;(计算精确值)s1m_:=NSumfa+i*(b-a)/m*(b-a)/m,i,0,m-1,k;(取小区间左端点的矩形公式)s2m_:=NSumfa+(i+1/2)*(b-a)/m*(b-a)/m,i,0,m-1,k(取小区间中点的矩形公式)s3m_:=NSumfa+i*(b-a)/m*(b-a)/m,i,1,m,k; (取小区间右端点的矩形公式)s4m_:=NSum(fa+i*(b-a)/m+fa+(i+1)*(b-a)/m)/2*(b- a)/m,i,0,m-1,k; (梯形公式) s5m_:=N(b-a)/m/6*(fa+fb)+2*Sumfa+i*(b-a)/m,。

3、i,1,m-1 +4*Sumfa+(i-1/2)*(b-a)/m,i,1,m),k;(辛普森公式) t=Tables1m,r1m,s2m,r2m,s3m,r3m,s4m,r4m,s5m,r5m, m,100,1000,1001.2可积的条件设f(x)=sinx,取a=0,b=1对于，矩形公式、梯形公式、辛普森公式的Mathematica程序为a=0;b=1;k=10;fx_:=Sinx;d=NIntegratefx,x,a,b,k;(计算精确值)s1m_:=NSumfa+i*(b-a)/m*(b-a)/m,i,0,m-1,k;(取小区间左端点的矩形公式)s2m_:=NSumfa+(i+1/2)。

4、*(b-a)/m*(b-a)/m,i,0,m-1,k(取小区间中点的矩形公式)s3m_:=NSumfa+i*(b-a)/m*(b-a)/m,i,1,m,k; (取小区间右端点的矩形公式)s4m_:=NSum(fa+i*(b-a)/m+fa+(i+1)*(b-a)/m)/2*(b- a)/m,i,0,m-1,k; (梯形公式) s5m_:=N(b-a)/m/6*(fa+fb)+2*Sumfa+i*(b-a)/m,i,1,m-1 +4*Sumfa+(i-1/2)*(b-a)/m,i,1,m),k;(辛普森公式) r1m_:=d-s1m;r2m_:=d-s2m;r3m_:=d-s3m;r4m_:= 。

5、t=Tables1m,r1m,s2m,r2m,s3m,r3m,s4m,r4m,s5m,r5m, m,100,1000,1001.3牛顿-莱布尼茨公式设函数在上连续，而且是的一个原函数，则有牛顿-莱布尼兹公式。函数在不连续、不存在原函数，但在上可积；函数在不连续，但在上可积。此外函数处处不连续、不存在原函数，在任意区间(长度大于0)上不可积。求原函数并验证牛顿-莱布尼兹公式的Mathematica程序fx_:=Sinx;Integratef(x),x(求不定积分)Fx_:=%(定义原函数)d=NIntegratef(x),x,a,b(求定积分)df=Fb-Fa (计算原函数的增量)三、实验所用软。

6、件及版本Mathematica5.0第二部分 实验计划(一)定积分的数值计算1.程序修改a=0;b=1;k=10;fx_:=Sinx;d=NIntegratefx,x,a,b,k;s1m_:=NSumfa+i*(b-a)/m*(b-a)/m,i,0,m-1,k; s2m_:=NSumfa+(i+1/2)*(b-a)/m*(b-a)/m,i,0,m-1,k s3m_:=NSumfa+i*(b-a)/m*(b-a)/m,i,1,m,k; s4m_:=NSum(fa+i*(b-a)/m+fa+(i+1)*(b-a)/m)/2*(b- a)/m,i,0,m-1,k; s5m_:=N(b-a)/m/6*。

7、(fa+fb)+2*Sumfa+i*(b-a)/m,i,1,m-1 +4*Sumfa+(i-1/2)*(b-a)/m,i,1,m),k; r1m_:=d-s1m;r2m_:=d-s2m;r3m_:=d-s3m;r4m_:=d-s4m;r5m_:=d-s5m; t=Tables1m,r1m,s2m,r2m,s3m,r3m,s4m,r4m,s5m,r5m, m,100,1000,100利用以上程序计算，并对几个公式比较。2.实验思路对以上程序，分别将sinx的x替换成1，x，x2，ex，In(1+x)(二)可积的条件1.实验思路：(1)如果函数f(x)在区间a,b上连续，则f(x)在区间a,b上可。

8、积，反之亦然。(2)设一连续函数，判断其是否可积。2.程序修改a=0;b=1;k=10;fx_:=Sinx;d=NIntegratefx,x,a,b,k; s1m_:=NSumfa+i*(b-a)/m*(b-a)/m,i,0,m-1,k; s2m_:=NSumfa+(i+1/2)*(b-a)/m*(b-a)/m,i,0,m-1,k s3m_:=NSumfa+i*(b-a)/m*(b-a)/m,i,1,m,k; s4m_:=NSum(fa+i*(b-a)/m+fa+(i+1)*(b-a)/m)/2*(b- a)/m,i,0,m-1,k; s5m_:=N(b-a)/m/6*(fa+fb)+2*Su。

9、mfa+i*(b-a)/m,i,1,m-1 +4*Sumfa+(i-1/2)*(b-a)/m,i,1,m),k; r1m_:=d-s1m;r2m_:=d-s2m;r3m_:=d-s3m;r4m_:=d-s4m;r5m_:=d-s5m; t=Tables1m,r1m,s2m,r2m,s3m,r3m,s4m,r4m,s5m,r5m, m,100,1000,100利用以上程序计算，并对几个公式比较。(三)牛顿-莱布尼茨公式1. 程序修改fx_:=Sinx;Integratef(x),x Fx_:=%d=NIntegratef(x),x,a,bdf=Fb-Fa r=d-df2实验思路(1)先对一个函数。

10、sinx在区间0,1时,运行程序计算。(2)在考虑其他函数，y=1,y=x，y=x2，y=ex，y=In(1+x)，y=sign(x)在0,1时，进行程序计算。第三部分 实验过程与结果实验一 定积分的实验计算1. 在mathmatica上输入以下程序a=0;b=1;k=10;fx_:=Sinx;d=NIntegratefx,x,a,b,k; s1m_:=NSumfa+i*(b-a)/m*(b-a)/m,i,0,m-1,k; s2m_:=NSumfa+(i+1/2)*(b-a)/m*(b-a)/m,i,0,m-1,k s3m_:=NSumfa+i*(b-a)/m*(b-a)/m,i,1,m,k;。

11、 s4m_:=NSum(fa+i*(b-a)/m+fa+(i+1)*(b-a)/m)/2*(b- a)/m,i,0,m-1,k; s5m_:=N(b-a)/m/6*(fa+fb)+2*Sumfa+i*(b-a)/m,i,1,m-1 +4*Sumfa+(i-1/2)*(b-a)/m,i,1,m),k; r1m_:=d-s1m;r2m_:=d-s2m;r3m_:=d-s3m;r4m_:=d-s4m;r5m_:=d-s5m; t=Tables1m,r1m,s2m,r2m,s3m,r3m,s4m,r4m,s5m,r5m, m,100,1000,100运行结果为：2. 在mathmatica上输入以下程。

12、序a=0;b=1;k=10;fx_:=1d=NIntegratefx,x,a,b,k; s1m_:=NSumfa+i*(b-a)/m*(b-a)/m,i,0,m-1,k; s2m_:=NSumfa+(i+1/2)*(b-a)/m*(b-a)/m,i,0,m-1,k s3m_:=NSumfa+i*(b-a)/m*(b-a)/m,i,1,m,k; s4m_:=NSum(fa+i*(b-a)/m+fa+(i+1)*(b-a)/m)/2*(b- a)/m,i,0,m-1,k; s5m_:=N(b-a)/m/6*(fa+fb)+2*Sumfa+i*(b-a)/m,i,1,m-1 +4*Sumfa+(i-。

13、1/2)*(b-a)/m,i,1,m),k; r1m_:=d-s1m;r2m_:=d-s2m;r3m_:=d-s3m;r4m_:=d-s4m;r5m_:=d-s5m; t=Tables1m,r1m,s2m,r2m,s3m,r3m,s4m,r4m,s5m,r5m, m,100,1000,100运行结果为：3. 在mathmatica上输入以下程序a=0;b=1;k=10;fx_:=xd=NIntegratefx,x,a,b,k; s1m_:=NSumfa+i*(b-a)/m*(b-a)/m,i,0,m-1,k; s2m_:=NSumfa+(i+1/2)*(b-a)/m*(b-a)/m,i,0,m。

14、-1,k s3m_:=NSumfa+i*(b-a)/m*(b-a)/m,i,1,m,k; s4m_:=NSum(fa+i*(b-a)/m+fa+(i+1)*(b-a)/m)/2*(b- a)/m,i,0,m-1,k; s5m_:=N(b-a)/m/6*(fa+fb)+2*Sumfa+i*(b-a)/m,i,1,m-1 +4*Sumfa+(i-1/2)*(b-a)/m,i,1,m),k; r1m_:=d-s1m;r2m_:=d-s2m;r3m_:=d-s3m;r4m_:=d-s4m;r5m_:=d-s5m; t=Tables1m,r1m,s2m,r2m,s3m,r3m,s4m,r4m,s5m,r。

15、5m, m,100,1000,100运行结果为：4. 在mathmatica上输入以下程序a=0;b=1;k=10;fx_:=x2d=NIntegratefx,x,a,b,k; s1m_:=NSumfa+i*(b-a)/m*(b-a)/m,i,0,m-1,k; s2m_:=NSumfa+(i+1/2)*(b-a)/m*(b-a)/m,i,0,m-1,k s3m_:=NSumfa+i*(b-a)/m*(b-a)/m,i,1,m,k; s4m_:=NSum(fa+i*(b-a)/m+fa+(i+1)*(b-a)/m)/2*(b-a)/m,i,0,m-1,k; s5m_:=N(b-a)/m/6*(f。

16、a+fb)+2*Sumfa+i*(b-a)/m,i,1,m-1 +4*Sumfa+(i-1/2)*(b-a)/m,i,1,m),k; r1m_:=d-s1m;r2m_:=d-s2m;r3m_:=d-s3m;r4m_:=d-s4m;r5m_:=d-s5m; t=Tables1m,r1m,s2m,r2m,s3m,r3m,s4m,r4m,s5m,r5m, m,100,1000,100运行结果为：5. 在mathmatica上输入以下程序a=0;b=1;k=10;fx_:=Expxd=NIntegratefx,x,a,b,k; s1m_:=NSumfa+i*(b-a)/m*(b-a)/m,i,0,m-。

17、1,k; s2m_:=NSumfa+(i+1/2)*(b-a)/m*(b-a)/m,i,0,m-1,k s3m_:=NSumfa+i*(b-a)/m*(b-a)/m,i,1,m,k; s4m_:=NSum(fa+i*(b-a)/m+fa+(i+1)*(b-a)/m)/2*(b- a)/m,i,0,m-1,k; s5m_:=N(b-a)/m/6*(fa+fb)+2*Sumfa+i*(b-a)/m,i,1,m-1 +4*Sumfa+(i-1/2)*(b-a)/m,i,1,m),k; r1m_:=d-s1m;r2m_:=d-s2m;r3m_:=d-s3m;r4m_:=d-s4m;r5m_:=d-s5。

18、m; t=Tables1m,r1m,s2m,r2m,s3m,r3m,s4m,r4m,s5m,r5m, m,100,1000,100运行结果为：6. 在mathmatica上输入以下程序a=0;b=1;k=10;fx_:= Log1+xd=NIntegratefx,x,a,b,k; s1m_:=NSumfa+i*(b-a)/m*(b-a)/m,i,0,m-1,k; s2m_:=NSumfa+(i+1/2)*(b-a)/m*(b-a)/m,i,0,m-1,k s3m_:=NSumfa+i*(b-a)/m*(b-a)/m,i,1,m,k; s4m_:=NSum(fa+i*(b-a)/m+fa+(i+。

19、1)*(b-a)/m)/2*(b- a)/m,i,0,m-1,k; s5m_:=N(b-a)/m/6*(fa+fb)+2*Sumfa+i*(b-a)/m,i,1,m-1 +4*Sumfa+(i-1/2)*(b-a)/m,i,1,m),k; r1m_:=d-s1m;r2m_:=d-s2m;r3m_:=d-s3m;r4m_:=d-s4m;r5m_:=d-s5m; t=Tables1m,r1m,s2m,r2m,s3m,r3m,s4m,r4m,s5m,r5m, m,100,1000,100运行结果为：7. 实验观察结果：随着n的增大，以及公式的更加精确性，求积分的误差d越来越小，结果越来越准确。实验二。

20、 可积条件1. 以实验一的6组数据为基础，再加一组数据；2. 在mathmatica上输入以下程序a=0;b=1;k=10;fx_:=Signxd=NIntegratefx,x,a,b,k; s1m_:=NSumfa+i*(b-a)/m*(b-a)/m,i,0,m-1,k; s2m_:=NSumfa+(i+1/2)*(b-a)/m*(b-a)/m,i,0,m-1,k s3m_:=NSumfa+i*(b-a)/m*(b-a)/m,i,1,m,k; s4m_:=NSum(fa+i*(b-a)/m+fa+(i+1)*(b-a)/m)/2*(b- a)/m,i,0,m-1,k; s5m_:=N(b-a。

21、)/m/6*(fa+fb)+2*Sumfa+i*(b-a)/m,i,1,m-1 +4*Sumfa+(i-1/2)*(b-a)/m,i,1,m),k; r1m_:=d-s1m;r2m_:=d-s2m;r3m_:=d-s3m;r4m_:=d-s4m;r5m_:=d-s5m; t=Tables1m,r1m,s2m,r2m,s3m,r3m,s4m,r4m,s5m,r5m, m,100,1000,100运行结果为：3. 实验观察结果：连续函数一定可积，但是非连续函数不一定不可积实验三 牛顿-莱布尼茨公式1. 在mathmatica上输入以下程序a=0;b=1;fx_:=x;Integratefx,xFx。

22、_=%d=NIntegratefx,x,a,bdf=Fb-Far=d-df运行结果为2. 在mathmatica上输入以下程序a=0;b=1fx_:=1;Integratefx,xFx_:=%d=NIntegratefx,x,a,bdf=Fb-Far=d-df运行结果为3. 在mathmatica上输入以下程序a=0;b=1fx_:=x;Integratefx,xFx_:=%a=0;b=1d=NIntegratefx,x,a,bdf=Fb-Far=d-df运行结果为4. 在mathmatica上输入以下程序a=0;b=1fx_:=x2;Integratefx,xFx_:=%d=NIntegra。

23、tefx,x,a,bdf=Fb-Far=d-df运行结果为5. 在mathmatica上输入以下程序a=0;b=1fx_:=Expx;Integratefx,xFx_:=%d=NIntegratefx,x,a,bdf=Fb-Far=d-df运行结果为6. 在mathmatica上输入以下程序a=0;b=1fx_:=Log1+x;Integratefx,xFx_:=%d=NIntegratefx,x,a,bdf=Fb-Far=d-df运行结果为7. 在mathmatica上输入以下程序a=0;b=1fx_:=SignxIntegratefx,xFx_:=%d=NIntegratefx,x,a,bdf=Fb-Far=d-df运行结果为8. 实验观察结果： 当f(x)连续时，牛顿-莱布尼茨公式成立第四部分 实验总结1 .将n无限细分时，也可用矩形公式，梯形公式，辛普森公式求定积分的近似值，其中他们的精确度依次增加。2. 如果f(x)是连续函数，则f(x)在区间a,b上可积，反之f(x)在区间a,b可积，但f(x)不一定连续3. 莱布尼茨公式成立，即若函数f(x)在a,b上连续，而且F(x)是f(x)的一个原函数，则有abf(x)dx=F(b)-F(a).21 / 21文档可自由编辑。

南通大学计算机类课程实验报告,南通大学实验报告相关推荐

  1. 计算机课程实验教学案例,联席会|第四届全国计算机类课程实验教学案例设计竞赛|...

    一.竞赛目的 为实现高校计算机类课程实验教学体系与教学内容的协同,推进探究性.创新性实验教学改革,提升高校教师的实验教学水平.以问题为导向,在教学中培育大学生创新意识和工程实践能力,促进高等教育质量工 ...

  2. 中职院校计算机专业课件视频,中职计算机类课程微视频的设计与制作浅谈

    李刚平 [摘 要] 随着网络教学在计算机类实验课程的广泛应用,微视频课程已成为实验类教学改革的研究热点.在如今的数字化"微"时代,微视频课将成为一种新的教学模式和学习方式.阐述计算 ...

  3. 计算机类高职院校课题,高职院校公共计算机类课程体系构建的教改研究

    龙源期刊网 http://www.360docs.net/doc/1b395d5a89eb172ded63b7b4.html 高职院校公共计算机类课程体系构建的教改研究 作者:刘永红 来源:<计 ...

  4. 斯坦福大学计算机类课程视频

    斯坦福大学计算机类课程都是以CS开头编号,可以在网址https://exploredegrees.stanford.edu/coursedescriptions/cs/查询,在网上可以登录查看课程的课 ...

  5. 南通大学计算机专业分数线2020,2020南通大学录取分数线_历年各专业分数线（2017-2019）_各省投档线_一品高考网...

    录取分数线,投档线 南通大学2019录取分数线是多少,各专业录取分数线是多少,是每个... 今年高考志愿填报已经结束,各批次录取工作陆续展开,各省各批次... 序号 省份 专业名称 招生数 省控线(文 ...

  6. 计算机类课程辅助工具设计,计算机辅助设计07075 .doc

    计算机辅助设计07075 湖北省高等教育自学考试课程考试大纲 课程名称:计算机辅助设计 课程代号:07075 第一部分 课程性质与目标 一.课程性质和特点 计算机辅助设计课程是高等教育自学考试计算机类 ...

  7. 2020全国高职计算机类就业率,2020中国大学生就业报告：这几个关键数据要注意...

    解读2020大学生就业报告:填报高考志愿可做重要参考 1.2019届(即2019年毕业的学生)本科毕业生平均月收入为5440元,(对比全国6亿人月薪1000元以下),课件读大学是农村子弟跳出阶层.帮助 ...

  8. 关于计算机类课程实验教学的思考

    由于计算机机房的限制,没有安装投影仪以及计算机控制管理软件,所以单纯的讲解效果非常不好.老师在上面讲,学生在下面干什么根本不知道.由于教授的实验课和理论课关系不是非常密切.试验的主要内容是网页设计,也 ...

  9. 计算机类课程嵌入式系统的特点及其应用,嵌入式系统的准确定义、特点及其重要性...

    嵌入式计算机系统的出现,是现代计算机发展史上的里程碑事件.嵌入式系统诞生于微型计算机时代,与通用计算机的发展道路完全不同,形成了独立的单芯片的技术发展道路.由于嵌入式系统的诞生,现代计算机领域中出现了 ...

最新文章

  1. device:GPU:0 but available devices are [ /job:localhost/replica:0/task:0/dev
  2. bread是可数还是不可数_雅思官方语法教程之——这个名词到底可不可数？
  3. sublime text3 使用SVN插件
  4. go语言代码规范详解
  5. 天津理工大学计算机工程学院院长,王法玉
  6. 深度学习面试的一些知识
  7. 路由技术——OSPF
  8. 使用Json.Net处理json序列化和反序列化接口或继承类
  9. windows，linux桌面系统管理
  10. c语言编程学习宝典,C语言学习宝典app
  11. win10专业版虚拟机配置服务器,win10专业版怎么运行虚拟机_win10专业版开启虚拟机的方法...
  12. php mysql 连接池_php 如何实现 数据库 连接池
  13. 计算机上的策略禁止用户安装,Win7系统如何使用组策略禁止安装软件？
  14. ajax发送put请求参数,请求参数和PUT方法(Request params and PUT method)
  15. 数字逻辑——时序逻辑电路
  16. 亚控科技工作中的编程知识小积累
  17. Julia教程：Julia语言入门
  18. Hunter’s Apprentice （猎人的学徒）——【Green 公式（ 判断多边形边界曲线顺/逆时针】
  19. prettier简单介绍
  20. SML中的local

热门文章

  1. 财务机器人实现高效赋能财务精细化运营
  2. cbac式_AABC,AAB,AABC,AABB,ABAB,ABAC,ABCA式的词语大全
  3. B1/gsm功放板下行测试流程
  4. typora+picgo+gitee一键搬家csdn
  5. 【OpenMLDB Meetup #1】会议纪要
  6. 关于python语句p p_关于Python语句P=-P,哪个选项的描述是正确的?_学小易找答案
  7. ESP8266呼吸灯亮度调节并且实时显示PWM数值
  8. WEB安全_基础入门_数据包拓展，安装Burp Suite，练习来源页伪造、浏览器信息伪造、HTTP动作练习、投票系统程序设计缺陷分析
  9. 利用PS把多张psd格式的图片转换为一张PDF格式
  10. 深划痕需要大面积补漆吗_什么程度的划痕才需补漆? 小心亏大了!