实验五_蒙特卡罗方法.ppt

1、，6 :蒙特卡洛方法实验，二维填充图绘制方法区域，体积计算问题冰激凌圆锥的体积计算事故问题和练习问题，fill(x，y，c)用c指定的颜色填充多边形。其中多边形的顶点由x，y确定。c指定的颜色不允许省略。由x，y确定的点必须形成闭合多边形。二维多边形填充fill()是格式，1.2维填充绘制方法，clfN=10dt=2 * pi/n:T=T=0: dt :2 * pi；T=t，t(1)；x=sin(t)；y=cos(t)；菲尔(x，y，c)；Axis off，rand生成介于0和1之间的均匀随机数。rand(m，n)生成Mn 0到1之间的均匀随机数，生成X=rand(10000，1)。Hist(

2、X)，(1)蒙特卡洛方法使用随机实验计算近似值，10，000个随机数在每个小区之间均匀分布，随机变量X下降到小区的概率仅与小区的长度相关，与小区的位置无关。2.面积，体积计算问题，示例1。计算以下两条曲线围成的区域的面积：End X=2*rand(N，1)-1；Y=2*rand(N，1)；II=find(Y=ABS(X)；M=长度(ii)；s=4 * m/N:X1=0:0.01:1；Y1=x1x2=1:-0.01:-1；y2=1 sqrt(1-x 2.2)；x3=-1:0.01:0；y3=-x3；Fill (x1、x2、x3、y1、y2、y3、c)，s=2.5460，范例2。计算半径为1的两个

3、直圆柱面的边界体积。x2 y2=1，x2 z2=；End P=rand(N，3)；X=P(:1)；Y=P(:2)；Z=P(:3)；Ii=find (x.2 y.2=1 meshz (x，y，ZZ)color map(0 0 1)axis off view(119，34)，v=；f是由积分函数名称组成的字符串。相关命令：dblquad()积分计算，示例3。积分计算，fun=inline(sqrt(x)-x . 2)；S=四元(fun，0.01，1)t=0:0.01:1；y=fun(t)；Fill (0，t，0，y，c) syms u，s0=rand (sqrt (u)-U2，0，1)，s=0.3

4、327 s0=1因此极坐标变换处理重点的方法，4 .试题和练习题，3。圆x2 (y 2)2=4内的小圆x2 (y 1)2=1后的图形，h=2 * pi/100；t=-pi/23360h 33603 * pi/2；x=cos(t)；y=1 sin(t)；x=2 * cos(t)；y=2 2 * sin(t)；N=长度(x)；x=x(n :-1:1)；y=y(n :-1:1)；Fill(x，X，Y，Y，c) axis off，大圆，小圆，4。蒙特卡洛方法在计算体积和计算面积方面有很多困难吗？5.介绍了用蒙特卡洛方法计算冰激凌圆锥体积和计算莫合盖问题的异同。6.在积分学中，叙述了两个茄子著名的问题：

5、莫合开问题和薇薇安问题。h=2 * pi/100；T=0:h:piR=0:0.0533601x=r * cos(t)；y=r * sin(t)；ZZ=sqrt(1-x . 2)；Meshz (x，y，ZZ)color map(001)axis off view(119，34)，7在区域1中投入2000个均匀随机点后，随机点落入抛物线包围区域的概率是请求面积与正方形面积的比率。functions=area1 (n) if nargin=0，n=2000End X=rand(N，1)；Y=rand(N，1)；II=find(Y=x . 2)；M=长度(ii)；s=m/N:X1=0:0.01:1；x

6、2=1:-0.01:0；y1=sqrt(x1)；Y2=x2.2使用Fill(x1，x2，y1，y2，c)、S=0.3333，(补例2) :随机投资点方法计算圆周率。座标原点为中心，第一个象限内的单位圆部分恰好为0，10，正方形Y=rand(N，1)；II=find(x . 2y . 2=1)；M=长度(ii)；P=4*m/N : Plot (x (ii)，y (ii)，b)，hold on theta=0:01: pi/2；x=cos(theta)；y=sin(theta)；Plot (x，y，r)，3.1560，3.1520，3.1400，积分计算，fun=inline(sqrt(1-x . 2)；S=