24小时热门版块排行榜    

查看: 1361  |  回复: 14
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

liao6417

新虫 (初入文坛)

[交流] 【求助】求一段monte carlo 程序 已有1人参与

求一段monte carlo 程序,不管是实现什么的,我主要是想看看这个运算法则的样子,以及写法

[ Last edited by wuli8 on 2009-11-13 at 21:40 ]
回复此楼

» 猜你喜欢

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

zhiping780

木虫 (著名写手)

天鹰帮帮主


小木虫(金币+0.5):给个红包,谢谢回帖交流
哥们   你搞monte carlo咋学的?我想搞    却一直不懂这东西杂用
目云天外,志在九霄中。。。。。
5楼2009-11-22 19:25:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 15 个回答

pfzhangphy

金虫 (小有名气)


lei0736(金币+1):谢谢 2010-03-15 11:07
建议找本MC的书看吧!看程序不如看书好!
2楼2009-11-16 14:28:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bay__gulf

金虫 (著名写手)

刘苏州

★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
zeoliters(金币+4,VIP+0):多谢回帖交流! 11-17 15:04
MC方法很好,只是形式复杂的力场不好加进去

当初学习MC时候写一个例子
用于逃离曲线上 local minimum, 找到的global minimum
实验函数y = cos (x) + cos (3x), 图像见图top
最低点-2.0, 收敛目标-1.99999
代码
CODE:
#include
#include
#include
int main()
{
  float a, d, t;
  float x0, x, y0, y;
  scanf("%f %f %f", &a, &d, &t);
  srand(0);
  x0=a;
  do {
    y0= cos(x0)+cos(3*x0);
    x= x0 + d*(rand()%10000/10000.0-0.5);
    y= cos(x) + cos(3*x);
    printf("%f\t%f\n", x, y);
    if (exp((y0-y)/t) > rand()%10000/10000.0) x0= x;
  } while (y > -1.99999);

  return 0;
}

需要手工输入三个参数,
a: x初值, d:尝试移动的最大范围, t: 相当于kT
图middle 和buttom 分别在(a,d,t)=(0.2, 0.1, 0.8) 时候, x 和y 随运行步数的演化
可见在有限的时间内, 通过简单的算法可以逃离曲线上 local minimum, 找到的global minimum
适合分子复杂, 不能一一列举的各个构象时候的构象分析
3楼2009-11-16 14:58:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wuchenwf

荣誉版主 (职业作家)

引用回帖:
Originally posted by bay__gulf at 2009-11-16 14:58:
MC方法很好,只是形式复杂的力场不好加进去

当初学习MC时候写一个例子
用于逃离曲线上 local minimum, 找到的global minimum
实验函数y = cos (x) + cos (3x), 图像见图top
最低点-2.0, 收敛目标-1.99999
代 ...

好资源!
4楼2009-11-16 21:52:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 南昌大学材料专硕311分求调剂 +6 77chaselx 2026-03-20 6/300 2026-03-21 07:24 by JourneyLucky
[考研] 265求调剂 +3 Jack?k?y 2026-03-17 3/150 2026-03-21 03:17 by JourneyLucky
[考研] 301求调剂 +10 yy要上岸呀 2026-03-17 10/500 2026-03-21 03:14 by JourneyLucky
[考研] 二本跨考郑大材料306英一数二 +3 z1z2z3879 2026-03-17 3/150 2026-03-21 02:29 by JourneyLucky
[考研] 307求调剂 +10 冷笙123 2026-03-17 10/500 2026-03-21 01:54 by JourneyLucky
[考研] 317求调剂 +8 申子申申 2026-03-19 13/650 2026-03-21 00:09 by 刘国森
[考研] 274求调剂 +10 S.H1 2026-03-18 10/500 2026-03-20 23:51 by JourneyLucky
[考研] 321求调剂 +9 何润采123 2026-03-18 11/550 2026-03-20 23:19 by JourneyLucky
[考研] 085600材料与化工 +8 安全上岸! 2026-03-16 8/400 2026-03-20 22:13 by luoyongfeng
[考研] 材料与化工 322求调剂 +4 然11 2026-03-19 4/200 2026-03-20 22:12 by luoyongfeng
[考研] 材料与化工专硕调剂 +7 heming3743 2026-03-16 7/350 2026-03-20 19:31 by zhukairuo
[考研] 0856调剂,是学校就去 +8 sllhht 2026-03-19 9/450 2026-03-20 14:25 by 无懈可击111
[考研] 320求调剂0856 +3 不想起名字112 2026-03-19 3/150 2026-03-19 22:53 by 学员8dgXkO
[考研] 生物学调剂招人!!! +3 山海天岚 2026-03-17 4/200 2026-03-19 21:34 by 怎么释怀
[考研] 0703化学调剂 +4 18889395102 2026-03-18 4/200 2026-03-19 16:13 by 30660438
[考研] 化学工程321分求调剂 +15 大米饭! 2026-03-15 18/900 2026-03-18 14:52 by haxia
[考研] 290求调剂 +3 p asserby. 2026-03-15 4/200 2026-03-17 16:35 by wangkm
[考研] 0854控制工程 359求调剂 可跨专业 +3 626776879 2026-03-14 9/450 2026-03-16 17:42 by 626776879
[考研] 0703化学调剂 290分有科研经历,论文在投 +7 腻腻gk 2026-03-14 7/350 2026-03-16 10:12 by houyaoxu
[考研] 070305求调剂 +3 mlpqaz03 2026-03-14 4/200 2026-03-15 11:04 by peike
信息提示
请填处理意见