复合辛普森求积_第1页
复合辛普森求积_第2页
全文预览已结束

下载本文档

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

文档简介

精品文档第三次实验实验名称复合辛普森求积计算积分实验时间2012.05.06姓名班级数应二班学号成绩一、实验目的，内容 1. 目的： 学习并理解复合辛普森求积计算积分的matlab实现。 2. 内容： 用matlab计算积分（ 精确值-0.4444），并求出达到。 二代码1.function Sn=ComSimpson(a,b,n)%复合辛普森求积%f表示被积函数，本题中由f.m调用%a,b分别表示积分上下限%n表示区间分割次数%sn表示该方法计算所返回的值 h=(b-a)/n;for k=0:n x(k+1)=a+k*h; x_k(k+1)=x(k+1)+1/2*h; if (x(k+1)=0)|(x_k(k+1)=0) x(k+1)=10(-10); x_k(k+1)=10(-10)；%误差精度 endendS1=h/6*(f(x(1)+f(x(n+1);%复合辛普森求积公式第一部分for i=2:n F_1(i)=h/3*f(x(i);endfor j=1:n F_2(j)=2*h/3*f(x_k(j);endS2=sum(F_1)+sum(F_2);%复合辛普森求积公式第二部分Sn=S1+S2;%算得值 f的表达式由f.m（见附）文件调用。 附：function y=f(x) y=sqrt(x).*log(x); 3. 数值结果: 在命令窗口输入如下指令时，出现如下结果 Sn=ComSimpson(0,1,300) Sn =-0.4438 Sn=ComSimpson(0,1,700) Sn =-0.4442四计算结果的分析: 当步长取值很小的时候，误差较大。当步长取得越细，即区间分割的越小时，实验计算值的精度越高，即越趋近精确值。五. 计算中出现的问题，解决方法及体会： 本实验过程中，了解了复合求积公式的几个构成，以及在matlab中的实现，深化了对于该问题的理解。一的操作中，本来个个很简单的问题，哪怕就是一个小小符号的不注意，也不会得到结果。实验的严谨性，细腻性有待进一步提升。对于实验，网上的程序有不少，关键是要看懂，弄懂，在实际的
人人文库> 全部分类> 行业资料 > 信息产业

温馨提示

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

最新文档

评论

0/150

提交评论