数值积分梯形算法.doc_第1页
数值积分梯形算法.doc_第2页
数值积分梯形算法.doc_第3页
全文预览已结束

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

数学实验报告实验序号：课件制作（1）日期：2014年6月20日班级姓名学号课件名称数值积分梯形算法问题背景描述：求某函数的定积分时，在多数情况下，被积函数的原函数很难用初等函数表达出来，因此能够借助微积分学的牛顿-莱布尼兹公式计算定积分的机会是不多的。另外，许多实际问题中的被积函数往往是列表函数或其他形式的非连续函数，对这类函数的定积分，也不能用不定积分方法求解。由于以上原因，数值积分的理论与方法一直是计算数学研究的基本课题。构造数值积分公式最通常的方法是用积分区间上的n 次插值多项式代替被积函数，由此导出的求积公式称为插值型求积公式。特别在节点分布等距的情形称为牛顿-柯茨公式，例如梯形公式与抛物线公式就是最基本的近似公式。但是它们的精度较差。而且高阶Newton-Cotes求积公式是不稳定的。因此，通常不用高阶求积公式得到比较精确的积分值，而是将整个积分区间分段，在每一小段上用低阶求积公式。这种方法称为复化求积方法。采用复化梯形求积算法时，将积分区间a,b划分为n等份，步长分点为 ，k=0,1,n。所谓复化求积法，就是先用低阶的Newton-Cotes公式求得每个子区间上的积分值，然后再求和，用 作为所求积分I的近似值公式实验目的： 了解数值积分梯形解法的基本思想，会用Mathematica软件实现完整的过程实验所用软件及版本：Mathematica7以上版本主要内容（要点）：数值积分梯形算法的公式，显示了精确值，近似值，误差及图形显示，其中等分点数用控件实现，很好的实现交互功能实验过程记录在Mathematica环境中输入以下代码： 接着按Shift+Enter就出现课件画面如下图，通过操纵界面上的控件，可得出下面几幅图像，从中分析可得出一些结论。思考与深入： 从交互式画面可以明显看到，复化求积公式在使用时，发现运用复化梯形公式求解积分效果较好，误差较小必须但是从理论上，我们分析，复合梯形公式的误差限的系数为-（b-a）/12，而且后面是h平方级还有f的二次导数，而辛普森的系数是-（b-a）/（180*16），后面h是四次方级的，f的导数为四次导数，显然辛普森的误差限更加小而高斯求积公式是对代数精度的方面有着更加好的结果。此课件还有改进之处，如“选择函数”可进行动态输入，并不只局限于列表中的几个函数。此课件有几个优点：1、代码开放可随时修改；2、运行方便，做数学课件易如反掌；3、可通过存为CDF格式放在网上，成为网络课件。当然，它也有使用上的难处，就是要对Mathematica要比较全面深入地了解，若能把程序设计的思想融入其中，则一定能设计出功能强大的数学
人人文库> 全部分类> 应用文书 > 技术指导

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论