实验五 MATLAB蒙特卡罗方法.pptVIP

*/10 ? ? ? 实验五:蒙特卡罗方法实验 面积、体积计算问题 冰淇淋锥的体积计算 思考题与练习题 蒙特卡罗方法——随机投点试验求近似解 引例. 给定曲线y =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.^2y.^3=x.^2); M=length(II); S=4*M/10000 plot(x(II),y(II),g.) S = 2.1136 例5.14 计算 其中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 例5.15 用蒙特卡罗方法计算 其中,积分区域是