实验五蒙特卡罗方法2概述.ppt

*/10 ? ? ? 实验五:蒙特卡罗方法实验 面积、体积计算问题 维维安尼的体积计算 思考题与练习题 例5.4 计算 其中D为y= x – 2与y2 = x 所围 D的边界曲线交点为：(–1，1)，(4，2)，被积函数在求积区域内的最大值为16。积分值是三维体积，该三维图形位于立方体区域 0≤ x ≤4，–1≤ y ≤2，0 ≤ z ≤16 内，立方体区域的体积为192。 data=rand(10000,3); x=4*data(:,1); y=-1+3*data(:,2); z=16*data(:,3); II=find(x=y.^2x=y+2z=x.*(y.^2)); M=length(II); V=192*M/10000 例7.求 | ln x | + | ln y | = 1所围平面图形的面积。 分析：由曲线方程，| ln x | ≤1，| ln y | ≤1， 1/e ≤ x ≤ e，1/e ≤ y ≤ e 正方形域：D ={(x，y) | 0.3≤x≤2.8，0.3≤y ≤2.8} S=(2.8-0.3)^2;Sk=[]; for k=1:5 data=0.3+2.5*rand(10000,2); x=data(:,1); y=data(:,2); II=find(abs(log(x))+abs(log(y))=1); M=length(II); s=S*M/N;Sk=[Sk,s]; end Sk 二维多边形填充图 fill() 使用格式 用c所指定的颜色对多边形填充.其中,多边形的顶点由x,y确定. c 指定颜色不允许省略;x,y确定的点必须形成封闭的多边形. fill(x,y,c) 填充图绘制方法: x1=0:.01:1;y1=sqrt(x1); x2=1:-.01:0;y2=x2.^2; fill([x1,x2],[y1,y2],r) y1=-1:.1:2;y2=2:-.1:-1; x11=y1.*y1;x22=y2+2; fill([x11,x22],[y1,y2],r) 例8 y =x2 , x = y 2 所围区域 例9 y2 = x与 y= x – 2 所围区域 例 10. | ln x | + | ln y | = 1所围区域 令 ? x1=linspace(-1,0,10);y1=x1+1; x2=linspace(0,1,10);y2=1-x2; x3=x2(10:-1:1);y3=-y2(10:-1:1); x4=x1(10:-1:1);y4=-y1(10:-1:1); xt=[x1,x2,x3,x4]; yt=[y1,y2,y3,y4]; X=exp(xt);Y=exp(yt); figure(3),fill(X,Y,r) 实验:维维安尼体体积计算 实验:蒙特卡罗方法计算体积 x=R*rand产生0到R之间的随机数 y=R*rand-R/2产生– R/2到R/2之间的随机数 z=-R+2*R*rand;产生0到R之间的随机数 维维安尼含于体积 = 2R3 的六面体 由于rand 产生0 到1之间的随机数,所以 N个点均匀分布于六面体中,维维安尼占有m个,则维维安尼与六面体体积之比近似为 m : N 实验任务一：记录第L次实验的实验数据及误差R=2,N=10000 实验任务二：计算L次实验数据均值及均值误差( mean 计算平均值 ) 序号 1 2 3 4 5 6 7 数据 误差 L 8 16 32 64 128 256 均值 误差 [a,b,c]=cylinder(1,100); c(1,:)=-2; c(2,:)=2; figure(2) mesh(a+1,b,c) hold on [u,v]=meshgrid(0:0.2:2*pi); x=2*cos(u).*cos(v); y=2*cos(u).*sin(v); z=2*sin(u); mesh(x,y,z) 维维安尼图形绘制程序 4. 说明L次蒙特卡罗实验平均值计算维维安尼的体积误差与实验次数之间关系 1.说明蒙特卡罗方法计算维维安尼的体积误差与哪些因素有关。 2.概率论中的贝努里大数定律在本次实验中体现如何 3.叙述概率论中著名独立同分布大数定律，并以这一大数定律解释实验的数据变换规律。 思考题与练习题 5.下面程序绘出二维图形填充图(右图)。分析每条语句功能,给程序中语句写注记 x1=-1:0.1:1; y1=x1.^2.