24小时热门版块排行榜    

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

liao6417

新虫 (初入文坛)

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

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

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

» 猜你喜欢

已阅   回复此楼   关注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的回帖
查看全部 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的回帖

zhiping780

木虫 (著名写手)

天鹰帮帮主


小木虫(金币+0.5):给个红包,谢谢回帖交流
哥们   你搞monte carlo咋学的?我想搞    却一直不懂这东西杂用
目云天外,志在九霄中。。。。。
5楼2009-11-22 19:25:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 广西大学材料导师推荐 +3 夏夏夏小正 2026-03-17 5/250 2026-03-21 22:20 by 金昊ML
[考研] 求调剂 +4 要好好无聊 2026-03-21 4/200 2026-03-21 18:57 by 学员8dgXkO
[考研] 工科0856求调剂 +3 沐析汀汀 2026-03-21 3/150 2026-03-21 18:30 by 学员8dgXkO
[考研] 0703化学297求调剂 +3 Daisy☆ 2026-03-20 3/150 2026-03-21 17:45 by ColorlessPI
[考研] 266求调剂 +3 哇呼哼呼哼 2026-03-20 3/150 2026-03-21 16:46 by barlinike
[考研] 279求调剂 +5 红衣隐官 2026-03-21 5/250 2026-03-21 14:59 by lature00
[考研] 一志愿天津大学化学工艺专业(081702)315分求调剂 +12 yangfz 2026-03-17 12/600 2026-03-21 03:30 by JourneyLucky
[考研] 一志愿中国石油大学(华东) 本科齐鲁工业大学 +3 石能伟 2026-03-17 3/150 2026-03-21 02:22 by JourneyLucky
[考研] 332求调剂 +4 ydfyh 2026-03-17 4/200 2026-03-21 02:20 by JourneyLucky
[考研] 324分 085600材料化工求调剂 +4 llllkkkhh 2026-03-18 4/200 2026-03-21 01:24 by JourneyLucky
[考研] 321求调剂 +9 何润采123 2026-03-18 11/550 2026-03-20 23:19 by JourneyLucky
[考研] 一志愿苏州大学材料求调剂,总分315(英一) +5 sbdksD 2026-03-19 5/250 2026-03-20 22:10 by luoyongfeng
[考研] 材料学求调剂 +4 Stella_Yao 2026-03-20 4/200 2026-03-20 20:28 by ms629
[考研] 材料工程专硕调剂 +5 204818@lcx 2026-03-17 6/300 2026-03-18 22:55 by 204818@lcx
[考研] 0703化学调剂 +3 妮妮ninicgb 2026-03-17 3/150 2026-03-18 10:29 by macy2011
[考研] 材料,纺织,生物(0856、0710),化学招生啦 +3 Eember. 2026-03-17 9/450 2026-03-18 10:28 by Eember.
[考研] 一志愿,福州大学材料专硕339分求调剂 +3 木子momo青争 2026-03-15 3/150 2026-03-17 07:52 by laoshidan
[考研] 11408 一志愿西电,277分求调剂 +3 zhouzhen654 2026-03-16 3/150 2026-03-17 07:03 by laoshidan
[考研] 304求调剂 +5 素年祭语 2026-03-15 5/250 2026-03-16 17:00 by 我的船我的海
[考研] 070303 总分349求调剂 +3 LJY9966 2026-03-15 5/250 2026-03-16 14:24 by xwxstudy
信息提示
请填处理意见