蒙特卡罗的方法计算三重积分.docVIP

蒙特卡罗的方法计算三重积分 　　【摘要】应用蒙特卡罗方法的随机投点法原理求解三重积分时，传统的方法是构造四维长方体，然后在该长方体内投点求解，本文创新提出根据被积函数的要求构造出包含积分区域的四维圆柱体并计算。数值结果表明：四维圆柱区域更贴近于原积分区域，此方法计算效果更好，在三重积分的计算上是个行之有效的方法。 　　【关键词】蒙特卡罗法；三重积分；随机投点法；四维圆柱体方法 　　Calculation of Triple Integral with Monte-Carlo Method 　　WANG Hong-taoLI Man-zhiSHEN You-jian 　　(School of Mathematics and Statistics,Hainan Normal University, Haikou Hainan 571158, China) 　　【Abstract】This paper briefly introduces the application of Monte Carlo(MC) method on computing triple integral.it gives one new method by structuring 4-dimension cylinder containing the integral domain. Results of the numerical calculations show that the method is faster and effective than the classic one. 　　【Key words】Monte Carlo method;Triple integral;Randomly cast-point method;4-dimention cylinder method 　　1蒙特卡罗法的基本原理 　　蒙特卡罗方法以随机模拟和统计试验为手段，从随机变量的概率分布中，通过选择随机数的方法产生一种符合该随机变量概率分布特性的随机数值序列，作为输入变量序列进行特定的模拟试验、求解的方法[1]。 在应用蒙特卡罗方法时，必须要产生非均匀的随机数，而产生特定的、非均匀的随机数序列，可行的方法是先产生一种均匀分布的随机数序列，然后设法转换成我们需要的随机数序列并以此作为数字模拟试验的输入变量序列进行模拟求解[2-6]。 　　根据三重积分的数学定义，三重积分相当于密度函数为f(x,y,z)的空间物体Ω的质量，可将质量假设为四维空间区域的体积。因为随机投点法需要构造一个能包含积分区域的四维空间区域，我们需要根据该空间区域的体积来计算，因此该方法又叫取体积法。 　　2蒙特卡罗法计算三重积分的原理及数值算法 　　2.1原理 　　定理[2,6-13]设f(x,y,z)为区域Ω上的有界函数，对于三重积分I=■f(x,y,z)dv,其中Ω∶a1≤x≤b1,a2(x)≤y≤b2(x),a3(x,y)≤z≤b3(x,y)，满足 　　1）a2(x),b2(x)在区间[a1,b1]上连续，a3(x,y),b3(x,y)在区域 　　D∶a1≤x≤b1,a2(x)≤y≤b2(x)上连续，且满足不等式 　　a2(x)≥c,b2(x)≤d,a3(x,y)≥e,b3(x,y)≤h 　　2）设M≥maxf(x,y,z)； 　　3）设四维超长方体?追由a1≤x≤b1,c≤y≤d,p≤z≤q,0≤g≤M所围成，体积为V?追； 　　4）区域?祝由a1≤x≤b1,a2(x)≤y≤b2(x),a3(x,y)≤z≤b3(x,y),0≤g≤f(x,y,z)所围； 　　5）在四维超长方体?追内产生N个均匀随机点(xi,yi,zi,gi),i=1,2,…,N，设(xi,yi,zi,gi),i=1,2,…,n为落在?祝中的n个随机数， 　　则当N充分大时，有 　　■ 　　2.2数值算法 　　１）赋初值：令落在区域?祝中的点的个数n=0;规定投点试验的总次数N; 　　2）产生四组相互独立的均匀随机数列ξ,η,ω,ψ~U(0,1)由x=a1+(b1-a1)ξ,y=c+(d-c)η,z=p+(q-p)ω,g=Mψ求xi,yi,zi,gi(i=1,2,…,N)； 　　3）选出满足a2(xi)≤yi≤b2(xi)的yi(i=1,2,…,N)，设它们是yki(i=1,2,…,m)； 　　4）选出满足a3(xi,yi)≤zki≤b3(xi,yi)且与yki(i=1,2,…,m)相对应的zki(i=1,2,…,m)，设它们是zki(i=1,2,…,l)； 　　5）选出满足gki≤f(xi,yi,zi)且与zki(i=1,2,…,l)相对应的gki(i=1,2,…