实验二-蒙特卡罗方法计算三维体积.docVIP

实验二-蒙特卡罗方法计算三维体积 实验二-蒙特卡罗方法计算三维体积 PAGE 实验二-蒙特卡罗方法计算三维体积 数学与应用数学系201 3 ～201 4 学年第 二 学期实验报告 班级： 学号： 姓名： 实验时间： 2014年 月 日 实验 项目 实验二 蒙特卡罗方法计算三维体积 所属 课程 数学实验 实 验 目 的 了解蒙特卡罗方法的原理，掌握随机数使用技术。 实 验 内 容 用两种不同的方法（蒙特卡罗法和重积分方法）计算下列三维体积和重积分，并比较计算分析结果。 1、画出由锥面 上方与球面 内部区域围成的图形（简称冰淇淋锥），并计算也该冰淇淋锥的体积 2、画出积分区域并计算，其中积分区域是由和所围成。 3、画出积分区域并计算，其中D为y= x – 2与y2 = x 所围，D的边界曲线交点为：(1，-1)，(4，2)。 实 验 结 果 ︵ 不 够 可 另 附 纸 ︶ (*1冰淇淋锥的体积*) t1=ParametricPlot3D[{r*Cos[t],r*Sin[t],r^2},{t,0,2Pi},{r,0,1},DisplayFunction?Identity]; t2=ParametricPlot3D[{Cos[u]*Sin[v],Sin[u]*Sin[v],1+Cos[v]},{u,0,2Pi},{v,0,Pi/2},DisplayFunction?Identity]; Show[t1,t2,DisplayFunction?$DisplayFunction] ?Graphics3D? (*1冰淇淋锥的体积*) n=1000; p=Complex{}; Do[m=0; Do[x=2*Random[Real,{0,1}]-1;y=2*Random[Real,{0,1}]-1;z=2*Random[Real,{0,1}]; R1=x^2+y^2; R2=Sqrt[R1]; If[z?R2 (z-1)^2?1-R1,m++],{k,1,n}]; AppendTo[p,N[8m/n]],{t,1,10}]; Print[p]; Sum[p[[t]],{t,1,10}]/10 {,,,,,,,,,} (*1冰淇淋锥的体积*) x=r*Sin[u]*Cos[v]; y=r*Sin[u]*Sin[v]; z=r*Cos[u]+1; s=Integrate[r^2*Sin[u],{v,0,2Pi},{u,0,Pi/4},{r,0,2Cos[u]}]; N[s] (*2体积*) s1=ParametricPlot3D[{r*Sin[u],r*Cos[u],r},{u,0,2Pi},{r,0,1},DisplayFunction?Identity]; t1=ParametricPlot3D[{r*Sin[u],r*Cos[u],1},{u,0,2Pi},{r,0,1},DisplayFunction?Identity]; Show[s1,t1,DisplayFunction?$DisplayFunction] ?Graphics3D? (*2体积*) n=1000; p=Complex{}; Do[m=0; Do[x=2*Random[Real,{0,1}]-1;y=2*Random[Real,{0,1}]-1;z=Random[Real,{0,1}]; u=2*Random[Real,{0,1}]; R1=x^2+y^2; R2=Sqrt[R1]; If[z?1 z?R2 u?R1+z^2,m++],{k,1,n}]; AppendTo[p,N[8m/n]],{t,1,10}]; Print[p]; Sum[p[[t]],{t,1,10}]/10 {,,,,,,,,,} (*2体积*) s=Integrate[x^2+y^2+z^2,{y,-1,1},{x,-Sqrt[1-y^2],Sqrt[1-y^2]},{z,Sqrt[x^2+y^2],1}]; N[s] (*3体积*) s=Plot[x-2,{x,1,4},DisplayFunction?Identity]; t=Plot[Sqrt[x],{x,1,4},DisplayFunction?Identity]; Show[s,t,DisplayFunction?$DisplayFunction] ?Graphics? n=10000;(*3体积*) p=Complex{}; Do[m=0; Do[x=4*Random[Real,{0,1}];y=3*Random[Real,{0,