MATLAB数值积分求值实验报告_第1页
MATLAB数值积分求值实验报告_第2页
MATLAB数值积分求值实验报告_第3页
全文预览已结束

下载本文档

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

文档简介

学生身份证班级统计1001名字教师易坤南实验主题用各种方法计算数值积分评估点1.设计目的(实习):1.理解MATLAB在实际问题中的应用2.通过实践加深对该语言中M文件的理解3.熟悉简单的程序结构，如循环结构(for循环，while循环)选择结构(if-else-if)、分支语句(switch-case-other)。2.实验内容：(1)。用左右矩形法、梯形法和辛普森公式计算了0，1上Y=x 2的定积分。(2)用蒙特卡罗随机投影法计算0，1上Y=1/(12)的定积分，得到pi的近似值；(3)用蒙特卡罗平均估计方法计算0，1上Y=x 2的定积分。3.详细设计：一、左右矩形法和梯形法：h=1/200；x=0:h:1。y=x.2；Z1=总和(y(1:200)*h%左矩形法Z2=总和(y(2:201)*h%右矩形法z=cum sum(y)；z11=z(200)* h；%等于z1z12=(z(201)-z(1)* h；%等于z2z3=梯形法，相当于z3=梯形法(y)*h2.重复辛普森公式方法：y=inline(x.2)；Z1=四(y，0，1，100) %simpleson公式Z2=四次方(y，0，1，100)%复数simpleson公式Z3=quad8 (y，0，1，100，轨迹(10)% simpleton 8阶公式法3.蒙特卡洛随机射击法；n=100000。k=0；对于i=1:nx=兰特；%为(0，1)间隔生成随机数y=兰特；如果y=1/(1 x2)；%至Y=1/(1 x 2)面积投掷点k=k1；目标目标z=k/nPi=4*k/n%由积分pi/4=k/n导出，前者是概率，后者是频率4.蒙特卡罗平均估计方法；n=10000。x=兰特(1，n)；y=x.2；z=总和(y)/n%均值估计方法4.实验结果：1.左右矩形法和梯形法：z1=0.3308z2=0.3358z3=0.33332.重复辛普森公式方法：z1=0.3333z2=0.3333警告： QUAD8已过时。我们用QUADL代替。在第八区35号在简单的418 0.00000000000 5.00000000 e-001 0.33333333z3=0.33333.蒙特卡洛随机射击法；z=0.7854pi=3.14164.蒙特卡罗平均估计方法；z=0.33655.实验总结有些规则与C语言的有些不同，但是结构形式和大多数语句是相同的，所以对于C语言来说也是可以接受的。通过各种求定积分的方法，我已经掌握了积分的计算，这对以后的学习很有帮助。然而，对于辛普森的复杂解决方案仍有一些疑问。这样的函数能直接调用内置函数吗？我觉得我已经掌握了一些基本的语言规则，这对我将来的学习非常有帮助，我感到非常
人人文库> 全部分类> 行业资料 > 管理策划

温馨提示

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

最新文档

评论

0/150

提交评论