[Alg]利用蒙特卡洛思想计算圆锥体积

祥知道

于 2017-12-03 15:36:20 发布

阅读量2.7k

点赞数

分类专栏：算法

本文链接：https://blog.csdn.net/humanking7/article/details/78701954

版权

算法专栏收录该内容

4 篇文章

订阅专栏

本文介绍了一种使用蒙特卡洛方法来估算圆锥体积的Matlab实现方式。通过大量的随机抽样点判断其是否位于给定圆锥内，以此计算圆锥体积并对比数学公式得出的精确值。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

原创文章，欢迎转载。转载请注明：转载自祥的博客

原文链接：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)%蒙特卡洛思想算的体积与公式计算的体积的误差

donate