数值分析+复合梯形公式实验报告_第1页
数值分析+复合梯形公式实验报告_第2页
数值分析+复合梯形公式实验报告_第3页
数值分析+复合梯形公式实验报告_第4页
全文预览已结束

下载本文档

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

文档简介

数学与信息工程学院实验报告课程名称：数值分析实验室：实验台号：班级：姓名：实验日期：2012年3月14日实验名称复合梯形公式实验目的和要求1.掌握复合梯形公式的基本思想。

2.编程实现用复合梯形公式求积分。

3.熟悉matlab软件的使用。实验内容和步骤：实验内容：根据matlab编写复合梯形公式求积分代码实验步骤：（1）开启软件平台——MATLAB，编程；（2）根据数值解法步骤编写M文件；①functiont=trapz(f,a,b,n)h=(b-a)/n;sum=0;fork=1:n-1x=a+k*h;sum=sum+feval(f,x);endt=h/2*(feval(f,a)+feval(f,b)+2*sum);②执行M文件trapz(f,a,b,n)（3）观察运行结果实验数据记录：数据1：trapz(inline('sqrt(x)'),0.5,1,2)ans=0.4299数据2：trapz(inline('sqrt(x)'),0.5,1,3)ans=0.4305数据3：trapz(inline('sqrt(x)'),0.5,1,4)ans=0.4307数据4：trapz(inline('sqrt(x)'),0.5,1,100)ans=0.4310数据5：trapz(inline('sqrt(x)'),0.5,1,1000)ans=0.4310实验结果分析：不同的数据代入运行结果不同，n越大数据越准确，n大到一定时，数据不再变化。M文件保存的时候要与函数文件名一致；M文件中f(a)要写成feval(f,a)；执行文件时函数前要用inline，否则运行不了；成绩评定

人人文库> 全部分类> 应用文书 > 作业报告

温馨提示

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

最新文档

评论

0/150

提交评论