数值积分与微分学习笔记.docVIP

数值积分与微分 实验目的：1）用matlab软件掌握梯形公式、辛普森公式和蒙特卡罗方法计算数值积分； 2）通过实例学习用数值积分和数值微分解决实际问题。 实验内容： 第一题：用梯形、辛普森和蒙特卡罗方法计算积分。改变步长（对梯形），改变精度要求（对辛普森），改变随机点数目（对蒙特卡罗），进行比较、分析。 y=e ,-2x2 解：用三种方法计算积分的源程序如下： h=4/50; x=-2:h:2; y=fun(x); z1=trapz(y)*h z2=quad(fun,-2,2) plot(x,y,g); n=1000; x1=rand(1,n); y1=fun(x1.*2); z3=sum(y1)*4/n 对梯形公式取h=4/50,4/100,4/10000;对辛普森分别取精度为10，,10；对蒙特卡罗方法分别取n=1000,10000,100000.得到的结果如下： 梯形公式 辛普森公式 蒙特卡罗方法 0.95438456767789 0.95449943824154 0.93999211059586 0.95447094168964 0.95449973610735 0.95794338594866 0.95449973322412 0.95449973610373 0.95427317381756 从得到的结果可以看到对梯