基于matlab的蒙特卡罗积分的实现.doc

概率论与数理统计论文基于matlab的蒙特卡罗定积分的实现 班级： 姓名： 基于matlab的蒙特卡罗定积分的实现摘 要：在对蒙特卡罗方法概念学习的基础上，利用计算机产生了随机数序列在此基础上，研究了蒙特卡罗积分，并应用matlab加以实现关键词：蒙特卡罗；MATLAB；定积分0 引言随着电子计算的出现和发展，近年来用概率模型来作近似计算的方法得到了很大的发展，即蒙特卡罗(MonteGarlo)方法它是一种采用统计抽样理论近似的求解数学与物理问题的方法，它既可以用来研究概率问题，也可以用来解决非概率问题蒙特卡罗法已被广泛地运用到各个领域中，如高维数学问题求解、医学技术中的诊断识别、大型系统的可靠性分析等一般蒙特卡罗方法在数学中最常见的应用就是蒙特卡罗积分对于那些由于计算过于复杂而难以得到解析解或者根本没有解析解的问题，蒙特卡罗方法能够提供一种有效而可行的解决方法1 蒙特卡罗方法的基本原理用蒙特卡罗方法解决数学分析问题时，基本思想是：首先建立与描述该问题有相似性的概率模型，利用这些相似把使该模型的一些数字特征(如概率或均值)与数学分析的解答(如积分值或微分方程的解)联系起来；然后对该模型进行随机模拟或统计抽样；最后利用所得的结果求出这些特征的统计估计值作为原来问题的近似解蒙特卡罗方法计算的结果收敛的理论依据来自于大数定律，且结果渐进地服从正态分布的理论依据是中心极限定理本文将对蒙特卡罗方法计算定积分，并用matlab加以实现2 蒙特卡罗方法求定积分的一个实例21 采用均匀随机数的蒙特卡罗方法计算和 matlab 实现 为了使用蒙特卡罗方法求该积分，首先选定一个矩形区域 D =0，1；O，1然后抽取n个随即点的坐标(i=1，2，n)，他们服从区域D上的二维均匀分布计算落在区域D内的随机点数设为m，即满足0则由强大数定律，要求的定积分近似的等于两个区域随机落点的比值 与矩形面积的乘积 matlab程序及其运行结果：22 matlab生成随机掷点效果图x1=rand(1,2000);y1=rand(1,length(x1);x=0:0.01:1;y2=0;y3=x.3;subplot(2,1,1),plot(x,y2,k,x,y3,k,x1,y1,.)xlabel(x轴),ylabel(y轴)m=1;for j=10:10:30000;a(m)=j;x0=rand(1,j);y0=rand(1,length(x0);s=0;n=length(x0);for i=1:n y4(i)=0;y5(i)=x0(i).3;if y0(i)y4(i) s=s+1;endendp(m)=s/n;m=m+1;endt=mean(p);txx=0;length(x0);yy=1/4;subplot(2,1,2),plot(a,p,b,xx,yy,b)t 0.249923 matlab均匀随机数法计算定积分由以上运行结果可知，均匀随机数法求的定积分 0.2499由定积分的几何意义可知，它是三次曲线与直线y = 0所围的平面图形的面积，由微积分基本公式可知，模拟值与精确值的误差率为0.43 小结Monte Carlo方法的实质是通过大量随机试验，利用概率论解决问题的一种数值方法，基本思想是基于概率和体积间的相似性，因此程序十分简单；蒙特卡罗方法在解决确定性问题时，具有较高的精确度；误差主要取决于样本的容量而与维数无关，随着抽取点增加，近似面积也将逼近真实面积；对问题求解的过程取决于所建立的概率模型参 考 文 献：1 概率论与数理统