24小时热门版块排行榜

返回列表

当前主题已经存档。

【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 anxt2006 的 5 个金币

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

anxt2006

木虫 (正式写手)

应助: 0 (幼儿园)
金币: 6086.7
散金: 38
红花: 1
帖子: 377
在线: 208.2小时
虫号: 578555
注册: 2008-07-04
专业: 凝聚态物性 II ：电子结构

[交流] [求助]fortran里如何产生从0到1的加权的随机数？

fortran里如何产生从0到1的加权的随机数？我想用蒙特卡罗积分，但是我的被积函数有尖峰，所以想找一个加权的随机数产生程序，在尖峰的地方多取一些点，在其他地方少取一些点。请高手指教。

[ Last edited by kuhailangyu on 2008-12-9 at 10:32 ]

回复此楼

» 猜你喜欢

1楼 2008-10-13 17:05:56

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

maomao1210

金虫 (正式写手)

程序强帖: 5
应助: 2 (幼儿园)
金币: 1431.3
散金: 242
红花: 16
沙发: 1
帖子: 991
在线: 441.5小时
虫号: 253215
注册: 2006-05-20
性别: MM
专业: 考古理论

SUBROUTINE(JIE,X)
IMPLICIT REAL*8 (A-H,O-Z)
   INTEGER*4 I(2)
C          READ(*,*) JIE
         IF (JIE.EQ.0) GO TO 1
   I(1)=0
   I(2)=JIE
  1    IX=1490*I(2)+3805*I(1)
   IE=3805*I(2)
   I(2)=IE-IE/32768*32768
   IX=IX-IX/32768*32768+IE/32768
   I(1)=IX-IX/32768*32768
   X=(I(1)*32768+I(2))*9.3132258D-10
   END
RETURN
END
这样就行了，你调用该程序的时候，给定一个种子JIE，然后这个程序就会给你一个数值X。

[ Last edited by maomao1210 on 2008-10-14 at 22:15 ]