原创文章,欢迎转载。转载请注明:转载自 祥的博客
原文链接:http://blog.csdn.net/humanking7/article/details/78701954
话不多说,直接上Matlab代码
%% 用蒙特卡洛思想算圆锥体积
clear
clc
%%
R = 0.5;
H = 0.1;
n = 10000000;%R和H越小,模拟次数越大,精确度越好
Pr = unifrnd(-R,R,2,n);
Ph = unifrnd(0,H,1,n);
cnt = 0;
for i = 1:n
r = sqrt(Pr(1,i)^2 + Pr(2,i)^2);
h = H-Ph(i);%坐标系不一样,也可以 h=Ph(i)
if (r/R)<(h/H)
cnt = cnt+1;%说明点在圆锥里面
end
end
V_f = 1/3*pi*R^2*H;%公式计算圆锥体积
V_rect = (2*R)^2*H;%公式计算长方体体积
scale = cnt/n;%圆锥体积占长方体的比例
V_m = scale*V_rect%利用蒙特卡洛思想计算圆锥体积
error_cmp = abs(V_m - V_f)%蒙特卡洛思想算的体积与公式计算的体积的误差