数学软件与实验(第四次作业).docVIP

数学软件与实验 第四次作业 姓名：江海强 班级：070921 学号 时间：2011-5-31 内容 一. 实验目的： 本实验的主要目的是学会运用MATLAB对积分计算和对求解积分的左点法、右点法、中点法、梯形法、辛普森公式法、蒙特卡罗方法的掌握以及熟悉对积分图形打印等操作。通过本次实验，还让我们学会一些基本的编程，了解积分以及概率论等数学知识，学会随机模拟，并且学会统计分析。 二. 实验过程： 第一题： 定积分计算的基本公式是牛顿－莱布尼兹公式。但当被积函数的原函数不知道时，这时就需要利用近似计算。特别是在许多实际应用中，被积函数甚至没有解析表达式，而是一条实验记录曲线，或一组离散的采样值，此时只能用近似方法计算定积分。对于本题，主要利用符号积分法、左点法、右点法、中点法、梯形法和辛普森公式和蒙特卡罗方法从看，通过不断增加随即点数，结果越与真实值相符syms u v x y　　z=int(int(exp(-(u^2+v^2)),u,0,y),v,0,x) A=rand(1,100);satisfaction_book=A，people_nember=geornd(0.1,[1 365]) 得出一年365天每天借书的人数。 　　在做下面的步骤之前先假设抽样是放回抽样，对此先介绍几何分布函数，其公式为，其中期望和方差分别为：，，则，即每天平均有10个人借书。 然后从第一天开始，有people_nember(1)个人借阅书本，根据书的热门程度来进行选择借阅的书本，即热门度越高，被借阅的机率就越大。由于大家每次借阅一本，看完后即放在书架的最左边，且每天大概只有10人借阅书本，故可假设一个人借阅完，另一个人才开始借阅。则若第一个人借阅书本a，放回后a在第一位，若第二个人借阅书本b，放回后b在第一位，365天之后，得出书本热门度的排放顺序。但是由于是随机模拟，故整个实验需要多次执行，可在程序前面家上for循环，本附录中的程序已经加上for循环语句了，并且实验的次数为100次，再求其均值并以图形方式显示出统计结果。 其输出图形如下： 由图像可以大概的知道，热门度越大，排列越在前面，即越热门的书被借阅的机会越大，这是一个很显然的结果！ 以下给出了流程图： 三. 实验结果： 　　首先先介绍误差函数erf。其公式如下： 　　它具有如下的三条性质： 　　（1）； （2）； （3）对于均值为，方差为的高斯分布，其取值落在内的概率为． 所以，与输出s = (pi^(1/2)*erf(1))/2相符合。 符号积分法、左点法、右点法、中点法、梯形法和辛普森公式和蒙特卡罗方法 算出的值 相对误差 符号 积分法 0.74682413281242702539946743613185 左点法 0.747140131778599 右点法 0.746508011219770 中点法 0.746824163469049 梯形法 0.746824071499185 辛普森 公式法 0.746824132854452 蒙特卡罗方法 0.0024 　　 其中相对误差的公式如下： 例如得到左点法的相对误差为： 先算出，其中。再对进行50等分划分，得出矩阵X,Y,Z，最后作图surfc(X,Y,Z)，得出的图形如下所示： 四. 实验总结： 对于运用多种方法求解积分得出的结果及相对误差，从中可以得到：不同的方法有不同的计算精度，且中点法、梯形法和辛普森公式法一般比左点法和右点法的计算精度高，其中辛普森公式法的计算精度最高。而对于蒙特卡罗方法来说，通过不断增加随即点数，结果越与真实值相符 clear;clc format long %-----------（1）符号积分法---------- syms x s=int(exp(-x^2),0,1) disp(精确值:); s_symbol=vpa(s) %------------（2）左点法------------ h=1/1000; x=h.*[0:1000]; y=exp(-x.^2); disp(左点法估计值:); s_left=h*sum(y)-h*y(1001) %输出左点法求得的结果 disp(左点法估计值与精确值的相对误差:); d_left=abs((s_symbol-s_left)/s_symbol) %输出其相对误差 %------------（2）右点法------------ disp(右点法估计值:); s_right=h*sum(y)-h*y(1)%输出右点法求得的结果 disp(右点法估计值与精确值的相对误差:); d_right=abs((s_symbol-s_right)/s_