作者ttvic ( 短 维)
看板Fortran
标题[问题] 圆锥体积 蒙地卡罗法
时间Mon Dec 23 23:32:10 2013
最近碰到需要使用蒙地卡罗法找出圆锥的体积
示意图:
https://www.dropbox.com/s/gemf3q3yqrmg1lt/cone.jpg
想法: 利用圆柱来抓出1/3的圆锥，但投出机率是错误的...
以下是我的程式
nh=0
do i=1, n
call random_number (rand)
r=rand
call random_number (rand)
a=2*3.1415926*rand !角度一 0~2*pi
call random_number (rand)
b=2*3.1415926*rand !角度二 0~2*pi
call random_number (rand)
z=rand !高度 0~1
x=r*cos(a)
y=r*sin(b)
if (z>=f1(x,y)) then
nh=nh+1
endif
enddo
pr=(real(nh)/real(n)) !←
主要是这里的机率一直抓不到
area1=(real(nh)/real(n))*3.1415926*1**2
write(*,*) 'con area',area1
!external function
real function f1(x,y)
f1=sqrt(x**2+y**2)
end function
不晓得边界的界定是错在哪边? 恳请大大帮忙~谢谢!
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 140.112.249.76
1F：推 iinon:可以有analytical 解为什麽不用? 12/26 15:05
2F：推 KAOKAOKAO:蒙地卡罗的练习啊 用解析解还用练习吗 01/18 21:26
3F：推 bonbbon:x=-1~1,y=-1~1,这个程式算出来的是角锥! 03/19 02:31
4F：→ bonbbon:y=r*sin(a) 才会是圆锥 03/19 02:32