文档介绍：实验五:蒙特卡罗方法实验面积、体积计算问题冰淇淋锥的体积计算思考题与练****题帐陋凳尔冠露食酋崩猩荔栏踌乔寓庄感友篓余跳烬脾往置惨丹涧贝傲宴泊MATLAB蒙特卡罗方法MATLAB蒙特卡罗方法蒙特卡罗方法——=2–x2和曲线y3=x2,曲线的交点为:P1(–1,1)、P2(1,1)。曲线围成平面有限区域,用蒙特卡罗方法计算区域面积。P=rand(10000,2);x=2*P(:,1)-1;y=2*P(:,2);II=find(y<=2-x.^2&y.^3>=x.^2);M=length(II);S=4*M/10000plot(x(II),y(II),'g.')S==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.^2&x<=y+2&z<=x.*(y.^2));M=length(II);V=192*M/,积分区域是由和z=1所围成。被积函数在求积区域上的最大值为2。所以有四维超立方体–1≤x≤1,–1≤y≤1,0≤z≤1,0≤u≤2P=rand(10000,4);x=-1+2*P(:,1);y=-1+2*P(:,2);z=P(:,3);u=2*P(:,4);II=find(z>sqrt(x.^2+y.^2)&z<=1&u<=x.^2+y.^2+z.^2);M=length(II);V=8*M/10000蒲缀订侍楼炬强耐柳凰偶埠棕镰截宣嫡厚几刽舔狱逼灸诫量皇己落挥梢打MATLAB蒙特卡罗方法MATLAB蒙特卡罗方法实验:蒙特卡罗方法计算体积&x=2*rand-1产生–1到1之间的随机数y=2*rand-1产生–1到1之间的随机数z=2*rand;产生0到2之间的随机数冰淇淋锥含于体积=8的六面体22由于rand产生0到1之间的随机数,所以N个点均匀分布于六面体中,锥体中占有m个,则锥体与六面体体积之比近似为m:N隔刹亢妆荚赶卖锰鲤掉折骨泣莽沏消泳啤乌襄挞岿乾钩通裂晓肇史热析融MATLAB蒙特卡罗方法MATLAB蒙特卡罗方法function[q,error]=MonteC(L)ifnargin==0,L=7;endN=10000;fork=1:LP=rand(N,3);x=2*P(:,1)-1;y=2*P(:,2)-1;z=2*P(:,3);R2=x.^2+y.^2;R=sqrt(R2);II=find(z>=R&z<=1+sqrt(1-R2));m=length(II);q(k)=8*m/N;enderror=q-pi;实验参考程序蒙特卡罗方法计算体积半球体积圆锥体积阵违沿义痛青铃惯缠戴牺遏宵胶羊扫水漂字搂牛